CN108459610B - Method for inhibiting liquid sloshing during lander power descent - Google Patents
Method for inhibiting liquid sloshing during lander power descent Download PDFInfo
- Publication number
- CN108459610B CN108459610B CN201810164238.6A CN201810164238A CN108459610B CN 108459610 B CN108459610 B CN 108459610B CN 201810164238 A CN201810164238 A CN 201810164238A CN 108459610 B CN108459610 B CN 108459610B
- Authority
- CN
- China
- Prior art keywords
- moment
- angular velocity
- disturbance
- control period
- difference
- 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 29
- 239000007788 liquid Substances 0.000 title claims abstract description 22
- 230000002401 inhibitory effect Effects 0.000 title claims description 4
- 230000001629 suppression Effects 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000005096 rolling process Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims 2
- 239000003380 propellant Substances 0.000 description 4
- 230000005764 inhibitory process Effects 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000004088 simulation Methods 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/244—Spacecraft control systems
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (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)
Abstract
The invention relates to a lander power descending attitude control method, in particular to a liquid shaking suppression attitude control method during lander power descending, and belongs to the technical field of spacecraft attitude control. Aiming at the problem of suppressing liquid shaking in a power descending section, the method adopts two observers with different filter coefficients to respectively estimate slowly-changing engine interference torque and the sum of the engine interference torque and the shaking interference torque. According to the historical values of the two estimated disturbance moments for a period of time, the state of the disturbance moment and the estimated state of the observer are judged, one estimated value of the disturbance moment is selected, attitude control feedforward compensation is carried out, and attitude control precision is improved.
Description
Technical Field
The invention relates to a lander power descending attitude control method, in particular to a liquid shaking suppression attitude control method during lander power descending, and belongs to the technical field of spacecraft attitude control.
Background
During the power descending of the detector, the guidance system gives a target attitude according to the position and speed information of the detector, the attitude control system determines the given current attitude according to the attitude, the detector is controlled to tend to the target attitude, and finally the detector is in soft landing on the surface of the moon.
When the lander propellant storage tank is a surface tension storage tank, in the landing process, besides the interference torque generated by the engine, the liquid propellant is easy to shake to generate shaking torque, so that the landing safety is endangered.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method overcomes the defects of the prior art, provides a method for inhibiting the liquid shaking in the power descent aiming at the attitude control of the lander with the liquid propellant in the power descent period, and improves the attitude control precision.
The technical solution of the invention is as follows:
a method for restraining the dynamic descent liquid shaking of a lander can estimate the angular velocity and the disturbance moment of the lander based on a Longberger observer method, and the angular velocity coefficient and the disturbance moment coefficient of the Longberger observer are both 0-1;
when the disturbance moment coefficient of the Robert observer is larger than zero and smaller than 0.5, the estimation result is slowly-changed engine disturbance moment and corresponding angular velocity, and is defined as slow-to-interference torque and slow-to-interference angular velocity, and the coefficient of the Robert observer at the moment is defined as slow-to-interference torque coefficient and slow-to-interference torque coefficient;
when the disturbance moment coefficient of the humper observer is more than or equal to 0.5 and less than 1, the estimation result is the sum of the engine disturbance moment and the shaking disturbance moment and the corresponding angular velocity, and the sum is defined as a fast-changing disturbance moment and a fast-changing angular velocity, and the coefficient of the humper observer at the moment is defined as a fast-changing angular velocity coefficient and a fast-changing disturbance moment coefficient;
according to the statistical characteristic of the difference value of the fast-changing disturbance moment and the slow-changing disturbance moment in a period of time, determining to adopt the fast-changing disturbance moment or the slow-changing disturbance moment as a feedforward compensation moment for attitude control;
the method is based on the estimation of the angular speed and the disturbance moment of the lander by a Romberg observer, and comprises the following specific steps:
(1) setting the initial state of the lander, the angular velocity ω is changed rapidly_e1(0)=[0 0 0]TQuickly changing disturbance moment Md_e1(0)=[0 0 0]T(ii) a Slowly changing angular velocity omega_e2(0)=[0 0 0]TSlow-to-break disturbance torque Md_e2(0)=[0 0 0]T(ii) a Initial attitude control moment Mc(0)=[0 0 0]T;
(2) Estimating the fast changing angular speed omega in the ith (i is more than or equal to 1) control period_e1(i),ω_e1(i)=ω_e1(i-1)+J-1·(Mc(i-1)+Md_e1(i-1)-cross(ω_e1(i-1),J·ω_e1(i-1)))·Δt+L11·(ω(i)-ω_e1(i-1)-J-1·(Mc(i-1)+Md_e1(i-1)-cross(ω_e1(i-1),J·ω_e1(i-1)))·Δt);
Wherein: omega_e1(i-1) fast upward-periodic angular velocity, Md_e1(i-1) is the upper-cycle fast-varying disturbance moment, L11For fast changing angular velocity coefficient, Mc(i-1) is an upper period attitude control moment, omega (i) is an ith control period, an attitude angular velocity is obtained through measurement, J is a lander inertia matrix, and delta t is a control period;
(3) estimating the fast-changing disturbance moment M in the ith control periodd_e1(i),Md_e1(i)=Md_e1(i-1)+L21·J·(ω(i)-ω_e1(i-1)-J-1·(Mc(i-1)+Md_e1(i-1)-cross(ω_e1(i-1),J·ω_e1(i-1)))·Δt)/Δt;
Wherein: l is21The interference moment coefficient is changed rapidly;
(4) estimating the slow change angular velocity omega in the ith control period_e2(i),ω_e2(i)=ω_e2(i-1)+J-1·(Mc(i-1)+Md_e2(i-1)-cross(ω_e2(i-1),J·ω_e2(i-1)))·Δt+L12·(ω(i)-ω_e2(i-1)-J-1·(Mc(i-1)+Md_e2(i-1)-cross(ω_e2(i-1),J·ω_e2(i-1)))·Δt);
Wherein: omega_e2(i-1) upper-period slow variation of angular velocity, Md_e2(i-1) Upper-cycle Slow Dry disturbance Torque, L12Is a slow varying angular velocity coefficient;
(5) estimating the slow-to-interference torque M in the ith control periodd_e2(i),Md_e2(i)=Md_e2(i-1)+L22·J·(ω(i)-ω_e2(i-1)-J-1·(Mc(i-1)+Md_e2(i-1)-cross(ω_e2(i-1),J·ω_e2(i-1)))·Δt)/Δt;
Wherein: l is22Is the slow-to-break disturbance torque coefficient;
(6) ith controlIn a period, obtaining a torque difference historical value matrix dMdeD matrix dMdeIs 3 lines NwMatrix of columns, NwCounting the number of cycles; the method comprises the step of obtaining M in the step (3)d_e1(i) Subtracting M of step (5)d_e2(i) Difference value M ofd_e1(i)-Md_e2(i) Assigning to a torque difference historical value matrix dMdeMod (i, N) ofw) Columns, wherein: the function mod () is a remainder function;
(7) in the ith control cycle, if i<NwAnd then: feedforward compensation moment DeltaMd(i)=[0 0 0]T;
(8) In the ith control period, if i is more than or equal to NwAnd then: calculating the mean value of the rolling moment difference Meanx (i),the rolling moment difference variance stdx (i) is calculated,wherein: j is 1 to Nw(ii) a If | dMde(1,i)-Meanx(i)|>stdx (i), then: will Md_e2(i) Is assigned to Δ Md(i) The 1 st element of (a); if | dMde(1, i) -Meanx (i) < stdx (i), then: will Md_e1(i) Is assigned to Δ Md(i) The 1 st element of (a);
(9) in the ith control period, if i is more than or equal to NwAnd then: calculating the mean value of the difference between the pitching moments mean (i),the difference variance stdy (i) of the pitching moment differences is calculated,wherein: j is 1 to Nw(ii) a If | dMde(2,i)-Meany(i)|>stdy (i), then: will Md_e2(i) Is assigned to Δ Md(i) The 2 nd element of (1); if | dMde(2, i) -means (i) < stdy (i), then: will Md_e1(i) Is assigned to Δ Md(i) The 2 nd element of (1);
(10) in the ith control period, if i is more than or equal to NwAnd then: calculating the mean value of the yaw moment difference Meanz (i),calculating the yaw moment difference variance stdz (i),wherein: j is 1 to Nw(ii) a If | dMde(3,i)-Meanz(i)|>stdz (i), then: will Md_e2(i) Is assigned to Δ Md(i) The 3 rd element of (a); if | dMde(3, i) -Meanz (i) ≦ stdz (i), then: will Md_e1(i) Is assigned to Δ Md(i) The 3 rd element of (a);
(11) in the ith control period, calculating attitude control moment Mc(i),Mc(i)=TPID(i)+ΔMd(i) Wherein: t isPID(i) The torque is controlled by PID.
Compared with the prior art, the invention has the beneficial effects that:
aiming at the problem of liquid shaking inhibition in a power descending section, the attitude control method adopts a humper observer with different filter coefficients to respectively estimate slowly-changing engine interference torque and the sum of the engine interference torque and shaking interference torque. According to the historical values of the two estimated disturbance moments for a period of time, the state of the disturbance moment and the estimated state of the observer are judged, one estimated value of the disturbance moment is selected, attitude control feedforward compensation is carried out, and attitude control precision is improved. The invention discloses a lander dynamic descent attitude control method, which aims at the problem of liquid sloshing inhibition in a dynamic descent section and adopts two observers with different filter coefficients to respectively estimate slowly-changing engine interference torque and the sum of the engine interference torque and the sloshing interference torque. According to the historical values of the two estimated disturbance moments for a period of time, the state of the disturbance moment and the estimated state of the observer are judged, one estimated value of the disturbance moment is selected, attitude control feedforward compensation is carried out, and attitude control precision is improved.
Detailed Description
The invention has the following implementation steps:
(1) setting the initial state of the lander, the angular velocity ω is changed rapidly_e1(0)=[0 0 0]TQuickly changing disturbance moment Md_e1(0)=[0 0 0]T(ii) a Slowly changing angular velocity omega_e2(0)=[0 0 0]TSlow-to-break disturbance torque Md_e2(0)=[0 0 0]T(ii) a Initial attitude control moment Mc(0)=[0 0 0]T;
(2) Estimating the fast changing angular speed omega in the ith (i is more than or equal to 1) control period_e1(i),ω_e1(i)=ω_e1(i-1)+J-1·(Mc(i-1)+Md_e1(i-1)-cross(ω_e1(i-1),J·ω_e1(i-1)))·Δt+L11·(ω(i)-ω_e1(i-1)-J-1·(Mc(i-1)+Md_e1(i-1)-cross(ω_e1(i-1),J·ω_e1(i-1)))·Δt);
Wherein: omega_e1(i-1) fast upward-periodic angular velocity, Md_e1(i-1) is the upper-cycle fast-varying disturbance moment, L11For fast changing angular velocity coefficient, Mc(i-1) is an upper period attitude control moment, omega (i) is an ith control period, an attitude angular velocity is obtained through measurement, J is a lander inertia matrix, and delta t is a control period;
(3) estimating the fast-changing disturbance moment M in the ith control periodd_e1(i),Md_e1(i)=Md_e1(i-1)+L21·J·(ω(i)-ω_e1(i-1)-J-1·(Mc(i-1)+Md_e1(i-1)-cross(ω_e1(i-1),J·ω_e1(i-1)))·Δt)/Δt;
Wherein: l is21The interference moment coefficient is changed rapidly;
(4) estimating the slow change angular velocity omega in the ith control period_e2(i),ω_e2(i)=ω_e2(i-1)+J-1·(Mc(i-1)+Md_e2(i-1)-cross(ω_e2(i-1),J·ω_e2(i-1)))·Δt+L12·(ω(i)-ω_e2(i-1)-J-1·(Mc(i-1)+Md_e2(i-1)-cross(ω_e2(i-1),J·ω_e2(i-1)))·Δt);
Wherein: omega_e2(i-1) upper-period slow variation of angular velocity, Md_e2(i-1) Upper-cycle Slow Dry disturbance Torque, L12Is a slow varying angular velocity coefficient;
(5) estimating the slow-to-interference torque M in the ith control periodd_e2(i),Md_e2(i)=Md_e2(i-1)+L22·J·(ω(i)-ω_e2(i-1)-J-1·(Mc(i-1)+Md_e2(i-1)-cross(ω_e2(i-1),J·ω_e2(i-1)))·Δt)/Δt;
Wherein: l is22Is the slow-to-break disturbance torque coefficient;
(6) in the ith control period, obtaining a torque difference historical value matrix dMdeD matrix dMdeIs 3 lines NwMatrix of columns, NwCounting the number of cycles; the method comprises the step of obtaining M in the step (3)d_e1(i) Subtracting M of step (5)d_e2(i) Difference value M ofd_e1(i)-Md_e2(i) Assigning to a torque difference historical value matrix dMdeMod (i, N) ofw) Columns, wherein: the function mod () is a remainder function;
(7) in the ith control cycle, if i<NwAnd then: feedforward compensation moment DeltaMd(i)=[0 0 0]T;
(8) In the ith control period, if i is more than or equal to NwAnd then: calculating the mean value of the rolling moment difference Meanx (i),the rolling moment difference variance stdx (i) is calculated,wherein: j is 1 to Nw(ii) a If | dMde(1,i)-Meanx(i)|>stdx (i), then: will Md_e2(i) Is assigned to Δ Md(i) The 1 st element of (a); if | dMde(1, i) -Meanx (i) < stdx (i), then: will Md_e1(i) Is assigned to Δ Md(i) The 1 st element of (a);
(9) the ithWhen controlling the period, if i is more than or equal to NwAnd then: calculating the mean value of the difference between the pitching moments mean (i),the difference variance stdy (i) of the pitching moment differences is calculated,wherein: j is 1 to Nw(ii) a If | dMde(2,i)-Meany(i)|>stdy (i), then: will Md_e2(i) Is assigned to Δ Md(i) The 2 nd element of (1); if | dMde(2, i) -means (i) < stdy (i), then: will Md_e1(i) Is assigned to Δ Md(i) The 2 nd element of (1);
(10) in the ith control period, if i is more than or equal to NwAnd then: calculating the mean value of the yaw moment difference Meanz (i),calculating the yaw moment difference variance stdz (i),wherein: j is 1 to Nw(ii) a If | dMde(3,i)-Meanz(i)|>stdz (i), then: will Md_e2(i) Is assigned to Δ Md(i) The 3 rd element of (a); if | dMde(3, i) -Meanz (i) ≦ stdz (i), then: will Md_e1(i) Is assigned to Δ Md(i) The 3 rd element of (a);
(11) in the ith control period, calculating attitude control moment Mc(i),Mc(i)=TPID(i)+ΔMd(i) Wherein: t isPID(i) The torque is controlled by PID.
Examples
(1) Lander inertia matrixThe control period Δ t is 0.128. Fast changing angular velocity coefficient L110.9, fast-varying disturbance moment coefficient L210.5, slow change angular velocity coefficient L12=0.9,Slow-dry disturbance torque coefficient L220.02. Counting the number of cycles Nw=50。
(2) Controlling the number of cycles to be less than NwTake the control cycle number equal to 1 as an example. Estimating fast varying angular velocity ω_e1(1)=[0.009 0.18 0.27]TQuick change interference moment Md_e1(1)=[0.17 2.62 3.93]TSlowly changing angular velocity omega_e2(1)=[0.009 0.18 0.27]TSlow-to-disturbance torque Md_e2(1)=[0.007 0.1 0.16]T(ii) a Moment difference historical value matrix dMdeColumn 1 of (2) is [ 0.1632.523.77]T(ii) a Calculating attitude control moment Mc(1)=[0.01 0.07 0.08]T。
(3) Controlling the number of cycles to be greater than or equal to NwTake the control cycle number equal to 50 as an example. Fast changing angular velocity omega_e1(50)=[0.05 0.22 0.32]TQuickly changing disturbance moment Md_e1(50)=[0.05 15 26]TSlowly changing angular velocity omega_e2(50)=[0.06 0.2 0.3]TSlow-to-break disturbance torque Md_e2(50)=[0.02 10 20]T(ii) a Moment difference historical value matrix dMdeColumn 50 of (2) is [ 0.0356]T(ii) a Rolling moment difference mean value Meanx (50) is 0.02, rolling moment difference variance stdx (50) is 0.06, pitching moment difference mean value Meany (50) is 0.2, pitching moment difference variance stdy (50) is 10, yawing moment difference mean value Meanz (50) is 0.3, yawing moment difference variance stdz (50) is 11, Δ Md(50)=[0.05 15 25]TPID control moment TPID(50)=[0.01 0.4 0.5]T. Attitude control moment Mc(50)=[0.06 15.4 25.5]T. Mathematical simulation shows that the attitude control moment acts on the lander to inhibit the liquid of the propellant of the lander from shaking, and high-precision attitude control in the process of power descent is realized. The angular velocity is given in degrees/s and the moment in Nm.
Those skilled in the art will appreciate that those matters not described in detail in the present specification are well known in the art.
Claims (10)
1. A method for inhibiting the liquid sloshing of lander power descending is characterized by comprising the following specific steps:
(1) setting the initial state of the lander, the angular velocity ω is changed rapidly_e1(0)=[0 0 0]TQuickly changing disturbance moment Md_e1(0)=[0 0 0]T(ii) a Slowly changing angular velocity omega_e2(0)=[0 0 0]TSlow-to-break disturbance torque Md_e2(0)=[0 0 0]T(ii) a Initial attitude control moment Mc(0)=[0 0 0]T;
(2) In the ith control period, a fast-changing angular velocity omega is estimated by using a Romberg observer_e1(i) A value of (d); i is more than or equal to 1;
(3) in the ith control period, a Romberg observer is used for estimating the fast-changing disturbance moment Md_e1(i) A value of (d);
(4) in the ith control period, a Romberg observer is used for estimating the slow-changing angular velocity omega_e2(i) A value of (d);
(5) in the ith control period, the slow-drying disturbance moment M is estimated by using a Romberg observerd_e2(i) A value of (d);
(6) in the ith control period, obtaining a torque difference historical value matrix dMde;
(7) In the ith control cycle, if i<NwAnd then: feedforward compensation moment DeltaMd(i)=[0 0 0]T;NwCounting the number of cycles;
(8) in the ith control period, if i is more than or equal to NwAnd then: calculating a rolling moment difference mean value Meanx (i); calculating a rolling moment difference variance stdx (i); wherein: j is 1 to Nw(ii) a If | dMde(1,i)-Meanx(i)|>stdx (i), then: will Md_e2(i) Is assigned to Δ Md(i) The 1 st element of (a); if | dMde(1, i) -Meanx (i) < stdx (i), then: will Md_e1(i) Is assigned to Δ Md(i) The 1 st element of (a);
(9) in the ith control period, if i is more than or equal to NwAnd then: calculating a mean value mean (i) of the differences of the pitching moments; calculating a pitching moment difference variance stdy (i); wherein: j is 1 to Nw(ii) a If | dMde(2,i)-Meany(i)|>stdy (i), then: will Md_e2(i) Is assigned to Δ Md(i) The 2 nd element of (1); if | dMde(2, i) -means (i) < stdy (i), then: will Md_e1(i) Is assigned to Δ Md(i) The 2 nd element of (1);
(10) in the ith control period, if i is more than or equal to NwAnd then: calculating a yaw moment difference mean value Meanz (i); calculating a yaw moment difference variance stdz (i); wherein: j is 1 to Nw(ii) a If | dMde(3,i)-Meanz(i)|>stdz (i), then: will Md_e2(i) Is assigned to Δ Md(i) The 3 rd element of (a); if | dMde(3, i) -Meanz (i) ≦ stdz (i), then: will Md_e1(i) Is assigned to Δ Md(i) The 3 rd element of (a);
(11) in the ith control period, calculating attitude control moment Mc(i)。
2. The landing gear power-down liquid sloshing suppression method according to claim 1, wherein: in the step (2), the angular velocity ω is rapidly changed_e1(i) The values of (A) are:
ω_e1(i)=ω_e1(i-1)+J-1·(Mc(i-1)+Md_e1(i-1)-cross(ω_e1(i-1),J·ω_e1(i-1)))·Δt+L11·(ω(i)-ω_e1(i-1)-J-1·(Mc(i-1)+Md_e1(i-1)-cross(ω_e1(i-1),J·ω_e1(i-1)))·Δt);
wherein: omega_e1(i-1) fast upward-periodic angular velocity, Md_e1(i-1) is the upper-cycle fast-varying disturbance moment, L11For fast changing angular velocity coefficient, McAnd (i-1) is an upper period attitude control moment, omega (i) is an ith control period, the attitude angular velocity is obtained through measurement, J is a lander inertia matrix, and delta t is a control period.
3. The landing gear power-down liquid sloshing suppression method according to claim 2, wherein: in the step (3), the disturbance moment M is changed rapidlyd_e1(i) The values of (A) are:
Md_e1(i)=Md_e1(i-1)+L21·J·(ω(i)-ω_e1(i-1)-J-1·(Mc(i-1)+Md_e1(i-1)-cross(ω_e1(i-1),J·ω_e1(i-1)))·Δt)/Δt;
wherein: l is21The fast-changing disturbance moment coefficient.
4. The landing gear power-down liquid sloshing suppression method according to claim 1, wherein: in the step (4), the angular velocity omega is slowly changed_e2(i) The values of (A) are: omega_e2(i)=ω_e2(i-1)+J-1·(Mc(i-1)+Md_e2(i-1)-cross(ω_e2(i-1),J·ω_e2(i-1)))·Δt+L12·(ω(i)-ω_e2(i-1)-J-1·(Mc(i-1)+Md_e2(i-1)-cross(ω_e2(i-1),J·ω_e2(i-1)))·Δt);
Wherein: omega_e2(i-1) upper-period slow variation of angular velocity, Md_e2(i-1) Upper-cycle Slow Dry disturbance Torque, L12Is a slowly varying angular velocity coefficient.
5. The landing gear power-down liquid sloshing suppression method according to claim 4, wherein: in the step (5), the slow-drying disturbance torque Md_e2(i) The values of (A) are:
Md_e2(i)=Md_e2(i-1)+L22·J·(ω(i)-ω_e2(i-1)-J-1·(Mc(i-1)+Md_e2(i-1)-cross(ω_e2(i-1),J·ω_e2(i-1)))·Δt)/Δt;
wherein: l is22Is the slow-to-disturb torque coefficient.
6. The landing gear power-down liquid sloshing suppression method according to claim 1, wherein: in the step (6), the historical value matrix dMdeIs 3 lines NwMatrix of columns, NwCounting the number of cycles; the method comprises the step of obtaining M in the step (3)d_e1(i) Subtracting M of step (5)d_e2(i) Difference value M ofd_e1(i)-Md_e2(i) Assigning to a torque difference historical value matrix dMdeMod (i, N) ofw) Columns, wherein: function mod () is a remainder function.
10. the landing gear power-down liquid sloshing suppression method according to claim 1, wherein: the steps are(11) In, the attitude control moment is: mc(i)=TPID(i)+ΔMd(i) Wherein: t isPID(i) The torque is controlled by PID.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810164238.6A CN108459610B (en) | 2018-02-27 | 2018-02-27 | Method for inhibiting liquid sloshing during lander power descent |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810164238.6A CN108459610B (en) | 2018-02-27 | 2018-02-27 | Method for inhibiting liquid sloshing during lander power descent |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108459610A CN108459610A (en) | 2018-08-28 |
CN108459610B true CN108459610B (en) | 2021-02-05 |
Family
ID=63217580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810164238.6A Active CN108459610B (en) | 2018-02-27 | 2018-02-27 | Method for inhibiting liquid sloshing during lander power descent |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108459610B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113311850B (en) * | 2021-04-25 | 2024-03-19 | 北京控制工程研究所 | Phase plane attitude control method for resisting slow-varying interference |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101074880A (en) * | 2007-07-24 | 2007-11-21 | 北京控制工程研究所 | Method for scanning entry safety landing area in moon detector suspension stage |
CN101122780A (en) * | 2007-09-04 | 2008-02-13 | 北京控制工程研究所 | Half-physical emulation test system for controlling and guiding, navigating and controlling soft landing for moon |
CN102981507A (en) * | 2012-11-29 | 2013-03-20 | 北京理工大学 | Control method of soft landing autonomic obstacle avoidance fixed-thruster |
CN103303495A (en) * | 2013-04-11 | 2013-09-18 | 北京控制工程研究所 | Method for estimating disturbance moment in power decreasing process |
CN103350758A (en) * | 2013-06-24 | 2013-10-16 | 北京航空航天大学 | Attitude-adjustable moon soft lander |
CN103466100A (en) * | 2013-08-23 | 2013-12-25 | 北京控制工程研究所 | Lander soft landing posture control method |
CN105346733A (en) * | 2015-08-07 | 2016-02-24 | 上海空间推进研究所 | Star lander storage box propellant management device and manufacturing method thereof |
CN105843233A (en) * | 2016-04-11 | 2016-08-10 | 哈尔滨工程大学 | Non-linear observer based autonomous underwater vehicle motion control method |
CN106628280A (en) * | 2016-11-23 | 2017-05-10 | 南京航空航天大学 | Soft landing simulation experimental device and simulation method for Spacecraft |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050242239A1 (en) * | 2004-05-03 | 2005-11-03 | Scott Mark W | Method for VTOL aircraft landing with external load |
-
2018
- 2018-02-27 CN CN201810164238.6A patent/CN108459610B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101074880A (en) * | 2007-07-24 | 2007-11-21 | 北京控制工程研究所 | Method for scanning entry safety landing area in moon detector suspension stage |
CN101122780A (en) * | 2007-09-04 | 2008-02-13 | 北京控制工程研究所 | Half-physical emulation test system for controlling and guiding, navigating and controlling soft landing for moon |
CN102981507A (en) * | 2012-11-29 | 2013-03-20 | 北京理工大学 | Control method of soft landing autonomic obstacle avoidance fixed-thruster |
CN103303495A (en) * | 2013-04-11 | 2013-09-18 | 北京控制工程研究所 | Method for estimating disturbance moment in power decreasing process |
CN103350758A (en) * | 2013-06-24 | 2013-10-16 | 北京航空航天大学 | Attitude-adjustable moon soft lander |
CN103466100A (en) * | 2013-08-23 | 2013-12-25 | 北京控制工程研究所 | Lander soft landing posture control method |
CN105346733A (en) * | 2015-08-07 | 2016-02-24 | 上海空间推进研究所 | Star lander storage box propellant management device and manufacturing method thereof |
CN105843233A (en) * | 2016-04-11 | 2016-08-10 | 哈尔滨工程大学 | Non-linear observer based autonomous underwater vehicle motion control method |
CN106628280A (en) * | 2016-11-23 | 2017-05-10 | 南京航空航天大学 | Soft landing simulation experimental device and simulation method for Spacecraft |
Non-Patent Citations (3)
Title |
---|
《In-Flight Identification of Magnetometer and Attitude Determination System for Land-survey Mini-satellite》;Yevgeny Somov 等;《IFAC-PapersOnLine》;20151231;第16-21页 * |
《月球着陆器最终下降段的制导与控制方法研究》;胡锦昌 等;《宇航学报》;20121130;第1577-1585页 * |
《空间飞行器动力学与控制研究综述》;刘付成 等;《上海航天》;20170228;第1-28页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108459610A (en) | 2018-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108710303B (en) | Spacecraft relative attitude control method containing multi-source disturbance and actuator saturation | |
CN108803649B (en) | A kind of VTOL Reusable Launch Vehicles active disturbance rejection sliding-mode control | |
CN111268176B (en) | Perturbation track four-pulse intersection rapid optimization method | |
US20180290719A1 (en) | Motion control device and motion control method for ship | |
CN105159084B (en) | A kind of manipulator nerve network control system and method with interference observer | |
CN107272639B (en) | Detection, estimation and its adjusting method of rigid spacecraft reaction wheel failure | |
CN105629734B (en) | A kind of Trajectory Tracking Control method of Near Space Flying Vehicles | |
CN111258221B (en) | Spacecraft fault-tolerant control method based on self-adaptive sliding mode theory | |
CN109739248B (en) | Ship-borne three-degree-of-freedom parallel stable platform stability control method based on ADRC | |
CN112445234A (en) | Attitude control method and device for spacecraft | |
CN108459610B (en) | Method for inhibiting liquid sloshing during lander power descent | |
WO2022095643A1 (en) | Self-adaptive iterative guidance method and device for aerospace vehicle | |
CN111061286A (en) | Method for realizing lateral overload control of aircraft by providing damping through filtering differentiation | |
CN113589689A (en) | Sliding mode controller design method based on multi-parameter adaptive neural network | |
CN108516107B (en) | Online engine thrust and specific impulse estimation method and system | |
CN113608534A (en) | Unmanned ship tracking control method and system | |
CN114326405B (en) | Neural network backstepping control method based on error training | |
CN111638643A (en) | Displacement mode drag-free control dynamics coordination condition determination method | |
CN113759718A (en) | Self-adaptive control method for airplane wing damage | |
CN107943097B (en) | Aircraft control method and device and aircraft | |
CN111007867B (en) | Hypersonic aircraft attitude control design method capable of presetting adjustment time | |
CN109445448B (en) | Self-adaptive integral sliding-mode attitude controller for wheel-controlled minisatellite | |
CN108415443B (en) | Control method for forced flight-around of non-cooperative target | |
CN111400902A (en) | Rocket debris landing area estimation method and device, electronic equipment and storage medium | |
CN114995163B (en) | Unmanned aerial vehicle immune control method |
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 |