Adaptive Fuzzy Modified Fixed-Time Fault-Tolerant Control on SE(3) for Coupled Spacecraft

This paper aims to solve the control problem of coupled spacecraft tracking maneuver in the case of actuator faults, inertia parametric uncertainties, and external disturbances. Firstly, the spacecraft attitude and position coupling kinematics and dynamics model are established on the Lie group SE(3), and the coupled relative motion tracking error model is derived by exponential coordinates. Then, considering the actuator faults, an adaptive fuzzy scheme is proposed to estimate the lumped disturbances in real time, and a novel modified fixed-time terminal sliding mode fault-tolerant control law is developed to deal with the actuator faults and compensate the lumped disturbances. Next, the Lyapunov method is used to prove the stability and convergence of the system. Finally, the proposed controller can achieve fast and high-precision fault-tolerant control goals, and its effectiveness and feasibility are verified by numerical simulation.


Introduction
In the context of the rapid development of space technology, new and higher requirements have been put forward for the mobility and accuracy of spacecraft. e modeling and control of the attitude and trajectory of relatively moving spacecraft has always been a hot research topic in the fields of space rendezvous and docking, spacecraft formation flying (SFF) [1,2]. Due to the strong coupling and nonlinearity of the relative motion of spacecraft's attitude and position motion, the conventional idea of dividing attitude and position motion into independent two-channel control ignores the influence of coupling between the two, although it satisfies the requirements of some space missions. However, for aerospace missions with high-precision requirements, the divide-and-conquer method will appear powerless [3]. erefore, seeking the integrated control of spacecraft attitude and position has theoretical guidance and is of great significance to engineering practice.
Due to long-term exposure to harsh space environments such as strong radiation and ultra-low temperature, the actuator will have various types of failures. erefore, the conventional control theory based on the normal operation of the actuator may be difficult to cope with the failure and may eventually cause the system to crash or fail. In addition, the spacecraft itself will also face the uncertainty of internal parameters and external disturbances, which brings huge challenges to the design of the control system. So it is particularly important to choose a suitable fault-tolerant control strategy for the aforementioned drawbacks, which also provides a strong guarantee for the long-term service operation of the spacecraft.
At present, many research results have been made on the problem of spacecraft attitude fault-tolerant control [4,5]. But for the spacecraft attitude and position coupling control system, when various faults occur in the relative attitude and position actuators at the same time, the related six-degreeof-freedom (6-DOF) fault-tolerant control algorithm design is not enough [6]. Dong et al. [7,8], studied the integrated fault-tolerant control of the spacecraft's position and attitude in the case of actuator failure based on the dual quaternion, and their numerical simulation results verified the effectiveness of the algorithm.
In recent years, the coupled modeling of rotation and translation relative motion based on different forms of rigid spacecraft has attracted widespread attention. Common spacecraft attitude and position coupling modeling and control forms mainly include dual quaternion [7,8], Lie group SE(3) [9,10], modified Rodrigues parameters (MRPs), and other forms [11,12]. Although the integrated modeling method of spacecraft attitude and position based on dual quaternion is widely used, dual quaternion also has its limitations. e model based on dual quaternion uses eight parameters to describe the three-dimensional motion, so it requires unitized constraints. Sometimes improper handling of this constraint will cause problems. Moreover, since the group function corresponding to the unit quaternion is left multiplication and right multiplication, the quaternion description rotation is not unique, which will cause ambiguity, and when it is serious, it will cause unwinding problem [3]. Moreover, describing the attitude based on MRPS is nonglobal and nonunique [13]. Compared with the traditional description method in Euclidean space, the geometric framework of Lie group SE(3) is more natural and concise, the analysis results are more realistic and credible, and the designed controller is more concise, so in recent years, it received attention gradually. Lie group SE(3) can describe the three-dimensional translation and rotation of a rigid body. Lee et al. [14], Sanyal et al. [15], and Bullo and Murray [16], conducted an in-depth study on the control of the 6-DOF motion of a spacecraft on SE (3). By using the relationship between the exponential mapping function and the logarithm mapping function of Lie group and Lie algebra, the motion spinor is transformed into the corresponding spacecraft attitude and position motion equation. On this basis, various simple controllers are designed to realize the pose tracking control target [17]. In this paper, the integrated model of spacecraft attitude and position coupling is established in the framework of Lie group SE (3), which is convenient for the design of fault-tolerant controllers in the following.
Fault-tolerant control mainly includes active fault-tolerant control and passive fault-tolerant control. Among them, considering the actuator failure belongs to the integrity design category of passive fault-tolerant control, it is also a hot research direction in the field of fault-tolerant control and has obtained rich research results [18][19][20][21][22]. Fuzzy approximation can make full use of the information ability of fuzzy logic systems; it is easier to construct and can approximate nonlinear functions with arbitrary accuracy. When the actuator fails, the uncertainty of the system increases. For the parameter uncertain system, the adaptive law can be constructed by the Lyapunov method, and the uncertain parameters in the model can be replaced by the adaptive control based on the principle of equivalence. Finally, the adaptive law is designed for the estimated parameters to make the closed-loop system stable. is mainstream adaptive control method has been widely used in the field of spacecraft control due to its simple design and easy to understand [23]. Recently, many major achievements in the engineering application of fuzzy approximation methods have been reported, such as application of adaptive fuzzy controller in industrial process [24,25]. At the same time, fuzzy control is also applied to robust faulttolerant control for fault detection and actuator faults [26][27][28][29]. In addition, the fuzzy control scheme to approximate the disturbance of the spacecraft has been successfully applied, and it is effective to combine the adaptive fuzzy controller of NFTSMC in [30,31] to reject the system uncertainty. Zhang et al. [32] applied fuzzy adaptive finite-time control to the 6-DOF SFF system and achieved success when considering the consensus control problem among the followers with signal transmission time delays.
Fixed-time control is developed on the basis of finitetime control. e difference between the two is only in the form of the sliding surface. e former can achieve fixedtime convergence without relying on the initial state, while the latter's convergence time is related to the initial state. Double-power fast terminal sliding mode control is a kind of fixed time control, which can be used to realize the fixed time stability of the system, which is more useful than the finitetime sliding mode control methods, such as terminal sliding mode (TSM) [33], fast terminal sliding mode (FTSM) [34], and nonsingular fast terminal sliding mode (NFTSM) [35]. It has faster convergence speed and better control effect. Shi et al. achieved attitude tracking control of rigid spacecraft on Lie group with fixed-time convergence [36] and global fixed time attitude tracking control for the rigid spacecraft with actuator saturation and faults [37]. Gao [3]. Mobayen et al. [39] proposed a new adaptive finite-time stabilization method based on global sliding mode to advance the steady state and transient performances of a class of chaotic flows in the presence of disturbances. Also, Mobayen and Pujol-Vázquez used robust LMI approach to deal with nonlinear feedback stabilization of continuous state-delay systems with lipschitzian nonlinearities and verified the effectiveness through experiments [40]. Jafari and Mobayen [41] combined LMI approach and second-order sliding set design for a class of uncertain nonlinear systems with disturbances.
Motivated by the facts mentioned above, this paper takes the leader-follower formation spacecraft as the research object. Firstly, a dynamic model of the relative tracking error of the spacecraft attitude and position coupling with model uncertainties, external disturbances, and actuator faults is derived on the Lie group SE(3). en, the adaptive fuzzy method is used to design the sliding mode controller to realize the fixed-time fault-tolerant control. e novelty of this paper is as follows: inspired by [3,32], the proposed model in this paper takes the actuator faults into consideration. Based on the established model, a modified double-power fast terminal sliding manifold is defined by the exponential coordinates and velocity tracking errors, and then adaptive fuzzy modified fixed-time faulttolerant control schemes is proposed, in which the adaptive fuzzy control technique is applied to reject the system lumped disturbances. Compared with finite-time stability and traditional fixed time stability, the control performance obtained in this paper has significant advantages in convergence accuracy and effectiveness. e structure of this paper is as follows: Section 2 introduces the mathematics preliminaries and the rigid body dynamics of the spacecraft on SE(3) with actuator faults. Section 3 adopts fuzzy adaptive method to design the sliding mode controller to realize the modified fixed-time faulttolerant control and uses the Lyapunov method to prove the stability of the system strictly. Section 4 verifies the effectiveness of this method through numerical simulation. Section 5 draws conclusions and summarizes this paper.

Mathematics Preliminaries
In order to facilitate the design and stability proof and analysis of the integrated attitude and position control system, the following will give some related definitions and stability theory and lemmas.

2.1.
Notations. For any column vector x � [x 1 , x 2 , . . . , x n ] T ∈ R n , define the following vector and operation: (2) ‖x‖ denotes the Euclidean norm or its induced norm.
x. (5) [·] ∧ represents an operator that maps a vector to a Lie algebra. For ςe (3), it maps a vector to an skewsymmetric matrix, that is, [·] ∨ represents the operator that maps the Lie algebra to a vector. For ςe (3), it maps the skew-symmetric matrix to a vector, which is

Relative Coupled Dynamics of Spacecraft on Lie Group
SE (3). Firstly, in order to describe the space orientation of the leader-follower spacecraft and establish the kinematics and dynamics model, we introduce three reference frames that are all orthogonal coordinate systems as shown in Figure 1. e Earth-centered inertial (ECI) reference frame with the origin at the center of the Earth is represented by I { } � x I , y I , z I , which is used to describe the absolute motion of the spacecraft relative to the Earth. e body-fixed frames of the leader spacecraft and the follower spacecraft can be expressed as Lb { } � x Lb , y Lb , z Lb and Fb { } � x Fb , y Fb , z Fb , respectively; their origin is at the center of mass of the spacecraft, and the axis coincides with the principal axis of inertia.
In nature, the configuration space of rigid body motion is SE(3), which can express translation and rotation of rigid body compactly. e special Euclidean group SE(3) is the semidirect product of the three-dimensional Euclidean space and the special orthogonal space, which can be expressed as SE(3) � R 3 ⋉SO(3). R 3 is used to describe the translation of the rigid body's center of mass; SO(3) � R ∈ R 3×3 |RR T � I, det(R) � 1 is a rotation group composed of a three-dimensional rotation matrix, which is used to represent the rotation of the rigid body around the center of mass. erefore, an element g in the Lie group SE(3) can express the configuration of the spacecraft [3]: where R ∈ SO(3) is the rotation matrix of the spacecraft from the body-fixed frame to the ECI reference frame and b ∈ R 3 is the position vector from the center of mass of the Earth to the center of mass of the spacecraft in the ECI coordinate system. e angular velocity and translational velocity of the spacecraft are defined as e above velocity vectors are defined in the body-fixed frame of the spacecraft. To describe the kinematics and dynamics equations below, it is necessary to introduce the Lie group SE(3) and its corresponding Lie algebra ςe(3) to meet the following mapping relations: e adjoint matrix of g � g(R, b) ∈ SE(3) can be expressed as e Lie algebra corresponding to ξ � v T ω T T can be expressed as e adjoint matrix of ξ ∧ � ξ ∧ (ω, v) ∈ ςe(3) can be expressed as e co-adjoint matrix of ξ ∧ � ξ ∧ (ω, v) ∈ ςe(3) can be expressed as e coupled kinematics of the spacecraft in the ECI coordinate system can be expressed as e above kinematics equation can be simplified as follows: e coupled dynamics of the spacecraft in the body-fixed frame can be expressed as where m and J are the mass and moment of inertia of the spacecraft, respectively; F g and M g are the gravity gradient force and the gravity gradient moment of the spacecraft, respectively; a J 2 is the perturbation caused by the Earth's oblateness; F c and M c are the control force and control torque of the spacecraft, respectively; and F d and M d are the unknown external force and external torque caused by radiation pressure, atmosphere drag, and other bounded uncertain disturbances, respectively. e specific forms of F g , M g , and a J 2 are as follows: where μ � 398, 600.44 km 3 s − 2 is the gravitational constant of the Earth, J 2 � 0.00108263 is the perturbation caused by the Earth's oblateness, and R e � 6378.14 km is the equatorial radius of the Earth. en, the coupling dynamics of the spacecraft can be expressed in a compact form as follows: is the unified input vector related to gravity. us, combining equations (9) and (12), the coupled kinematics and dynamics of the spacecraft can be expressed compactly as Next, based on the above equations, the coupled relative motion tracking error dynamics will be derived. Let g o be the actual pose configuration of the leader spacecraft, which the leader spacecraft can be real or virtual; g b be the actual pose configuration of the follower spacecraft.
en the actual relative pose configuration between the leader-follower spacecraft is as follows: Let g d be the desired pose configuration of the leader spacecraft, then the desired relative pose configuration between the leader-follower spacecraft is us the pose configuration tracking error is as follows: In general, the desired relative pose configuration is a constant value, and the desired relative linear velocity and angular velocity are zero, which means that the follower spacecraft and the leader spacecraft keep relatively stationary in a fixed configuration. e configuration tracking error of the follower spacecraft can be expressed by exponential coordinates as e velocity tracking error of the follower spacecraft can be expressed by exponential coordinates as en h e can be expressed on SE (3) as rough the logarithm mapping between the Lie group SE(3) and the Lie algebra ςe(3), we can get the following results: Taking equation (12) into equation (22) yields However, in actual space missions, the inertia matrix Ξ is uncertain due to fuel consumption and external disturbances, so the actual inertia matrix Ξ can be expressed as where Ξ 0 is the nominal part and ΔΞ is the uncertainty part. en the inverse of the inertia matrix can be expressed as erefore, (26) can be rewritten as where H is a known deterministic term of the system and Δd is the lumped disturbances, including uncertainties and external disturbances. en the coupling model of relative motion spacecraft can be expressed as follows:

