CN105628176A - Rotation machinery torsional vibration signal collection analysis method - Google Patents

Rotation machinery torsional vibration signal collection analysis method Download PDF

Info

Publication number
CN105628176A
CN105628176A CN201610011172.8A CN201610011172A CN105628176A CN 105628176 A CN105628176 A CN 105628176A CN 201610011172 A CN201610011172 A CN 201610011172A CN 105628176 A CN105628176 A CN 105628176A
Authority
CN
China
Prior art keywords
sampling
torsional angle
instantaneous
signal
tor
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.)
Pending
Application number
CN201610011172.8A
Other languages
Chinese (zh)
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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Henan Electric Power Co Ltd
Henan Enpai High Tech Group Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Henan Electric Power Co Ltd
Henan Enpai High Tech Group Co Ltd
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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Henan Electric Power Co Ltd, Henan Enpai High Tech Group Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201610011172.8A priority Critical patent/CN105628176A/en
Publication of CN105628176A publication Critical patent/CN105628176A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H11/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties

Abstract

The invention discloses a rotation machinery torsional vibration signal collection analysis method, comprising steps of installing a marking device at a rotation shaft measurement cross section of a device to be tested, installing a matching detection device, removing low frequency interference in an outputted pulse signal, converting the pulse signal to a standard TTL pulse signal, utilizing a high frequency counter to perform signal collection to obtain a series counting value, using an obtained series counting value to solve a series transient angular speed value, solving the series uniform angle collection transient torsional angle speed value, solving a series equal-time-interval sampling transient torsional angle speed value, solving the real torsional angle speed signal, performing integral transformation on a real torsional angle speed signal after the trend item is eliminated to obtain a torsional vibration angular displacement signal. The invention maximally improves the accuracy and reliability of rotation machinery torsional vibration signal analysis and has high reliability and high measurement accuracy.

Description

A kind of rotating machinery torsion vibration signal capturing analysis method
Technical field
The present invention relates to a kind of rotating machinery torsion vibration signal capturing analysis method, particularly relate to a kind of rotating machinery torsion vibration signal capturing analysis method suitable in all kinds of rotating machineries such as steam turbine, gas turbine, electromotor, compressor, blower fan, pump and electromotors.
Background technology
Twisting vibration is the recurrent class vibration faults of rotating machinery such as steam turbine, electromotor, blower fan. Under significantly twisting vibration impact or long-term twisting vibration effect, rotating shaft surface there will be " Flos Cannabis formula " crackle with cross section angle at 45 ��, it is easy to causes rotating shaft fatigue damage. With shaft bending vibration-testing the difference is that, the test of twisting vibration is increasingly complex, and the influence factor suffered by test result also compares many, and the measuring technology bending vibration measuring technology that is far from is so ripe. And twisting vibration is safe and stable to rotating machinery and reliability service has tremendous influence, therefore, twisting vibration Testing Technology Study is carried out, raising torsional oscillation test signal accuracy and reliability is significant.
Twisting vibration test is divided into contact and contactless two big classes. Contact measurement method includes: Geiger that method, dynamic strain method and acceleration transducer method etc. Contactless measurement includes: laser method and impulse method etc. For contact measurement method, the impulse method principle in non-contact measurement method is simple, easily realizes at the scene, is a kind of measuring method more common at present. The present invention is built upon on contactless impulsive measurement method basis, to alleviate the impact on impulse method torsional vibration measurement result of every interference factor for target, and proposed a kind of torsion vibration signal analytic process.
Summary of the invention
It is an object of the invention to provide a kind of rotating machinery torsion vibration signal capturing analysis method, contactless acquisition geometry is adopted to extract torsional vibration signals from the pulse signal detected, accuracy and reliability, the advantage with high reliability and high measurement accuracy that rotating machinery torsion vibration signal is analyzed can be greatly enhanced.
The present invention adopts following technical proposals:
A kind of rotating machinery torsion vibration signal capturing analysis method, comprises the following steps successively:
A: section is measured in the rotating shaft that the indication device being evenly arranged with sampled point is arranged on Devices to test, and installs supporting detection equipment, adjusts the distance between detection equipment and rotating shaft and angle, enables detection equipment stably output pulse signal;
B: remove the low-frequency disturbance in the pulse signal that detection equipment exports;
C: the pulse signal after removing low-frequency disturbance is converted into the output of standard TTL pulse signal;
D: utilize high-frequency counter that standard TTL pulse signal is carried out signal sampling, obtain the time that on indication device, each sampled point turns over, i.e. series count value;
E: utilize the serial count value obtained, according to the number of divisions on indication device, try to achieve the instantaneous angular velocity of each sampled point on indication device respectively, the set of the instantaneous angular velocity of each sampled point on indication device, constitute series instantaneous angular velocity value;
F: deduct rotating shaft mean angular velocity value from indication device in the instantaneous angular velocity of each sampled point, try to achieve the instantaneous torsional angle velocity amplitude of the equiangular sampling of each sampled point on indication device respectively, the set of the instantaneous torsional angle velocity amplitude of the equiangular sampling of each sampled point on indication device, constitutes the series instantaneous torsional angle velocity amplitude of equiangular sampling;
G: the instantaneous torsional angle velocity amplitude of equiangular sampling of each sampled point is carried out resampling, instantaneous for equiangular sampling torsional angle velocity amplitude is converted into the constant duration instantaneous torsional angle velocity amplitude of sampling of each sampled point obtained with constant duration sample mode, the set of the constant duration instantaneous torsional angle velocity amplitude of sampling of each sampled point obtained, constitutes the series constant duration instantaneous torsional angle velocity amplitude of sampling;
H: utilize least square method to eliminate the trend term in the series constant duration instantaneous torsional angle velocity amplitude of sampling, try to achieve true torsional angle rate signal;
I: the true torsional angle rate signal after eliminating trend term is integrated conversion, obtains torsional oscillation angular displacement signal.
In described step B, high pass filter is utilized to remove the low-frequency disturbance in the pulse signal that detection equipment exports.
In described step C, the pulse signal after removal low-frequency disturbance is converted into the standard TTL pulse signal output of 0-5V.
In described step D, utilizing high-frequency counter that standard TTL pulse signal is carried out signal sampling, obtain the time that on indication device, each sampled point turns over, namely series count value, is designated as �� t respectively1,��t2,...,��tN, wherein, N represents total number of sample points, i.e. total number of teeth of toothed disc.
In described step E, following formula is utilized to calculate at tnThe instantaneous angular velocity �� of ith sample point on moment indication devicen(i),
, N represents total number of sample points, i.e. total number of teeth of toothed disc; I represents ith sample point on indication device, i.e. the subscript of the i-th tooth on toothed disc, �� tiFor the time that ith sample point turns over.
In described step F, following formula is utilized to calculate the instantaneous torsional angle velocity amplitude �� of equiangular sampling of ith sample point on indication devicetor(i),
Subscript tor represents torsional angle speed, and subscript n represents that instantaneous angular velocity, i represent ith sample point on indication device,For rotating shaft mean angular velocity value,Calculated by below equation:
In described step G, the method utilizing linear interpolation, the instantaneous torsional angle velocity amplitude �� of equiangular sampling to each sampled pointtorI () carries out resampling, be translated into the constant duration instantaneous torsional angle velocity amplitude �� of sampling of each sampled pointtor' (t), in variable, subscript �� represents constant duration sampling;
If sampling time interval is �� t, the sampling time corresponding to sample point m is tm, the time respectively t of 2 of the front and back moment adjacent with the sample point m equiangularly sampled point that sample mode obtainsk,tk+1, then the sample point m obtained is sampled at t=t with constant durationmThe constant duration instantaneous torsional angle velocity amplitude of sampling in momentFor:
ω t o r ′ ( t ) | t = t m = ω t o r ( t ) t = t k + t m - t k t k + 1 - t k [ ω t o r ( t ) | t = t k + 1 - ω t o r ( t ) | t = t k ] ,
Wherein,For sampling the sample point m obtained at t=t with constant durationmThe constant duration instantaneous torsional angle velocity amplitude of sampling in moment,WithThe sample point that respectively equiangularly sample mode the obtains instantaneous angular velocity value at k and k+1 sample point place.
In described step H, least square method is utilized to eliminate the constant duration instantaneous torsional angle velocity amplitude �� of samplingtorTrend term in ' (t);
The instantaneous torsional angle velocity amplitude �� if constant duration is sampledtorTrend term �� in ' (t)1T () is the quadratic polynomial about time variable t:
��1(t)=at2+ bt+c,
Coefficient a, b, the c in multinomial is solved, from instantaneous angular velocity value �� with method of least squaretor' (t) eliminates trend term ��1T (), tries to achieve true torsional angle rate signal �� "tor(t), �� "tor(t)=�� 'tor(t)-��1(t)��
In described step I, to true torsional angle rate signal �� "torT () is integrated, obtain torsional oscillation angular displacement signal �� (t), �� (t)=�� �� "tor(t)dt��
The described indication device being provided with sampled point adopts the uniform toothed disc of calibration, and supporting detection equipment adopts eddy current sensor or magnetoresistive transducer.
There is advantages that
(1) it is configured with signal conditioner, pulse signal conditioning detection equipment Inspection arrived is the output of 0-5V standard TTL pulse signal, high-frequency counter can be accessed easily, to improve counting precision and the accuracy of detection unit count device, reduce the phenomenons such as leakage tooth;
(2) signal that equiangularly sampling configuration obtains be converted into sample the signal that obtains with constant duration, decrease signal analysis error, it is to avoid when data post processes (during such as spectrum analysis) produce deviation effects test result;
(3) reject the trend term in signal with method of least square, decrease and obtained the distorted signals in angular displacement signal process by angular velocity signal integration.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the present invention;
Fig. 2 is the torsional vibration signals figure of detection equipment output;
The torsional vibration signals that Fig. 3 is the output of detection equipment modulates interval pulse signal graph heterogeneous;
Fig. 4 is equiangular sampling and constant duration sampling difference schematic diagram;
Fig. 5 is the torsional angle displacement oscillogram before rejecting trend term;
Fig. 6 is the torsional angle shift spectrum figure before rejecting trend term;
Fig. 5 is carried out rejecting the torsional angle displacement oscillogram after trend term processes by Fig. 7;
Fig. 6 is carried out rejecting the torsional angle shift spectrum figure after trend term processes by Fig. 8.
Detailed description of the invention
Below in conjunction with drawings and Examples, the present invention is made with detailed description:
As it is shown in figure 1, rotating machinery torsion vibration signal capturing analysis method of the present invention, comprise the following steps successively:
A: section is measured in the rotating shaft that the indication device being evenly arranged with sampled point is arranged on Devices to test, and installs supporting detection equipment, adjusts the distance between detection equipment and rotating shaft and angle, enables detection equipment stably output pulse signal.
In the present invention, the indication device being provided with sampled point adopts the uniform toothed disc of calibration, supporting detection equipment to adopt eddy current sensor or magnetoresistive transducer, eddy current sensor or magnetoresistive transducer that the pulse signal detected can be converted into signal of telecommunication output.
When the rotating shaft of Devices to test is without torsional oscillation, the pulse interval that eddy current sensor or magnetoresistive transducer collect is uniform; And when the rotating shaft of Devices to test exists torsional oscillation, the pulse interval that eddy current sensor or magnetoresistive transducer collect is uneven, wherein just contain torsional vibration signals, as shown in Figures 2 and 3.
B: remove the low-frequency disturbance in the pulse signal that detection equipment exports. In the present invention, available high pass filter removes the low-frequency disturbance in the pulse signal that detection equipment exports.
C: the pulse signal after removing low-frequency disturbance is converted into the output of standard TTL pulse signal. In the present invention, the pulse signal after removing low-frequency disturbance can be converted into the standard TTL pulse signal output of 0-5V.
D: utilize high-frequency counter that standard TTL pulse signal is carried out signal sampling, obtain the time that on indication device, each sampled point turns over, namely series count value, is designated as �� t respectively1,��t2,...,��tN, wherein, N represents total number of sample points, i.e. total number of teeth of toothed disc.
E: utilize the serial count value obtained, according to the number of divisions on indication device, try to achieve the instantaneous angular velocity of each sampled point on indication device respectively, i.e. series instantaneous angular velocity value;
The instantaneous angular velocity �� of each sampled point i on indication devicenI the set of (), constitutes series instantaneous angular velocity value.
In the present invention, following formula is utilized to calculate the instantaneous angular velocity �� of ith sample point on tn moment indication devicen(i),
, N represents total number of sample points, i.e. total number of teeth of toothed disc; I represents ith sample point on indication device, i.e. the subscript of the i-th tooth on toothed disc, �� tiFor the time that ith sample point turns over.
F: deduction rotating shaft mean angular velocity value from series instantaneous angular velocity value, obtains the series instantaneous torsional angle velocity amplitude of equiangular sampling;
The instantaneous torsional angle velocity amplitude �� of the equiangular sampling of each sampled point on indication devicetorI the set of (), constitutes the series instantaneous torsional angle velocity amplitude of equiangular sampling.
In the present invention, following formula is utilized to calculate the instantaneous torsional angle velocity amplitude �� of equiangular sampling of ith sample point on indication devicetor(i),
ω t o r ( i ) = ω n ( i ) - ω ‾ ;
Wherein, subscript tor represents torsional angle speed, and subscript n represents that instantaneous angular velocity, i represent ith sample point on indication device,For rotating shaft mean angular velocity value,Calculated by below equation:
ω ‾ = 1 N Σ i = 1 N ω n ( i ) ;
When sampling according to step A to F, each tooth on each sampled point and toothed disc on indication device, 1 pulse signal all can be being exported through eddy current sensor or magnetoresistive transducer, therefore obtained sampled signal equiangularly mode obtains, and the interval sampled unequal, this sample mode will cause producing great deviation (during such as spectrum analysis) when data post processes. And during follow-up signal analysis, if sampled signal obtains in constant duration mode, then can greatly reduce deviation when data post processes, it is ensured that certainty of measurement. Therefore, the present invention also has additional step G.
G: the instantaneous torsional angle velocity amplitude �� of equiangular sampling of each sampled point that equiangularly sample mode is obtainedtorI () carries out resampling, be translated into the constant duration instantaneous torsional angle velocity amplitude �� of sampling of each sampled point obtained with constant duration sample modetor' (t), in variable, subscript �� represents constant duration sampling.
The constant duration instantaneous torsional angle velocity amplitude �� of sampling of each sampled point on indication devicetorThe set of ' (t), constitutes the series constant duration instantaneous torsional angle velocity amplitude of sampling. Fig. 4 is equiangular sampling and constant duration sampling difference schematic diagram;
In the present invention, instantaneous for series equiangular sampling torsional angle velocity amplitude is converted into the series constant duration instantaneous torsional angle velocity amplitude of sampling by the method for available linear interpolation;
If sampling time interval is �� t, the sampling time corresponding to sample point m is tm, the time respectively t of 2 of the front and back moment adjacent with the sample point m equiangularly sampled point that sample mode obtainsk,tk+1, then the sample point m obtained is sampled at t=t with constant durationmThe constant duration instantaneous torsional angle velocity amplitude of sampling in momentFor:
ω t o r ′ ( t ) | t = t m = ω t o r ( t ) t = t k + t m - t k t k + 1 - t k [ ω t o r ( t ) | t = t k + 1 - ω t o r ( t ) | t = t k ] ,
Wherein,For sampling the sample point m obtained at t=t with constant durationmThe constant duration instantaneous torsional angle velocity amplitude of sampling in moment,WithThe sample point that respectively equiangularly sample mode the obtains instantaneous angular velocity value at k and k+1 sample point place.
H: in order to obtain instantaneous torsional angle displacement signal, it is necessary to sample instantaneous torsional angle velocity amplitude �� to constant durationtor' (t) is integrated. Due to the constant duration instantaneous torsional angle velocity amplitude �� of samplingtor' (t) exists trend term, i.e. the instantaneous torsional angle velocity amplitude �� of constant duration sampling in a period of timetorThe meansigma methods of ' (t) is not 0 so that sampled instantaneous torsional angle velocity amplitude �� by constant durationtorTo there is significantly drift phenomenon in the instantaneous torsional angle displacement signal that ' (t) integration obtains, cause distorted signals, and as shown in Figure 5 and Figure 6, Fig. 5 is the torsional angle displacement oscillogram before rejecting trend term, and Fig. 6 is the torsional angle shift spectrum figure before rejecting trend term.
Therefore, in the present invention, utilize least square method to eliminate the constant duration instantaneous torsional angle velocity amplitude �� of samplingtorTrend term in ' (t); The instantaneous torsional angle velocity amplitude �� if constant duration is sampledtorTrend term �� in ' (t)1T () is the quadratic polynomial about time variable t:
��1(t)=at2+ bt+c, solves coefficient a, b, the c in multinomial with method of least square, from the constant duration instantaneous torsional angle velocity amplitude �� of samplingtor' (t) eliminates trend term ��1T (), tries to achieve true torsional angle rate signal �� "tor(t), �� "tor(t)=�� 'tor(t)-��1(t)��
Fig. 5 is carried out rejecting the torsional angle displacement oscillogram after trend term processes by Fig. 7, and Fig. 6 is carried out rejecting the torsional angle shift spectrum figure after trend term processes by Fig. 8.
I: to eliminating the true torsional angle rate signal �� after trend term "torT () is integrated, obtain torsional oscillation angular displacement signal.
In step I, to true torsional angle rate signal �� "torT () is integrated, obtain torsional oscillation angular displacement signal �� (t),
�� (t)=�� �� "tor(t)dt��
The present invention adopts contactless acquisition geometry to extract torsional vibration signals from the pulse signal detected; accuracy and reliability that rotating machinery torsion vibration signal is analyzed can be greatly enhanced; present invention could apply to the monitoring of all kinds of rotating machinery twisting vibration, analyzing and diagnosing and protection, improve safety and reliability during unit operation
The advantage with high reliability and high measurement accuracy.
Compared with prior art, there is advantages that
(1) it is configured with signal conditioner, pulse signal conditioning detection equipment Inspection arrived is the output of 0-5V standard TTL pulse signal, high-frequency counter can be accessed easily, to improve counting precision and the accuracy of detection unit count device, reduce the phenomenons such as leakage tooth;
(2) signal that equiangularly sampling configuration obtains be converted into sample the signal that obtains with constant duration, decrease signal analysis error, it is to avoid when data post processes (during such as spectrum analysis) produce deviation effects test result;
(3) reject the trend term in signal with method of least square, decrease and obtained the distorted signals in angular displacement signal process by angular velocity signal integration.

