CN103424225B - A kind of method of testing rotatable parts sound amount of unbalance - Google Patents

A kind of method of testing rotatable parts sound amount of unbalance Download PDF

Info

Publication number
CN103424225B
CN103424225B CN201310319876.8A CN201310319876A CN103424225B CN 103424225 B CN103424225 B CN 103424225B CN 201310319876 A CN201310319876 A CN 201310319876A CN 103424225 B CN103424225 B CN 103424225B
Authority
CN
China
Prior art keywords
air floating
floating table
unbalance
axis
rotating part
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
CN201310319876.8A
Other languages
Chinese (zh)
Other versions
CN103424225A (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.)
Beijing Institute of Control Engineering
Original Assignee
Beijing Institute of Control Engineering
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 Beijing Institute of Control Engineering filed Critical Beijing Institute of Control Engineering
Priority to CN201310319876.8A priority Critical patent/CN103424225B/en
Publication of CN103424225A publication Critical patent/CN103424225A/en
Application granted granted Critical
Publication of CN103424225B publication Critical patent/CN103424225B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Of Balance (AREA)

Abstract

The invention discloses a kind of method of testing rotatable parts sound amount of unbalance, the present invention is by the test to spacecraft rotating part sound amount of unbalance, obtain quiet unbalancing value during spacecraft rotatable parts operation on orbit, on the basis of test result by counterweight to reduce the size of quiet unbalancing value, during making spacecraft rotatable parts operation on orbit, the interference of celestial body is met and allow requirement, the present invention has directly applied to many satellites with flexible rotatable parts, show that the method successfully solves the quiet unbalance dynamic test problem of flexible rotatable parts by data in-orbit.The method is simple, and measuring accuracy is high.

Description