Actuator Configuration with Faults.
In this paper, the actuator of attitude control is reaction flywheel, and the actuator of orbit control is thruster.4 reaction flywheels and 4 pairs of thrusters are used. en the control vector can be expressed as where D ∈ R 6×12 is the configuration matrix, u � [u 1 , u 2 , u 3 , . . . , u 12 ] T is the actual control vector, and u i (i � 1, 2, 3, . . . , 12) is the torque or force that each flywheel or thruster can provide. e four reaction flywheels adopt the traditional installation method of three orthogonal and one inclined Mathematical Problems in Engineering installation, and the configuration structure is shown in Figure 2. e control torque distribution matrix of the reaction flywheel is Eight thrusters are installed symmetrically at the middle point of each edge of the cube in pairs. e installation mode of thrust passing through the center of mass is adopted. e configuration structure is shown in Figure 3. e control force distribution matrix of the thruster is en the attitude and position coupling integrated control distribution matrix D is According to the cause of the faults of the actuator, the failure of the actuator can be divided into the following categories: stuck, loose, saturated, damaged, or invalid. e abovementioned faults types can be unified as follows: where u c (t) is the control command of the actuator; u is the stuck fault of the actuator with bounded value and satisfies the constraint is the actuator effectiveness matrix and satisfies the constraint 0 ≤ σ i ≤ 1. σ i � 0 denotes that the i th actuator does not supply any control output; σ i � 1 means there is no fault for the i th actuator; and 0 < σ i < 1 implies the i th actuator has partially lost its effectiveness [38].
Taking equations (31) and (35) into equation (30), the integrated dynamic equation of relative motion spacecraft considering actuator fault can be expressed as follows: where

Adaptive Fuzzy Modified Fixed-Time Fault-Tolerant Controller Design and Stability Analysis
In this part, our goal is to design a fault-tolerant controller on the relative coupled dynamics so that the configuration of the spacecraft can converge to the desired state in the presence of model uncertainties and external disturbances and actuator faults in fixed time.

Introduction of Fuzzy Approximation Technique.
e spacecraft exhibits strong nonlinearity due to the influence of lumped disturbances, which will affect the performance of the controller. So it is critical to approximate the lumped disturbances with high accuracy, and fuzzy logic system (FLS) is an effective way to realize this objective. e fuzzy approximation method can make full use of fuzzy linguistic information to approximate any nonlinear continuous function. It has a good effect in fitting nonlinear function. It can approach nonlinear continuous function with arbitrary precision. e structure and basic theory of fuzzy approximation system are given below [32].
X � (x 1 , x 2 , . . . , x n ) T ∈ R n is the input variable of the FLS, and M fuzzy rules are designed for each component of the input variable, then the whole system has n M fuzzy rules, and the specific expression of each fuzzy rule is where l i � 1, . . . , M is the number of fuzzy rules for each input variable x k , z is the output of the fuzzy system, A l k is the fuzzy set of system input variables, and B l is the fuzzy set of system output.
If the fuzzy system adopts singleton fuzzifier, centeraverage defuzzifier, and product inference engine, the output of fuzzy approximation system can be obtained as follows: where μ A l k (x k ) is the membership function corresponding to the input variable x k ; in this paper, Gauss membership function is used with the form where a l k ,x l k , and b l k are all positive real parameters with 0 < a l k ≤ 1. x l k is the abscissa corresponding to the membership function μ A l k (x k ) when the maximum value is 1. Let W � (z 1 , z 2 , . . . , z M ), then equation (38) can be rewritten as follows: where β is the basis function, which can be expressed as Based on the above introduction, the total external disturbances of the follower spacecraft can be estimated by the fuzzy approximation as where W * is the optimal weight matrix and ε is the bounded approximation error of FLS. Let W be the estimation of W * . e sliding surface S is the input variable of β(X), and then the estimated value Δd of the total external disturbances Δd of the spacecraft can be expressed as en the estimation error of the optimal weight matrix is In order to design and analyze the controller, some assumptions are given below.

Assumption 1.
e output of FLS is bounded, and the estimated value of the external total disturbances is bounded such that where d m is a positive constant.

Assumption 2.
e approximation error of FLS is bounded such that where ε m is a positive constant.
Assumption 3. e optimal weight matrix of FLS is bounded such that where W m is a positive constant.

Assumption 4.
e faults of the actuators satisfy the constraint rank(DE) � 6, and this means that the redundant actuators can still combine enough control output to complete the given goal.

Remark 1.
Because the mass, moment of inertia, fault amplitude, input variables of FLS, and external disturbance of the system are bounded, so Assumption 1 is reasonable; Assumptions 2 and 3 have the property that fuzzy approximation system can fit any nonlinear continuous function, and Assumption 4 does not consider underactuated system, so it is also reasonable.

Controller Design.
In order to achieve the control goal of modified fixed time stability, 3 sliding surface forms are proposed as follows.
Firstly, the finite-time terminal sliding mode is denoted as S � ξ e + C 1 η e + C 2 sig α η e .
From the above analysis and discussion, we can see that the modified fixed-time terminal sliding mode is a combination of the finite-time terminal sliding mode and the fixedtime terminal sliding mode under different conditions. To guarantee that the relative motion can spacecraft converge to the desired state in the expected time even in the case of actuator faults, the corresponding three kinds of adaptive fuzzy sliding mode controllers are designed as follows: (1) Corresponding to (48), adaptive fuzzy finite-time fault-tolerant controller (AF-Finite) is given as (2) Corresponding to (49), adaptive fuzzy fixed-time fault-tolerant controller (AF-Fixed) is given as where ( en the adaptive update law of the optimal weight matrix W is given by where c > 0 is an auxiliary parameter independent of control.

Stability Analysis.
In this part, we will take the stability proof of AF-finite fault-tolerant controller as an example for stability analysis, and the other two (AF-Fixed and AF-MFixed) stability analysis methods are the same as is. Some lemmas are given before the stability analysis.
Next, the Lyapunov method will be used to prove the reachability of sliding mode variables and the convergence of the system states. Proof. When equation (48) reaches the sliding mode surface S � 0, such that then we have A candidate Lyapunov function is selected as follows: Mathematical Problems in Engineering 9 Taking the derivative of V with respect to time yields: where a 1 � 2λ min (G(η e )C 1 ), a 2 � 2 (1+α/2) λ min (G(η e )C 2 ), by using Lemma 4 we know that the eigenvalues of matrix G(η e ) are all positive; in addition, C 1 , C 2 ∈ R 6×6 are both positive definite diagonal matrices, and then we have a 1 , a 2 > 0. By using Lemma 1, we can conclude that V will converge to equilibrium point in finite time T, such that So, η e can also converge to equilibrium point in finite time; according to equation (62), ξ e will converge to equilibrium point in finite time too. □ Remark 3. Since αC 2 diag(|η e | α− 1 )G(η e )ξ e is included in equation (51), when η e reaches the equilibrium point before ξ e , the control output will become infinite. e singularity can be avoided by selecting α ∈ ((1/2), 1) in this paper. Because when η e � 0, the following equation holds: Remark 4. When we choose the other two sliding surfaces in equations (49) and (50) and use the same proof method as eorem 1, we can also get that the system state will converge to the equilibrium point in fixed time. By using equation (49) and Lemma 2, η e , ξ e will converge to the equilibrium point in fixed time T f : where a 1f � 2 (1+α 1 /2) λ min (G(η e )C 1 ) and a 2f � 2 (1+α 2 /2) λ min (G(η e )C 2 ). By using equation (50) and Lemma 3, η e , ξ e will converge to the equilibrium point in modified fixed time T m : where a 1m � 2 (1+(1/2)+(1/2)α 1 +((1/2)α 1 − (1/2))sgn(|η e |− 1)/2) λ min (G(η e )C 1 ) and a 2m � 2 (1+α 2 /2) λ min (G(η e )C 2 ).
To guarantee that the sliding surface can reach S � 0 in finite time, eorem 2 is proposed.

Theorem 2.
For the relative motion spacecraft system with actuator faults in equation (36), the sliding surface S of the system can converge to a small region containing zero in finite time when using sliding surface in equation (48) and fuzzy adaptive control law in equations (51) and (54).
Proof. Another candidate Lyapunov function is selected as follows: Taking the derivative of V 1 yields Substituting the adaptive law of equation (54) into the above equation yields where Δ is defined as From Assumptions 2 and 3, we know that the following inequalities are satisfied: where χ 1 and χ 2 are defined to satisfy the following equations: (74) en equation (71) can be simplified as e above equation can be rewritten as where , by using Lemma 1, V 1 will converge to the equilibrium point in finite time.
When χ 1 > 0, that is, V 1 > (Δ ′ /χ 1 ), V 1 will converge to the region Δ 1 containing zero in finite time T 1 : When χ 2 > 0, that is, V 1 > (Δ ′ /χ 2 ) (2/1+α) , V 1 will converge to the region Δ 2 containing zero in finite time T 2 : ln According to equations (77) and (78), we can conclude that V 1 will converge to the region Δ containing zero in finite time T ′ : Since the following inequality holds the sliding surface S of the system can converge to a small region Δ S � ������������ (2/λ min (Ξ 0 ))Δ containing zero in finite time. e estimated value of the optimal weight matrix can also converge to the true value. □ Remark 5. When we choose the other two controllers in equations (52) and (53) and use the same adaptive update law in equation (54), we can also draw the conclusion that the sliding surface converges in fixed time. By using equation (52) and Lemma 2, V 1 will converge to the region Δ f containing zero in fixed time T f ′ . where e derivation process of χ 1f , χ 1f , χ 2f , and χ 2f can be referenced with eorem 2.
By using equation (53) and Lemma 3, V 1 will converge to the region Δ m containing zero in fixed time T m ′ :

Numerical Simulation Analysis
In this part, three kinds of fuzzy adaptive finite-time and fixed-time fault-tolerant control algorithms proposed in this paper are simulated to verify the effectiveness of the algorithms. Before the simulation, the input expression of the fuzzy approximation system is defined as Seven fuzzy membership functions are selected as follows: In the simulation, the mass and moment of inertia of the follower spacecraft and the leader spacecraft are chosen the same as e leader spacecraft moves on a Molniya orbit, and its initial orbital elements are given in Table 1 [3].
Assuming that the leader spacecraft moves along an ideal orbit, its orbit is generated by offline calculation. At the initial moment, the body-fixed frame of the leader spacecraft coincides with the orbital coordinate system, and its initial pose configuration and initial velocity are e position vector is expressed in the ECI system, and its unit iskm, and the velocity vector is expressed in the body-fixed coordinate system, and the units are rad/s and km/s. e definition of the initial pose configuration and initial velocity parameters of the follower spacecraft relative to the leader spacecraft are shown in Table 2.
e desired pose configuration and desired velocity of the follower spacecraft relative to the leader spacecraft are shown in Table 3.
In other words, the control goal is to keep the attitude synchronization between the follower spacecraft and the leader spacecraft, hover under it, and the relative velocity of the two is zero. e uncertain part of the mass and inertia matrix and the external disturbances are selected as follows: During the simulation, the maximum output of the reaction flywheel and thruster are 1 N.m and 10 N, respectively.
at is, the boundary of control force is [− 10, 10]N, and the control torque is limited to [− 1, 1]N · m. e parameters of the controllers are chosen as in Table 4. e specific fault types of each flywheel and each thruster are shown in Table 5. Figures 4-7 show the output of AF-MFixed and its comparison with AF-Finite and AF-Fixed under the normal condition of the actuator; Figures 8-11 show the output with the actuator fault. e above output results verify the stability analysis of the proposed control scheme. Figure 4 illustrates the pose configuration and the velocity tracking error of AF-MFixed without actuator fault. It can be seen that the attitude and angular velocity tracking errors quickly converge to the equilibrium state within 18 s, and the convergence accuracy is finally maintained within 1 × 10 − 4 deg and 2 × 10 − 5 deg/s, respectively; the position and translational velocity tracking errors quickly converge to the equilibrium state within 60 s, and the convergence accuracy is finally maintained within 8 × 10 − 5 m and 1 × 10 − 6 m/s, respectively. Figure 5(a) shows the output torque of each flywheel and Figure 5(b) shows the output force of each thruster under normal conditions. It is evident that all the actuator outputs are bounded. Figure 6 shows the comparison of the pose configuration and the velocity tracking error norms among AF-Finite, AF-Fixed, and AF-MFixed without actuator fault. It can be seen from Figures 6(a) and 6(b) that the convergence rates and convergence accuracy of the three methods are almost the same for attitude and angular velocity tracking errors. From Figures 6(c) and 6(d), as for position and translational velocity tracking errors, AF-MFixed and AF-Fixed are superior to AF-Finite in convergence rates, AF-Finite have the highest convergence accuracy of position tracking error, and the minimum overshoot of translational velocity tracking error, but it has the lowest convergence accuracy of translational velocity tracking error, and the control performance of AF-Fixed is between the other two. In summary, AF-MFixed has more obvious advantages in terms of rapidity than AF-Finite and accuracy of control performance than AF-Fixed, which also confirms the analysis and discussion of AF-MFixed in Remark 2. Figure 7 shows the comparison of the integration of control force and torque of the three methods without actuator fault. e integral of the control output often represents the control energy. From Figure 7(a), we can know that AF-Finite has the largest energy consumption of control force (integration of control force) and AF-Fixed and AF-MFixed have the similar control force energy consumption. However, it can be seen from Figure 7(b) that there is no significant difference in the energy consumption of control torque (integration of control torque) among the three. erefore, AF-MFixed also has great advantages in reducing control energy consumption. Figure 8 illustrates the pose configuration and the velocity tracking error of AF-MFixed with actuator fault. It can be seen that the attitude and angular velocity tracking errors quickly converge to the equilibrium state within 20 s, and the convergence accuracy is finally maintained within 3 × 10 − 4 deg and 2 × 10 − 5 deg/s, respectively. e position tracking errors quickly converge to the equilibrium state within 60 s, and the convergence accuracy is finally maintained within 4 × 10 − 3 m and tends to decrease; the translational velocity tracking errors quickly converge to the equilibrium state within 75 s, and the convergence accuracy is finally maintained within 2 × 10 − 5 m/s. Figure 9(a) shows the output torque of each flywheel, and Figure 9(b) shows the output force of each thruster under fault conditions. It is evident that all the actuator outputs are bounded, and the output curves well reflect the types of the fault. It is worth noting that the control torques and forces do not vanish completely when the control goal is     achieved because they also need to compensate for the total disturbances and actuator faults to keep the relative pose configuration between the follower spacecraft and the leader spacecraft. Figure 10 shows the comparison of the pose configuration and the velocity tracking errors norms among AF-Finite, AF-Fixed, and AF-MFixed with actuator fault. It can be seen from Figures 10(a) and 10(b) that the convergence rates of the three methods are almost the same for attitude and angular velocity tracking errors, but AF-MFixed has a significant advantage over AF-Fixed in terms of convergence accuracy for attitude tracking error. We can know from  accurately than AF-Fixed at the same convergence speed. In summary, AF-MFixed has more obvious advantages in terms of rapidity and accuracy of control performance, which also confirms the analysis and discussion of AF-MFixed in Remark 2. Figure 11 shows the comparison of the integration of control force and torque of the three methods with actuator fault. From Figure 11(a), we can know that AF-Finite has the lowest energy consumption of control force, and AF-MFixed is slightly higher than that of AF-Fixed in order to achieve fast convergence performance. However, it can be seen from Figure 11(b) that AF-Fixed has the lowest energy consumption of control torque, and AF-Finite is slightly higher than that of AF-Fixed, but the difference in the   16 Mathematical Problems in Engineering energy consumption of control torque is not obvious among the three compared with energy consumption of control force.

Conclusions
In this paper, adaptive fuzzy modified Fixed-time faulttolerant control schemes on SE(3) for coupled spacecraft were proposed to solve the attitude and position tracking problem with external disturbances, model uncertainties, and actuator faults simultaneously. From the comparative analysis of the three control strategies, we can see that AF-MFixed can achieve the control goals of fast convergence and higher tracking accuracy; the settling time of the closedloop tacking system can be independent of the initial states. e integrated attitude and position modeling method based on Lie group SE(3) is simple and can be applied to solve the problem of 6-DOF in practical aerospace engineering. e   fuzzy adaptive control scheme can estimate the total disturbances and fault information with high accuracy. e parameter tuning of the proposed algorithm is simple, avoiding the tedious adjustment of too many parameters. Moreover, the algorithm is suitable for both actuator failure and normal condition, and the control performance under fault condition will be slightly lower than that under normal condition, which also shows that the robustness proposed in this paper has a strong advantage, and it has potential engineering application value. However, the practical problem is that the actuator fault information is difficult to obtain in real time.
e establishment of a fault observer for fault diagnosis will be the work of the author in the future, and the estimated fault information will be applied to the design of fault-tolerant controllers. In addition, the next research content of this paper will consider the influence of sensor measurement noise and fault, and design a more robust fault-tolerant controller.

Data Availability
e data used to support the findings of this study are included in this paper.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.