Claims (10)

1. a rotating machinery torsion vibration signal capturing analysis method, it is characterised in that comprise the following steps successively:
A: section is measured in the rotating shaft that the indication device being evenly arranged with sampled point is arranged on Devices to test, and installs supporting detection equipment, adjusts the distance between detection equipment and rotating shaft and angle, enables detection equipment stably output pulse signal;
B: remove the low-frequency disturbance in the pulse signal that detection equipment exports;
C: the pulse signal after removing low-frequency disturbance is converted into the output of standard TTL pulse signal;
D: utilize high-frequency counter that standard TTL pulse signal is carried out signal sampling, obtain the time that on indication device, each sampled point turns over, i.e. series count value;
E: utilize the serial count value obtained, according to the number of divisions on indication device, try to achieve the instantaneous angular velocity of each sampled point on indication device respectively, the set of the instantaneous angular velocity of each sampled point on indication device, constitute series instantaneous angular velocity value;
F: deduct rotating shaft mean angular velocity value from indication device in the instantaneous angular velocity of each sampled point, try to achieve the instantaneous torsional angle velocity amplitude of the equiangular sampling of each sampled point on indication device respectively, the set of the instantaneous torsional angle velocity amplitude of the equiangular sampling of each sampled point on indication device, constitutes the series instantaneous torsional angle velocity amplitude of equiangular sampling;
G: the instantaneous torsional angle velocity amplitude of equiangular sampling of each sampled point is carried out resampling, instantaneous for equiangular sampling torsional angle velocity amplitude is converted into the constant duration instantaneous torsional angle velocity amplitude of sampling of each sampled point obtained with constant duration sample mode, the set of the constant duration instantaneous torsional angle velocity amplitude of sampling of each sampled point obtained, constitutes the series constant duration instantaneous torsional angle velocity amplitude of sampling;
H: utilize least square method to eliminate the trend term in the series constant duration instantaneous torsional angle velocity amplitude of sampling, try to achieve true torsional angle rate signal;
I: the true torsional angle rate signal after eliminating trend term is integrated conversion, obtains torsional oscillation angular displacement signal.
2. rotating machinery torsion vibration signal capturing analysis method according to claim 1, it is characterised in that: in described step B, utilize high pass filter to remove the low-frequency disturbance in the pulse signal that detection equipment exports.
3. rotating machinery torsion vibration signal capturing analysis method according to claim 1, it is characterised in that: in described step C, the pulse signal after removal low-frequency disturbance is converted into the standard TTL pulse signal output of 0-5V.
4. rotating machinery torsion vibration signal capturing analysis method according to claim 1, it is characterized in that: in described step D, utilize high-frequency counter that standard TTL pulse signal is carried out signal sampling, obtain the time that on indication device, each sampled point turns over, namely series count value, is designated as �� t respectively1,��t2,...,��tN, wherein, N represents total number of sample points, i.e. total number of teeth of toothed disc.
5. rotating machinery torsion vibration signal capturing analysis method according to claim 4, it is characterised in that: in described step E, utilize following formula to calculate at tnThe instantaneous angular velocity �� of ith sample point on moment indication devicen(i),
N represents total number of sample points, i.e. total number of teeth of toothed disc; I represents ith sample point on indication device, i.e. the subscript of the i-th tooth on toothed disc, �� tiFor the time that ith sample point turns over.
6. rotating machinery torsion vibration signal capturing analysis method according to claim 5, it is characterised in that: in described step F, utilize following formula to calculate the instantaneous torsional angle velocity amplitude �� of equiangular sampling of ith sample point on indication devicetor(i),
Subscript tor represents torsional angle speed, and subscript n represents that instantaneous angular velocity, i represent ith sample point on indication device,For rotating shaft mean angular velocity value,Calculated by below equation:
ω ‾ = 1 N Σ i = 1 N ω n ( i ) .
7. rotating machinery torsion vibration signal capturing analysis method according to claim 6, it is characterised in that: in described step G, the method utilizing linear interpolation, the instantaneous torsional angle velocity amplitude �� of equiangular sampling to each sampled pointtorI () carries out resampling, be translated into the constant duration instantaneous torsional angle velocity amplitude �� of sampling of each sampled pointtor' (t), in variable, subscript �� represents constant duration sampling;
If sampling time interval is �� t, the sampling time corresponding to sample point m is tm, the time respectively t of 2 of the front and back moment adjacent with the sample point m equiangularly sampled point that sample mode obtainsk,tk+1, then the sample point m obtained is sampled at t=t with constant durationmThe constant duration instantaneous torsional angle velocity amplitude of sampling in momentFor:
ω t o r ′ ( t ) | t = t m = ω t o r ( t ) t = t k + t m - t k t k + 1 - t k [ ω t o r ( t ) | t = t k + 1 - ω t o r ( t ) | t = t k ] ,
Wherein,For sampling the sample point m obtained at t=t with constant durationmThe constant duration instantaneous torsional angle velocity amplitude of sampling in moment,WithThe sample point that respectively equiangularly sample mode the obtains instantaneous angular velocity value at k and k+1 sample point place.
8. the rotating machinery torsion vibration signal capturing analysis method stated according to claim 7, it is characterised in that: in described step H, utilize least square method to eliminate the constant duration instantaneous torsional angle velocity amplitude �� of samplingtorTrend term in ' (t);
The instantaneous torsional angle velocity amplitude �� if constant duration is sampledtorTrend term �� in ' (t)1T () is the quadratic polynomial about time variable t:
��1(t)=at2+ bt+c,
Coefficient a, b, the c in multinomial is solved, from instantaneous angular velocity value �� with method of least squaretor' (t) eliminates trend term ��1T (), tries to achieve true torsional angle rate signal �� "tor(t), �� "tor(t)=�� 'tor(t)-��1(t)��
9. rotating machinery torsion vibration signal capturing analysis method according to claim 8, it is characterised in that: in described step I, to true torsional angle rate signal �� "torT () is integrated, obtain torsional oscillation angular displacement signal �� (t), �� (t)=�� �� "tor(t)dt��
10. rotating machinery torsion vibration signal capturing analysis method as claimed in any of claims 1 to 9, it is characterized in that: the described indication device being provided with sampled point adopts the uniform toothed disc of calibration, and supporting detection equipment adopts eddy current sensor or magnetoresistive transducer.
CN201610011172.8A 2016-01-08 2016-01-08 Rotation machinery torsional vibration signal collection analysis method Pending CN105628176A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610011172.8A CN105628176A (en) 2016-01-08 2016-01-08 Rotation machinery torsional vibration signal collection analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610011172.8A CN105628176A (en) 2016-01-08 2016-01-08 Rotation machinery torsional vibration signal collection analysis method

