CN103808323A - Cosine transition acceleration path method for satellite attitude tracking maneuver - Google Patents
Cosine transition acceleration path method for satellite attitude tracking maneuver Download PDFInfo
- Publication number
- CN103808323A CN103808323A CN201210439266.7A CN201210439266A CN103808323A CN 103808323 A CN103808323 A CN 103808323A CN 201210439266 A CN201210439266 A CN 201210439266A CN 103808323 A CN103808323 A CN 103808323A
- Authority
- CN
- China
- Prior art keywords
- path
- satellite
- attitude
- motor
- acceleration
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
The invention provides a cosine transition acceleration path method for satellite attitude tracking maneuver. Off-centered half periodic cosine function transition is employed at a mutation position of angular acceleration of a BCB path; according to the angular acceleration path, predetermined maneuvering direction vector is projected to the satellite system three shafts; integration is conducted to obtain an angular velocity path; and then kinematics is utilized for calculation of attitude angle or attitude quaternion route. Compared with the prior art, the vibration of flexible attachment can be reduced greatly to improve the stability and rapidity, and realize good flexibility inhibitory effect. The design method is simple, easy for engineering implement, and suitable for satellite attitude maneuver with large flexible coupling.
Description
Technical field
The present invention relates to attitude of satellite control techniques, especially the attitude of satellite is followed the tracks of motor-driven trajectory planning, a kind of method that uses cosine transition acceleration path.
Background technology
In order to increase satellite imagery fabric width, instant observation is realized in accident area, or realize three-dimensional imaging along forward sight and the backsight of trajectory direction via satellite, the fast reserve ability of satellite has been proposed to strict requirement.It is motor-driven motor-driven with closed loop that the motor-driven method of the attitude of satellite is divided into open loop.Open loop is motor-driven, and to model, qualitative requirement is higher really, general only motor-driven for single shaft; Motor-driven step instruction mode and the path planning pattern of being divided into of closed loop, step instruction mode often has overshoot, and path planning pattern can be avoided overshoot, and mobile process is reliable gently.So, use path planning pattern to realize the attitude of satellite motor-driven, especially three-axis attitude is motor-driven is a kind of trend.
A new generation's satellite structure often needs to carry large-scale solar array or large-scale deployable antenna, the mutation process of angular acceleration is equivalent to the step response process of second-order system, easily evoke the vibration of flexible appendage, especially windsurfing, its damping ratios is very little, vibration modal frequency is low, the effect of step exciting force will cause windsurfing violent oscillatory motion last very long, have a strong impact on the attitude stability of satellite body, increase the difficulty of attitude control.For the maneuver model in trace command path, the planning quality in path has a strong impact on the vibration suppressioning effect of flexible appendage.
Prior art is checked in to published paths planning method to be had:
1) BCB path: be the more motor-driven path of current practical application, according to bang-bang control principle, according to the path of the mode planning of " accelerate-at the uniform velocity-Chang value of normal value is slowed down ", this path can make this system respond fast, but the sudden change of angular acceleration (being uncontinuity) often causes the judder of flexible windsurfing, and control accuracy is not high.As shown in Figure 1, for applying maximum BCB path schematic diagram, this path planning is simple.Nearly all open loop is motor-driven is all to adopt BCB mode, but, stability requirement more and more higher satellite more and more stronger for flexible coupling must be made improvement to instruction mode.
2) S type path: permed for 2010 the 25th in " the motor-driven path planning of satellite based on multi-objective Evolutionary Algorithm " shown at aviation power journal by Shen Xiaoning etc., proposed the motor-driven path based on S type rate curve for sudden change of acceleration problem, the method application multi-objective optimization algorithm is found the parameter of S function.But its path complex forms, and speed is larger, and topworks is easily saturated, is not suitable for Project Realization.As shown in Figure 2, for parabolic type acceleration path schematic diagram, can find out,
,
,
,
moment angular acceleration differential is larger, and compares with BCB path, does not make full use of accessible maximum angular acceleration.
3) parabolic type acceleration path: permed for 2011 the 39th in " the Spacecraft Large Angle Attitude Maneuver path planning " shown at Central China University of Science and Technology's journal by Zheng Lijun etc., propose the angular acceleration based on the parabolic function design motor-driven path of syllogic, the constant in BCB path is replaced by parabolic function.There is certain flexible inhibition, but owing to rationally not utilizing the maximum output torque of topworks, make the path time of planning longer, and the differential of segment section angular acceleration is larger, is further improved.
The problem that above-mentioned each path exists is: BCB path planning is simple, but flexible inhibition is not made to consideration; Flexible inhibition is considered in S type path, but complex forms is not easy to Project Realization; Consideration has been made to flexible inhibition in parabolic type acceleration path, but to sacrifice planning time as cost.Therefore need the path planning that can consider flexible inhibition and energy Practical.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind ofly follows the tracks of motor-driven cosine transition acceleration Path Method for the attitude of satellite, can suppress flexible appendage vibration, and form is simple, and engineering is easy to the motor-driven path of realizing.
For solving the problems of the technologies described above, of the present inventionly a kind ofly follow the tracks of motor-driven cosine transition acceleration Path Method for the attitude of satellite, it is the improvement to BCB path, adopt the half period cosine function transition of biasing in angular acceleration sudden change place in BCB path, projecting to satellite body according to angular acceleration path according to predetermined motor-driven direction vector is three axles, after integration, obtain angular velocity path, then use kinematics to resolve attitude angle or attitude quaternion path.
Described cosine transition acceleration Path Method, specifically comprises the steps:
If when the flexible coupling coefficient of satellite flexible appendage is large, choose transit time
value is large, otherwise chooses transit time
be worth little;
Step 3, choose maximum angular acceleration amplitude according to actuator stem force square restriction, satellite inertia size and motor-driven direction
Consider that topworks's moment is saturated, choose
, in formula
for topworks's maximum output torque on an axle,
for the maximum principal moments of satellite body three axles;
Step 4, determine at the uniform velocity section angular velocity size according to the limit of range of speed measuring device
Desirable in order to reserve surplus
, in formula
for the measurable maximum angular rate of speed measuring device;
Step 7, to planning angular acceleration path be assigned on body series three axles and obtain by direction vector definite in step 1
;
Cosine transition acceleration Path Method of the present invention, compared with prior art, its advantage and beneficial effect are:
1, make full use of the maximum output of topworks, on the basis of as far as possible reducing planning time, utilize cosine function transition to reduce the flexible vibration bringing in mobile process;
2, the cosine transition function adopting, owing to being zero at initial and end differential, does not have angular acceleration saltus step and the saltus step of angular acceleration differential, and flexible inhibition is good;
3, cosine transition section only need to be selected transit time and just can plan migration path, and Project Realization is simple;
4, acceleration path and known inertia battle array multiply each other, and obtain the torque command of softening, for open loop maneuver model, use such torque command can suppress the vibration of flexible appendage, and then improve degree of stability and rapidity.
Accompanying drawing explanation
Fig. 1 is existing BCB path schematic diagram;
Fig. 2 is existing parabolic type acceleration path schematic diagram;
Fig. 3 is cosine transition acceleration of the present invention path schematic diagram;
Fig. 4 is cosine transition acceleration path design process flow diagram of the present invention;
Fig. 5 is that cosine transition acceleration path contrasts figure with BCB path to flexible inhibition;
Fig. 6 is that cosine transition acceleration path and parabolic type acceleration path are to flexible inhibition
Effect comparison figure.
Embodiment
Below with reference to drawings and Examples, the present invention is described in further detail:
As shown in Figure 3, be cosine transition acceleration path schematic diagram.Adopt the half period cosine function transition of biasing in angular acceleration sudden change place in existing BCB path, projecting to satellite body according to angular acceleration path according to predetermined motor-driven direction vector is three axles, after integration, obtain angular velocity path, then use kinematics to resolve attitude angle or attitude quaternion path.
In figure
for maximum angular acceleration amplitude;
for transit time;
the constant accelerating sections time.Four transition sections adopt cosine function to describe, as first paragraph transition cosine curve is
.Needing explanation, do not provide S type path schematic diagram, be because S type path is take angular velocity as plan objects, and this method comparability is not strong, and S type path is Project Realization difficulty.
As shown in Figure 4, be the design flow diagram in cosine transition acceleration path, to adopt three axles of quaternion feedback motor-driven as example, specifically comprise the steps:
The hypercomplex number initial value of supposing the relative targeted attitude of satellite initial attitude is
, satellite need to complete the axle around Euler
rotation
;
Choose transit time by engineering is actual
if when the flexible coupling coefficient of satellite flexible appendage is large, choose transit time
value is large, otherwise chooses transit time
be worth little;
Choose maximum angular acceleration amplitude according to the restriction of actuator stem force square, satellite inertia size and motor-driven direction
, consider saturated can the choosing of topworks's moment
, in formula
for topworks's maximum output torque on an axle,
for the maximum principal moments of satellite body three axles;
Determine at the uniform velocity section angular velocity size according to the limit of range of speed measuring device
, desirable in order to reserve surplus
, in formula
for the measurable maximum angular rate of speed measuring device;
If determined constant angular acceleration
, the constant acceleration time
, transit time
afterwards, time at the uniform velocity
size need to determine according to needed motor-driven angle, at the uniform velocity motor-driven angle is
If
, illustrate and do not need at the uniform velocity section, at this moment can determine to revise according to engineering problem
,
or
, or three parameters revise simultaneously, only need to make
just can, in step 6, provide be revise
.
Step 7, to the angular acceleration path of planning by the Central European pulling shaft of step 1
direction vector be assigned on body series three axles and obtain angular acceleration vector
;
As shown in Figure 5, for using BCB path and the cosine transition acceleration path result contrast figure with identical control parameters simulation test, planning was realized in BCB path at 42 seconds, planning was realized in cosine transition acceleration path at 46 seconds, in figure, provide the 45-80 front quadravalence flexible mode figure of second, can find out residual oscillation situation, although cosine transition acceleration path has extended planning time a little, to flexible inhibiting effect successful.
As shown in Figure 6, for being used parabolic type acceleration path, same maneuver simulation test contrasts figure with the emulation that uses cosine transition acceleration path, planning was realized in parabolic type acceleration path at 55 seconds, than BCB path length 13 seconds, but but do not use cosine transition acceleration path good to the flexible vibration suppressioning effect of remnants.In figure, provide the front quadravalence flexible mode figure of 55-80 second.
Claims (2)
1. follow the tracks of motor-driven cosine transition acceleration Path Method for the attitude of satellite for one kind, it is characterized in that: to the improvement in BCB path, adopt the half period cosine function transition of biasing in angular acceleration sudden change place in BCB path, projecting to satellite body according to angular acceleration path according to predetermined motor-driven direction vector is three axles, after integration, obtain angular velocity path, then use kinematics to resolve attitude angle or attitude quaternion path.
2. according to claim 1ly follow the tracks of motor-driven cosine transition acceleration Path Method for the attitude of satellite, it is characterized in that: it specifically comprises the steps:
Step 1, calculating attitude of satellite anglec of rotation size, determine motor-driven direction;
Step 3, choose maximum angular acceleration amplitude according to actuator stem force square restriction, satellite inertia size and motor-driven direction
;
Step 4, determine at the uniform velocity section angular velocity size according to the limit of range of speed measuring device
;
Step 7, the angular acceleration path of target is assigned on satellite system three axles and is obtained by direction vector definite in step 1
;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210439266.7A CN103808323B (en) | 2012-11-07 | 2012-11-07 | A kind of attitude of satellite follows the tracks of motor-driven cosine transition angular acceleration Path Method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210439266.7A CN103808323B (en) | 2012-11-07 | 2012-11-07 | A kind of attitude of satellite follows the tracks of motor-driven cosine transition angular acceleration Path Method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103808323A true CN103808323A (en) | 2014-05-21 |
CN103808323B CN103808323B (en) | 2016-09-07 |
Family
ID=50705434
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210439266.7A Active CN103808323B (en) | 2012-11-07 | 2012-11-07 | A kind of attitude of satellite follows the tracks of motor-driven cosine transition angular acceleration Path Method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103808323B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005312A (en) * | 2015-06-29 | 2015-10-28 | 哈尔滨工业大学 | Satellite planning trajectory method based on maximum angular acceleration and maximum angular velocity |
CN106527471A (en) * | 2017-01-25 | 2017-03-22 | 上海航天控制技术研究所 | Trigonometric function track planning method used for restraining flexible vibration in attitude maneuver process and system thereof |
CN107515611A (en) * | 2017-07-28 | 2017-12-26 | 北京控制工程研究所 | A kind of sinusoidal motor-driven paths planning method of superimposed type mixing |
CN107608213A (en) * | 2017-10-12 | 2018-01-19 | 上海航天控制技术研究所 | A kind of Parameters design of the motor-driven path planning of the attitude of satellite |
CN107807657A (en) * | 2017-11-29 | 2018-03-16 | 南京理工大学 | A kind of Flexible Spacecraft self-adaptation control method based on path planning |
CN107831521A (en) * | 2017-10-16 | 2018-03-23 | 中国西安卫星测控中心 | Low orbit satellite tracks the window calculation method of non-orbital flight high dynamic target |
CN111798701A (en) * | 2020-07-07 | 2020-10-20 | 中国船舶工业系统工程研究院 | Unmanned ship path tracking control method, system, storage medium and terminal |
CN112422184A (en) * | 2020-09-28 | 2021-02-26 | 东方红卫星移动通信有限公司 | Rotation control method and device of coarse pointing device for space optical communication |
CN113830330A (en) * | 2021-09-30 | 2021-12-24 | 北京控制工程研究所 | Satellite attitude pointing method and system based on relay satellite measurement and control |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2022718A1 (en) * | 2006-05-12 | 2009-02-11 | NEC TOSHIBA Space Systems, Ltd. | Attitude control data creating method, and attitude control device applying the method |
EP2224307A2 (en) * | 2009-02-03 | 2010-09-01 | The Boeing Company | Spacecraft Acquisiton Maneuvers Using Position-Based Gyroless Control |
CN102298390A (en) * | 2011-06-24 | 2011-12-28 | 北京航空航天大学 | Anti-disturbance flexible spacecraft attitude and vibration composite control method |
CN102645773A (en) * | 2007-04-27 | 2012-08-22 | 三星电子株式会社 | Gate driving circuit and liquid crystal display having same |
-
2012
- 2012-11-07 CN CN201210439266.7A patent/CN103808323B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2022718A1 (en) * | 2006-05-12 | 2009-02-11 | NEC TOSHIBA Space Systems, Ltd. | Attitude control data creating method, and attitude control device applying the method |
CN102645773A (en) * | 2007-04-27 | 2012-08-22 | 三星电子株式会社 | Gate driving circuit and liquid crystal display having same |
EP2224307A2 (en) * | 2009-02-03 | 2010-09-01 | The Boeing Company | Spacecraft Acquisiton Maneuvers Using Position-Based Gyroless Control |
CN102298390A (en) * | 2011-06-24 | 2011-12-28 | 北京航空航天大学 | Anti-disturbance flexible spacecraft attitude and vibration composite control method |
Non-Patent Citations (3)
Title |
---|
S.PARMAN ET.AL: "Controlling the attitude maneuvers of flexible spacecraft by using time-optimal/fuel-efficient shaped inputs", 《JOURNAL OR SOUND AND VIBRATION》 * |
郑立君: "挠性卫星姿态大角度机动路径的设计与优化", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
雷拥军等: "一种航天器姿态快速机动及稳定控制方法", 《中国空间科学技术》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005312B (en) * | 2015-06-29 | 2017-11-03 | 哈尔滨工业大学 | One kind is based on maximum angular acceleration and maximum angular rate satellite planned trajectory method |
CN105005312A (en) * | 2015-06-29 | 2015-10-28 | 哈尔滨工业大学 | Satellite planning trajectory method based on maximum angular acceleration and maximum angular velocity |
CN106527471A (en) * | 2017-01-25 | 2017-03-22 | 上海航天控制技术研究所 | Trigonometric function track planning method used for restraining flexible vibration in attitude maneuver process and system thereof |
CN106527471B (en) * | 2017-01-25 | 2019-10-01 | 上海航天控制技术研究所 | Inhibit the method for planning track and system of flexible vibration during attitude maneuver |
CN107515611A (en) * | 2017-07-28 | 2017-12-26 | 北京控制工程研究所 | A kind of sinusoidal motor-driven paths planning method of superimposed type mixing |
CN107515611B (en) * | 2017-07-28 | 2020-11-10 | 北京控制工程研究所 | Superposition type hybrid sine maneuvering path planning method |
CN107608213B (en) * | 2017-10-12 | 2020-11-03 | 上海航天控制技术研究所 | Parameter design method for satellite attitude maneuver path planning |
CN107608213A (en) * | 2017-10-12 | 2018-01-19 | 上海航天控制技术研究所 | A kind of Parameters design of the motor-driven path planning of the attitude of satellite |
CN107831521A (en) * | 2017-10-16 | 2018-03-23 | 中国西安卫星测控中心 | Low orbit satellite tracks the window calculation method of non-orbital flight high dynamic target |
CN107831521B (en) * | 2017-10-16 | 2020-10-23 | 中国西安卫星测控中心 | Window calculation method for low-orbit satellite to track non-orbit flying high-dynamic target |
CN107807657B (en) * | 2017-11-29 | 2021-01-26 | 南京理工大学 | Flexible spacecraft attitude self-adaptive control method based on path planning |
CN107807657A (en) * | 2017-11-29 | 2018-03-16 | 南京理工大学 | A kind of Flexible Spacecraft self-adaptation control method based on path planning |
CN111798701A (en) * | 2020-07-07 | 2020-10-20 | 中国船舶工业系统工程研究院 | Unmanned ship path tracking control method, system, storage medium and terminal |
CN111798701B (en) * | 2020-07-07 | 2022-07-26 | 中国船舶工业系统工程研究院 | Unmanned ship path tracking control method, system, storage medium and terminal |
CN112422184A (en) * | 2020-09-28 | 2021-02-26 | 东方红卫星移动通信有限公司 | Rotation control method and device of coarse pointing device for space optical communication |
CN113830330A (en) * | 2021-09-30 | 2021-12-24 | 北京控制工程研究所 | Satellite attitude pointing method and system based on relay satellite measurement and control |
CN113830330B (en) * | 2021-09-30 | 2023-08-29 | 北京控制工程研究所 | Satellite attitude pointing method and system based on relay satellite measurement and control |
Also Published As
Publication number | Publication date |
---|---|
CN103808323B (en) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103808323A (en) | Cosine transition acceleration path method for satellite attitude tracking maneuver | |
Gao et al. | Multi-objective robust $ H_ {\infty} $ control of spacecraft rendezvous | |
CN103412491B (en) | A kind of Spacecraft feature axis attitude maneuver index time-varying sliding-mode control | |
CN106697333B (en) | A kind of robust analysis method of spacecraft orbit control strategy | |
CN105911866B (en) | The finite time full-order sliding mode control method of quadrotor unmanned vehicle | |
CN110347170B (en) | Reusable carrier reentry segment robust fault-tolerant guidance control system and working method | |
CN104570742B (en) | Feedforward PID (proportion, integration and differentiation) control based rapid high-precision relative pointing control method of noncoplanar rendezvous orbit | |
Zhao et al. | Output feedback continuous terminal sliding mode guidance law for missile-target interception with autopilot dynamics | |
CN105786024A (en) | Airborne photoelectric platform high precision tracking controller based on model error compensation and tracking control method thereof | |
CN103955223A (en) | Posture and path coupling control method for deep space exploration soft landing process | |
CN107479566A (en) | Flexible satellite attitude control method based on three-stage path planning | |
CN104309822B (en) | A kind of spacecraft single impulse water-drop-shaped based on parameter optimization is diversion track Hovering control method | |
Gui et al. | Attitude tracking of a rigid spacecraft using two internal torques | |
CN107807657A (en) | A kind of Flexible Spacecraft self-adaptation control method based on path planning | |
CN109283841A (en) | Flexibility of joint and armed lever space manipulator control method flexible | |
CN105607485A (en) | Adaptive fault tolerance control method for flexible liquid-filled satellite attitude based on fault characteristic model | |
CN103235509A (en) | Rotating member disturbance compensation method based on momentum wheel | |
CN109426147A (en) | The adaptive gain that spacecraft is combined after capture satellite adjusts control method | |
CN107861386A (en) | A kind of anti-interference attitude control ground checking system and its control method based on angular speed observer | |
Zhang et al. | Combined control of fast attitude maneuver and stabilization for large complex spacecraft | |
CN102508492B (en) | Method for realizing great circle flight of aircraft in constant height between isometric waypoints | |
Xiaojian et al. | Impact angle control over composite guidance law based on feedback linearization and finite time control | |
Chen et al. | Decoupling attitude control of a hypersonic glide vehicle based on a nonlinear extended state observer | |
CN106005483A (en) | Active attitude control method of modular mobile phone satellite | |
Gao et al. | Attitude tracking control of a quadrotor based on linear active disturbance rejective control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |