CN112874818B - Finite time state feedback control method of spacecraft rendezvous system - Google Patents
Finite time state feedback control method of spacecraft rendezvous system Download PDFInfo
- Publication number
- CN112874818B CN112874818B CN202110068888.2A CN202110068888A CN112874818B CN 112874818 B CN112874818 B CN 112874818B CN 202110068888 A CN202110068888 A CN 202110068888A CN 112874818 B CN112874818 B CN 112874818B
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- time
- following
- space model
- finite time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 11
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 7
- 230000005484 gravity Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000014759 maintenance of location Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/242—Orbits and trajectories
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a finite time state feedback control method of a spacecraft rendezvous system. The invention designs a finite time state feedback controller aiming at a spacecraft rendezvous system with saturated actuators based on low-gain feedback and event trigger control, saves the computing resources of the spacecraft rendezvous system and enables two spacecrafts to complete rendezvous tasks in finite time.
Description
Technical Field
The invention belongs to the technical field of space control, and designs a finite time state feedback control method based on event triggering aiming at a spacecraft rendezvous system. The spacecraft rendezvous system with the saturated actuator is controlled within limited time, so that the purpose that the spacecraft rendezvous in limited time is achieved.
Background
Spacecraft rendezvous is an important technology for current and future space flight missions. Therefore, it is very important to design an effective control method for a spacecraft rendezvous system so as to realize that two spacecrafts complete rendezvous tasks in a limited time.
The spacecraft orbit rendezvous is that a target spacecraft running on a space orbit meets another tracking spacecraft which tracks the target spacecraft, and the tracking spacecraft finally rendezvouss with the target spacecraft by adjusting the running orbit of the tracking spacecraft. In consideration of the real situation, since the acceleration generated by the thruster of the spacecraft is limited, if the actually designed controller does not consider the upper limit of the acceleration, the instability of the orbit intersection system can be caused, and finally the intersection failure of the spacecraft can be caused. Moreover, due to the development of the spacecraft technology, the realization that two spacecrafts complete the rendezvous mission in limited time also becomes an important index in spacecraft rendezvous, so that the limited time control of a spacecraft rendezvous system is particularly important.
With the development of modern networked control systems, in order to reduce sampling update of the system in data transmission and reduce the data transmission amount and data calculation amount of the system, an event trigger mechanism is introduced to achieve the purpose of saving system resources. Therefore, it is necessary to design an effective method for meeting two spacecrafts in a limited time and saving computing resources for a spacecraft meeting system with saturated actuators.
Disclosure of Invention
The invention provides a finite time state feedback control method based on event triggering to realize finite time control of a spacecraft rendezvous system.
The invention designs a finite time state feedback controller based on low gain feedback and event trigger control in consideration of the influence of the saturation of an actuator of a spacecraft rendezvous system. The invention establishes a relative motion equation of the spacecraft orbit intersection system with actuator saturation, and the designed controller realizes the effective control of the spacecraft intersection system.
The method comprises the following specific steps:
step 1, establishing a relative motion equation of a spacecraft rendezvous system
The target spacecraft is assumed to be running on a circular orbit with radius R and the relative distance between the two is R. And establishing a coordinate system o-xyz of the orbit of the target spacecraft, wherein a coordinate axis x is the direction of the radius R of the circular orbit, a coordinate axis y is the running direction of the tracking spacecraft, the coordinate axis z is perpendicular to the plane of the target spacecraft moving relative to the earth mass center, the direction and the coordinate axes x and y form a right-hand coordinate system, and the origin o is the mass center of the target spacecraft.
Let the gravity constant μ be GM, where M is the mass of the surrounded planet and G is the gravity constant. Can calculate the orbit angular velocity of the target spacecraftThe relative motion equation between the target spacecraft and the tracking spacecraft can be deduced according to the Newton's theory of motion:
wherein,and x, y and z respectively represent the relative distance between the tracking spacecraft and the target spacecraft in the directions of the x axis, the y axis and the z axis. a is x ,a y ,a z Acceleration components in the directions of coordinate axes x, y, and z, respectively. Omega x ,ω y ,ω z The maximum accelerations in the directions of the three coordinate axes are respectively. sat (-) represents a unit saturation function.
Step 2, establishing a spacecraft rendezvous system state space model
Taylor expansion of γ at the origin and retention to the first order differential term yields a linearized differential equation:
by selecting a state vector
Obtaining a state space model
step 3, designing time-varying parameters
The design time varying parameter ξ (t) is as follows,
θ c =θ c (ξ 0 ) ≧ 1 is a constant. And a scalar theta c (ξ 0 ) Can be obtained by the following formula
θ c (ξ 0 )=6ξ 0 λ max (U(ξ 0 )W(ξ 0 ) -1 ),
Wherein W (ξ) 0 ) And U (xi) 0 ) Are unique symmetric positive solutions of the following Lyapunov equation, respectively
Step 4, designing event trigger conditions
The event triggering conditions are designed as follows:
wherein,
alpha is an event triggering parameter. 0 < delta < 1/6 is an alternative positive scalar quantity such that 0 < alpha < 1
The smaller δ, the higher the sampling frequency of the system. And time t e [ t ∈ [ [ t ] k ,t k+1 ) K ∈ N is the time of event trigger, and N represents a natural number set. Xi 0 =ξ(t 0 ),ξ Δ =ξ(T Δ ),t 0 Indicating the initial time of the system, T Δ Indicating the moment when the system is stable.Is the variable of the error, and is,
represents the current state X (t) and the last sampling state X (t) of the system k ) The error between.
Step 5, designing a finite time state feedback controller
The finite time state feedback controller is designed as follows
U(t)=-B T P(ξ(t))X(t k ),t∈[t k ,t k+1 ),
P(ξ(t))∈R 6×6 Is a solution of the following parametric Lyapunov equation
A T P(ξ(t))+P(ξ(t))A-P(ξ(t))BB T P(ξ(t))=-ξ(t)P(ξ(t))
Step 6, designing an ellipsoid set
Defining two sets
ε(t)={X:6ξ(t)X T P(ξ(t))X≤1},
And
|' represents a 2-norm of the matrix or vector, and epsilon (t) is a set of ellipsoids. When X (t) belongs to the setWhen the actuator is not saturated.
The calculation results in that,
‖B T P(ξ(t))X‖ 2 =X T P(ξ(t))BB T P(ξ(t))X≤6ξ(t)X T P(ξ(t))X,
that is, when X (t) k ) E ε (t), the actuator will not saturate, i.e.
sat(B T P(ξ(t))X(t k ))=B T P(ξ(t))X(t k )
Step 7, establishing a closed loop system state space model
Substituting the designed finite time state feedback controller into a state space model of a spacecraft rendezvous system to obtain the following closed-loop system state space model
Considering for arbitraryThe actuator does not saturate. Further simplified to obtain the following closed-loop system state space model
Step 8, stability analysis of closed loop system
According to the Lyapunov stability theory, the following Lyapunov function is selected
V(X,t)=6ξ(t)X T P(ξ(t))X
V (X, t) vs. time t e [ t ] k ,t k+1 ) Is a derivative of
Substituting the event trigger condition can further obtain
That is to say that the temperature of the molten steel,
substituting the designed time-varying parameter xi (t) into the above formula to obtain
This also indicates that the closed loop system is stable for a finite time T.
Aiming at a spacecraft rendezvous system with actuator saturation, the invention designs a finite time state feedback controller based on low-gain feedback and event triggering conditions, thereby avoiding the occurrence of actuator saturation, saving the computing resources of the system and realizing that two spacecrafts complete rendezvous tasks in finite time.
Detailed Description
Step 1, establishing a relative motion equation of a spacecraft rendezvous system
The target spacecraft is assumed to run on a circular orbit with radius R, and the relative distance between the two orbits is R. Establishing a coordinate system o-xyz of the orbit of the target spacecraft, wherein the x direction of a coordinate axis is the direction of the radius R of a circular orbit, the y direction of the coordinate axis is the direction of tracking the running of the spacecraft, the z direction of the coordinate axis is perpendicular to the plane of the target spacecraft moving relative to the earth mass center, the direction and the coordinate axes x and y form a right-hand coordinate system, and an origin o is the mass center of the target spacecraft.
Let the gravity constant μ be GM, where M is the mass of the surrounded planet and G is the gravity constant. The orbit angular velocity of the target spacecraft can be calculatedThe relative motion equation between the target spacecraft and the tracking spacecraft can be deduced according to the Newton's theory of motion:
wherein,x, y, z represent the relative distance of the tracking spacecraft from the target spacecraft in the x, y, z directions, respectively. a is a x ,a y ,a z Acceleration components in the directions of coordinate axes x, y, and z, respectively. Omega x ,ω y ,ω z The maximum accelerations in the directions of the three coordinate axes are respectively. sat (-) represents a unit saturation function.
Step 2, establishing a spacecraft rendezvous system state space model
Taylor expansion of γ at the origin and retention to the first order derivative term can be obtained
Linearized differential equation:
by selecting a state vector
Obtaining a state space model
step 3, designing time-varying parameters
The design time varying parameter ξ (t) is as follows,
θ c =θ c (ξ 0 ) ≧ 1 is a constant. And a scalar theta c (ξ 0 ) Can be obtained by the following formula
θ c (ξ 0 )=6ξ 0 λ max (U(ξ 0 )W(ξ 0 ) -1 ),
Wherein W (ξ) 0 ) And U (xi) 0 ) Are unique symmetric positive solutions of the following Lyapunov equation, respectively
Step 4, designing event trigger conditions
The event triggering conditions are designed as follows:
wherein,
alpha is an event triggering parameter. 0 < delta < 1/6 is an alternative positive scalar quantity such that 0 < alpha < 1
The smaller δ, the higher the sampling frequency of the system. And time t e [ t ∈ [ [ t ] k ,t k+1 ) And k ∈ N is the time of the event trigger. Xi 0 =ξ(t 0 ),ξ Δ =ξ(T Δ ),t 0 Indicating the initial time of the system, T Δ Indicating the moment when the system is stable.Is the variable of the error, and is,
represents the current state X (t) and the last sampling state X (t) of the system k ) The error between.
Step 5, designing a finite time state feedback controller
A finite time state feedback controller is designed,
U(t)=-B T P(ξ(t))X(t k ),t∈[t k ,t k+1 ),
P(ξ(t))∈R 6×6 is a solution of the following parametric Lyapunov equation
A T P(ξ(t))+P(ξ(t))A-P(ξ(t))BB T P(ξ(t))=-ξ(t)P(ξ(t))
Step 6, designing an ellipsoid set
Defining two sets
ε(t)={X:6ξ(t)X T P(ξ(t))X≤1},
And
II denotes a 2-norm matrix or vectorAnd ε (t) is a set of ellipsoids. When X (t) belongs to the setWhen this happens, the actuator is not saturated.
The calculation results show that the method has the advantages of high efficiency,
‖B T P(ξ(t))X‖ 2 =X T P(ξ(t))BB T P(ξ(t))X≤6ξ(t)X T P(ξ(t))X,
that is, when X (t) k ) E ε (t), the actuator will not saturate, i.e.
sat(B T P(ξ(t))X(t k ))=B T P(ξ(t))X(t k )
Step 7, establishing a closed loop system state space model
Substituting the designed finite time state feedback controller into a state space model of a spacecraft rendezvous system to obtain the following closed-loop system state space model
Considering for arbitraryThe actuator does not saturate. Further simplified to obtain the following closed-loop system state space model
Step 8, stability analysis of closed loop system
According to the Lyapunov stability theory, the following Lyapunov function is selected
V(X,t)=6ξ(t)X T P(ξ(t))X
V (X, t) versus time t ∈ [ t ] k ,t k+1 ) Is a derivative of
Substituting the event trigger condition can further obtain
substituting the designed time-varying parameter xi (t) into the above formula to obtain
This also indicates that the closed loop system is stable for a finite time T.
Claims (1)
1. A finite time state feedback control method of a spacecraft rendezvous system is characterized by comprising the following steps:
the method comprises the following steps: establishing a relative motion equation of a spacecraft rendezvous system
Assuming that the target spacecraft runs on a circular orbit with the radius of R, and the relative distance between the target spacecraft and the circular orbit is R; establishing a coordinate system o-xyz of the orbit of the target spacecraft, wherein the x direction of a coordinate axis is the direction of the radius R of a circular orbit, the y direction of the coordinate axis is the direction of tracking the running of the spacecraft, the z direction of the coordinate axis is perpendicular to the plane of the target spacecraft moving relative to the earth mass center, the direction and the coordinate axes x and y form a right-hand coordinate system, and an origin o is the mass center of the target spacecraft;
making the gravity constant mu equal to GM, wherein M is the mass of a surrounded planet, and G is a universal gravity constant; determining the orbital angular velocity of a target spacecraftThe relative motion equation between the target spacecraft and the tracking spacecraft can be deduced according to the Newton's theory of motion:
wherein,x, y and z respectively represent the relative distance between the tracking spacecraft and the target spacecraft in the x, y and z directions; a is x ,a y ,a z Acceleration components in the directions of coordinate axes x, y and z are respectively; omega x ,ω y ,ω z The maximum acceleration in the directions of three coordinate axes are respectively; sat (·) represents a unit saturation function;
step two: establishing a spacecraft rendezvous system state space model
Taylor expansion of γ at the origin and retention to the first order differential term yields a linearized differential equation:
by selecting a state vector
Obtaining a state space model
step three: time-varying parametric design
The design time-varying parameter ξ (t) is as follows,
θ c =θ c (ξ 0 ) 1 or more is a constant; and a scalar theta c (ξ 0 ) Can be obtained by the following formula
θ c (ξ 0 )=6ξ 0 λ max (U(ξ 0 )W(ξ 0 ) -1 ),
Wherein W (ξ) 0 ) And U (xi) 0 ) Are unique symmetric positive solutions of the following Lyapunov equation, respectively
Step four: event trigger condition design
The event triggering conditions are designed as follows:
wherein,
α is an event triggering parameter; 0 < delta < 1/6 is an alternative positive scalar such that 0 < alpha < 1
The smaller the delta, the higher the sampling frequency of the system; and time t e [ t ∈ [ [ t ] k ,t k+1 ) K belongs to N and is the time of event triggering, and N represents a natural number set; xi 0 =ξ(t 0 ),ξ Δ =ξ(T Δ ),t 0 Indicating the initial time of the system, T Δ Indicating the moment when the system is stable;is the variable of the error, and is,
represents the current state X (t) and the last sampling state X (t) of the system k ) An error therebetween;
step five: finite time state feedback controller design
A finite-time state feedback controller for controlling the output of the inverter,
U(t)=-B T P(ξ(t))X(t k ),t∈[t k ,t k+1 ),
P(ξ(t))∈R 6×6 is a solution of the following parametric Lyapunov equation
A T P(ξ(t))+P(ξ(t))A-P(ξ(t))BB T P(ξ(t))=-ξ(t)P(ξ(t))
Step six: design ellipsoid set
Defining two sets
ε(t)={X:6ξ(t)X T P(ξ(t))X≤1},
And
|' represents a 2-norm of the matrix or vector, epsilon (t) is a set of ellipsoids; when X (t) belongs to the setWhen the actuator is not saturated;
the calculation results in that,
‖B T P(ξ(t))X‖ 2 =X T P(ξ(t))BB T P(ξ(t))X≤6ξ(t)X T P(ξ(t))X,
that is, when X (t) k ) E ε (t), the actuator will not saturate, i.e.
sat(B T P(ξ(t))X(t k ))=B T P(ξ(t))X(t k )
Step seven: establishing a closed-loop system state space model
Substituting the designed finite time state feedback controller into a state space model of a spacecraft rendezvous system to obtain the following closed-loop system state space model
Considering for arbitraryThe actuator is not saturated; further simplified to obtain the following closed-loop system state space model
Step eight: stability analysis of closed loop systems
According to the Lyapunov stability theory, the following Lyapunov function is selected
V(X,t)=6ξ(t)X T P(ξ(t))X
V (X, t) versus time t ∈ [ t ] k ,t k+1 ) Is a derivative of
Substituting the event trigger condition can further obtain
substituting the designed time-varying parameter xi (t) into the above formula to obtain
This also indicates that the closed loop system is stable for a finite time T.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110068888.2A CN112874818B (en) | 2021-01-19 | 2021-01-19 | Finite time state feedback control method of spacecraft rendezvous system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110068888.2A CN112874818B (en) | 2021-01-19 | 2021-01-19 | Finite time state feedback control method of spacecraft rendezvous system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112874818A CN112874818A (en) | 2021-06-01 |
CN112874818B true CN112874818B (en) | 2022-08-26 |
Family
ID=76049742
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110068888.2A Active CN112874818B (en) | 2021-01-19 | 2021-01-19 | Finite time state feedback control method of spacecraft rendezvous system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112874818B (en) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105843077B (en) * | 2016-04-07 | 2019-03-29 | 杭州电子科技大学 | A kind of gain switching method of Spacecraft Rendezvous system and maximum domain of attraction estimation |
CN105974822B (en) * | 2016-06-13 | 2019-02-22 | 北京航空航天大学 | A kind of spacecraft, which is independently diversion, intersects the verification method of control system ground validation device |
CN106407619B (en) * | 2016-11-16 | 2019-08-02 | 哈尔滨工业大学 | A kind of linear feedback Stabilization method of Control constraints Spacecraft Rendezvous control system |
CN110262280B (en) * | 2019-02-26 | 2020-07-14 | 北京控制工程研究所 | Intelligent data analysis and decision support system for spacecraft rendezvous and docking flight control |
CN110727199B (en) * | 2019-11-18 | 2022-09-30 | 哈尔滨工业大学 | Time-varying feedback limited time stabilization method for controlling restricted spacecraft rendezvous control system |
-
2021
- 2021-01-19 CN CN202110068888.2A patent/CN112874818B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN112874818A (en) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Adaptive backstepping control of spacecraft rendezvous and proximity operations with input saturation and full-state constraint | |
CN104527994B (en) | Multi-polar cross-over becomes the track set time soon and holds position sensing tracking and controlling method | |
McFarland et al. | Adaptive nonlinear control of agile antiair missiles using neural networks | |
Qinglei et al. | Adaptive integral-type sliding mode control for spacecraft attitude maneuvering under actuator stuck failures | |
JPH0789499A (en) | Method and device for controlling inclined orbit attitude of momentum bias space ship | |
CN103076807B (en) | A kind of control method of drive lacking flexible spacecraft attitude stabilization | |
Iyer et al. | Variable structure slewing control and vibration damping of flexible spacecraft | |
Zha et al. | Quaternion-based nonlinear trajectory tracking control of a quadrotor unmanned aerial vehicle | |
Jin et al. | Underactuated spacecraft angular velocity stabilization and three-axis attitude stabilization using two single gimbal control moment gyros | |
CN112286058A (en) | Time-varying feedback limited time control method for controlling restricted satellite formation flight system | |
Kojima et al. | Model predictive steering control law for double gimbal scissored-pair control moment gyros | |
Fawzy Ahmed et al. | Design of gain schedule fractional PID control for nonlinear thrust vector control missile with uncertainty | |
Bhattacharjee et al. | Robust control strategy for quadcopters using sliding mode control and model predictive control | |
Leomanni et al. | An MPC-based attitude control system for all-electric spacecraft with on/off actuators | |
Özdoğan et al. | Design, modeling, and control allocation of a heavy-lift aerial vehicle consisting of large fixed rotors and small tiltrotors | |
Kojima et al. | Steering control law for double-gimbal scissored-pair CMG | |
CN112874818B (en) | Finite time state feedback control method of spacecraft rendezvous system | |
Inumoh et al. | Bounded gain-scheduled LQR satellite control using a tilted wheel | |
CN112882390B (en) | Finite time output feedback control method for spacecraft rendezvous system | |
Kumar et al. | Exponential reaching law based robust trajectory tracking for unmanned aerial vehicles | |
Xin et al. | Spacecraft position and attitude control with Theta-D technique | |
Wang et al. | Review on guidance and control of aerospace vehicles: recent progress and prospect | |
Dodds et al. | Sliding-mode control system for the three-axis attitude control of rigid-body spacecraft with unknown dynamics parameters | |
Guo et al. | Integrated power, attitude, and vibration control of gyroelastic body | |
Loubar et al. | Altitude back-stepping control of quadcopter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |