Articles | Volume 12, issue 1
https://doi.org/10.5194/ms-12-59-2021
https://doi.org/10.5194/ms-12-59-2021
Research article
 | 
04 Feb 2021
Research article |  | 04 Feb 2021

Analysis and control of state jump in space deployable structures under alternating temperature loads

Congcong Chen, Tuanjie Li, Yaqiong Tang, and Zuowei Wang
Abstract

State jump has been experimentally observed in space deployable structures working in alternating temperature environments. State jump is a phenomenon in which the geometric shape of the structure changes after the temperature loading and unloading process, which makes the working accuracy of the space deployable structure intrinsically unpredictable. This paper aims to investigate the causes of this state jump phenomenon and seek measures to reduce its effect. Firstly, the static multiple-stable-state phenomenon resulting in state jump is analyzed for clearance joints in deployable structures. Then, an equivalent model consisting of a variable stiffness spring and a contact element for state jump analysis is proposed, which is verified by a finite element simulation. Influence factors and control methods of state jump are further explored. Finally, numerical results of a space deployable structure of an umbrella-shaped antenna show the effectiveness of the developed analytical method.

1 Introduction

As the carrier of receiving and reflecting electromagnetic signals, the satellite antenna is a key piece of equipment for deep space explorations, satellite navigation, satellite communications and electronic reconnaissance. In order to achieve the ability of long-distance information transmission and capture weak signals, the satellite antenna is required to have a large aperture (Alfred and Joseph, 1998; Arnol and Naderi, 1989; Rahmat-Samii and Densmore, 2015). At the same time, due to the limitation of launch capability and space, large satellite antennas must be light and deployable. For this reason, space deployable structures have been widely used in constructing large satellite antennas.

The space deployable structure (Gantes, 2001) is a special kind of deployable mechanism, which can transform from a folding configuration with a small volume to a stable load-bearing structure configuration with a large volume or surface area. It is found in the ground test that even if temperature loads applied on space structures are identical, the structural state is inconsistent. Moreover, affected by solar radiation and earth shadow, space temperature is alternating from −180 to 150 C (Hu et al., 2013). These reasons make the working accuracy of space deployable structures intrinsically unpredictable.

So far, researches about space deployable structures under temperature loads are mainly focused on analysis and control of thermal deformations and thermally induced vibrations. For example, Tang et al. (2014) studied thermal deformation and surface accuracy of large deployable antennas in orbit based on the stochastic finite element method. Lu et al. (2019) studied thermal deformation and the shape adjustment method of a planar-phased array antenna structure. Thornton and Kim (1993) studied thermally induced vibrations of the cantilever beam of the Hubble Space Telescope and established the stability criterion of thermally induced vibrations. Shen et al. (2013) established thermal-structural analysis model of flexible beams using the absolute nodal coordinate formulation and studied its thermally induced vibrations. Azadi et al. (2017) studied thermally induced vibrations of smart solar panels and analyzed the influence of heat radiation and orbit parameters on the amplitude of thermal vibrations. Fazelzadeh and Azadi (2017) modeled the orbiting smart satellite panels as a functionally graded material beam and studied the control strategy of thermally induced vibrations.

While most of the above studies are mainly based on ideal structures, some other studies have shown that the thermal response of the structure is strongly affected by joint clearances in structure, which is inevitable in space deployable structures due to manufacturing errors and assembly requirements. Kim et al. (1999) studied the thermal-creak-induced dynamics caused by external temperature changes of space structures with clearances, but the influence of temperature changes on position accuracy and stability of the structure is not considered. With the urgent need for high-accuracy space deployable structures, the influence of clearance on the accuracy and stability of deployable structures cannot be neglected. To fill this research gap and meet the needs of high-precision and high-stability deployable structures in aerospace industry, further study of the state jump phenomenon is needed for space deployable structures with clearance joints.

This paper is organized as follows: first, kinematic joints of space deployable structures are classified and the static multiple-stable-state phenomenon in clearance joints is analyzed in Sect. 2. Then, the state jump analysis model is established and verified, and interfering factors and control methods of the state jump phenomenon are studied in Sect. 3. In Sect. 4, the above contents are verified by an example of an umbrella antenna. And conclusions are summarized in Sect. 5.

2 Static multiple-stable-state phenomenon in clearance joints

2.1 Classification of clearance joints

Ideally, the degree of freedom of a deployable structure is 0 when it is in deployment state. However, the existence of clearances in kinematic joints makes the structure movable to some extent.

Kinematic joints in space deployable structures are divided into three types: one-dimensional clearance joint, two-dimensional clearance joint and three-dimensional clearance joint. The one-dimensional, two-dimensional and three-dimensional clearance joints allows small linear relative displacement, small planar displacement and small spatial relative displacement between the two components connected by the joint, respectively. Figure 1 shows these three types of clearance joints, in which I and II are two components connected by the clearance joint, and c is the size of the clearance. For illustration, the clearance in Fig. 1 is enlarged.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f01

Figure 1Classification of clearance joints.

Download

2.2 Static multiple-stable-state phenomenon

Without losing the generality, it is assumed that the pairing element of component II is a second-order continuous surface S=Γx,y,z. When component II is fixed and a spatial external force F=(Fx,Fy,Fz) at the centroid of component I is applied, component I will eventually be balanced under the combined action of the external force F and reaction force Fr provided by component II, that is

(1) F + F r = F + F N + F T = 0 ,

where FN, FT are normal and tangential components of reaction force Fr, respectively. According to the force equilibrium condition, FN, FT and F are coplanar.

For an arbitrary point (x,y,z) on S, Eq. (1) has an exclusive resolution (FN0,FT0). If and only if FN is positive (pointing from the component II to the component I) and |FT|≤μ|FN|, the component I can be balanced at this point, where μ is the friction coefficient of the contact surface. When the above conditions are not satisfied, component I cannot be balanced at this point.

By solving Eq. (1), the following can be found:

  1. When there is no friction in the clearance joint (μ=0), component I under a given external force F can be balanced at only one point, as shown in Fig. 2. OI and OII are the centers of component I and component II, respectively.

  2. When μ≠0, component I could be balanced in multiple positions, as shown in Fig. 3. The filled circles in Fig. 3 represent possible equilibrium positions of component I. For two components connected by a one-dimensional, two-dimensional or three-dimensional clearance joint, their relative equilibrium positions can be expressed as a straight line, a spatial curve or a spatial surface, respectively.

In summary, due to the existence of clearance and friction in deployable structure joints, the static equilibrium position of space deployable structure components under external force is not unique. In other words, there is a static multiple-stable-state phenomenon. This phenomenon makes the state jump occur in space deployable structures.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f02

Figure 2Static equilibrium position of frictionless clearance joints.

Download

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f03

Figure 3Static equilibrium positions of friction clearance joints.

Download

3 State jump phenomenon in clearance joints

3.1 Analysis model of the state jump phenomenon

The analysis model of deployable structure component used in this paper is shown in Fig. 4. The structure component is composed of a massless spring and a lumped mass block, and it may deform axially under the temperature loads. In order to distinguish intrinsic factors related to the state jump phenomenon, the quasi-static assumption is introduced; that is, the accumulation and release of elastic energy in the component is instantly completed.

When the temperature change is ΔT, the length change of the component is

(2) Δ l = α L Δ T l ,

where Δl, l and αL denote the thermal deformation, original length and coefficient of linear thermal expansion of the component, respectively.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f04

Figure 4State jump analysis model.

Download

According to the Hooke's law, the relationship between the axial force and the section stress σ and deformation Δl of the flexible member is

(3) F = σ A = E A l Δ l ,

where E and A are the elastic modulus and cross-sectional area of the component, respectively.

When the axial force in the component is insufficient to overcome the maximum friction force provided by the contact surface, the length of the component will remain the same and the energy of the external temperature load will be accumulated in the form of elastic energy. If the difference between the actual deformation and thermal deformation of the component is Δx, the axial force in the component can be calculated as

(4) F = E A l t Δ x = k t Δ x ,

where lt is the stress-free length of the component at time t. Therefore, the non-linear spring is used to establish the equivalent mechanical model of the structure component.

3.2 State jump analysis

The temperature applied on the structure is assumed to be T0 at the beginning, and after a certain loading process, the temperature returns to T0 for the first time.

When there is no friction in the clearance joint and the clearance size is large enough, the component will expand and contract freely under external temperature loads. When the temperatures before and after loading are the same, the structural states will also be the same; in other words, the lumped mass block will return to its original position at the end of the loading process.

