# A New Adaptive System for Suppression of Transient Wind Disturbance in Large Antennas

Authors: Wei Liang; Jin Huang; Jie Zhang
Journal: International Journal of Antennas and Propagation
 (2019)

---

## Abstract

Under transient wind disturbance, vibration deformation of the large antenna surface profile causes deterioration of pointing performance. This paper presents a new adaptive system to suppress unknown transient wind disturbance. First, to monitor the vibration, based on the acceleration measurement and a low-order flexible model considering equivalent identification of forces, the real-time estimation of the vibration state is obtained in an unbiased minimum-variance way. Next, a novel four-cable-actuator mechanism with a circular slide track is proposed for suppressing the vibration, in which the locations of the cable drivers on the slide track are determined according to the attitude of the antenna, and the expected tension distribution of the cables is found by the vibration state and the optimal gain of the linear quadratic regulator (LQR). In the end, the simulation implementation of a 7.3 m antenna under the transient wind disturbance is used to verify the effectiveness of the proposed method.


---

## Body

## 1. Introduction

Large antennas are widely used in deep space exploration, satellite communications, radio astronomy, and other fields. With the development trend of the large-aperture and high-frequency bands for reflector antennas, a higher pointing accuracy is required under such conditions. For the European Space Agency designed 35 m deep space antenna, a pointing accuracy of 0.006° is required in order to operate it at 35 GHz [1]. As large reflector antennas have a low natural frequency and primarily operate in open air without radome, stochastic wind disturbance obviously increases structural vibration deformation (SVD) which results in pointing performance degradation [2]. To reduce the pointing error (PE), the first step should be to measure the SVD and then to control and compensate for it.There are many measurement methods to get the SVD. But it is very difficult for a large antenna because conventional contact-type transducers require fixed reference points which are rarely found in real situations [3]. Furthermore, the attitude of a large antenna is different during the working process. Although noncontact photogrammetry and laser-based transducers can sense relative displacement [4], these instruments are also incapable of the whole working condition measurement because of light obstruction. On the contrary, accelerometers can satisfy the requirements of the measurements of multiple points without fixed reference points and light obstruction. For reconstructing displacement from acceleration, the frequency-domain integration approach (FDIA) and two integrators with a high-pass digital filter are most commonly used [5]. In references [6, 7], some accelerometers were placed on the antenna reflector to study the dynamic behaviour of the antenna based on the FDIA under the transient wind disturbance. However, the conventional integration approaches need to remove the pseudostatic and low-frequency signals to overcome the integral drift. Moreover, if the FDIA driven by FFT and IFFT is executed in a relatively short time window, the method causes serious discretization errors [5]. Therefore, they are the major drawbacks for complete displacement reconstruction in real time.And to control and compensate for the PE caused by wind disturbance, several control schemes have been developed in large antennas. These methodologies can be roughly divided into two categories. The first category involves using the fixed compensation with the lookup table. For example, Ukita proposed an auxiliary pointing correction system using a lookup table, in which data are compiled from all sky-pointing measurements under different wind conditions. This method was able to reduce the RMS PE from 2.4 to 1.2 arcsec at a static wind speed of 9.3 m/s [6]. However, the method is inconvenient for application because a large amount of data need to be collected, and it only focuses on static wind. The second category is to consider the wind spectrum characteristics in the controller design of the hinge servo system. Using three wind models, Gawronski applied various control theories such as PID, LQG, and H∞ and discussed their basic properties, tracking precision, and limitations [8]. Qiu proposed an extended state observer using the Davenport wind spectrum [9]. These methods have a certain amelioration in suppressing transient wind disturbance. However, the category only adjusts the PE obtained by the angle encoders, the serious SVD-induced PE, that is, the dominant source of the total PE in large antennas is ignored [10]. So a pointing error analytical model was proposed by Zhang to estimate the pointing error caused by the antenna structure deformation [11]. Based on the model, Zhang used different controller design methods such as the model predictive controller [12], the linear-quadratic-Gaussian controller, and the sliding mode controller [13] to compensate for the SVD-induced PE via a rotating shaft servo system. However, the SVD of the antenna is not suppressed. As a result, there is no effective control method to reduce SVD of a large antenna under transient wind disturbance.In response to the above problems, a new adaptive system is proposed to suppress transient wind disturbance in large antennas. This paper is organized as follows: In Section2, the adaptive system is outlined. Next, in Sections 3, a low-order flexible antenna model under transient wind forces and cable tensions in a three-dimensional space is introduced. Section 4 presents a real-time estimator of the vibration state based on acceleration measurements and the low-order flexible model considering equivalent identification of force. The cable tension controller to suppress the vibration based on the LQR is proposed in Section 5. In Section 6, the simulation example is undertaken. Finally, in Section 7, the work is summarized.


## 2. Synthetic Design Scheme with Integration of Mechanical, Sensing, and Control Technologies

The adaptive system is composed of a vibration monitoring and control system, as shown in Figure1. The accelerometers placed on the reflector are used to obtain the vibration state of SVD caused by wind forces and tensions of the cables. And for the convenience of analysis and implementation, the four-cable-actuator mechanism with a circular slide track is designed to suppress unknown transient wind forces; the actual tension is measured by the tension sensor. The state estimator of vibration and tension controller of cables are discussed in Sections 4 and 5, respectively.Figure 1
Schematic diagram of a large antenna under wind disturbance. (a) Original structure. (b) Original structure equipped with the adaptive system.
(a)
(b)As shown in Figures1 and 2, the center of the slide track with radius R2 is located at the center of the antenna base. The drivers can move on the slide track with the change of the azimuth and elevation angles of the antenna. When the azimuth angle is fixed, the “1” and “3” drivers are located at the intersection of the X-axis and the slide track, and the distance between the “2” or “4” driver and the Y-axis is L. Then, the coordinate of the drivers is found, as shown in Figure 2:(1)L=dcosEm+l,where Em is the elevation angle, d is the distance between the elevation axis and the antenna aperture surface, and l is the distance between the azimuth axis and the elevation axis.Figure 2
Geometric relations of various elements.The four tensions act on the edge of the antenna aperture surface with radiusR1. The connection line between action points of T1 and T3 is perpendicular to the elevation axis, and the vectors of T1 and T3 are located in the XZ plane; the connection line between action points of T2 and T4 is parallel to the elevation axis, and the vectors of T2 and T4 are located in a plane perpendicular to the XY plane.The heights of action points ofTi (i=1,2,3,4) are(2)H1=H2−R1cosEm,H2=H4=dsinEm+h,H3=H2+R1cosEm,where h is the height of the elevation axis.The angles betweenTi and negative Z-axis are(3)α1=arctanR2−R1sinEm+LH1,α2=α4=arctanR22−L2−R1H2,α3=arctanR2+L−R1sinEmH3.Components ofTi along X, Y, and Z directions are(4)T1x=T1sinα1;T1y=0;T1z=−T1cosα1,(5)T2x=0;T2y=T2sinα2;T2z=−T2cosα2,(6)T3x=−T3sinα3;T3y=0;T3z=−T3cosα3,(7)T4x=0;T4y=−T4sinα4;T4z=−T4cosα4.


## 3. A Low-Order Flexible Antenna Model under Transient Wind Forces and Cable Tensions in a Three-Dimensional Space

A flexible antenna structural model consisting of numerous nodes and elements has hundreds of thousands of degrees of freedom. The model can be characterized by matrices representing the mass, stiffness, and damping of each component. The inputs of the flexible model are transient wind forces and cable tensions in a three-dimensional space. The outputs of the model are accelerations of measurement points for obtaining the state of the vibration in Section4 and displacements of all reflector nodes for assessing the SVD and PE in Section 6.For describing the flexible vibration of the antenna reflector, consider the continuous-time dynamic model under moving forces in a three-dimensional space:(8)Mp¨xp¨yp¨z+Dp˙xp˙yp˙z+Kpxpypz=Bx00Fx+BTxT1T2T3T4+0By0Fy+BTyT1T2T3T4+00BzFz+BTzT1T2T3T4,(9)axayaz=C0p¨xp¨yp¨z,dxdydz=pxpypz.Here, the antenna reflector hast nodes, and the number of measuring points is m. p=pxTpyTpzTT, p˙=p˙xTp˙yTp˙zTT, and p¨=p¨xTp¨yTp¨zTT are the 3t×1 nodal displacement, velocity, and acceleration vectors of the reflector in the XYZ coordinate system. Also, Fx, Fy, and Fz are all the t×1 input vectors of wind forces which can act on any node of the reflector in X, Y, and Z directions, respectively, as shown in Figure 1; the 4×1 input vector of cable tensions is T=T1T2T3T4T; BTx, BTy, and BTz are the tension decomposition matrices which depend on equations (4)–(7). a=axTayTazTT and d=dxTdyTdzTT are the 3m×1 output vector of measuring points and 3t×1 output vector of nodes. And the 3t×3t matrices M, D, and K are for the masses, damping factors, and stiffness parameters, respectively. The input matrices Bx, By, and Bz all have the dimension t×t. The output matrix C0 has the dimension 3m×3t.For getting a real-time measurement model, a low-order antenna model based on the modal superposition method is introduced. We choose then order modes. To express the structure in modal coordinates, the n×1 modal displacement pm is introduced, which satisfies(10)p=Φpm.Here, the3t×n matrix of mode shapes Φ in X, Y, and Z directions is defined as(11)Φ=ΦxΦyΦz=φx_11...φx_1n.........φx_t1...φx_tnφy_11...φy_1n.........φy_t1...φy_tnφz_11...φz_1n.........φz_t1...φz_tn,where φx_ij is the ith node modal displacement of the jth mode in the X direction.Substituting (11) into (8) and (9), left multiplying (8) by ΦT, and assuming proportional damping, we obtain(12)Mmp¨m+Dmp˙m+Kmpm=ΦxTBxFx+BTxT+ΦyTByFy+BTyT+ΦzTBzFz+BTzT,(13)a=C0Φp¨m,d=Φpm.In this expression, now the modal mass matrixMm, modal stiffness matrix Km, and modal damping matrix Dm satisfy(14)Mm=ΦTMΦ,Km=ΦTKΦ,Dm=ΦTDΦ.Next, left multiplying (12) by Mm−1 yields(15)p¨m+Γp˙m+Λpm=BiF,(16)Bi=Mm−1ΦxTBx,Mm−1ΦyTBy,Mm−1ΦzTBz,(17)F=Fx+BTxTFy+BTyTFz+BTzT,Γ=Mm−1Dm=diag2ξ1ω1,2ξ2ω2,…,2ξnωn,Λ=Mm−1Km=diagω12,ω22,…,ωn2,where ξi is the damping ratio of the ith mode and ωi is the natural frequency of the ith mode.For the sake of convenience, (15) and (13) can be transformed into state-space equations:(18)x˙t=Acxt+BcFt(19)at=Caxt+DcFt,d=Cdxt,(20)Αc=0I−Λ−Γ,Bc=0Bi,(21)Ca=−CoΦΛ−CoΦΓ,Cd=Φ0,Dc=CoΦBi,(22)xt=pmtp˙mt.


## 4. Real-Time Estimation of the Vibration State by Acceleration Measurements and the Low-Order Flexible Model considering Equivalent Identification of Force

First, by discretizing the continuous system of (18) and (19) via a zero-order holder, we obtain(23)xk+1=Axk+BFk+wk,(24)ak=Caxk+DcFk+vk,(25)dk=Cdxk,where A=eAcΔt, in which Δt is the sampling time; B=A−IAc−1Bc; wk is the model error caused by model reduction, and its covariance matrices Q=EwkwlT≥0; vk is the measurement noise, and its covariance matrices R=EvkvlT≥0.However,Fk and xk are unknown in real, so Fk and xk should be identified successively. Gillijns presented an unbiased minimum-variance input and state estimation method for linear discrete-time systems which can solve this problem [14]. The method has been used for input-response estimation by acceleration measurements for structural systems [15, 16]. However, these papers only loaded and identified one or two node forces in one direction. Identifying multisource moving force in a three-dimensional space not only leads to poor computation efficiency but also does not satisfy the direct invertibility requirement in which the number of forces to be identified cannot exceed the number of modes [17]. Since our goal is to estimate xk, it is not necessary to identify the actual multisource moving forces (wind forces and cable tensions) in a three-dimensional space. Therefore, under the premise of stability, observability, controllability, and direct invertibility criteria, the complex problem of identifying unknown spatial distribution forces can be simplified into identifying the equivalent forces of some nodes in a certain direction.The state forcesFkx are the same before and after equivalent:(26)BFk=Be1Fk′=Fkx.For the estimator design, the new state equation of (23) is(27)xk+1=Axk+BFk+wk=Axk+Be1Fk′+wk,Be1=B□BF,B=BX,BY,BZ,BF=B□11B□12...B□1qB□21B□22...B□2q............B□t1B□t1...B□tq,where □ can be a certain direction: X, Y, or Z, and B□∈BX,BY,BZ. BF is the input matrix of force positions. If the jth identified equivalent force is in the X direction at the ith node, then BXij=1 in BF or otherwise BXij=0.The multisource moving forces in a three-dimensional space are(28)Fk=Fk1XFk2X⋯FktXFk1YFk2Y⋯FktYFk1ZFk2Z...FktZT,where FkiX, FkiY, and FkiZ are the ith node forces at step k in X, Y, and Z directions.The theoretical equivalent forces onq nodes are(29)Fk′=Be1+BFk,where Be1+ is the pseudoinverse (pinv) of Be1. And the pinv formed by singular value decomposition not only solves the inverse of the nonsquare matrix but also avoids numerical instabilities due to rank deficiency.The new equation of (24) is(30)ak=Caxk+DcFk+vk=Caxk+De1Fk′+vk,De1=DcB+Be1.Then, based on the acceleration measurement and the new system model comprisingA, Be1, Ca, and De1, the real-time estimation of the identified equivalent forces and vibration state xkk is obtained in an unbiased minimum-variance way; for details on the procedures, readers can refer to the original works [14, 15], and the equations used in the paper are given in Appendix.


## 5. Cable Tension Controller to Suppress the Vibration Based on LQR

For the controller design, the new state equation is(31)xk+1=Axk+BFk+wk=Axk+Fkx+wk,Be2=I.Based on the LQR, the state force that suppresses vibration is obtained from the estimated value of the vibration statexkk:(32)Fk−x=−KLQRxkk.The optimal gainKLQR is derived from minimization of the quadratic cost function:(33)min:JF−x=∑k=0∞xkkTQ1xkk+Fk−xTR1Fk−x,KLQR=Be2TSBe2+R1−1Be2TSA.Here, the weighting parametersQ1 and R1 affected the system output and control input, respectively. Therefore, they certainly satisfied the following conditions: Q1≥0 and R1>0.The infinite horizon solutionS of the associated discrete-time Riccati equation satisfies(34)ATSA−S−ATSBe2Be2TSBe2+R1−1Be2TSA+Q1=0.F
k
−
x should be generated by the expected cable tension of the next time step:(35)Fk−x=BBTxTk+1eBTyTk+1eBTzTk+1e=BBTTk+1e.Then, the expected tension distribution of the cables is found by (32) and (35):(36)Tk+1e=−BBT+KLQRxkk.And a cable-actuator control system (CACS) can be approximately considered a first order of inertia [18]:(37)T˙i=−1τTi+1τTie,where Ti and Tie represent the actual output tension and the expected tension of the ith cable. τ denotes the time constant on inertial constraint units subjected to the element in the whole closed-loop actuator system. Therefore, the servo bandwidth of the CACS is 1/τ.Then, with the sampling timeΔt, by discretizing the continuous cable-actuator control system via a zero-order holder and considering the saturation constraint, we obtain(38)Tk+1=S10000S20000S30000S4T1_k+1T2_k+1T3_k+1T4_k+1=S10000S20000S30000S4e−1/τΔtT1_k+1−e−1/τΔtT1_kee−1/τΔtT2_k+1−e−1/τΔtT2_kee−1/τΔtT3_k+1−e−1/τΔtT3_kee−1/τΔtT4_k+1−e−1/τΔtT4_ke=STe−1/τΔtI4×4Tk+1−e−1/τΔtI4×4Tke,where the unit matrix I4×4 has dimensions of 4×4. The saturation constraint ST used to limit output tension to the upper and lower saturation values is(39)Si=1,forTmax≥Ti≥0,TmaxTi,forTi>Tmax,0,otherwise,where Tmax cannot exceed the allowable stress of the antenna structure and the cable.


## 6. Simulation Examples

The transient wind disturbance rejection simulation of a 7.3 m antenna is used to illustrate the adaptive system. Table1 shows the parameters used for simulation under the largest windward area. The material of the rope is nylon whose Young’s modulus is denoted by E, and the cross-sectional area with a diameter of 1.8 cm is denoted by Ar [19]. Tmax is set to 1000 N. To ensure real-time monitoring and control, the runtime of each recursive step in the adaptive system should be smaller than the sampling time. We use MATLAB commands “tic” and “toc” to test each recursive step runtime. In an ordinary computer (with the Intel Core i5-6500 CPU and 16 GB RAM), each recursive time step takes about 0.7 ms. So the sampling time is set to 1 ms which is larger than 0.7 ms.Table 1
Parameters used for simulation.
Value
Parameters of the adaptive system
E
m
∘
0
h
m
3.8915
d
m
2.4445
l
m
0.4150
R
1
m
3.65
R
2
m
10
E
Gpa
2.5
A
r
m
2
2
.
54
×
10
−
4
L
m
2.8595
H
1
m
0.2415
H
2
m
3.8915
H
3
m
7.5415
α
1
∘
88.0629
α
2
∘
56.7364
α
3
∘
59.6104
α
4
∘
56.7364
Other parameters
Wind speed (m/s)
v
Wind direction
Front of the antenna
Sampling time (s)
0.001
Simulation time (s)
4


### 6.1. Transient Wind Model of a 7.3 m Antenna

Wind load can be classified into mean wind and gust wind. The wind speed at any point can be expressed as the sum of mean wind speed and gust wind speed [20, 21]:(40)v=vZbZZbαs+Δv,where mean wind satisfies the exponential law relation of the wind profile, vZb is the average wind speed at the reference height of Zb in the XYZ coordinate system and set to 10 m/s, Zb is set to 1 m, and αs is the ground roughness set to 0.16 for open terrain. And Δv can be obtained by the Davenport filter from white noise with an average value of zero.We set up the wind field of the antenna in the CFX software, import the wind speed data, and perform transient analysis, as shown in Figure3. After the analysis, the transient wind forces can be obtained.Figure 3
Wind field of the antenna.


### 6.2. A Low-Order Flexible Antenna Model for Estimator and Controller Design

In order to get an approximate full-order model (AFM) for real-time monitoring and control simulation in the next section and get the low-order flexible model for estimator and controller design, the modal superposition method solved by MATLAB is introduced as described in Section3. Since the characteristics of the wind spectrum are mainly associated with low frequency, resulting in the PE in low frequency [7], we choose the AFM consisting of the first 26 modes within 50 Hz. By loading the same wind forces, setting the same sampling rate and simulation time, the displacements of the reflector surface are obtained by transient analysis of the full-order model using ANSYS and the first 26-order antenna model with MATLAB. Then, according to the best-fitting paraboloid theory [10, 22], we get the PE of the reflector surface in the UVW coordinate system:(41)θ=θi,θiiT=fdW,where θi and θii are the PEs in the UW plane and VW plane, respectively, as shown in Figure 4. And if the reflector does not have deformation, then θT=0. dW obtained through coordinate transformation is the displacement vector of t nodes in the W direction.Figure 4
PEs in theUVW coordinate system.The amplitude ofθ is(42)θT=θi2+θii2.As shown in Figure5, the PE of the first 26 modes is basically the same as the PE of the full-order antenna model, with a relative error of 0.79%. Therefore, the AFM can be used to approximate the full-order antenna model. Moreover, transient wind mainly excites the fundamental frequency mode of the elevation direction and the azimuth direction for the antenna, as shown in Figures 5 and 6. As a result, the first 2-mode antenna model within 10 Hz is used as the low-order flexible model.Figure 5
PEs under wind forces in the (a) time domain and (b) frequency domain.
(a)
(b)Figure 6
The first two modes of the antenna in the (a) elevation direction and (b) azimuth direction.
(a)
(b)


### 6.3. Adaptive System to Suppress Transient Wind Disturbance

As shown in Figure7, the antenna structure consists of the original antenna (the AFM is used in simulation) and the four-cable-actuator mechanism with a circular slide track. For the calculation formula of Mk and Kk, see Appendix. The initial state of the system without considering the effects of temperature and gravity is assumed to be zero, and the covariance matrices Q and R are assigned values 1e−10 and 1e−6 on the diagonal, respectively. The data of acceleration measurements are polluted by white noise with an average value of zero and a maximum amplitude of 1.2×10−3m/s2.Figure 7
Adaptive system to suppress transient wind disturbance.A network of structure nodes, sensors, and action points of identified equivalent forces (IEFs) and tensions is shown in Figure8. The sensor placement issue is not discussed here (refer to [20, 23] for sensor placement strategies). For the choice of equivalent force positions, results remain sensitive to correct assumptions on the force locations [24]; they should be determined according to the actual identification effect.Figure 8
A network of structure nodes, sensors, and action points of IEFs and tensions.For the value ofτ, the servo bandwidth of the CACS should be at least larger than the main frequency of structural vibration. Hence, a value of τ is set to 1/16π; the servo bandwidth (8 Hz) corresponding to the τ is higher than the second mode frequency (6.6 Hz). And we think that if the servo bandwidth can cover the frequency range of larger amplitude in Ti under the absence of inertial constraint in the CACS, the better control effect will be achieved. Therefore, after the adaptive system is operated without considering the inertial constraint in the CACS, the spectral characteristics of the ideal Ti can be obtained, as shown in Figure 9. Since the amplitude spectrum of the ideal Ti is mainly distributed within 80 Hz, another value of τ is set to 1/160π. Figure 10 shows the Bode diagram of the CACS under two different τ; the frequency corresponding to −3dB is the servo bandwidth of the CACS.Figure 9
FFT ofTi without considering the inertial constraint in the CACS.Figure 10
Bode diagram of the CACS under two differentτ.
(a)
(b)As shown in Figure11, the IEFs of the two nodes in the X direction under τ=1/160π are in close agreement with the theoretical equivalent forces which are solved by (29). Similar good recognition results are also obtained under the absence of inertial constraint and τ=1/16π.Figure 11
Equivalent forcesF1x′ (a) and F2x′ (b) under τ=1/160π.
(a)
(b)To suppress the SVD under the transient wind disturbance, tensions of the cables can be adaptively controlled, as shown in Figure12. With the higher servo bandwidth of the CACS under τ=1/160π, the Ti is very close to the ideal Ti; with the lower servo bandwidth of the CACS under τ=1/16π, the phase of the Ti has a certain lag with respect to the ideal Ti. The longitudinal deformation of the cable, as shown in Figure 12, is obtained by the following formula in which the cable is assumed as a linear spring [25]:(43)ΔLi=Tikci,kci=EArLi,where Li is the length of the cable.Figure 12
Tensions and longitudinal deformations of the cables. The red curve is the idealTi without considering the inertial constraint in the CACS; the blue curve is the Ti under τ=1/160π; the green curve is the Ti under τ=1/16π.
(a)
(b)
(c)
(d)Therefore, theTi under τ=1/160π has a better control effect on SVD and PE compared to the Ti under τ=1/16π, as shown in Figures 13 and 14. In the end, the RMS PEs with the adaptive system under the absence of inertial constraint, τ=1/160π and τ=1/16π, are reduced by 80.26%, 79.44%, and 66.45% compared with the uncontrolled situation, respectively.Figure 13
SVD inW direction at 3 s.
(a)
(b)
(c)Figure 14
PEs without control and with control.


## 6.1. Transient Wind Model of a 7.3 m Antenna

Wind load can be classified into mean wind and gust wind. The wind speed at any point can be expressed as the sum of mean wind speed and gust wind speed [20, 21]:(40)v=vZbZZbαs+Δv,where mean wind satisfies the exponential law relation of the wind profile, vZb is the average wind speed at the reference height of Zb in the XYZ coordinate system and set to 10 m/s, Zb is set to 1 m, and αs is the ground roughness set to 0.16 for open terrain. And Δv can be obtained by the Davenport filter from white noise with an average value of zero.We set up the wind field of the antenna in the CFX software, import the wind speed data, and perform transient analysis, as shown in Figure3. After the analysis, the transient wind forces can be obtained.Figure 3
Wind field of the antenna.


## 6.2. A Low-Order Flexible Antenna Model for Estimator and Controller Design

In order to get an approximate full-order model (AFM) for real-time monitoring and control simulation in the next section and get the low-order flexible model for estimator and controller design, the modal superposition method solved by MATLAB is introduced as described in Section3. Since the characteristics of the wind spectrum are mainly associated with low frequency, resulting in the PE in low frequency [7], we choose the AFM consisting of the first 26 modes within 50 Hz. By loading the same wind forces, setting the same sampling rate and simulation time, the displacements of the reflector surface are obtained by transient analysis of the full-order model using ANSYS and the first 26-order antenna model with MATLAB. Then, according to the best-fitting paraboloid theory [10, 22], we get the PE of the reflector surface in the UVW coordinate system:(41)θ=θi,θiiT=fdW,where θi and θii are the PEs in the UW plane and VW plane, respectively, as shown in Figure 4. And if the reflector does not have deformation, then θT=0. dW obtained through coordinate transformation is the displacement vector of t nodes in the W direction.Figure 4
PEs in theUVW coordinate system.The amplitude ofθ is(42)θT=θi2+θii2.As shown in Figure5, the PE of the first 26 modes is basically the same as the PE of the full-order antenna model, with a relative error of 0.79%. Therefore, the AFM can be used to approximate the full-order antenna model. Moreover, transient wind mainly excites the fundamental frequency mode of the elevation direction and the azimuth direction for the antenna, as shown in Figures 5 and 6. As a result, the first 2-mode antenna model within 10 Hz is used as the low-order flexible model.Figure 5
PEs under wind forces in the (a) time domain and (b) frequency domain.
(a)
(b)Figure 6
The first two modes of the antenna in the (a) elevation direction and (b) azimuth direction.
(a)
(b)


## 6.3. Adaptive System to Suppress Transient Wind Disturbance

