Next Article in Journal
Low-Profile High-Efficiency Transmitarray Antenna for Beamforming Applications
Previous Article in Journal
SP-YOLO-Lite: A Lightweight Violation Detection Algorithm Based on SP Attention Mechanism
Previous Article in Special Issue
Sign Language Gesture Recognition and Classification Based on Event Camera with Spiking Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Collision Avoidance Second Order Sliding Mode Control of Satellite Formation with Air-Floated Platform Semi-Physical Simulation

1
Shanghai Electro-Mechanical Engineering Institute, Shanghai 201109, China
2
School of Astronautics, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Electronics 2023, 12(14), 3179; https://doi.org/10.3390/electronics12143179
Submission received: 24 April 2023 / Revised: 7 July 2023 / Accepted: 17 July 2023 / Published: 21 July 2023
(This article belongs to the Special Issue Design, Dynamics and Control of Robots)

Abstract

:
As the number of satellites in orbit increases, the issue of flight safety in spacecraft formation orbit control has become increasingly prominent. With this in mind, this paper designs a second-order terminal sliding mode controller for spacecraft formation obstacle avoidance based on an artificial potential function (APF). To demonstrate the effectiveness of the controller, this paper first constructs a Lyapunov function to prove its stability and then verifies its theoretical validity through numerical simulation. Finally, a satellite simulator is used for semi-physical simulation to verify the practical effectiveness of the controller proposed in this paper.

1. Introduction

Spacecraft formation is a relatively new mode of flight [1]. The members of the formation are physically separated but functionally integrated as a collective. Compared to traditional large spacecraft, spacecraft formations boast superior performance and more integrated functions. Through coordinated and cooperative efforts, members of the formation can accomplish tasks that cannot be completed by traditional single spacecraft, such as distributed aperture radar, stereo-imaging for earth observation, and deep space exploration [2,3,4].
Numerous studies have been conducted on the cooperative control of spacecraft formation flight. In the early stages, Lawton [5] categorized control strategies for spacecraft formation flight into three types: master–slave mode, virtual structure mode, and behavior mode. Scharf et al. [6] summarized the control strategy based on subsequent research results and added multi-input multi-output mode and cyclic mode to Lawton’s classification. In addition to these methods, other control methods include reference trajectory planning, passive decomposition, and consensus theory-based control.
In fact, the master–slave mode is the most extensively studied control structure for spacecraft formation. It employs a hierarchical arrangement of each spacecraft controller and simplifies formation control into a tracking problem for each spacecraft [7]. Research on master–slave mode formation control initially focused on the robot formation problem [8]. Wang et al. [9] were the first to extend the master–slave mode strategy to the spacecraft formation problem and designed the attitude cooperative control rate for formation keeping and neighboring spacecraft. Subsequently, the master–slave strategy has been widely used in research on spacecraft formations.
Wang et al. [10] applied the master–slave mode to the autonomous rendezvous and docking process of formations and discussed potential problems within it. Xu et al. [11] designed a sliding mode tracking control algorithm for the attitude–orbit coupling 6-DOF satellite. Stansbery et al. [12] used a state-dependent Riccati equation to design a full-state feedback nonlinear controller to simultaneously realize translational and rotational six-DOF control. Park et al. [13] used a state-dependent Riccati equation as a nonlinear controller to study optimal reconfiguration and the maintenance of satellites flying in formation. Massari et al. [14] applied SDRE technology to the attitude and orbit coupling control of the Proba-3 formation mission. Kumar et al. [15] studied a hybrid linear/nonlinear controller composed of linear control and time-optimal bang-bang control method to achieve the effective maneuvering of formation spacecraft.
Considering non-ideal situations such as faults or disturbances in spacecraft formations, some scholars have also conducted research on robust control methods. The design process of the sliding mode control method is simple and has strong robustness, making it widely used in the aerospace field. However, traditional first-order sliding mode has problems such as chattering. Therefore, Feng et al. [16] proposed using a second-order sliding mode controller to solve these problems while improving accuracy and maintaining the strong robustness of the controller.
Mondal S et al. [17] designed a second-order terminal sliding mode controller while using an adaptive rate to estimate the uncertainty upper bound online and update controller parameters. Ran et al. [18] aimed at the finite-time maneuvering problem of spacecraft and designed an adaptive second-order terminal sliding mode controller based on the adaptive idea. Liu et al. [19] designed a second-order nonsingular terminal sliding mode controller to realize coupled spacecraft attitude–orbit control. Chen [20] and Geng [21] combined adaptive control with a super-twisting controller to estimate and update controller gain online, effectively enhancing the anti-interference ability of the controller for the on-orbit service problem of out-of-control spacecraft.
Spacecraft formations operate in a complex orbital environment with external disturbances such as Earth’s flat-rate perturbations, atmospheric perturbations, and solar pressure perturbations, as well as non-ideal situations such as uncertain model parameters. It is necessary to design a suitable control law to make the configuration stable and maintainable while meeting accuracy requirements. Now that the second-order sliding mode has been introduced in the aerospace field, it needs to be more deeply combined with corresponding mission requirements to solve practical problems. Therefore, this paper designs an adaptive second-order sliding mode controller considering system uncertainty to support subsequent experimental verification.
With the gradual development of the aerospace field, the number of on-orbit spacecraft is gradually increasing, and obstacles such as scrapped spacecraft and space debris pose a growing threat to the development of the aerospace industry [22]. Based on the above problems, Nair R et al. [23] designed a robust controller based on the potential function and sliding mode method to enable spacecraft formation to avoid obstacles. Actually, the idea of the artificial potential function was first proposed by Utkin et al. [24], and since then this method has been used in [25] for formation control. Palacios L et al. [26] proposed two algorithms for close-range maneuver control based on the hybrid linear quadratic regulator and artificial potential function, respectively.
The concept of the potential function was introduced by Oussama Khatib into the research on robot anti-collision motion control, such as the research by the authors of [27,28], who are committed to the researching sensor perception ability, and the control algorithm mainly adopts the typical simple form of artificial potential function. Starting from this simple algorithm, more complex ones can be developed. In addition, Dusan Stipanovic also contributed important results in the control of flying systems. The authors of [29] propose theoretical and experimental results on the control of multi-agent nonholonomic systems, design and implement a novel decentralized control scheme that achieves dynamic formation control and collision avoidance for a group of non-holonomic robots, derive a feedback law using Lyapunov-type analysis that guarantees collision avoidance and tracking of reference trajectory for a single robot, extend this result to the case of multiple nonholonomic robots, and show how different multi-agent problems are. The authors of [30] deal with the formation control problem without collisions for second-order multi-agent systems. The avoidance complementary component is formed by applying repulsive vector fields with an unstable focus structure. Using the well-known input-to-state stability property, the control law for second-order agents is derived in a constructive manner starting from the first-order case.
By combining optimal sliding mode control technology with the artificial potential function, a new spacecraft rendezvous control method was proposed, and the stability of the system was proved by Lyapunov theory. Boyarko G et al. [31] proposed and solved the minimum energy optimal rendezvous problem using Pontryagin’s minimum principle while considering collision constraints [32,33], all considered accidental collisions between formation members and designed controllers to enable formation members to avoid external obstacles. Qi et al. [34] proposed a traditional control sliding mode method with anti-collision navigation for the attitude and orbit control process of spacecraft autonomous rendezvous and docking, considering obstacle-avoidance situations.
With the development of one-rocket multi-satellite missions, the number of satellite members in formation ranges from a few to more than a dozen. When performing tasks such as tracking, configuration initialization, and configuration maintenance, the orbital control law must ensure that there will be no accidental collisions between formation members. On the other hand, it also requires formation members to have the ability to avoid external obstacles. Therefore, a series of flight safety issues such as obstacle avoidance and anti-collision during formation flight must be considered in the design of the control law. Most obstacle avoidance algorithms in the current literature on aircraft formation control use traditional sliding mode control combined with artificial potential functions. In this paper, the artificial potential function is further improved, and the proposed second-order sliding mode method is used for obstacle avoidance control.
It should be pointed out that although the above methods have achieved some achievements, proof of engineering effectiveness is lacking. Chen L [35] proposed an adaptive integral sliding mode-based aircraft formation controller. Based on this, a second-order terminal sliding mode obstacle avoidance controller based on the artificial potential function is proposed, and the stability of the system is proved by the Lyapunov method. Numerical simulations are carried out. Finally, the controller is used to control the satellite simulator to prove the engineering effectiveness of the proposed controller.
The chapters of this paper are arranged as follows: Section 2 models spacecraft formation dynamics. In Section 3, a second-order terminal sliding mode obstacle avoidance controller based on artificial potential function is proposed, and its stability is proved. In Section 4, numerical simulations are presented to demonstrate the effectiveness of the controller. In Section 5, a semi-physical simulation using a satellite simulator is performed to demonstrate the engineering effectiveness of the controller. Section 6 is the conclusion, which summarizes the whole paper.