Publications (1)

Publication Number Publication Date
CN105628176A true CN105628176A (en) 2016-06-01

Family

ID=56043322

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610011172.8A Pending CN105628176A (en) 2016-01-08 2016-01-08 Rotation machinery torsional vibration signal collection analysis method

Country Status (1)

Country Link
CN (1) CN105628176A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108593955A (en) * 2018-04-04 2018-09-28 天津大学 For moment rotatation speed test method in the case of cyclic irregularity
CN109708872A (en) * 2017-10-20 2019-05-03 株洲中车时代电气股份有限公司 A kind of train gear-box shaft coupling method for diagnosing faults, apparatus and system
CN111256971A (en) * 2020-02-19 2020-06-09 东南大学 Turbine small shaft bending vibration and torsional vibration parallel measurement method based on overspeed protection signal
CN111523231A (en) * 2020-04-22 2020-08-11 中国华能集团清洁能源技术研究院有限公司 Subsynchronous oscillation analysis method based on EEMD and Prony method
CN112945551A (en) * 2021-01-27 2021-06-11 重庆大学 Gear ring dynamic deformation detection system and evaluation method
CN113340631A (en) * 2021-05-17 2021-09-03 西安交通大学 Torsional vibration testing device and signal analysis method
CN113358351A (en) * 2021-06-17 2021-09-07 武汉理工大学 Rotating shaft end face torsional vibration extraction method and device based on photogrammetry
CN114088325A (en) * 2021-11-24 2022-02-25 大连海事大学 Shafting multi-measuring-point synchronous torsional vibration high-precision monitoring method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57179628A (en) * 1981-04-28 1982-11-05 Ishikawajima Harima Heavy Ind Co Ltd Torsional vibration measuring apparatus for rotating shaft
US5390545A (en) * 1993-01-11 1995-02-21 Caterpillar Inc. Apparatus for measuring torsional vibrations of rotating machinery
CN101241045A (en) * 2008-03-11 2008-08-13 四方电气(集团)有限公司 Secondary sampling method for speed-measuring fluted disc teeth number of self-adapting large steam turbine generator
CN102183290A (en) * 2011-01-24 2011-09-14 哈尔滨工程大学 Method for improving time frequency conversion accuracy of shaft torsional vibration signals
CN104034407A (en) * 2014-06-10 2014-09-10 国家电网公司 Method for reducing periodic error in rotating machinery torsional vibration signal pulse measuring method
CN104089693A (en) * 2014-04-19 2014-10-08 云南电力试验研究院（集团）有限公司电力研究院 Rotation shaft torsional vibration modality online measurement method for turbo generator set

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57179628A (en) * 1981-04-28 1982-11-05 Ishikawajima Harima Heavy Ind Co Ltd Torsional vibration measuring apparatus for rotating shaft
US5390545A (en) * 1993-01-11 1995-02-21 Caterpillar Inc. Apparatus for measuring torsional vibrations of rotating machinery
CN101241045A (en) * 2008-03-11 2008-08-13 四方电气(集团)有限公司 Secondary sampling method for speed-measuring fluted disc teeth number of self-adapting large steam turbine generator
CN102183290A (en) * 2011-01-24 2011-09-14 哈尔滨工程大学 Method for improving time frequency conversion accuracy of shaft torsional vibration signals
CN104089693A (en) * 2014-04-19 2014-10-08 云南电力试验研究院（集团）有限公司电力研究院 Rotation shaft torsional vibration modality online measurement method for turbo generator set
CN104034407A (en) * 2014-06-10 2014-09-10 国家电网公司 Method for reducing periodic error in rotating machinery torsional vibration signal pulse measuring method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
丁康 等: "提高扭转振动瞬态信号处理精度的方法研究", 《振动工程学报》 *
孙凤芝: "《数值计算方法与实验》", 31 January 2013, 黑龙江大学出版社 *
孙月明等: "《机械振动学（测试与分析）》", 30 November 1991, 浙江大学出版社 *
胡伟: "双质量飞轮动态扭转角度信号测试方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109708872A (en) * 2017-10-20 2019-05-03 株洲中车时代电气股份有限公司 A kind of train gear-box shaft coupling method for diagnosing faults, apparatus and system
CN108593955A (en) * 2018-04-04 2018-09-28 天津大学 For moment rotatation speed test method in the case of cyclic irregularity
CN111256971A (en) * 2020-02-19 2020-06-09 东南大学 Turbine small shaft bending vibration and torsional vibration parallel measurement method based on overspeed protection signal
CN111256971B (en) * 2020-02-19 2021-12-07 东南大学 Turbine small shaft bending vibration and torsional vibration parallel measurement method based on overspeed protection signal
CN111523231A (en) * 2020-04-22 2020-08-11 中国华能集团清洁能源技术研究院有限公司 Subsynchronous oscillation analysis method based on EEMD and Prony method
CN112945551A (en) * 2021-01-27 2021-06-11 重庆大学 Gear ring dynamic deformation detection system and evaluation method
CN113340631A (en) * 2021-05-17 2021-09-03 西安交通大学 Torsional vibration testing device and signal analysis method
CN113358351A (en) * 2021-06-17 2021-09-07 武汉理工大学 Rotating shaft end face torsional vibration extraction method and device based on photogrammetry
CN113358351B (en) * 2021-06-17 2022-10-28 武汉理工大学 Rotating shaft end face torsional vibration extraction method and device based on photogrammetry
CN114088325A (en) * 2021-11-24 2022-02-25 大连海事大学 Shafting multi-measuring-point synchronous torsional vibration high-precision monitoring method

