# Dynamics and control of a tethered satellite system based on the SDRE method.

1. IntroductionA tethered satellite system (TSS) usually consists of several spacecraft, a long tether which connects these spacecraft and a tether deployment/retrieval system. In 1895, Tsiolkovsky conceived the concept of a TSS [1]. In 1974, Colombo et al. proposed a "Shuttle-borne Skyhook" for low-orbital altitude research, which marked the advent of tethered satellite systems (TSS) [2]. Since then, there are several potential applications in space technologies, such as scientific experiments, deployment or retrieval of satellites, power generation, measurement of aerodynamic forces.

A classical TSS consists of two satellites connected by a tether, and one of the simplest models is a one-piece dumbbell, where the two satellites are assumed to be point masses and the tether is considered to be a straight massless line [3]. He et al. studied the stability of a TSS by utilizing Taylor expansion formula and stability criterion, and a range-rate control algorithm is developed to achieve the stable control of tether's deploying, keeping, and retrieving [4]. Nakanishi et al. investigated the in-plane periodic solutions for a dumbbell model in elliptic orbits using bifurcation, and their trajectories are projected on the van der Pol planes to predict when the delayed feedback control will need to act to maintain the periodic motions [5]. Kojima et al. conducted an experimental study to investigate the effect of gravity for the libration of tethered satellite system by delayed feedback control [6]. Zhang et al. presented several criteria on the existence of periodic solutions for a TSS in an elliptical orbit, and the uniqueness of periodic solutions for the TSS in a circular orbit is presented on the basis of coincidence degree theory [7]. Hu et al. presented the theoretical and experimental studies of the deployment and retrieval control of a TSS, and the ideas of online optimization and receding horizon control were applied to design a feedback controller for the TSS [8]. Jung et al. presented a dynamic analysis of a TSS with a moving mass based on a two-piece dumbbell model, where the moving mass is conveyed between two satellites along a straight massless tether [9]. The same authors analyzed the nonlinear dynamic behavior of deployment and retrieval of a three-body TSS on a variable-radius orbit, and the system is modeled as a two-piece dumbbell model with six degrees of freedom, which consists of three point masses and two straight massless tethers [10]. Hong and Varatharajoo discussed the development of mathematical models for a flexible TSS in both planar and coplanar states, where the TSS consists of three rigid bodies with two flexible tethers, each connecting two rigid bodies with one located in the center and serving as the mothership [11]. Hong et al. developed the analytic solutions for a TSS subjected to internal tether tension moment and external aerodynamic torque for spin-up and spin-down maneuvers by the approximation of Euler's equations of motion via Fresnel integrals and sine and cosine integrals [12]. Furthermore, there are some review papers addressing the technologies and applications of a TSS [13-17].

For nonlinear optimal control, Pearson presented a linear time and state-dependent approximation to a nonlinear and nonstationary system optimized with respect to a quadratic performance index by treating it as an instantaneously linear stationary system in 1962 [18]. This approach was further investigated by Coultier et al., who extended the LQR approach to propose a so-called state-dependent Riccati equation (SDRE) method for the nonlinear quadratic regulation problems [19]. Zhao and Deng investigated the relationship between the technique by SDRE and Hamilton-Jacobi-Isaacs equations for nonlinear [H.sup.[infinity]] control design by establishing Lyapunov matrix equations for partial derivatives of the solution of the SDREs and introducing symmetry measure for some related matrices [20]. Nekoo investigated a general case which has control nonlinearities and time-varying weighting matrix Q is solved with three approaches: exact solution, online control update, and power series approximation [21]. Chang and Bentsman extended the continuous time SDRE technique to discrete time under input and state constraints, yielding constrained discrete time SDRE, referred to as CD-SDRE [22]. Babaei and Salamci presented a model reference adaptive control algorithm for a class of nonlinear systems by using a SDRE technique [23]. Nekoo and Geranmehr used a combination of an SDRE controller and SDRE observer to control a class of nonlinear nonaffine control systems [24]. Wang et al. presented a generalized SDRE approach for continuous time nonlinear systems to achieve a mixed nonlinear quadratic regulator and [H.sup.[infinity]] control performance criteria [25]. Korayem and Nekoo used the finite-time horizon SDRE for controlling a class of nonlinear systems, and the derivations of the SDRE for two classes of systems were presented in the finite-time horizon [26]. Korayem and Nekoo investigated the finite-time optimal and suboptimal controls for time-varying systems with state and control nonlinearities, and a finite-time constraint imposed on the equation leads to a state-dependent differential Riccati equation (SDDRE) [27]. Furthermore, there are some review papers addressing the SDRE method [28-32].

This paper investigated a TSS modeled by a two-piece dumbbell model, which leads to five nonlinear coupled ordinary differential equations. There are two sets of equations of motion formulated. One is based on the first satellite and the other is based on the mass center of the TSS. There are two reasons to formulate the two sets of equations. One is to facilitate their mutual comparison due to the complex formulations. The other is to provide them for different application situations. In order to suppress the librations of the system, the SDRE control is applied, which aims to control a nonlinear system and provides the flexibility of controller design. Besides, to increase the convergence rates of the librations, the prescribed degree of stability is incorporated into the SDRE control. This paper is organized as follows. Section 2 formulates the two sets of equations of motion of the TSS. Section 3 presents the SDRE control with the prescribed degree of stability. Section 4 analyzes the nonlinear dynamics of the TSS. Section 5 demonstrates the libration suppressions by the SDRE control. Section 6 summarizes some significant conclusions.

2. Dynamic Modeling

A tethered satellite system consists of three satellites connected by a tether, where each satellite is considered as a point mass. Besides, the tether is assumed to be massless, and its length is constant. Thus, the system can be simplified as a two-piece dumbbell model shown in Figure 1, where the three points A, B, and D represent three satellites, and they are connected to each other by a constant-length massless tether. These satellites are assumed to be point masses and to be expressed as [m.sub.A], [m.sub.B], and [m.sub.D], which refer to the points A, B, and D, respectively. Besides, there is a fixed coordinate system X-Y defined on the center of Earth O. In this paper, one neglects the environmental disturbance torques due to the aerodynamics, the solar radiation, the Earth magnetic filed, and so forth.

This section presents two dynamic models to describe the motions of the tethered satellite system. The formulation of the first one is based on satellite A, and the other is based on the mass center of the system. Since both formulations describe the same system, they must have the same degrees of freedom, which is equal to five in the system. However, both formulations are expressed in terms of different variables. Each formulation can be used to perform the analysis for a given specific prescribed motion. For instance, if a circular motion is imposed on satellite A, the first set of formulations should be adapted. If the same circular motion is imposed on the mass center of the system, the second set of formulations should be selected.

2.1. Formulations Based on Point A. Figure 2(a) shows the coordinate system defined at point A, and the position of point A is expressed as the orbital radius vector [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] with the true anomaly [[psi].sub.A]. Besides, another coordinate system fixed on point A is defined by the two unit vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The relative position vectors from points A to D and from D to B are, respectively, expressed as [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], which can show the configuration of the tether. Referring to Figure 2(a), angles [[phi].sub.1] and [[phi].sub.2], respectively, represent the directions of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] based on the unit vector [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The two angles are introduced to investigate the libration of the tethered satellite system based on different coordinate systems. If they have an equal value, the system can be treated as a one-piece dumbbell model.

Referring to Figure 2(a), the position vectors of the three point masses from the origin of the X-Y coordinate system are given as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

where [R.sub.A], [R.sub.B], [R.sub.D], [L.sub.1], and [L.sub.2] represent the length of vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], respectively; [c.sub.1] and [c.sub.2] refer to cos [[phi].sub.1] and cos [[phi].sub.2], respectively. Note that the three position vectors are expressed in terms of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

By differentiating (1) with respect to time, the velocities of the three vectors are given as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

Therefore, the kinetic energy and the potential energy due to the Earth's gravitational field are, respectively, written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

where G and [M.sub.e] are the universal gravitational constant and the Earth mass, respectively.

The nonlinear dynamic equations of motion of the tethered satellite system can be obtained by applying Euler-Lagrange's equation as

d/dt ([partial derivative][L.sub.g]/[partial derivative][[??].sub.i]) - [partial derivative][L.sub.h]/[partial derivative][q.sub.i] = [Q.sub.i], (5)

where [q.sub.i] is the generalized coordinate, [Q.sub.i] is the nonconservative forces corresponding to the generalized coordinate, and [L.sub.g] is the Lagrangian as

[L.sub.g] = T - V (6)

in which T and V are the kinetic energy and the potential energy, respectively. They are expressed in (3) and (4).

Based on the two-piece dumbbell model, the generalized coordinates are [R.sub.A], [[psi].sub.A], [[phi].sub.1], [[phi].sub.2], and [L.sub.1]. Besides, the nonconservative forces [Q.sub.i] refer to the driving forces or torques applied to some points and the tether. Substituting (1) to (4) and (6) into (5), one has the nonlinear dynamic equations of motion as follows.

For [q.sub.i] = [R.sub.A], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)

where [F.sub.A] is the external force applied to point A. One has

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

Note that (9) to (11) are expressed based on the coordinate system, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

For [q.sub.i] = [[psi].sub.A], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)

where [T.sub.O] is the external torque applied to point O and [c.sub.12] and [s.sub.12] refer to cos([[phi].sub.1] - [[phi].sub.2]) and sin([[phi].sub.1] - [[phi].sub.2]), respectively.

For [q.sub.i] = [[phi].sub.1], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (13)

where [T.sub.A] is the external torque applied to point A.

For [q.sub.i] = [[phi].sub.2], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (14)

where [T.sub.D] is the external torque applied to point D.

For [q.sub.i] = [L.sub.1], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

where [F.sub.L] is the external force applied to the tether.

