Modelling and Simulation in Engineering

Modelling and Simulation in Engineering / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 497916 |

V. S. Khoroshilov, A. E. Zakrzhevskii, "Modelling of Spacecraft Dynamics at Deployment of Large Elastic Structure", Modelling and Simulation in Engineering, vol. 2012, Article ID 497916, 13 pages, 2012.

Modelling of Spacecraft Dynamics at Deployment of Large Elastic Structure

Academic Editor: Zeki Ayag
Received20 Apr 2012
Revised29 Jun 2012
Accepted02 Jul 2012
Published22 Aug 2012


In this paper, a new approach to the modelling of the deployment dynamics of a flexible multi-body system with the time dependent configurations is demonstrated in the frame of the study the dynamics of a spacecraft with the gyro-gravitational system of stabilization. Primarily the gravitational stabilizer that is made as a pantograph structure is in a compact form. The deployment of a flexible pantograph structure is performed after placing the spacecraft into orbit and completion of the preliminary damping by a special jet-propelled system, and after uncaging the gyros. After its deployment, the pantograph turns into an elongated structure that serves as a gravitational stabilizer and carrier of solar batteries. The objective of the study is the creation of the generalized mathematical model and the conducting of the computational modelling of the spacecraft dynamics. The equations of motion are derived with the use of the Eulerian-LaGrangian formalism and symbolic computing. Numerical simulations of the typical operational mode of the system are conducted taking into account various control profiles for the deployment. Numerical results indicate that the system used for attitude stabilization ensures the shape of the deployed design and prescribed accuracy of the orientation.

1. Introduction

Structures created from elements delivered into orbit in a compact condition are one of the basic components of modern space systems. Spacecraft (SC) with long flexible appendages for the exploration of the Earth’s physical fields and gravitationally stabilized SC that deploys the long boom with the mass concentrated at the tip have been used since the end of the fifties of the previous century. The deployment of the space structures brings in considerable disturbances into the dynamics of the SC about its mass centre. At present, the study of such deployable structures is actually for minimization of deployment duration and power resources, for analysis of the effect of attached systems of a changeable configuration on the SC attitude. The problem of weight reduction of such systems is constantly actual. Therefore, elements of such systems should be considered as elastic bodies.

There is a large number of studies in the literature dedicated to the deployment of elastic appendages from the fixed basis as well as from the rotating SC [16], including the deployment of the gravity-gradient boom [7, 8]. The maximum bending moments and the deflection of the booms at different simplifying hypothesis were evaluated in these papers; the optimization of the boom deployment timing was also fulfilled. In a relevant paper, Cherchas [9] studied the dynamics of the satellites stabilized by the rotation in the process of extension of the long flexible booms. He determined the maximum angles of nutation precession after deployment and the maximum bending moments and deflections of the rod. The equations of motion were derived by the method of Lagrange, and the elastic degrees of freedom were discredited based on the modal analysis. Lips and Modi [10], Tsuchiya [11], and Lips et al. [12] also described the dynamics of the deployed elastic elements in terms of the modal analysis with the mode shapes of natural vibrations dependent on time. In a different approach to the problem under consideration, Banerjee and Kane [6] presented a unique method for the description of beam dynamics, as it is being deployed from or retrieved into a rigid rotating body. The essence of the method was the modelling of the beam as a chain of elastically connected rigid links, then the method worked with the equations of motion linearized in the modal coordinates for the links located at the considered time out of the main body. Banerjee [13] and Banerjee and Nagarajan [14] extended the analysis to cases of large deflections with the order-n formulation for a variable number of bodies. The simulation results showed that the deflections of the boom tip were very sensitive to the rate of deployment/retrieval, the forcing into turned out to be less persistent process out of two ones. It was established that the character of the angular velocity of the rigid body from/into which the beam is being deployed/retrieved influences the beam behaviour.

The special position among the large structures to be deployed in space belongs to tethered spacecraft, which is connected by long thin cables, and lightweight cable nets deployable in space to serve as platforms for very large structures. One of the most important operations during a tethered satellite system mission is the deployment and retrieval of a subsatellite from or to a spacecraft. There are many studies devoted to analytical and numerical analysis of the equations of tethers motion, focused on estimation or creation of control laws of their deployment. Levin [15] and Williams [16] give the detailed analysis of basic aspects of this problem.

Only a few of the studies, which deal with the dynamics of the SC of variable configuration, take into account the dynamics of the deployment actuating mechanism. Zakrzhevskii and Khoroshilov [17] studied the dynamics of the orbital deployment of an elastic ring-shaped antenna taking into account dynamics of actuating drive. Creamer [7] studied the transverse motion of the truss boom deployed along an axis, in order to compare the theory and experiment. The boom was deployed from a canister in the horizontal plane and was supported by a bungee cord tied with the slider at the ceiling of the building in order to compensate the gravitation. Dranovskii et al. [18] stop at the development of general mathematical model of the SC during the deployment of the flexible gravitation stabilizer according to the specified programme.

Dranovskii et al. [19] studied the dynamics of SC with regard to an elastic pantograph structure deployment, but transversal displacements of the pantograph were not taken into account here. Works, except for work [19], which would investigate the dynamics of deployment in orbit of pantograph designs, are unknown to the authors. This research generalises work [19] by the additional account of transversal displacements of the deployed pantograph design with objective to study their effect on the dynamics of the SC with gyro-gravy stabilization and its elements.

2. Problem Formulation

2.1. Physical Model of System

Here, the SC that includes gyro-gravitational system of stabilization is studied during the deployment of a flexible pantograph structure according to a specified program into the elongated flexible gravitational stabilizer (GS) that serves also as a carrier of solar batteries. SC includes two gyro-dampers (GD) which are installed in order to reduce amplitudes of angular oscillations.

The SC is placed in a circular orbit about the Earth of radius 6600 km. Basic elements of the SC under consideration are shown in Figure 1. Here, element 1 is the SC main module, 2 is the spatial structure that consist of two coupled plain pantographs, and 3 and 4 are the GD.

Figure 1 shows symbolically only three tiers of the design during deployment. Here, the structures up to 20 tiers were studied. The deployment starts when the points of joining each a plain pantograph to the main module begin to approach synchronously towards each other. The actuating drive dynamics are not taken into consideration since it is very simple and may be realized as a device that has no effect on the SC dynamics. The law of changing of the distance between the points of joining of the structure to the main module is prescribed in time.

Such deployable structures possess a considerable flexibility. Because of requirements to maintenance of the shape of the gravitational stabilizer, so that, for instance, to stabilize the spacecraft in the orbital reference frame with the predefined accuracy, some restrictions must be imposed on the deployment and design parameters. They can be determined in the process of numerical simulations. The GD role in the process of deployment and after its termination is also studied further.

2.2. Mechanical Model of System

