Available online at www.sciencedirect.com ScienceDirect Procedia Engineering00 (2014) 000 000 www.elsevier.com/locate/procedia APISAT2014, 2014 Asia-Pacific International Symposium on Aerospace Technology, APISAT2014 Study on Impact Dynamics of Development for Solar Panel with Cylindrical Clearance Joint Wang Xu-peng a, b, *, Liu Geng a a School of Mechatronic Engineering, Northwestern Polytechnical University, Xi an Shanxi 710072, China b Aviation Equipment Research Institute, Qing An Group Co, Xi an Shanxi 710077, China Abstract In this study, an improved contact model in revolute joint with clearance is established by considering the axial dimension of the joint, energy dissipation and the nonlinear power exponent from material. The friction effect is described using the modified Coulomb friction model. Then simulation is carried out to investigate the influence of clearance size, input driving torque, and friction coefficient on impact dynamics of development for solar panel with cylindrical clearance joint. The results show that the effect of the clearance size, input driving torque and friction coefficient on angular displacement is slight, but obvious on angular velocity, and extraordinarily obvious on the joint reaction force, which may effects the accuracy and dynamic stability of the system, even induce the dynamic behavior of system tends to be nonlinear. 2014 The Authors. Published by Elsevier Ltd. Peer-review under responsibility of Chinese Society of Aeronautics and Astronautics (CSAA). Keywords: impact dynamic; satellite solar panel; cylindrical clearance joints; contact model; dynamic behavior Nomenclature α(deg) ω(deg/sec) F (Newton) angular displacement angular velocity contact force * Corresponding author. Tel.: +86 29 846 364 88; fax: +86 29 846 364 88. E-mail address:wangxupeng@avic.com 1877-7058 2014 The Authors. Published by Elsevier Ltd. Peer-review under responsibility of Chinese Society of Aeronautics and Astronautics (CSAA).
2 Wang Xupeng,Liu Geng/ Procedia Engineering00 (2014) 000 000 1. Introduction In general dynamic analysis of spacecraft system, such as dynamics of development for solar panel, it is assumed that the kinematic joints are ideal or perfect, that is, clearance, deformation, wear and lubrication effects are neglected. However, all of the mechanical systems do not have perfect kinematic joints due to the manufacture and assemblage errors, and the influence of wear [1-5]. Such clearance is necessary to allow the relative motion between the connected bodies and to permit the components assemblage as well, but it also leads to degradation of the performance of mechanical systems. By dynamic simulation of spacecraft, Chen and Xu [6] studied the influence of deploying flexible appendages on the attitude of main body in ADAMS and Matlab/Simulink circumstance. Bai et al. [7] presented an approach for the effect analysis of the spacecraft attitude motion undergone the flexible solar panel development and locking. Wang et al. [8] analyzed several important problems such as deployment procedure of flexible appendages on complex spacecraft, and corresponding program and application software are introduced. However, these studies have less focus on the effects of clearances on the dynamics performance of the spacecraft system. Nurre G S et al. [9] and Olivier A et al. [10] indicated that contact force induced by the clearance in joints have a negligible influence on the dynamic characteristic of spacecraft. The paper focus on investigating the joint clearance effects on the dynamics characteristics of solar panel with cylindrical clearance. The contact dynamics model in the clearance is established based on the improved contact model, which considers the axial dimension of the clearance joint, energy dissipation and the nonlinear power exponent from material, and the friction effect also is considered using the modified Coulomb friction model. Then the dynamics model of small monoplane satellite system with cylindrical clearance is established and the large numbers of dynamics simulations are presented. Finally, the effects of clearance size, input driving torque and friction coefficient on the dynamics performance of development for solar panel are analyzed in detail. 2. Modeling revolute clearance joint 2.1. Modeling a joint clearance As described above, the existence of clearance in the joint of small monoplane satellite system is inevitable. The key to modeling revolute clearance joint lies in effectively describing relative motion between the joint elements, and making an accurate judgment whether contact appearance. Fig. 1(a) depicts a revolute joint with clearance, and the clearance c is defined as follows: c R R (1) i j Where R i and R j are the radii of the bearing and journal. Fig.1. (a) Revolute joint with clearance; (b) Types of journal motion inside the bearing boundaries; (c) Revolution joint clearance with contact.
Wang Xupeng,Liu Geng / Procedia Engineering 00 (2014) 000 000 3 There are three main modeling methods for revolute clearance joints, namely, the massless link approach, the spring-damping approach and the momentum exchange approach [11]. The third method is more widely used, as it truly depicts the contact effect by simulating the elasticity properties of contacting surfaces, and the dissipation of energy during the impact process, it also classifies the journal motion inside the bearing into three different models, the free flight mode, the contact or following mode, and the impact mode, which are illustrated in Fig. 1(b). According to Fig. 1(c), the XOY coordinate frame represents the global coordinate system. R io and R jo are the positional vectors of the bearing and journal. Point O i indicates the center of the bearing and the center of the journal is denoted by the point O j. The clearance vector is given by [12] e= o o ij Ri Rj (2) Where e ij is the eccentricity vector that defines the distance between journal and bearing center, which can be stated as e= e e 2 2 ij x y (3) The contact penetration depth can be shown by e ij c (4) When e ij is bigger than c, the contact comes into between bearing and journal. In order to estimate the dissipation of energy during the impact process, the velocity between bearing and journal must be calculated, the relative velocities in the normal and tangential directions can be computed as vn n (5) vt t (6) Where t represents the unit tangential vector, it can be obtained by reversing the unit normal vector n for 90. 2.2. Modeling of contact force It is crucial to evaluate the contact force efficiently and exactly, which is the basic for studying on the impact dynamics of development for solar panel with clearance joint. The expression of the contact force is generally as follows [13] F F F C N T (7) where F N and F T are normal and tangential force component, respectively. So far, based on Hertz s law, a variety of F N models have been proposed by Brändlein and his co-authors [14], Yang and Sun [15], Dubowsky and Freudenstein [16], and Goldsmith [17] et al, it is worth noting that all these F N models do not consider the energy dissipation during the contact process. In order to overcome the limitation above, Kelvin and Voigt [18], Hunt and Crossley [19], Lankarani and Nikravesh [20], Zhiying and Qishao [21], and Flores et al [11] researchers extended the Hertz s law to accommodate energy dissipation in the form of internal damping. In particular, the model
4 Wang Xupeng,Liu Geng/ Procedia Engineering00 (2014) 000 000 developed by Lankarani and Nikravesh is widely used in dynamic studies of mechanical multi-body systems with clearance joints, and this model is expressed as [20] n K D 0 FN FK FD 0 0 (8) Where F K represents the elastic force component, F D explains the energy dissipation. K and δ are the generalized stiffness parameter and the nonlinear power exponent from material and the geometric properties of contacting objects. D is the damping coefficient, and is the relative impact velocity. The stiffness coefficient K is given by 4 RR i j K 3 ( ) R-R i j i j (9) In which the material parameters σ i and σ j are given by 1 l E 2 l l,( l i, j) (10) And υ i and υ j are the Poisson s ratio and Young s modulus associated with each body, respectively. The damping coefficient is given as [17] D n 2 3 (1 ) K c r ( ) 4 (11) ( ) Where c r is the restitution coefficient and is the initial impact velocity. However, the contact model of Lankarani and Nikravesh is only effective in the case of large clearance size with a small normal load [22]. Liu et al. [23] extended the Hertz contact law and presented a compliant force model, which was compared and validated with the results obtained with FEM analysis, and it also can be used to a wider range of impact conditions to meet different clearance and load. In this contact force model, the relationship between load and rigid body displacement can be expressed as F * E 1/2 ( ) 0 ' K 2 2( c ) 0 0 (12) Where δ is the relative indentation, c is the radial clearance size and the E * is the composite modulus of the colliding bodies, which can be expressed as 2 2 * 1l 1l 1 E ( ),( l i, j) E E l l (13) But the contact model established by Liu et al. does not consider the nonlinear power exponent from material, and energy dissipation during the contact process. In order to overcome these limitations, based on the contact
Wang Xupeng,Liu Geng / Procedia Engineering 00 (2014) 000 000 5 model of Liu, the paper introduces an improved contact model, which considers energy dissipation, the nonlinear power exponent from material and the axial dimension of the joint, the improved contact model can be evaluated using the following expression * EL 1 FNimp FKimp FDimp 2 2( c ) 0 0 n 1/2 ( ) Dimp 0 (14) Similar to the contact model of Lankarani and Nikravesh, the first term represents the elastic force and the second term explains the energy dissipation. D imp is the improved damping coefficient and can be written as [24] 0 0 Dimp Step(,0,0,, D ) D ( ) (3 2 ) 0 D 2 max max max max max max max max (15) Usually, the value of the δ max is 0.1mm, and the D max is the max damping coefficient relative to improved stiffness K N. The D max can be expressed as D K max N The value of the η is %[24], and the improved stiffness K N can be defined as (16) K N * n ( Fkimp ) E L n 1 1/2 1 1 1/2 ( ( ) ( ) ) 2 2( c ) 4 2( c ) (17) Fig.2. (a) Coefficient of friction varying with slips velocity; (b) Sketch map of small monoplane satellite system. The tangential friction F T, known as the Coulomb, is also crucial to depict the contact force. Classical Coulomb model does not refer the influence of relative slip velocity on tangential friction, so it does not accurately describe the conversion of friction between different states [25]. In recent years, Dubowsky [16], Rooney and Deravi [26], Threlfall [27], and Bai [28] et al. researchers modified the Coulomb model, established various new friction models. The paper used the Bai s friction model, which is written as F =- ( v ) F T T Nimp (18)
6 Wang Xupeng,Liu Geng/ Procedia Engineering00 (2014) 000 000 Where μ (v T ) is the friction coefficient relative to relative slip velocity v T that can be written as dsign( v ) vt > vd ( v) step( v, v,, v, ) sign( v) v v v step( v, v,, v, ) v < v T d d s s s T d T s s s s T s (19) Where, μ d and μ s are the dynamic friction coefficient and the max static friction coefficient, v d and v s are friction transition velocity and the stiction transition velocity. The Fig. 2(a) below shows the coefficient of friction varies with slip velocity [28] 3. Numerical example of solar panel with cylindrical clearance joints Fig. 2(b) shows a typical small satellite system with monoplane which consists of satellite body and solar panel. In the model, solar panel and satellite body are connected by a revolute cylindrical joint with clearance. In this work, the satellite body and solar panel are considered as rigid body. The solar panel is driven by input torque, which is supplied by torsion spring. The input torque is expressed as T= T 0 k (20) Where T 0 is the initial torque and k is the stiffness coefficient of the torsion spring, θ is the angle displacement of the solar panel. The initial coefficient of the torsion spring as follow, T 0 is 0.135N.m, and k is 0.0015N.m/deg. The mass and inertia properties of the small monoplane satellite system are shown as Table 1. Table 2 presents the parameter used in the dynamic simulation. Table 3 presents the friction coefficient of lubricated and non-lubricated. Table 1 Mass and inertia properties of small monoplane satellite system. Description Satellite body Solar panel Mass (kg) 768.0 7.68 J X (kg/m 2 ) 104.96 1.6137 J Y (kg/m 2 ) 81.92 0.4156 J Z (kg/m 2 ) 104.96 2.0293 J XY (kg/m 2 ) 0 0 J YZ (kg/m 2 ) 0 0 J ZX (kg/m 2 ) 0 0 Table 2 Parameters used in the dynamic simulations. Parameter Value Coefficient η 1 Yong s modulus 207GPa Poisson s ratio 0.3 Diameter of bearing Width of bearing 20mm 16mm Integration step size 0.0001s Integration tolerance 0.00001s
Wang Xupeng,Liu Geng / Procedia Engineering 00 (2014) 000 000 7 Table 3. Friction coefficient under lubricated and no lubricated. Description lubricated No lubricated Vs (mm/sec) 0.1 0.1 Vd (mm/sec) 10 10 ms 0.08 0.3 md 0.05 0.25 4. Numerical example of solar panel with cylindrical clearance joints This section contains extensive results obtained from computational simulations, which are based on the main functional parameters of the small monoplane satellite system, chiefly, clearance size, input driving torque, and friction coefficient. 4.1. Influence of the clearance size The influence of the clearance size on the dynamic behavior of the solar panel is investigated in this section. Fig.3. Simulation results under different clearance: (a) Angular displacement of solar panel; (b) Angular velocity of solar panel; (c) Angular displacement of satellite; (d) Angular velocity of satellite
8 Wang Xupeng,Liu Geng/ Procedia Engineering00 (2014) 000 000 In Figs 3, 4, and 5, the solar panel angle displacement, angle velocity, and the satellite angle displacement, angle velocity, as well as the joint reaction force and the journal center trajectories are used to illustrate the dynamic behaviour of the solar panel when different clearance sizes are considered. The values for the clearance sizes are chosen to be 0.10, 0.25, and 0.50mm. Figure 3(a), (c) shows that the angle displacement curves of solar panel and satellite are almost coincident to the ideal curves when clearance sizes are chosen to be 0.10 and 0.25mm, and the error between actual curve and ideal curve becomes larger when clearance size is 0.5mm. It indicates that the clearance has slight effect on angle displacement. Figure 3(b), (d) shows that the effect of clearance on angle velocity is obvious. The angle velocity of solar panel and satellite are shaky around the ideal value, and the higher clearance size, the more obvious shake, which will have effect on the movement accuracy and dynamic stability. Figure 4 shows the joint contact force is strongly affected by the clearance. The contact force in clearance joint is obviously shaken and the amplitude is increased from the ideal joint, contact force becomes to impulse type. And the larger clearance size, the contact force convergence slower, resulting in more time spent in the panel to reach steady movement state. Fig.4. The joint reaction force under different clearance: (a) Ideal joint; (b) clearance joint Figure 5 shows that when the clearance size is increased, there are higher numbers of impacts take place between journal and bearing, which indicates the dynamic behaviour of system tends to be nonlinear when the clearance size is increased. Fig.5. Journal center trajectories under different clearance size: (a) c=0.25; (b) c=0.5
4.2. Influence of input driving torque Wang Xupeng,Liu Geng / Procedia Engineering 00 (2014) 000 000 9 The influence of the driving torque on the dynamic behaviour of the solar panel is investigated in this section. Define the clearance size is 0.5mm, the angle displacement of solar panel is 56.5, and the friction coefficient is consistent with 4.1. Fig.6. Simulation results under different driving torque: (a) Angular displacement of solar panel, (b) Angular velocity of solar panel, (c) Angular displacement of solar panel, (d) Angular velocity of solar panel Fig.7. The joint reaction force under different driving torque: (a) Ideal joint; (b) clearance joint
10 Wang Xupeng,Liu Geng/ Procedia Engineering00 (2014) 000 000 In order to study the effect of the driving torque on the dynamic behaviour of the solar panel, three different values for input driving torque were selected, T, 3T and 5T. The simulation results are shown in Figs 6, 7, and 8. By and large, the influence of the clearance joint on the output results is similar for three different driving torque simulations, that is, as the driving torque increases, the solar panel can quickly reach predetermined angle, but it also cause angle velocity curve and the joint reaction force curve shake significantly. In addition, for higher values for driving torque, the system s behaviour tends to be more obvious nonlinear or chaotic. This conclusion can be confirmed in the journal center trajectories maps displayed in Fig.8. Fig.8. Journal center trajectories under different driving torque: (a) 3 times driving torque; (b) 5 times driving torque 4.3..Inflence of the friction coefficient The influence of the friction coefficient on the dynamic behaviour of the solar panel is investigated in this section. Again, define the clearance size is 0.5mm, simulation time is 5s, and the input driving torque is T. In order to study the effect of the friction coefficient on the dynamic behavior of the solar panel, the simulation was performed with two types: non lubrication and lubrication. Figs 9(a), (c) shows that friction coefficient has less effect on angle displacement of the solar panel and satellite, and the simulation curve of lubricated and nonlubricated is consistent. Comparing to angle velocity without lubrication, the simulation curve of lubricated has obvious shake, and this is due to without lubrication, the friction damping is lager, so the energy of the system, impact force and impact numbers are inhibited. This observation can be drawn by examination of Figs. 9(b), (d), 10 and 11, respectively.
Wang Xupeng,Liu Geng / Procedia Engineering 00 (2014) 000 000 11 Fig.9. Simulation results under different friction coefficient: (a) Angular displacement of solar panel, (b) Angular velocity of solar panel, (c) Angular displacement of satellite, (d) Angular velocity of satellite Fig.10. The joint reaction force under different friction coefficient: (a) Ideal joint; (b) clearance joint
12 Wang Xupeng,Liu Geng/ Procedia Engineering00 (2014) 000 000 Concluding remarks Fig.11. Journal center trajectories under different driving torque: (a) 3 times driving torque; (b) 5 times driving torque In order to study on impact dynamics for development of solar panel with cylindrical clearance joint, a general methodology for modelling and simulating clearance joint in small monoplane satellite system was presented and discussed throughout this work. Considering the axial dimension of the joint, the improved contact dynamics model in the joint clearance was established and the friction force also was considered using the modified Coulomb friction model. The main functional parameters of the small monoplane satellite system were used during the dynamics simulation, namely, clearance size, input driving torque, and friction coefficient. From the numerical simulations performed, it can be conclude that the clearance has obvious effects on the dynamics characteristics. (1)The effect of the clearance on angle displacement is slight, but obvious on angle velocity, and extraordinarily obvious on the joint reaction force, the joint reaction force is obviously higher than that without clearance. The angle velocity curve and joint reaction force curve of solar panel are shaky around the ideal value, and the higher clearance size, the more obvious shake. (2)Under the same clearance size and friction coefficient, the higher input driving torque can make the solar panel more quickly reach predetermined angle, but it also causes angle velocity curve and the joint reaction force curve shake significantly. (3)Under the same clearance and input driving torque, the higher friction force can effectively inhibit the joint reaction force, which can make the system reach stabilization more quickly. (4)Higher clearance size and input driving torque, and lower friction coefficient may induce obvious shake, and larger shake amplitudes of the angular velocity and the joint reaction force, which may effects the accuracy and dynamic stability of the system, even induce the dynamic behaviour of system tends to be nonlinear. Finally, it should be highlighted that the result presented in this paper does not consider the effect of the elasticity of the solar panel, in fact, it may tend to reduce the value of the angular displacement, and angular velocity as well as the joint reaction force, even change the entire behaviour of the system. Acknowledgements The support of the National Science Foundation of China (51275423) and the Foundation of the Discipline innovative talent recruitment program (B13044) is gratefully acknowledged. References [1] P. Flores, Modeling and simulation of wear in revolute clearance joints in multibody systems, Mechanism and Machine Theory 37(2009)1211-1222. [2] J.C.G. Oeden, Analysis of clearance in multibody system, Multi-body System Dynamics 13(2005)401-420.
Wang Xupeng,Liu Geng / Procedia Engineering 00 (2014) 000 000 13 [3] I. Khemili, L. Romdhane. Dynamic Analysis of a Flexible Slider-Crank Mechanism with Clearance, European Journal of Mechanics A/Solids 27(2008) 882-898. [4] P. Flores, A parametric study on the dynamic response of planar multibody systems with multiple clearance joints, Nonlinear Dynamics 61(2010) 633-653. [5]P. Flores, H.M. Lankarani. Spatial rigid-multibody systems with lubricated spherical clearance joints: modeling and simulation, Nonlinear Dynamics 60(2010)99-114. [6] T. Chen, S.J. Xu, Dynamic Simulation of Spacecraft Appendages Deploying, Aerospace Control 23(1) (2005) 79-83. [7] Z.F. Bai, Y. Zhang, H. Tian, Dynamics Simulation of Deployment of Solar Panels in Fault Modes, Journal of System Simulation 19(13) (2007) 3067-3072. [8] X.R. Ma, T.S. Wang, Dynamics Analysis of Lager Spacecraft With Deploying Flexible Appendages, Chinese Space Science and Technology 8(4) (2000) 1-8. [9] G.S. Nurre, J.P. Sharkey, J.D. Nelson, Preserving mission, on-orbit modifications to Hubble space telescope pointing control system, Journal of Guidance Control and Dynamics 18(2) (1995) 222-229. [10] A. Olivier, R. Jesus, Modeling of joints with clearance in flexible multibody system, International Journal of Solids and Structures 39(2002) 41-63. [11] P. Flores, J. Ambrósio, Revolute joints with clearance in multibody systems, Computers and Structures 82 (2004) 1359-1369. [12] T. Qiang, Y. Zhang, L. Chen, P. Flores, Dynamics of spatial flexible multibody systems with clearance and lubricated spherical joints, Computers and Structures 87 (13-14) (2009)913-929. [13] Selcuk Erkaya, Investigation of joint clearance effects on welding robot manipulators, Robotics and Computer-Integrated Manufacturing 28 (2012) 449-457. [14] J. Brändlein, P. Eschamann, L. Hasbargen, Die Wälzlagerpraxis, Handbuch für die Berechnung and Gestaltung von Lagerungen, Vereinigte Fachverlage, Germany, 1998. [15] D.C.H. Yang, Z.S. Sun, A rotary model for spur gear dynamics, Journal of Mechanisms, Transmissions and Automation in Design107(1985)529-535. [16] S. Dubowsky, F. Freudenstein, Dynamic analysis of mechanical systems with clearances, Part 1: formulation of dynamic model, Journal of Engineering for industry 93(1971)305-309. [17] W. Goldsmith, Impact-The theory and Physical Behaviors of Colliding Solids, Edward Arnold Ltd., London, England, 1960. [18] G.Gilardi, I. Sharf, Literature survey of contact dynamics modeling, Mechanism and Machine Theory 37(2002)1213-1239. [19] K.H. Hunt, F.R.E. Crossley, Coefficient of restitution interpreted as damping in vibroimpact, Journal of Applied Mechanics 7(1975)440-445. [20] Q. Zhiying, L. Qishao, Analysis of impact process based on restitution coefficient, Journal of Dynamics and Control 4(2006)294-298. [21] H.M. Lankarani, P.E. Nikravesh, A contact force model with hysteresis damping for impact analysis of multibody systems, Journal of Mechanical Design112(1990)369-376. [22] LOVE AEH, A Treatise on the Mathematical Theory of Elasticity (Fourth Edition), New York: Dover Publications, 1944. [23] C.S. Liu, K. Zhang, Y. Lei, The compliance contact model of cylindrical joints with clearance, Acta Mech Sinica 21(2005)451-458. [24] MDI. Building models in ADAMS/View.2002. [25] P. Flores, Dynamic Analysis of Mechanical Systems with Imperfect Kinematic Joints. Portugal: Universidad Do Minho PhD thesis, 2004. [26] G.T. Rooney, P. Deravi, Coulomb Friction in Mechanism Sliding Joints. Mechanism and Machine Theory 17(1982) 207-211. [27] D.C. Threlfall, The Inclusion of Coulomb Friction in Mechanisms Programs with Particular Reference to Dram. Mechanism and Machine Theory13 (1978)475-483. [28] Z.F. Bai, Y. Zhang, H. Tian, Dynamics characteristics of mechanisms with joint clearance, Journal of Vibration and Shock 30(11) (2011) 17-20.