Equations (7) and (12) to (15) are the nonlinear dynamic equations of motion for the tethered satellite system based on the coordinate system [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

2.2. Formulations Based on the Mass Center. Figure 2(b) shows the coordinate system defined at the mass center, where point C represents the mass center of the tethered satellite system, and the position vector from points O to C is expressed as [[??].sub.C] with the true anomaly [[psi].sub.C]. The relative position vectors from points C to the three point masses are [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Referring

to Figure 2(b), another coordinate system fixed on point C is defined by the two unit vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Referring to Figure 2(b), angles [[theta].sub.1] and [[theta].sub.2], respectively, represent the directions of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] based on the unit vector [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Since the total length of the tether is a constant value, it can be written as

L = [L.sub.1] + [L.sub.2], (16)

where [L.sub.1] and [L.sub.2] are the lengths of vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], respectively.

To formulate the dynamic equations of the two-piece dumbbell model, one will apply Euler-Lagrange's equation, and the position and velocity of each point mass will be expressed in terms of the relative vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Furthermore, they will be written based on the coordinate system ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]).

Referring to Figure 2(b), one has these two equations as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17)

Besides, since point C is the mass center of the tethered satellite system, one has the equation as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18)

where the mass ratios are written as

[[rho].sub.A] = [m.sub.A]/M, [[rho].sub.B] = [m.sub.B]/M, [[rho].sub.D] = [m.sub.D]/M, (19)

in which the total mass of the system is given by

M = [m.sub.A] + [m.sub.B] + [m.sub.D]. (20)

Simultaneously solving (17) and (18) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

Note that (21) expresses the relative vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], in terms of the relative vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

By taking the time derivatives of (21), the velocities of the relative vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], are given as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22)

Referring to Figure 2(b) again, the relative vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] can be expressed in the coordinate system ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]) as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

where [c.sub.1], [c.sub.2], [s.sub.1], and [s.sub.2] refer to cos [[theta].sub.1], cos [[theta].sub.2], sin [[theta].sub.1], and sin [[theta].sub.2], respectively.

By taking the time derivatives of (23), the velocity of the relative vectors [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is expressed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

To obtain the dynamic equations of the system by applying Euler-Lagrange's equation, the kinetic energy is written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (26)

the potential energy is written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (27)

Using Euler-Lagrange's equation shown in (5), one obtains the equations of motion as follows.

For [q.sub.i] = [R.sub.C], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (28)

where [F.sub.C] is the external force applied to the center of the system and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (29)

For [q.sub.i] = [[psi].sub.C], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (30)

where [T.sub.O] is the external torque at point O, [c.sub.12] = cos([[theta].sub.1] - [[theta].sub.2]), and [s.sub.12] = sin([[theta].sub.1] - [[theta].sub.2]).

For [q.sub.i] = [[theta].sub.1], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (31)

where [T.sub.A] is the external torque applied at point A.

For [q.sub.i] = [[theta].sub.2], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (32)

where [T.sub.D] is the external torque applied at point D.

For [q.sub.i] = [L.sub.1], (5) leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (33)

where [F.sub.L] is the external force applied at the tether.

Equations (28) and (30) to (33) are the nonlinear dynamic equations of motion for the tethered satellite system based on the coordinate system ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]). Note that the term related to [[??].sub.C] only appears in (28) but not in (30) to (33). This property is not shown in the equations of motion based on point A.

3. SDRE Control with a Specified Degree of Stability

A nonlinear and observable system is written as

[??](t) = f(x) + g(x)u(t), x(Q) = [x.sub.0], (34)

where x(t) and u(t) are the state and input vectors, respectively.

To perform a regulation problem, an infinite-time performance index is written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (35)

where Q(t) [greater than or equal to] 0 and R(t) > 0 are state-dependent weighting matrices.

To apply the LQR, (34) can be factorized in a linear-like structure as

[??](t) = A(x) x(t) + B(x)u(t), (36)

where f(x) = A(x)x(t) and g(x) = B(x). There is an infinite number of forms to factorize f(x) as A(x)x(t), so this property increases the flexibility of the SDRE.

Solving (35) and (36) by following the formulation of the LQR leads to

u(t) = -K(x) x(t) = -[R.sup.-1] (x) [B.sup.T](x) P(x) x(t), (37)

where P(x) is the symmetric and positive-definite solution of the algebraic state-dependent Riccati equation as

P(x) A(x) + [A.sup.T] (x) P(x) -P(x) B(x) [R.sup.-1] (x) [B.sup.T] (x) P(x) + Q(x) = 0. (38)

To ensure the closed-loop system with a specified degree of stability, the performance index is rewritten as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (39)

where [alpha] is a positive parameter. Note that the pair [A + [alpha]I, B] is completely stabilizable, where I is an identity matrix.

Define new state variables and controls as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (40)

Using (40), the system equations shown in (36) can be modified as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (41)

Note that the initial condition is the same as [??](0) = x(0) = [x.sub.0]. Also, the performance index is transformed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (42)

Solving (41) and (42) by following the formulation of the LQR leads to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (43)

where [bar.P](x) is the symmetric and positive-definite solution of the algebraic state-dependent Riccati equation as

[bar.P](x) [A (x) + [alpha]I] + [[A(x) + [alpha]I].sup.T] [bar.P](x) -[bar.P](x) B(x) [R.sup.-1](x) [B.sup.T](x) [bar.P](x) + Q(x) = 0. (44)

Substituting (43) into (41) leads to the optimal closed-loop system as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (45)

Roughly speaking, the system seemingly has eigenvalues with real parts less than -[alpha], so the system is treated with a degree of stability of at least [alpha].

Based on (43), the optimal control can be transformed to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (46)

Note that the original control (see (37)) has the same structure as the optimal control (see (46)) of the modified system.

4. Nonlinear Dynamic Analysis

This section presents the nonlinear dynamic analysis without applying any control theory.

4.1. Free Libration. To demonstrate the consistence of the proposed formulations, the numerical simulations of free libration are performed. The initial configuration is that three satellites and the mass center of the Earth are collinear, and satellite D is at the middle point between satellites A and B. Figure 3 shows the initial configuration, where the distance from the Earth center to satellite A is 42,522 km and the length of the tether is 10,000 km. The three satellites have the same initial orbital rate [[??].sub.A] = 7.2923 x [10.sup.-5] rad/s. Hence, the initial velocities of these satellites can be determined and are expressed as [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], shown in Figure 3. The masses of satellites are identical as 10,000 kg, and the mass of satellite is 100 kg. Based on the aforementioned initial state, the initial conditions are given as [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] for the formulations based on point A. The initial conditions are [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] for the formulations based on point C.

Since both formulations represent the same system and the same initial state is applied, the numerical simulations must be identical. Figure 4 shows distances [R.sub.A] and [R.sub.C] as functions to time. The results show that the distances periodically change and one period is around 2.389 x [10.sup.4] seconds. Besides, the mean values of [R.sub.A] and [R.sub.C] are 7.7083 x [10.sup.4] and 8.1791 x [10.sup.4] km, respectively. Also, their amplitude is 3.4561 x [10.sup.4] and 3.4269 x [10.sup.4] km, respectively. Figure 5 shows the true anomalies [[psi].sub.A] and [[psi].sub.C] as functions of time. The results show that they do not linearly change with time. Figure 6 shows the libration angles [[phi].sub.1] and [[phi].sub.2] and the tether length [L.sub.1] as functions of time. Since the total length of the tether is 10,000 km, the tether length [L.sub.1] must be less than the total length. This figure shows that length [L.sub.1] changes between 0 and 10,000 km. Besides, while length [L.sub.1] reaches maximum value 10,000 km, the libration angle 02 will have the change of 360[degrees]. Similarly, while length [L.sub.1] reaches the minimum value 0, the libration angle [[phi].sub.1] also has the change of 360[degrees]. Figure 7 shows the libration angles [[theta].sub.1] and [[theta].sub.2] and tether length [L.sub.1] as functions of time, and this figure is similar to Figure 6.

4.2. Free Librations with a Prescribed Function [L.sub.1](t). One considers that satellite D is a climber, which will move along the tether based on a prescribed function. The climber is initially located at [L.sub.10] = 0.1 km. It ascends while accelerating until reaching [L.sub.11] = 0.1L. It cruises at a constant velocity [v.sub.c] until [L.sub.12] = [L.sub.1f] - 0.1L, where [L.sub.1f] is the final location of the climber. During the acceleration and deceleration periods, the velocity varies as a sinusoidal function. Hence, the motion of the climber is described by the following piecewise function as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (47)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (48)

In this section, the initial conditions are the same as those in Section 4.1, except for the fact that [L.sub.1] = 5000 km and [L.sub.1] = 0 are removed. They are replaced by the prescribed function shown in (47) to (48), where the constant velocity [v.sub.c] and the final position [L.sub.1f] are specified as 0.5 km/s and 4,628 km, respectively. Since both formulations based on points A and C represent the same system, and the same initial state is considered, the numerical simulations by both formulations must be identical. Figure 8 shows distances [R.sub.A] and [R.sub.C] as functions to time. The results show that the distances periodically change, and one period is around 2.372 x [10.sup.4] seconds. Besides, the mean values of [R.sub.A] and [R.sub.C] are 7.6706 x [10.sup.4] and 8.1393 x [10.sup.4] km, respectively. Also, their amplitude is 3.4184 x [10.sup.4] and 3.3896 x [10.sup.4] km, respectively. Figure 9 shows the true anomalies [[psi].sub.A] and [[psi].sub.C] as functions of time. The results show that they do not linearly change with time. To compare the results with those in Section 4.1, all the values are similar, so the prescribed function does not affect [R.sub.A] and [R.sub.C] very much. Figure 10 shows libration angles [[phi].sub.1] and [[phi].sub.2] and tether length [L.sub.1] as functions of time. The results show that this prescribed function [L.sub.1](t) yields the librations angles varying around between [+ or -] 50[degrees]. Figure 11 shows libration angles [[theta].sub.1] and [[theta].sub.2] as well as tether length [L.sub.1] as functions of time, and the results are similar to those shown in Figure 10.

