Planning Allocation for GTO-GEO Transfer Spacecraft with Triple Orthogonal Gimbaled Thruster Boom

: This paper proposes an electric propulsion platform based on a triple orthogonal gimbaled thruster boom to realize the GTO-GEO transfer process. The adjustment mechanism of the gimbaled thruster boom signiﬁcantly improves the range of thrust vector variation enhances the efﬁciency of thrust vector adjustment, and reduces the spacecraft burn-up. Additionally, to achieve the application performance, a planning allocation method based on the model prediction algorithm is proposed and veriﬁed through numerical simulation.


Introduction
The geostationary orbit (GEO) is a special type where spacecraft have the same Earth's rotation period around its axis [1].The GEO is significant in promoting the development of the economy and space technology.Currently, there are several ways to deliver satellites into GEO orbit.Some launch vehicles can deliver payloads directly to GEO, but this requires the rocket to carry a huge amount of fuel, which also imposes limitations on the size of the payload.Therefore, for most missions, the launch vehicle launches the satellite into a Geostationary Transfer Orbit (GTO) and separates from the satellite.The orbiting process from GTO-GEO is then accomplished by an onboard chemical propulsion system [2].However, it requires GEO spacecraft to carry a significant amount of fuel, which limits the size and mass of the payload [3].GEO spacecraft has recently developed towards a high load-bearing ratio and long life.Therefore, many scholars have extensively researched electric propulsion technology to replace traditional chemical propulsion technology to realize the orbit transfer task of GEO orbit spacecraft [4].
During space propulsion, the energy consumption of various propulsion systems depends heavily on the type of mission.GEO is a typical high orbit, and for this type of space launch mission, achieving fast orbit entry through high thrust is not optimal.The real concern is the mass of the load that can be carried by the satellite platform, which is why electric propulsion systems have been applied and validated during GTO-GEO transfers in recent years.Compared with the traditional chemical propulsion system, the electric propulsion system has the advantages of high specific impulse and adjustable thrust.The launch mass of the spacecraft can be reduced by about 50% under the condition of carrying the same payload mass.The dry mass ratio of the propellant to the spacecraft can be reduced from about 1.5 to about 0.3, which effectively reduces the cost and significantly improves the carrying capacity of the spacecraft platform and reduces the fuel consumption of GEO spacecraft during orbital transfer [5].
On the other hand, despite the low thrust of the electric propulsion system, it can still provide high economic efficiency during high precision attitude and orbit control of GEO spacecraft [6][7][8].However, the current layout of the electric propulsion system commonly used in GEO spacecraft is the traditional quadrilateral layout in a fixed direction [9].Since this layout does not allow for real-time adjustment of thrust direction during reorbiting, there is an inevitable problem of thrust loss.This will consume a large amount of fuel and reduce the propulsion efficiency.Meanwhile, the change of thrust vector is entirely dependent on the attitude adjustment of the satellite platform, which intensifies the coupling of attitude control and orbit control of the satellite platform and increases the unnecessary fuel consumption for attitude adjustment [10,11].
To achieve real-time adjustment of thrust direction during GTO-GEO transfer and further improve the comprehensive performance of the propulsion platform, it is necessary to change the installation of thrusters and propulsion platforms so that the thrust angle can be flexibly adjusted according to the mission requirements.With the characteristics of a high-precision controlled drive, the space extension arm has been widely used in various space in-orbit missions.It is the core maneuvering equipment for spacecraft inorbit service [12].Thanks to its high flexibility, the space robotic arm can assist or replace astronauts in the harsh space environment to complete various in-orbit operational tasks, significantly improving the safety and effectiveness of space operations and applications.And with the increasing complexity of space science missions and the development of related technologies to promote the configuration of space robotic arms from a single arm to multi-arm development [13].For instance, Space Shuttle Columbia was equipped with a six-degree-of-freedom single-arm configuration [14], while the ISS used a seven-degree-offreedom redundant single-arm configuration [15,16].ISS and robotic astronaut in a doublearm configuration, ISS Japanese arm in a tandem configuration with macro and microrobotic arms [17].In the future, multi-arm configuration space robots may be developed to accomplish more complex operations through multi-arm cooperative movements.
In recent years, teams have applied the gimbaled thruster boom to the thrust vectoring process of electric propulsion platforms [18].For example, Boeing proposed the 702SP spacecraft platform [19], equipped with four XIPS-25 xenon ion thrusters with thrust vectoring mechanisms.ESA proposed the Electra spacecraft platform [20], which uses two sets of three electric thrusters with a space robotic arm, called gimbaled thruster boom, to perform the orbit holding.The movement of the gimbaled thruster boom can significantly improve the coverage of the thrust vector and enhance the control performance of the spacecraft platform [21].From the above study and in-orbit validation, it is clear that the GTO-GEO transfer strategy based on the gimbaled thruster boom has some engineering realizability.It also provides a more efficient transfer scheme for future GEO orbiting satellites with significant engineering application value.However, most current propulsion schemes employ multi-degree-of-freedom space extension arms, which increases the mechanical design difficulty of the propulsion platform and places higher requirements on the overall control system performance [22].At the same time, the coupling between the adjustment motion of the gimbaled thruster boom and the attitude control and orbit control of the spacecraft platform will be more severe under environmental disturbances [23].If the complex multidegree-of-freedom gimbaled thruster boom is directly applied to thrust vector regulation, there would be a trade-off between control system complexity and regulation performance.
Motivated by the above-mentioned research, Considering the complexity of mechanical structure and control system and engineering realizability, this paper proposes a triple orthogonal gimbaled thruster boom to realize the adjustment of thrust vector during GTO-GEO transfer.Within a certain range, the three rotating joints can realize the real-time change of thrust vector pointing according to the mission requirements.This fast response feature can significantly improve the efficiency of the propulsion platform.In addition, if the target angle is beyond the range of motion of the joint angle, the thrust vector adjustment can be combined with the traditional attitude adjustment method, i.e., the movement of the extension arm, to compensate for the required attitude angle change of the satellite.
The gimbaled thruster boom adjustment mechanism significantly increases the range of thrust vector variation, further improving the efficiency of thrust vector adjustment and reducing the total fuel consumption of the satellite.On the other hand, the triple orthogonal design can significantly simplify the mechanical system and the control system, resulting in better engineering realizability and higher reliability than the existing complex multidegree-of-freedom gimbaled thruster boom.At the same time, a drive planning method based on a model prediction algorithm is proposed and simulated to achieve high-precision drive control of the triple orthogonal gimbaled thruster boom.
This paper is organized as follows.Section 2 presents the mathematical preliminaries used in this study.Section 3 details the planning allocation design, while Section 4 presents the numerical simulation to demonstrate the feasibility and effectiveness of the proposed method.Finally, the paper concludes the main conclusions in Section 5.This paper proposes a triple orthogonal configuration of the gimbaled thruster boom to achieve the directional control of the thrust vector for the GTO-GEO orbit transfer process, as illustrated in Figure 1.