Similar Documents

Publication Publication Date Title
CN105628176A (en) Rotation machinery torsional vibration signal collection analysis method
CN101639395B (en) Improved holographic dynamic balancing method of high-speed main shaft
CN105547698A (en) Fault diagnosis method and apparatus for rolling bearing
CN102507205B (en) Method for checking vibration fault of fan blade of aerial engine
CN105973448A (en) Rotating blade vibration measuring method and system
CN104535323A (en) Locomotive wheelset bearing fault diagnosis method based on angular domain-time domain-frequency domain
EP2259033A2 (en) Lateral, angular and torsional vibration monitoring of rotordynamic systems
CN102183366A (en) Device and method for vibration measurement and failure analysis of rolling bearing
CN105466550B (en) Non-homogeneous lack sampling blade tip-timing vibration signal reconstruction method and its device
CN104215456B (en) Plane clustering and frequency-domain compressed sensing reconstruction based mechanical fault diagnosis method
CN104819841B (en) Built-in-coding-information-based single sensing flexible angle-domain averaging method
CN108844733B (en) Gear state monitoring index extraction method based on KL divergence and root mean square value
CN105675113A (en) Rotating machinery angular domain vibration signal acquisition device and method based on microsensor
CN102590543A (en) Rotary speed measuring device and rotary speed measuring method
Lin et al. A review and strategy for the diagnosis of speed-varying machinery
CN104459187A (en) Device and method for measuring rotating speed of large rotating equipment
Luo et al. Effective and accurate approaches for wind turbine gearbox condition monitoring
CN105372029A (en) Engine rotor blade amplitude data recognition method based on tip timing technology
CN105277362A (en) Gear fault detection method on the basis of multi-position turning angle signals of encoders
Wang et al. A two-stage method using spline-kernelled chirplet transform and angle synchronous averaging to detect faults at variable speed
CN105987809A (en) Centrifugal-compressor semi-open-type impeller crack detection method based on random resonance
CN104090126B (en) A kind of method of testing of accelerometer bandwidth
CN108317052B (en) Method and device for detecting damage factor of gear and wind generating set
CN103821749A (en) On-line diagnosis method of stall and surge of axial fan
CN106644040B (en) Rotating shaft torsional vibration detection method and device based on multiple sensors

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20160601

RJ01 Rejection of invention patent application after publication