2. Spacecraft Formation Orbit Dynamics Model

This chapter analyzes the uncertain factors such as the relative motion model of the elliptical orbit of the short-range spacecraft formation and the perturbation factors.
Considering the control input u and external disturbance d from onboard, the nonlinear dynamic model of the elliptical orbit spacecraft formation under local vertical-local horizontal (LVLH) is as follows:
M ρ ¨ + C ρ ˙ + G = u + d ,
where ρ = [ x , y , z ] T is the position vector of the secondary satellite relative to the primary satellite in the local vertical-local horizontal (LVLH); ω is the angular velocity of the primary satellite; M is the mass matrix; and the definition of C , G is as follows:
C = 2 ω M 0 1 0 1 0 0 0 0 0 G = M μ x R + ρ 3 , ω 2 μ x y + R R + ρ 3 1 R , ω 2 μ y z R + ρ 3 T ,
ρ ¨ = x ¨ 2 n y ˙ n y n 2 x i + y ¨ + 2 n x ˙ + n ˙ x n 2 y j + z ¨ k .
Considering the uncertainty of parameters, modeling errors, and external disturbances in the system, the relative motion model of the spacecraft formation in the elliptical orbit can be expressed as:
M 0 ρ ¨ + C 0 ρ ˙ + G 0 = u + d + Δ d + D ( ρ , ρ ˙ , ρ ¨ ) ,
where M 0 , C 0 , G 0 R 3 × 3 is the known nominal coefficient matrix of the system; Δ M , Δ C R 3 × 3 is the coefficient uncertainty caused by fuel consumption, liquid sloshing, etc.; u R 3 is the control input of the system; d R 3 is the external disturbance; and Δ d R 3 includes the actuator installation deviation and output deviation, etc. The uncertainty item D ρ , ρ ˙ , ρ ¨ = Δ M ρ ¨ Δ C ρ ˙ Δ G R 3 is recorded as the uncertain item for which the system does not contain an external interference.
Define D ¯ ρ , ρ ˙ , ρ ¨ = d + Δ d + D ρ , ρ ˙ , ρ ¨ as a compound term that includes uncertainties in system parameters, modeling errors, and external disturbances.
Equation (4) is finally rewritten as follows:
M 0 ρ ¨ + C 0 ρ ˙ + G 0 = u + D ¯ ( ρ , ρ ˙ , ρ ¨ ) .
Suppose D ¯ ρ , ρ ˙ , ρ ¨ and its derivative term D ¯ ˙ ( ρ , ρ ˙ , ρ ¨ ) are bounded. Let φ = M 0 1 D ¯ ˙ ( ρ , ρ ˙ , ρ ¨ ) , assuming it satisfies the following inequality:
φ = M 0 1 D ¯ ˙ < d m ,
where d m is an unknown positive scalar.
Let x = ρ , f ( x ) = C 0 ρ ˙ G 0 / M ; the orbital dynamics model of the spacecraft formation can be expressed as follows:
x ¨ t = f x + u + D ¯ .

3. Second-Order Terminal Sliding Mode Controller for Spacecraft Formation Based on Modified Artificial Potential Function

The concept of the artificial potential function (APF) is a virtual force method in nature. It is premised on the idea that a tracking spacecraft navigates within an imaginary potential field, and a scalar function is established to represent the relative position within the state space. The desired position corresponds to the global minimum value of this scalar function, while regions with higher potential functions indicate the hazardous areas for potential collisions during the tracking spacecraft’s movement. Theoretically, the negative gradient of the specifically designed APF offers a secure route for the spacecraft to reach its target position while evading obstacles.
In summary, the potential function is defined as follows:
U = U a t t + U r e p ,
where U a t t is the gravitational potential function and U r e p is the repulsive potential function.
The following quadratic function is chosen as the gravitational potential function U a t t .
U a t t = 1 2 ρ ρ d T P ρ ρ d ,
where P is a positive definite symmetric matrix and ρ d is the desired position.
Depending on the chosen gravitational potential function, when ρ ρ d is, U a t t 0 , its gradient is as follows:
r U a t t = ρ ρ d T P .
According to (10), the negative gradient vector of the gravitational potential function points to the desired position.
The value of the repulsive potential function significantly increases when approaching the hazardous area of the obstacle, whereas it decreases or even approaches zero when moving away from the hazardous area. In other words, the overall potential function U is primarily influenced by the repulsive potential function within the danger zone, and the attractive potential function takes effect once the spacecraft escapes the danger zone.
The hazardous area is considered as a sphere with the obstacle’s center of mass as the sphere’s center. The repulsive potential function can be defined as follows:
U r e p = 1 2 ψ i exp ρ r o , i T N i ρ r o , i σ i ,
where N i is a positive definite symmetric matrix; r o , i represents the centroid position of the ith obstacl; and ψ i and σ i , respectively, represent the height and radius of the repulsive potential function of the ith obstacle. According to the obstacles of different sizes and the safety factor, the selection of the scope of the danger area can be obtained by adjusting N i , ψ i and σ i .
Depending on the chosen repulsive potential function, U r e p max when ρ r o , i . Its gradient is as follows:
r U r e p = i = 1 n ψ i σ i exp ρ r o , i T N i ρ r o , i σ i ρ r o , i T N i .
According to Equation (12), the negative gradient vector of the gravitational potential function deviates from the obstacle.
When considering multiple obstacles, the artificial potential function model of spacecraft formation is as follows:
U = U r e p + U a t t = 1 2 ρ ρ d T P ρ ρ d + i = 1 n 1 2 ψ i exp ρ r o , i T N i ρ r o , i σ i .
In order to make the value of the artificial potential function U at the desired position ρ = ρ d zero and improve the efficiency of the obstacle avoidance process, the artificial potential function is modified as follows:
U = U r e p + U a t t = 1 2 ρ ρ d T P ρ ρ d + 1 2 ρ ρ d T M ρ ρ d · i = 1 n 1 2 ψ i exp ρ r o , i T N i ρ r o , i σ i .
According to the definition of the artificial potential function, it works by driving the tracking spacecraft from a place with a higher potential function value to a desired position with a lower potential function, and the control force generated at a certain position is determined according to the gradient at that position. If the repulsive force generated by the region with large potential function value is limited in magnitude, and the relative initial velocity of the tracking satellite and the obstacle is large, the tracking spacecraft cannot slow down to zero outside the danger zone in time, resulting in the danger of collision. Therefore, the relative position and the relative velocity of the tracking spacecraft and the obstacle should be considered simultaneously when designing the control law of obstacle avoidance navigation.
The sliding surface is defined as follows:
s = k ρ U + ρ ˙ ρ ˙ d = k ρ U + e ˙ ,
where:
ρ U = ρ ρ d T P + ρ ρ d T M · i = 1 n ψ i σ i exp ρ r o , i T N i ρ r o , i σ i ρ ρ d T M ρ ρ d · i = 1 n ψ i σ i exp ρ r o , i T N i ρ r o , i σ i ρ r o , i T N i .
Taking the first derivative of the sliding mode surface s with respect to time gives s ˙ as follows:
s ˙ = k ρ U ρ ρ ˙ + k ρ U ρ d ρ ˙ d + k i = 1 n ρ U ρ i r o , i r ˙ o , i + ρ ¨ ρ ¨ d ,
where:
ρ U ρ = P + M i = 1 n ψ i σ i exp ρ r o , i T N i ρ r o , i σ i 2 M ρ ρ d i = 1 n ψ i σ i exp ρ r o , i T N i ρ r o , i σ i ρ r o , i T N i + ρ ρ d T M ρ ρ d i = 1 n ψ i σ i 2 exp ρ r o , i T N i ρ r o , i σ i 2 N i ρ r o , i ρ r o , i T N i σ i N i ,
ρ U ρ d = P M i = 1 n ψ i exp ρ r o , i T N i ρ r o , i σ i + 2 M ρ ρ d i = 1 n ψ i σ i exp ρ r o , i T N i ρ r o , i σ i . ρ r o , i T N i ,
i = 1 n ρ U ρ i r o , i r ˙ o , i = 2 M ρ ρ d · i = 1 n ψ i σ i exp ρ r o , i T N i ρ r o , i σ i ρ r o , i T N i r ˙ o , i ρ ρ d T M ρ ρ d i = 1 n ψ i σ i 2 exp ρ r o , i T N i ρ r o , i σ i 2 N i ρ r o , i ρ r o , i T N i σ i N i r ˙ o , i .
Let us define:
Ω ρ φ ρ ρ ˙ + ρ U ρ d ρ ˙ d + i = 1 n ρ U ρ i r o , i r ˙ o , i .
At this point, the derivative of the sliding mode s is as follows:
s ˙ = k Ω + e ¨ .
According to Equations (7) and (8), the second derivative of the sliding mode can be expressed as follows:
s ¨ = k Ω ˙ + d e ¨ d t = k Ω ˙ + M 0 1 u ˙ d C 0 ρ ˙ + G 0 d t d ρ ¨ d d t + φ ,
where M 0 is the reversible nominal mass matrix defined above. Based on the idea of dynamic sliding mode, combined with the second-order sliding mode, the traditional sliding mode is introduced into the new switching function through differential processing. Select the switching function as follows:
σ = s + β 1 s ˙ p / q ,
where σ = σ 1 , σ 2 , σ 3 T , β = diag β 1 , β 2 , β 3 , β i > 0 , the definition of s p q is sgn ( s ) | s | p q , and p, q is positive odd and satisfies 1 < p / q < 2 . The derivative of Equation (25) is taken as follows:
σ ˙ = s ˙ + p / q β 1 s ˙ p / q 1 s ¨ = p / q β 1 s ˙ p / q 1 s ¨ + q / p β s ˙ 2 p / q .
Before designing the controller, the following lemma is given:
Lemma 1.
If there exists a continuous function V that satisfies the following conditions:
(1)
V is a positive definite function.
(2)
There are positive real numbers β > 0 , γ ( 0 , 1 ) and a region near the origin that satisfies:
V ˙ ( t ) + β V γ ( t ) 0 .
Then, the corresponding system can converge in finite time, and the convergence time t f is:
t f 1 β ( 1 γ ) V 1 γ ( 0 ) .
Lemma 2.
If the real numbers a 0 and b 0 form the two right sides of a right triangle, then the inequality a + b a 2 + b 2 holds according to the relationship of the side lengths of the right triangle.
After giving the above two lemmas, the form of the controller is given below.
In the process of configuration initialization of spacecraft formation, there are obstacles around the formation members that can be measured in relative position, relative velocity, and size. The second-order terminal sliding mode control law and adaptive update law in the following form are adopted, and the appropriate control parameters are selected to ensure that the formation members tend to the desired position while maintaining a safe anti-collision distance with respect to the obstacles.
u = 0 t u ˙ d t u ˙ = u ˙ nom + u ˙ com u ˙ nom = M 0 d ρ ¨ d d t + d C 0 ρ ˙ + G 0 d t k M 0 Ω ˙ q / p β M 0 s ˙ 2 p / q u ˙ com = k M 0 σ M 0 d ^ m + ε sgn ( σ ) d ^ ˙ m = η p q β 1 s ˙ p / q 1 σ .
Proof. 
Define the composite interference upper bound error as follows:
d ˜ m = d ^ m d m .
The Lyapunov function is chosen as follows:
V 1 = 1 2 σ T σ + γ 0 d ˜ m 2 ,
where γ > 0 is the controller parameter to be designed. We take the derivative of V 1 :
V ˙ 1 = σ T σ ˙ + γ 0 d ˜ m d ^ ˙ m = p q β 1 s ˙ p / q 1 σ T s ¨ + q p β s ˙ 2 p / q + γ 0 d ^ m d m d ^ ˙ m = p q β 1 s ˙ p / q 1 σ T M 0 1 u ˙ d C 0 ρ ˙ + G 0 d t d ρ ¨ d d t + φ + k Ω ¨ + q p β s ˙ 2 p / q + γ 0 η p q β 1 s ˙ p / q 1 σ d ^ m d m .
According to Equations (28) and (31):
V ˙ 1 p q β 1 s ˙ p / q 1 φ σ d ^ m σ + d m σ d m σ + γ 0 η d ^ m d m σ k σ 2 ε σ p q β 1 s ˙ p / q 1 d m φ σ p q β 1 s ˙ p / q 1 σ γ 0 η σ d ^ m d m = c 1 σ c 2 d ˜ m ,
where:
c 1 p q β 1 s ˙ p / q 1 d m φ c 2 p q β 1 s ˙ p / q 1 σ γ 0 η σ .
Equation (32) is further rewritten into the following form:
V ˙ 1 c 1 2 σ 2 c 2 2 γ γ 2 d ˜ m 2 min c 1 2 , c 2 2 γ × σ 2 + γ 2 d ˜ m 2 .
According to Lemma 2, σ 2 0 , γ 0 2 d ˜ m 2 0 , σ 2 + γ 0 2 d ˜ m 2 1 2 σ T σ + γ 0 d ˜ m 2 , which is:
V ˙ 1 c V 1 1 2 ,
where c = min c 1 2 , c 2 2 γ > 0 ; according to Lemma 1, for any initial state σ ( 0 ) 0 , the system can converge to σ = 0 in a finite time.
Further, we need to prove that sliding mode s ( t ) can also converge to s = 0 in a finite time at σ = 0 .
The Lyapunov function is chosen as follows:
V 2 = 1 2 s T s .
We take the derivative of V 2 , which is associated with Equation (25), and the result is as follows:
V ˙ 2 = s T s ˙ = s T β q / p s q / p β min q / p 2 s 2 ( p + q ) / p = c s V α ( t ) ,
where β min is the smallest absolute value component of β ; c s = β min q / p 2 ( p + q ) / ( 2 p ) ; α = ( p + q ) / 2 p ( 0 , 1 ) ; and, according, to Lemma 1, the sliding mode s ( t ) can converge in finite time, which is:
s = k r U + ρ ˙ ρ ˙ d = 0 ,
and
U ˙ = ρ U · ρ ˙ + ρ d U · ρ ˙ d + i = 1 n r o , i U · r ˙ o , i .
In fact, the relative velocity of the obstacle is relatively small compared to the relative velocity of the spacecraft. When the desired position is a constant, only the first term can be considered on the right-hand side of Equation (40), which is:
U ˙ ρ U · ρ ˙ .
Simultaneous Equations (38) and (40) can be obtained as follows:
U ˙ = k ρ U 2 .
The above equation shows that t , U ˙ 0 , ρ U 0 ; since r d is the only stable equilibrium point of U, we have U 0 , and we can further obtain that the system state converges to the desired state, which is:
lim t ρ ρ d , lim t v = 0 .
In summary, by choosing the sliding mode of Equation (15) and the controller of Equation (28), the system can converge in finite time and can complete the obstacle avoidance task. □

4. Simulation and Results

The primary satellite orbit parameters are a = 7000 km , μ = 3.986 × 10 5 km 3 · s 2 , e = 0.3 , θ 0 = 0 rad , while n = μ a 1 e 2 / r 2 is the instantaneous angular velocity of the primary satellite and r is the distance between the primary satellite and the center of the Earth.
In order to verify the obstacle avoidance ability of the control algorithm, the following numerical simulation is designed. For three slave satellites with initial position deviation, a master–slave control strategy is adopted to realize configuration initialization and the configuration maintenance of an equilateral triangle. Among them, the appropriate position is selected, and two obstacles are set between the initial and desired positions of the two slave satellites. The initialization of the simulation is shown in the Table 1. The simulation verifies whether the two slave satellites can successfully avoid the obstacles at a safe distance during the configuration initialization process.
The composite interference is defined as in Equation (44), and the upper bound of the composite interference is estimated by the adaptive rate. Through numerical simulation, the robustness of the control method under uncertain interference is verified.
d = 1.2 × 10 3 × sin ( 0.01 t ) cos ( 0.01 t ) sin ( 0.01 t ) .
The danger area of the designed obstacle in the simulation is a sphere with radius R = 50 m , and the artificial potential function parameters are set as shown in the Table 2.
The controller parameters are set as shown in the Table 3.
The simulation results using the second order sliding mode terminal controller are shown as follows. Figure 1 and Figure 2 show the three-dimensional motion trajectory of the spacecraft relative to the obstacle from two perspectives, respectively. The blue shadow sphere in the figure indicates the spatial influence range of the obstacle, and the radius of each sphere is set at 50 m . The spacecraft approaches from the starting point to the end point and finally reaches the desired position to form a configuration and maintain it.
Figure 3 and Figure 4 represent the relative distance between the spacecraft and each obstacle. When the spacecraft approaches the obstacle at nearly 50 m , the distance expands rapidly, keeps away from the obstacle or close to the edge of the danger area of the obstacle, and approaches the expected position. The distance in the whole process is greater than 50 m , indicating that the spacecraft avoids the collision danger area of the obstacle. It can be understood that there is no collision between the spacecraft and the two obstacles.
Taking Sat1 as an example, it can be seen that the spacecraft reaches the desired state after avoiding obstacles to achieve the purpose of final control. It can be seen from Figure 5 that the obstacle avoidance time is about 1000 s. Then, the steady-state error of configuration maintenance is within 0.005 m. It can be seen from Figure 6 that the controller output increases when the spacecraft is close to the obstacle.
In order to facilitate the comparison of the effects of this method, the following traditional sliding mode controller designed by the exponential approach law is used for numerical simulation of obstacle avoidance:
s ˙ = k Ω + e ¨ u = k ρ φ ρ ρ ˙ k ρ U ρ d ρ ˙ d k i = 1 n ρ U ρ i r o , i r ˙ o , i λ 1 s λ 2 sgn ( s ) α f ( ρ , v ) + r d .
The initial parameters of the simulation are consistent with the above. The parameters of the traditional sliding mode controller (SMC) are designed as Table 4:
The simulation results are as Figure 7, Figure 8, Figure 9 and Figure 10:
It can be seen that when the spacecraft approaches an obstacle, the output of the controller increases significantly, and the entire process of avoiding the obstacle lasts for a long time. Compared with the traditional sliding mode control, the upper limit of the output of the control amount of the second-order sliding mode control is smaller, the duration of the control output is shorter, and the second-order sliding mode control weakens the vibration phenomenon of the controller.
The simulation results in this section show that the second-order sliding mode obstacle avoidance controller proposed by Section 3 can realize the spacecraft formation formation reconstruction task in the presence of obstacles.

5. Semi-Physical Simulation Experiment of Satellite Simulator

In order to prevent the advanced control methods of spacecraft formations from being disconnected from practical applications, it is necessary to develop control hardware and software related to simulating operation in the space environment, establish an experimental platform that can simulate satellite formation control, and study new types of practical control methods that are conducive to engineering practice. The air flotation platform is used to simulate the flight of satellites in orbit and to verify the accuracy of satellite formation modeling and the performance of the formation controller. It has flourished due to its low cost, high reliability, and wide range of applications.
In this chapter, a semi-physical simulation experiment is carried out with a satellite simulator to verify the effectiveness of the second-order terminal sliding mode spacecraft obstacle avoidance controller proposed by Section 3, and the results are analyzed.

5.1. Establishment of Satellite Simulator Simulation Environment

The simulation environment of the satellite simulator mainly includes the satellite simulator subsystem, measurement subsystem, control subsystem, and marble air floating platform. The overall schematic of the semi-physical simulation environment is shown in Figure 11, and the real scene is shown in Figure 12.
Through the coordinated networking work of the above subsystems, an experimental platform that can simulate the control of space networked satellite formations is formed. Semi-physical simulation will be carried out through the coordinated control method proposed; in theory, combined with the problems encountered in the actual experiment process and analysis, the control method will be further improved, and the advanced control method will truly be integrated with the actual control application to achieve the desired control effect.
The satellite simulator subsystem is the most important part of the test, and the components of this subsystem are shown in Table 5.
The satellite simulator is connected with a rigid reflecting ball to establish a rigid body. The vision camera captures the shape of the rigid body and calculates the position and attitude information of the satellite simulator. The host computer communicates with the communication component of the satellite simulator through the WiFi network to obtain the real-time position and attitude information of the satellite simulator and display it.
The control subsystem uses the control law to calculate the reference trajectory and transmits the desired trajectory to the satellite simulator through the Wifi signal. The lower computer of the satellite simulator controls the jet nozzle through the internal feedback loop to air injection so that the satellite simulator can track the reference trajectory in real time, to complete the control goal. The control flow of the satellite simulator is shown in Figure 13.
It is worth stating that the reference trajectory is a set of time-discrete position data, and the time interval is the jet period of the satellite simulator jet nozzle.
The jet nozzle distribution of the satellite simulator is shown in Figure 14. The satellite simulator has a total of 8 jet nozzles. The air injection of the jet nozzle is controlled by controlling the switch of the solenoid valve, so that the force in both x and y directions and the torque in z direction can be provided.
Using the assembly method of the above-mentioned control actuator, the three-degree-of-freedom movement can be easily and effectively simulated, including two shifts in the horizontal plane and rotation around the z-axis. This satellite simulator can verify the coupled motion process of attitude and trajectory. This experiment is mainly used for obstacle avoidance verification of the two-position plane. The rotational movement around the z-axis only needs to maintain the angular orientation.

5.2. Validation of the Controller

In order to enhance the credibility of the test, the parameters of the obstacle avoidance controller are selected in accordance with Table 3. The artificial potential function parameters are shown in Table 6. The initial position, desired position, and radius settings are shown in Table 7.
The marble air floating platform consists of two 6 m × 3 m marble pieces stitted into a complete 6 m × 6 m platform; the satellite simulator, through air injection, makes the air foot and marble air floating platform produce a layer of stable air film, making the satellite simulator in the floating microgravity produce a frictionless state and a simulation of the space weightlessness environment, to achieve three degrees of freedom of free movement without external force.
Through the coordination of the above subsystems, an experimental platform is constructed to simulate the formation control of space network satellites, so that the effectiveness of the proposed control algorithm can be experimentally verified.
The controller parameters are set as shown in the following table.
The obstacle avoidance motion trajectory of the satellite simulator and the relative motion error are shown in Figure 15 and Figure 16. Figure 17 is a video screenshot during the obstacle avoidance process of the satellite simulator. The tracking error between the satellite simulator and the desired position at the next time instant is shown in Figure 18. It can be seen that under the action of the second-order terminal sliding mode obstacle avoidance controller proposed in Section 3, the satellite simulator can complete the control task in the predetermined situation.
In the process of carrying out control experiments, due to limited engineering technologies such as manufacturing, processing, and assembly, the thrust output of the simulator has large installation errors and output errors, and the quality and moment of inertia of the simulator are subject to uncertain changes with the consumption of working fluid. At the same time, there are various external disturbances such as friction in the experimental environment, resulting in inevitable differences between the experimental results and the numerical simulation results. For the numerical simulation results, it is necessary to further improve the accuracy of the modeling based on the experimental results, fully consider the uncertainty factors such as changes in internal parameters of the system and external interference, and study a robust control algorithm similar to this method to reduce the sensitivity to interference. At the same time, we must also consider the limitations of hardware computing power and practical application. For example, improper handling of the integration link in the control algorithm will cause the experimental control process to diverge and produce unpredictable hazards.
In summary, the second-order terminal sliding mode obstacle avoidance controller proposed in this paper is effectively verified by the test of the satellite simulator.

6. Results

This paper first combines adaptive control with second-order terminal sliding mode control and designs an adaptive second-order terminal sliding mode controller (SOTSMC) to solve the uncertainty control problem of the system.
Then, flight safety issues such as obstacle avoidance and anti-collision of satellite formations during configuration initialization or reconstruction are considered; the artificial potential function guidance is combined with the adaptive second-order terminal sliding mode controller; and a navigation controller is designed so that the formation can achieve high-precision target tracking, configuration initialization, configuration reconstruction, and other tasks at the same time, to enable the avoidance of obstacles with measurable position and speed. Numerical simulation analysis shows that the second-order sliding mode control proposed in this paper has the ability to avoid obstacles after combining with the navigation algorithm. Compared with the traditional sliding mode control, the control accuracy is improved and the vibration phenomenon of the controller output is significantly weakened.
Finally, using the three-axis satellite simulator, the application of the control algorithm proposed in this paper in obstacle avoidance scenarios was physically tested. Based on the test results, combined with numerical simulation and actual situation, a comparative analysis is made to verify the effectiveness of the controller and the engineering application value. Although limited by the allocation of hardware resources and the influence of some non-ideal test environments, some control processes in the actual application results do not achieve the ideal results, but the physical tests have also further combined the theory of the controller with the actual engineering application, highlighting some practical problems that were missed in the theoretical analysis, so that in the future controller design process, theoretical research can make targeted progress towards solving actual control problems.
On the basis of the research of this thesis, in-depth research is required from the following aspects in the follow-up.
The integration process in the second-order sliding mode control weakens the controller vibration problem caused by the symbol term and causes significant harm to actual engineering applications. Experimental verification and analysis show that when the manufacturing accuracy of the controller is not high, there is an error between the actual output of the controller and the nominal output. The error will gradually accumulate due to integration during the control process, affecting the subsequent control process.Therefore, it is necessary to further study the integration saturation and controller fault tolerance of the second-order sliding mode control in practical applications in the future.
In the process of obstacle avoidance and navigation, if the output amplitude of the controller is limited, and the initial speed of the formation members is high, the speed cannot be reduced to zero in time before reaching the dangerous area, so there will be a risk of collision. Therefore, the output saturation of the controller needs to be further considered in the design process of the controller. At the same time, it is necessary to further consider the research on the control algorithm of the controlled object under drive conditions.

Author Contributions

Methodology, J.Z.; Validation, Y.W.; Investigation, J.J.; Writing—original draft, Y.W.; Writing—review & editing, J.Z.; Supervision, C.C. and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

For the consideration of data privacy, the experimental data of the semi-physical simulation will not be disclosed in this paper.

Conflicts of Interest

The authors hereby declare that there is no conflict of interest in this article.

