CN109612667B - Method and system for identifying displacement of satellite flexible accessory at rail end by utilizing gyroscope data - Google Patents
Method and system for identifying displacement of satellite flexible accessory at rail end by utilizing gyroscope data Download PDFInfo
- Publication number
- CN109612667B CN109612667B CN201910020517.XA CN201910020517A CN109612667B CN 109612667 B CN109612667 B CN 109612667B CN 201910020517 A CN201910020517 A CN 201910020517A CN 109612667 B CN109612667 B CN 109612667B
- Authority
- CN
- China
- Prior art keywords
- vibration
- satellite
- axis
- displacement
- orbit
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Gyroscopes (AREA)
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
The invention provides a method for identifying the displacement of an in-orbit vibration end part of a satellite flexible accessory by utilizing gyro data, which comprises the following steps of: an attitude angular velocity obtaining step: obtaining attitude angular velocity measurement corresponding to i-axis of under-damped free vibration region after satellite on-orbit jet closed-loop controlQuantity data omegai(t); a filtering step: for omegai(t) filtering to obtain filtered i-axis attitude angular velocity data omega'i(t); time series calculation step: according to ω'i(t) calculating η time series of acquired modal variablesi(t) end displacement calculating step according to ηi(t) and end node matrix data diCalculating to obtain the on-orbit vibration end displacement m of the satellite flexible accessoryi(t) of (d). Correspondingly, the invention also provides a system for identifying the displacement of the in-orbit vibration end part of the satellite flexible accessory by utilizing the gyro data. According to the method, other sensors do not need to be installed on the flexible accessory, and only the measurement data of the existing attitude sensor of the satellite platform is utilized for analysis and processing.
Description
Technical Field
The invention relates to the field of satellite monitoring, in particular to a method and a system for identifying the displacement of an in-orbit vibration end part of a satellite flexible accessory by utilizing gyroscope data. In particular to an in-orbit identification method and system for vibration state information of a flexible satellite, and more particularly to an identification method and system for end displacement of the flexible satellite during in-orbit vibration.
Background
With the continuous development of aerospace technology, a large space structure is an important development direction in the aerospace field and also is a necessary infrastructure for space development (in-orbit spacecraft dynamics parameter identification technology research-cloud-climbing). In order to obtain the dynamic characteristics of the flexible satellite structure, the vibration mode parameters are usually calculated by methods such as finite element modeling, ground test, in-orbit identification and the like. Model simplification, condition assumption and the like in the finite element modeling process influence model precision, and particularly accurate modeling of contact mechanisms such as hinges and the like is difficult. Due to the influence of factors such as gravity and atmospheric resistance, large flexible structures are difficult to assemble on the ground and perform full-scale kinetic parameter identification tests. Therefore, for spacecraft with large flexible appendages, it is difficult to obtain precise structural dynamics through finite element modeling or ground testing. Based on the factors, the on-orbit modal parameter identification research on large flexible structures such as solar arrays, space spread antennas and the like is urgent and necessary, and meanwhile, the on-orbit modal parameter identification research method has higher theoretical significance and practical application value. The flexible satellite structure dynamic parameters comprise modal frequency, modal damping, array type, coupling coefficient and the like, have important physical significance, and can provide necessary support for structural design, structural health monitoring, structural fault diagnosis, structural vibration control and other aspects of space flexible parts [ CN102982196A ].
At present, modal parameter identification methods in the field of structural dynamics mainly comprise a frequency domain method, a time domain method and a recently-developed time-frequency domain method. Most methods require on-orbit excitation and arrangement of sensors on a flexible accessory to acquire information, such as [ large solar panel mode parameters on-orbit identification research _ li xiao-master paper ] research on minimum configuration number and optimal layout scheme of sensors on a flexible accessory; CN105486474A introduces a system and method for realizing in-orbit modal identification of satellite flexible components, which needs to perform pulse excitation on a flexible accessory, receive and monitor pulse response signals of each measuring point, and acquire acceleration response signals generated by the satellite flexible components in the steady-state operation process of the satellite in orbit; two methods of satellite solar array sensor layout are introduced [ CN106557633A ] and [ CN107609296A ]. The binocular vision measurement-based on-orbit identification method for the dynamic characteristics of the unfolding structure of the solar sail is characterized in that structural vibration displacement information is directly extracted from an image, and then the dynamic characteristics of the structure are obtained in real time through a working mode analysis technology, so that the on-orbit identification of the dynamic characteristics of the structure is realized. There are also two patents that propose on-orbit identification methods for flexible satellite modal parameters by using attitude angular velocity data measured by a gyroscope: CN103970964A introduces an in-orbit identification method for modal parameters of flexible satellites, which needs to acquire the moment applied to the flexible satellite body by the actuator and the angular velocity information of the flexible satellite body relative to the inertial coordinate system, and obtain the transfer function from the modal parameters and the moment to the angular velocity by using a subspace identification algorithm. CN105157728A proposes a flexible satellite modal parameter identification method capable of suppressing gyro noise influence, which is also to identify the modal frequency and modal damping ratio parameter of the whole satellite by using the measurement data of the satellite body angular velocity when the satellite is in orbit flight, and perform differential processing on the gyro data to suppress the identification error caused by the two noise parts.
The in-orbit kinetic parameter identification method for the flexible satellite has two constraints in the aspect of engineering application: firstly, the on-orbit spacecraft is difficult to apply known excitation required by kinetic parameter identification, and excitation sources are only generated for excitation by unfolding and folding of the on-orbit spacecraft structure, butt joint and separation of the structure, ignition of an engine and the like, and signals of the excitation sources are difficult to measure [ research on the kinetic parameter identification technology of the on-orbit spacecraft-cloud acquisition ]. Secondly, the number of the vibration sensors arranged on the flexible accessory is limited by the implementation of engineering. General flexible accessories need to be unfolded on the rail, and various speed, acceleration and displacement sensors are installed and fixed, cables are laid and the like to have adverse effects on an unfolding mechanism, so that the design difficulty is increased. In addition, the number of slip ring signal channels of the driving mechanism needs to be increased for rotating flexible parts such as a solar array.
Disclosure of Invention
In view of the defects in the prior art, the present invention provides a method and a system for identifying the displacement of the in-orbit vibration end of a satellite flexible accessory by using gyro data.
The method for identifying the in-orbit vibration end displacement of the satellite flexible accessory by utilizing the gyro data comprises the following steps of:
an attitude angular velocity obtaining step: obtaining attitude angular velocity measurement data omega corresponding to i axis of under-damped free vibration area after satellite on-orbit air injection closed-loop controli(t), wherein the i-axis is any one of an X-axis, a Y-axis and a Z-axis in the space rectangular coordinate system;
a filtering step: for omegai(t) filtering to obtain filtered i-axis attitude angular velocity data omegai′(t)；
Time series calculation step: according to omegai' (t) calculating a time series η of derived modal variablesi(t)；
End displacement calculating step according to ηi(t) and end node matrix data diCalculating to obtain the on-orbit vibration end displacement m of the satellite flexible accessoryi(t)。
Preferably, in the attitude angular velocity acquisition step, ωiAnd (t) is obtained by measuring through a satellite gyroscope during the period of under-damped free vibration after the on-orbit air injection control of the satellite.
Preferably, the filtering step comprises the steps of:
a focus vibration band setting step: setting a frequency band of interest in the i-axis direction to [ f1,f2]Hz；
A stripping step: let the band-pass filter pair omegai(t) line filtering, stripping omegaiThe periodic motion and the high-frequency vibration component exceeding the set value in (t) only reserve the vibration information in the vibration injection frequency band to obtain omegai′(t)。
Preferably, the method further comprises a band adjusting step, wherein the band adjusting step comprises the following steps:
overall displacement obtaining step of ηi(t) and vibrational mode matrix data DiMultiplying to obtain the integral displacement condition of each node on the flexible accessory in a set period;
judging a main vibration mode: judging whether the flexible accessory has only one main vibration mode in each vibration direction according to the overall displacement condition; if yes, continuing to execute the end displacement calculation step; if not, returning to the step of executing the filtering and adjusting the concerned vibration frequency band.
Preferably, in the step of setting the vibration frequency band of interest, the vibration frequency band of interest in the i-axis direction is set according to a calculation result of a ground whole star structure dynamics analysis.
Preferably, the time series calculating step comprises the steps of:
a first derivative transformation step: according to the following formula for ωi' (t) treatment:
in the formula:
the first derivative of the modal variable corresponding to the i axis;
Jithe moment of inertia of the satellite in the direction of the i axis;
Bithe rotational coupling coefficient of the vibration mode of the flexible accessory corresponding to the i axis;
an integration step: to pair
Integral calculation is carried out to obtain ηi(t)。
Preferably, mi(t) is obtained by the following calculation formula:
mi(t)＝diηi(t)。
the invention also provides a system for identifying the displacement of the on-orbit vibration end part of the satellite flexible accessory by utilizing the gyro data, which comprises the following modules:
an attitude angular velocity acquisition module: obtaining attitude angular velocity measurement data omega corresponding to i axis of under-damped free vibration area after satellite on-orbit air injection closed-loop controli(t), wherein the i-axis is any one of an X-axis, a Y-axis and a Z-axis in the space rectangular coordinate system;
a filtering module: for omegai(t) filtering to obtain filtered i-axis attitude angular velocity data omegai′(t)；
A time series calculation module: according to omegai' (t) calculating a time series η of derived modal variablesi(t)；
End displacement calculation Module according to ηi(t) and end node matrix data diThe satellite flexible accessory is obtained through calculationDisplacement m of rail vibration endi(t)。
Preferably, the filtering module comprises the following modules:
focus vibration band setting module: setting a frequency band of interest in the i-axis direction to [ f1,f2]Hz；
A stripping module: let the band-pass filter pair omegai(t) line filtering, stripping omegaiThe periodic motion and the high-frequency vibration component exceeding the set value in (t) only reserve the vibration information in the vibration injection frequency band to obtain omegai′(t)；
The system for identifying the on-orbit vibration end displacement of the satellite flexible accessory by utilizing the gyro data also comprises a frequency band adjusting module, wherein the frequency band adjusting module comprises the following modules:
overall displacement acquisition Module ηi(t) and vibrational mode matrix data DiMultiplying to obtain the integral displacement condition of each node on the flexible accessory in a set period;
the main vibration mode judging module: judging whether the flexible accessory has only one main vibration mode in each vibration direction according to the overall displacement condition; if yes, the end displacement calculation module is continuously executed; if not, returning to the execution of the filtering module and adjusting the concerned vibration frequency band.
Preferably, mi(t) is obtained by the following calculation formula:
mi(t)＝diηi(t)。
compared with the prior art, the invention has the following beneficial effects:
1. the method meets the requirement of the on-orbit structure dynamics characteristic identification of the flexible satellite, and extracts the on-orbit vibration end displacement information of the flexible satellite.
2. Compared with other methods for calculating the dynamic parameters of the in-orbit structure of the satellite, the method does not need to install other sensors on the flexible accessory, and only utilizes the measurement data of the existing attitude sensor of the satellite platform to carry out analysis processing.
3. The increased satellite design, manufacturing difficulty and risk of on-orbit operation due to the installation of other vibration sensors are avoided.
4. The identified vibration state information of the flexible accessory can be used for dynamic analysis, fault diagnosis, performance evaluation and the like of the satellite on-orbit structure.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a schematic view of a configuration of a single-wing solar array satellite;
FIG. 2 is a change curve of attitude angular velocity of a satellite during an on-orbit jet closed-loop control period;
FIG. 3 is a time domain comparison graph before and after satellite attitude angular velocity filtering;
FIG. 4 is a frequency domain comparison before and after filtering of satellite attitude angular velocity;
FIG. 5 is a graph of the change in flexural modal variables;
FIG. 6 is a schematic illustration of an array of flexural attachment vibrational modes;
FIG. 7 is a graph of end displacement variation of vibration of a flexible attachment;
FIG. 8 is a flow chart of a method for identifying in-orbit vibration tip displacements of a satellite flexure attachment using gyroscope data.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that variations and modifications can be made by persons skilled in the art without departing from the spirit of the invention. All falling within the scope of the present invention.
In the description of the present invention, it is to be understood that the terms "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc., indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of description and simplicity of description, but do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, are not to be construed as limiting the present invention.
As shown in fig. 8, the method for identifying the displacement of the orbiting vibration end of the satellite flexible accessory by using the gyro data provided by the present invention comprises the following steps: an attitude angular velocity obtaining step: obtaining attitude angular velocity measurement data omega corresponding to i axis of under-damped free vibration area after satellite on-orbit air injection closed-loop controli(t), wherein the i-axis is any one of an X-axis, a Y-axis and a Z-axis in the space rectangular coordinate system; a filtering step: for omegai(t) filtering to obtain filtered i-axis attitude angular velocity data omegai' (t); time series calculation step: according to omegai' (t) calculating a time series η of derived modal variablesi(t) end displacement calculating step according to ηi(t) and end node matrix data diCalculating to obtain the on-orbit vibration end displacement m of the satellite flexible accessoryi(t) of (d). Preferably, in the attitude angular velocity acquisition step, ωiAnd (t) is obtained by measuring through a satellite gyroscope during the period of under-damped free vibration after the on-orbit air injection control of the satellite. That is, the displacement of the vibration end of the flexible attachment in orbit is identified only by using the gyro measurement information loaded on the satellite platform, and displacement, speed and acceleration vibration sensors mounted on the flexible attachment are not used.
The filtering step comprises the steps of: a focus vibration band setting step: setting a frequency band of interest in the i-axis direction to [ f1,f2]Hz; a stripping step: let the band-pass filter pair omegai(t) line filtering, stripping omegaiThe periodic motion and the high-frequency vibration component exceeding the set value in (t) only reserve the vibration information in the vibration injection frequency band to obtain omegai' (t). That is, in the stripping step, the attitude angular velocity measurement data is subjected to filtering processing using a band-pass filter, and long-period motion and other dither error components of the satellite body are stripped. Preferably, in the step of setting the vibration frequency band of interest, the vibration frequency band of interest in the i-axis direction is set according to a calculation result of a ground whole star structure dynamics analysis.
In the embodiment, the number of gyros is usedThe method for identifying the displacement of the on-orbit vibration end of the satellite flexible accessory further comprises a frequency band adjusting step, wherein the frequency band adjusting step comprises the step of acquiring the overall displacement of ηi(t) and vibrational mode matrix data DiMultiplying to obtain the integral displacement condition of each node on the flexible accessory in a set period; judging a main vibration mode: judging whether the flexible accessory has only one main vibration mode in each vibration direction according to the overall displacement condition; if yes, continuing to execute the end displacement calculation step; if not, returning to the step of executing the filtering and adjusting the concerned vibration frequency band.
The time series calculation step comprises a first derivative conversion step and an integration step. Wherein, the first derivative conversion step: according to the following formula for ωi' (t) treatment:
in the formula:
the first derivative of the modal variable corresponding to the i axis; j. the design is a squareiThe moment of inertia of the satellite in the direction of the i axis; b isiThe rotational coupling coefficient of the vibration mode of the flexible attachment corresponding to the i axis. An integration step: to pair
Integral calculation is carried out to obtain ηi(t) of (d). Preferably, mi(t) is obtained by the following calculation formula: m isi(t)＝diηi(t)。
Correspondingly, the invention also provides a system for identifying the in-orbit vibration end displacement of the satellite flexible accessory by utilizing gyro data, which comprises the following modules: an attitude angular velocity acquisition module: obtaining attitude angular velocity measurement data omega corresponding to i axis of under-damped free vibration area after satellite on-orbit air injection closed-loop controli(t), wherein the i-axis is any one of an X-axis, a Y-axis and a Z-axis in the space rectangular coordinate system; a filtering module: for omegai(t) filtering to obtainObtaining filtered i-axis attitude angular velocity data omegai' (t); a time series calculation module: according to omegai' (t) calculating a time series η of derived modal variablesi(t) end displacement calculation Module according to ηi(t) and end node matrix data diCalculating to obtain the on-orbit vibration end displacement m of the satellite flexible accessoryi(t) of (d). Preferably, in the attitude angular velocity acquisition module, ωiAnd (t) is obtained by measuring through a satellite gyroscope during the period of under-damped free vibration after the on-orbit air injection control of the satellite. That is, the displacement of the vibration end of the flexible attachment in orbit is identified only by using the gyro measurement information loaded on the satellite platform, and displacement, speed and acceleration vibration sensors mounted on the flexible attachment are not used.
The filtering module comprises the following modules: focus vibration band setting module: setting a frequency band of interest in the i-axis direction to [ f1,f2]Hz; a stripping module: let the band-pass filter pair omegai(t) line filtering, stripping omegaiThe periodic motion and the high-frequency vibration component exceeding the set value in (t) only reserve the vibration information in the vibration injection frequency band to obtain omegai' (t). That is, in the stripping module, the attitude angular velocity measurement data is filtered by using a band-pass filter, and long-period motion and other high-frequency vibration error components of the satellite body are stripped. Preferably, in the attention vibration frequency band setting module, the attention vibration frequency band in the i-axis direction is set according to a ground whole star structure dynamics analysis calculation result.
In an embodiment, the system for identifying the in-orbit vibration end displacement of the satellite flexible accessory by utilizing the gyro data further comprises a frequency band adjusting module, wherein the frequency band adjusting module comprises a whole displacement acquisition module which is ηi(t) and vibrational mode matrix data DiMultiplying to obtain the integral displacement condition of each node on the flexible accessory in a set period; the main vibration mode judging module: judging whether the flexible accessory has only one main vibration mode in each vibration direction according to the overall displacement condition; if yes, the end displacement calculation module is continuously executed; if not, returning to the execution of the filtering module and adjusting the concerned vibration frequency band.
The time sequence calculation module comprises a first derivative conversion module and an integration module. Wherein the first derivative conversion module: according to the following formula for ωi' (t) treatment:
in the formula:
the first derivative of the modal variable corresponding to the i axis; j. the design is a squareiThe moment of inertia of the satellite in the direction of the i axis; b isiThe rotational coupling coefficient of the vibration mode of the flexible attachment corresponding to the i axis. An integration module: to pair
Integral calculation is carried out to obtain ηi(t) of (d). Preferably, mi(t) is obtained by the following calculation formula: m isi(t)＝diηi(t)。
Preferred embodiments:
a remote sensing satellite is provided with a single-wing solar cell array, and the configuration of the single-wing solar cell array is shown in figure 1. The rolling X-axis attitude angular velocity measurement data during jet control after launch into orbit is shown in fig. 2.
Step 1: selecting rolling X-axis attitude angular velocity measurement data omega of satellite in orbit running periodX(t) analysis was performed (FIG. 2).
Step 2: according to the dynamic analysis and calculation result of the ground whole star structure, the frequency of the main vibration mode in the X direction is about 0.37Hz, and the concerned vibration frequency band is set to be [0.259,0.481]]Hz, using 5 th order Butterworth band-pass filter to measure data omega of satellite attitude angular velocity in orbit running periodX(t) carrying out filtering processing, stripping long-period motion and high-frequency vibration components in the rolling X-axis attitude angular velocity signal, only retaining vibration information in a frequency band of interest, and obtaining the filtered rolling X-axis attitude angular velocity omegaX′(t)。
The time domain contrast and the power spectral density contrast before and after the filtering of the attitude angular velocity measurement data are respectively shown in fig. 3 and fig. 4. As is clear from both time domain contrast and frequency domain contrast, signals within the [0.259,0.481] Hz band are preserved, and signals of the remaining frequency bands are significantly attenuated.
And step 3: according to the formula
The filtered single-axis attitude angular velocity omegaX' (t) conversion to the first derivative of the modal variable
Wherein, JX＝6872.13kg·m2The moment of inertia of the satellite in the direction of the rolling X axis, BX＝33.48762m·kg1/2The coefficient of rotational coupling for the vibrational modes of the flexure attachment. Second derivative of modal variables
Integration is performed to obtain a time series η of modal variablesX(t), see FIG. 5.
And 4, step 4: array data D of vibration modeX(FIG. 6) and modal variables ηXAnd (t) multiplying, and obtaining the vibration displacement condition of each node on the flexible accessory in the period. The relative displacement of any two nodes on the flexible accessory can be obtained by utilizing the difference of the vibration displacement of any two nodes on the flexible accessory.
And 5: matrix d with end nodesX(fig. 6) the end displacement of the flexible attachment during the in-orbit operation of the satellite can be obtained: m isX(t)＝dXηX(t) of (d). In this example, the pattern of end nodes is dX＝0.277kg-1/2The end displacement results obtained are shown in figure 7.
Those skilled in the art will appreciate that, in addition to implementing the systems, apparatus, and various modules thereof provided by the present invention in purely computer readable program code, the same procedures can be implemented entirely by logically programming method steps such that the systems, apparatus, and various modules thereof are provided in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Therefore, the system, the device and the modules thereof provided by the present invention can be considered as a hardware component, and the modules included in the system, the device and the modules thereof for implementing various programs can also be considered as structures in the hardware component; modules for performing various functions may also be considered to be both software programs for performing the methods and structures within hardware components.
The foregoing description of specific embodiments of the present invention has been presented. It is to be understood that the present invention is not limited to the specific embodiments described above, and that various changes and modifications may be made by one skilled in the art within the scope of the appended claims without departing from the spirit of the invention. The embodiments and features of the embodiments of the present application may be combined with each other arbitrarily without conflict.
Claims (6)
1. A method for identifying the displacement of the in-orbit vibration end part of a satellite flexible accessory by utilizing gyro data is characterized by comprising the following steps:
an attitude angular velocity obtaining step: obtaining attitude angular velocity measurement data omega corresponding to i axis of under-damped free vibration area after satellite on-orbit air injection closed-loop controli(t), wherein the i-axis is any one of an X-axis, a Y-axis and a Z-axis in the space rectangular coordinate system;
a filtering step: for omegai(t) filtering to obtain filtered i-axis attitude angular velocity data omega'i(t); time series calculation step: according to ω'i(t) calculating η time series of acquired modal variablesi(t)；
End displacement calculating step according to ηi(t) and end node matrix data diCalculating to obtain the on-orbit vibration end displacement m of the satellite flexible accessoryi(t)；
The filtering step comprises the steps of:
a focus vibration band setting step: setting a frequency band of interest in the i-axis direction to [ f1,f2]Hz；
A stripping step: let the band-pass filter pair omegai(t) performing a filtering process to strip off ωi(t) inThe periodic motion and the dither component exceeding the set value are obtained by retaining only the dither information in the dither band of interest to obtain ω'i(t)；
The time-series calculating step includes the steps of:
a first derivative transformation step: according to the following formula to ω'i(t) carrying out treatment:
in the formula:
the first derivative of the modal variable corresponding to the i axis;
Jithe moment of inertia of the satellite in the direction of the i axis;
Bithe rotational coupling coefficient of the vibration mode of the flexible accessory corresponding to the i axis;
an integration step: to pair
Integral calculation is carried out to obtain ηi(t)；
M in the tip displacement calculating stepi(t) is obtained by the following calculation formula:
mi(t)＝diηi(t)。
2. the method for identifying the in-orbit vibration end displacement of the satellite flexible accessory according to claim 1, wherein in the attitude angular velocity obtaining step, ω isiAnd (t) is obtained by satellite gyro measurement in an underdamped free vibration period after the satellite is subjected to the on-orbit air injection closed-loop control.
3. The method for identifying the end displacement of the in-orbit vibration of the satellite flexible accessory according to the claim 1, further comprising the step of adjusting the frequency band, wherein the step of adjusting the frequency band comprises the steps of:
overall displacement obtaining step of ηi(t) and vibrational mode matrix data DiMultiplying to obtain the integral displacement condition of each node on the flexible accessory in a set period;
judging a main vibration mode: judging whether the flexible accessory has only one main vibration mode in each vibration direction according to the overall displacement condition; if yes, continuing to execute the end displacement calculation step; if not, returning to the step of executing the filtering and adjusting the concerned vibration frequency band.
4. The method for identifying the in-orbit vibration end displacement of the satellite flexible accessory by using the gyro data as claimed in claim 1, wherein in the interested vibration frequency band setting step, the interested vibration frequency band in the i-axis direction is set according to the calculation result of the ground whole satellite structure dynamics analysis.
5. A system for identifying the displacement of the in-orbit vibration end of a satellite flexible accessory by using gyro data is characterized by comprising the following modules:
an attitude angular velocity acquisition module: obtaining attitude angular velocity measurement data omega corresponding to i axis of under-damped free vibration area after satellite on-orbit air injection closed-loop controli(t), wherein the i-axis is any one of an X-axis, a Y-axis and a Z-axis in the space rectangular coordinate system;
a filtering module: for omegai(t) filtering to obtain filtered i-axis attitude angular velocity data omega'i(t); a time series calculation module: according to ω'i(t) calculating η time series of acquired modal variablesi(t) end displacement calculation Module according to ηi(t) and end node matrix data diCalculating to obtain the on-orbit vibration end displacement m of the satellite flexible accessoryi(t)；
The filtering module comprises the following modules:
focus vibration band setting module: setting a frequency band of interest in the i-axis direction to [ f1,f2]Hz；
A stripping module: let the band-pass filter pair omegai(t) filteringTreatment, stripping of omegaiIn the (t), the periodic motion and the dither component exceeding the set value are obtained by retaining only the dither information in the dither band of interest to obtain ω'i(t)；
The time series calculation module includes: a first derivative conversion module and an integration module;
the first derivative conversion module: according to the following formula to ω'i(t) carrying out treatment:
in the formula:
the first derivative of the modal variable corresponding to the i axis;
Jithe moment of inertia of the satellite in the direction of the i axis;
Bithe rotational coupling coefficient of the vibration mode of the flexible accessory corresponding to the i axis;
an integration module: to pair
Integral calculation is carried out to obtain ηi(t)；
M in the end displacement calculation modulei(t) is obtained by the following calculation formula:
mi(t)＝diηi(t)。
6. the system for identifying the in-orbit vibration end displacement of the satellite flexible accessory according to claim 5, further comprising a frequency band adjusting module, wherein the frequency band adjusting module comprises the following modules:
overall displacement acquisition Module ηi(t) and vibrational mode matrix data DiMultiplying to obtain the integral displacement condition of each node on the flexible accessory in a set period;
the main vibration mode judging module: judging whether the flexible accessory has only one main vibration mode in each vibration direction according to the overall displacement condition; if yes, the end displacement calculation module is continuously executed; if not, returning to the execution of the filtering module and adjusting the concerned vibration frequency band.
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 CN109612667A (en) | 2019-04-12 |
| CN109612667B true 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) |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2558724C2 (en) * | 2013-12-17 | 2015-08-10 | Открытое акционерное общество "Акционерная компания по транспорту нефти "Транснефть" (ОАО "АК "Транснефть") | Diagnostic complex for determination of pipeline position, and method for determining relative displacement of pipeline as per results of two and more inspection passes of diagnostic complex for determination of pipelines position |
| US9605964B2 (en) * | 2014-01-03 | 2017-03-28 | 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 |
| CN107194077B (en) * | 2017-05-25 | 2021-02-09 | 北京空间飞行器总体设计部 | Calculation method for obtaining vibration suppression response of antenna arm |
| CN107792393B (en) * | 2017-09-25 | 2020-07-14 | 上海卫星工程研究所 | Master-slave non-contact embedded satellite ground verification system and verification method thereof |
| CN107505846B (en) * | 2017-10-12 | 2018-09-21 | 北京航空航天大学 | A kind of anti-interference attitude harmony verification device of Space Manipulator System and control method |
-
2019
- 2019-01-09 CN CN201910020517.XA patent/CN109612667B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109612667A (en) | 2019-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Reaction wheel disturbance modeling, jitter analysis, and validation tests for solar dynamics observatory | |
| CN110929386A (en) | Method for identifying whole satellite flexible vibration modal parameters by using satellite gyroscope data | |
| CN109612665B (en) | Method and system for identifying flexible vibration modal parameters of whole satellite | |
| JP5474938B2 (en) | Accelerometer and error compensation method | |
| CN109655218B (en) | Method and system for identifying whole-satellite flexible vibration modal frequency by using satellite gyroscope data | |
| CN108827299A (en) | A kind of attitude of flight vehicle calculation method based on improvement quaternary number second order complementary filter | |
| CN106989761B (en) | A kind of spacecraft Guidance instrumentation on-orbit calibration method based on adaptive-filtering | |
| CN106289246A (en) | A kind of rods arm measure method based on position and orientation measurement system | |
| CN108168548B (en) | Pedestrian inertial navigation system and method assisted by machine learning algorithm and model | |
| CN106767846B (en) | Three axis stabilized satellite does not have to the attitude acquisition method and system of gyro | |
| CN109612664B (en) | Method and system for identifying on-orbit vibration state of satellite flexible accessory by utilizing gyroscope data | |
| Hoang et al. | Noise attenuation on IMU measurement for drone balance by sensor fusion | |
| Carratù et al. | IMU self-alignment in suspensions control system | |
| CN109612667B (en) | Method and system for identifying displacement of satellite flexible accessory at rail end by utilizing gyroscope data | |
| Bazin et al. | Feasibility of in-flight quadrotor individual motor thrust measurements | |
| CN109916399A (en) | A kind of attitude of carrier estimation method under shade | |
| CN109612666B (en) | Method and system for identifying on-orbit relative displacement of satellite flexible accessory by utilizing gyroscope data | |
| CN114370930A (en) | Method for measuring satellite structure micro-angle vibration by using MHD micro-angle vibration sensor | |
| CN112326162B (en) | Wing elastic deformation measuring method for airborne distributed POS | |
| Kaswekar et al. | Sensor fusion based vibration estimation using inertial sensors for a complex lightweight structure | |
| Catelani et al. | Reliability and Functional Analysis of IMU systems under temperature-based stress tests | |
| CN109917797A (en) | Utilize the whole star flexible vibration modal damping method and system of satellite gyroscope data identification | |
| Wei et al. | MEMS-based gyro-stellar inertial attitude estimate for NSPO micro-sat program | |
| CN110501029A (en) | The online calibration method of MEMS gyroscope load is carried towards cluster aircraft | |
| CN113984042B (en) | Series combined navigation method applicable to high-dynamic aircraft |
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 |