CN101762284B - Determining method of single-frame control moment gyro dynamic unbalance disturbance moment - Google Patents
Determining method of single-frame control moment gyro dynamic unbalance disturbance moment Download PDFInfo
- Publication number
- CN101762284B CN101762284B CN2009102432691A CN200910243269A CN101762284B CN 101762284 B CN101762284 B CN 101762284B CN 2009102432691 A CN2009102432691 A CN 2009102432691A CN 200910243269 A CN200910243269 A CN 200910243269A CN 101762284 B CN101762284 B CN 101762284B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- gyro
- rotor
- sin
- cos
- 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
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
The present invention relates to a determining method of positive cone-shaped single-frame control moment gyro group motion unbalance disturbance moment; a rotor coordinate system, a frame coordinate system and a gyro base coordinate system are built, the gyro rotor disturbance moment in the rotor coordinate system determines the disturbance moment of the gyro rotor in the star coordinate system by continuous coordinate conversion; the determining method fills up the technology blank, attitude disturbance analysis and flutter analysis to quick moonlet are carried out possibly by utilizing disturbance moment continuously, so as to provide basis for judging satellite imagery quality and improving image quality in earlier stage; the determining method can adjust the disturbance moment of the moment gyro by adjusting the installation mode of the gyro and the processing precision of the gyro, so as to achieve the purpose of improving the image quality finally; in the method, the disturbance moment generated by the rotor dynamic unbalance is determined from the point of the disturbance moment generated by the dynamic unbalance of a fixed shaft rotating part by combining the two-freedom-degree rotation of the single-frame control moment gyro rotor through continuous coordinate conversion, the method is simple and easy.
Description
Technical field
The present invention relates to a kind of definite method of control-moment gyro unbalance dynamic disturbance torque, particularly the disturbance torque of single frame frame control moment gyro dynamic unbalance generation is determined method, belong to the satellite dynamics technical field, can be used for attitude of satellite interference, flutter analysis.
Background technology
People such as Yin Qiuyan 2007 in " disturbance restraining method research in the counteraction flyback " literary composition from the counteraction flyback movement characteristic, labor the wheel sound is uneven disturbs, set up corresponding mathematics model, adopted corresponding inhibition method and carried out simulation calculation according to the effect characteristics of disturbing.The disturbance torque that produces for the reaction wheel unbalance dynamic determines that research is more at present, and is also comparatively ripe, but less for the research of control-moment gyro disturbance torque.
Development along with quick moonlet, control-moment gyro becomes the preferred configuration of attitude control actuator so that bigger angular momentum ability can be provided, yet the small imbalance that the control-moment gyro rotor quality distributes equally can be in rotary course produces direction and size all in time and the disturbance that changes causes vibratory response, attitude disturbance response and the flutter response of satellite to satellite.At present, also do not have a kind of special method of determining the unbalance dynamic disturbance torque of control-moment gyro, cannot carry out attitude disturbance analysis, flutter analysis etc. to the quick moonlet that adopts control-moment gyro, test the resolution of satellite and whether adhere to specification thereby also just can't pass judgment on the moonlet image quality according to the flutter analysis result.
The present invention is directed to present technological gap, a kind of definite method of disturbance torque of single frame control-moment gyro unbalance dynamic generation has been proposed first, the disturbance torque that utilizes the present invention to determine, can carry out attitude disturbance analysis, flutter analysis to quick moonlet, thereby pass judgment on the satellite imagery quality, the resolution of test satellite can also finally make the resolution of satellite reach designing requirement by the disturbance torque of adjusting control-moment gyro.
The disturbance torque that relevant controlling moment gyro unbalance dynamic at present causes determines that method do not see disclosed pertinent literature report as yet.
Summary of the invention
Technology of the present invention is dealt with problems and is: overcome the deficiencies in the prior art, a kind of definite method of disturbance torque of single frame control-moment gyro unbalance dynamic generation is provided, fill up the prior art blank, the inventive method provides the input foundation for satellite flutter analysis and flutter to the impact analysis and the final satellite imagery quality of improving of satellite imagery quality.
Technical solution of the present invention is: a kind of definite method of single-frame control moment gyro dynamic unbalance disturbance moment, realize by following steps:
The first step is utilized the rotor quality characteristic test method, obtains the amount J of the principal axis of inertia relative rotation axi line inclined degree of each gyrorotor of tolerance control-moment gyro group
Xz, J
Yz
Second step obtained the dynamic unbalance value and the starting phase angle of gyrorotor according to formula group (1),
Wherein, I
0Be the dynamic unbalance value of gyrorotor, φ is the starting phase angle of gyrorotor dynamic unbalance;
The 3rd step, satellite control system topworks measure portion record each gyrorotor of control-moment gyro group relatively the corner of self framework be angle of rotor γ
iThe corner of relative gyro pedestal with framework is framework corner ζ
i, and the angular velocity that rotates around turning axle of rotor
I=1,2...N, N are total numbers of control-moment gyro;
In the 4th step, utilized for second step obtained I
0Go on foot the angular velocity that obtains with φ and the 3rd
Obtain gyrorotor at rotor coordinate system ox according to formula (2)
my
mz
mMiddle disturbance torque,
T wherein
DmFor gyrorotor at rotor coordinate system ox
my
mz
mThe disturbance torque of following generation,
Be rotor coordinate system ox
my
mz
mX
mTo, y
mTo unit vector;
In the 5th step, utilize the 3rd to go on foot the angle of rotor γ that obtains
iWith framework corner ζ
i, obtaining frame coordinates according to formula (3), (4) is ox
ry
rz
rRelative gyro base coordinate system ox
sy
sz
sTransition matrix A
RsiWith rotor coordinate system ox
my
mz
mFrame coordinates is ox relatively
ry
rz
rTransition matrix A
Mri,
A
RsiThe frame coordinates that is i gyro is the transition matrix of relative gyro base coordinate system, A
MriIt is the transition matrix of the relative frame coordinates of the rotor coordinate system system of i gyro;
In the 6th step, obtain each gyro base coordinate system ox of positive taper control-moment gyro group by formula (5)
sy
sz
sTo celestial body coordinate system o
bx
by
bz
bTransition matrix, each gyro gimbal axle and satellite celestial body coordinate system o
bx
by
bz
bZ
bAxle coplane, and and z
bThe angle of axle is β, and gimbal axis is at satellite celestial body coordinate plane x
by
bProjection and x
bThe axle clamp angle is α
i,
Wherein, M
SbiBe the transition matrix of i gyro base coordinate system to the celestial body coordinate system;
In the 7th step, utilize the transition matrix A that obtains in the 5th step
Rsi, A
MriWith the 6th transition matrix M that obtains of step
Sbi, formula (2) through three coordinate conversion, is transformed into satellite celestial body coordinate system from the rotor coordinate system, obtain the disturbance torque T that every i gyrorotor unbalance dynamic produces in satellite celestial body coordinate system
Di,
In the 8th step, obtain the disturbance torque T that single frame control-moment gyro group produces in satellite celestial body coordinate system according to formula (6)
d,
Design concept of the present invention:
1, sets up the coordinate system of using in the analysis
Definition rotor coordinate system ox
my
mz
m, be the rotor body coordinate system, be connected with rotor;
Definition frame coordinate system ox
ry
rz
r, coordinate axis along gimbal axis, axis of torque and rotor angular momentum axle, is connected with framework respectively, initial time ox
ry
rz
rAnd ox
my
mz
mOverlap;
Definition gyro base coordinate system ox
sy
sz
s, when frame corners was zero, it and frame coordinates were ox
ry
rz
rOverlap.Arrange according to the configuration of gyro group and the direction of each gyro gimbal axle, can determine the installation matrix M of each gyro pedestal.Each coordinate system as shown in Figure 1.
2, the conversion between each coordinate system
The corner of the relative framework of control-moment gyro rotor is designated as γ, and the corner of the relative gyro pedestal of framework is designated as ζ.The gyro base coordinate system can be obtained by twice rotation to the conversion of rotor coordinate system, promptly earlier around x
rAxle turns over ζ, again around z
mAxle turns over γ.
Frame coordinates is ox
ry
rz
rRelative gyro base coordinate system ox
sy
sz
sTransition matrix be
Rotor coordinate system ox
my
mz
mFrame coordinates is ox relatively
ry
rz
rTransition matrix be
For the single frame control-moment gyro group of taper configuration, a certain axon of the gimbal axis of every gyro and celestial body (for example z axle) coplane, with this axle clamp angle β, gimbal axis is at celestial body coordinate plane x
by
bProjection and x
bAxle clamp angle α, the angular momentum of frame corners zero-bit is parallel to x
by
bThe plane can determine that each gyro base coordinate system gets transition matrix to the celestial body coordinate system, and promptly the installation matrix of gyro pedestal is
3, determine disturbance torque that the rotor unbalance dynamic produces according to d'Alembert principle under the rotor coordinate system, the moment of inertia that during rotation of dead axle rotatable parts rotation axis is produced is:
Particularly when constant speed is rotated
Wherein
Be the angular velocity that rotatable parts rotate around dead axle, U
dBe unbalance dynamic, do not overlap by rotor principal axis of inertia and rotation and cause J
Xz, J
YzBe the amount of the principal axis of inertia relative rotation axi line inclined degree of tolerance rotor, obtain by the rotor quality characteristic test.
For the single frame control-moment gyro, its rotor is done uniform rotation around the axis of angular momentum, so the disturbance torque T of rotor unbalance dynamic generation
DmAt rotor coordinate system o
mx
my
mz
mIn be
4, determine the disturbance torque that the rotor unbalance dynamic produces under the celestial body coordinate system
Disturbance torque expression formula according to the rotor unbalance dynamic produces under the rotor coordinate system is transformed into it under celestial body coordinate system by continuous coordinate transform, finally determines the disturbance torque that the rotor unbalance dynamic produces under the celestial body coordinate system.
At celestial body o
bx
by
bz
bIn the coordinate system, T
dFor
With each transition matrix substitution above, must the rotor unbalance dynamic be to the disturbance torque that satellite body produces
The present invention compared with prior art beneficial effect is:
(1) the present invention proposes definite method of the disturbance torque of single frame control-moment gyro rotor unbalance dynamic generation first, filled up technological gap, make the later use disturbance torque to quick moonlet carry out the attitude disturbance analysis, flutter analysis becomes possibility, thereby for passing judgment on the satellite imagery quality, improve picture quality, foundation in earlier stage is provided;
(2) machining precision of mounting means that definite method of the present invention can also be by regulating gyro and gyro is adjusted the disturbance torque of control-moment gyro, finally reaches the purpose that improves picture quality;
(3) the present invention is set out by the disturbance torque that the unbalance dynamic of dead axle rotatable parts produces, and in conjunction with the two-freedom rotation of single frame control-moment gyro rotor, by continuous coordinate conversion, determines the disturbance torque that its rotor unbalance dynamic produces, and method is simple.
Description of drawings
Fig. 1 is a single frame control-moment gyro coordinate system synoptic diagram of the present invention;
Fig. 2 is the positive five face cone shape single frame control-moment gyro group configuration synoptic diagram of the present invention.
Embodiment
Be configured as example with five face cone shape single frame gyro groups and specify performing step of the present invention.
One, the gimbal axis of each gyro in the gyro group and satellite celestial body coordinate system o
bx
by
bz
bZ
bAxle coplane, and and z
bThe angle of axle is β, and gimbal axis is at satellite celestial body coordinate plane x
by
bProjection and x
bAxle clamp angle α
i, i=1,2 ... 5, each gyro base coordinate system ox then
sy
sz
sTo celestial body coordinate system o
bx
by
bz
bTransition matrix M
SbiFor
As shown in Figure 2, five face cone shape single frame gyro groups, the gimbal axis of each gyro in the gyro group is at satellite celestial body coordinate plane x
by
bProjection and x
bThe axle clamp angle is respectively α
1=90 °, α
2=π-18 °, α
3=-(π-54 °), α
4=-54 °, α
5=18 °, can get the installation matrix of five gyro pedestals in the substitution formula (5).
Two, utilize the rotor quality characteristic test method, obtain the amount J of the principal axis of inertia relative rotation axi line inclined degree of each gyrorotor of tolerance
Xz, J
Yz
Three, obtain the dynamic unbalance value and the starting phase angle of gyrorotor according to formula group (1),
Wherein, I
0Be the dynamic unbalance value of gyrorotor, φ is the starting phase angle of gyrorotor dynamic unbalance.
Four, satellite control system topworks measure portion record each gyrorotor relatively the corner of self framework be angle of rotor γ
iThe corner of relative gyro pedestal with framework is framework corner ζ
i, and the angular velocity that rotates around turning axle of rotor
Five, obtain each gyrorotor at rotor coordinate system ox according to formula (2)
my
mz
mMiddle disturbance torque,
T wherein
DmFor gyrorotor at rotor coordinate system ox
my
mz
mThe disturbance torque of following generation,
Be rotor coordinate system ox
my
mz
mx
mTo, y
mTo unit vector.
Six, obtaining frame coordinates according to formula (3), (4) is ox
ry
rz
rRelative gyro base coordinate system ox
sy
sz
sTransition matrix A
RsiWith rotor coordinate system ox
my
mz
mFrame coordinates is ox relatively
ry
rz
rTransition matrix A
Mri
A
RsiThe frame coordinates that is i gyro is the transition matrix of relative gyro base coordinate system, A
MriIt is the transition matrix of the relative frame coordinates of the rotor coordinate system system of i gyro.
Seven, utilize the transition matrix M of formula (5)
SbiAnd the transition matrix A of formula (3), (4)
Rsi, A
Mri, formula (2) through three coordinate conversion, is transformed into satellite celestial body coordinate system from the rotor coordinate system, obtain the disturbance torque T of i gyrorotor in satellite celestial body coordinate system
Di,
Eight, obtain the disturbance torque T of gyro group in satellite celestial body coordinate system according to formula (6)
d,
For the single frame moment gyro group of five face cone shape configurations, ignore the starting phase angle φ of rotor dynamic unbalance
i, can get its rotor unbalance dynamic and be total disturbance torque that celestial body produces
In the formula, ζ
i(i=1,2 ..., 5) and be respectively the corner of five frameworks; γ is the corner of rotor, and five control-moment gyro angle of rotor are identical, is provided by satellite control system topworks measure portion.
So far, determined the uneven total disturbance torque that produces of positive five face cone shape single frame control-moment gyro group motions, for the satellite flutter is determined also and then carry out evaluation of imaging quality the anticipated import foundation is provided.
Above-mentioned definite total disturbance torque substitution based on the flexible spacecraft dynamic modeling method, in conjunction with the stiffness characteristics of moonlet parts, is mainly considered the flexibility of the sun wing, the moonlet flexible spacecraft trembling vibration mechanical equation of foundation
In, utilize the derivation algorithm of ordinary differential equation initial-value problem to resolve, obtain the attitude response of satellite under the effect of flutter disturbance torque, the attitude response of satellite under the effect of flutter disturbance torque deducted utilize the controlled attitude part of attitude of satellite control rule, obtain the flutter response that satellite causes under the effect of flutter disturbance torque, and then carry out satellite flutter response the influence of image quality is passed judgment on.
The unspecified part of the present invention belongs to general knowledge as well known to those skilled in the art.
Claims (1)
1. definite method of the uneven disturbance torque of positive taper single frame control-moment gyro group motion is characterized in that realizing by following steps:
The first step is utilized the rotor quality characteristic test method, obtains the amount J of the principal axis of inertia relative rotation axi line inclined degree of described each gyrorotor of gyro group of tolerance
Xz, J
Yz
Second step obtained the dynamic unbalance value and the starting phase angle of gyrorotor according to formula group (1),
Wherein, I
0Be the dynamic unbalance value of gyrorotor, φ is the starting phase angle of gyrorotor dynamic unbalance;
The 3rd step, satellite control system topworks measure portion record each gyrorotor of described gyro group relatively the corner of self framework be angle of rotor γ
iThe corner of relative gyro pedestal with framework is framework corner ζ
i, and the angular velocity that rotates around turning axle of rotor
I=1,2...N, N are total numbers of the gyro in the described gyro group;
In the 4th step, utilized for second step obtained I
0Go on foot the angular velocity that obtains with φ and the 3rd
Obtain gyrorotor at rotor coordinate system ox according to formula (2)
my
mz
mMiddle disturbance torque,
T wherein
DmFor gyrorotor at rotor coordinate system ox
my
mz
mThe disturbance torque of following generation,
Be rotor coordinate system ox
my
mz
mX
mTo, y
mTo unit vector;
In the 5th step, utilize the 3rd to go on foot the angle of rotor γ that obtains
iWith framework corner ζ
i, obtaining frame coordinates according to formula (3), (4) is ox
ry
rz
rRelative gyro base coordinate system ox
sy
sz
sTransition matrix A
RsiWith rotor coordinate system ox
my
mz
mFrame coordinates is ox relatively
ry
rz
rTransition matrix A
Mri,
A
RsiThe frame coordinates that is i gyro is the transition matrix of relative gyro base coordinate system, A
MriIt is the transition matrix of the relative frame coordinates of the rotor coordinate system system of i gyro;
In the 6th step, obtain each gyro base coordinate system ox of described gyro group by formula (5)
sy
sz
sTo satellite celestial body coordinate system o
bx
by
bz
bTransition matrix, each gyro gimbal axle and satellite celestial body coordinate system o
bx
by
bz
bZ
bAxle coplane, and and z
bThe angle of axle is β, and gimbal axis is at satellite celestial body coordinate plane x
by
bProjection and x
bThe axle clamp angle is α
i,
Wherein, M
SbiBe the transition matrix of i gyro base coordinate system to satellite celestial body coordinate system;
In the 7th step, utilize the transition matrix A that obtains in the 5th step
Rsi, A
MriWith the 6th transition matrix M that obtains of step
Sbi, formula (2) through three coordinate conversion, is transformed into satellite celestial body coordinate system from the rotor coordinate system, obtain the disturbance torque T that every i gyrorotor unbalance dynamic produces in satellite celestial body coordinate system
Di,
In the 8th step, obtain the disturbance torque T that described gyro group produces in satellite celestial body coordinate system according to formula (6)
d,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102432691A CN101762284B (en) | 2009-12-30 | 2009-12-30 | Determining method of single-frame control moment gyro dynamic unbalance disturbance moment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102432691A CN101762284B (en) | 2009-12-30 | 2009-12-30 | Determining method of single-frame control moment gyro dynamic unbalance disturbance moment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101762284A CN101762284A (en) | 2010-06-30 |
CN101762284B true CN101762284B (en) | 2011-11-16 |
Family
ID=42493585
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102432691A Active CN101762284B (en) | 2009-12-30 | 2009-12-30 | Determining method of single-frame control moment gyro dynamic unbalance disturbance moment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101762284B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102063521B (en) * | 2010-10-12 | 2012-07-25 | 北京理工大学 | Design method for configuration-adjustable single-framework control moment gyro system |
CN102176159B (en) * | 2011-02-28 | 2013-05-01 | 哈尔滨工业大学 | Satellite attitude control system failure diagnosis device and method based on state observer and equivalent space |
CN102778891B (en) * | 2012-08-03 | 2014-04-16 | 北京理工大学 | Parameter selection method adopting onboard control moment gyroscope group vibration-isolating platform |
CN106933241B (en) * | 2017-03-30 | 2019-11-29 | 北京航空航天大学 | Single-gimbal control momentum gyro spacecraft fault tolerant control method based on fault de couple |
CN107796546B (en) * | 2017-09-22 | 2020-07-14 | 上海卫星工程研究所 | Dynamic measurement method for on-orbit output torque of satellite torque gyro set |
-
2009
- 2009-12-30 CN CN2009102432691A patent/CN101762284B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN101762284A (en) | 2010-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101750200A (en) | Method for determining flutter response of high-resolution minisatellites | |
CN101762284B (en) | Determining method of single-frame control moment gyro dynamic unbalance disturbance moment | |
US9536446B2 (en) | Motion simulation system controller and associated methods | |
CN104527994B (en) | Multi-polar cross-over becomes the track set time soon and holds position sensing tracking and controlling method | |
CN106672266B (en) | Balance point Halo orbit phase modulation orbit transfer method considering time constraint | |
CN103019091B (en) | Flexible spacecraft fault-tolerant attitude control method based on linear extended state observer | |
JP4808354B2 (en) | Posture estimation of tiltable objects using modified quaternion data notation | |
CN107628272B (en) | Satellite rotary part Wheel static and dynamic imbalance torque self-compensating device | |
CN106896821B (en) | A kind of angular momentum management method of variable speed control moment gyro | |
CN104570742B (en) | Feedforward PID (proportion, integration and differentiation) control based rapid high-precision relative pointing control method of noncoplanar rendezvous orbit | |
CN107600464A (en) | Utilize the flywheel control capture sun and Direct to the sun method of star sensor information | |
CN106586034A (en) | Self-compensating method for dynamic and still unbalancing disturbance moment of satellite rotating part | |
CN102323825B (en) | Torque compensation control method of DGMSCMG (double-gimbal magnetically suspended control moment gyroscope) system for spacecraft maneuver | |
CN109033604B (en) | Satellite dynamics modeling method containing rotation load and determination method of stress at bearing | |
CN107402516B (en) | Rank saturation the fuzzy PD attitude control method is passed based on joint executing agency | |
CN108128484B (en) | Double-star system track keeping method based on linear quadratic regulator | |
CN104166345A (en) | Decoupling and disturbance-attenuation method for rotor system of magnetic suspension control moment gyroscope | |
CN106933241A (en) | Single-gimbal control momentum gyro spacecraft fault tolerant control method based on fault de couple | |
CN105259906A (en) | Apparatus and method of improving spacecraft attitude stability | |
CN103955138B (en) | A kind of dynamic middle imaging satellite attitude control method based on increment type drift angle | |
CN203966431U (en) | The target detection of a kind of dexterous ammunition device and control analogue system | |
CN112857348A (en) | Angular velocity measuring method using magnetic suspension bearing | |
CN109018442B (en) | Novel low-cost satellite three-axis attitude time-sharing decoupling high-multiplexing air injection control method | |
Saberi et al. | Design and analysis of gimbal thruster configurations for 3-axis satellite attitude control | |
JPWO2021024591A5 (en) |
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 |