References

  1. Beard, R.W.; Lawton, J.; Hadaegh, F.Y. A coordination architecture for spacecraft formation control. IEEE Trans. Control. Syst. Technol. 2001, 9, 777–790. [Google Scholar] [CrossRef] [Green Version]
  2. Di Mauro, G.; Lawn, M.; Bevilacqua, R. Survey on guidance navigation and control requirements for spacecraft formation-flying missions. J. Guid. Control. Dyn. 2018, 41, 581–602. [Google Scholar] [CrossRef]
  3. Lee, D.; Sanyal, A.K.; Butcher, E.A. Asymptotic tracking control for spacecraft formation flying with decentralized collision avoidance. J. Guid. Control. Dyn. 2015, 38, 587–600. [Google Scholar] [CrossRef] [Green Version]
  4. Li, P.; Zhu, Z.H.; Meguid, S.A. State dependent model predictive control for orbital rendezvous using pulse-width pulse-frequency modulated thrusters. Adv. Space Res. 2016, 58, 64–73. [Google Scholar] [CrossRef]
  5. Lawton, J.R.T. A Behavior-Based Approach to Multiple Spacecraft Formation Flying; Brigham Young University: Provo, UT, USA, 2000. [Google Scholar]
  6. Scharf, D.P.; Hadaegh, F.Y.; Ploen, S.R. A survey of spacecraft formation flying guidance and control. Part ii: Control. In Proceedings of the 2004 American Control Conference, Boston, MA, USA, 30 June–2 July 2004; Volume 4, pp. 2976–2985. [Google Scholar]
  7. Jun, S.; Jing, H.; Zhang, X.; Huang, T. Dynamics and control of spacecraft formation flying in earth orbit. Mech. Eng. 2019, 41, 117–136. [Google Scholar]
  8. Wang, P.K.C. Navigation strategies for multiple autonomous mobile robots moving in formation. J. Robot. Syst. 1991, 8, 177–195. [Google Scholar] [CrossRef]
  9. Wang, P.K.C.; Hadaegh, F.Y. Coordination and control of multiple micro-spacecraft moving in formation. J. Astronaut. Sci. 1996, 44, 315–355. [Google Scholar]
  10. Wang, P.K.C.; Hadaegh, F.Y. Formation flying of multiple spacecraft with autonomous rendezvous and docking capability. IET Control. Theory Appl. 2007, 1, 494–504. [Google Scholar] [CrossRef]
  11. Xu, Y.; Tatasch, A.; Fitz-Coy, N. Chattering free sliding mode control for a 6 DOF formation flying mission. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, San Francisco, CA, USA, 15–18 August 2005; p. 6464. [Google Scholar]
  12. Stansbery, D.; Cloutier, J. Nonlinear control of satellite formation flight. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Denver, CO, USA, 14–17 August 2000; p. 4436. [Google Scholar]
  13. Park, H.-E.; Park, S.; Choi, K. Satellite formation reconfiguration and station-keeping using state-dependent Riccati equation technique. Aerosp. Sci. Technol. 2011, 15, 440–452. [Google Scholar] [CrossRef]
  14. Massari, M.; Zamaro, M. Application of SDRE technique to orbital and attitude control of spacecraft formation flying. Acta Astronaut. 2014, 94, 409–420. [Google Scholar] [CrossRef]
  15. Kumar, B.S.; Ng, A. Time-optimal low-thrust formation maneuvering using a hybrid linear/nonlinear controller. J. Guid. Control. Dyn. 2009, 32, 343–347. [Google Scholar] [CrossRef]
  16. Feng, Y.; Han, F.; Yu, X. Chattering free full-order sliding-mode control. Automatica 2014, 50, 1310–1314. [Google Scholar] [CrossRef]
  17. Mondal, S.; Mahanta, C. Adaptive second order terminal sliding mode controller for robotic manipulators. J. Frankl. Inst. 2014, 351, 2356–2377. [Google Scholar] [CrossRef]
  18. Ran, D.; Xiao, B.; Sheng, T.; Chen, X. Adaptive finite time control for spacecraft attitude maneuver based on second-order terminal sliding mode. Proc. Inst. Mech. Eng. Part J. Aerosp. Eng. 2017, 231, 1415–1427. [Google Scholar] [CrossRef]
  19. Liu, H.; Shi, X.; Zhang, J. Adaptive second-order terminal sliding mode control for approaching non-cooperative targets. Syst. Eng. Electron. Technol. 2016, 38, 2353–2360. [Google Scholar]
  20. Chen, B.; Geng, Y. Adaptive sliding mode controller design for in-orbit service of runaway spacecraft. J. Aeronaut. 2015, 36, 1639–1649. (In Chinese) [Google Scholar]
  21. Chen, B.; Geng, Y. Super twisting controller for on-orbit servicing to non-cooperative target. Chin. J. Aeronaut. 2015, 28, 285–293. [Google Scholar] [CrossRef] [Green Version]
  22. Bai, X.; Chen, L.; Zhang, Y.; Tang, G. Survey on space target collision warning technology. J. Astronaut. 2013, 34, 1027–1039. [Google Scholar]
  23. Nair, R.R.; Behera, L.; Kumar, V.; Jamshidi, M. Multisatellite formation control for remote sensing applications using artificial potential field and adaptive fuzzy sliding mode control. IEEE Syst. J. 2014, 9, 508–518. [Google Scholar] [CrossRef]
  24. Utkin, V.I.; Drakunov, S.V.; Hashimoto, H.; Harashima, F. Robot path obstacle avoidance control via sliding mode approach. In Proceedings of the IROS’91: IEEE/RSJ International Workshop on Intelligent Robots and Systems’ 91, Osaka, Japan, 3–5 November 1991; pp. 1287–1290. [Google Scholar]
  25. Hu, Q.; Dong, H.; Zhang, Y.; Ma, G. Tracking control of spacecraft formation flying with collision avoidance. Aerosp. Sci. Technol. 2015, 42, 353–364. [Google Scholar] [CrossRef]
  26. Palacios, L.; Ceriotti, M.; Radice, G. Close proximity formation flying via linear quadratic tracking controller and artificial potential function. Adv. Space Res. 2015, 56, 2167–2176. [Google Scholar] [CrossRef] [Green Version]
  27. Flacco, F.; Kröger, T.; Luca, A.D.; Khatib, O. A depth space approach to human-robot collision avoidance. In Proceedings of the 2012 IEEE International Conference on Robotics and Automation, St Paul, MN, USA, 14–18 May 2012; pp. 338–345. [Google Scholar]
  28. Minguez, J.; Montano, L.; Khatib, O. Reactive collision avoidance for navigation with dynamic constraints. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Lausanne, Switzerland, 30 September–4 October 2002; Volume 1, pp. 588–594. [Google Scholar]
  29. Silvia, M.; Stipanović, D.M.; Graunke, C.R.; Intlekofer, K.A.; Spong, M.W. Formation control and collision avoidance for multi-agent non-holonomic systems: Theory and experiments. Int. J. Robot. Res. 2008, 27, 107–126. [Google Scholar]
  30. Francisco, F.J.; Avilés, D.; Aranda-Bricaire, E. Formation Control for Second-Order Multi-Agent Systems with Collision Avoidance. Machines 2023, 11, 208. [Google Scholar]
  31. Boyarko, G.; Yakimenko, O.; Romano, M. Optimal rendezvous trajectories of a controlled spacecraft and a tumbling object. J. Guid. Control. Dyn. 2011, 34, 1239–1252. [Google Scholar] [CrossRef]
  32. Ning, Z.; Xia, Y. Coordination control design for formation reconfiguration of multiple spacecraft. IET Control. Theory Appl. 2015, 9, 2222–2231. [Google Scholar]
  33. Feng, L.C.; Bai, Y.Z.; Chen, X.Q. Finite-time sliding mode control for spacecraft obstacle avoidance rendezvous. J. Astronaut. 2016, 37, 1342–1348. [Google Scholar]
  34. Li, Q.; Yuan, J.; Wang, H. Sliding mode control for autonomous spacecraft rendezvous with collision avoidance. Acta Astronaut. 2018, 151, 743–751. [Google Scholar] [CrossRef]
  35. Chen, L.; Baodi, X. Spacecraft formation controller design based on adaptive intergated sliding model control. Air Space Def. 2021, 4, 65–70. [Google Scholar]
Figure 1. The trajectory of the spacecraft relative to the obstacle.
Figure 1. The trajectory of the spacecraft relative to the obstacle.
Electronics 12 03179 g001
Figure 2. The trajectory of spacecraft relative to obstacles (X-Y projection).
Figure 2. The trajectory of spacecraft relative to obstacles (X-Y projection).
Electronics 12 03179 g002
Figure 3. The distance of Sat1 with respect to the obstacle.
Figure 3. The distance of Sat1 with respect to the obstacle.
Electronics 12 03179 g003
Figure 4. The distance of Sat2 with respect to the obstacle.
Figure 4. The distance of Sat2 with respect to the obstacle.
Electronics 12 03179 g004
Figure 5. The relative position error of Sat1.
Figure 5. The relative position error of Sat1.
Electronics 12 03179 g005
Figure 6. The Sat1 controller output.
Figure 6. The Sat1 controller output.
Electronics 12 03179 g006
Figure 7. The three-dimensional trajectory relative to obstacles (SMC).
Figure 7. The three-dimensional trajectory relative to obstacles (SMC).
Electronics 12 03179 g007
Figure 8. The distance of satellite with respect to the obstacle.
Figure 8. The distance of satellite with respect to the obstacle.
Electronics 12 03179 g008
Figure 9. The relative position error of satellite.
Figure 9. The relative position error of satellite.
Electronics 12 03179 g009
Figure 10. The SMC controller output.
Figure 10. The SMC controller output.
Electronics 12 03179 g010
Figure 11. The overall composition of the simulation system.
Figure 11. The overall composition of the simulation system.
Electronics 12 03179 g011
Figure 12. The real scene of semi-physical simulation environment.
Figure 12. The real scene of semi-physical simulation environment.
Electronics 12 03179 g012
Figure 13. The control flow chart of the satellite simulator.
Figure 13. The control flow chart of the satellite simulator.
Electronics 12 03179 g013
Figure 14. The satellite simulator jet nozzle distribution.
Figure 14. The satellite simulator jet nozzle distribution.
Electronics 12 03179 g014
Figure 15. The obstacle avoidance trajectory of the satellite simulator.
Figure 15. The obstacle avoidance trajectory of the satellite simulator.
Electronics 12 03179 g015
Figure 16. The relative position error of satellite simulators.
Figure 16. The relative position error of satellite simulators.
Electronics 12 03179 g016
Figure 17. The obstacle avoidance trajectory of the satellite simulators.
Figure 17. The obstacle avoidance trajectory of the satellite simulators.
Electronics 12 03179 g017
Figure 18. The tracking error of the satellite simulator.
Figure 18. The tracking error of the satellite simulator.
Electronics 12 03179 g018
Table 1. Selection of initial/desired position. ( m or m / s ).
Table 1. Selection of initial/desired position. ( m or m / s ).
Initial StatePosition of ObstaclesDesired State
ρ 10 = [ 200 , 200 , 100 ] T r o , 11 = [ 130 , 180 , 75 ] T ρ d 1 = [ 0 , 50 , 0 ] T
ρ 20 = [ 243.3 , 250 , 110 ] T r o , 12 = [ 80 , 70 , 30 ] T ρ d 2 = 50 [ cos π 6 , sin π 6 , 0 ] T
ρ 30 = [ 163.3 , 105 , 100 ] T r o , 21 = [ 163.3 , 185 , 70 ] T ρ d 3 = 50 [ cos π 6 , sin π 6 , 0 ] T
ρ ˙ 10 = ρ ˙ 20 = ρ ˙ 30 = 0 r o , 22 = [ 123.3 , 75 , 40 ] T ρ ˙ d 1 = ρ ˙ d 2 = ρ ˙ d 3 = 0
r ˙ o , 1 i = r ˙ o , 2 i = 0
Table 2. Artificial potential function parameters.
Table 2. Artificial potential function parameters.
PM N i σ i ψ i
0.012 × diag ( 1 , 1 , 1 ) diag ( 0.8 , 0.6 , 0.4 ) diag ( 15 , 15 , 15 ) 10 5 × 10 3
Table 3. The controller parameters.
Table 3. The controller parameters.
p / q ε η λ β k
5 / 3 0.1 10 5 diag ( 1 , 1 , 1 ) diag ( 1 , 1 , 1 ) 10 3 × diag ( 1 , 1 , 1 )
Table 4. SMC controller parameters.
Table 4. SMC controller parameters.
k λ 1 λ 2 α
diag ( 0.5 , 0.3 , 0.4 ) diag ( 0.3 , 0.3 , 0.3 ) diag ( 0.3 , 0.3 , 0.3 ) 0.5
Table 5. Satellite simulator subsystem components.
Table 5. Satellite simulator subsystem components.
Type of ComponentComponent Name
Industrial control computer
Digital output card
Wireless network card
Control and communication componentsWireless router
Solenoid valve
Air jet nozzle
Rigid body reflecting ball
High-pressure composite gas cylinder
Cylinder terminal connector
Filters
High-pressure regulating valve
Air supply componentsHigh-pressure barometer
High-pressure gas pipeline fittings
Low-pressure regulating valve
Low-pressure barometer
Low-pressure gas pipeline fittings
Lithium ion battery
DC-DC cable
Power supply componentsPlug for connection
Air switch
Electric wire
Table 6. Artificial potential function parameters.
Table 6. Artificial potential function parameters.
PM N i σ i ψ i
0.012 × diag ( 1 , 1 , 1 ) diag ( 0.8 , 0.6 , 0.4 ) diag ( 15 , 15 , 15 ) 10 5 × 10 3
Table 7. Initial/desired position selection (mm).
Table 7. Initial/desired position selection (mm).
Initial StatePosition of ObstaclesDesired State
ρ 1 = [ 1400 , 1400 ] T r o , 1 = [ 1000 , 100 ] T ρ d 1 = [ 1818 , 1915 ] T
ρ 2 = [ 100 , 1400 ] T r o , 2 = [ 900 , 1000 ] T ρ d 2 = [ 170 , 1887 ] T
ρ 3 = [ 1100 , 1400 ] T r 1 = r 2 = 800 ρ d 3 = [ 1247 , 1836 ] T
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, J.; Wang, Y.; Jia, J.; Chi, C.; Li, H. Collision Avoidance Second Order Sliding Mode Control of Satellite Formation with Air-Floated Platform Semi-Physical Simulation. Electronics 2023, 12, 3179. https://doi.org/10.3390/electronics12143179

AMA Style

Zhang J, Wang Y, Jia J, Chi C, Li H. Collision Avoidance Second Order Sliding Mode Control of Satellite Formation with Air-Floated Platform Semi-Physical Simulation. Electronics. 2023; 12(14):3179. https://doi.org/10.3390/electronics12143179

Chicago/Turabian Style

Zhang, Ji, Yili Wang, Jun Jia, Chuanguo Chi, and Huayi Li. 2023. "Collision Avoidance Second Order Sliding Mode Control of Satellite Formation with Air-Floated Platform Semi-Physical Simulation" Electronics 12, no. 14: 3179. https://doi.org/10.3390/electronics12143179

APA Style

Zhang, J., Wang, Y., Jia, J., Chi, C., & Li, H. (2023). Collision Avoidance Second Order Sliding Mode Control of Satellite Formation with Air-Floated Platform Semi-Physical Simulation. Electronics, 12(14), 3179. https://doi.org/10.3390/electronics12143179

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop