Research Article  Open Access
Dynamics Analysis of Deployable Structures considering a TwoDimensional Coupled ThermoStructural Effect
Abstract
The deployment accuracy of deployable structures is affected by temperature and flexibility. To obtain the higher accuracy, various measures such as the optimization design and the control process are employed, and they are all based on deployment dynamics characteristics of deployable structures. So a precise coupled thermostructural deployment dynamics analysis is important and necessary. However, until now, only a onedimensional thermal effect is considered in the literatures because of simplicity, which reduces the accuracy of the model. Therefore, in this paper, a new model coupling mechanical field with a temperature field is presented to analyze the deployment dynamics of a deployable structure with scissorlike elements (SLEs). The model is based on the absolute nodal coordinate formulation (ANCF) and is established via a new lockingfree beam element whose formulation is extended to account for the twodimensional thermally induced stresses due to the heat expansion for the first time. Namely, in the formulation, the thermal influences are along twodimensional directions, the axial direction and the transverse direction, rather than along a onedimensional direction. The validity and precision of the proposed model are verified using a flexible pendulum example. Finally, the dynamics of a linear deployable structure with three SLEs modeled by the element is simulated under a temperature effect.
1. Introduction
Over the past three decades, deployable structures have been extensively used in space missions because they own the properties of high stiffness, low mass, and small folding volume [1–4]. And with the rapid development of aerospace industry, higher demands for these structures have been proposed. One of the most important is higher deployment accuracy, which is affected by thermal and flexible effects [5].
In order to achieve higher accuracy, the precise deployment dynamics characteristic of deployable structures is essential because it is fundamental to the optimization design [6] and the control process [7]. However, due to the high cost of the imitated space environment on the ground, it is difficult to predict the dynamic capability of these structures by experiment. So some practical methods are developed for dynamic analysis of these flexible structures. The most classical approach is the floating frame of reference formulation [8], but it is deficient in dynamic stiffening and will lead to an imprecise simulation result. To overcome this drawback, Shabana presented a new approach called the absolute nodal coordinate formulation (ANCF) [9], which can satisfy the dynamic stiffening automatically [10, 11] and leads to a constant mass matrix and zero centrifugal and Coriolis forces [12]. Because of these advantages, several beam elements [11, 13–16], plate elements [17], and brick elements [18] applied to analyze various objects of study have been developed within this formulation. A twodimensional shear beam element proposed by Omar and Shabana [11] to consider the shear deformation effect is among them. However, in a later research, it was found out that this beam element deriving elastic forces on the basis of the nonlinear elastic theory suffers from volumetric Poisson’s locking, thickness locking, and shear locking, which influence the accuracy and computational efficiency of ANCF greatly [19, 20]. In order to avoid the Poisson’s locking, an exact description of elastic forces is presented by Sopanen and Mikkola [21], but the simplest way is setting the Poisson ratio to zero directly. Possible alternatives to solve other locking problems are selective reduced integration [22–24] and the redefinition of elastic forces based on the HellingerReissner principle [20, 25].
As the absolute nodal coordinate formulation grows and is verified in many fields [26, 27], some researchers start to apply it in space deployable structures. Li et al. [28] use ANCF to simulate the deployment dynamics of the deployable antenna reflector. Tian et al. [29] simulate the deployment dynamics of two types of deployable structures based on ANCF and find a good agreement between simulations and experiments. In these simulations, all of researchers do not take into account the temperature effect and only investigate the deployable systems under the normal temperature. Li and Wang [30] fill the gap by discussing the deployable driving torque at the lowest and highest temperature within ANCF, but their model is based on a simple Euler–Bernoulli beam element and only the axial thermal effect is pondered. As mentioned above, however, a precise dynamics analysis is important and necessary for deployable structures. Besides, several highprofile examples have verified that the temperature effect has a great influence on deployable structures [31, 32]. So the multidimensional coupled thermostructural effect on the deployment process of deployable structures must be considered.
To this end, a new model coupling mechanical field with a temperature field of a deployable structure with scissorlike elements (SLEs) based on ANCF is presented. The model is established by a new planar lockingfree shear deformable beam element inspired from the gradientdeficient element which alleviates the thickness locking. In order to simulate the deployment dynamics characteristic accurately, the formulation of the shear deformable beam element is extended to contain the thermal stress due to heat expansions. Note that the temperature field is uniform and the heat expansions are along twodimensional directions; that is to say, the thermal loads are applied in an axial direction and a transverse direction, which is a great improvement compared with that of Li and Wang [30]. Also, the heat expansion coefficients are equal in twodimensional directions because of the isotropic material in the model. This formulation is verified by ANSYS with a simple pendulum driven by gravity. In the end, on the basis of the new model, the simulation results of deployment dynamics of the deployable structure considering the coupled thermostructural effect are calculated out and compared with results of only flexible analysis of the same structure to discuss the temperature effects.
This paper is organized as follows. The coupled thermostructural formulation of the lockingfree shear deformable beam element with twodimensional thermal loads is discussed in Section 2. In Section 3, the formulation is demonstrated and a model of the deployable structure consisting of SLEs is established with the formulation. Also, the coupled thermostructural analysis of other deployable structures comprising rods can be created by using this formulation. The simulation results of the deployable structure with SLEs are illustrated in Section 4. Finally, some valuable conclusions are drawn in Section 5.
2. A LockingFree Shear Deformable Beam Element under the Thermo Effect
In this paper, using a twodimensional shear deformable beam element [11] to model the rods is appropriate because motions and forces of the deployable structure with SLEs are in plane and no torsion occurs. However, the beam element suffers from the locking phenomena. So a new lockingfree shear deformable beam element is presented in Figure 1. In the improved beam kinematics, an additional slope vector is added in the middle of the gradientdeficient element [19].
An arbitrary point on the beam element can be expressed as where is the additional shape function compared with [19]
The matrix is a identity matrix, the elements are given as follows, and their detailed derivation process is shown in the appendix. where , , is the length of the beam element, and are the element coordinates in the element local coordinate system. So the nodal coordinate vector can be written as where , , and are the global positions of nodes and , , , and are the global slopes of the element nodes.
The velocity at any material point can be obtained:
Thus, the kinetic energy of the beam element is expressed as follows: where is the mass matrix and is constant. where is the volume of the beam element and is the material density of the beam.
Before deducing the strain energy, in order to bring the thermal effect into the formulation of the beam element, some assumptions are introduced. The material is isotropous, and the temperature field is homogeneous along the beam element. So the strain vector due to the change of temperature can be written as [33] where is the coefficient of thermal expansion and is the change in the temperature field. According to continuum mechanics theory, the strain vector due to the deformation can be derived as follows. First, the deformation gradient is defined as where is the th row of the element shape function and is the absolute nodal coordinate vector of the beam element in the initial configuration. Hence, the Lagrangian strain tensor due to flexibility can be written as where
The strain tensor can be written as a vector form because of its symmetry: where and are normal strains along the direction, respectively, and is the shear strain.
So the strain vector considering the coupled thermostructural effect can be obtained:
In terms of constitutive relation, the stress vector related to the strain vector is given by where indicates the coefficient matrix of material elasticity, and it can be written as follows because of the isotropic homogenous material: where where and indicate the elastic modulus and Poisson’s ratio of the beam material, respectively, and indicates the shear modulus. So the strain energy is
The elastic forces can be deduced by using the strain energy . where is the stiffness matrix. where
Finally, according to Lagrange’s equation [9], the equation of motion for the beam element affected by temperature is given by where is the generalized external forces applied in the element.
3. Validation
A simple example that a rectangular pendulum moves from the horizontal rest position under the effect of gravity (), as shown in Figure 2, is examined to validate the proposed model with a twodimensional thermal load based on ANCF. And the thermal load is applied by changing the temperature from 0°C to 200°C suddenly. And other simulation parameters such as the material property and geometry size of the pendulum are given in Table 1.

Figure 3 shows the curves of displacement of point a which are evaluated by ANSYS, the 1D thermostructural model [30], and the proposed model, respectively. The two models are in good agreement with ANSYS within a large overall motion. It demonstrates that the proposed beam element is free of locking phenomena. Further, in order to observe the difference among three results clearly, a margin calculation is carried out. The difference between ANSYS and the flexible model, the difference between the 1D thermostructural model and the flexible model, and the difference between the proposed model and the flexible model are compared in Figure 4. Here, the flexible model implies no temperature effect. It is clear from the difference value present in Figure 4 that the proposed model is in good agreement with ANSYS, whereas the 1D thermostructural model has a huge gap with ANSYS. So the proposed model has higher precision. Therefore, the development of the new twodimensional thermostructural model is significant.
So the proposed model is used to assemble the deployable structure with SLEs under the thermal effect. The dynamics equations with the constraint conditions can be expressed as follows: where is the generalized coordinate vector of the deployable structure, and are the mass matrix and the stiffness matrix, respectively, is the generalized force vector, is the Lagrange multipliers, is the constraint equation, and is the derivative matrix of with respect to . Note that the stiffness matrix is the function of the temperature besides the generalized coordinate vector .
4. Numerical Results
In this section, a deployable structure with three scissorlike elements considering the coupled thermostructural effect is discussed. As shown in Figure 5, a scissorlike element is assembled by two rods which can be treated as beams with a pin, and the three scissorlike elements are connected to a deployable system with pins. A driving force equal to 50 N is applied at point a_{1}. The material property and geometry size of the structure are given in Table 2. And the motion simulation is carried out within the time interval , where the system is uniformly revelational to the temperature profile presented in Figure 6 with a maximum temperature of 150°C and a minimum temperature of −150°C [30]. In order to demonstrate the temperature influence on the dynamics of the mechanism, the simulation results are compared with the same system only considering the flexible effect.

Figure 7 shows the deflection curves of rod b_{1}a_{2}, rod b_{2}a_{3}, and rod b_{3}a_{4} in the changing temperature field as shown in Figure 6 and in the constant temperature field. Referring to Figure 7, it can be seen that the timedependent temperature leads to many small vibrations which can be observed 2.7 s after clearly when the two curves differ greatly at the moment. These small vibrations which will influence the deployment accuracy greatly arise because thermally induced transverse strain results in varying moment of inertia related to the crosssectional area besides the change of the length of rods resulting from thermally induced axial stress. Another notable point from the pictures is that the temperature influence on the deflection of the rods seems to weaken in this time period of 2 s to 2.5 s. A probable reason is that the lowest temperature effect on the deployable structure partly offsets the highest temperature effect in the earlier period.
(a) Rod b1a2 deflection with time
(b) Rod b2a3 deflection with time
(c) Rod b3a4 deflection with time
The motion diagrams of a representative point a_{4} selected to stand for the dynamics characteristics of the deployable structure are shown in Figure 8. Referring to Figure 8, the temperature impacting on deployment displacement and velocity is not obvious, but there is a huge influence on deployment acceleration of point a_{4}. The thermal effect results in additional vibrations emerging in acceleration, which can be observed by comparing the thermoflexible curve with the flexible curve. Besides, referring to Figures 6 and 8(c), it is found that acceleration vibrations become larger perspicuously at inflection points of the temperature profile such as 0.25 s, 1.25 s, 1.75 s, and 2.75 s. These vibrations may lead to unexpected collision and impact which might result in the invalidation of the deployable structure or even mission failure.
(a) Displacement with time
(b) Velocity with time
(c) Acceleration with time
5. Conclusions
In this paper, the dynamics analysis of the deployable structure with three SLEs considering the coupled thermostructural effect based on ANCF is presented. (1)In order to establish the precise model coupling of the mechanical field with the temperature field of the mechanism, the formulation of a new lockingfree beam element with the additional slope vector is extended to account for the thermally induced stresses in twodimensional directions, that is, the axial direction and transverse direction(2)It shows that the new beam element is free of locking phenomena as a result of a good agreement between the calculating results of the simple pendulum and the simulation results of the ANSYS within a large overall motion. Then, the validity and the precision of the proposed thermostructural model are demonstrated through the further margin calculation: the difference of the proposed model is approximately identical to the difference of ANSYS, whereas the difference of the 1D thermostructural model has a relatively large gap with the difference of ANSYS. Therefore, the proposed model is more precise than the 1D thermalstructural model. On the other hand, it also demonstrates that the development of the new twodimensional thermostructural model is significant(3)In the final simulations of the deployable structure, by contrasting the thermoflexible deflection with the flexible deflection, it is found that the main influence of temperature is that the change of temperature leads to many vibrations: for the deflection of rods, many small vibrations arise because thermally induced transverse strain results in varying moment of inertia related to the crosssectional area besides the change of the length of rods resulting from thermally induced axial stress. On the other hand, for the deployment dynamics characteristics, temperature has a huge impact on deployment acceleration. The thermal effect will result in additional vibrations emerging in acceleration of an arbitrary point on the deployable structure, and these vibrations will become larger perspicuously at inflection points of temperature change
Appendix
For the new beam element, the displacement field is defined as where are the polynomial coefficients.
The nodal displacement conditions shown in Figure 1 can be expressed as
Thus, the polynomial coefficients can be derived from (A.2):
By substituting (A.3) into (A.1), the displacement field is rewritten with nodal displacement: where
Data Availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Conflicts of Interest
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (grant number 51175422).
References
 Z. You and S. Pellegrino, “Deployable mesh reflector,” The International Symposium on Spatial, Lattice, and Tension Structures, pp. 103–112, 1994. View at: Google Scholar
 A. G. Cherniavsky, V. I. Gulyayev, V. V. Gaidaichuk, and A. I. Fedoseev, “Large deployable space antennas based on usage of polygonal pantograph,” Journal of Aerospace Engineering, vol. 18, no. 3, pp. 139–145, 2005. View at: Publisher Site  Google Scholar
 Y. Wang, R. Liu, H. Yang, Q. Cong, and H. Guo, “Design and deployment analysis of modular deployable structure for large antennas,” Journal of Spacecraft and Rockets, vol. 52, no. 4, pp. 1101–1111, 2015. View at: Publisher Site  Google Scholar
 B. Li, S. M. Wang, R. Yuan, X. Z. Xue, and C. J. Zhi, “Dynamic characteristics of planar linear array deployable structure based on scissorlike element with joint clearance using a new mixed contact force model,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 230, no. 18, pp. 3161–3174, 2016. View at: Publisher Site  Google Scholar
 Y. Tang, T. Li, Z. Wang, and H. Deng, “Surface accuracy analysis of large deployable antennas,” Acta Astronautica, vol. 104, no. 1, pp. 125–133, 2014. View at: Publisher Site  Google Scholar
 J. Sun, Q. Tian, and H. Hu, “Structural optimization of flexible components in a flexible multibody system modeled via ANCF,” Mechanism and Machine Theory, vol. 104, pp. 59–80, 2016. View at: Publisher Site  Google Scholar
 Y. Zhang, B. Duan, and T. Li, “A controlled deployment method for flexible deployable space antennas,” Acta Astronautica, vol. 81, no. 1, pp. 19–29, 2012. View at: Publisher Site  Google Scholar
 J. García de Jalón and E. Bayo, Kinematic and Dynamic Simulation of Multibody Systems, Mechanical Engineering Series, Springer, New York, 1994. View at: Publisher Site
 A. A. Shabana, Dynamics of Multibody Systems, Cambridge University Press, 2005. View at: Publisher Site
 M. Berzeri and A. A. Shabana, “Study of the centrifugal stiffening effect using the finite element absolute nodal coordinate formulation,” Multibody System Dynamics, vol. 7, no. 4, pp. 357–387, 2002. View at: Publisher Site  Google Scholar
 M. A. Omar and A. A. Shabana, “A twodimensional shear deformable beam for large rotation and deformation problems,” Journal of Sound and Vibration, vol. 243, no. 3, pp. 565–576, 2001. View at: Publisher Site  Google Scholar
 A. A. Shabana, Computational Continuum Mechanics, Cambridge University Press, 2012.
 M. Berzeri and A. A. Shabana, “Development of simple models for the elastic forces in the absolute nodal coordinate formulation,” Journal of Sound and Vibration, vol. 235, no. 4, pp. 539–565, 2000. View at: Publisher Site  Google Scholar
 A. A. Shabana and R. Y. Yakoub, “Three dimensional absolute nodal coordinate formulation for beam elements: theory,” Journal of Mechanical Design, vol. 123, no. 4, pp. 606–613, 2001. View at: Publisher Site  Google Scholar
 R. Y. Yakoub and A. A. Shabana, “Three dimensional absolute nodal coordinate formulation for beam elements: implementation and applications,” Journal of Mechanical Design, vol. 123, no. 4, pp. 614–621, 2001. View at: Publisher Site  Google Scholar
 S. von Dombrowski, “Analysis of large flexible body deformation in multibody systems using absolute coordinates,” Multibody System Dynamics, vol. 8, no. 4, pp. 409–432, 2002. View at: Publisher Site  Google Scholar
 O. N. Dmitrochenko and D. Y. Pogorelov, “Generalization of plate finite elements for absolute nodal coordinate formulation,” Multibody System Dynamics, vol. 10, no. 1, pp. 17–43, 2003. View at: Publisher Site  Google Scholar
 L. Kübler, P. Eberhard, and J. Geisler, “Flexible multibody systems with large deformations using absolute nodal coordinates for isoparametric solid brick elements,” in ASME 2003 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, pp. 31–52, Chicago, Illinois, USA, September 2003. View at: Publisher Site  Google Scholar
 D. GarcíaVallejo, A. M. Mikkola, and J. L. Escalona, “A new lockingfree shear deformable finite element based on absolute nodal coordinates,” Nonlinear Dynamics, vol. 50, no. 12, pp. 249–264, 2007. View at: Publisher Site  Google Scholar
 B. A. Hussein, H. Sugiyama, and A. A. Shabana, “Coupled deformation modes in the large deformation finiteelement analysis: problem definition,” Journal of Computational and Nonlinear Dynamics, vol. 2, no. 2, p. 146, 2007. View at: Publisher Site  Google Scholar
 J. T. Sopanen and A. M. Mikkola, “Description of elastic forces in absolute nodal coordinate formulation,” Nonlinear Dynamics, vol. 34, no. 1/2, pp. 53–74, 2003. View at: Publisher Site  Google Scholar
 K. S. Kerkkänen, J. T. Sopanen, and A. M. Mikkola, “A linear beam finite element based on the absolute nodal coordinate formulation,” Journal of Mechanical Design, vol. 127, no. 4, pp. 621–630, 2005. View at: Publisher Site  Google Scholar
 J. Gerstmayr, M. K. Matikainen, and A. M. Mikkola, “A geometrically exact beam element based on the absolute nodal coordinate formulation,” Multibody System Dynamics, vol. 20, no. 4, pp. 359–384, 2008. View at: Publisher Site  Google Scholar
 K. Nachbagauer, A. S. Pechstein, H. Irschik, and J. Gerstmayr, “A new lockingfree formulation for planar, shear deformable, linear and quadratic beam finite elements based on the absolute nodal coordinate formulation,” Multibody System Dynamics, vol. 26, no. 3, pp. 245–263, 2011. View at: Publisher Site  Google Scholar
 H. Sugiyama and Y. Suda, “A curved beam element in the analysis of flexible multibody systems using the absolute nodal coordinates,” Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multibody Dynamics, vol. 221, no. 2, pp. 219–231, 2007. View at: Publisher Site  Google Scholar
 W. S. Yoo, J. H. Lee, S. J. Park, J. H. Sohn, D. Pogorelov, and O. Dmitrochenko, “Large deflection analysis of a thin plate: computer simulations and experiments,” Multibody System Dynamics, vol. 11, no. 2, pp. 185–208, 2004. View at: Publisher Site  Google Scholar
 W. S. Yoo, M. S. Kim, S. H. Mun, and J. H. Sohn, “Large displacement of beam with base motion: flexible multibody simulations and experiments,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 5051, pp. 7036–7051, 2006. View at: Publisher Site  Google Scholar
 P. Li, C. Liu, Q. Tian, H. Hu, and Y. Song, “Dynamics of a deployable mesh reflector of satellite antenna: parallel computation and deployment simulation 1,” Journal of Computational and Nonlinear Dynamics, vol. 11, no. 6, article 061005, 2016. View at: Publisher Site  Google Scholar
 Q. Tian, J. Zhao, C. Liu, C. Zhou, and H. Hu, “Dynamics of space deployable structures,” in ASME 2015 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, ASME, Boston, USA, 2015. View at: Google Scholar
 T. Li and Y. Wang, “Deployment dynamic analysis of deployable antennas considering thermal effect,” Aerospace Science and Technology, vol. 13, no. 45, pp. 210–215, 2009. View at: Publisher Site  Google Scholar
 E. Thornton and R. Foster, “Dynamic response of rapidly heated space structures,” in Computational Nonlinear Mechanics in Aerospace Engineering, pp. 451–477, 33rd Structures, Structural Dynamics and Materials Conference, Dallas, TX, USA, 2013. View at: Google Scholar
 C. L. Foster, M. L. Tinker, G. S. Nurre, and W. A. Till, “Solararrayinduced disturbance of the Hubble Space Telescope pointing system,” Journal of Spacecraft and Rockets, vol. 32, no. 4, pp. 634–644, 1995. View at: Publisher Site  Google Scholar
 B. A. Boley, J. H. Weiner, and E. H. Dill, “Theory of thermal stresses,” Physics Today, vol. 14, no. 3, pp. 58–60, 1961. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Qi’an Peng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.