As shown in Figure7, the antenna structure consists of the original antenna (the AFM is used in simulation) and the four-cable-actuator mechanism with a circular slide track. For the calculation formula of Mk and Kk, see Appendix. The initial state of the system without considering the effects of temperature and gravity is assumed to be zero, and the covariance matrices Q and R are assigned values 1e−10 and 1e−6 on the diagonal, respectively. The data of acceleration measurements are polluted by white noise with an average value of zero and a maximum amplitude of 1.2×10−3m/s2.Figure 7
Adaptive system to suppress transient wind disturbance.A network of structure nodes, sensors, and action points of identified equivalent forces (IEFs) and tensions is shown in Figure8. The sensor placement issue is not discussed here (refer to [20, 23] for sensor placement strategies). For the choice of equivalent force positions, results remain sensitive to correct assumptions on the force locations [24]; they should be determined according to the actual identification effect.Figure 8
A network of structure nodes, sensors, and action points of IEFs and tensions.For the value ofτ, the servo bandwidth of the CACS should be at least larger than the main frequency of structural vibration. Hence, a value of τ is set to 1/16π; the servo bandwidth (8 Hz) corresponding to the τ is higher than the second mode frequency (6.6 Hz). And we think that if the servo bandwidth can cover the frequency range of larger amplitude in Ti under the absence of inertial constraint in the CACS, the better control effect will be achieved. Therefore, after the adaptive system is operated without considering the inertial constraint in the CACS, the spectral characteristics of the ideal Ti can be obtained, as shown in Figure 9. Since the amplitude spectrum of the ideal Ti is mainly distributed within 80 Hz, another value of τ is set to 1/160π. Figure 10 shows the Bode diagram of the CACS under two different τ; the frequency corresponding to −3dB is the servo bandwidth of the CACS.Figure 9
FFT ofTi without considering the inertial constraint in the CACS.Figure 10
Bode diagram of the CACS under two differentτ.
(a)
(b)As shown in Figure11, the IEFs of the two nodes in the X direction under τ=1/160π are in close agreement with the theoretical equivalent forces which are solved by (29). Similar good recognition results are also obtained under the absence of inertial constraint and τ=1/16π.Figure 11
Equivalent forcesF1x′ (a) and F2x′ (b) under τ=1/160π.
(a)
(b)To suppress the SVD under the transient wind disturbance, tensions of the cables can be adaptively controlled, as shown in Figure12. With the higher servo bandwidth of the CACS under τ=1/160π, the Ti is very close to the ideal Ti; with the lower servo bandwidth of the CACS under τ=1/16π, the phase of the Ti has a certain lag with respect to the ideal Ti. The longitudinal deformation of the cable, as shown in Figure 12, is obtained by the following formula in which the cable is assumed as a linear spring [25]:(43)ΔLi=Tikci,kci=EArLi,where Li is the length of the cable.Figure 12
Tensions and longitudinal deformations of the cables. The red curve is the idealTi without considering the inertial constraint in the CACS; the blue curve is the Ti under τ=1/160π; the green curve is the Ti under τ=1/16π.
(a)
(b)
(c)
(d)Therefore, theTi under τ=1/160π has a better control effect on SVD and PE compared to the Ti under τ=1/16π, as shown in Figures 13 and 14. In the end, the RMS PEs with the adaptive system under the absence of inertial constraint, τ=1/160π and τ=1/16π, are reduced by 80.26%, 79.44%, and 66.45% compared with the uncontrolled situation, respectively.Figure 13
SVD inW direction at 3 s.
(a)
(b)
(c)Figure 14
PEs without control and with control.


## 7. Conclusion

In this paper, a novel antenna monitoring and control system is proposed to suppress transient wind disturbance. From the initial study, the following conclusions are drawn:(1)
In the new adaptive system, the accelerometers placed on the reflector are used to obtain the vibration state of SVD caused by wind forces and tensions of the cables; the four-cable-actuator mechanism with a circular slide track is designed to suppress unknown transient wind forces.(2)
To monitor the vibration, the first equivalent transformation is carried out to simplify the complex problem of identifying unknown spatial distribution forces including wind forces and cable tensions into identifying the equivalent forces of some nodes; based on the acceleration measurement and a low-order model, the identified equivalent forces and the vibration state are estimated successively in an unbiased minimum-variance way.(3)
To control the vibration, the second equivalent transformation is carried out to get cable tensions from the state forces for suppressing the vibration; the state forces are determined by the vibration state and the optimal gain of the LQR.(4)
Since the wind conditions are ever changing in the real environment, the loads acting on the antenna are different under different wind directions or speeds. To further improve control performance under random wind disturbance, how to optimize the system parameters and enhance the system robustness are key problems to be addressed in the future work.(5)
This method is currently applicable to the situation where rotating shaft of antenna are locked, such as the communication between the antenna and the geostationary spacecraft. In order to achieve an accurate tracking motion of the antenna, flexible multibody dynamics research of the system still needs to be done.



---

*Source: 2015341-2019-09-04.xml*