A kind of method of testing rotatable parts sound amount of unbalance
Technical field
The present invention relates to a kind of method of testing rotatable parts sound amount of unbalance, particularly relate to the method for the spacecraft sound amount of unbalance particularly tested with flexible rotatable parts, belong to field of measuring technique.
Background technology
Generally, with the spacecraft of rotatable parts, when the quiet unbalancing value of rotatable parts is larger, the disturbance torque produced will affect the attitude in-orbit of spacecraft, the work of other useful load spaceborne is caused to be affected, need the size measuring the quiet unbalancing value of rotatable parts on ground comparatively accurately for this reason, and reduce amount of unbalance by counterweight as far as possible.In order to reduce the quiet unbalancing value of rotatable parts, key is the test to quiet unbalancing value, and this test generally adopts dynamic balancing machine to carry out on ground.In order to improve the precision of test, tested equipment must High Rotation Speed, and this and rotatable parts are real work operating mode in-orbit inconsistent (during the operation on orbit of rotatable parts, velocity of rotation is often not high) after launching.Therefore, said method exists following not enough: said method is only adapted to the good tested parts of rigidity, for the equipment that there is more obvious flexible deformation's characteristic, because distortion is larger with the relation of rotating speed, the state that high rotating speed during test makes the size of tested parts generation distortion and operation on orbit be is inconsistent, what finally cause test result to reflect is the quiet unbalancing value of tested parts under larger rotating speed, and can not reflect the amount of unbalance under tested parts actual condition.
Summary of the invention
The technical matters that the present invention solves is: overcome the deficiencies in the prior art, propose a kind of method of testing rotatable parts sound amount of unbalance, and the method is simple, and measuring accuracy is high.
Technical scheme of the present invention is: a kind of method of testing rotatable parts sound amount of unbalance, and step is as follows:
(1) step of air floating table X, Y-axis moment of inertia is measured: float air floating table, to air floating table coarse balance, the counterweight of known quality is placed in a certain position on air floating table, the gravitational torque utilizing counterweight to produce makes air floating table produce angular acceleration, adopt the stage body angular velocity of gyro to measure air floating table, the stage body angular velocity of gyro to measure is carried out curve fitting and obtains the stage body angular acceleration of air floating table, the angular acceleration that the interference of deduction air floating table self produces, determines the moment of inertia of air floating table X, Y-axis according to theorem of angular momentum;
(2) step of tested parts rotating part static-unbalance is measured: tested parts rotating part is placed in respectively the x-axis of air floating table and negative x-axis direction, float air floating table, to air floating table coarse balance tested parts are not rotated, the gravitational torque utilizing tested parts to produce makes air floating table produce angular acceleration, adopt the stage body angular velocity of gyro to measure air floating table, the stage body angular velocity of gyro to measure is carried out curve fitting and obtains the stage body angular acceleration of air floating table, according to the air floating table X that theorem of angular momentum integrating step (1) is determined, the moment of inertia of Y-axis calculates the static-unbalance of tested parts rotating part,
(3) step of the disturbance torque that tested parts unbalance dynamic produces is measured: be placed on air floating table by tested parts rotating part, float air floating table, to air floating table coarse balance, the rotating part of tested parts is rotated, under rotation condition, adopt stage body attitude angle and the angular velocity of gyro to measure air floating table, the air floating table stage body attitude angle utilizing wave filter to go out gyro to measure and angular velocity carry out filtering and spectrum analysis and obtain the disturbance torque that tested parts unbalance dynamic produces;
(4) counterweight of certain mass is installed on tested parts rotating part, the disturbance torque of the tested parts unbalance dynamic generation after installing counterweight is obtained according to the measuring process of step (3), the disturbance torque relatively installing tested parts unbalance dynamic generation before and after counterweight obtains the disturbance torque phase place of counterweight generation, calibrate the tested parts rotation phase that disturbance torque data initial time is corresponding, thus obtain the disturbance torque of calibrated test component unbalance dynamic generation;
(5) disturbance torque that the test component unbalance dynamic that the tested parts static-unbalance utilizing step (2) to measure to obtain, step (4) calibrate produces calculates the couple-unbalance of tested parts rotating part in conjunction with the installation site of tested parts on air floating table.
The described method utilizing theorem of angular momentum to determine the moment of inertia of air floating table X, Y-axis:
Wherein J is the moment of inertia of air floating table X-axis or Y-axis, and m is counterbalance mass, and g is test ground local gravitational acceleration, and l is the horizontal range of counterweight center and air floating table center of rotation, the air floating table angular acceleration obtained is measured after increasing counterweight, for the air floating table angular acceleration that measurement when not increasing counterweight obtains, i.e. the angular acceleration that produces of air floating table self interference.
The method of described calculating tested parts rotating part static-unbalance is:
mr x = J x ( ω . x 0 - ω . x 180 ) 2 g
mr y = J y ( ω . y 0 - ω . y 180 ) 2 g
mr = mr x 2 + mr y 2
α = ctg ( mr y mr x )
Wherein:
Mr xfor the X-axis component of static-unbalance under tested part coordinates system;
Mr yfor the Y-axis component of static-unbalance under tested part coordinates system;
Mr is tested parts rotating part static-unbalance;
J xfor the moment of inertia of air floating table X-axis;
J yfor the moment of inertia of air floating table Y-axis;
for the air floating table stage body angular acceleration that measurement when tested parts rotating part is placed in air floating table x-axis obtains;
for tested parts rotating part is placed in the air floating table stage body angular acceleration that when air floating table bears x-axis, measurement obtains;
α is static-unbalance phasing degree;
G is test ground local gravitational acceleration.
The described method obtaining the disturbance torque that calibrated test component unbalance dynamic produces is:
T x 10 = a 10 cos ( ω 0 t + α 10 ) = J x ω . x 10
T x 11 = a 11 cos ( ω 0 t + α 11 ) = J x ω . x 11
T xp11=a xp11cos(ω 0t+α xp11)=T x11-T x10
Δα=α xpxp11
T xn10=a 10cos(ω 0t+α 10+Δα)
Wherein:
T x10for X-axis disturbance torque when tested parts rotating part rotates before increase counterweight;
A 10for X-axis disturbance torque amplitude when tested parts rotating part rotates before increase counterweight;
ω 0for the velocity of rotation of tested parts rotating part;
α 10for X-axis disturbance torque phase place when tested parts rotating part rotates before increase counterweight;
J xfor the moment of inertia of air floating table X-axis;
for the derivative of X-axis angular velocity Sine-Fitting result when tested parts rotating part rotates before increase counterweight;
T x11for X-axis disturbance torque when tested parts rotating part rotates after increase counterweight;
A 11for X-axis disturbance torque amplitude when tested parts rotating part rotates after increase counterweight;
α 11for X-axis disturbance torque phase place when tested parts rotating part rotates after increase counterweight;
for the derivative of X-axis angular velocity Sine-Fitting result when tested parts rotating part rotates after increase counterweight;
T xp11for the X-axis disturbance torque that counterweight when tested parts rotating part rotates produces;
A xp11for the X-axis disturbance torque amplitude that counterweight when tested parts rotating part rotates produces;
α xp11for the X-axis disturbance torque phase place that counterweight when tested parts rotating part rotates produces;
α xpfor counterweight produces the notional phase of disturbance torque;
Δ α is phase differential;
T xn10for the disturbance torque that the test component unbalance dynamic after phase calibration produces.
The method of described calculating tested parts rotating part couple-unbalance is:
T d = T xn 10 - mr ( g + L z ω 0 2 ) cos ( ω 0 t + α - π 2 )
C = | | T d | | ω 0 2
Wherein:
T dfor the couple that couple-unbalance produces;
T xn10for the disturbance torque that the test component unbalance dynamic after phase calibration produces;
Mr is tested parts rotating part static-unbalance;
α is static unbalance phasing degree;
ω 0for the velocity of rotation of tested parts rotating part;
L zfor tested parts rotating part barycenter is to the vertical range of air floating table center of rotation;
G is test ground local gravitational acceleration;
C is the couple-unbalance of tested parts rotating part.
The present invention's beneficial effect is compared with prior art: the present invention is by the test to spacecraft rotating part sound amount of unbalance, obtain quiet unbalancing value during spacecraft rotatable parts operation on orbit, on the basis of test result by counterweight to reduce the size of unbalancing value, during making spacecraft rotatable parts operation on orbit, the interference of celestial body is met and allow requirement, the present invention has directly applied to many satellites with flexible rotatable parts, shows that the method successfully solves the quiet unbalance dynamic test problem of flexible rotatable parts by data in-orbit.The method is simple, and measuring accuracy is high.
Accompanying drawing explanation
Fig. 1 is realization flow figure of the present invention;
Fig. 2 is the composition structural drawing of testing apparatus of the present invention;
Fig. 3 is the angular velocity curve of air floating table X-axis;
Fig. 4 disturbs the air floating table X-axis angular velocity change curve caused;
Air floating table X-axis attitude angle curve when Fig. 5 is static unbalance test;
Air floating table X-axis attitude angular velocity curve when Fig. 6 is static unbalance test;
Air floating table X-axis attitude angle curve when Fig. 7 is unbalance dynamic test
Air floating table X-axis attitude angular velocity curve when Fig. 8 is unbalance dynamic test;
Fig. 9 is air floating table X-axis attitude angular velocity FFT transformation curve;
Figure 10 is filtered air floating table X-axis attitude angular velocity curve;
Figure 11 is the matched curve of air floating table X-axis angular velocity;
Air floating table X-axis attitude angle curve when Figure 12 is unbalance dynamic test after counterweight;
Air floating table X-axis attitude angular velocity curve when Figure 13 is unbalance dynamic test after counterweight;
Figure 14 is air floating table X-axis angular velocity FFT transformation curve after counterweight;
Figure 15 is the air floating table X-axis attitude angular velocity curve after counterweight post filtering;
Figure 16 is air floating table X-axis angular velocity matched curve after counterweight.
Embodiment
Below in conjunction with the drawings and specific embodiments, further detailed description is done to the present invention:
The unbalance dynamic of rotatable parts determines primarily of static-unbalance and couple-unbalance, and couple-unbalance directly cannot be tested and obtains, and needs to be obtained by static unbalance and unbalance dynamic test result calculations.After obtaining static-unbalance and couple-unbalance, can eliminate by quality counterweight the disturbance torque that quiet unbalance dynamic brings greatly, thus reduce rotatable parts and rotate impact on the attitude of satellite.The attitude angular velocity that this test utilizes three floating gyros to carry out three-axis air-bearing table carries out high-acruracy survey, and then the angular acceleration of stage body is obtained by curve, also need to know the moment of inertia of air floating table to obtain disturbance torque, therefore main testing procedure, as shown in Figure 1:
(1) step of air floating table X, Y-axis moment of inertia is measured: float air floating table, to air floating table coarse balance, the counterweight of known quality is placed in a certain position on air floating table, the gravitational torque utilizing counterweight to produce makes air floating table produce angular acceleration, adopt the stage body angular velocity of gyro to measure air floating table, the stage body angular velocity of gyro to measure is carried out curve fitting and obtains the stage body angular acceleration of air floating table, the angular acceleration that the interference of deduction air floating table self produces, determines the moment of inertia of air floating table X, Y-axis according to theorem of angular momentum;
(2) step of tested parts rotating part static-unbalance is measured: tested parts rotating part is placed in respectively the x-axis of air floating table and negative x-axis direction, float air floating table, to air floating table coarse balance tested parts are not rotated, the gravitational torque utilizing tested parts to produce makes air floating table produce angular acceleration, adopt the stage body angular velocity of gyro to measure air floating table, the stage body angular velocity of gyro to measure is carried out curve fitting and obtains the stage body angular acceleration of air floating table, according to the air floating table X that theorem of angular momentum integrating step (1) is determined, the moment of inertia of Y-axis calculates the static-unbalance of tested parts rotating part,
(3) step of the disturbance torque that tested parts unbalance dynamic produces is measured: be placed on air floating table by tested parts rotating part, float air floating table, to air floating table coarse balance, the rotating part of tested parts is rotated, under rotation condition, adopt stage body attitude angle and the angular velocity of gyro to measure air floating table, the air floating table stage body attitude angle utilizing wave filter to go out gyro to measure and angular velocity carry out filtering and spectrum analysis and obtain the disturbance torque that tested parts unbalance dynamic produces;
(4) counterweight of certain mass is installed on tested parts rotating part, the disturbance torque of the tested parts unbalance dynamic generation after installing counterweight is obtained according to the measuring process of step (3), the disturbance torque relatively installing tested parts unbalance dynamic generation before and after counterweight obtains the disturbance torque phase place of counterweight generation, calibrate the tested parts rotation phase that disturbance torque data initial time is corresponding, thus obtain the disturbance torque of calibrated test component unbalance dynamic generation;
(5) disturbance torque that the test component unbalance dynamic that the tested parts static-unbalance utilizing step (2) to measure to obtain, step (4) calibrate produces calculates the couple-unbalance of tested parts rotating part in conjunction with the installation site of tested parts on air floating table.
Test system as shown in Figure 2, is made up of simulation control subsystem and uphole equipment on three-axis air-bearing table, platform.During test, the rotating part of tested parts is positioned at the top of three-axis air-bearing table, and the rotation axis of rotating part is vertical with horizontal direction, overlaps with the Z axis of air floating table.The angular velocity of air floating table is obtained by the gyro to measure be positioned on air floating table, and test data and test instruction realize transmission by the wireless communication module on platform between industrial computer and ground control cabinet.
The method of the moment of inertia of measurement air floating table X, Y-axis in step (1): the moment utilizing counterweight to produce and the angular acceleration of stage body, utilize theorem of angular momentum to calculate air floating table moment of inertia, computing formula is as follows:
J = mgl ω · 1 - ω · 2 - - - ( 1 )
Wherein J is the moment of inertia of air floating table X-axis or Y-axis, and m is counterbalance mass, and g is test ground local gravitational acceleration, and l is the horizontal range of counterweight center and air floating table center of rotation, the air floating table angular acceleration obtained is measured after increasing counterweight, for the air floating table angular acceleration that measurement when not increasing counterweight obtains, i.e. the angular acceleration that produces of air floating table self interference.
Fig. 3 is the angular velocity curve of air floating table X-axis under counterweight effect, and Fig. 4 disturbs the air floating table X-axis angular velocity change curve caused.After repetitive measurement, calculate air floating table X by formula (1), Y-axis moment of inertia is as shown in table 1.
Table 1 moment of inertia reckoner
The method calculating tested parts rotating part static-unbalance is:
mr x = J x ( ω · x 0 - ω · x 180 ) 2 g
mr y = J y ( ω · y 0 - ω · y 180 ) 2 g
mr = mr x 2 + mr y 2
α = ctg ( mr y mr x )
Wherein:
Mr xfor the X-axis component of static-unbalance under tested part coordinates system;
Mr yfor the Y-axis component of static-unbalance under tested part coordinates system;
Mr is tested parts rotating part static-unbalance;
J xfor the moment of inertia of air floating table X-axis;
J yfor the moment of inertia of air floating table Y-axis;
for the air floating table stage body angular acceleration that measurement when tested parts rotating part is placed in air floating table x-axis obtains;
for tested parts rotating part is placed in the air floating table stage body angular acceleration that when air floating table bears x-axis, measurement obtains;
α is static-unbalance phasing degree;
G is test ground local gravitational acceleration.
The attitude angle curve of air floating table X-axis when Fig. 5 is static unbalance test, the angular velocity change curve of air floating table X-axis when Fig. 6 is static unbalance test, calculate static-unbalance when tested parts rotating part benchmark is placed in zero-bit (being 0 degree), result of calculation is as shown in table 2.
Table 2 static-unbalance test chart
Static-unbalance mr(kgmm) mr x(kgmm) mr y(kgmm) Phase angle [alpha] (°)
58.4 -25.3 52.6 115.68
Unbalance dynamic is tested: during test, and tested parts rotate, and by air floating table coarse balance under rotation condition, make the extended testing system time as much as possible, improves measuring accuracy; The attitude angle using gyro to measure three-axis air-bearing table to export and angular velocity, extract the frequency band signals needed by wave filter, analyze the disturbance torque providing tested parts unbalance dynamic and produce.
The cycle indicator signal (rotating part zero cross signal) provided with tested parts is starting point, test curve when Fig. 7-Figure 10 is the rotation of tested parts, wherein Fig. 7 is X-axis attitude angle curve, Fig. 8 is X-axis attitude angular velocity curve, Fig. 9 is X-axis angular velocity signal FFT transformation curve, and Figure 10 is the X-axis angular velocity curve of stable state after filtering.Filtered X-axis attitude angular velocity signal is carried out Sine-Fitting, as shown in figure 11, and calculates corresponding disturbance torque according to air floating table moment of inertia.
X-axis angular velocity Sine-Fitting result: ω x10=0.01239sin (3.697t-2.132) (°/s) (6)
X-axis disturbance torque is: T x 10 = J x ω · x 10 = 1.9805 cos ( 3.697 t - 2.132 ) ( Nm ) - - - ( 7 )
Unbalance dynamic test start-phase is demarcated: the counterweight of installing certain mass on tested parts rotating part, by comparing the disturbance torque change of installing before and after counterweight, obtain the disturbance torque phase place that counterweight produces, because counterweight is known in the position of tested parts rotating part, therefore the rotation phase of tested parts corresponding to disturbance torque data initial time can be marked.
The method obtaining the disturbance torque that calibrated test component unbalance dynamic produces is:
T x 10 = a 10 cos ( ω 0 t + α 10 ) = J x ω · x 10
T x 11 = a 11 cos ( ω 0 t + α 11 ) = J x ω · x 11
T xp11=a xp11cos(ω 0t+α xp11)=T x11-T x10
Δα=α xpxp11
T xn10=a 10cos(ω 0t+α 10+Δα)
Wherein:
T x10for X-axis disturbance torque when tested parts rotating part rotates before increase counterweight;
A 10for X-axis disturbance torque amplitude when tested parts rotating part rotates before increase counterweight;
ω 0for the velocity of rotation of tested parts rotating part;
α 10for X-axis disturbance torque phase place when tested parts rotating part rotates before increase counterweight;
J xfor the moment of inertia of air floating table X-axis;
for the derivative of X-axis angular velocity Sine-Fitting result when tested parts rotating part rotates before increase counterweight;
T x11for X-axis disturbance torque when tested parts rotating part rotates after increase counterweight;
A 11for X-axis disturbance torque amplitude when tested parts rotating part rotates after increase counterweight;
α 11for X-axis disturbance torque phase place when tested parts rotating part rotates after increase counterweight;
for the derivative of X-axis angular velocity Sine-Fitting result when tested parts rotating part rotates after increase counterweight;
T xp11for the X-axis disturbance torque that counterweight when tested parts rotating part rotates produces;
A xp11for the X-axis disturbance torque amplitude that counterweight when tested parts rotating part rotates produces;
α xp11for the X-axis disturbance torque phase place that counterweight when tested parts rotating part rotates produces;
α xpfor counterweight produces the notional phase of disturbance torque;
Δ α is phase differential;
T xn10for the disturbance torque that the test component unbalance dynamic after phase calibration produces.
After increasing counterweight, tested parts rotate, and the cycle indicator signal provided with tested parts is starting point, and filtered X-axis attitude angular velocity signal is carried out Sine-Fitting, and as shown in figure 16, then, after increasing counterweight, X-axis angular velocity Sine-Fitting result is:
ω x11=0.006708sin(3.697t-2.867)(°/s)(8)
Calculate corresponding disturbance torque according to air floating table moment of inertia, then, after increasing counterweight, when tested parts rotate, X-axis disturbance torque is:
T x 11 = J x ω · x 11 = 1.0723 cos ( 3.697 t - 2.867 ) ( Nm ) - - - ( 9 )
Disturbance torque according to can obtain counterweight generation without disturbance torque during counterweight and the disturbance torque after increasing counterweight is: T xp11=T x11-T x10=1.386cos (3.696t+1.555) (Nm) (10)
The benchmark defining tested parts rotating part moment that overlaps with tested parts X-axis (i.e. zero-bit) is nominal initial time.From the installation of counterweight, the phase place that counterweight relatively rotates reference part on tested parts is 326.25 degree, when tested parts rotate (during at nominal initial time) when reference part overlaps with tested parts X-axis, the gravitational torque that counterweight produces and the phase place of centrifugal moment under tested part coordinates system are 236.25 degree.Calculate known according to formula (10): when the cycle indicator signal chosen tested parts and provide is starting point, the moment that counterweight produces is 89.095 degree relative to the phase place of tested part coordinates system, then the now phase place of data more advanced than nominal initial time 147.155 degree (2.568 radians).Be then starting point with the nominal time, can be obtained by according to formula (6), (7), test component unbalance dynamic test X-axis angular velocity and corresponding disturbance torque are:
ω xn10=0.01239sin(3.697t+0.436)(°/s)(11)
T xn10=1.9805cos(3.697t+0.436)(Nm)(12)
Couple-unbalance calculates: when tested parts rotate, couple being expressed as (taking the nominal time as starting point) under test component coordinate system that couple-unbalance produces:
T d = T xn 10 - mr ( g + L z ω 1 2 ) cos ( ω 1 t + α - π 2 ) = 1.9805 cos ( 3.696 t + 0.436 ) - 0.0584 * ( 9.8015 + 1.93 * 3.696 2 ) cos ( 3.696 t + 0.4482 ) = 0.1339 cos ( 3.696 t - 2.512 ) ( Nm ) - - - ( 13 )
Wherein T xn10for the disturbance torque that the test component unbalance dynamic after phase calibration produces, mr is static-unbalance, and α is static unbalance phasing degree, ω 1for the velocity of rotation of rotating part, L zthe rotating part barycenter that measurement obtains is to the vertical range of air floating table center of rotation, and g is test ground local gravitational acceleration.
Then tested parts rotating part couple-unbalance size is:
C = | | T d | | ω 1 2 = 0.1339 ( 2 π / 1.7 ) 2 = 0.009802 kgm 2 - - - ( 14 )
Couple initial phase is under tested part coordinates system :-2.512rad
Test result gathers as shown in table 4.
Table 4 test result gathers
The content be not described in detail in instructions belongs to the known technology of professional and technical personnel in the field.

Claims (4)

1. test a method for rotatable parts sound amount of unbalance, it is characterized in that step is as follows:
(1) step of air floating table X, Y-axis moment of inertia is measured: float air floating table, to air floating table coarse balance, the counterweight of known quality is placed in a certain position on air floating table, the gravitational torque utilizing counterweight to produce makes air floating table produce angular acceleration, adopt the stage body angular velocity of gyro to measure air floating table, the stage body angular velocity of gyro to measure is carried out curve fitting and obtains the stage body angular acceleration of air floating table, the angular acceleration that the interference of deduction air floating table self produces, determines the moment of inertia of air floating table X, Y-axis according to theorem of angular momentum;
(2) step of tested parts rotating part static-unbalance is measured: tested parts rotating part is placed in respectively the x-axis of air floating table and negative x-axis direction, float air floating table, to air floating table coarse balance tested parts are not rotated, the gravitational torque utilizing tested parts to produce makes air floating table produce angular acceleration, adopt the stage body angular velocity of gyro to measure air floating table, the stage body angular velocity of gyro to measure is carried out curve fitting and obtains the stage body angular acceleration of air floating table, according to the air floating table X that theorem of angular momentum integrating step (1) is determined, the moment of inertia of Y-axis calculates the static-unbalance of tested parts rotating part,
Wherein, the method calculating tested parts rotating part static-unbalance is:
mr x = J x ( ω · x 0 - ω · x 180 ) 2 g
mr y = J y ( ω · y 0 - ω · y 180 ) 2 g
mr = mr x 2 + mr y 2
α = ctg ( mr y mr x )
Wherein:
Mr xfor the X-axis component of static-unbalance under tested part coordinates system;
Mr yfor the Y-axis component of static-unbalance under tested part coordinates system;
Mr is tested parts rotating part static-unbalance;
J xfor the moment of inertia of air floating table X-axis;
J yfor the moment of inertia of air floating table Y-axis;
for the air floating table stage body angular acceleration that measurement when tested parts rotating part is placed in air floating table x-axis obtains;
for tested parts rotating part is placed in the air floating table stage body angular acceleration that when air floating table bears x-axis, measurement obtains;
α is static-unbalance phasing degree;
G is test ground local gravitational acceleration;
(3) step of the disturbance torque that tested parts unbalance dynamic produces is measured: be placed on air floating table by tested parts rotating part, float air floating table, to air floating table coarse balance, the rotating part of tested parts is rotated, under rotation condition, adopt stage body attitude angle and the angular velocity of gyro to measure air floating table, the air floating table stage body attitude angle utilizing wave filter to go out gyro to measure and angular velocity carry out filtering and spectrum analysis and obtain the disturbance torque that tested parts unbalance dynamic produces;
(4) counterweight of certain mass is installed on tested parts rotating part, the disturbance torque of the tested parts unbalance dynamic generation after installing counterweight is obtained according to the measuring process of step (3), the disturbance torque relatively installing tested parts unbalance dynamic generation before and after counterweight obtains the disturbance torque phase place of counterweight generation, calibrate the tested parts rotation phase that disturbance torque data initial time is corresponding, thus obtain the disturbance torque of calibrated test component unbalance dynamic generation;
(5) disturbance torque that the test component unbalance dynamic that the tested parts static-unbalance utilizing step (2) to measure to obtain, step (4) calibrate produces calculates the couple-unbalance of tested parts rotating part in conjunction with the installation site of tested parts on air floating table.
2. a kind of method of testing rotatable parts sound amount of unbalance according to claim 1, is characterized in that: the described method utilizing theorem of angular momentum to determine the moment of inertia of air floating table X, Y-axis:
J = mgl ω · 1 - ω · 2
Wherein J is the moment of inertia of air floating table X-axis or Y-axis, and m is counterbalance mass, and g is test ground local gravitational acceleration, and l is the horizontal range of counterweight center and air floating table center of rotation, the air floating table angular acceleration obtained is measured after increasing counterweight, for the air floating table angular acceleration that measurement when not increasing counterweight obtains, i.e. the angular acceleration that produces of air floating table self interference.
3. a kind of method of testing rotatable parts sound amount of unbalance according to claim 1, is characterized in that: described in obtain the method for disturbance torque that calibrated test component unbalance dynamic produces and be:
T x 10 = a 10 cos ( ω 0 t + α 10 ) = J x ω · x 10
T x 11 = a 11 cos ( ω 0 t + α 11 ) = J x ω · x 11
T xp11=a xp11cos(ω 0t+α xp11)=T x11-T x10
Δα=α xpxp11
T xn10=a 10cos(ω 0t+α 10+Δα)
Wherein:
T x10for X-axis disturbance torque when tested parts rotating part rotates before increase counterweight;
A 10for X-axis disturbance torque amplitude when tested parts rotating part rotates before increase counterweight;
ω 0for the velocity of rotation of tested parts rotating part;
α 10for X-axis disturbance torque phase place when tested parts rotating part rotates before increase counterweight;
J xfor the moment of inertia of air floating table X-axis;
for the derivative of X-axis angular velocity Sine-Fitting result when tested parts rotating part rotates before increase counterweight;
T x11for X-axis disturbance torque when tested parts rotating part rotates after increase counterweight;
A 11for X-axis disturbance torque amplitude when tested parts rotating part rotates after increase counterweight;
α 11for X-axis disturbance torque phase place when tested parts rotating part rotates after increase counterweight;
for the derivative of X-axis angular velocity Sine-Fitting result when tested parts rotating part rotates after increase counterweight;
T xp11for the X-axis disturbance torque that counterweight when tested parts rotating part rotates produces;
A xp11for the X-axis disturbance torque amplitude that counterweight when tested parts rotating part rotates produces;
α xp11for the X-axis disturbance torque phase place that counterweight when tested parts rotating part rotates produces;
α xpfor counterweight produces the notional phase of disturbance torque;
Δ α is phase differential;
T xn10for the disturbance torque that the test component unbalance dynamic after phase calibration produces.
4. a kind of method of testing rotatable parts sound amount of unbalance according to claim 1, is characterized in that: the method for described calculating tested parts rotating part couple-unbalance is:
T d = T xn 10 - mr ( g + L z ω 0 2 ) cos ( ω 0 t + α - π 2 )
C = | | T d | | ω 0 2
Wherein:
T dfor the couple that couple-unbalance produces;
T xn10for the disturbance torque that the test component unbalance dynamic after phase calibration produces;
Mr is tested parts rotating part static-unbalance;
α is static unbalance phasing degree;
ω 0for the velocity of rotation of tested parts rotating part;
L zfor tested parts rotating part barycenter is to the vertical range of air floating table center of rotation;
G is test ground local gravitational acceleration;
C is the couple-unbalance of tested parts rotating part.
CN201310319876.8A 2013-07-26 2013-07-26 A kind of method of testing rotatable parts sound amount of unbalance Active CN103424225B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310319876.8A CN103424225B (en) 2013-07-26 2013-07-26 A kind of method of testing rotatable parts sound amount of unbalance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310319876.8A CN103424225B (en) 2013-07-26 2013-07-26 A kind of method of testing rotatable parts sound amount of unbalance

Publications (2)

Publication Number Publication Date
CN103424225A CN103424225A (en) 2013-12-04
CN103424225B true CN103424225B (en) 2015-11-25

Family

ID=49649308

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310319876.8A Active CN103424225B (en) 2013-07-26 2013-07-26 A kind of method of testing rotatable parts sound amount of unbalance

Country Status (1)

Country Link
CN (1) CN103424225B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109632185A (en) * 2018-12-13 2019-04-16 上海航天控制技术研究所 A kind of in-orbit estimation method of large inertia rotatable parts unbalance dynamic characteristic

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104713680A (en) * 2013-12-11 2015-06-17 中国航空工业第六一八研究所 Method for static imbalance test for inertial platform framework
CN104763632B (en) * 2015-04-08 2017-03-29 上海磁浮交通发展有限公司 A kind of oil-free vortex air compressor balance method
CN104748913A (en) * 2015-04-09 2015-07-01 哈尔滨工业大学 Double freedom degree precision centrifuge spindle dynamic unbalance scale identification method based on micro displacement sensor
CN107631738B (en) * 2017-09-20 2021-02-09 哈尔滨工业大学 Method for identifying imbalance of gyro flywheel couple
CN108519181B (en) * 2018-03-09 2020-06-09 北京航天控制仪器研究所 Modeling and testing method for mass unbalance moment of platform body
CN108760118B (en) * 2018-03-09 2020-07-14 北京航天控制仪器研究所 Device and method for measuring mass unbalance moment of platform body of inertia platform
CN109632186A (en) * 2018-12-21 2019-04-16 上海航天控制技术研究所 The in-orbit estimation method of unbalance dynamic characteristic and equipment of large inertia rotatable parts
CN110285922B (en) * 2019-01-31 2020-11-13 上海卫星工程研究所 Static unbalance test modeling method for two-dimensional rotary table
CN110196139B (en) * 2019-01-31 2020-12-11 上海卫星工程研究所 Static unbalance testing method for two-dimensional rotary table
CN110702315A (en) * 2019-06-06 2020-01-17 上海航天控制技术研究所 On-orbit estimation method and device for static imbalance characteristic of large-inertia rotating part
CN111982402B (en) * 2020-04-09 2021-11-05 西北工业大学 Aircraft engine rotor part matching optimization method considering initial unbalance amount

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102620892A (en) * 2011-12-15 2012-08-01 上海卫星工程研究所 Dynamic balance testing method for rotatable part
CN102778333A (en) * 2011-12-15 2012-11-14 上海卫星工程研究所 Method for performing dynamic balance test on large rotating part

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102620892A (en) * 2011-12-15 2012-08-01 上海卫星工程研究所 Dynamic balance testing method for rotatable part
CN102778333A (en) * 2011-12-15 2012-11-14 上海卫星工程研究所 Method for performing dynamic balance test on large rotating part

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Air Bearing Based Satellite Attitude Dynamics Simulator for Control Software Research and Development;B. N. Agrawal ert al.;《Proceedings of SPIE》;20011231;第4366卷;第204-214页 *
一种利用单轴气浮平台测量角动量的方法;许峰 等;《科学技术与工程》;20071231;第6479-6482页 *
测试单轴气浮台转动惯量的一种方法;李季苏;《航天控制》;19941231;第69-73页 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109632185A (en) * 2018-12-13 2019-04-16 上海航天控制技术研究所 A kind of in-orbit estimation method of large inertia rotatable parts unbalance dynamic characteristic

Also Published As

Publication number Publication date
CN103424225A (en) 2013-12-04

Similar Documents

Publication Publication Date Title
CN103424225B (en) A kind of method of testing rotatable parts sound amount of unbalance
CN102620892B (en) Dynamic balance testing method for rotatable part
CN103575299B (en) Utilize dual-axis rotation inertial navigation system alignment and the error correcting method of External Observation information
CN104296908B (en) Three freedom degree air floating platform disturbance torque composition measuring apparatus
CN101290326B (en) Parameter identification calibration method for rock quartz flexibility accelerometer measuring component
CN101975872B (en) Method for calibrating zero offset of quartz flexible accelerometer component
CN102879014B (en) Optical imaging autonomous navigation semi-physical simulation testing system for deep space exploration proximity process
CN101105503B (en) Acceleration meter assembling error scalar rectification method for strapdown type inertia navigation measurement combination
CN103743413B (en) Heeling condition modulated is sought northern instrument alignment error On-line Estimation and is sought northern error compensating method
CN103090866B (en) Method for restraining speed errors of single-shaft rotation optical fiber gyro strapdown inertial navigation system
CN102735267B (en) Measuring method for inertial measurement device in sled testing
CN102636183B (en) Quadratic overload term test method for flexible gyroscope based on optical fiber monitoring and centrifuge with two-axis turntable
CN106844880A (en) A kind of satellite momentum wheel disturbs shake experiment and data interpretation method
CN102072797B (en) Method for measuring unbalance of spindle in measurement of dynamic balance of tire and tire balancing machine
CN103234560B (en) A kind of method of strapdown inertial navitation system (SINS) Zero positioning
CN105300597B (en) Three-axis air-bearing table barycenter balancing method and device
CN103323625B (en) Error calibration compensation method of accelerometers in MEMS-IMU under dynamic environment
CN102768043B (en) Integrated attitude determination method without external observed quantity for modulated strapdown system
CN106246163A (en) Nearly drill bit dynamic directional survey method and device
CN105807780B (en) A kind of anti-interference attitude control method and checking device based on flywheel output bias
CN105628025A (en) Constant-rate offset frequency/mechanically dithered laser gyro inertial navigation system navigation method
CN102778333B (en) Method for performing dynamic balance test on large rotating part
CN106017452A (en) Dual gyro anti-disturbance north-seeking method
CN101738203A (en) Optimal position calibration method of static drifting zero and primary acceleration related term error model of flexible gyroscope
CN103591919B (en) For the method and apparatus that precision centrifuge static radius is measured

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant