CN109612667A - The in-orbit end displacement method and system of satellite flexible appendage is recognized using gyro data - Google Patents
The in-orbit end displacement method and system of satellite flexible appendage is recognized using gyro data Download PDFInfo
- Publication number
- CN109612667A CN109612667A CN201910020517.XA CN201910020517A CN109612667A CN 109612667 A CN109612667 A CN 109612667A CN 201910020517 A CN201910020517 A CN 201910020517A CN 109612667 A CN109612667 A CN 109612667A
- Authority
- CN
- China
- Prior art keywords
- vibration
- axis
- orbit
- satellite
- flexible appendage
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
- G01M7/025—Measuring arrangements
Abstract
The present invention provides a kind of methods using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage comprising the steps of: attitude angular velocity obtaining step: the corresponding attitude angular velocity measurement data ω of i axis in the underdamping free vibration area after obtaining the in-orbit jet closed-loop control of satellitei(t)；Filter step: to ωi(t) it is filtered and obtains filtered i axis attitude angular velocity data ω 'i(t)；Time series calculates step: according to ω 'i(t) the time series η for obtaining mode variable is calculatedi(t)；End displacement calculates step: according to ηi(t) with end node formation data diThe in-orbit vibration end displacement m of satellite flexible appendage is calculatedi(t).Correspondingly, the present invention also provides a kind of systems using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage.This method does not need other installation sensors on flexible appendage, is analyzed and processed merely with the measurement data of the existing attitude sensor of satellite platform.
Description
Technical field
The present invention relates to satellite monitoring fields, and in particular, to a kind of to be existed using gyro data identification satellite flexible appendage
The method and system of rail vibration end displacement.A kind of more particularly to in-orbit identification method for flexible satellite vibrational state information
With system, the discrimination method and system of end displacement during the specifically a kind of pair of in-orbit vibration of flexible satellite.
Background technique
With the continuous development of space technology, large space structure has been an important development side in space industry
To also will being the infrastructure [in-orbit spacecraft dynamics parameter identification technique research-Yu Dengyun] of space development indispensability.For
Flexible satellite structural dynamic characteristics are obtained, calculates and vibrates frequently with the methods of finite element modeling, ground experiment and in-orbit identification
Modal parameter.Model simplification, condition hypothesis in finite element modeling process etc. affect model accuracy, are especially difficult to hinge
Etc. contacts mechanism Accurate Model.Due to the influence of the factors such as gravity and atmospheric drag, Large Flexible Structure be difficult ground into
Row assembles and implements full-scale dynamic parameters identification test.Therefore, for the spacecraft with large-scale flexible appendage, it is difficult to
Accurate structural dynamic parameter is obtained by finite element modeling or ground experiment.Based on above-mentioned factor, to sun battle array, space exhibition
It is very urgent and necessary while also with higher that the big flexible structure such as Kaitian's line, which carries out in-orbit Modal Parameter Identification research,
Theory significance and practical application value.Flexible satellite structural dynamic parameter includes modal frequency, modal damping, formation, coupling
Coefficient etc., these parameters have important physical significance, can for the structure design of space flexible part, monitoring structural health conditions,
The application of structure failure diagnosis, structural vibration control etc. provides necessary support [CN102982196A].
The Modal Parameters Identification in Structural Dynamics field mainly includes frequency domain method, time domain approach and in recent years at present
Three kinds of the time-frequency domain method of rise.Most methods are needed using in-orbit excitation and placement sensor obtains letter on flexible appendage
Breath, such as [large-sized solar windsurfing modal parameter in-orbit identification research _ Li Xiaoran-Master's thesis] has studied senses on flexible appendage
The minimal configuration quantity and optimal location scheme of device；[CN105486474A] describes a kind of in-orbit mode of satellite flexible part
The realization system and method for identification needs to carry out pulse excitation to flexible appendage, receives and monitor the impulse response of each measuring point
Signal, and acquire acceleration responsive signal caused by the in-orbit steady state operation Satellite flexible part of satellite；
[CN106557633A] and [CN107609296A] describes the two methods of satellite sun array sensor placement.
The solar sail deployed configuration dynamic characteristic in-orbit identification method of [CN106408570A] based on Binocular vision photogrammetry, from image directly
Structural vibration displacement information is extracted, and then obtains the dynamic characteristic of structure in real time by operational modal analysis technology, realizes structure
The in-orbit identification of dynamic characteristic.Also two patents propose the attitude angular velocity data using gyro to measure to flexible satellite mode
The in-orbit identification method of parameter: [CN103970964A] describes a kind of flexible satellite modal parameter in-orbit identification method, needs
The angular velocity information of torque and flexible satellite body relative inertness coordinate system that executing agency is applied on flexible satellite body is collected,
The transmission function of modal parameter with torque to angular speed is obtained using Subspace Identification algorithm.[CN105157728A] proposes one
It kind can inhibit the flexible satellite Modal Parameters Identification of gyro noise influence, satellite body when flight in-orbit also with satellite
The measurement data of angular speed recognizes the modal frequency and damping ratios parameter of whole star, carries out difference processing inhibition to gyro data
Identification Errors caused by this two parts noise.
There are two constraints in terms of engineer application for the above in-orbit dynamic parameters identification method of flexible satellite: first is that in-orbit boat
Its device is difficult to required for applying dynamic parameters identification known excitation, can only utilize the expansion and gathering of in-orbit space structure,
The docking and separation of structure, the igniting of engine etc. generate driving source and are motivated, and these excitation source signals are all difficult to measure
[in-orbit spacecraft dynamics parameter identification technique research-Yu Dengyun].Second is that installing the number of vibrating sensor on flexible appendage
Amount receives the limitation of engineering construction.General flexible appendage is required to in-orbit expansion, all kinds of speed, accelerations, displacement sensor
Installation fixation, cable laying etc. have an adverse effect to unfolding mechanism, increase design difficulty.In addition, for sun battle array
Equal rotational flexibilities component also needs to increase driving mechanism slip ring signal number of channels.
Summary of the invention
For the defects in the prior art, it is attached using gyro data identification satellite flexibility that the object of the present invention is to provide a kind of
The method and system of the in-orbit vibration end displacement of part.
The method using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage provided according to the present invention, packet
Containing following steps:
Attitude angular velocity obtaining step: the i axis in the underdamping free vibration area after obtaining the in-orbit jet closed-loop control of satellite
Corresponding attitude angular velocity measurement data ωi(t), wherein i axis is X-axis, Y-axis, any in Z axis in rectangular coordinate system in space
A axis；
Filter step: to ωi(t) it is filtered and obtains filtered i axis attitude angular velocity data ωi′(t)；
Time series calculates step: according to ωi' (t) calculates the time series η for obtaining mode variablei(t)；
End displacement calculates step: according to ηi(t) with end node formation data diIt is in-orbit that satellite flexible appendage is calculated
Vibrate end displacement mi(t)。
Preferably, in attitude angular velocity obtaining step, ωiIt (t) is that underdamping is freely shaken after the in-orbit gas puff Z-pinch of satellite
It in dynamic period, is obtained by satellite gyroscope measurement.
Preferably, the filter step comprises the steps of:
Concern vibration frequency band setting steps: the concern vibration frequency band of setting i axis direction is [f1,f2]Hz；
Strip step: enable bandpass filter to ωi(t) ω is removed in row filtering processingiIt (t) is more than the period of setting value in
Movement and high-frequency vibration component only retain the vibration information closed in vibration injection band, obtain ωi′(t)。
It preferably, also include frequency band set-up procedure, the frequency band set-up procedure comprises the steps of:
Global displacement obtaining step: by ηi(t) with mode of oscillation formation data DiIt is multiplied, obtains each node on flexible appendage
In the global displacement situation of setting period；
Principal oscillation mode judgment step: according to the global displacement situation, judge flexible appendage on each direction of vibration
Whether a principal oscillation mode is only had；If so, continuing to execute end displacement calculates step；If it is not, then returning to execution filtering step
Suddenly and concern vibration frequency band is adjusted.
Preferably, in concern vibration frequency band setting steps, according to the whole star structural dynamical model calculated result in ground, setting
Frequency band is vibrated in the concern of i axis direction.
Preferably, time series calculates step and comprises the steps of:
First derivative step of converting: according to the following formula to ωi' (t) is handled:
In formula:For the first derivative of the corresponding mode variable of i axis；
JiFor satellite i axis direction rotary inertia；
BiFor the rotation coefficient of coup of the corresponding flexible appendage mode of oscillation of i axis；
Integration step: rightIntegral calculation is carried out, η is obtainedi(t)。
Preferably, mi(t) it is obtained by following calculation formula:
mi(t)=diηi(t)。
The in-orbit system for vibrating end displacement of satellite flexible appendage is recognized using gyro data the present invention also provides a kind of,
Comprising with lower module:
Attitude angular velocity obtains module: the i axis in the underdamping free vibration area after obtaining the in-orbit jet closed-loop control of satellite
Corresponding attitude angular velocity measurement data ωi(t), wherein i axis is X-axis, Y-axis, any in Z axis in rectangular coordinate system in space
A axis；
Filter module: to ωi(t) it is filtered and obtains filtered i axis attitude angular velocity data ωi′(t)；
Time series computing module: according to ωi' (t) calculates the time series η for obtaining mode variablei(t)；
End displacement computing module: according to ηi(t) with end node formation data diIt is in-orbit that satellite flexible appendage is calculated
Vibrate end displacement mi(t)。
Preferably, the filter module includes with lower module:
Concern vibration frequency band setup module: the concern vibration frequency band of setting i axis direction is [f1,f2]Hz；
Strip module: enable bandpass filter to ωi(t) ω is removed in row filtering processingiIt (t) is more than the period of setting value in
Movement and high-frequency vibration component only retain the vibration information closed in vibration injection band, obtain ωi′(t)；
System using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage also includes frequency band adjustment module,
The frequency band adjustment module includes with lower module:
Global displacement obtains module: by ηi(t) with mode of oscillation formation data DiIt is multiplied, obtains each node on flexible appendage
In the global displacement situation of setting period；
Principal oscillation mode judgment module: according to the global displacement situation, judge flexible appendage on each direction of vibration
Whether a principal oscillation mode is only had；If so, continuing to execute end displacement computing module；If it is not, then returning to execution filtering mould
Block simultaneously adjusts concern vibration frequency band.
Preferably, mi(t) it is obtained by following calculation formula:
mi(t)=diηi(t)。
Compared with prior art, the present invention have it is following the utility model has the advantages that
1, the present invention meets the needs of in-orbit structural dynamic characteristics of flexible satellite recognize, and extracts the in-orbit vibration of flexible satellite
End displacement information.
2, compared with other calculate the method for the in-orbit structural dynamic parameter of satellite, this method is not needed on flexible appendage
Other installation sensors, are analyzed and processed merely with the measurement data of the existing attitude sensor of satellite platform.
3, increased design of satellites, manufacture difficulty and wind in orbit are avoided because installing other vibrating sensors
Danger.
4, the flexible appendage vibrational state information picked out can be used for the in-orbit structural dynamical model of satellite, fault diagnosis,
Performance Evaluation etc..
Detailed description of the invention
Upon reading the detailed description of non-limiting embodiments with reference to the following drawings, other feature of the invention,
Objects and advantages will become more apparent upon:
Fig. 1 is certain single-blade sun battle array satellite configuration schematic diagram；
Fig. 2 is the attitude angular velocity change curve in satellite in-orbit jet closed-loop control period；
Fig. 3 is the time domain comparison diagram of measuring satellite angular velocities filtering front and back；
Fig. 4 is the frequency domain comparison diagram of measuring satellite angular velocities filtering front and back；
Fig. 5 is the change curve of flexible mode variable；
Fig. 6 is the formation schematic diagram of flexible appendage mode of oscillation；
Fig. 7 is the end displacement change curve of flexible appendage vibration；
Fig. 8 is the method flow diagram using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage.
Specific embodiment
The present invention is described in detail combined with specific embodiments below.Following embodiment will be helpful to the technology of this field
Personnel further understand the present invention, but the invention is not limited in any way.It should be pointed out that the ordinary skill of this field
For personnel, without departing from the inventive concept of the premise, various modifications and improvements can be made.These belong to the present invention
Protection scope.
In the description of the present invention, it is to be understood that, term " on ", "lower", "front", "rear", "left", "right", " perpendicular
Directly ", the orientation or positional relationship of the instructions such as "horizontal", "top", "bottom", "inner", "outside" is orientation based on the figure or position
Relationship is set, is merely for convenience of description of the present invention and simplification of the description, rather than the device or element of indication or suggestion meaning are necessary
It with specific orientation, is constructed and operated in a specific orientation, therefore is not considered as limiting the invention.
As shown in figure 8, provided by the invention recognize the in-orbit vibration end displacement of satellite flexible appendage using gyro data
Method comprising the steps of: attitude angular velocity obtaining step: the underdamping after obtaining the in-orbit jet closed-loop control of satellite is freely shaken
The corresponding attitude angular velocity measurement data ω of i axis in dynamic areai(t), wherein i axis is X-axis, the Y-axis, Z in rectangular coordinate system in space
Any one axis in axis；Filter step: to ωi(t) it is filtered and obtains filtered i axis attitude angular velocity data ωi′(t)；
Time series calculates step: according to ωi' (t) calculates the time series η for obtaining mode variablei(t)；End displacement calculates step:
According to ηi(t) with end node formation data diThe in-orbit vibration end displacement m of satellite flexible appendage is calculatedi(t).Preferably,
In attitude angular velocity obtaining step, ωi(t) it is underdamping free vibration period after the in-orbit gas puff Z-pinch of satellite, passes through satellite
What gyro to measure obtained.That is, the gyro that the in-orbit vibration end displacement of identification flexible appendage is loaded merely with satellite platform
Metrical information, without using installed on flexible appendage displacement, speed, acceleration vibrating sensor.
The filter step comprises the steps of: concern vibration frequency band setting steps: the concern vibration frequency of setting i axis direction
Band is [f1,f2]Hz；Strip step: enable bandpass filter to ωi(t) ω is removed in row filtering processingiIt (t) is more than setting value in
Periodic motion and high-frequency vibration component, only retain close vibration injection band in vibration information, obtain ωi′(t).That is,
In strip step, attitude angular velocity measurement data is filtered using bandpass filter, removes the length of satellite body
Periodic motion and other high-frequency vibration error components.Preferably, in concern vibration frequency band setting steps, according to the whole star structure in ground
Frequency band is vibrated in the concern of dynamic analysis calculated result, setting i axis direction.
In embodiment, the method using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage also includes frequency band
Set-up procedure, the frequency band set-up procedure comprise the steps of: global displacement obtaining step: by ηi(t) with mode of oscillation formation
Data DiIt is multiplied, global displacement situation of each node in setting period on acquisition flexible appendage；Principal oscillation mode judgment step: root
According to the global displacement situation, judge whether flexible appendage only has a principal oscillation mode on each direction of vibration；If so,
It continues to execute end displacement and calculates step；It executes filter step if it is not, then returning and adjusts concern vibration frequency band.
It includes first derivative step of converting and integration step that time series, which calculates step,.Wherein, first derivative step of converting:
According to the following formula to ωi' (t) is handled:
In formula:For the first derivative of the corresponding mode variable of i axis；JiFor satellite i axis direction rotary inertia；Bi
For the rotation coefficient of coup of the corresponding flexible appendage mode of oscillation of i axis.Integration step: rightIntegral calculation is carried out, η is obtainedi
(t).Preferably, mi(t) it is obtained by following calculation formula: mi(t)=diηi(t)。
Correspondingly, the in-orbit vibration end displacement of satellite flexible appendage is recognized using gyro data the present invention also provides a kind of
System, comprising with lower module: attitude angular velocity obtains module: the underdamping after obtaining the in-orbit jet closed-loop control of satellite is free
The corresponding attitude angular velocity measurement data ω of the i axis of vibrating areai(t), wherein i axis be rectangular coordinate system in space in X-axis, Y-axis,
Any one axis in Z axis；Filter module: to ωi(t) it is filtered and obtains filtered i axis attitude angular velocity data ωi′(t)；
Time series computing module: according to ωi' (t) calculates the time series η for obtaining mode variablei(t)；End displacement computing module:
According to ηi(t) with end node formation data diThe in-orbit vibration end displacement m of satellite flexible appendage is calculatedi(t).Preferably,
It is obtained in module in attitude angular velocity, ωi(t) it is underdamping free vibration period after the in-orbit gas puff Z-pinch of satellite, passes through satellite
What gyro to measure obtained.That is, the gyro that the in-orbit vibration end displacement of identification flexible appendage is loaded merely with satellite platform
Metrical information, without using installed on flexible appendage displacement, speed, acceleration vibrating sensor.
The filter module includes with lower module: concern vibration frequency band setup module: frequency is vibrated in the concern of setting i axis direction
Band is [f1,f2]Hz；Strip module: enable bandpass filter to ωi(t) ω is removed in row filtering processingiIt (t) is more than setting value in
Periodic motion and high-frequency vibration component, only retain close vibration injection band in vibration information, obtain ωi′(t).That is,
In strip module, attitude angular velocity measurement data is filtered using bandpass filter, removes the length of satellite body
Periodic motion and other high-frequency vibration error components.Preferably, in concern vibration frequency band setup module, according to the whole star structure in ground
Frequency band is vibrated in the concern of dynamic analysis calculated result, setting i axis direction.
In embodiment, the system using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage also includes frequency band
Module is adjusted, the frequency band adjustment module includes with lower module: global displacement obtains module: by ηi(t) with mode of oscillation formation
Data DiIt is multiplied, global displacement situation of each node in setting period on acquisition flexible appendage；Principal oscillation mode judgment module: root
According to the global displacement situation, judge whether flexible appendage only has a principal oscillation mode on each direction of vibration；If so,
Continue to execute end displacement computing module；It executes filter module if it is not, then returning and adjusts concern vibration frequency band.
Time series computing module includes first derivative conversion module and integration module.Wherein, first derivative conversion module:
According to the following formula to ωi' (t) is handled:
In formula:For the first derivative of the corresponding mode variable of i axis；JiFor satellite i axis direction rotary inertia；Bi
For the rotation coefficient of coup of the corresponding flexible appendage mode of oscillation of i axis.Integration module: rightIntegral calculation is carried out, η is obtainedi
(t).Preferably, mi(t) it is obtained by following calculation formula: mi(t)=diηi(t)。
Preferred embodiment:
Certain remote sensing satellite is equipped with single-blade solar battery array, and configuration is shown in Fig. 1.After transmitting is entered the orbit, during gas puff Z-pinch
It is as shown in Figure 2 to roll X-axis attitude angular velocity measurement data.
Step 1: choosing the rolling X-axis attitude angular velocity measurement data ω in satellite period in orbitX(t) it is analyzed
(Fig. 2).
Step 2: according to the whole star structural dynamical model calculated result in ground, the frequency of X-direction principal oscillation mode is about
0.37Hz, the vibration frequency band that concern is arranged is [0.259,0.481] Hz, is existed using the fertile hereby bandpass filter of 5 rank Barts to satellite
Attitude angular velocity measurement data ω run time railX(t) it is filtered, removing rolls in X-axis attitude angular velocity signal
The phugoid mode of motion and high-frequency vibration component only retain the vibration information in concern frequency band, obtain filtered rolling X-axis attitude angle
Speed omegaX′(t)。
Fig. 3 and Fig. 4 are shown in the time domain comparison and power spectral density comparison of attitude angular velocity measurement data filtering front and back respectively.From
Time domain comparison and frequency domain comparison can be clearly seen, and the signal within [0.259,0.481] Hz frequency band is retained, remaining frequency range
Signal is substantially weakened.
Step 3: according to formulaBy filtered single-axis attitude angular velocity omegaX' (t) is converted into mode
The first derivative of variableWherein, JX=6872.13kgm2It is satellite in the rotary inertia for rolling X-direction, BX=
33.48762m·kg1/2For the rotation coefficient of coup of flexible appendage mode of oscillation.Again to the first derivative of mode variableInto
Row integral, obtains the time series η of mode variableX(t), see Fig. 5.
Step 4: by the formation data D of mode of oscillationX(Fig. 6) and mode variable ηX(t) it is multiplied, can be obtained on flexible appendage
Vibration displacement situation of each node in this period.It is made the difference using the vibration displacement of two node any on flexible appendage available
The relative displacement of any two node on flexible appendage.
Step 5: bringing the formation d of end node intoXThe flexible appendage in (Fig. 6) this section of satellite can be obtained period in orbit
End displacement: mX(t)=dXηX(t).In this example, the formation of end node is dX=0.277kg-1/2, obtained end displacement
As a result see Fig. 7.
One skilled in the art will appreciate that in addition to realizing system provided by the invention in a manner of pure computer readable program code
It, completely can be by the way that method and step be carried out programming in logic come so that provided by the invention other than system, device and its modules
System, device and its modules are declined with logic gate, switch, specific integrated circuit, programmable logic controller (PLC) and insertion
The form of controller etc. realizes identical program.So system provided by the invention, device and its modules may be considered that
It is a kind of hardware component, and the knot that the module for realizing various programs for including in it can also be considered as in hardware component
Structure；It can also will be considered as realizing the module of various functions either the software program of implementation method can be Hardware Subdivision again
Structure in part.
Specific embodiments of the present invention are described above.It is to be appreciated that the invention is not limited to above-mentioned
Particular implementation, those skilled in the art can make various deformations or amendments within the scope of the claims, this not shadow
Ring substantive content of the invention.In the absence of conflict, the feature in embodiments herein and embodiment can any phase
Mutually combination.
Claims (10)
1. a kind of method using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage, which is characterized in that include
Following steps:
Attitude angular velocity obtaining step: the i axis in the underdamping free vibration area after obtaining the in-orbit jet closed-loop control of satellite is corresponding
Attitude angular velocity measurement data ωi(t), wherein i axis is X-axis, Y-axis, any of Z axis in rectangular coordinate system in space
Axis；
Filter step: to ωi(t) it is filtered and obtains filtered i axis attitude angular velocity data ω 'i(t)；
Time series calculates step: according to ω 'i(t) the time series η for obtaining mode variable is calculatedi(t)；
End displacement calculates step: according to ηi(t) with end node formation data diThe in-orbit vibration of satellite flexible appendage is calculated
End displacement mi(t)。
2. the method according to claim 1 using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage,
It is characterized in that, in attitude angular velocity obtaining step, ωi(t) when being underdamping free vibration after the in-orbit gas puff Z-pinch of satellite
Phase is obtained by satellite gyroscope measurement.
3. the method according to claim 1 using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage,
It is characterized in that, the filter step comprises the steps of:
Concern vibration frequency band setting steps: the concern vibration frequency band of setting i axis direction is [f1,f2]Hz；
Strip step: enable bandpass filter to ωi(t) ω is removed in row filtering processingiIt (t) is more than the periodic motion of setting value in
With high-frequency vibration component, only retains the vibration information closed in vibration injection band, obtain ω 'i(t)。
4. the method according to claim 3 using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage,
It is characterized in that, also including frequency band set-up procedure, the frequency band set-up procedure is comprised the steps of:
Global displacement obtaining step: by ηi(t) with mode of oscillation formation data DiIt is multiplied, obtains each node on flexible appendage and setting
The global displacement situation of timing phase；
Principal oscillation mode judgment step: according to the global displacement situation, judge flexible appendage on each direction of vibration whether
An only principal oscillation mode；If so, continuing to execute end displacement calculates step；Filter step is executed simultaneously if it is not, then returning
Adjustment concern vibration frequency band.
5. the method according to claim 3 using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage,
It is characterized in that, according to the whole star structural dynamical model calculated result in ground, i axis is arranged in concern vibration frequency band setting steps
Frequency band is vibrated in the concern in direction.
6. the method according to claim 1 using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage,
It is comprised the steps of: it is characterized in that, time series calculates step
First derivative step of converting: according to the following formula to ω 'i(t) it is handled:
In formula:For the first derivative of the corresponding mode variable of i axis；
JiFor satellite i axis direction rotary inertia；
BiFor the rotation coefficient of coup of the corresponding flexible appendage mode of oscillation of i axis；
Integration step: rightIntegral calculation is carried out, η is obtainedi(t)。
7. the method according to claim 1 using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage,
It is characterized in that, mi(t) it is obtained by following calculation formula:
mi(t)=diηi(t)。
8. a kind of system using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage, which is characterized in that include
With lower module:
Attitude angular velocity obtains module: the i axis in the underdamping free vibration area after obtaining the in-orbit jet closed-loop control of satellite is corresponding
Attitude angular velocity measurement data ωi(t), wherein i axis is X-axis, Y-axis, any of Z axis in rectangular coordinate system in space
Axis；
Filter module: to ωi(t) it is filtered and obtains filtered i axis attitude angular velocity data ω 'i(t)；
Time series computing module: according to ω 'i(t) the time series η for obtaining mode variable is calculatedi(t)；
End displacement computing module: according to ηi(t) with end node formation data diThe in-orbit vibration of satellite flexible appendage is calculated
End displacement mi(t)。
9. the system according to claim 8 using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage,
It is characterized in that, the filter module includes with lower module:
Concern vibration frequency band setup module: the concern vibration frequency band of setting i axis direction is [f1,f2]Hz；
Strip module: enable bandpass filter to ωi(t) ω is removed in row filtering processingiIt (t) is more than the periodic motion of setting value in
With high-frequency vibration component, only retains the vibration information closed in vibration injection band, obtain ω 'i(t)；
System using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage also includes frequency band adjustment module, described
It includes with lower module that frequency band, which adjusts module:
Global displacement obtains module: by ηi(t) with mode of oscillation formation data DiIt is multiplied, obtains each node on flexible appendage and setting
The global displacement situation of timing phase；
Principal oscillation mode judgment module: according to the global displacement situation, judge flexible appendage on each direction of vibration whether
An only principal oscillation mode；If so, continuing to execute end displacement computing module；Filter module is executed simultaneously if it is not, then returning
Adjustment concern vibration frequency band.
10. the system according to claim 8 using the in-orbit vibration end displacement of gyro data identification satellite flexible appendage,
It is characterized in that, mi(t) it is obtained by following calculation formula:
mi(t)=diηi(t)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910020517.XA CN109612667B (en) | 2019-01-09 | 2019-01-09 | Method and system for identifying displacement of satellite flexible accessory at rail end by utilizing gyroscope data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910020517.XA CN109612667B (en) | 2019-01-09 | 2019-01-09 | Method and system for identifying displacement of satellite flexible accessory at rail end by utilizing gyroscope data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109612667A true CN109612667A (en) | 2019-04-12 |
| CN109612667B CN109612667B (en) | 2020-10-09 |
Family
ID=66015701
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910020517.XA Active CN109612667B (en) | 2019-01-09 | 2019-01-09 | Method and system for identifying displacement of satellite flexible accessory at rail end by utilizing gyroscope data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109612667B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204214627U (en) * | 2014-11-26 | 2015-03-18 | 武汉理工大学 | Marine propulsion shafting whirling vibration test unit |
| US20150192415A1 (en) * | 2014-01-03 | 2015-07-09 | The Boeing Company | Gyro quadrature stabalization with demodulation phase error nulling |
| CN105092156A (en) * | 2015-09-21 | 2015-11-25 | 上海卫星工程研究所 | Triaxial air bearing table high-precision mass center balance adjusting device and method |
| EP2985509A1 (en) * | 2013-12-17 | 2016-02-17 | Oil Transporting Joint Stock Company "Transneft" (JSC "Transneft") | Device and method for determining position of pipeline |
| CN105446348A (en) * | 2015-12-25 | 2016-03-30 | 北京理工大学 | Distributed control method capable of improving control precision of flexible spacecraft |
| CN106672272A (en) * | 2015-11-10 | 2017-05-17 | 上海新跃仪表厂 | Flexibility parameter in-orbit identification ground testing system for spacecraft with flexible structure |
| CN107194077A (en) * | 2017-05-25 | 2017-09-22 | 北京空间飞行器总体设计部 | A kind of computational methods for obtaining antenna arm vibration suppression response |
| CN107505846A (en) * | 2017-10-12 | 2017-12-22 | 北京航空航天大学 | A kind of anti-interference attitude harmony of Space Manipulator System verifies device and control method |
| CN107792393A (en) * | 2017-09-25 | 2018-03-13 | 上海卫星工程研究所 | The non-contact internal satellite ground checking system of principal and subordinate and its verification method |
| RU2683871C2 (en) * | 2016-02-25 | 2019-04-02 | Алексей Владимирович Молчанов | Method and system of non-destructive monitoring constructions of structures and method and system for monitoring the location of structures based on such a method |
-
2019
- 2019-01-09 CN CN201910020517.XA patent/CN109612667B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2985509A1 (en) * | 2013-12-17 | 2016-02-17 | Oil Transporting Joint Stock Company "Transneft" (JSC "Transneft") | Device and method for determining position of pipeline |
| US20150192415A1 (en) * | 2014-01-03 | 2015-07-09 | The Boeing Company | Gyro quadrature stabalization with demodulation phase error nulling |
| CN204214627U (en) * | 2014-11-26 | 2015-03-18 | 武汉理工大学 | Marine propulsion shafting whirling vibration test unit |
| CN105092156A (en) * | 2015-09-21 | 2015-11-25 | 上海卫星工程研究所 | Triaxial air bearing table high-precision mass center balance adjusting device and method |
| CN106672272A (en) * | 2015-11-10 | 2017-05-17 | 上海新跃仪表厂 | Flexibility parameter in-orbit identification ground testing system for spacecraft with flexible structure |
| CN105446348A (en) * | 2015-12-25 | 2016-03-30 | 北京理工大学 | Distributed control method capable of improving control precision of flexible spacecraft |
| RU2683871C2 (en) * | 2016-02-25 | 2019-04-02 | Алексей Владимирович Молчанов | Method and system of non-destructive monitoring constructions of structures and method and system for monitoring the location of structures based on such a method |
| CN107194077A (en) * | 2017-05-25 | 2017-09-22 | 北京空间飞行器总体设计部 | A kind of computational methods for obtaining antenna arm vibration suppression response |
| CN107792393A (en) * | 2017-09-25 | 2018-03-13 | 上海卫星工程研究所 | The non-contact internal satellite ground checking system of principal and subordinate and its verification method |
| CN107505846A (en) * | 2017-10-12 | 2017-12-22 | 北京航空航天大学 | A kind of anti-interference attitude harmony of Space Manipulator System verifies device and control method |
Non-Patent Citations (2)
| Title |
|---|
| 向明江: ""太阳翼不同转角的卫星在轨模态频率计算方法"", 《航天器工程》 * |
| 吕旺 等: ""风云四号卫星在轨动力学特性分析"", 《第十届全国多体动力学与控制暨第五届全国航天动力学与控制学术会议》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109612667B (en) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105698765B (en) | Object pose method under double IMU monocular visions measurement in a closed series noninertial systems | |
| Lai et al. | Attitude estimation based on fusion of gyroscopes and single antenna GPS for small UAVs under the influence of vibration | |
| CN109682377B (en) | Attitude estimation method based on dynamic step gradient descent | |
| CN107747953B (en) | Multi-sensor data and track information time synchronization method | |
| Fedorov et al. | Using of measuring system MPU6050 for the determination of the angular velocities and linear accelerations | |
| SE464431B (en) | SET AND DEVICE FOR DETERMINING THE LOCATION FOR A BODY | |
| CN109612665A (en) | Recognize the method and system of whole star flexible vibration modal parameter | |
| CN109655218A (en) | With the method and system of the whole star flexible vibration modal frequency of satellite gyroscope data identification | |
| Hoang et al. | A new technique on vibration optimization of industrial inclinometer for MEMS accelerometer without sensor fusion | |
| CN110929386A (en) | Method for identifying whole satellite flexible vibration modal parameters by using satellite gyroscope data | |
| CN112665574A (en) | Underwater robot attitude acquisition method based on momentum gradient descent method | |
| JPS621880B2 (en) | ||
| Carratù et al. | IMU self-alignment in suspensions control system | |
| Hoang et al. | Noise attenuation on IMU measurement for drone balance by sensor fusion | |
| CN109612664A (en) | The in-orbit vibrational state method and system of satellite flexible appendage is recognized using gyro data | |
| CN112345199B (en) | Method for correcting impact of vibration of attack angle sensor of temporary-impulse high-speed wind tunnel | |
| CN109612667A (en) | The in-orbit end displacement method and system of satellite flexible appendage is recognized using gyro data | |
| CN109917797A (en) | Utilize the whole star flexible vibration modal damping method and system of satellite gyroscope data identification | |
| CN109612666A (en) | The in-orbit relative displacement method and system of satellite flexible appendage is recognized using gyro data | |
| CN109916399B (en) | Carrier attitude estimation method under shadow | |
| Hong et al. | Application of EKF for missile attitude estimation based on “SINS/CNS” integrated guidance system | |
| Burdette et al. | Using an array of in-canopy sensors for guided airdrop system state estimation | |
| CN109781103A (en) | Attitude angular velocity measurement method and system based on double super platforms | |
| Yang et al. | Model-free integrated navigation of small fixed-wing UAVs full state estimation in wind disturbance | |
| JP2012198057A (en) | Attitude estimation apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |