CN107631738B - Method for identifying imbalance of gyro flywheel couple - Google Patents

Method for identifying imbalance of gyro flywheel couple Download PDF

Info

Publication number
CN107631738B
CN107631738B CN201710851087.7A CN201710851087A CN107631738B CN 107631738 B CN107631738 B CN 107631738B CN 201710851087 A CN201710851087 A CN 201710851087A CN 107631738 B CN107631738 B CN 107631738B
Authority
CN
China
Prior art keywords
rotor
amplitude
unbalance
roll angle
phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710851087.7A
Other languages
Chinese (zh)
Other versions
CN107631738A (en
Inventor
霍鑫
刘海媛
皮永江
郭照升
冯四兆
赵辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201710851087.7A priority Critical patent/CN107631738B/en
Publication of CN107631738A publication Critical patent/CN107631738A/en
Application granted granted Critical
Publication of CN107631738B publication Critical patent/CN107631738B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Balance (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

The invention discloses a method for identifying unbalance of a gyro flywheel couple, and relates to a method for identifying unbalance of a gyro flywheel couple. The method aims to solve the problems that the existing classical influence coefficient method can generate partial non-solution conditions and the difference between different solutions is large. The invention comprises the following steps: firstly, adding even unbalance type trial weights to N phases of a rotor, acquiring rotor roll angle data under the condition of identifying the rotating speed, and obtaining roll angle-frequency multiplication amplitude phi corresponding to the N phases of the trial weights added according to the roll angle datak(ii) a Secondly, the roll angle-frequency multiplication amplitude value corresponding to the N phase adding trial weights obtained in the step one is phikIdentifying original even unbalance of the gyrowheel rotor by adopting a least square estimation algorithm based on an amplitude influence coefficient method or a least square estimation algorithm based on an amplitude and phase influence coefficient method; and thirdly, correcting the even unbalance of the gyro flywheel according to the identification result of the original even unbalance of the gyro flywheel rotor obtained in the second step. The invention is used for the field of dynamic balance of rotating machinery.

Description

Method for identifying imbalance of gyro flywheel couple
Technical Field
The invention relates to the field of dynamic balance of rotating machinery, in particular to a gyro flywheel pair unbalance identification method.
Background
The miniaturization of the spacecraft has become a main hotspot and inevitable trend of space technology development by combining the characteristics of light weight, small volume, low cost, short development period, easy technology updating, high functional density, rapidness, flexibility, good maneuverability and the like. The attitude control system is one of the most important subsystems of the spacecraft, is a main contributor to the mass, the volume, the cost, the power consumption and the like of the tiny spacecraft, and the light weight realization and the multifunctional integration of the attitude control system are important means for promoting the development of the tiny spacecraft.
A gyro flywheel is a novel attitude control and attitude sensing device, integrates the functions of a variable-speed double-frame control moment gyro and a dynamic tuning gyro, and simultaneously realizes three-dimensional moment output and two-dimensional attitude angular rate measurement by controlling the angular momentum amplitude and the direction of a mechanical rotating body. The device improves the integration level and efficiency of the spacecraft attitude control system, reduces the volume, the quality, the power consumption and the cost of the spacecraft attitude control system, and is an important way for realizing the research of the micro spacecraft attitude control technology.
Due to the limitations of the processing technology and the assembly level, the mass distribution of the rotor of the gyro flywheel is not uniform, and vibration can be caused during high-speed rotation. The rotor tilts along the equatorial axis at the same frequency as the rotating speed due to the even unbalance of the gyro flywheel, so that the moment output precision and the attitude sensitive precision of the gyro flywheel are reduced, the normal work of a gyro flywheel system is influenced, and the system is unstable or even damaged in severe cases. Therefore, realizing high-precision dynamic unbalance identification and adjustment is a key technical problem which needs to be solved in gyro flywheel research.
When the gyro flywheel even unbalance is identified by using a classical influence coefficient method (the patent name is: an even unbalance identification method based on gyro flywheel moment coil current, the patent number is 201710556832.5), under the condition of a closed loop of a gyro flywheel tilting control loop, the even unbalance is solved by solving a simultaneous equation set by using a moment coil current signal, but the situation of no solution may exist according to the difference selected by the simultaneous equation set.
Disclosure of Invention
The invention aims to solve the problems that the existing classical influence coefficient method can generate partial non-solution conditions and the difference between different solutions is large, and provides a gyro-flywheel even imbalance identification method.
A gyro flywheel pair unbalance identification method comprises the following steps:
step one, adding even unbalance type trial weights to N phases of a rotor, acquiring rotor roll angle data under the condition of identifying the rotating speed, and obtaining roll angle-frequency multiplication amplitude phi corresponding to the N phases of the trial weights added according to the roll angle datak
Step two, obtaining the product by utilizing the step oneThe tilt angle-frequency multiplication amplitude corresponding to the N phase adding trial weights is phikIdentifying original even unbalance of the gyrowheel rotor by adopting a least square estimation algorithm based on an amplitude influence coefficient method or a least square estimation algorithm based on an amplitude and phase influence coefficient method;
and step three, correcting the even unbalance of the gyro flywheel according to the original even unbalance identification result of the gyro flywheel rotor obtained in the step two.
The invention has the beneficial effects that:
the method can realize identification and correction of even unbalance of the gyro flywheel, effectively improve the mass distribution of the gyro flywheel rotor, reduce the first frequency multiplication amplitude of the tilt angle of the gyro flywheel rotor from 0.0232 degrees to 0.0022 degrees, and has obvious correction effect when the even unbalance correction amplitude of the gyro flywheel is 90.4 percent, as shown in figure 4. Therefore, moment noise of the gyro flywheel can be effectively reduced, and the running stability of the rotor is improved, so that the moment output precision and the attitude measurement precision of the rotor are improved.
Aiming at the defects of a classical influence coefficient method (the patent name is: an even unbalance identification method based on gyro flywheel moment coil current), the invention adopts a roll angle signal under the open-loop condition of a roll control loop to directly reflect the vibration caused by even unbalance, but not the moment coil current under the closed-loop condition of the roll control loop to indirectly reflect the vibration caused by even unbalance; meanwhile, the existence and reasonableness of the knowledge are guaranteed by using an improved influence coefficient method based on least square.
According to the invention, even unbalance calculation is carried out by linear least square fitting by utilizing a rotor tilt angle signal and adopting an improved influence coefficient method based on least square under the open loop condition of a gyro flywheel tilt control loop, so that the existence of the solution can be ensured. The effect of identifying the original even unbalance of the gyro flywheel by using the improved influence coefficient method based on least square is shown in the table I, and the table I can know that the identification can be realized by using more than three groups of data and the identification result has little difference.
The effect pair of the original even unbalance of the gyrorotor is identified by using an improved influence coefficient method (the invention) based on least square and a classical influence coefficient method (the patent name is: an even unbalance identification method based on gyrorotor moment coil current), for example, as shown in the second table, the second table shows that the identification result of the improved influence coefficient method based on least square is stable, and no solution exists.
In conclusion, the improved influence coefficient method based on the least square is more suitable for identifying the even unbalance of the gyro flywheel.
TABLE I improved influence coefficient method identification result (3600r/min) based on least squares
Figure BDA0001413588690000021
Figure BDA0001413588690000031
TABLE II comparison of identification results (3600r/min) using two algorithms
Figure BDA0001413588690000032
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a diagram showing the repeatability of the gyro flywheel tilt angle-frequency doubling amplitude varying with the rotation speed;
FIG. 3 is a graph showing the variation of the roll angle-doubling amplitude with the trial weight added phase for different trial weight phases;
FIG. 4 is a comparison graph of a gyro flywheel rotor tilt angle-frequency multiplication amplitude before and after identification and correction by using the identification algorithm of the present invention under the open loop condition of the tilt control loop.
Detailed Description
The first embodiment is as follows: as shown in fig. 1, a method for identifying an imbalance of a gyro-flywheel couple includes the following steps:
the specific process of determining the identification rotating speed of the even unbalance amount by the early preparation work is as follows:
1) selecting a test rotating speed of omega within the range of the working rotating speed of the gyro flywheeli,i=1,2,…,n;
2) When testing the rotation speed omegaiAfter the operation is stable, the gyro flywheel works in an open loop state of a tilting control loop, the data of the tilting angle of a rotor of the gyro flywheel within 10 seconds are recorded, and the step is repeatedly executed for 5 times;
3) calculating the roll angle data acquired in the step 2) by adopting fast Fourier transform to obtain a test rotating speed omega for 5 timesiCorresponding roll angle-frequency multiplication amplitude is phiijJ-1, 2, …,5, as shown in fig. 2.
4) According to 5 times of tilting angle, namely frequency multiplication amplitude phiijMethod for determining identification rotating speed omega of gyro flywheel rotor along with repeatability of rotating speed change curveT
Step one, adding even unbalance type trial weights to N phases of a rotor, acquiring rotor roll angle data under the condition of identifying the rotating speed, and obtaining roll angle-frequency multiplication amplitude phi corresponding to the N phases of the trial weights added according to the roll angle datak
Step two, the roll angle-frequency multiplication amplitude value corresponding to the N phase adding trial weights obtained in the step one is phikIdentifying original even unbalance of the gyrowheel rotor by adopting a least square estimation algorithm based on an amplitude influence coefficient method or a least square estimation algorithm based on an amplitude and phase influence coefficient method;
and step three, correcting the even unbalance of the gyro flywheel according to the original even unbalance identification result of the gyro flywheel rotor obtained in the step two.
The method can realize identification and correction of even unbalance of the gyro flywheel, effectively improve the mass distribution of the gyro flywheel rotor, reduce the first frequency multiplication amplitude of the tilt angle of the gyro flywheel rotor from 0.0232 degrees to 0.0022 degrees, and has obvious correction effect when the even unbalance correction amplitude of the gyro flywheel is 90.4 percent, as shown in figure 4. Therefore, moment noise of the gyro flywheel can be effectively reduced, and the running stability of the rotor is improved, so that the moment output precision and the attitude measurement precision of the rotor are improved.
Aiming at the defects of a classical influence coefficient method (the patent name is: an even unbalance identification method based on gyro flywheel moment coil current), the invention adopts a roll angle signal under the open-loop condition of a roll control loop to directly reflect the vibration caused by even unbalance, but not the moment coil current under the closed-loop condition of the roll control loop to indirectly reflect the vibration caused by even unbalance; meanwhile, the existence and reasonableness of the knowledge are guaranteed by using an improved influence coefficient method based on least square.
The second embodiment is as follows: the first difference between the present embodiment and the specific embodiment is: in the first step, even unbalance type trial weights are added to N phases of the rotor, rotor roll angle data are collected under the condition of identifying the rotating speed, and roll angle-frequency multiplication amplitude phi corresponding to the N phases of the phase adding trial weights is obtained according to the roll angle datakThe specific process comprises the following steps:
one step, when no test weight is added, the gyro flywheel rotor is enabled to identify the rotating speed omegaTStably running, recording gyroscopic flywheel rotor roll angle data within 10 seconds under the open loop condition of a roll control loop, and calculating the rotor roll angle by adopting fast Fourier transform to obtain an initial roll angle-frequency multiplication amplitude phi0
Step one, two, in the phase theta of the gyro flywheel rotorkWhere k is 1,2, …, N, N is not less than 3, N is the number of positions to which trial weights are added, and the added mass is mkIs in even unbalance form (namely, the mass m is respectively added at the symmetrical positions of the upper surface and the lower surface of the rotorkTrial weight of/2);
step three, when the gyro flywheel rotor identifies the rotating speed omegaTWhen the gyroscope flywheel rotor stably runs, recording gyroscopic flywheel rotor roll angle data within 10 seconds under the open-loop condition of a roll control loop;
step four, calculating the inclination angle of the rotor by adopting fast Fourier transform to obtain the inclination angle-frequency multiplication amplitude phi corresponding to N phase adding trial weightsk
Other steps and parameters are the same as those in the first embodiment.
The third concrete implementation mode: the present embodiment differs from the first or second embodiment in that: what is needed isIn the second step, the N phase adding trial weights obtained in the first step are used, and the corresponding roll angle-frequency multiplication amplitude is phikThe specific process of identifying the original even unbalance of the gyrorotor rotor by adopting a least square estimation algorithm based on an amplitude influence coefficient method or a least square estimation algorithm based on an amplitude and phase influence coefficient method is as follows:
step two, firstly: the phase and amplitude of the original even unbalance of the gyro flywheel are assumed to be theta0And U0The equivalent mass of the original even unbalance amount on the test weight surface (the plane added with the test weight) is m0. Based on an influence coefficient method, the relation between the roll angle and the even unbalance when the N phases of the rotor are subjected to the test weight in the form of adding the even unbalance is expressed in a complex form as follows:
Figure BDA0001413588690000051
wherein C is1,C2For the influence coefficient, theta, of the gyro flywheel prototype to be correctedC1C2Are respectively C1And C2A corresponding phase; u shapekFor added even unbalance, phikAnd thetaφkRepresenting the roll angle-frequency multiplication amplitude and phase when the test weight in the form of even unbalance is added to the N phases; e denotes the base of the natural logarithm, i denotes the unit of an imaginary number, eixX is θ, cosx + isinxφkC1UC20Or thetak
Rewriting formula (1) to the following form:
Figure BDA0001413588690000052
the formula (2) is arranged as Axk=bkIn the form of (a), wherein,
Figure BDA0001413588690000053
Figure BDA0001413588690000054
Figure BDA0001413588690000055
wherein A is a parameter matrix, xk、bkIs an intermediate variable;
according to least squares, by sum of squares of residuals
Figure BDA0001413588690000061
The minimum is the standard, the parameter matrix A can be estimated, and the inclination angle can meet the requirement when no test weight is added
Figure BDA0001413588690000062
An optimal solution of the original even imbalance can be estimated with amplitude and phase m respectively0And theta0
Step two: the specific process of the least square estimation algorithm based on the amplitude influence coefficient method or the least square estimation algorithm based on the amplitude and phase influence coefficient method is as follows:
(1) the least square estimation algorithm based on the amplitude influence coefficient method comprises the following steps: that is, equation (1) is expressed in scalar form as:
Figure BDA0001413588690000063
the sum of the squared residuals is then expressed as:
Figure BDA0001413588690000064
according to the least square method, more than or equal to three groups of experimental data mk,θkAnd phikCarry-in (7) to identify m with minimum sum of squared residuals as the criterion0And theta0
(2) The least square estimation algorithm based on the amplitude and phase influence coefficient method is as follows: that is, equation (1) is expressed in scalar form as:
Figure BDA0001413588690000065
wherein theta ismTo add the mth phase, θ, of the trial weight in even-unbalanced formnTo add the nth phase of the trial weight in even unbalance form, k, m, N is 1,2, …, N ≧ 3, the sum of the squares of the residuals is expressed as:
Figure BDA0001413588690000066
wherein the content of the first and second substances,
Figure BDA0001413588690000067
wherein EkError of the system of equations, Ek TIs EkTransposing;
according to the least square method, the experimental data phi when no test weight is added0And thetaφ0And experimental data m of two groups or morek,θk,φk,θφkWith equations (9) and (10), the parameters are estimated using the residual sum of squares, J, as a criterion
Figure BDA0001413588690000071
According to
Figure BDA0001413588690000072
Identification m0And theta0(estimate the optimal solution of the original even imbalance amplitude and phase, i.e., m0And theta0)。
Other steps and parameters are the same as those in the first or second embodiment.
The fourth concrete implementation mode: the difference between this embodiment mode and one of the first to third embodiment modes is: the specific process of correcting the even unbalance of the gyro flywheel according to the identification result of the original even unbalance of the gyro flywheel rotor obtained in the step two in the step three is as follows:
phase θ on trial weight surface0The added mass at the +180 DEG position is m0And balancing the balance weight in an even unbalance mode to finish the correction. Enabling the gyro flywheel rotor to identify the rotation speed omegaTAnd (4) stably running, recording a tilting angle-frequency multiplication amplitude of the gyro flywheel within 10 seconds under the condition of open loop of a tilting control loop, and verifying the even unbalance identification and correction effect by comparing the tilting angle-frequency multiplication amplitude with the correction forward tilting angle-frequency multiplication amplitude.
Other steps and parameters are the same as those in one of the first to third embodiments.
The following examples were used to demonstrate the beneficial effects of the present invention:
the first embodiment is as follows:
the method comprises the following steps: even unbalance type test weights with the mass of 0.5g are added at the positions of 0 degree, 60 degrees, 120 degrees, 180 degrees, 240 degrees, 300 degrees and 360 degrees of a test weight plane of the gyrorotor flywheel rotor respectively. And when the gyroduower rotor stably runs at the identification rotating speed of 3600r/min, recording the roll angle data of the gyroduower rotor within 10 seconds under the open-loop condition of the roll control loop. Calculating the tilt angle of the rotor by adopting fast Fourier transform to obtain tilt angle-frequency multiplication amplitudes phi corresponding to different trial weight phases1234567. The variation of the rotor roll angle-doubled frequency amplitude with the phase of the added trial weight is shown in fig. 3.
Step two: according to the least square estimation algorithm based on the amplitude, the relationship between the tilt angle of the gyro flywheel rotor and the even imbalance under different test weight conditions expressed in a scalar form is shown as a formula (6), the sum of squares of the residual errors is shown as a formula (7), the minimum sum of squares of the residual errors is used as a standard, seven groups of experimental data are used for obtaining that the phase of the original even imbalance is 265 degrees, and the equivalent mass of the original even imbalance on the test weight surface is 0.518 g.
Step three: the correction is completed by adding a weight in the form of an even unbalance with a mass of 0.518g at phase 85 °. After correction, when the gyrowheel rotor stably runs at the rotating speed of 3600r/min, the data of the roll angle of the gyrowheel rotor within 10 seconds under the open-loop condition of the roll control loop are recorded, the roll angle of the rotor is calculated by adopting fast Fourier transform, the frequency doubling amplitude of the roll angle before and after correction is compared, and the even unbalance identification and correction effect is verified, as shown in fig. 4, the frequency doubling amplitude of the roll angle is reduced from 0.0232 degrees to 0.0022 degrees.
The present invention is capable of other embodiments and its several details are capable of modifications in various obvious respects, all without departing from the spirit and scope of the present invention.

Claims (2)

1. A gyro flywheel pair imbalance identification method is characterized by comprising the following steps: the gyro flywheel pair unbalance identification method comprises the following steps:
step one, adding even unbalance type trial weights to N phases of a rotor, acquiring rotor roll angle data under the condition of identifying the rotating speed, and obtaining roll angle-frequency multiplication amplitude phi corresponding to the N phases of the trial weights added according to the roll angle datak(ii) a The specific process is as follows:
one step, when no test weight is added, the gyro flywheel rotor is enabled to identify the rotating speed omegaTStably running, recording gyroscopic flywheel rotor roll angle data within 10 seconds under the open loop condition of a roll control loop, and calculating the rotor roll angle by adopting fast Fourier transform to obtain an initial roll angle-frequency multiplication amplitude phi0
Step one, two, in the phase theta of the gyro flywheel rotorkWhere k is 1,2, …, N, N is not less than 3, and the added mass is mkA trial weight in the form of an even imbalance of (c);
step three, when the gyro flywheel rotor identifies the rotating speed omegaTWhen the gyroscope flywheel rotor stably runs, recording gyroscopic flywheel rotor roll angle data within 10 seconds under the open-loop condition of a roll control loop;
step four, calculating the inclination angle of the rotor by adopting fast Fourier transform to obtain the inclination angle-frequency multiplication amplitude phi corresponding to N phase adding trial weightsk
Step two, the roll angle-frequency multiplication amplitude value corresponding to the N phase adding trial weights obtained in the step one is phikUsing a baseIdentifying original even unbalance of the gyro flywheel rotor by a least square estimation algorithm of an amplitude influence coefficient method or a least square estimation algorithm based on the amplitude and phase influence coefficient method; the specific process is as follows:
step two, firstly: the phase and amplitude of the original even unbalance of the gyro flywheel are assumed to be theta0And U0The equivalent mass of the original even unbalance on the test weight surface is m0(ii) a Based on an influence coefficient method, the relation between the roll angle and the even unbalance when the N phases of the rotor are subjected to the test weight in the form of adding the even unbalance is expressed in a complex form as follows:
Figure FDA0002767048180000011
wherein C is1,C2For the influence coefficient, theta, of the gyro flywheel prototype to be correctedC1C2Are respectively C1And C2A corresponding phase; u shapekFor added even unbalance, phikAnd thetaφkRepresenting the roll angle-frequency multiplication amplitude and phase when the test weight in the form of even unbalance is added to the N phases;
rewriting formula (1) to the following form:
Figure FDA0002767048180000021
the formula (2) is arranged as Axk=bkIn the form of (a), wherein,
Figure FDA0002767048180000022
Figure FDA0002767048180000023
Figure FDA0002767048180000024
wherein A is a parameter matrix, xk、bkIs an intermediate variable;
step two: the specific process of the least square estimation algorithm based on the amplitude influence coefficient method or the least square estimation algorithm based on the amplitude and phase influence coefficient method is as follows:
(1) the least square estimation algorithm based on the amplitude influence coefficient method comprises the following steps: that is, equation (1) is expressed in scalar form as:
Figure FDA0002767048180000025
the sum of the squared residuals is then expressed as:
Figure FDA0002767048180000026
more than or equal to three groups of experimental data mk,θkAnd phikSubstituting into equation (7), and identifying m with minimum sum of squared residuals as the standard0And theta0
(2) The least square estimation algorithm based on the amplitude and phase influence coefficient method is as follows: that is, equation (1) is expressed in scalar form as:
Figure FDA0002767048180000027
wherein theta ismTo add the mth phase, θ, of the trial weight in even-unbalanced formnTo add the nth phase of the trial weight in even imbalance, the sum of the squared residuals is expressed as:
Figure FDA0002767048180000028
wherein the content of the first and second substances,
Figure FDA0002767048180000031
wherein EkError of the system of equations, Ek TIs EkTransposing;
the experimental data phi when no test weight is added0And thetaφ0And more than or equal to two sets of experimental data mk,θk,φk,θφkWith equations (9) and (10), the parameters are estimated using the residual sum of squares, J, as a criterion
Figure FDA0002767048180000032
According to
Figure FDA0002767048180000033
Identification m0And theta0
And step three, correcting the even unbalance of the gyro flywheel according to the original even unbalance identification result of the gyro flywheel rotor obtained in the step two.
2. The method for identifying the imbalance of the gyrometer according to claim 1, wherein: the specific process of correcting the even unbalance of the gyro flywheel according to the identification result of the original even unbalance of the gyro flywheel rotor obtained in the step two in the step three is as follows:
phase θ on trial weight surface0The added mass at the +180 DEG position is m0And balancing the balance weight in an even unbalance mode to finish the correction.
CN201710851087.7A 2017-09-20 2017-09-20 Method for identifying imbalance of gyro flywheel couple Active CN107631738B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710851087.7A CN107631738B (en) 2017-09-20 2017-09-20 Method for identifying imbalance of gyro flywheel couple

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710851087.7A CN107631738B (en) 2017-09-20 2017-09-20 Method for identifying imbalance of gyro flywheel couple