For the case of deployment of the elastic pantograph structure in orbit, the generalized mechanical model may be represented as a main rigid body 𝑆1 and body 𝑆2 of variable configuration attached to it, as it is shown in Figure 2. Here, 𝐶𝑋𝑌𝑍 is an earth-centred inertial reference frame with axis 𝐶𝑌 along the earth’s spin axis, axis 𝐶𝑍 pointing to the vernal equinox, and axis 𝐶𝑋 completing the orthogonal triad. 𝑂𝑥𝑦𝑧 is a main body fixed reference frame with origin at the some point 𝑂 of the body 𝑆1 with axis 𝑂𝑧 along the design position of the GS axis. The body 𝑆1 is the gyro-static part; it includes the damping gyros, which do not change the rotational body inertia. The body 𝑆2 is the deployed part with 𝐶2 as the centre of mass and does not comprise any gyrostatic components. The body 𝑆1 motion is defined by the velocity vector 𝐯𝑂 of the point 𝑂 and vector of absolute angular velocity 𝝎. The position vector 𝐫 defines the location of the arbitrary point 𝑃 with respect to the inertial reference frame 𝐶𝑋𝑌𝑍, and the position vector 𝐫—with respect to the reference frame 𝑂𝑥𝑦𝑧, fixed to the main body. In contrast to the problem of dynamics of the relative motion of the attached bodies described by Lurie [20], here one has the more general case when the expression for 𝐫 depends on time 𝑡 explicitly, not only through the following generalized coordinates: 𝐫=𝐫𝑞1,,𝑞𝑛,𝑡,(1) since the deployment of the pantograph takes place in accordance with the prescribed in time law. As a result, 𝐫 changes in time during deployment even when the relative elastic motion of the structure is absent.

The additional frame of reference (Figure 3) can be useful also at the problem setting. The frame of reference 𝐶𝑥or𝑦or𝑧or is fixed in the mass centre of the SC. This frame is introduced in such a way as in [21] with 𝐶𝑥or along the SC velocity vector, 𝐶𝑦or along the orbit binormal, and 𝐶𝑧or along the local vertical.

In Figure 3, Ω is the longitude of the ascending node from a point of the vernal equinox, 𝑖 is the orbit inclination, 𝑢 is the argument of latitude, and 𝝎or is the vector of the orbital angular velocity.

Neglecting the orbital precession, one can consider the absolute angular velocity of the orbital frame of reference to be equal to 𝝎or.

For both pantographs, the elastic rods have their length 2𝑎 and mass 2𝑚𝑎(𝑖), where index 𝑖 is used for the numbering of tiers that form the full rhombuses. These rods are connected by the cylindrical hinges at the points 𝐸(𝑅)(𝑖),𝐸(𝐿)(𝑖),(𝑖=0,𝑁), where 𝑁 is the number of the tiers, and at points 𝐻(𝑅)(𝑖),𝐻(𝐿)(𝑖),𝐹(𝑅)(𝑖),𝐹(𝐿)(𝑖),(𝑖=1,𝑁), where the superscripts identify the pantograph according to Figure 4.

These two pantographs are denoted by 𝐿 and 𝑅 accordingly. Thus, each pantograph consists of the zero semirhombus and the subsequent 𝑁 rhombuses with the sides 𝑎. The pantographs are connected between themselves and also out of their plane by rigid rods of the length ̃𝑏 and mass 𝑚𝑏 at points 𝐻(𝑅)(𝑖),𝐻(𝐿)(𝑖) and 𝐹(𝑅)(𝑖),𝐹(𝐿)(𝑖), (𝑖=1,𝑁), on which solar batteries are mounted. Rods of length 𝐻(𝑅)(𝑖)𝐻(𝐿)(𝑖)=𝐹(𝑅)(𝑖)𝐹(𝐿)(𝑖)=̃𝑏, (𝑖=1,𝑁) may be considered as the third system, which is denoted further by the superscript 𝑀.

Change of the configuration of this structure occurs when the distance 2𝑏(𝑡) between points 𝐻0(𝑅),𝐹0(𝑅) and 𝐻0(𝐿),𝐹0(𝐿) of joining each a plain pantograph with the deployment actuating drive that is established on the main body begin to decrease under the given law. The actuating drive dynamics are not taken into consideration here since this device is very simple and may be realized as a device that has no effect on the SC dynamics.

At the initial instant, the pantograph structure is in a compact form (a transport condition). The inclination angles of the rods that combine in all tiers (including the zero tier) with respect to axis Ox are equal to 5. After the deployment termination, values of the specified angles reach 75. Here, 𝑂𝑥𝑦𝑧 is the frame of reference that is fixed in the main body with the origin in point 𝑂, which is located on the plane 𝐻(𝑅)(0)𝐻(𝐿)(0)𝐹(𝑅)(0)𝐹(𝐿)(0) at the intersection of diagonals of the given rectangle. Axes 𝑂(𝑅)𝑥1,𝑂𝑥,𝑂(𝐿)𝑥2 and 𝑂(𝑅)𝑧1,𝑂𝑧,𝑂(𝐿)𝑧2 are parallel.

Elastic rods of the physical model have been replaced in the mechanical model by equivalent constructions of two rigid rods connected by the spring-bias cylindrical hinges with damping. Damping is used in order to approach the dynamics to reality at least in a qualitative sense. Stiffness of the springs in the hinges is defined from the condition of equal deflections of either of the two constructions (Figure 5) for equal loads. The equivalent construction (in Figure 5 below) has the same deflection when spring resistance 𝑐joint in the hinge is equal to 3𝐸𝐽/(2𝑎). Such a replacement is completely justified, as the configuration of the design is defined by mutual positions of the middle and end points of each rod. Besides, let the centre of each hinge have certain position and velocity that do not contradict to connections in an arbitrary state of a pantograph design. Let this state be identical for a real design and for its model. One can define for such a state of both systems their tensor of inertia, moment of momentum, and kinetic and potential energies. One can see that these values are almost identical for the real design and for the model. Differences are defined only by the small items connected with difference of local dynamic characteristics of a real bent rod in length 2𝑎 and its model of two direct rods. These differences are negligible quantities as compared to these values for all design. The further complication of the model will not bring new effects. Introduction of the additional internal degrees of freedom concerned with elastic deviations in a direction of axis 𝑂𝑦 is not justified, since external forces are absent in this direction.

2.3. Mathematical Model of System

The equations of motion of the system under consideration become the most compact and convenient for numerical integration (both for angular and relative motions), if one chooses the instant position of the mass centre 𝐶 as an origin. In this case, kinetic energy can be presented in the following form: 𝑇𝑆1+𝑆2=12𝑚𝐯2𝑖1𝑑𝑚=2𝑚𝐯𝐶+𝝆+𝝎×𝝆2=1𝑑𝑚2𝑚𝐯2𝐶+𝝆2+(𝝎×𝝆)2+2𝐯𝐶𝝆+2𝐯𝐶(𝝎×𝝆)+2𝝆(𝝎×𝝆)𝑑𝑚.(2)

Here, 𝝆 is the position vector of the system arbitrary point with respect to the instant position of the centre of mass 𝐶.

After integration, the expression for the kinetic energy of the whole system is given by 𝑀𝑇=2𝐯2𝐶+𝑇𝑂𝑟𝑀2𝐫𝐶2+12𝝎×Θ𝐶𝝎+𝝎×𝐊𝐶𝑟.(3)

The following notations are used here: Θ𝐶 is the inertia tensor of the whole system with respect to point 𝐶; 𝐊𝐶𝑟=𝑚𝐫×𝐫𝑑𝑚𝑀𝐫𝐶×𝐫𝐶 is the relative moment of momentum of the deployed part with respect to point 𝐶; 𝐫𝐶 is position vector of the instant position of the centre of mass 𝐶 in the main body fixed frame of reference; the symbol * indicates the time differentiation in the main body fixed reference frame; 𝑀 is the total mass of the system; 𝐯𝐶 is the absolute velocity of the centre of mass; 𝑇𝑂𝑟 is the kinetic energy of the relative motion of the carried bodies calculated under condition of definition of relative velocities of their points with respect to point 𝑂. It is easy to show that expression (3) is transformed to the more simple form if to calculate relative velocities of the specified points with respect to the instant position of the mass centre 𝐶 as follows: 𝑀𝑇=2𝐯2𝐶+𝑇𝐶𝑟+12𝝎×Θ𝐶𝝎+𝝎×𝐊𝐶𝑟.(4)

Note that though kinetic energy of the relative motion is a scalar, it follows from the previously mentioned that it is necessary to specify a chosen origin for determination of the relative velocities for its correct calculation.

Application of the Euler’s operator 𝐸𝑗()=(𝑑/𝑑𝑡)(𝜕()/𝜕̇𝑞𝑗)𝜕()/𝜕𝑞𝑗 to the expression (3) leads to the following Lagrange’s equations of the second kind for the generalized coordinates 𝑞𝑠: 𝐸𝑠𝑇𝑂𝑟𝑀𝐫𝐶𝜕𝐫𝐶𝜕𝑞𝑠12𝝎𝜕Θ𝐶𝜕𝑞𝑠+̇𝝎𝝎𝜕𝐊𝐶𝑟𝜕̇𝑞𝑠+𝝎𝐸𝑠𝐊𝐶𝑟=𝑄𝑠,(5) where 𝐸𝑗()=(𝜕/𝜕𝑡)(𝜕()/𝜕̇𝑞𝑗)𝜕()/𝜕𝑞𝑗 is also the Euler’s operator, but the time differentiation is performed here in the reference frame 𝑂𝑥𝑦𝑧; 𝑄𝑠 are the generalized forces that take into account the elastic and dissipative characteristics of the structure; 𝐦𝐶 is the gravitational torque.

It is important to mention here that the term 𝐸𝑠(𝐊𝐶𝑟) does not reduce to (2)𝜕𝐊𝐶𝑟/𝜕𝑞𝑠, as it is shown in [20] for 𝐸𝑠(𝐊𝑂𝑟), since the case to be considered here is the more general case of nonstationary connections (see expression (1)).

The following notations are used here: Θ𝐶 is the inertia tensor of the whole system with respect to the mass centre 𝐶, 𝐊𝐶𝑟=𝑚𝐫×𝐫𝑑𝑚𝑀𝐫𝐶×𝐫𝐶 is the relative moment of momentum of the deployed part with respect to point 𝐶, 𝐫𝐶 is the position vector of the instant position of the mass centre 𝐶 in the main body fixed frame of reference, and the symbol * indicates the time differentiation in the main body fixed reference frame.

Euler-Lagrange equation in the vector form can be obtained simple enough if to use the moment of momentum theorem of the whole system. The total moment of momentum of the system looks like 𝐊𝐶=Θ𝐶𝝎+𝐊𝐶𝑟.(6)

Because of application of the specified theorem, the SC attitude dynamics equation may be obtained in the following form: Θ𝐶̇𝝎+Θ𝐶Θ𝝎+𝝎×𝐶𝝎+𝝎×𝐊𝐶𝑟=𝐦𝐶,(7) which coincides with known equation [20].

If to supplement equations (5) and (7) by kinematic equations, one obtains a closed system of equations of motion. Now, one can pose the Cauchy problem and study the dynamics of the spacecraft in time, when its geometry changes. The parameters of Rodrigues-Hamilton were chosen as kinematic parameters, which determine the attitude of the main body in the orbital reference frame in order to anticipate possible singularities in the kinematical equations of the spacecraft while conducting the simulation. Reference [17] gives details of mutual transformations between quaternion and conventional attitude angles.

Further, it is necessary to choose proper generalized co-ordinates 𝑞𝑠 and determine expressions for 𝐫𝑐(𝑞1,𝑞2,,𝑞𝑁,𝑡), Θ𝑂(𝑞1,𝑞2,,𝑞𝑁,𝑡), 𝑇𝑂𝑟(𝑞1,𝑞2,,𝑞𝑁,̇𝑞1,̇𝑞2,,̇𝑞𝑁,𝑡), 𝐊𝑂𝑟(𝑞1,𝑞2,,𝑞𝑁,̇𝑞1,̇𝑞2,,̇𝑞𝑁,𝑡), Π(𝑞1,𝑞2,,𝑞𝑁,𝑡), and all their derivatives with respect to time and generalized co-ordinates, which appear in the expressions for coefficients of the governing equations (5) and (7).

Choice of the generalized co-ordinates plays the extremely important role at studying of the dynamics of complex mechanical systems. The successful choice allows to write the motion equations in the compact and transparent form and to conduct their effective integration. On the contrary, the irrational choice leads to bulky expressions for coefficients of the equations that stipulate ineffective integration (large time of calculations, necessity of significant capacities of the operative memory, accumulation of errors at evaluations, and “stiffness” of the system).

The SC state after its placing into orbit and completion of the preliminary damping in the orbital frame of reference, when the gyroscopes are uncaged, is considered as the initial system state.

In this case, it is impossible to consider that the design centreline remains rectilinear, as Coriolis forces act on its points during the deployment in conditions of the rotational-translational motion. These forces can produce deviations of an elastic design in the transversal direction in the orbital plane. At the same time, transversal forces are absent in the direction of the orbit binormal. Hence, the deviations in the direction of axis 𝑂𝑦 can be neglected.

Values of co-ordinates 𝑥𝑖,𝑧𝑖(𝑖=1,𝑁) of points 𝐸𝑖 of each tier, lying on axis 𝑂𝑧 on the straight lines connecting points 𝐸(𝐿)(𝑖) and 𝐸(𝑅)(𝑖) (Figure 4), and angles ̃𝛽𝑖(𝑖=1,2) (Figure 1) have been chosen as the generalized co-ordinates in this study. Note that 𝑧𝑖=𝑧𝑖,𝑝(𝑡)+𝑧𝑖,𝑒, where 𝑧𝑖,𝑝(𝑡) are functions of time, which are determined by the deployment time schedule, and 𝑧𝑖,𝑒 are independent variables.

Using Figure 4, it is simple to show the position vector of the mass centre of 𝑖th tier of the pantograph design is 𝐫𝑂(𝑖)=𝑥(1+𝑖)+𝑥(𝑖)2,0,(𝑧(1+𝑖)+𝑧(𝑖))2.(8)

As a result, for SC as a whole, 𝐫𝑂=𝑎𝑥,𝑎𝑦,𝑎𝑧𝑚𝐾𝐴+0,0,̃𝑧0𝑚0+𝑚1𝑀+𝑁𝑖=1𝑥𝑖,0,𝑧𝑖𝑚𝑖𝑀.(9)

Here, {𝑎𝑥,𝑎𝑦,𝑎𝑧} is the position vector of the mass centre of the main body in the frame of reference 𝑂𝑥𝑦𝑧; 𝑚𝐾𝐴 is mass of the main body; 𝑚𝑖(𝑖=0,𝑁) is mass of 𝑖th tier; {𝑥𝑖,0,𝑧𝑖} is the position vectors of points 𝐸𝑖; 𝑚𝑖=(𝑚𝑖+𝑚𝑖+1)/2, 𝑚𝑁=𝑚𝑁/2, ̃𝑧0=𝑧0𝑘0; 𝑘0 is the coefficient defining the position of the mass centre of the zero tier on axis 𝑂𝑧; 𝑏(𝑡)=𝐻0𝑂=𝑂𝐹0 is known function of time which determines the law of the design deployment.

2.4. Main Dynamical Values

The inertia tensor of the pantograph design can be written in the form of the sum of the local tensors of inertia of single rods in the frame of reference that is parallel to axis 𝑂𝑥𝑦𝑧 with the origin in the mass centre of corresponding rods, and the tensor of inertia of all point masses disposed in the centres of mass of the rods that are equal to masses of these rods. It gives the following expressions for 𝑖th tier: Θ𝑂1,1,𝑖=2𝐽𝑀𝑖+𝑎2𝑚𝑖sin2𝛼𝑖3+𝑎2𝑚𝑖sin2𝛽𝑖3+𝑚𝑀𝑖𝑑𝑧2𝑖2+2𝑚𝑀𝑖𝑑𝑥2𝑖𝑢2𝑖𝐿2𝑖2𝑚𝑀𝑖𝑑𝑧𝑖𝑧𝑖+2𝑚𝑀𝑖𝑧2𝑖+𝑚𝑖5𝑑𝑧𝑖22+2𝑑𝑥2𝑖𝑢2𝑖𝑒𝐿12𝑖8𝑑𝑧𝑖𝑧𝑖+8𝑏𝑏2+𝑧2𝑖,Θ𝑂1,3,𝑖=𝑎2𝑚𝑖sin2𝛼𝑖+sin2𝛽𝑖6+𝑚𝑀𝑖𝑥𝑖𝑑𝑧𝑖2𝑧𝑖+𝑑𝑥𝑖𝑑𝑧𝑖0.5+2𝑢2𝑖𝐿2𝑖+𝑧𝑖+𝑚𝑖4𝑥𝑖𝑑𝑧𝑖2𝑧𝑖+𝑑𝑥𝑖𝑑𝑧𝑖2.5+2𝑢2𝑖𝐿2𝑖+4𝑧𝑖,Θ𝑂2,2,𝑖=16𝐽𝑖𝐿2𝑖+𝑚𝑀𝑖𝑑𝑥2𝑖𝐿2𝑖+4𝑢2𝑖+𝑑𝑧2𝑖𝐿2𝑖+4𝑢2𝑖4𝑑𝑥𝑖𝐿2𝑖𝑥𝑖4𝑑𝑧𝑖𝐿2𝑖𝑧𝑖+4𝐿2𝑖𝑥2𝑖+𝑧2𝑖+𝑚𝑖𝑑𝑥2𝑖5𝐿2𝑖+4𝑢2𝑖+𝑑𝑧2𝑖5𝐿2𝑖+4𝑢2𝑖16𝑑𝑥𝑖𝐿𝑖2𝑥𝑖16𝑑𝑧𝑖𝐿2𝑖𝑧𝑖+16𝐿2𝑖𝑥2𝑖+𝑧2𝑖/2𝐿2𝑖,Θ𝑂3,3,𝑖=2𝐽𝑀𝑖+𝑎2𝑚𝑖cos2𝛼𝑖3+𝑎2𝑚𝑖cos2𝛽𝑖3+𝑚𝑀𝑖𝑑𝑥2𝑖2+2𝑚𝑀𝑖𝑑𝑧2𝑖𝑢2𝑖𝐿2𝑖2𝑚𝑀𝑖𝑑𝑥𝑖𝑥𝑖+2𝑚𝑀𝑖𝑥2𝑖+𝑚𝑖5𝑑𝑥2𝑖2+2𝑑𝑧2𝑖𝑢2𝑖𝐿2𝑖8𝑑𝑥𝑖𝑥𝑖+8𝑏𝑏2+𝑥2𝑖,Θ𝑂1,1,0=2𝐽𝑀0+4𝑚0𝑏𝑏2𝑚0𝑏2+𝑚0𝑎2+2𝑚𝑀0𝑧20+𝑚0𝑎2sin2𝛼06+𝑚0𝑎2sin2𝛽06,Θ𝑂1,3,0=𝑚0𝑎2sin2𝛼0+sin2𝛽0,Θ12𝑂2,2,0=4𝐽0+2𝑚𝑀0𝑏2+𝑚0𝑎2+2𝑚𝑀0𝑧20,Θ𝑂3,3,0=𝑚0+2𝑚𝑀0𝑏2+12𝐽𝑀0+24𝑚0𝑏𝑏2+𝑚0𝑎2cos2𝛼0+𝑚0𝑎2cos2𝛽06,(10) where Θ𝑂𝑘,𝑗,𝑖,(𝑘,𝑗=1,2,3) are (𝑘,𝑗) component of the inertia tensor for 𝑖th tier; 𝑚𝑖 is mass of one rod in the rhombus; 𝑚𝑀𝑖 is mass of one connecting rod; 𝐽𝑖 is proper moment of inertia of the rod with respect to its principal central axis that is parallel to axis 𝑂𝑦; 𝐽𝑀𝑖 is the characteristic moment of inertia of the connecting rod with respect to its principal central axis that is parallel to axis 𝑂𝑥; 𝑑𝑥𝑖=𝑥𝑖𝑥𝑖1; 𝑑𝑧𝑖=𝑧𝑖𝑧𝑖1; 𝑢𝑖 is length of the line segment 𝐻𝑖𝑂𝑖; 𝐿𝑖 is length of 𝐸𝑖𝐸𝑖1; angles of inclination of rods to axis 𝑂𝑥𝛼𝑖 and 𝛽𝑖 are shown in Figure 6 for the first tier; index 𝑖 shows that the given value concerns to 𝑖th tier.

The expression for the kinetic energy of relative motion may be presented as follows: 𝑇𝑟𝑂𝑖=128𝐽𝑖𝛼2𝑖+𝛽2𝑖𝐿6𝑖𝑢2𝑖+𝑚𝑀𝑖16𝑎2𝑑𝑥2𝑖𝑑𝑥𝑧1𝑡2𝑖+16𝑎2𝑑𝑥𝑧1𝑡2𝑖𝑑𝑧2𝑖64𝑎2𝑑𝑥𝑖𝑑𝑥1𝑡𝑖𝑑𝑥𝑧1𝑡𝑖𝐿2𝑖𝑢2𝑖64𝑎2𝑑𝑥𝑧1𝑡𝑖𝑑𝑧𝑖𝑑𝑧1𝑡𝑖𝐿2𝑖𝑢2𝑖+16𝐿4𝑖𝑢2𝑖4𝑎2𝑑𝑥1𝑡2𝑖+4𝑎2𝑑𝑧1𝑡2𝑖4𝑑𝑥1𝑡𝑖𝐿2𝑖𝑥1𝑡𝑖4𝑑𝑧1𝑡𝑖𝐿2𝑖𝑧1𝑡𝑖+4𝐿2𝑖𝑥1𝑡2𝑖+𝑧1𝑡2𝑖+𝑚𝑖16𝑎2𝑑𝑥2𝑖𝑑𝑥𝑧1𝑡2𝑖+16𝑎2𝑑𝑥𝑧1𝑡2𝑖𝑑𝑧2𝑖64𝑎2𝑑𝑥𝑖𝑑𝑥1𝑡𝑖𝑑𝑥𝑧1𝑡𝑖𝐿2𝑖𝑢2𝑖64𝑎2𝑑𝑥𝑧1𝑡𝑖𝑑𝑧𝑖𝑑𝑧1𝑡𝑖𝐿2𝑖𝑢2𝑖+16𝐿4𝑖𝑢2𝑖×𝑑𝑥1𝑡2𝑖5𝐿2𝑖+4𝑢2𝑖+𝑑𝑧1𝑡2𝑖5𝐿2𝑖+4𝑢2𝑖16𝑑𝑥1𝑡𝑖𝐿2𝑖𝑥1𝑡𝑖16𝑑𝑧1𝑡𝑖𝐿2𝑖𝑧1𝑡𝑖+16𝐿2𝑖𝑥1𝑡2𝑖+𝑧1𝑡2𝑖/64𝐿6𝑖𝑢2𝑖.(11) For the gyro-dampers: 𝑇(gd𝑖)𝑟=𝜔rot𝑖Θ𝑂𝑖rot𝑖𝜔rot𝑖2+𝜔𝑐𝑠𝑖Θ𝑂𝑖𝑐𝑠𝑖𝜔𝑐𝑠𝑖2,(𝑖=1,2).(12) Here, the following notations are used additionally: 𝑇𝑟𝑂𝑖 is the kinetic energy of the relative motion of 𝑖th tier when the relative velocities are evaluated in the frame of reference 𝑂𝑥𝑦𝑧; 𝑥1𝑡𝑖=𝜕𝑥𝑖/𝜕𝑡; 𝑧1𝑡𝑖=𝜕𝑧𝑖/𝜕𝑡; 𝑑𝑥1𝑡𝑖=𝜕(𝑑𝑥𝑖)/𝜕𝑡; 𝑑𝑧1𝑡𝑖=𝜕(𝑑𝑧𝑖)/𝜕𝑡;𝑑𝑥𝑧1𝑡𝑖=2(𝑑𝑥𝑖𝑑𝑥1𝑡𝑖+𝑑𝑧𝑖𝑑𝑧1𝑡𝑖); 𝑖 is the tier number; 𝜔rot𝑖̇𝛽={𝑖,Ωrot𝑖,0};𝜔𝑐𝑠𝑖̇𝛽={𝑖,0,0} in the GD reference frames.

Substituting all used notations into (11) and (12), one may apply the Eulerian operator to such expressions with respect to particular generalized co-ordinates taking into consideration the fact that the generalized co-ordinate and a velocity with index 𝑖 are included into the full expression for 𝑇𝑂𝑟 for the tiers 𝑖 and tier 𝑖1 (except the first tier since motion of elements of the zero tier depends explicitly on time but not on generalized co-ordinates).

Expression for projections of the angular momentum into axis 𝑂𝑦 for 𝑖th tier has been reduced to the following form: 𝐾𝑟𝑂𝑖=8𝐽𝑖𝑑𝑥1𝑡𝑖𝑑𝑧𝑖𝑑𝑥𝑖𝑑𝑧1𝑡𝑖𝐿2𝑖+𝑚𝑖𝑥1𝑡𝑖𝑑𝑧3𝑖+4𝑎2+4𝑑𝑥2𝑖+3𝑑𝑧2𝑖𝑧𝑖1𝑎+42+𝑑𝑥2𝑖5𝑑𝑧2𝑖𝑧𝑖+𝑥1𝑡𝑖1𝑑𝑧3𝑖+4𝑎2+4𝑑𝑥2𝑖+5𝑑𝑧2𝑖𝑧𝑖1+4𝑎2+4𝑑𝑥2𝑖+3𝑑𝑧2𝑖𝑧𝑖𝑑𝑥3𝑖+4𝑎2+5𝑑𝑥2𝑖+4𝑑𝑧2𝑖𝑥𝑖1+4𝑎2+3𝑑𝑥2𝑖+4𝑑𝑧2𝑖𝑥𝑖𝑧1𝑡𝑖1+𝑑𝑥3𝑖+4𝑎23𝑑𝑥2𝑖4𝑑𝑧2𝑖𝑥𝑖14𝑎2+5𝑑𝑥2𝑖+4𝑑𝑧2𝑖𝑥𝑖𝑧1𝑡𝑖/2𝐿2𝑖+𝑚𝑀𝑖𝑥1𝑡𝑖𝑑𝑧3𝑖+4𝑎2+𝑑𝑥2𝑖+𝐿2𝑖𝑧𝑖1+4𝑎2𝑑𝑥2𝑖+𝐿2𝑖𝑧𝑖+𝑥1𝑡𝑖1𝑑𝑧3𝑖+4𝑎2𝑑𝑥2𝑖+𝐿2𝑖𝑧𝑖1+4𝑎2+𝑑𝑥2𝑖+𝐿2𝑖𝑧𝑖𝑑𝑥3𝑖+4𝑎2𝑑𝑧2𝑖+𝐿2𝑖𝑥𝑖1+4𝑎2+𝑑𝑧2𝑖+𝐿2𝑖𝑥𝑖𝑧1𝑡𝑖1+𝑑𝑥3𝑖+4𝑎2𝑑𝑧2𝑖𝐿2𝑖𝑥𝑖1+4𝑎2+𝑑𝑧2𝑖𝐿2𝑖𝑥𝑖𝑧1𝑡𝑖/2𝐿2𝑖.(13)

Two remaining components of 𝐊𝑂𝑟 are equal to zero in the considered problem as all relative motions happen in the planes that are parallel to the plane 𝑂𝑥𝑧.

For 𝑖th gyro-damper 𝐾𝑂𝑟(𝑖)=𝐓𝑐_gd(𝑖){𝐽𝑐𝑠̇̃𝛽𝑖,𝐻𝑔,0}(𝑖=1,2). The transformation matrix 𝐓𝑐_gd(𝑖) has the following form: |||||̃𝛽cos𝜎sin𝑖+̃𝛽𝑖0̃𝛽sin𝜎cos𝑖+̃𝛽𝑖0̃𝛽sin𝜎0cos𝑖+̃𝛽𝑖0̃𝛽sin𝑖+̃𝛽𝑖0̃𝛽sin𝜎cos𝜎sin𝑖+̃𝛽𝑖0̃𝛽cos𝑖+̃𝛽𝑖0|||||cos𝜎.(14) It is defined by the sequence of rotations of the frames of reference fixed to the GD cases around the axes 𝑂𝑥𝑦𝑧 (see Figure 7).

The potential energy of elastic displacements of the pantograph structure for 𝑖th tier looks like Π𝑖=𝑐joint𝛼𝑖𝛼𝑖12+𝛼𝑖+1𝛼𝑖2+𝛼𝑖+2𝛼𝑖+12+𝛽𝑖𝛽𝑖12+𝛽𝑖+1𝛽𝑖2+𝛽𝑖+2𝛽𝑖+12.(15) Here, 𝛼 and 𝛽 with subscripts are angles that are shown in Figure 6. Naturally, this expression must be corrected for two last tiers.

The generalized forces of flexible displacements may be determined as 𝑄𝑥𝑖=𝜕Π𝑖/𝜕𝑥𝑖, 𝑄𝑧𝑖=𝜕Π𝑖/𝜕𝑧𝑖.

Except the generalized forces of elastic displacements, the design is subjected to the forces of Coulomb friction in its real hinges and damping torques in the fictitious hinges, which have appeared because of finite-element modelling of elastic rods according to Figure 5. There are six real hinges (four side hinges and two central ones) in each tier. For simplicity, it is possible to consider that moments of a Coulomb friction in such hinges are identical, as exact quantitative estimations cannot be determined. Actually, these moments are dependent on time at longitudinal oscillations of the design. If to write the virtual work for all real hinges of the design as 𝛿𝐴𝑖=6𝑀Coul(𝛿𝛼𝑖+𝛿𝛽𝑖)sign(̇𝛼𝑖+̇𝛽𝑖) and select coefficients at variations of generalized co-ordinates in the expression of the total virtual work, one can obtain the expressions for the generalized forces corresponding to the Coulomb friction.

The information on structural damping hardly can be enough exacted, nevertheless a qualitative evaluation needs to be spent. Considering that the fictitious hinge in the equivalent system (in Figure 5 below) includes not only a spring but also a viscous damper that corresponds to Kelvin-Voigt model, it is possible to write expressions for a dissipative function of four fictitious hinges of each a tier as follows: 𝑅𝑖=𝑟joinṫ𝛼𝑖̇𝛼𝑖12+̇𝛼𝑖+1̇𝛼𝑖2+̇𝛼𝑖+2̇𝛼𝑖+12+̇𝛽𝑖̇𝛽𝑖12+̇𝛽𝑖+1̇𝛽𝑖2+̇𝛽𝑖+2̇𝛽𝑖+12.(16)

The approximate value of the dissipative parameter 𝑟joint can be defined from the Lagrange equation for a separate equivalent design, if to reduce it to its classical form ̈̇Δ+2𝑛Δ+𝜔20Δ=0, where 𝑛=3𝑟joint/(𝑚𝑎2). Taking into account the expression for the defined above spring stiffness in the hinge, one can obtain the expression 𝑟joint=𝛿𝑎2𝑐joint𝑚/3(4𝜋2𝛿2)(𝛿𝑎/2𝜋)2𝑐joint𝑚/3. Here 𝛿 is a decrement of oscillations. Its value can be accepted such, as for steel designs. As a result, it is possible to obtain the expressions for the generalized forces of the structural damping in the fictitious hinges in the following form: 𝑄𝑅𝑥𝑖=𝜕𝑗𝑅𝑗/𝜕𝑥𝑖,𝑄𝑅𝑧𝑖=𝜕𝑗𝑅𝑗/𝜕𝑧𝑖. Corresponding analytical expressions for all main dynamic values have been obtained in Mathematica 5.

The potential energy of the GD springs has the following form: ΠGD=̃𝛽1̃𝛽0122𝑘4𝛽1+̃𝛽2̃𝛽0222𝑘4𝛽2,(17) where 𝑘4𝛽𝑖(𝑖=1,2) is the GD springs resistance, and ̃𝛽0𝑖(𝑖=1,2) are the values that ̃𝛽𝑖 possess, if the frame of reference 𝑂𝑥𝑦𝑧 does not rotate in the inertial basis.

For definition of coefficients of the motion equations (5) and (7), the obtained vector and tensor expressions and all their derivatives should be transformed to the frame of reference that is fixed in an instant position of the mass centre with its axes parallel to the 𝑂𝑥𝑦𝑧 axes with the same names.

According to [20], it is possible to write Θ𝐶=Θ𝑂𝑀𝐸𝑟𝑐2𝑟𝑐𝑟𝑐,𝐾𝐶𝑟=𝐾𝑂𝑟𝑀𝑟𝑐×𝑟𝑐,Θ𝐶=Θ𝑂𝑀2𝐸𝑟𝑐𝑟𝑐𝑟𝑐𝑟𝑟𝑟𝑐,𝜕Θ𝐶𝜕𝑞𝑠=𝜕Θ𝑂𝜕𝑞𝑠2𝑀𝐸𝑟𝑐𝜕𝑟𝑐𝜕𝑞𝑠+2𝑀𝑟𝑐𝜕𝑟𝑐𝜕𝑞𝑠,𝐾𝐶𝑟=𝐾𝑂𝑟𝑀𝑟𝑐×𝑟𝑐.(18)

Here, 𝐸 is the identity matrix; the expressions 𝑎𝑏 and 𝑎×𝑏 mean scalar and vector products of vectors accordingly; the expression 𝑎𝑏 means a dyadic product.

Having substituted the obtained expressions in (5) and (7) after completing them with suitable kinematic equations, one can formulate the Cauchy problem for study the of dynamics of the system under consideration during its deployment.

3. Simulation

For numerical integration, the system is written in the following implicit form: 𝑁eq𝑗=1𝐴𝑖𝑗𝑦1,𝑦2,,𝑦𝑁eq,𝑡̇𝑦𝑗=𝑓𝑖𝑦1,𝑦2,,𝑦𝑁eq,𝑡𝑖=1,𝑁eq.(19)

In the case under consideration, 𝑦𝑖=𝜔𝑖(𝑖=1,2,3). If the problem is written in the quaternions, 𝑦𝑖(𝑖=4,7) define the quaternion components. If the problem is written in the Bryan angles, then the same variables define the angles (in this case 𝑦70). Remaining variables are the generalized co-ordinates and their first derivatives with respect to time, that is, 𝑦4+4𝑖=𝑥𝑖,𝑦5+4𝑖=𝑥1𝑡𝑖,𝑦6+4𝑖=𝑧𝑖,𝑦7+4𝑖=𝑧1𝑡𝑖𝑖=,𝑦1,𝑁8+4𝑁=̃𝛽1,𝑦9+4𝑁=̇̃𝛽1,𝑦10+4𝑁=̃𝛽2,𝑦11+4𝑁=̇̃𝛽2.(20)

So, the total number of variables is 𝑁eq=11+4𝑁. The original computation package is developed for numerical integration of the ordinary differential equations in the frame of the Cauchy problem. The majority of operators of the program is obtained in the form of Fortran expressions in Mathematica 5 on the basis of realization of some packages written especially for the studied system, and after replacement of the obtained bulky expressions with simple enough variables.