When there is no friction in the clearance joint and the clearance size is small, elastic energy will be accumulated in the component when the length variation of the component caused by temperature is greater than the clearance size. Once the temperature changes inversely, the energy accumulated in the component will be released, and there is no energy dissipation in the process when the complete elastic deformation is considered. Therefore, the lumped mass block will still return to its original position after unloading.

When there is friction in the clearance joint, energy dissipation and transfer will occur during the loading process, which makes it impossible for the lumped mass block to return to its original position. There exists the state jump phenomenon.

As shown in Fig. 5, when the clearance size is large enough, the system will go through five typical stages. In the figure, xt and xT are actual displacement (deformation) and theoretical thermal displacement (deformation) of the component, respectively.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f05

Figure 5Typical thermal loading process.

Download

3.2.1 Stage S1: static stage 1

In this stage, the stress caused by temperature change is not enough to overcome the friction force, and the actual length of the component remains unchanged. The component is subjected to an axial force Ft(|Ft|≤f), where f is the maximum static friction provided by the contact surface. The actual deformation of the component in this stage is xt=0.

3.2.2 Stage S2: sliding stage 1

In this stage, the insufficient friction cannot keep the component maintaining its current length, and the component begins to deform. The axial force in the component is |Ft|=f. The actual deformation of the component is xt=xT-sign(dxTdt)ltfEA, where sign(A) is the sign function.

3.2.3 Stage S3: static stage 2

In this stage, the external environment temperature begins to change reversely, and the elastic energy stored in the component is released. The axial force in the component is |Ft|≤f. The actual deformation of the component in this stage is xt=xt2.

3.2.4 Stage S4: static stage 3

Similar to S1, in this stage, the actual length of the component remains unchanged, and the axial force in the component is |Ft|≤f. The actual deformation of the component in this stage is xt=xt2.

3.2.5 Stage S5: sliding stage 2

Similar to S2, the friction is insufficient and the component begins to deform reversely. The axial force in the component is |Ft|=f. The actual deformation of component I is xt=xT+sign(dxTdt)ltfEA.

3.3 Validation of the analysis model

Geometric and material parameters of component I are shown in Table 1, and the friction force is 30 N. The state jump value Det can be calculate as Det=ltfEA=0.5×30709×9-6×1000mm=0.024mm, as discussed in Sect. 3.2.

Table 1Geometric and material parameters of component I.

Download Print Version | Download XLSX

The finite element simulation model is established and analyzed using the ANSYS software to verify the effectiveness of the above analysis model. As shown in Fig. 6, a BEAM188 element and a SOLID185 element are chosen to simulate component I and component II, respectively. The contact pairs between the lower surface of component I and the upper surface of component II are established using a CONTA170 element and a TARGET169 element, respectively. Two CONTA178 elements with an initial gap size of c are used to simulate the horizontal contact between the components. A vertical force Fn is applied to keep the components in contact state.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f06

Figure 6Schematic diagram of finite element model.

Download

In the simulation, component II is regarded as the grid, the clearance size is 1 mm, the friction coefficient between the components is 0.3, the vertical force is −100 N (the corresponding friction is 30 N), the variation of the environmental temperature is 100 K and the temperature changes in a sinusoidal law.

The displacement and force curves of component I by simulation analysis are shown in Fig. 7.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f07

Figure 7Displacement curve of component I.

Download

The color cloud picture of displacement of component I is shown in Fig. 8. The maximum deformation at the right of the component is 0.0244 mm, which is consistent with the result of the equivalent analysis model.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f08

Figure 8Color cloud picture of displacement.

Download

3.4 Control of the state jump phenomenon

It can be seen from Sect. 3.2 that, for two components connected by a one-dimensional clearance joint, the deformation of the component after being loaded is ltfEA, which means the state jump value is independent of the motion path of the component. The main factors affecting the state jump value are the magnitude of friction and the axial stiffness of the component, while the maximum thermal deformation and clearance mainly affect the shape of the structure during the temperature loading process. Therefore, the following measures may be taken to control the state jump phenomenon.

  1. Control the friction between components. The main factors affecting the friction between components are the contact pressure, friction type, material type and roughness of the contact surface. Therefore, the control of the state jump phenomenon can be realized by choosing and controlling component materials, processing technology and components matching types.

  2. Control the deformation stiffness of components. The main factor affecting the deformation stiffness of components is the material of the component. Besides, the deformation stiffness of the component can also be increased by setting preload springs in the clearance joint. Therefore, the control of the state jump phenomenon can be realized by reasonably selecting the materials of the components and adding preload springs in kinematic joints.

The following should be noted:

  1. Applications of above control measures may affect the deployment process or other performances of the structure, which should be considered comprehensively.

  2. With the increase in the dimension and number of clearance joints, the affecting factors may be different and the control of the state jump phenomenon may become much more complicated.

4 Simulation analysis

4.1 Modeling of deployable structure of umbrella antenna

The equivalent modeling method proposed in Sect. 3 is used to analyze the deployable structure of an umbrella-shaped antenna. The deployment principle of the umbrella-shaped deployable antenna is shown in Fig. 9. The motor drives the slider D to move to realize the closure and expansion of the antenna.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f09

Figure 9Deployment principle of umbrella antenna.

Download

In Fig. 9, link AB is fixed to the antenna rib, A is a joint at the root of antenna rib, B is located at an eccentric position of the rib root, BC is the connecting rod and CD is the sliding disk. The dotted line in Fig. 9 represents the folded state of the antenna while the solid line represents the unfolded state. After being deployed and locked in place, the mechanism degenerates into a structure and AC is equivalent to the frame. The equivalent deployable structure is shown in Fig. 10. XG1, YG1, XG2 and YG2 are the centroid coordinates of links AB and BC. θ2 and θ4 are the angles between the links and the horizontal direction.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f10

Figure 10Deployable structure with clearance joints.

Download

Dynamic equations of the structure are established using the Lagrange method,

(5) d d t ( T - V ) q j - ( T - V ) q j = Q j ( j = 1 , 2 , , 6 ) ,

where q=(XG1,YG1,XG2,YG2,θ2,θ4) denotes the generalized coordinates, T is the kinetic energy, V is the potential energy of the system and Qj is the generalized force corresponding to qj. The generalized force Qj is a function of contact forces and friction forces in clearance joints. The classical Hertz contact force model (Lankarani and Nikravesh, 1990) and Ambrósio friction model (Ambrósio, 2003) are introduced to calculate Qj.

The classical Hertz contact force model is

(6)FN=Kδn,(7)K=43πh1+h2R1R2R1+R21/2,(8)h1=1-ρ12πE1,h2=1-ρ22πE2,

where K is the stiffness coefficient, δ is the penetration depth, and n is the force index which depends on material and geometric properties. For metal contact, n is generally 1.5. R1 and R2 are the radii of two contact bodies, and the material coefficient hi depends on Young's modulus Ei and Poisson's ratio ρi of the contact body.

The Coulomb friction model is the most classical friction model but it does not take the tangential velocity between the two contact components into account. When the tangential velocity is in the vicinity of zero, the value of friction force is changed from Ft to Ft, leading to numerical solving difficulties. To address these problems, Ambrósio proposed a modified Coulomb friction model (Ambrósio, 2003):

(9)FT=μfμdFN,(10)μd=0vtv0vt-v0v1-v0v0vtv11vtv1,

where μf is sliding friction coefficient, μd is dynamic correction coefficient, vt is the tangential velocity, and v0 and v1 are given tolerances for the tangential velocity.

In order to obtain the contact force and friction force in clearance joints, the penetration depth and tangential contact velocity should be calculated. The penetration in a clearance joint is shown in Fig. 11.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f11

Figure 11Penetration in clearance joint.

Download

In Fig. 11, C and D are the centers of the journal and the bearing, respectively, and rC and rD are coordinate vectors of points C and D. The eccentricity vector dDC which connects the centers of the bearing and the journal is expressed as

(11) d DC = r C - r D ,

Then the penetration depth in the clearance joint can be expressed as

(12) δ DC = d DC - ( R 2 - R 1 ) ,

where R1 and R2 are the radii of the journal and the bearing, respectively.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f12

Figure 12Flow chart of simulation solution.

Download

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f13

Figure 13Center trajectories of the journal and bearing.

Download

The unit normal vector at the contact point is expressed as

(13) n DC = d DC d DC ,

Then the coordinate vectors of contact points C1 and D1 are calculated by

(14) r C 1 = r C + R 1 n DC r D 1 = r D + R 2 n DC ,

The relative penetration velocity between the journal and bearing can be expressed as

(15) v D 1 C 1 = r ˙ C 1 - r ˙ D 1 .

The normal contact velocity and the tangential contact velocity between the journal and the bearing are expressed as

(16) v n = v D 1 C 1 n DC n DC v t = v D 1 C 1 - v n .

Contact force FN and friction force FT now can be obtained by using Eqs. (6) and (9).

The change of temperature results in the change of component stiffness, and correspondingly; the stress-free length of the component changes. When the component length is varied with time or temperature, the center positions of the journal and the bearing change, resulting in the variation of the contact force and friction force.

4.2 Computational strategy of dynamic equations

The fourth-order Runge–Kutta method is used to solve dynamic equations, and the flow chart of the solving process is shown in Fig. 12.

  1. Give the total number of iterations as N and set the current iteration number as n=0.

  2. Update the generalized displacement as qn, the generalized velocity q˙n and the component length ln.

  3. Calculate the penetration depth δ using Eq. (12).

  4. If the journal and bearing are in contact (δ>0), calculate FN and FT using Eqs. (6) and (9), otherwise, set FN=FT=0.

  5. Calculate the generalized displacement qn+1 and generalized velocity q˙n+1 using Eq. (5).

  6. If n<N, set n=n+1 and return to step (2), otherwise, end the simulation.

4.3 State jump and control in deployable structure

For the above-mentioned clearance structure, geometric and material parameters are shown in Table 2 and simulation parameters are shown in Table 3. At the initial moment, the center of the journal and the bearing are coincident, and the generalized velocity and acceleration are set as 0. The joints between the links and the frame are assumed to be ideal and the clearance size of the joint between AB and BC is 0.1 mm. The temperature increases linearly from 0 to 100 K in the first half of loading and then decreases linearly to 0 K.

Table 2Geometric and material parameters of deployable structure.

Download Print Version | Download XLSX

Table 3Simulation parameters of dynamic solution process.

Download Print Version | Download XLSX

Simulation results are shown in Fig. 13. Curves in Fig. 13a and b are the center trajectories of the journal and the bearing, respectively, and Fig. 13c shows the trajectory of the bearing center relative to the journal center. As can be seen from Fig. 13, relative positions of deployable structural members will change while the temperatures before and after loading are the same; that is, the state jump phenomenon occurs.

Deflection angle curves of AB under different clearance sizes are shown in Fig. 14. It is shown that the deflection angle of AB varies with the clearance size. The smaller the clearance is, the earlier the deflection occurs. When the clearance size is large enough (greater than the maximal thermal deformation), AB and BC do not come into contact during the whole loading process.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f14

Figure 14Deflection angle of AB under different clearance sizes.

Download

Deflection angle curves of link AB under different friction coefficients are shown in Fig. 15. It is shown that the deflection angle of AB varies with the friction coefficients. Generally speaking, the larger the friction coefficient, the larger the range of the deflection angle.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f15

Figure 15Deflection angle of AB under different friction coefficients.

Download

One torsion spring is added to the clearance joint of the deployable structure, and its effect on the state jump phenomenon is analyzed. The torque M provided by the torsion spring with an angle α is

(17) M = K M Δ σ = K M ( α - α 0 ) ,

where α is the flare angle, α0 is the flare angle of torsion spring in the free state and KM is the elastic coefficient of torsional spring. The effect of torsional spring on the structure is equivalent to the external force varying with the flare angle.

The torsion spring is assumed to be in the stress-free state when the structure is in ideal position (the center of the journal coincides with that of the bearing), and the deflection angle of AB under different elastic coefficients is shown in Fig. 16.

https://ms.copernicus.org/articles/12/59/2021/ms-12-59-2021-f16

Figure 16Deflection angle of AB under different elastic coefficients.

Download

It can be seen from Fig. 16 that the state jump value will change when the torsion spring is installed in the clearance joint. The deflection angle decreases with the increase in the elastic coefficient of the torsion spring, in a certain range. Therefore, by using reasonable methods, the state jump phenomenon in deployable structures can be controlled.

5 Conclusion

In this paper, the static multiple-stable-state phenomenon in clearance joints of deployable structures is clarified through force analysis. An equivalent model consisting of a variable stiffness spring and a contact element is proposed, which can capture the state jump phenomenon in deployable structures. The control methods have been presented for reducing the effect of state jump on structural responses. Numerical results of the deployable structure of an umbrella-shaped antenna have verified the proposed analysis method. Some conclusions have been summarized in the following.

  • a.

    Clearance and friction in kinematic joints are the basic reasons for the generation of the state jump phenomenon.

  • b.

    The magnitude of friction and axial stiffness of components are the main factors affecting the state jump value.

  • c.

    The state jump phenomenon can be controlled by reasonable measures, such as reducing the friction force, increasing the stiffness of structure components and adding pre-stressed springs to the structure.

6 Future work

Future work will focus on improving the performance of the space deployable structure under complex temperature loads.

Code availability

The code is available upon request by contact with the corresponding author.

Data availability

All data included in this study are available upon request by contact with the corresponding author.

Author contributions

CC did all of the work for this research under the supervision of TL and YT, and ZW improved the overall writing and language.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Robotics and advanced manufacturing”. It is not associated with a conference.

Financial support

This research has been supported by the National Natural Science Foundation of China (grant nos. 51775403 and 51905401) and the Natural Science Foundation of Shaanxi Province (grant no. 2019JQ-505).

Review statement

This paper was edited by Guimin Chen and reviewed by two anonymous referees.

References

Alfred, U. M. R. and Joseph, N. P.: Global satellite communications technology and systems – An overview, Space Commun., 16, 55–69, https://doi.org/10.1109/RADAR.2000.851948, 1998. 

Ambrósio, J. A. C.: Impact of Rigid and Flexible Multibody Systems: Deformation Description and Contact Models, in: Virtual Nonlinear Multibody Systems, edited by: Schiehlen, W. and Valášek, M., Springer, Dordrecht, 57–81, https://doi.org/10.1007/978-94-010-0203-5_4, 2003. 

Arnol, R. and Naderi, F. M.: Advanced architectures and the required technologies for next-generation communications satellite systems, Acta Astronaut., 20, 185–195, https://doi.org/10.1016/0094-5765(89)90068-4, 1989. 

Azadi, E., Fazelzadeh, A., and Ahmad, M.: Thermally induced vibrations of smart solar panel in a low-orbit satellite, Adv. Space Res., 59, 1502–1513, https://doi.org/10.1016/j.asr.2016.12.034, 2017.  

Fazelzadeh, S. A. and Azadi E.: Thermoelastic vibration and maneuver control of smart satellites, Aircr. Eng. Aerosp. Tec., 89, 477–490, https://doi.org/10.1108/AEAT-11-2015-0241, 2017. 

Gantes, C. J.: Deployable structures – analysis and design, Southampton: WIT Press, England, 2001. 

Hu, H. Y., Tian, Q., Zhang, W., Jin, D. P., Hu, G. K., and Song, Y. P.: Nonlinear dynamics and control of large deployable space structures composed of trusses and meshes, Adv. Mech., 43, 390–414, https://doi.org/10.6052/1000-0992-13-045, 2013 (in Chinese). 

Kim, Y. A.: Thermal creak induced dynamics of space structures, PhD thesis, Massachusetts Institute of Technology, Boston, 1999. 

Lankarani, H. M. and Nikravesh, P. E.: A contact force model with hysteresis damping for impact analysis of multibody systems, J. Mech. Des., 112, 369–376, https://doi.org/10.1115/1.2912617, 1990. 

Lu, G. Y., Zhou, J. Y., Cai, G. P., Fang, G. Q., Lv, L. L., and Peng, F. J.: Studies of thermal deformation and shape control of a space planar phased array antenna, Aerosp. Sci. Technol., 93, 105311, https://doi.org/10.1016/j.ast.2019.105311, 2019. 

Rahmat-Samii, Y. and Densmore, A. C.: Technology Trends and Challenges of Antennas for Satellite Communication Systems, IEEE T. Antenn. Propag., 63, 1191–1204, 2015. 

Shen, Z. X., Tian, Q., Liu, X. N., and Hu, G. K.: Thermally induced vibrations of flexible beams using Absolute Nodal Coordinate Formulation, Aerosp. Sci. Technol., 29, 386–393, https://doi.org/10.1016/j.ast.2013.04.009, 2013. 

Tang, Y. Q., Li, T. J., Wang, Z. W., and Deng, H. Q.: Surface accuracy analysis of large deployable antennas, Acta Astronaut., 104, 125–133, https://doi.org/10.1016/j.actaastro.2014.07.029, 2014. 

Thornton, E. A. and Kim, Y. A.: Thermal induced bending vibrations of a flexible rolled-up solar array, J. Spacecraft Rockets., 30, 438–448, https://doi.org/10.2514/3.25550, 1993.