Publications (2)

Publication Number Publication Date
CN107631738A CN107631738A (en) 2018-01-26
CN107631738B true CN107631738B (en) 2021-02-09

Family

ID=61103028

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710851087.7A Active CN107631738B (en) 2017-09-20 2017-09-20 Method for identifying imbalance of gyro flywheel couple

Country Status (1)

Country Link
CN (1) CN107631738B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109143861B (en) * 2018-09-07 2021-06-29 哈尔滨工业大学 Active vibration suppression method for variable speed tilting momentum wheel system based on torquer
CN109211215B (en) * 2018-10-26 2022-03-18 哈尔滨工业大学 Three-degree-of-freedom flexible support rotor tilting vibration control method
CN110702315A (en) * 2019-06-06 2020-01-17 上海航天控制技术研究所 On-orbit estimation method and device for static imbalance characteristic of large-inertia rotating part
CN111578966B (en) * 2020-04-09 2021-07-16 哈尔滨工程大学 Hemisphere harmonic oscillator characteristic parameter identification method based on LMS algorithm
CN113029436B (en) * 2021-03-01 2023-10-20 大唐东北电力试验研究院有限公司 Method for dynamically balancing long shaft of three-cylinder two-exhaust 200MW steam turbine generator assembly
CN113156987B (en) * 2021-04-15 2022-05-31 哈尔滨工业大学 Spacecraft actuating mechanism combining double-frame shear type moment gyroscope and flywheel and control method thereof

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103424225A (en) * 2013-07-26 2013-12-04 北京控制工程研究所 Method for measuring dynamic and static unbalance of rotating component
CN105478245A (en) * 2015-12-30 2016-04-13 哈尔滨工业大学 Identification method on basis of vibration detection of main shaft and for dynamic unbalance value of auxiliary shaft of double-freedom degree precision centrifugal machine

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103424225A (en) * 2013-07-26 2013-12-04 北京控制工程研究所 Method for measuring dynamic and static unbalance of rotating component
CN105478245A (en) * 2015-12-30 2016-04-13 哈尔滨工业大学 Identification method on basis of vibration detection of main shaft and for dynamic unbalance value of auxiliary shaft of double-freedom degree precision centrifugal machine

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
三自由度陀螺飞轮转子动不平衡辨识方法研究;郭照升;《中国优秀硕士学位论文全文数据库 工程科技第II辑2017年第2期》;20170215;第25-44页 *