4.3. Free Libration with the Prescribed Function [L.sub.1](t) on a Circular Orbit with a Constant Angular Rate. In this subsection, one considers that the tethered satellite system moves on a circular orbit, so the radius and the true anomaly become known parameters, and the system is governed by two ordinary differential equations. The tether length function [L.sub.1](t) is the same as (47).

For the formulations based on point A, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The initial conditions are set as zeros. Figure 12 shows the libration angles [[phi].sub.1](t) and [[phi].sub.2](t) as functions of time based on the prescribed function [L.sub.1](t) on a circular orbit. The results show that there is a transient response during the first 1.5 x [10.sup.5] seconds, and the steady state shows that the libration angles oscillate around zero degrees with the amplitude of 0.01 degrees. For the formulations based on point [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The initial conditions are also set as zeros. Note that the initial state is different from the former case, because the circular orbit is assigned to different points. Figure 13 shows the libration angles [[phi].sub.1](t) and [[phi].sub.2](t) as functions of time based on the prescribed function [L.sub.1](t) on a circular orbit. The results show that there is a transient response during the first 1.5 x [10.sup.5] seconds, and the steady state shows that the libration angles oscillate around zero degrees with the amplitude of 0.01 degrees.

The following simulations are similar to the former ones, except for the fact that the initial conditions are replaced by [[phi].sub.1] = [[phi].sub.2] = [[theta].sub.1] = [[theta].sub.2] = 0.1 degrees. Based on the formulations based on point A, Figure 14 shows the libration angles [[phi].sub.1](t) and [[phi].sub.2](t) as functions of time based on the prescribed function [L.sub.1](t) on a circular orbit. The results show that there is a transient response during the first 1.5 x [10.sup.5] seconds, and the steady state shows that the libration angles oscillate around zero degrees with the amplitude of 0.1 degrees. Based on the formulations based on point C, Figure 15 shows the libration angles [[phi].sub.1](t) and [[phi].sub.2](t) as functions of time based on the prescribed function [L.sub.1](t) on a circular orbit. The results show that there is a transient response during the first 1.5 x [10.sup.5] seconds, and the steady state shows that the libration angles oscillate around zero degrees with the amplitude of0.1 degrees.

5. Libration Suppression

This section demonstrates the libration suppression by applying the SDRE control with a prescribed degree of stability, which is presented in Section 3.

5.1. Libration Control for the Equations of Motion Based on Point A. One considers that the tethered satellite system moves on a circular orbit with a constant rotating rate and satellite D is fixedly located at the middle point of the tether. Thus, there are only two ordinary differential equations taken into account, and each equation is considered with a control torque, which is obtained based on the SDRE control method.

Based on the formulations of the SDRE control, the two equations of motion are factorized as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (49)

where [A.sup.([phi]), [B.sup.([phi]), and [C.sup.([phi]) are 2 x 2 matrices, D is a 2 x 2 identity matrix, [phi] = [[[[phi].sub.1] [phi].sub.2]].sup.T], and u = [[[T.sub.A] [T.sub.D]].sup.T]. Besides, the components of the matrices and the vector are written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (50)

The simulations are performed based on the following parameters: [R.sub.A] = 42,522 km, [[??].sub.A] = 7.2923 x [10.sup.-5] rad/s, and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Besides, the initial conditions are given as [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. For simplicity, the weighting matrices in the SDRE control are selected as Q = diag(1, 1, 0,0) and R = diag(1, 1). In order to demonstrate the performance of the prescribed degree of stability in the SDRE control, there are two cases tested as follows.

Case 1 ([alpha] = 0). Figure 16 shows the libration angles as functions of time, and the results show that the SDRE control can suppress the librations around within 10 orbits. Although the time responses seemingly have high-frequency oscillations in the transient states (see Figure 16(a)), the oscillation period is around 1,930 seconds (see Figure 16(b)). Figure 17 shows the control torques as function of time, and the results show that the maximum control torque is around 0.16 MN-m.

Case 2 ([alpha] = [10.sup.-6]). Seemingly, the convergence rate for [alpha] = 0 is little slow, so this test case is based on [alpha] = [10.sup.-6]. Figure 18 shows the libration angles as functions of time, and the results show that the SDRE control can suppress the librations around within 8 orbits, which indicates that the proposed SDRE control with a prescribed degree of stability can effectively reduce the converge time. This figure also shows a high-frequency oscillation in the transient state, and the oscillation period is the same as 1,930 seconds. Figure 19 shows the control torques as function of time, and the results show that the maximum control torque is around 0.28 MN-m, which is larger than those presented in previous case. This is the penalty to reduce the convergence time.

5.2. Libration Control for the Equations of Motion Based on Point C. The same orbit and orbital rate are considered in this subsection, and satellite D is also fixedly located at the middle point of the tether. Thus, there are still only two ordinary differential equations taken into account, and each equation is considered with a control torque, which will be determined by applying the SDRE control method.

Based on the formulations of the SDRE control, the two equations of motion are factorized as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (51)

where [A.sup.([theta])], [B.sup.([theta])], and [C.sup.([theta])] are 2 x 2 matrices and [theta] refers to [theta] = [[[[theta].sub.1] [[theta].sub.2]].sup.T]. Besides, the components of the matrices and the vector are written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (52)

In order to compare the simulations with those in previous subsection, the initial conditions are the same, so the initial libration angles are given as [[theta].sub.1](0) = 0.9649[degrees] and [[theta].sub.2](0) = -1.0351[degrees]. Based on the initial condition, the orbital radius is given as [R.sub.C] = 47,521km, and the orbital rate is the same as [[psi].sub.C] = 7.2923 x [10.sup.-5] rad/s. Besides, the parameters are given as [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. In the SDRE control, the weighting matrices are the same as Q = diag(1, 1,0,0) and R = diag(1, 1). Note that the weighting matrices have different meaning due to different state variables. Thus, although most of parameters are the same as those in previous subsection, the simulation results will be different. Similarly, two cases are tested to demonstrate the prescribed degree of stability.

Case 3 ([alpha] = 0). Figure 20 shows the libration angles as functions of time, and the results show that the SDRE control can suppress the librations around within 10 orbits. Although the time responses seemingly have high-frequency oscillations in the transient states (see Figure 20(a)), the oscillation period is around 3,954.8 seconds (see Figure 20(b)), which is almost twice time in Figure 16(b). Figure 21 shows the control torques as function of time, and the results show that the maximum control torque is around 0.35 MN-m, which is also larger than the torque in Figure 16.

Case 4 ([alpha] = [10.sup.-6]). Seemingly, the convergence rate for [alpha] = 0 is little slow, so this test case is based on [alpha] = [10.sup.-6]. Figure 22 shows the libration angles as functions of time, and the results show that the SDRE control can suppress the librations around within 8 orbits, which indicates that the proposed SDRE control with a prescribed degree of stability can effectively reduce the converge time. This figure also shows a high-frequency oscillation in the transient state, and the oscillation period is the same as 3,954.8 seconds. Figure 23 shows the control torques as function of time, and the results show that the profile is similar to that in Case 1, and the difference is that the torque converges faster.

6. Conclusions

This paper investigates the nonlinear dynamic analysis and control of a tethered satellite system. First, the system is modeled by a two-piece dumbbell. There are three satellites in the system, and each satellite is assumed to be a point mass. A constant-length massless tether connects the three satellites. Based on the model, two sets of equations of motion are derived, which are based on the first satellite and the mass center of the system. Secondly, the SDRE with a prescribed degree of stability is proposed, which can reduce the convergence time of a nonlinear system. Thirdly, the nonlinear dynamic analysis is performed by numerical simulations. The results show that the librations of the system affected the motion of the middle satellite. Also, its motion affects the orbital radius and orbital rate. Finally, the proposed SDRE approach is developed and is applied to the system, and the numerical examples demonstrate the efficiency of the proposed approach in the system.

http://dx.doi.org/10.1155/2016/3610507

Competing Interests

The author declares that they have no competing interests.

Acknowledgments

The study was sponsored by a Grant, MOST 103-2221-E-011-028, from the Ministry of Science and Technology, Taiwan.

References

[1] K. E. Tsiolkovsky, Speculations Between Earth and Sky, Isd-vo AN-SSSR, 1895, (reprinted in 1959).

[2] G. Colombo, E. M. Goposchkin, M. D. Grossi, and G. C. Weiffenbach, "Shuttle-borne skyhook: a new tool for low-orbital-altitude research," Contract NAS8-02138, Smithsonian Institution Astrophysical Observatory Report, Cambridge, Mass, USA, 1974.

[3] M. L. Cosmo and E. C. Lorenzini, Tethers in Space Handbook, NASA Marshall Space Flight Center, Huntsville, Ala, USA, 3rd edition, 1997

[4] Y. He, B. Liang, and W. Xu, "Study on the stability of Tethered Satellite System," Acta Astronautica, vol. 68, no. 11-12, pp. 1964-1972, 2011.

[5] K. Nakanishi, H. Kojima, and T. Watanabe, "Trajectories of in-plane periodic solutions of tethered satellite system projected on van der Pol planes," Acta Astronautica, vol. 68, no. 7-8, pp. 1024-1030, 2011.

[6] H. Kojima, Y. Furukawa, and P. M. Trivailo, "Experimental study on delayed feedback control for libration of tethered satellite system," Journal of Guidance, Control, and Dynamics, vol. 35, no. 3, pp. 998-1002, 2012.

[7] W. Zhang, F. B. Gao, and M. H. Yao, "Periodic solutions and stability of a tethered satellite system," Mechanics Research Communications, vol. 44, pp. 24-29, 2012.

[8] H. Hu, H. Wen, and D. P. Jin, "Dynamic control and ground-based experiments of a tethered satellite system," in IUTAM Symposium on Nonlinear Dynamics for Advanced Technologies and Engineering Design: Proceedings of the IUTAM Symposium on Nonlinear Dynamics for Advanced Technologies and Engineering Design, held Aberdeen, UK, 27-30 July 2010, vol. 32 of IUTAM Bookseries, pp. 375-386, Springer, Amsterdam, The Netherlands, 2013.

[9] W. Jung, A. P. Mazzoleni, and J. Chung, "Dynamic analysis of a tethered satellite system with a moving mass," Nonlinear Dynamics, vol. 75, no. 1-2, pp. 267-281, 2014.

[10] W. Jung, A. P. Mazzoleni, and J. Chung, "Nonlinear dynamic analysis of a three-body tethered satellite system with deployment/retrieval," Nonlinear Dynamics, vol. 82, no. 3, pp. 11271144, 2015.

[11] A. A. T. Hong and R. Varatharajoo, "Dynamics of a flexible tethered satellite system utilising various materials for coplanar and non-coplanar models," Advances in Space Research, vol. 56, no. 4, pp. 648-663, 2015.

[12] A. A. T. Hong, R. Varatharajoo, and F. Romli, "Analytic solutions for Tethered Satellite System (TSS) subject to internal and external torques undergoing a spin-up and spin-down maneuvers," Aircraft Engineering and Aerospace Technology, vol. 87, no. 6, pp. 514-520, 2015.

[13] A. K. Misra and V. J. Modi, "Dynamics and control of tether connected two-body systems--a brief review," in Proceedings of the Paris International Astronautical Federation Congress, vol. 1, Paris, France, 1982.

[14] M. P. Cartmell and D. J. McKenzie, "A review of space tether research," Progress in Aerospace Sciences, vol. 44, no. 1, pp. 1-21, 2008.

[15] P. Williams, "A review of space tether technology," Recent Patents on Space Technology, vol. 2, no. 1, pp. 22-36, 2012.

[16] Y. Chen, R. Huang, X. Ren, L. He, and Y. He, "History of the tether concept and tether missions: a review," ISRN Astronomy and Astrophysics, vol. 2013, Article ID 502973, 7 pages, 2013.

[17] Y. Chen, R. Huang, L. He, X. Ren, and B. Zheng, "Dynamical modelling and control of space tethers: a review of space tether research," Nonlinear Dynamics, vol. 77, no. 4, pp. 1077-1099, 2014.

[18] J. D. Pearson, "Approximation methods in optimal control I. Sub-optimal control," Journal of Electronics and Control, vol. 13, no. 5, pp. 453-469, 1962.

[19] J. R. Coultier, C. N. Souza, and C. P. Mracek, "Nonlinear regulation and nonlinear H, control via the state-dependent Riccati equation technique. Part 1: theory; Part 2: examples," in Proceedings of the International Conference on Nonlinear Problems in Aviation and Aerospace, pp. 117-141, Daytona Beach, Fla, USA, May 1996.

[20] X. Zhao and F. Deng, "Solution of the HJI equations for nonlinear HM control design by state-dependent Riccati equations approach," Journal of Systems Engineering and Electronics, vol. 22, no. 4, pp. 654-660, 2011.

[21] S. R. Nekoo, "Nonlinear closed loop optimal control: a modified state-dependent Riccati equation," ISA Transactions, vol. 52, no. 2, pp. 285-290, 2013.

[22] I. Chang and J. Bentsman, "Constrained discrete-time state-dependent Riccati equation technique: a model predictive control approach," in Proceedings of the 52nd IEEE Conference on Decision and Control (CDC 2013), pp. 5125-5130, Firenze, Italy, December 2013.

[23] N. Babaei and M. U. Salamci, "State dependent riccati equation based model reference adaptive control design for nonlinear systems," in Proceedings of the 24th International Conference on Information, Communication and Automation Technologies (ICAT '13), pp. 1-8, IEEE, Sarajevo, Bosnia and Herzegovina, November 2013.

[24] S. R. Nekoo and B. Geranmehr, "Nonlinear observer-based optimal control using the state-dependent Riccati equation for a class of non-affine control systems," Control Engineering and Applied Informatics, vol. 16, no. 2, pp. 5-13, 2014.

[25] X. Wang, J. Long, W. Sun, and W. Lian, "The generalized state dependent riccati equation control of continuous time nonlinear systems," in Proceedings of the International Mechanical Engineering Congress and Exposition (ASME '14), American Society of Mechanical Engineers, 2014.

[26] M. H. Korayem and S. R. Nekoo, "Nonlinear optimal control via finite time horizon state-dependent Riccati equation," in Proceedings of the 2nd RSI/ISM International Conference on Robotics and Mechatronics (ICRoM2014), pp. 878-883, Tehran, Iran, October 2014.

[27] M. H. Korayem and S. R. Nekoo, "Finite-time state-dependent Riccati equation for time-varying nonaffine systems: rigid and flexible joint manipulator control," ISA Transactions, vol. 54, pp. 125-144, 2015.

[28] J. R. Cloutier, "State-dependent Riccati equation techniques: an overview," in Proceedings of the American Control Conference, vol. 2, pp. 932-936, Albuquerque, NM, USA, June 1997.

[29] J. R. Cloutier and D. T. Stansbery, "The capabilities and art of State-dependent Riccati equation-based design," in Proceedings of the American Control Conference, vol. 1, pp. 86-91, May 2002.

[30] T. Cimen, "State-dependent Riccati equation (SDRE) control: a survey," in Proceedings of the 17th World Congress of the International Federation of Automatic Control (IFAC '08), pp. 611, Seoul, Republic of Korea, July 2008.

[31] T. Cimen, "Systematic and effective design of nonlinear feedback controllers via the state-dependent Riccati equation (SDRE) method," Annual Reviews in Control, vol. 34, no. 1, pp. 32-51, 2010.

[32] T. Cimen, "Survey of state-dependent Riccati equation in nonlinear optimal feedback control synthesis," Journal of Guidance, Control, and Dynamics, vol. 35, no. 4, pp. 1025-1047, 2012.

Yong-Lin Kuo

Graduate Institute of Automation and Control, National Taiwan University of Science and Technology, No. 43, Sec. 4, Keelung Road, Taipei City 10607, Taiwan

Correspondence should be addressed to Yong-Lin Kuo; yLkuo@yahoo.com

Received 1 April 2016; Accepted 25 July 2016

Academic Editor: Anton V. Doroshin

Caption: Figure 1: A tethered satellite system.

Caption: Figure 2: Simplified tethered satellite systems based on (a) point A and (b) point C.

Caption: Figure 3: Initial configuration for free libration.

Caption: Figure 4: The distances from Earth center as functions of time for free libration.

Caption: Figure 5: True anomaly as functions of time for free libration.

Caption: Figure 6: Libration angles and tether length [L.sub.1] by the formulations based on point A.

Caption: Figure 7: Libration angles and tether length [L.sub.1] by the formulations based on point C.

Caption: Figure 8: The distances from Earth center as function of time for a prescribed function [L.sub.1](t).

Caption: Figure 9: True anomaly as function of time for a prescribed function [L.sub.1](t).

Caption: Figure 10: Libration angles and tether length [L.sub.1] by the formulations based on point A for a prescribed function [L.sub.1](t).

Caption: Figure 11: Libration angles and tether length [L.sub.1] by the formulations based on point C with a prescribed function [L.sub.1](t).

Caption: Figure 12: Libration angles [[phi].sub.1](t) and [[phi].sub.2](t) as functions of time based on the prescribed function [L.sub.1](t) on a circular orbit with zero initial conditions.

Caption: Figure 13: Libration angles [[theta].sub.1](t) and [[theta].sub.2](t) as functions oftime based on the prescribed function [L.sub.1](t) on a circular orbit with zero initial conditions.

Caption: Figure 14: Libration angles [[phi].sub.1](t) and [[phi].sub.2](t) as functions of time based on the prescribed function [L.sub.1](t) on a circular orbit.

Caption: Figure 15: Libration angles [[theta].sub.1](t) and [[theta].sub.2](t) as functions oftime based on the prescribed function [L.sub.1] (t) on a circular orbit.

Caption: Figure 16: Libration angles from the formulations based on point A by applying the SDRE control with [alpha] = 0 for (a) the first 10 orbits and (b) the first 0.1 orbits.

Caption: Figure 17: Control torques from the formulations based on point A by applying the SDRE control with [alpha] = 0 for (a) the first 10 orbits and (b) the first 0.1 orbits.

Caption: Figure 18: Libration angles from the formulations based on point A by applying the SDRE control with [alpha] = [10.sup.-6] for (a) the first 10 orbits and (b) the first 0.1 orbits.

Caption: Figure 19: Control torques from the formulations based on point A by applying the SDRE control with [alpha] = [10.sup.-6] for (a) the first 10 orbits and (b) the first 0.1 orbits.

Caption: Figure 20: Libration angles from the formulations based on point C by applying the SDRE control with [alpha] = 0 for (a) the first 10 orbits and (b) the first 0.1 orbits.

Caption: Figure 21: Control torques from the formulations based on point C by applying the SDRE control with [alpha] = 0 for (a) the first 10 orbits and (b) the first 0.1 orbits.

Caption: Figure 22: Libration angles from the formulations based on point C by applying the SDRE control with [alpha] = [10.sup.-6] for (a) the first 10 orbits and (b) the first 0.1 orbits.

Caption: Figure 23: Control torques from the formulations based on point C by applying the SDRE control with [alpha] = [10.sup.-6] for (a) the first 10 orbits and (b) the first 0.1 orbits.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Kuo, Yong-Lin |

Publication: | International Journal of Aerospace Engineering |

Date: | Jan 1, 2016 |

Words: | 7565 |

Previous Article: | Preliminary design of a small unmanned battery powered tailsitter. |

Next Article: | Aerodynamic optimization based on continuous adjoint method for a flexible wing. |

Topics: |