The pantograph structures having 5, 10, and 20 tiers were studied. The Cauchy problem integration step was chosen from a condition of stability of the computation process on a studied interval of time. Its magnitude varied in the range of 0.01–0.001 s. Initial values of components of 𝝎 differ from the components of orbital angular velocity on ±0.0001 s−1.

It is obvious that the dynamics of a complex flexible structure depend on the time history of the deployment (low of the deployment). Two such time histories were used here. The deployment of space designs often used such a law 𝑏(𝑡) that the velocity and acceleration time history of an actuating motor ̇𝑏(𝑡) looks like the line with the mark “stand” in Figure 8. Such a law is used conventionally at the deployment of such space structures as solar batteries. Such a function ̇𝑏(𝑡) stand has two angular points, and function ̈𝑏(𝑡) stand (and a function of an actuating force or torque) has four points of discontinuities. Action of such a law on an oscillating system brings additional perturbations in its dynamics.

Therefore, the alternative law of deployment of the design was considered. For this law, ̈𝑏(𝑡) does not contain points of discontinuities. Such a law can be constructed as a solution of the optimal control problem of rest-to-rest transfer of a mechanical system with internal oscillating degrees of freedom at an integral minimisation of relative accelerations of the generalised co-ordinates (see [22], Special Case II). The velocity time history of an actuating motor ̇𝑏(𝑡) opt and ̈𝑏(𝑡) opt for the second law looks like the line 2 in Figure 8.

The orbit inclination 𝑖 is assumed to be 1, the longitude Ω of the ascending node is taken to be zero. The moment of the central Newtonian field, corresponding to a circular orbit of 600 km altitude, was considered as the external perturbing moment. Though the SC movement along an orbit is not considered here, the orbit parameters are required to calculate the gravitational torque and projections of the total SC moment of momentum to the absolute frame of reference. It is necessary for the monitoring of errors of the numerical integration of the initial value problem.

Key system parameters are shown in the Table 1.


Mass of main body 𝑚 1 350 kg
Rods mass 𝑚 𝑙 1 kg
Rods bending stiffness 𝐸 𝐽 20–80 N m2
Decrement of oscillations 𝜗 0.001
Nonzero components of the main body inertia tensor 𝐽 𝑥 𝑥 4000 kg m2
𝐽 𝑦 𝑦 5000 kg m2
𝐽 𝑧 𝑧 2000 kg m2
𝑎 𝑥 0 m
Shift of main body centre of mass 𝑎 𝑦 0 m
𝑎 𝑧 −1 m
Angular momentum of one GD rotor 𝑟 𝑜 𝑡 20 kg m2/s
GD damping coefficient 𝑘 3 𝛽 1 , 𝛽 2 40 N m/s2
Duration of deployment 𝑇 𝑓 100–1000 s

A series of evaluations have been conducted at different values of system parameters in the process of the numerical simulation of the deployment dynamics of the pantograph design. The attitude modification of the SC main body under the effect of the structure deployment and dynamics of a modification of the shape of the construction depending on the adopted parameters were studied.

Design deployment essentially increases the components of the inertia tensor Θ𝐶1,1 and Θ𝐶2,2 and decreases slightly the component Θ𝐶3,3. Figure 9 shows the time history of the diagonal components of the total inertia tensor of SC for both standard and optimal laws of deployment. Here, one can see that Θ1,1 and Θ2,2, which differ at the beginning, then behave practically synchronously with oscillations of the pantograph length. Component Θ3,3 decreases as the distances of pantograph’s points to axis Oz decrease noticeably during the deployment. One can see that oscillation amplitudes of the pantograph length are essentially less when optimum laws are used. Generally speaking, the inertia tensor is not diagonal one in the presence of transversal design deviations in the direction of axis 𝑂𝑥. The nonzero component Θ1,3 appears at that, but in view of smallness of these deviations as compared to the longitudinal displacements, it is negligibly small also as compared to the diagonal components of the inertia tensor and has no essential effect on the SC dynamics.

The increase in the inertia tensor components leads to sharp decrease in the component 𝜔2 of the SC angular velocity vector (Figure 10) and naturally to loss of the SC orientation in the orbital frame of reference. At the same time, generalized co-ordinates 𝑧𝑖(𝑖=4,8,12,16,20) behave according to Figure 11. Such behaviour speaks about appreciable longitudinal oscillations of the design at the deployment stage. These oscillating components are induced by excitation of elastic oscillations of design rods. Their amplitudes grow with increase of number of the tier. Longitudinal oscillations have noticeable influence upon the components of the inertia tensor (Figure 9). Behaviour of angles of the rods tilt 𝛼𝑖 with the same values of the subscripts is shown in Figure 12. Here, one can see that in the beginning of the deployment, magnitudes of the angles become less than their initial values under the effect of forces of rest inertia. As far as the relative acceleration, which the deploying mechanism tries to impose to the upper tiers, is greater than that for the lower ones, their “draw” is more appreciable.

The transversal relative deviations of the design longitudinal axis defined by the generalized co-ordinates 𝑥𝑖(𝑖=4,8,12,16,20) are shown in Figure 13. Practically, all these deviations are directed opposite to the vector of the SC velocity under the influence of the Coriolis forces. Figures 9, 10, 11, 12, and 13 correspond to the standard law of deployment. If to use the optimal law mentioned above, at the same values of the remaining parameters, the amplitudes of oscillations 𝑧𝑖 are essentially smaller than what one can see in Figure 14.

One of the considered variants has allowed to analyze the dynamics of the system at the same flexural rigidity, but at twice as long deployment, that is, at 𝑇𝑓 =1000 s. The behaviour of the generalized co-ordinates becomes smooth enough in this case, amplitudes of oscillations of the design length do not exceed 0.15 m. The transversal deviations have a smooth mode. The design deviates like a cantilever beam. The oscillating components are superimposed on these deviations. These deviations tend to reduction of their velocity with time more appreciable than in the previous case as duration of deployment is twice as long here.

Study of SC behaviour with the five-tier pantograph on long intervals of time has shown that SC with switched-off GD under the influence of the gravitational moment enters into a condition of simple harmonic pitch oscillations. Taking into account the Coulomb friction in design, hinges have no effect on damping of SC attitude oscillations practically as elastic longitudinal oscillations of the design damp quickly enough as compared with duration of deployment. Design transverse vibrations, on the contrary, damp very slowly even under the influence of forces of structural damping as their amplitudes and velocities are very small as one can see in Figure 15.

At deployment of this design from very heavy ideally stabilised space station (when the mass of the design can be neglected in comparison with the total mass of the station, and elastic oscillations have no effect on its attitude) at deployment within 100 seconds, the amplitudes and frequencies of longitudinal and transversal oscillations visibly decrease.

The analysis of the SC dynamics with the switched-on GD has shown that the pitch SC oscillations damp with sufficient for the practice speed at the accepted values of parameters of the GD system (Figure 16). The fact is noteworthy that design transverse vibrations interact with the pitch SC oscillations and damp together with them even if the structural damping is not taken into consideration (Figure 17). The behaviour of angles of GD deviation ̃𝛽𝑖 in the process of stabilization of the pitch oscillations is shown in Figure 18.

Obtained results were compared with the case of increasing rods rigidity up to the rigid rods, the behaviour of which is described in [19]. Here asymptotic coincidence of the results was obtained. In general, the control of the total moment of momentum magnitude was realized to confirm the results. According to the known theorem in [20], 𝐊𝐶(𝑡)=𝐊𝐶(𝑡0)+𝑡𝑡0𝐦𝐶(𝜏)𝑑𝜏. The term in the left part is calculated on each step of evaluations as the total current moment of momentum of the system under consideration. The first term in the right part is calculated using the initial conditions, and the second one is calculated on each step of the integrations. Good coincidence of the left and right part values indicates the absence of errors in the motion equations and in the computing program. Such a check has been included in the computing programme. During monitoring, 8 significant figures coincided for each projection, as a rule.

4. Conclusions

The present paper deals with the exploration of the dynamics of the gyro-gravitational stabilized spacecraft in the mode of the initial deployment of the elastic structure after the completion of the preliminary stabilization of the disordered rotation of the SC around its mass centre. Novelty of the approach consists in taking into account the additional internal degrees of freedom of the pantograph design in comparison with known earlier settings of the problem; in using the method of derivation of the dynamic equations of mechanical systems with internal degrees of freedom and nonstationary connections developed by author; in using optimum with respect to damping of elastic oscillations control profile for deployment of flexible designs. The mathematical model developed for this case may be regarded as the generalization of the theory of a flexible multibody system with the time-dependent configurations. The approach may be successfully extended to the modelling of the deployment dynamics of other space structures with the significant change of the configuration in the process of their functioning. The computational Fortran-package developed for the numeric simulation and the code for Mathematica 5© developed for the symbolic calculation have general characteristics, which may be easily adopted for other deployed systems. The obtained data illustrate the method of mathematical formulation of the system dynamics with the time-dependent configuration as well as the behaviour of the SC in the process of deployment of the flexible structures. The information to be obtained in the process of the numerical simulation has allowed drawing the conclusion that deployment of the studied design may be comprehensible to practice only at a rational choice of the parameters of the SC and its deployment law. Obtained data permit to select the most appropriate deployment and gyro-dampers parameters.


  1. W. E. Lang and G. N. Honeycutt, “Simulation of deployment dynamics of spinning spacecraft,” NASA TN D-4074, 1967. View at: Google Scholar
  2. G. J. Cloutier, “Dynamics of deployment of extendible booms from spinning space vehicles,” Journal of Spacecraft and Rockets, vol. 5, pp. 547–552, 1968. View at: Publisher Site | Google Scholar
  3. E. J. Bowers Jr. and C. E. Williams, “Optimization of RAE satellite boom deployment timing,” Journal of Spacecraft and Rockets, vol. 7, no. 9, pp. 1057–1062, 1970. View at: Google Scholar
  4. R. Barakat, “Transverse vibrations of a moving thin rod,” Journal of the Acoustical Society of America, vol. 43, pp. 533–539, 1968. View at: Publisher Site | Google Scholar
  5. B. Tabarrok, C. M. Leech, and Y. I. Kim, “On the dynamics of an axially moving beam,” Journal of the Franklin Institute, vol. 297, no. 3, pp. 201–220, 1974. View at: Google Scholar
  6. A. K. Banerjee and T. R. Kane, “Extrusion of a beam from a rotating base,” Journal of Guidance, Control, and Dynamics, vol. 12, no. 2, pp. 140–146, 1989. View at: Google Scholar
  7. N. G. Creamer, “Deployment of a flexible beam from an oscillating base,” Journal of Guidance, Control, and Dynamics, vol. 15, no. 2, pp. 527–529, 1992. View at: Google Scholar
  8. J. A. Lewis and E. E. Zajac, “A two-gyro, gravity-gradient satellite attitude control system,” Bell System Technical Journal, vol. 43, pp. 2705–2765, 1964. View at: Google Scholar
  9. D. B. Cherchas, “Dynamics of spin-stabilized satellites during extension of long flexible booms,” Journal of Spacecraft and Rockets, vol. 8, no. 7, pp. 802–804, 1971. View at: Google Scholar
  10. K. W. Lips and V. J. Modi, “Three-dimensional response characteristics for spacecraft with deploying flexible appendages,” Journal of guidance and control, vol. 4, no. 6, pp. 650–656, 1981. View at: Google Scholar
  11. K. Tsuchiya, “Dynamics of a spacecraft during extension of flexible appendages,” Journal of Guidance, Control, and Dynamics, vol. 6, no. 2, pp. 100–103, 1983. View at: Google Scholar
  12. K. W. Lips, W. B. Graham, F. R. Vigneron, and D. G. Hunter, “Dynamics and control characteristics for the WISP 300 m dipole antenna/shuttle configuration,” AAS Paperno. 8, pp. 85–365, 1985. View at: Google Scholar
  13. A. K. Banerjee, “Order-n formulation of extrusion of a beam with large bending and rotation,” Journal of Guidance, Control, and Dynamics, vol. 15, no. 1, pp. 121–127, 1992. View at: Google Scholar
  14. A. K. Banerjee and S. Nagarajan, “Efficient simulation of large overall motion of beams undergoing large deflection,” Multibody System Dynamics, vol. 1, no. 1, pp. 113–126, 1997. View at: Google Scholar
  15. E. M. Levin, Dynamic Analysis of Space Tether Missions, American Astronautical Society by Univelt, 2007.
  16. P. Williams, “Predictive Control of Tethered Satellite Systems, Model Predictive Control, Tao Zheng (Ed.),” InTech, 2010, View at: Google Scholar
  17. A. E. Zakrzhevskii and V. S. Khoroshilov, “Dynamics of the orbital deployment of an elastic ring-shaped antenna,” International Applied Mechanics, vol. 46, no. 6, pp. 718–729, 2010. View at: Publisher Site | Google Scholar
  18. V. I. Dranovskii, V. S. Khoroshylov, and A. E. Zakrzhevskii, “Spacecraft dynamics with regard to elastic gravitational stabilizer deployment,” Acta Astronautica, vol. 64, no. 5-6, pp. 501–513, 2009. View at: Publisher Site | Google Scholar
  19. V. I. Dranovskii, A. E. Zakrzhevskii, A. P. Kovalenko, and V. S. Khoroshilov, “On the dynamics of deployment of an orbital structure with elasticelements,” International Applied Mechanics, vol. 42, no. 8, pp. 959–965, 2006. View at: Publisher Site | Google Scholar
  20. A. I. Lurie, Analytical Mechanics, Springer, 2002.
  21. V. V. Beletsky, “Motion of an artificial satellite about its center of mass,” Israel Program For Scientific Translations, Jerusalem NASA-TT-F-429, TT-67-51366, 1966. View at: Google Scholar
  22. A. E. Zakrzhevskii, “Slewing of flexible spacecraft with minimal relative flexible acceleration,” Journal of Guidance, Control, and Dynamics, vol. 31, no. 3, pp. 563–570, 2008. View at: Publisher Site | Google Scholar

Copyright © 2012 V. S. Khoroshilov and A. E. Zakrzhevskii. 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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.