Also Published As

Publication number Publication date
CN107631738A (en) 2018-01-26

Similar Documents

Publication Publication Date Title
CN107631738B (en) Method for identifying imbalance of gyro flywheel couple
JP5084419B2 (en) Compass sensor calibration method, apparatus and medium in consideration of magnetic environment, and azimuth measuring method, apparatus and medium using the same
CN107356265B (en) Method for identifying even unbalance amount based on gyro flywheel torque coil current
CN109163734B (en) Autonomous calibration method based on dual-axis optical fiber rotation modulation integrated navigation system
CN107992063B (en) Variable-parameter nutation damping-based precession control method for variable-speed tilting momentum wheel
CN106441357B (en) A kind of single-shaft-rotation SINS axial direction gyroscopic drift bearing calibration based on damping network
CN115876182B (en) Electrode error modeling method of hemispherical resonator gyroscope
CN106017452B (en) Double tops disturbance rejection north finding method
CN115574798B (en) Method for improving comprehensive performance of hemispherical resonator gyroscope
CN104697553B (en) Fiber-optic gyroscope strapdown inertial navigation system accelerometer interior bar arm calibration method
CN114166242B (en) Calibration method and system for hemispherical resonator gyro detection signal non-uniformity
CN110806220B (en) Inertial navigation system initial alignment method and device
CN101655371A (en) Position signal damping method of inertial navigation system based on variable damping coefficient
Lv et al. A method of low-cost IMU calibration and alignment
CN113503894A (en) Inertial navigation system error calibration method based on gyroscope reference coordinate system
CN106802143A (en) A kind of hull deformation angle measuring method based on inertial instruments and Iterative-Filtering Scheme
WO2022222938A1 (en) Motion state monitoring-based adaptive horizontal attitude measurement method
US10310132B2 (en) Absolute vector gravimeter and methods of measuring an absolute gravity vector
CN107270899A (en) Long endurance inertial navigation system damping switching method based on switching control
EP2667155B1 (en) Inertial unit with vibratory gyroscopes mounted on a carousel and angle measuring method
CN109143861B (en) Active vibration suppression method for variable speed tilting momentum wheel system based on torquer
CN115979311B (en) PIGA cross quadratic term coefficient calibration method, system, equipment and medium
CN112747772A (en) Request-based inertial odometer moving base coarse alignment method
CN116698084A (en) Real-time zero-bias drift suppression and zero-bias stability improvement method for symmetrical gyroscope
CN111609869A (en) Positive and negative multi-position fiber-optic gyroscope orientation effect judgment method based on hypothesis testing

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