Mathematical Preliminaries
requirements.This fast response feature can significantly improve the efficiency of the propulsion platform.In addition, if the target angle is beyond the range of motion of the joint angle, the thrust vector adjustment can be combined with the traditional attitude adjustment method, i.e., the movement of the extension arm, to compensate for the required attitude angle change of the satellite.The gimbaled thruster boom adjustment mechanism significantly increases the range of thrust vector variation, further improving the efficiency of thrust vector adjustment and reducing the total fuel consumption of the satellite.On the other hand, the triple orthogonal design can significantly simplify the mechanical system and the control system, resulting in better engineering realizability and higher reliability than the existing complex multi-degree-of-freedom gimbaled thruster boom.At the same time, a drive planning method based on a model prediction algorithm is proposed and simulated to achieve high-precision drive control of the triple orthogonal gimbaled thruster boom.
This paper is organized as follows.Section 2 presents the mathematical preliminaries used in this study.Section 3 details the planning allocation design, while Section 4 presents the numerical simulation to demonstrate the feasibility and effectiveness of the proposed method.Finally, the paper concludes the main conclusions in Section 5.

Preparation
This paper proposes a triple orthogonal configuration of the gimbaled thruster boom to achieve the directional control of the thrust vector for the GTO-GEO orbit transfer process, as illustrated in Figure 1.The coordinate system OoXoYoZo is the spacecraft orbit coordinate system, and the spacecraft body coordinates system ObXbYbZb origin coincides with the spacecraft orbit coordinate system OoXoYoZo in this scheme.The coordinate system O0X0Y0Z0 is the base coordinate system of the gimbaled thruster boom, and the relative position vector between its origin and the origin of the spacecraft body coordinate system ObXbYbZb is l0.OiXiYiZi (i = 1, 2, 3) represent the three linkages coordinate systems of the gimbaled thruster boom, respectively.The thruster is mounted along the X-axis direction of the spacecraft orbit coordinate system OoXoYoZo.The joint lengths of joint 1 and joint 2 are d1 and d2, respectively, and the joint 3 linkage length is a1.The position vector of the origin of the wrist coordinate system concerning the origin of the base coordinate system.
) represent the three linkages coordinate systems of the gimbaled thruster boom, respectively.The thruster is mounted along the X-axis direction of the spacecraft orbit coordinate system O o X o Y o Z o .The joint lengths of joint 1 and joint 2 are d 1 and d 2 , respectively, and the joint 3 linkage length is a 1 .The position vector of the origin of the wrist coordinate system concerning the origin of the base coordinate system.Remark 1.To realize the re-orbit task, the projection of the required thrust vector direction in the spacecraft orbit coordinate system needs to be determined.If the angle is within the motion range, the joint can be controlled to adjust the thrust vector pointing.Otherwise, the needs to be adjusted the spacecraft attitude in conjunction with the gimbaled thruster boom.Also, it should be noted that in the proposed orthogonal configuration, the initial rotation angle of joint 2 is set to 90 deg.
The GTO-GEO orbit transfer process consists of two main phases: in-orbit plane adjustment and out-of-orbit plane adjustment.In the first phase, the spacecraft's orbital plane should be adjusted to coincide with the GEO orbital plane.In the second phase, the spacecraft is accelerated along the orbital velocity direction at the apogee, which reduces the orbital eccentricity and increases the orbital semi-major axis until the orbit becomes circular.This circularization process helps align the current orbit with the GEO orbit.
The GTO-GEO orbital transfer thrust control method used in this paper is shown in Figure 2. Remark 1.To realize the re-orbit task, the projection of the required thrust vector direction in the spacecraft orbit coordinate system needs to be determined.If the angle is within the motion range, the joint can be controlled to adjust the thrust vector pointing.Otherwise, the needs to be adjusted the spacecraft attitude in conjunction with the gimbaled thruster boom.Also, it should be noted that in the proposed orthogonal configuration, the initial rotation angle of joint 2 is set to 90 deg.
The GTO-GEO orbit transfer process consists of two main phases: in-orbit plane adjustment and out-of-orbit plane adjustment.In the first phase, the spacecraft's orbital plane should be adjusted to coincide with the GEO orbital plane.In the second phase, the spacecraft is accelerated along the orbital velocity direction at the apogee, which reduces the orbital eccentricity and increases the orbital semi-major axis until the orbit becomes circular.This circularization process helps align the current orbit with the GEO orbit.
The GTO-GEO orbital transfer thrust control method used in this paper is shown in Figure 2.During the thruster start-up process, the thrust direction is divided into two directional components: the orbital and the out-of-orbit planes.The thrust control angle α is defined as the angle between the projection of the ions vector in the orbit plane and the direction of the spacecraft geocentric vector in the orbit coordinate system, with the ions vector pointing in the radial direction being positive.The thrust control angle β defined as the angle between the ions vector and the orbit plane, and the ions vector pointing to the direction of the spacecraft angular momentum is positive.To facilitate the study, the coupling between the inside and outside of the orbit plane is ignored in this paper, and the analytical solution of the thrust control angle is defined as follows: During the thruster start-up process, the thrust direction is divided into two directional components: the orbital and the out-of-orbit planes.The thrust control angle α is defined as the angle between the projection of the ions vector in the orbit plane and the direction of the spacecraft geocentric vector in the orbit coordinate system, with the ions vector pointing in the radial direction being positive.The thrust control angle β defined as the angle between the ions vector and the orbit plane, and the ions vector pointing to the direction of the spacecraft angular momentum is positive.To facilitate the study, the coupling between the inside and outside of the orbit plane is ignored in this paper, and the analytical solution of the thrust control angle is defined as follows: where e is the spacecraft orbit eccentricity, J orbit is the spacecraft orbit true perigee angle, and u is the spacecraft orbit latitude amplitude angle.
For the scheme proposed in this paper, the mapping between the required thrust control angle of the spacecraft and the attitude angle of the spacecraft body coordinate system concerning the orbital coordinate system is given as follows: where Θ y denotes the pitch angle of the spacecraft body and coordinates system with respect to the spacecraft orbit coordinate system, and Θ z indicates the yaw angle of the spacecraft body coordinate system concerning the spacecraft orbit coordinate system.

GTO-GEO Transfer Spacecraft Kinematics and Dynamics
The GTO-GEO transfer spacecraft kinematics and dynamics are given as follows: . . . .
where r i is the position vector of the spacecraft with respect to the center of the Earth in the ECI, v i is the velocity vector of the spacecraft with respect to the center of the Earth in the geocentric inertial coordinate system, µ is the earth's gravitational parameter, a d is the perturbation acceleration of the spacecraft, m b is the mass of the spacecraft, C ib is the conversion matrix of the spacecraft body coordinate system to the ECI, f thrust,b is the projection of the thrust vector of the thrusters in the spacecraft body coordinate system, q bo,b is the projection of the vector part of the attitude quaternion of the spacecraft body coordinate system with respect to the spacecraft orbit coordinate system under the spacecraft body coordinate system, q bo,b1 is the scalar part of the attitude quaternion of the spacecraft body coordinate system with respect to the spacecraft orbit coordinate system, ω bo,b is the projection of the angular velocity vector of the spacecraft body coordinate system with respect to the spacecraft orbit coordinate system under the spacecraft body coordinate system, J body is the spacecraft inertia matrix, C bo is the conversion matrix from the spacecraft orbit coordinate system to the spacecraft body coordinate system matrix, ω oi,o is the projection of the angular velocity vector of the spacecraft orbit coordinate system with respect to the geocentric inertial coordinate system under the spacecraft orbital coordinate system, h f is the angular momentum of the spacecraft mounted flywheel, τ d is the disturbance moment of the spacecraft which takes into account the non-spherical uptake of the Earth, atmospheric drag and solar pressure, τ c is the control moment of the spacecraft, and τ thrust,b is the projection of the coupling moment generated by the thruster to the spacecraft attitude control process under the spacecraft body coordinate system.The direction of the thrust vector is controlled by three rotational joints.And θ = θ 1 θ 2 θ 3 T denotes the angle of rotation of each joint.In the spacecraft body coordinate system, the thrust vector, as well as the coupling moment, can be written as: where f is the thrust magnitude, C Wb is the transformation matrix from the spacecraft body coordinate system to the gimbaled thruster boom wrist coordinate system, C bz and C bx represent the attitude transformation matrix for an angle rotation about the j-axis (j = x, z), respectively.C W0 and r W0 represent the transformation matrix and the relative position vector from the base coordinate system to the wrist coordinate system of the gimbaled thruster boom.The combined transformation matrix T 0W can be written as follows: where s i (i = 1, 2, 3) and c i (i = 1, 2, 3) represent sin (θ i ) and cos (θ i ), respectively.
2.2.Gimbaled Thruster Boom Dynamics 2.2.1.Kinematics of the Thrust Vector-Regulated Gimbaled Thruster Boom The coordinate system L i is established at the center of mass of each connecting rod and the direction of the axes is taken to be the same as the direction of the main axis of inertia of the connecting rod to obtain the coordinate system L i concerning the base coordinate system of the gimbaled thruster boom.
where g 0L 1 (i = 1, 2, 3) represents the transformation matrix of the joint relative to the base, I 3 denotes 3 × 3 identity matrix.
In such a coordinate system, the linkage inertia matrix Ψ i has the following general form.
where m i is the mass of the linkage, I xi , I yi and I zi are the moments of inertia of the i linkage, respectively.From this, the bit shape of the base coordinate system and the wrist coordinate system is obtained as follows: To construct the kinematic rotation of the rotating joints, note that the unit vector of each joint axis of rotation ω i (i = 1, 2, 3) is: Taking the points q i (i = 1, 2, 3) on the axes respectively, thus: The motion spiral of each joint ξ i (i = 1, 2, 3) is generated by: Then, the kinematic positive solution mapping of the gimbaled thruster boom has the following form.
The exponential product formula expands the equations to obtain the kinematic positive solution mapping g 0W (θ) of the thrust vector regulating the gimbaled thruster boom.

Thrust Vector Adjustment Dynamics
To calculate the inertia matrix of the gimbaled thruster boom, it is easy to see that: where V b 0L i (i = 1, 2, 3) represents the projection of the velocity of the joint concerning the base, ξ + i (i = 1, 2, 3) is used to calculate the Jacobi matrix J i (i = 1, 2, 3) for each connecting rod. Thus: The inertia matrix of the system is M(θ) = 3 ∑ i=1 J T i ψ i J i , and the Gauche and centrifugal forces can be obtained from the inertia matrix by the following equation.
Finally, the forces acting on the gimbaled thruster boom can be expressed as: where V represents the potential energy, which can be expressed as: where h i is the height of the center of mass of the ith rod.
The transformation matrix of each linkage centre-of-mass coordinate system with respect to the base coordinate system can be obtained using the positive kinematic positive solution mapping formula.
Therefore, the height of the center of mass of each connecting rod can be obtained as follows: Taking it into the potential energy formula and differentiating it, yields: So far, the kinetic equation of the gimbaled thruster boom can be obtained as follows: where τ ARM is the driving moment of each rotating joint, M(θ) is the inertia matrix of the gimbaled thruster boom, C θ, .θ is the Gauche and centrifugal forces on the system, and θ is the potential energy equation of the system.

Controller Description
In the above section, we have determined the relationship between the thrust angle and the joint angle of the gimbal thruster boom during the GTO-GEO transfer.The next step is to perform drive control for the gimbal thruster boom.Considering the thrust angle adjustment demand of fast response and high accuracy, MPC has been applied in various research fields with its convenient modeling and high robustness [24,25].Therefore, this paper uses MPC to implement the drive control of the gimbaled thruster boom during GTO-GEO transfer.
Considering the gimbaled thruster boom dynamics as shown in Equation ( 25), writing ..
θ to the left side of the equation yields: ..
Defining the state of the system x = θ .θ T , the nonlinear state space equation of the system can be expressed as follows: . x = f(x) + g(x)u, (27) where f (x) and g(x) can be obtained from equation Equation (26).Defining x 0 as a certain state, u denotes the control input, then the above equation can be linearized around x 0 by using Taylor's Expanded Form as follows: .
Therefore, the linear state space equation can be expressed as follows. . where Considering the different initialization points, the linear state space equation can be discretized as follows.
where A k , B k and C k represent the different linearization points, respectively.Since the model predictive control is a rolling finite-time-domain optimization strategy, the planning allocation can be determined in each sampling unit.This kind of online iteration can be used to guarantee optimal allocation.Therefore, the model predictive control is adopted for planning the allocation of the gimbaled thruster boom in this paper.
Taking the linear state space equation expressed by iteration yields: Therefore, the planning allocation of the gimbaled thruster boom can be transformed into tracking optimization based on model predictive control as follows: where the function l e (x(j), u(j)) represents each tracking step's objective, the set of linear time-varying systems Equation (30), and the sets U X denote the set of constraints on the inputs and states.Let the tracking problem have an equilibrium point, set as x s , and use x s to describe the stability conveniently, which can be obtained from the following equation.
where X e is the robust control invariant set [24], i.e., for any ∀x ∈ X e , ∃u ∈ U, such that h(x, u) ∈ X e .The computation of the robust control invariant set can be referred to in the literature [25].According to the description above, the planning allocation of the gimbaled thruster boom can be demonstrated in Figure 3.In Section 2, it has been determined how the thrust angle changes.Therefore, we need to adjust the thrust angle by gimbaled thruster boom in this section.

Stability Analysis
To demonstrate the stability of the planning allocation of the gimbaled thruster boom, the following assumptions are made.

Stability Analysis
To demonstrate the stability of the planning allocation of the gimbaled thruster boom, the following assumptions are made.Assumption 1.The robust control invariant set X e has been obtained.Assumption 2. The nominal system x(j + 1) =h(x(j), u(j)) is strongly dissipative concerning the supply rate s(x, u), i.e., there exists a continuous storage function λ(•) and a K ∞ function α(•) such that the following equation holds for any x ∈ X and u ∈ U.
Assumption 3. Defining the set of N-steps reachable concerning the equilibrium point x s .
i.e., the set of all states that can reach x s in N steps.Suppose X N is a tight set and x s ∈ int(X N ), where int(•) denotes the set interior.
Assumption 4.There exists K ∞ class of functions γ(•), such that for any x ∈ X N , there exists a feasible solution to the optimization problem.Therefore, the following equation holds.
Lemma 1 gives the stability of the nominal system concerning the equilibrium point x s [27].
Lemma 1.Consider the problem for the nominal system and assume that Assumption 1-4 holds, then there exists a function λ(•) for all, such that the condition Equation ( 35) holds.Defining the rolling cost function as follows: Then the following control problem and optimization problem are equivalent.
Furthermore, the optimal value of the problem, set to V 0 N (•), is a Lyapunov function of the nominal closed-loop system with respect to the equilibrium point x s after solving the control optimization problem.
Here, the actual system is the original nonlinear system, set as We show that the optimal control rule obtained by solving the control optimization problem is also stable for the closed-loop system constituted by imposing it on the system when the deviation between the nominal and the actual system is bounded.Assumption 5. Let the deviation of the nominal system x(k + 1) =h(x(k),u(k)) and the actual system x(k + 1) =h(x(k),u(k)) be bounded, i.e., Theorem 1.Let Assumptions 1-5 hold and the partial derivative has an upper bound, i.e., for all x ∈ X, there is Assume that there exists a real number ε s > 0 and ρ s > 0, such that the following equation holds: where α(•) is the function mentioned in Assumption 2 and H is a constant that expands with where ρ min is defined as: And Ω ρ is the level set of the Lyapunov function The set Ω ρ max is the maximum level set within the set X N .At this point for any initial value x(0) ∈ Ω ρ max , the closed-loop system can converge to the robust control invariant set X e in a finite number of steps and remains within this set.
Proof of Theorem 1.According to Assumption 5, thus a Taylor expansion of the Lyapunov function V 0 N (x) can be obtained as: For x ∈ X e , one can find a positive real number H, such that the higher order expansion term of the above equation satisfies that: Since the initial conditions are equal, we have the following: If Equation (42) holds, there exists ρ s , such that the last term of the above equation is less than −ε s .Since ρ max ≥ ρ s , for all x(k) ∈ Ω ρ max , we have: The above equation illustrates that when x(k) − x s ≥ ρ s , the Lyapunov function V 0 N (x(k)) continues to fall and eventually, the system state x(k) will arrive within the set x − x s ≤ ρ s in a finite number of steps.
Given the definition of the set Ω ρ min , the system will always be inside the set Ω ρ min when the system satisfies x ∈ Ω ρ min .Therefore, the actual closed-loop system will reach the set X e in a finite number of steps and stay within X e all the time, and the bounded stability is proved.

Simulation Condition Configuration
To demonstrate the feasibility of the proposed GTO-GEO orbit transfer scheme based on the gimbaled thruster boom, this section will compare the orbit transfer scheme with the traditional quadrilateral layout in fixed direction.The parameters of the orbit and the gimbaled thruster boom involved in the simulation are shown in Table 1.The numerical simulation result of the GTO-GEO orbital transfer process are shown in Figures 4 and 5.

GTO-GEO Orbital Transfer Simulation Result and Discussion
The numerical simulation result of the GTO-GEO orbital transfer process are shown in Figures 4 and 5.As can be seen from Equation ( 1), the values of thrust control angle α and thrust control angle β are related to the satellite's true proximity angle, and the thrust control As can be seen from Equation ( 1), the values of thrust control angle α and thrust control angle β are related to the satellite's true proximity angle, and the thrust control angle in a single cycle is given in this paper considering that the trend of their changes is the same in each orbital cycle (Figure 6).angle in a single cycle is given in this paper considering that the trend of their changes is the same in each orbital cycle (Figure 6).Remark 2. From Equation ( 2), it can be seen that the thrust control angle α angle and β angle correspond to the negative values of pitch angle and yaw angle, respectively.Considering the motion angle of the Triple orthogonal Gimbaled Thruster Boom used in this paper, the thrust control angle requirement of 15deg,+15deg , the thrust control angle can be changed by driving joint 1 and joint 3, respectively, during the orbiting mission and then combined with the satellite attitude adjustment to complete the final thrust pointing requirement.The triple orthogonal gimbaled thruster boom can compensate for the required thrust vector pointing and save fuel consumption.At the same time, it can also significantly improve the efficiency of thrust vector adjustment due to its fast response time.Remark 2. From Equation (2), it can be seen that the thrust control angle α angle and β angle correspond to the negative values of pitch angle and yaw angle, respectively.Considering the motion angle of the Triple orthogonal Gimbaled Thruster Boom used in this paper, the thrust control angle requirement of Θ y ∈ [−90deg, 0deg), Θ z ∈ [−15deg, +15deg), the thrust control angle can be changed by driving joint 1 and joint 3, respectively, during the orbiting mission and then combined with the satellite attitude adjustment to complete the final thrust pointing requirement.The triple orthogonal gimbaled thruster boom can compensate for the required thrust vector pointing and save fuel consumption.At the same time, it can also significantly improve the efficiency of thrust vector adjustment due to its fast response time.
The change process of thrust control angle has been given in the GTO-GEO transfer process above, and we can make the thrust control angle meet the mission requirements by adjusting the spacecraft attitude and driving the gimbaled thruster boom joint motion.From the simulation result, it can be seen the thrusters are fixed to the spacecraft platform in the conventional scheme, so the full range of thrust control angle adjustment needs to be done through the attitude adjustment of the spacecraft.The proposed scheme can complete the above process by combining attitude adjustment and gimbaled thruster boom joint rotation.It saves the fuel required for the adjustment process and improves the efficiency of thrust vector adjustment, which has a high engineering application value.

Gimbaled Thruster Boom Simulation Result and Discussion
In the previous section, we have determined the variation requirement of the thrust control angle by the GTO-GEO orbit transfer process, and based on this, we have obtained the target angles of each joint of the gimbaled thruster boom in this process, which denote , as shown in Figure 11.From the simulation result, it can be seen the thrusters are fixed to the spacecraft platform in the conventional scheme, so the full range of thrust control angle adjustment needs to be done through the attitude adjustment of the spacecraft.The proposed scheme can complete the above process by combining attitude adjustment and gimbaled thruster boom joint rotation.It saves the fuel required for the adjustment process and improves the efficiency of thrust vector adjustment, which has a high engineering application value.

Gimbaled Thruster Boom Simulation Result and Discussion
In the previous section, we have determined the variation requirement of the thrust control angle by the GTO-GEO orbit transfer process, and based on this, we have obtained the target angles of each joint of the gimbaled thruster boom in this process, which denote From the simulation result, it can be seen the thrusters are fixed to the spacecraft platform in the conventional scheme, so the full range of thrust control angle adjustment needs to be done through the attitude adjustment of the spacecraft.The proposed scheme can complete the above process by combining attitude adjustment and gimbaled thruster boom joint rotation.It saves the fuel required for the adjustment process and improves the efficiency of thrust vector adjustment, which has a high engineering application value.

Gimbaled Thruster Boom Simulation Result and Discussion
In the previous section, we have determined the variation requirement of the thrust control angle by the GTO-GEO orbit transfer process, and based on this, we have obtained the target angles of each joint of the gimbaled thruster boom in this process, which denote θ = θ 1 θ 2 θ 3 T , as shown in Figure 11.As shown in Figure 11, it is joint 1 and joint 3 that need to be controlled for planning, and joint 2 remains unchanged.From Equation (2), we can see that the joint 1 turning angle corresponds to the -y axis and the joint 3 turning angle corresponds to the -z-axis, and the thrust vector during the GTO-GEO transfer is changed by adjustment.Also, we can notice that their changing trends are consistent with Figures 8 and 10.Next, in this section, the planning allocation simulation analysis of the gimbaled thruster boom is performed by the planning algorithm designed in Section 3.
The results of the planning control of the joint angle and angular velocity of the triple orthogonal gimbaled thruster boom are shown in Figures 12 and 13.As shown in Figure 11, it is joint 1 and joint 3 that need to be controlled for planning, and joint 2 remains unchanged.From Equation (2), we can see that the joint 1 turning angle corresponds to the -y axis and the joint 3 turning angle corresponds to the -z-axis, and the thrust vector during the GTO-GEO transfer is changed by adjustment.Also, we can notice that their changing are consistent with Figures 8 and 10.Next, in this section, the planning allocation simulation analysis of the gimbaled thruster boom is performed by the planning algorithm designed in Section 3.
The results of the planning control of the joint angle and angular velocity of the triple orthogonal gimbaled thruster boom are shown in Figures 12 and 13.As shown in Figure 11, it is joint 1 and joint 3 that need to be controlled for planning, and joint 2 remains unchanged.From Equation (2), we can see that the joint 1 turning angle corresponds to the -y axis and the joint 3 turning angle corresponds to the -z-axis, and the thrust vector during the GTO-GEO transfer is changed by adjustment.Also, we can notice that their changing trends are consistent with Figures 8 and 10.Next, in this section, the planning allocation simulation analysis of the gimbaled thruster boom is performed by the planning algorithm designed in Section 3.
The results of the planning control of the joint angle and angular velocity of the triple orthogonal gimbaled thruster boom are shown in Figures 12 and 13.As can be seen from Figures 12 and 13, the three-degree-of-freedom gimbaled thruster boom can accomplish the required angle change requirement by the proposed MPC planning control method.It is noteworthy that joints 1 and 3 reach the peak angle change simultaneously during the 0.2-0.25 days period, after which the planning direction becomes reversed.Since the MPC planning process is a simultaneous three-axis planning, joint 2 is affected by a certain coupling effect at that moment, but the change angle error is only 10 −5 , which is acceptable.
The results of the triple orthogonal gimbaled thruster boom drive planning error are shown in Figure 14.From Figure 14, it can be seen that with the drive planning algorithm proposed in this paper, the rotation angles of the three joints of the gimbaled thruster boom can accurately track the target curve with the angular tracking errors of 0.017 deg, 0.006 deg, As can be seen from Figures 12 and 13, the three-degree-of-freedom gimbaled thruster boom can accomplish the required angle change requirement by the proposed MPC planning control method.It is noteworthy that joints 1 and 3 reach the peak angle change simultaneously during the 0.2-0.25 days period, after which the planning direction becomes reversed.Since the MPC planning process is a simultaneous three-axis planning, 2 is affected by a certain coupling effect at that moment, but the angle error is only 10 −5 , which is acceptable.
The results of the triple orthogonal gimbaled thruster boom drive planning error are shown in Figure 14.As can be seen from Figures 12 and 13, the three-degree-of-freedom gimbaled thruster boom can accomplish the required angle change requirement by the proposed MPC planning control method.It is noteworthy that joints 1 and 3 reach the peak angle change simultaneously during the 0.2-0.25 days period, after which the planning direction becomes reversed.Since the MPC planning process is a simultaneous three-axis planning, joint 2 is affected by a certain coupling effect at that moment, but the change angle error is only 10 −5 , which is acceptable.
The results of the triple orthogonal gimbaled thruster boom drive planning error are shown in Figure 14.From Figure 14, it can be seen that with the drive planning algorithm proposed in this paper, the rotation angles of the three joints of the gimbaled thruster boom can accurately track the target curve with the angular tracking errors of 0.017 deg, 0.006 deg, From Figure 14, it can be seen that with the drive planning algorithm proposed in this paper, the rotation angles of the three joints of the gimbaled thruster boom can accurately track the target curve with the angular tracking errors of 0.017 deg, 0.006 deg, and 0.011 deg, respectively, which illustrates the effectiveness of the drive planning algorithm.Meanwhile, it is proved that the driving planning of the gimbaled thruster boom and the attitude adjustment of the satellite platform can meet the thrust vector pointing requirement during the GTO-GEO transfer process.The feasibility of the GTO-GEO transfer scheme proposed in this paper is further verified.
Based on the above simulation results, it can be seen that the conventional scheme has a fixed thruster arrangement on the spacecraft platform.This limits the thrust control angle adjustment, which requires adjusting the spacecraft's attitude to achieve the full range of thrust control angle adjustment.In contrast, the proposed thrust vector adjustment scheme based on a tri-orthogonal thrust vector adjustment gimbaled thruster boom can accomplish this through attitude adjustment and joint rotation.The fuel required in the adjustment process is reduced, and the efficiency of thrust vector adjustment is improved.Therefore, the proposed scheme has high value for engineering applications.

Conclusions
To reduce mission cycle time, decrease fuel consumption during GTO-GEO transfer, and enhance overall spacecraft efficiency, this paper proposes a GTO-GEO transfer scheme using a triple orthogonal gimbaled thruster boom for thrust vector adjustment.By rotating the joints, the spatial pointing of the thrust vector can be changed, allowing for compensation of the required attitude angle change of the spacecraft through the motion of the gimbaled thruster boom.A planning allocation method based on the model prediction algorithm is also proposed for high-precision control.The stability of the closed-loop system is rigorously proven using the Lyapunov method, and the simulation results demonstrate the feasibility and effectiveness of the proposed transfer scheme.Compared to the quadrilateral layout in a fixed direction, the proposed scheme achieves 95% thrust efficiency, approximately 22.9% reduction in fuel consumption, and 29.7% reduction in mission time.Furthermore, the proposed algorithm effectively completes gimbaled thruster boom control during the GTO-GEO transfer process and achieves a higher control performance.

Figure 1 .
Figure 1.The triple orthogonal configuration of the gimbaled thruster boom.

Figure 1 .
Figure 1.The triple orthogonal configuration of the gimbaled thruster boom.The coordinate system O o X o Y o Z o is the spacecraft orbit coordinate system, and the spacecraft body coordinates system O b X b Y b Z b origin coincides with the spacecraft orbitcoordinate system O o X o Y o Z o in this scheme.The coordinate system O 0 X 0 Y 0 Z 0 isthe base coordinate system of the gimbaled thruster boom, and the relative position vector between its origin and the origin of the spacecraft bodycoordinate system O b X b Y b Z b is l 0 .O i X i Y i Z i (i = 1, 2, 3) represent the three linkages coordinate systems of the gimbaled thruster boom, respectively.The thruster is mounted along the X-axis direction of the spacecraft orbit coordinate system O o X o Y o Z o .The joint lengths of joint 1 and joint 2 are d 1 and d 2 , respectively, and the joint 3 linkage length is a 1 .The position vector of the origin of the wrist coordinate system concerning the origin of the base coordinate system.

Assumption 1 .Assumption 2 .
The robust control invariant set e has been obtained.The nominal system ) ( (j+1)= (j), (j) x h x u is strongly dissipative concerning the supply rate ( , ) s x u , i.e., there exists a continuous storage function ()  λ and a  K function ()  α such that the following equation holds for any  x and  u .

Figure 6 .
Figure 6.Thrust control angle changes during the orbital transfer.

Figure 6 .
Figure 6.Thrust control angle changes during the orbital transfer.

Figure 8 .
Figure 8. Pitch angle change of proposed scheme.

Figure 8 .
Figure 8. Pitch angle change of proposed scheme.Figure 8. Pitch angle change of proposed scheme.

Figure 8 .
Figure 8. Pitch angle change of proposed scheme.Figure 8. Pitch angle change of proposed scheme.

Figure 9 .
Figure 9. Yaw angle change of traditional scheme.

Figure 10 .
Figure 10.Yaw angle change of proposed scheme.

Figure 10 .
Figure 10.Yaw angle change of proposed scheme.

Figure 10 .
Figure 10.Yaw angle change of proposed scheme.

Figure 11 .
Figure 11.Target angle of each joint.

Figure 11 .
Figure 11.Target angle of each joint.