CN105786008B - A kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior - Google Patents

A kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior Download PDF

Info

Publication number
CN105786008B
CN105786008B CN201610196190.8A CN201610196190A CN105786008B CN 105786008 B CN105786008 B CN 105786008B CN 201610196190 A CN201610196190 A CN 201610196190A CN 105786008 B CN105786008 B CN 105786008B
Authority
CN
China
Prior art keywords
mrow
msub
mtd
mover
mtr
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
Application number
CN201610196190.8A
Other languages
Chinese (zh)
Other versions
CN105786008A (en
Inventor
郭雷
张培喜
乔建忠
李文硕
徐健伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201610196190.8A priority Critical patent/CN105786008B/en
Publication of CN105786008A publication Critical patent/CN105786008A/en
Application granted granted Critical
Publication of CN105786008B publication Critical patent/CN105786008B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • G05D1/0816Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
    • G05D1/0825Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

A kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior, for flywheel saturation and frictional behavior influence Spacecraft Control precision, the spacecraft the coupled dynamical equation containing flywheel saturation and frictional behavior is built first, secondly, set up the friction interference estimator for wheel friction characteristic, again, for flexible appendage shuttle belt come interference, design flexible vibration observer;Then, for friction Interference Estimation error and flexible vibration disturbance-observer error, design anti-saturation controller is suppressed;Finally by anti-saturation controller, friction interference estimator and the control gain of flexible vibration observer is asked for, composite layered anti-interference controller is designed, the anti-interference gesture stability of spacecraft under multi-source interference effect is completed;This method can significantly improve the Spacecraft Attitude Control precision using flywheel executing agency, available for the high-precision earth observation satellite of aerospace field, the high-precision attitude control of the spacecraft such as space telescope.

Description

A kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior
Technical field
The present invention relates to a kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior, it is adaptable to needs Carry executive capability limited and need to realize the Flexible Spacecraft control system of high-precision control, belong to spacecraft attitude Control field.
Background technology
Space action is explored recently as the mankind to enrich constantly, task complexity is increasing, and demand is also increasing, It is required that spacecraft, which carries bigger solar cell tabula rasa, can provide more energy.In addition, with spacecraft Mission Operations It is more remote, requirements at the higher level also are proposed to spacecraft communication antenna, it is necessary to which the structure more stronger antenna of big, power is in order to remote Distance receives faint signal;The demand such as power supply and communication all causes spacecraft to need to carry increasing annex, from transmitting For cost and technology enforcement difficulty, carrying the quality volume of spacecraft has strict limitation, therefore above-mentioned annex is generally used Low quality, the flexible structure design of Low rigidity, substantial amounts of flexible appendage are used, when spacecraft body carries out motor-driven, flexibility knot Structure can be produced and is directed to, so that Spacecraft Attitude Control precision can be had a strong impact on, or even the last task of influence.
In addition the reliability and long-term working stability of spacecraft attitude control system are always the pass during spacecraft is developed Key technology.Flywheel is one of most important execution unit in satellite attitude control system, the long-life launched in recent years, high-precision Degree, multifunctional triaxial stabilized satellite, are not almost exceptionally used as main execution unit by the use of flywheel.But flywheel has very Salient feature, is limited to machining accuracy, can there is a certain degree of friction, so that bring flywheel to perform error, it is another Aspect, the output torque size of the flywheel in kind of actual physical system is strictly limited, therefore needs also exist for further consideration The problem of saturation and frictional behavior.Another aspect moment of friction is transferred to spacecraft body by flywheel wheel body, causes spacecraft Body is shivered, so that very big trouble can be brought to spacecraft attitude control system.Therefore, boat is completed in order to more accurate Its device gesture stability, must pull against the influence that above-mentioned two class is mainly disturbed during Spacecraft guidance and control.
Number of patent application be 201510294341.9 in propose a kind of anti-interference appearance based on reaction wheel frictional behavior State control method, but there are problems that two:(1) this article does not consider the vibration that the subsidiary flexible appendage of Spacecraft is brought Interference, can be impacted to spacecraft precision;(2) do not have the saturated characteristic for considering flywheel, can be by certain in actual use Limitation.Number of patent application be 201510303102.5 in propose a kind of Flexible Spacecraft control for flywheel low speed friction System and method for processed, wherein the method proposed equally exists Similar Problems:(1) do not have the saturated characteristic for considering that flywheel is present; (2) what is be related in method is completely different with this patent form for wheel friction observer, does not ensure that finite time is fast Rub interference in speed tracking, therefore can be inferior in precision and rapidity context of methods.
The content of the invention
The technology of the present invention solves problem:Overcome it is existing it is not enough there is provided a kind of for flywheel saturation and frictional behavior Flexible Spacecraft control method, can be to be used as main executing agency's Spacecraft system using flywheel by the use of this method Anti- flywheel saturation and friction, capability are provided, the high-precision attitude control of wheel control Spacecraft system can be realized.
The present invention technical solution be:A kind of Flexible Spacecraft control for flywheel saturation and frictional behavior Method, implementation step is as follows:
The first step, sets up the Spacecraft the coupled dynamical equation containing flywheel saturation and frictional behavior:
For saturation and frictional behavior common in flywheel executing agency, and consideration flexible appendage kinetics equation simultaneously, The flexible spacecraft dynamics system model with flywheel saturation and frictional behavior is set up, is expressed as follows:
In formula, t represents the time, and J is the rotary inertia of spacecraft;For spacecraft attitude angular acceleration, F is spacecraft Coupling matrix between posture and flexible structure, η (t) is the mode of oscillation of flexible appendage, and ω is flexible appendage mode of oscillation pair The vibration frequency answered, ξ is the damping of flexible appendage mode, Tc(t) control that the attitude controller of Space Vehicle System is resolved is represented Torque, sat (Tc(t) it is) to consider the saturation control moment after flywheel saturated characteristic, Mf(t) after for consideration wheel friction characteristic The wheel friction interference of introducing;d1(t) outer space environmental disturbances torque being subject to for spacecraft;
It can further obtain:
In formulaRepresent the vibration interference that flexible appendage shuttle belt comes, d1(t) it is space flight The outer space environmental disturbances torque that device is subject to;Flywheel saturation and frictional behavior are already had accounted in this equation, and it is flexible attached Flexible spacecraft dynamics system model with flywheel saturation and frictional behavior, is further transformed into by the interference that part shuttle belt comes State space form, the new system under state space form is represented is as follows:
X (t) is system mode in formula,θ (t) is the attitude angle of spacecraft,For the posture of spacecraft Angular speed,For sytem matrix,For control input matrix;
Second step, is characterized for vibration interference caused by flexible appendage by following interference models:
In formula, V, w (t) and W are expressed as the output matrix of vibration interference caused by flexible appendage, state variable, are System matrix, H represents the gain matrix of saturated controller output, wheel friction interference and environmental disturbances, wherein vibration annex is done The state variable disturbedWherein intermediate variable R=(1- FTI-1F)-1
Following flexible vibration observer is constructed to flexible appendage vibration interference:
Wherein v (t) is an auxiliary State Variable of flexible vibration observer,It is that flexible appendage is directed to interference d0 (t) estimate, L is the gain of interference observer to be asked;
For wheel friction interference, the moment of friction M that fly wheel system is subject tof(t) mainly the solid friction power of bearing is included The gentle dynamic resistance square of viscous friction torque that square and lubrication belt come, torque and Speed of Reaction Wheels phase that above-mentioned two class is come due to lubrication belt Close, it is comparatively small in the case of Speed of Reaction Wheels is relatively low, therefore the low speed friction characteristic of flywheel mostlys come from solid friction, flywheel Kinetic model state-space representation formula it is as follows:
In formulaFor the angular acceleration of flywheel, Ω (t) is Speed of Reaction Wheels,For the moment of friction rate of change of flywheel, Wherein D represents the damped coefficient of flywheel, Rotary Inertia of Flywheel Jw, bearing static ramp parameter β, Coulomb friction torque Mf0
For wheel friction, following wheel friction interference estimator is designed:
WhereinWithIt is the estimate of Speed of Reaction Wheels and moment of friction, parameter k1With k2, positive parameter alpha1With α2So thatWithCan be in Finite-time convergence in Ω (t) and Mf(t), it is ensured that wheel friction Tracking error can be cut down in finite time;
3rd step, design anti-saturation controller:
With reference to friction interference estimator, flexible appendage Vibration device, and consider flywheel damp constraint characteristic, enter one The following anti-saturation controller of step design:
In formula, u (t)=sat (Tc(t)) exported for anti-saturation controller, that is, consider the saturation after flywheel saturated characteristic Control moment, K is that anti-saturation controller controls gain, and due to considering flywheel saturation, the maximum output torque for setting flywheel is umax, and umax>0 so as to have:
The control moment instruction that u (t) in formula receives for flywheel, not over the maximum of flywheel, can also be protected in addition Demonstrate,prove Space Vehicle System stable, realization has spacecraft high-precision attitude under saturation and friction condition and controlled.
The advantage of the present invention compared with prior art is:
A kind of Flexible Spacecraft control method for saturation and frictional behavior being related in the present invention is mainly directed towards The Flexible Spacecraft of flywheel low speed friction is controlled to control as the Flexible Spacecraft of main executing agency using flywheel System, while considering saturation and frictional behavior that fly wheel system is present so that for the more comprehensive of executing agency's error analysis, With widely using face.In addition the friction interference observer with finite time convergence control, energy have been used for wheel friction Enough quickly and accurately tracking wheel friction interference, add the capability of fast response of Flexible Spacecraft control system, significantly Improve speed, precision and the stability of Spacecraft Attitude Control method.
Brief description of the drawings
Fig. 1 is a kind of design cycle for the Flexible Satellite Attitude control method for being directed to flywheel saturation and frictional behavior of the present invention Figure.
Embodiment
Illustrate implementing for system and method, satellite operation by taking the satellite system of a class generalized ribbon flexible appendage as an example In earth observation pattern, have high requirements to attitude control accuracy and stability tool;
As shown in figure 1, specific implementation step of the present invention is as follows:
1st, the flexible satellite the coupled dynamical equation containing flywheel saturation and frictional behavior is set up
For saturation and frictional behavior common in flywheel executing agency, and consideration flexible appendage kinetics equation simultaneously, The flexible satellite dynamic system model with flywheel saturation and frictional behavior is set up, is expressed as follows:
In formula, t represents the time, and J is the rotary inertia of satellite;For attitude of satellite angular acceleration, F be the attitude of satellite with Coupling matrix between flexible structure, η (t) is the mode of oscillation of flexible appendage, and ω shakes for flexible appendage mode of oscillation is corresponding Dynamic frequency, ξ is the damping of flexible appendage mode, Tc(t) control moment that the attitude controller of satellite system is resolved, sat are represented (Tc(t) it is) to consider the saturation control moment after flywheel saturated characteristic, Mf(t) it is flying for being introduced after consideration wheel friction characteristic Wheel friction interference;d1(t) outer space environmental disturbances torque being subject to for satellite;
It can further obtain:
In formulaRepresent the vibration interference that flexible appendage shuttle belt comes, d1(t) it is space flight The outer space environmental disturbances torque that device is subject to;Flywheel saturation and frictional behavior are already had accounted in this equation, and it is flexible attached The interference that part shuttle belt comes, is further transformed into shape by the flexible satellite dynamic system model with flywheel saturation and frictional behavior State space form, the new system under state space form is represented is as follows:
In formulaθ (t) is the attitude angle of satellite,For the attitude angular velocity of satellite, x (t) is system shape State,For sytem matrix,For control input matrix;
2nd, characterized for vibration interference caused by flexible appendage by following interference models:
In formula, V, w (t) and W are expressed as the output matrix of vibration interference caused by flexible appendage, state variable, are System matrix, H represents the gain matrix of saturated controller output, wheel friction interference and environmental disturbances, wherein vibration annex is done The state variable disturbedWherein intermediate variable R=(1- FTI-1F)-1
Following flexible vibration observer is constructed to flexible appendage vibration interference:
Wherein v (t) is an auxiliary State Variable of flexible vibration observer,It is that flexible appendage is directed to interference d0 (t) estimate, L is the gain of interference observer to be asked;
For wheel friction interference, the moment of friction M that fly wheel system is subject tof(t) mainly the solid friction power of bearing is included The gentle dynamic resistance square of viscous friction torque that square and lubrication belt come, torque and Speed of Reaction Wheels phase that above-mentioned two class is come due to lubrication belt Close, it is comparatively small in the case of Speed of Reaction Wheels is relatively low, therefore the low speed friction characteristic of flywheel mostlys come from solid friction, flywheel Kinetic model state-space representation formula it is as follows:
In formulaFor the angular acceleration of flywheel, Ω (t) is Speed of Reaction Wheels,For the moment of friction rate of change of flywheel, Wherein D represents the damped coefficient of flywheel, Rotary Inertia of Flywheel Jw, bearing static ramp parameter β, Coulomb friction torque Mf0
For wheel friction, following wheel friction interference estimator is designed:
WhereinWithIt is the estimate of Speed of Reaction Wheels and moment of friction, parameter k1With k2, positive parameter alpha1With α2So thatWithCan be in finite time convergence control in Ω (t) and Mf(t), it is ensured that wheel friction with Track error can be cut down in finite time;
3rd, anti-saturation controller is designed
With reference to friction interference estimator, flexible appendage Vibration device, and consider flywheel damp constraint characteristic, enter one The following anti-saturation controller of step design:
In formula, u (t)=sat (Tc(t)) exported for anti-saturation controller, that is, consider the saturation after flywheel saturated characteristic Control moment, K is that anti-saturation controller controls gain, and due to considering flywheel saturation, the maximum output torque for setting flywheel is umax, and umax>0 so as to have:
The control moment instruction that u (t) in formula receives for flywheel, not over the maximum of flywheel, can also be protected in addition Demonstrate,prove satellite system stable, realization has satellite high-precision gesture stability under saturation and friction condition.
The content not being described in detail in description of the invention belongs to prior art known to professional and technical personnel in the field.

Claims (1)

1. a kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior, it is characterised in that:Including following Step:The Spacecraft the coupled dynamical equation containing flywheel saturation and frictional behavior is built first, secondly, is set up for flying Take turns frictional behavior friction interference estimator, again, in Spacecraft due to for flexible appendage shuttle belt come interference, Design flexible vibration observer;Then, for friction Interference Estimation error and flexible vibration disturbance-observer error, anti-saturation is designed Controller is suppressed;Increase finally by anti-saturation controller, friction interference estimator and the control of flexible vibration observer is asked for Benefit, designs composite layered anti-interference controller, completes the anti-interference gesture stability of spacecraft under multi-source interference effect;Specific steps It is as follows:
The first step, sets up the Spacecraft the coupled dynamical equation containing flywheel saturation and frictional behavior
Spacecraft Mission Operations it is more remote, power supply and communication need all cause spacecraft to need to carry increasing attached Part, above-mentioned annex is generally using low quality, the flexible structure design of Low rigidity, and substantial amounts of flexible appendage is used, in spacecraft sheet When body carries out motor-driven, so that Spacecraft Attitude Control precision can be had a strong impact on, or even the last task of influence;Machinery is limited to add Work precision, can have a certain degree of friction, so as to bring flywheel to perform error, on the other hand, the material object of actual physical system The output torque size of flywheel is strictly limited, therefore the problem of need also exist for further considering saturation and frictional behavior;Separately One side moment of friction is transferred to spacecraft body by flywheel wheel body, causes spacecraft body to be shivered, so that can be to boat Its device attitude control system brings very big trouble;Therefore, Spacecraft Attitude Control is completed in order to more accurate, is set in spacecraft The influence that above-mentioned two class is mainly disturbed is must pull against during meter;
For saturation and frictional behavior common in flywheel executing agency, and flexible appendage kinetics equation is considered simultaneously, set up Flexible spacecraft dynamics system model with flywheel saturation and frictional behavior, is expressed as follows:
<mrow> <msub> <mo>&amp;Sigma;</mo> <mn>1</mn> </msub> <mo>:</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>J</mi> <mover> <mi>&amp;theta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>F</mi> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>s</mi> <mi>a</mi> <mi>t</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>M</mi> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mn>2</mn> <mi>&amp;xi;</mi> <mi>&amp;omega;</mi> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msup> <mi>&amp;omega;</mi> <mn>2</mn> </msup> <mi>&amp;eta;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msup> <mi>F</mi> <mi>T</mi> </msup> <mover> <mi>&amp;theta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
In formula, t represents the time, and J is the rotary inertia of spacecraft;For spacecraft attitude angular acceleration, F is spacecraft attitude Coupling matrix between flexible structure, η (t) is the mode of oscillation of flexible appendage,Mould is vibrated for the single order of flexible appendage State,For the second order vibration mode of flexible appendage, ω is the corresponding vibration frequency of flexible appendage mode of oscillation, and ξ is flexible appendage The damping of mode, Tc(t) control moment that the attitude controller of Space Vehicle System is resolved, sat (T are representedc(t) it is) winged to consider Take turns the saturation control moment after saturated characteristic, Mf(t) the wheel friction interference to be introduced after consideration wheel friction characteristic;d1(t) it is The outer space environmental disturbances torque that spacecraft is subject to;
It can further obtain:
<mrow> <mo>(</mo> <mi>J</mi> <mo>-</mo> <msup> <mi>FF</mi> <mi>T</mi> </msup> <mo>)</mo> <mover> <mi>&amp;theta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>=</mo> <msub> <mi>d</mi> <mn>0</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <mi>s</mi> <mi>a</mi> <mi>t</mi> <mo>(</mo> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> <mo>+</mo> <msub> <mi>M</mi> <mi>f</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow>
In formulaRepresent the vibration interference that flexible appendage shuttle belt comes, d1(t) for spacecraft by The outer space environmental disturbances torque arrived;Flywheel saturation and frictional behavior are already had accounted in this equation, and flexible appendage shakes The dynamic interference brought, is further transformed into state by the flexible spacecraft dynamics system model with flywheel saturation and frictional behavior Space form, the new system under state space form is represented is as follows:
<mrow> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>A</mi> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>B</mi> <mrow> <mo>(</mo> <mi>s</mi> <mi>a</mi> <mi>t</mi> <mo>(</mo> <mrow> <msub> <mi>T</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>M</mi> <mi>f</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>d</mi> <mn>0</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
X (t) is system mode in formula,θ (t) is the attitude angle of spacecraft,For the attitude angle speed of spacecraft Degree,For sytem matrix,For control input matrix;
Second step, is characterized for vibration interference caused by flexible appendage by following interference models:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>d</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>V</mi> <mi>w</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mover> <mi>w</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>W</mi> <mi>w</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>H</mi> <mrow> <mo>(</mo> <mi>s</mi> <mi>a</mi> <mi>t</mi> <mo>(</mo> <mrow> <msub> <mi>T</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>M</mi> <mi>f</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
In formula, V, w (t) and W are expressed as output matrix, state variable, the system square of vibration interference caused by flexible appendage Battle array, H represents the gain matrix of saturated controller output, wheel friction interference and environmental disturbances, wherein vibration annex interference State variableWherein intermediate variable R=(1-FTI- 1F)-1
Substantial amounts of flexible appendage is used, when spacecraft body carries out motor-driven, so that Spacecraft Attitude Control essence can be had a strong impact on Degree, or even the last task of influence, following flexible vibration observer is constructed to flexible appendage vibration interference:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>d</mi> <mo>^</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>V</mi> <mover> <mi>&amp;omega;</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mover> <mi>&amp;omega;</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>v</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>L</mi> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mover> <mi>v</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mrow> <mo>(</mo> <mi>W</mi> <mo>+</mo> <mi>L</mi> <mi>B</mi> <mi>V</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>v</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>-</mo> <mi>L</mi> <mi>x</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>+</mo> <mi>L</mi> <mi>A</mi> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mi>L</mi> <mi>B</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <msub> <mover> <mi>M</mi> <mo>^</mo> </mover> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mi>L</mi> <mi>B</mi> <mo>+</mo> <mi>H</mi> <mo>)</mo> </mrow> <mi>u</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
Wherein v (t) is an auxiliary State Variable of flexible vibration observer,It is that flexible appendage is directed to interference d0(t) estimate Evaluation, L is the gain of interference observer to be asked;
For wheel friction interference, the moment of friction M that fly wheel system is subject tof(t) mainly the solid friction torque and profit of bearing are included The gentle dynamic resistance square of viscous friction torque that slide strips are come, above-mentioned two class is related to Speed of Reaction Wheels due to the torque that lubrication belt comes, It is comparatively small in the case of Speed of Reaction Wheels is relatively low, therefore the low speed friction characteristic of flywheel mostlys come from solid friction, flywheel it is dynamic Mechanical model state-space representation formula is as follows:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mover> <mi>&amp;Omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>M</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msup> <msub> <mi>J</mi> <mi>w</mi> </msub> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>M</mi> <mi>f</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <mi>D</mi> <mi>&amp;Omega;</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;beta;</mi> <mi>&amp;Omega;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>M</mi> <mi>f</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mo>(</mo> <mrow> <mi>&amp;Omega;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>M</mi> <mrow> <mi>f</mi> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
In formulaFor the angular acceleration of flywheel, Ω (t) is Speed of Reaction Wheels,For the moment of friction rate of change of flywheel, wherein D Represent the damped coefficient of flywheel, Rotary Inertia of Flywheel Jw, bearing static ramp parameter β, Coulomb friction torque Mf0;Moment of friction is led to Cross flywheel wheel body and be transferred to spacecraft body, cause spacecraft body to be shivered,
For wheel friction, the following wheel friction interference estimator with finite time convergence control of design:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mover> <mover> <mi>&amp;Omega;</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mover> <mi>M</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msup> <msub> <mi>J</mi> <mi>w</mi> </msub> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>+</mo> <mi>D</mi> <mover> <mi>&amp;Omega;</mi> <mo>^</mo> </mover> <mo>+</mo> <msub> <mover> <mi>M</mi> <mo>^</mo> </mover> <mi>f</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>e</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <msub> <mi>sgne</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>k</mi> <mn>2</mn> </msub> <msub> <mi>e</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mi>sgn</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <msub> <mi>J</mi> <mi>w</mi> </msub> <mi>sgn</mi> <mo>(</mo> <mrow> <msub> <mi>e</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
WhereinWithIt is the estimate of Speed of Reaction Wheels and moment of friction, parameter k1And k2, positive ginseng Number α1With α2So thatWithCan be in finite time convergence control in Ω (t) and Mf(t), it is ensured that wheel friction tracking error It can be cut down in finite time;Wheel friction interference can be quickly and accurately tracked, Flexible Spacecraft control is added The capability of fast response of system processed;
3rd step, design anti-saturation controller
The output torque size of the flywheel in kind of actual physical system is strictly limited, with reference to friction interference estimator, flexibility Accessory vibration observer, and consider flywheel damp constraint characteristic, further design following anti-saturation controller:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>u</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>s</mi> <mi>a</mi> <mi>t</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>2</mn> <mi>K</mi> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>d</mi> <mo>^</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>M</mi> <mo>^</mo> </mover> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
In formula, u (t)=sat (Tc(t)) exported for anti-saturation controller, that is, consider the saturation controling power after flywheel saturated characteristic Square, K is that anti-saturation controller controls gain, due to considering flywheel saturation, and the maximum output torque for setting flywheel is umax, and umax>0 so as to have:
<mrow> <mi>u</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>&amp;le;</mo> <msub> <mi>T</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>u</mi> <mi>max</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>T</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </mtd> <mtd> <mrow> <msub> <mi>u</mi> <mi>max</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>T</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
The control moment instruction that u (t) in formula receives for flywheel, not over the maximum of flywheel, also ensure that boat in addition Its device system is stable, and realization has spacecraft high-precision attitude under saturation and friction condition and controlled.
CN201610196190.8A 2016-03-31 2016-03-31 A kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior Active CN105786008B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610196190.8A CN105786008B (en) 2016-03-31 2016-03-31 A kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610196190.8A CN105786008B (en) 2016-03-31 2016-03-31 A kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior

Publications (2)

Publication Number Publication Date
CN105786008A CN105786008A (en) 2016-07-20
CN105786008B true CN105786008B (en) 2017-09-22

Family

ID=56395443

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610196190.8A Active CN105786008B (en) 2016-03-31 2016-03-31 A kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior

Country Status (1)

Country Link
CN (1) CN105786008B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107450588B (en) * 2017-09-29 2018-06-22 北京航空航天大学 A kind of strong disturbance rejection control method of Flexible Spacecraft control system
CN109141471B (en) * 2018-06-28 2020-07-10 上海卫星工程研究所 Implementation method for screening high-stability satellite flywheel based on micro-interference torque
CN108983801B (en) * 2018-08-14 2021-05-28 北京航空航天大学 Anti-interference attitude control method for spacecraft based on dynamic characteristics of reaction flywheel
CN109613827B (en) * 2018-12-29 2021-04-02 西北工业大学 Track intersection control method for translation points with unknown relative speeds
CN109991848B (en) * 2019-03-28 2022-02-01 鲁东大学 Anti-interference control method and device
CN110502025B (en) * 2019-08-01 2020-11-10 北京航空航天大学 Spacecraft attitude control method considering reactive flywheel characteristics and power limitation
CN113031668B (en) * 2021-02-20 2022-06-03 北京控制工程研究所 Method and system for estimating disturbance moment caused by asymmetric rotation scanning load
CN113419550B (en) * 2021-06-23 2023-04-11 北京遥感设备研究所 Output torque saturation anti-interference attitude control method and system for aircraft actuating mechanism

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5201833A (en) * 1991-07-19 1993-04-13 General Electric Company Attitude control system with reaction wheel friction compensation
CN104155969A (en) * 2014-07-28 2014-11-19 北京航空航天大学 Friction inhibiting and offsetting verification system for zero momentum wheel
CN104590588A (en) * 2014-12-04 2015-05-06 哈尔滨工业大学 Flexible satellite attitude orbit coupling control method based on isolation allowance method and pulse width fusion strategy
CN104932510A (en) * 2015-06-04 2015-09-23 北京航空航天大学 Flexible spacecraft attitude control system and flexible spacecraft attitude control method in allusion to flywheel low-speed friction

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102109856B (en) * 2010-12-24 2012-07-25 中国科学院长春光学精密机械与物理研究所 System for identifying parameters of reaction flywheel
CN102122180B (en) * 2011-02-16 2012-10-17 哈尔滨工业大学 Flywheel simulator based on FPGA
CN103116357B (en) * 2013-03-14 2016-05-11 北京航空航天大学 A kind of sliding-mode control with anti-interference fault freedom

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5201833A (en) * 1991-07-19 1993-04-13 General Electric Company Attitude control system with reaction wheel friction compensation
CN104155969A (en) * 2014-07-28 2014-11-19 北京航空航天大学 Friction inhibiting and offsetting verification system for zero momentum wheel
CN104590588A (en) * 2014-12-04 2015-05-06 哈尔滨工业大学 Flexible satellite attitude orbit coupling control method based on isolation allowance method and pulse width fusion strategy
CN104932510A (en) * 2015-06-04 2015-09-23 北京航空航天大学 Flexible spacecraft attitude control system and flexible spacecraft attitude control method in allusion to flywheel low-speed friction

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李冬柏.考虑未知飞轮摩擦力矩的航天器姿态跟踪鲁棒自适应控制.《宇航学报》.2016,第37卷(第2期), *

Also Published As

Publication number Publication date
CN105786008A (en) 2016-07-20

Similar Documents

Publication Publication Date Title
CN105786008B (en) A kind of Flexible Spacecraft control method for being directed to flywheel saturation and frictional behavior
CN104932510B (en) A kind of Flexible Spacecraft control system for flywheel low speed friction and method
CN103412491B (en) A kind of Spacecraft feature axis attitude maneuver index time-varying sliding-mode control
CN104267732B (en) Flexible satellite high stability attitude control method based on frequency-domain analysis
CN106406086A (en) Large flexible spacecraft interference compensation method based on sliding mode disturbance observer
CN104898686B (en) A kind of anti-interference attitude control method based on reaction wheel rubbing characteristics
CN104656666A (en) Relative orbit design and high-precision posture pointing control method aiming at space non-cooperative target
CN106094528A (en) A kind of spatial flexible robot arm vibration suppression algorithm
CN103092208A (en) Spacecraft high-accuracy speediness attitude maneuver method based on single gimbal control moment gyro (SGCMG) and reaction wheel (RW)
CN103235597A (en) Rapid stable joint control method for attitude maneuver of spacecraft
CN109062042A (en) A kind of finite time Track In Track control method of rotor craft
CN104656447A (en) Differential geometry nonlinear control method for aircraft anti-interference attitude tracking
CN102736518A (en) Composite anti-interference controller comprising measurement and input time delay for flexible spacecraft
CN105843237A (en) Spacecraft attitude reference instruction generation method for suppressing flexible vibration
CN105607485A (en) Adaptive fault tolerance control method for flexible liquid-filled satellite attitude based on fault characteristic model
CN102663520A (en) Suboptimal solution solving method based on optimal control problem pseudo-spectral method solving framework
CN106354147A (en) Control system and control method for posture rail coupling under shifting state of big barycenter based on engine pre-swinging
CN108427272A (en) Flexible Spacecraft control based on flexible mode observation and vibration suppressing method
CN107861386A (en) A kind of anti-interference attitude control ground checking system and its control method based on angular speed observer
CN103863578B (en) Mars landing device jet thrust device and control moment gyroscope combined control system
CN105938368A (en) Flexible spacecraft attitude control method based on limited time friction estimation
Keim et al. Spacecraft inertia estimation via constrained least squares
Xu et al. Robust finite-time attitude tracking control of a CMG-based AUV with unknown disturbances and input saturation
Lu et al. Active disturbance rejection sliding mode altitude and attitude control of a quadrotor with uncertainties
Wang et al. Control system design for multi-rotor mav

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant