CN101451898B - Steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method - Google Patents
Steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method Download PDFInfo
- Publication number
- CN101451898B CN101451898B CN2009100766594A CN200910076659A CN101451898B CN 101451898 B CN101451898 B CN 101451898B CN 2009100766594 A CN2009100766594 A CN 2009100766594A CN 200910076659 A CN200910076659 A CN 200910076659A CN 101451898 B CN101451898 B CN 101451898B
- Authority
- CN
- China
- Prior art keywords
- axle
- fundamental vibration
- data
- constantly
- vibration amplitude
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a real time diagnosing method for turbo unit rotor thermal bending malfunctions. The method processes the calculation, analysis and judgments for the vibrating data via collecting the rotor vibrating signals of the turbo unit; and the method processes the real-time calculation of storage on both sides of the rotor shaft relative vibration frequency of vibration amplitude and phase values, and real-time axis vibration frequency of vibration amplitude in real-time authentication, determine any side of the rotor shaft vibration amplitude of the vibration frequency is greater than the vibration amplitude threshold. Then the method is combined with least square method, based on the axial vibration frequency of vibration data axis vibration frequency of the vibration amplitude of the trend of flat-based authentication and shaft vibration frequency vibration phase of verification, such as the trend of moderate real-time quantitative calculation and analysis. Calculated in real-time quantitative analysis based on the combined results of the verification, automatic real-time diagnosis of whether units imbalance rotor thermal bending failure. The invention has the advantages of scientific method; reliable conclusions can automatically real-time online monitoring, fault diagnosis and so on.
Description
Technical field
The invention belongs to rotating machinery vibrating condition monitoring and fault diagnosis field, particularly a kind of steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method of the automatic monitoring of large turbo-type generator group vibrational state real-time online.
Background technology
The flexural deformation that occurs after turbine rotor is heated becomes thermal flexure, and thermal flexure will cause the variation of rotor balancing state, so thermal flexure also becomes thermal unbalance.Thermal flexure is a kind of comparatively common oscillation phenomenon, and the reason that causes thermal flexure is diversified.The defective of rotor material can cause thermal flexure.The inhomogeneous after expansion that causes being heated of material is inhomogeneous, causes rotor thermal bending; The unrelieved stress of rotor in manufacture process, the back stress relief of being heated causes thermal flexure.Usually up and down there is the temperature difference in cylinder, and the upper temp of shutting down the back rotor makes rotor generation thermal flexure than bottom height, needs could recover behind the jiggering of long period.The time of the continuous jiggering of large turbo-type generator group can not be less than the stipulated time, and the jiggering time is too short, and thermal flexure can't recover, and causes that vibration is excessive.Water or the lower steam of temperature from jet chimney, oxygen-eliminating device and sealing system cause the reduction of rotor local temperature, cause rotor thermal bending.The oil that enters in the rotor bore forms the circulation of vaporizing-condense, and rotor is produced asymmetric cooling or heating, forms the thermal flexure of rotor.Loosening rotor and the impeller contact site non-uniform temperature of causing of shrunk-on disc forms thermal flexure.Generator amature causes rotor thermal bending because inhomogeneous cooling is even, internal friction, turn-to-turn short circuit.The work of rotor of steam turbo generator thermal flexure Fault Diagnosis all is to be finished by veteran expert, diagnosis financial cost height, and the cycle is long.Therefore, propose a kind of rotor of steam turbo generator thermal flexure On-line Fault real-time diagnosis method and just seem very important.
Steam-electric generating set rotor thermal bending fault real-time diagnosis method provided by the invention carries out real-time automatic on-line monitoring, analyzes, diagnoses unit rotor thermal bending fault, improves fault diagnosis efficient and accuracy.
Summary of the invention
The purpose of this invention is to provide can automatic on-line monitoring, accurate a kind of steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method of tracing trouble.This method is calculated realization based on the relative vibration amplitude of axle and the phase data of steam turbine operation rotor in conjunction with computer program.
The technical solution used in the present invention is: a kind of steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method is characterized in that it comprises:
(1) data acquisition, record near the radial journal bearing of harvester group rotor both sides in real time axle vibration data and key signal relatively;
(2) shake data in real time computing and storage of axle at the axle of machine group rotor both sides vibration data relatively, utilizes the FFT frequency spectrum analysis method, calculates rotor A, the relative fundamental vibration frequency vibration amplitude of B two side shafts A in real time synchronously
Ra, A
RbWith phase place P
Ra, P
Rb, and storage relative fundamental vibration frequency vibration amplitude of rotor two side shafts and phase value, wherein FFT is a fast fourier transform;
(3) the axle fundamental vibration amplitude real-time verification that shakes carries out current axle shake fundamental vibration amplitude and amplitude thresholds A in real time
T1Comparison, if the axle of any side of current rotor shakes the amplitude of fundamental vibration greater than A
T1, write down this T constantly so
1And carry out follow-up computational analysis.
(4) the axle mild property of the fundamental vibration amplitude trend checking of shaking is from T
1Constantly be truncated to T0 fundamental vibration amplitude (μ m) data constantly forward, calculate T0 constantly to T
1The linear model best-fit slope of fundamental vibration amplitude data is constantly carried out the quadravalence polynomial expression best fit analysis of fundamental vibration amplitude data simultaneously.According to axle shake the slope a of fundamental vibration amplitude linear model best-fit and the square error e1 of quadravalence polynomial expression best-fit, judge whether the axle mild property of the fundamental vibration amplitude trend checking of shaking is passed through.
(5) the axle mild property of the fundamental vibration phase place trend checking of shaking is from T
1Constantly be truncated to the T0 axle constantly fundamental vibration phase data (phase unit for °) of shaking forward, calculate T0 constantly to T
1Fundamental vibration phase data is constantly carried out the best linear fit analysis.According to axle shake the slope m of linear model best-fit of fundamental vibration phase data and the square error e2 of linear model best-fit, judge whether the axle mild property of the fundamental vibration phase place trend checking of shaking is passed through.
(6) rotor thermal bending unbalance fault is judged, in conjunction with axle fundamental vibration amplitude real-time verification, axle the shake result of 3 real-time analysiss of fundamental vibration phase place trend mild property checking of the mild property checking of fundamental vibration amplitude trend and axle that shakes that shakes, judge whether drawing the large turbo-type generator group, rotor thermal bending unbalance fault takes place.
Steam Turbine rotor thermal bending unbalance fault diagnostic method of the present invention utilizes relative vibration amplitude of the axle of unit operation rotor and phase data, obtain the fault diagnosis conclusion through the computational analysis judgement, has methodological science, conclusion is reliable, can realize advantages such as automatic real time on-line monitoring, tracing trouble.
Description of drawings
Fig. 1 is the monitoring synoptic diagram of steam-electric generating set rotor thermal bending unbalance fault.
Fig. 2 is the rotor thermal bending unbalance fault real time diagnostic functional flow diagram.
Fig. 3 is the axle mild property of the fundamental vibration amplitude trend checking process flow diagram that shakes.
Fig. 4 is the axle mild property of the fundamental vibration phase place trend checking process flow diagram that shakes.
Embodiment
The large steam turbine-generator set rotor thermal bending unbalance fault real time diagnostic method that the present invention proposes mainly is made up of data acquisition, axle data in real time computing and storage, axle fundamental vibration amplitude real-time verification, the axle checking of the mild property of fundamental vibration amplitude trend, the axle links such as the checking of the mild property of fundamental vibration phase place trend, rotor thermal bending unbalance fault judgement of shaking of shaking of shaking of shaking, and its functional flow diagram as shown in Figure 2.In the real-time diagnosis process, axle the shake mild property checking of fundamental vibration amplitude trend and axle 3 links of the mild property of fundamental vibration phase place trend checking of shaking of fundamental vibration amplitude real-time verification, axle of shaking are carried out in real time, and in the fault verification link simultaneously according to the diagnostic result of 3 links, guaranteed the reliability of failure diagnostic process and the accuracy of diagnostic result thus.Further specify concrete implementation step and diagnostic method below in conjunction with accompanying drawing.
1. data acquisition
In the monitoring synoptic diagram of steam-electric generating set rotor thermal bending unbalance fault as shown in Figure 1, near the A of rotor side, B side radial journal bearing and key phase gear, place the vibrating data collection sensor, the vibration data signal of vibrating data collection sensor acquisition is transported to by special vibrating data acquisition conditioning device and inserts vibration at high speed data collecting card in the slot that industrial microcomputer (IPC) provides, record near the radial journal bearing of harvester group rotor both sides in real time axle vibration data and key signal relatively.Each passage technology parameter of vibrating data collection card is 50ks/s, 24bit.
2. shake data in real time computing and storage of axle
The relative vibration data of axle at machine group rotor both sides utilizes the FFT frequency spectrum analysis method, carries out real-time synchronometer calculation and Analysis.Calculate rotor A, the relative fundamental vibration frequency vibration amplitude of B two side shafts A in real time synchronously
Ra, A
RbWith phase place P
Ra, P
RbAxle the shake frequency of fundamental frequency working speed correspondence when being meant rotor stable state operate as normal, i.e. 50Hz.The logical frequency vibration amplitude of the relative fundamental vibration frequency vibration amplitude of storage rotor two side shafts and phase value and the relative vibration of axle, data are every storage in 0.1 second once.
3. the axle fundamental vibration amplitude real-time verification that shakes
The setting shaft fundamental vibration amplitude thresholds A that shakes
T1, A
T1=62 μ m carry out current axle shake fundamental vibration amplitude and A in real time
T1Comparison, if the axle of any side of current rotor shakes the amplitude of fundamental vibration greater than A
T1, write down this T constantly so
1, and carry out follow-up computational analysis.The amplitude of fundamental vibration all is less than or equal to A if the axle of rotor both sides shakes
T1, so fault diagnostic program reenter data acquisition, the axle shake the data in real time computing and the storage link.Axle the shake frequency of fundamental frequency working speed correspondence when being meant rotor stable state operate as normal, i.e. 50Hz.
4. the axle mild property of the fundamental vibration amplitude trend checking of shaking
According to the relative fundamental vibration amplitude data of vibrating of T1 moment front axle of storage, from T
1Constantly be truncated to T0 fundamental vibration amplitude data (vibration amplitude unit is μ m) constantly forward, | T
1-T0|=P
T01, P
T01Be Preset Time segment length, P
T01=3600 seconds.Axle the shake frequency of fundamental frequency working speed correspondence when being meant rotor stable state operate as normal, i.e. 50Hz.The axle shake fundamental vibration amplitude data be every 0.1 second the storage once.At T0 constantly to T
1Fundamental vibration amplitude data are constantly carried out the best linear fit analysis.Calculate T0 constantly to T
1The best straight line slope of representative input data on the least square method meaning of fundamental vibration amplitude data constantly.The axle fundamental vibration amplitude data of shaking are fitted form into formula 1.
f＝ax+b ……(1)
Wherein, f is the axle fundamental vibration amplitude best linear fit value of shaking, the sequence X that x is made of constantly axle vibrational data acquisition, and a is a slope, b is an intercept.Adopt least square method to carry out match, be about to according to formula
The square error of computational data minimizes square error, obtains the slope a of linear model.Wherein, N is the shake data number of fundamental vibration amplitude data Y of axle, f
iBe i element of best linear fit, y
iBe shake i element of fundamental vibration amplitude data Y of axle.
T according to storage
1The fundamental vibration amplitude data of the relative vibration of front axle constantly are from T
1Constantly be truncated to T forward
0Fundamental vibration amplitude data (vibration amplitude unit is μ m) constantly, | T
1-T
0|=P
T01, P
T01Be Preset Time segment length, P
T01=3600 seconds.The axle shake fundamental vibration amplitude data be every 0.1 second the storage once.At T
0Constantly to T
1Fundamental vibration amplitude data are constantly carried out quadravalence polynomial expression best fit analysis on the least square method meaning.The axle fundamental vibration amplitude of shaking is fitted form into formula 2.
h＝a
0+a
1·x
1+a
2·x
2+a
3·x
3+a
4·x
4 ……(2)
Wherein, h is the axle fundamental vibration amplitude quadravalence polynomial expression best-fit values of shaking, the sequence X that x is made of constantly axle vibrational data acquisition, a
0, a
1, a
2, a
3, a
4It is polynomial each the rank coefficient of quadravalence.Adopt least square method to carry out match, be about to according to formula
The square error e1 of computational data minimizes square error e1.Wherein, N is the shake data number of fundamental vibration amplitude data Y of axle, h
iBe i element of quadravalence polynomial expression best-fit, y
iBe shake i element of fundamental vibration amplitude data Y of axle.
If satisfy following two conditions simultaneously, judge that so the axle mild property of the fundamental vibration amplitude trend checking of shaking passes through.Two conditions comprise: the shake slope a of fundamental vibration amplitude linear model of (1) axle falls into the slope range interval [a of setting
Min, a
Max] in, i.e. a
Min≤ a≤a
Max, a wherein
Min=0.005, a
Max=0.03; (2) the shake square error e1 of fundamental vibration amplitude quadravalence polynomial expression best-fit of axle falls into the scope interval [e1 of setting
Min, e1
Max] in, i.e. e1
Min≤ e1≤e1
Max, e1 wherein
Min=0, e1
Max=20.Its FB(flow block) as shown in Figure 3.
5. the axle mild property of the fundamental vibration phase place trend checking of shaking
According to the T1 of the storage front axle fundamental vibration phase data of vibration relatively constantly, be truncated to constantly fundamental vibration phase data of T0 (phase unit for °) constantly forward from T1, | T
1-T
0|=P
T01, P
T01Be Preset Time segment length, P
T01=3600 seconds.Axle the shake frequency of fundamental frequency working speed correspondence when being meant rotor stable state operate as normal, i.e. 50Hz.The axle shake the fundamental vibration phase data be every 0.1 second the storage once.At T
0Constantly to T
1Fundamental vibration phase data is constantly carried out the best linear fit analysis.Calculate T
0Constantly to T
1The best straight line slope of representative input data on the least square method meaning of fundamental vibration phase data constantly.The axle fundamental vibration phase place g data of shaking are fitted form into formula 3.
g＝mx+n ……(3)
Wherein, g is the axle fundamental vibration phase place best linear fit value of shaking, the sequence X that x is made of constantly axle vibrational data acquisition, and m is a slope, n is an intercept.Adopt least square method to carry out match, be about to according to formula
The square error e2 of computational data minimizes square error e2, obtains the slope m of linear model.Wherein, N is the shake data number of fundamental vibration phase data P of axle, g
iBe i element of best linear fit, p
iBe shake i the element of fundamental vibration phase data P of axle.
If satisfy following two conditions simultaneously, judge that so the axle mild property of the fundamental vibration phase place trend checking of shaking passes through.Two conditions comprise: the shake slope m of fundamental vibration phase linearity model of (1) axle falls into the slope range interval [m of setting
Min, m
Max] in, i.e. m
Min≤ m≤m
Max, m wherein
Min=-0.003, m
Max=0.003; (2) the shake square error e2 of fundamental vibration phase linearity model best-fit of axle falls into the scope interval [e2 of setting
Min, e2
Max] in, i.e. e2
Min≤ e2≤e2
Max, e2 wherein
Min=0, e2
Max=20.Its FB(flow block) as shown in Figure 4.
6. rotor thermal bending unbalance fault is judged
According to above-mentioned axle fundamental vibration amplitude real-time verification, axle the shake result of 3 real-time analysiss of fundamental vibration phase place trend mild property checking of the mild property checking of fundamental vibration amplitude trend and axle that shakes that shakes, judge whether rotor thermal bending unbalance fault takes place.If satisfy above-mentioned 3 checkings simultaneously, can failure judgement take place so.
Embodiment
Utilize this method can realize real-time monitoring, analysis, diagnosis to steam-electric generating set rotor thermal bending unbalance fault.Relative vibration signal of turbine generator unit shaft that diagnostic method needs and analysis of vibration signal are handled the key signal that needs and can be obtained from the supervisory instrument (TSI) of configuration Turbo-generator Set or can obtain from professional vibrating data collection conditioning device.In the present embodiment, relative vibration signal of turbine generator unit shaft and analysis of vibration signal are handled the key signal that needs and are obtained from the professional vibrating data collection conditioning device that links to each other with vibration transducer.The monitoring synoptic diagram of steam-electric generating set rotor thermal bending unbalance fault as shown in Figure 1, high-speed data acquisition card inserts in the slot that industrial microcomputer (IPC) provides.Requirement according to high-speed data acquisition card, specialty vibrating data collection conditioning device handles the relative vibration signal of turbine generator unit shaft and analysis of vibration signal is handled the key signal that needs, and relative vibration signal of turbine generator unit shaft after treatment and analysis of vibration signal are handled the high-speed data acquisition card in the key signal input IPC that needs.According to the concrete unit rotor thermal bending unbalance fault computer diagnosis program of this method design, fault diagnostic program is installed in the industrial microcomputer (IPC).Once diagnosis cyclic process in the unit rotor thermal bending unbalance fault real time diagnostic program comprises the data acquisition that relates in the diagnostic method, axle data in real time computing and storage, axle fundamental vibration amplitude real-time verification, the axle checking of the mild property of fundamental vibration amplitude trend, the axle series of computation analysis verification links such as mild property checking of fundamental vibration phase place trend and rotor thermal bending unbalance fault judgement of shaking of shaking of shaking of shaking.
At first, industrial microcomputer (IPC) is gathered the key signal of relative vibration signal of turbine generator unit shaft and analysis of vibration signal processing needs in real time by high-speed data acquisition card.
Whether suppose program monitoring, diagnosing low pressure rotor thermal bending unbalance fault takes place.At the low pressure rotor A of unit, the relative vibration data of axle of B both sides, utilize FFT (fast fourier transform) frequency spectrum analysis method, calculate the relative fundamental vibration frequency vibration amplitude of rotor two side shafts A in real time synchronously
Ra, A
RbWith phase place P
Ra, P
RbAxle the shake frequency of fundamental frequency working speed correspondence when being meant rotor stable state operate as normal, i.e. 50Hz.Storage relative fundamental vibration frequency vibration amplitude of rotor two side shafts and phase value, data are every storage in 0.1 second once.
Fault diagnostic program carries out current axle shake fundamental vibration amplitude and amplitude thresholds A in real time
T1Comparison, the setting shaft fundamental vibration amplitude thresholds A that shakes
T1, A
T1=62 μ m, the amplitude of fundamental vibration is greater than A if the axle of any side of current rotor shakes
T1, write down this T constantly so
1And carry out follow-up computational analysis.The amplitude of fundamental vibration all is less than or equal to A if the axle of rotor both sides shakes
T1, fault diagnostic program can not enter follow-up analyzing and diagnosing link so, reenter data acquisition, the axle shake the data in real time computing and the storage link.Suppose that the current axle fundamental vibration amplitude that shakes of low pressure rotor A side is 70 μ m, write down this T constantly
1And carry out follow-up computational analysis.
Fault diagnostic program is at the relative vibration signal of a side shaft of low pressure rotor, carries out the axle mild property checking of fundamental vibration amplitude trend and axle 2 checkings of the mild property of fundamental vibration phase place trend checking of shaking of shaking, and 2 proof procedures are to carry out synchronously in real time.Any one checking in 2 checkings was lost efficacy, and all can cause program to enter next diagnostic analysis circulation.
In axle shakes the checking of the mild property of fundamental vibration amplitude trend, from T
1Constantly be truncated to T forward
0Fundamental vibration amplitude data (vibration amplitude unit is μ m) constantly, | T
1-T
0|=P
T01, P
T01Be Preset Time segment length, P
T01=3600 seconds.At T
0Constantly to T
1Fundamental vibration amplitude data are constantly carried out the best linear fit analysis.Calculate T
0Constantly to T
1The best straight line slope of representative input data on the least square method meaning of fundamental vibration amplitude data constantly.The axle fundamental vibration amplitude f data of shaking are fitted form into formula 1.
f＝ax+b ……(1)
Wherein, f is the axle fundamental vibration amplitude best linear fit value of shaking, the sequence X that x is made of constantly axle vibrational data acquisition, and a is a slope, b is an intercept.Adopt least square method to carry out match, be about to according to formula
The square error of computational data minimizes square error, obtains the slope a of linear model.Wherein, N is the shake data number of fundamental vibration amplitude data Y of axle, f
iBe i element of best linear fit, y
iBe shake i element of fundamental vibration amplitude data Y of axle.
Simultaneously, at T
0Constantly to T
1Fundamental vibration amplitude data are constantly carried out quadravalence polynomial expression best fit analysis on the least square method meaning.The axle fundamental vibration amplitude data of shaking are fitted form into formula 2.
h＝a
0+a
1·x
1+a
2·x
2+a
3·x
3+a
4·x
4 ……(2)
Wherein, h is the axle fundamental vibration amplitude quadravalence polynomial expression best-fit values of shaking, the sequence X that x is made of constantly axle vibrational data acquisition, a
0, a
1, a
2, a
3, a
4It is polynomial each the rank coefficient of quadravalence.Adopt least square method to carry out match, be about to according to formula
The square error e1 of computational data minimizes square error e1.Wherein, N is the shake data number of fundamental vibration amplitude data Y of axle, f
iBe i element of quadravalence polynomial expression best-fit, y
iBe shake i element of fundamental vibration amplitude data Y of axle.
If satisfy following two conditions simultaneously, judge that so the axle mild property of the fundamental vibration amplitude trend checking of shaking passes through.Two conditions comprise: the shake slope a of fundamental vibration amplitude linear model of (1) axle falls into the slope range interval [a of setting
Min, a
Max] in, i.e. a
Min≤ a≤a
Max, a wherein
Min=0.005, a
Max=0.03; (2) the shake square error e1 of fundamental vibration amplitude quadravalence polynomial expression best-fit of axle falls into the scope interval [e1 of setting
Min, e1
Max] in, i.e. e1
Min≤ e1≤e1
Max, e1 wherein
Min=0, e1
Max=20.
Suppose that the shake slope of fundamental vibration amplitude linear model of the current axle of low pressure rotor A side is 0.015 to fall into the slope range interval [a of setting
Min, a
Max]; Simultaneously, the shake square error of fundamental vibration amplitude quadravalence polynomial expression best-fit of axle is 15, falls into the scope interval [e1 of setting
Min, e1
Max], judge that the axle mild property of the fundamental vibration amplitude trend checking of shaking passes through.
In axle shakes the checking of the mild property of fundamental vibration phase place trend, from T
1Constantly be truncated to T forward
0Fundamental vibration phase data constantly (phase unit be °), | T
1-T
0|=P
T01, P
T01Be Preset Time segment length, P
T01=3600 seconds, at T
0Constantly to T
1Fundamental vibration phase data is constantly carried out the best linear fit analysis; Calculate T
0Constantly to T
1The best straight line slope of representative input data on the least square method meaning of fundamental vibration phase data constantly; The axle fundamental vibration phase place of shaking is fitted form into formula 3,
g＝mx+n ……(3)
Wherein, g is the axle fundamental vibration phase place best linear fit value of shaking, the sequence X that x is made of constantly axle vibrational data acquisition, and m is a slope, n is an intercept.Adopt least square method to carry out match, be about to according to formula
The square error e2 of computational data minimizes square error e2, obtains the slope m of linear model.Wherein, N is the shake data number of fundamental vibration phase data P of axle, g
iBe i element of best linear fit, p
iBe shake i the element of fundamental vibration phase data P of axle.
If satisfy following two conditions simultaneously, judge that so the axle mild property of the fundamental vibration phase place trend checking of shaking passes through.Two conditions comprise: the shake slope m of fundamental vibration phase linearity model of (1) axle falls into the slope range interval [m of setting
Min, m
Max] in, i.e. m
Min≤ m≤m
Max, m wherein
Min=-0.003, m
Max=0.003; (2) the shake square error e2 of fundamental vibration phase linearity model best-fit of axle falls into the scope interval [e2 of setting
Min, e2
Max] in, i.e. e2
Min≤ e2≤e2
Max, e2 wherein
Min=0, e2
Max=20.
Suppose that the shake slope of fundamental vibration phase linearity model of the current axle of low pressure rotor A side is 0.001 to fall into the slope range interval [m of setting
Min, m
Max]; Simultaneously, the shake square error of fundamental vibration amplitude quadravalence polynomial expression best-fit of axle is 15, falls into the scope interval [e2 of setting
Min, e2
Max].Judge that the axle mild property of the fundamental vibration phase place trend checking of shaking passes through.
At last, fault diagnostic program judges whether to take place rotor thermal bending unbalance fault according to axle the shake result of the mild property of fundamental vibration phase place trend checking of the mild property checking of fundamental vibration amplitude trend and axle that shakes.If satisfy above-mentioned 2 checkings simultaneously, can judge the generation rotor thermal bending unbalance fault so.According to current supposed situation, low pressure rotor A side satisfies the axle mild property checking of fundamental vibration amplitude trend and axle 2 checkings of the mild checking of fundamental vibration phase place trend of shaking of shaking simultaneously, therefore can judge low pressure rotor generation rotor thermal bending unbalance fault.Above-mentioned series of computation analysis verification link is carried out in the diagnostic routine circulation, judges whether current unit rotor thermal bending unbalance fault takes place, and realizes the real-time diagnosis of rotor thermal bending unbalance fault.
Claims (6)
1. steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method is characterized in that it comprises:
(1) data acquisition, near the A of rotor side, B side radial journal bearing and key phase gear, place the vibrating data collection sensor, the vibration data signal of vibrating data collection sensor acquisition is transported to the vibration at high speed data collecting card that inserts in the industrial microcomputer slot by special vibrating data acquisition conditioning device, record near the radial journal bearing of harvester group rotor both sides in real time axle vibration data and key signal relatively;
(2) shake data in real time computing and storage of axle at the axle of machine group rotor both sides vibration data relatively, utilizes the FFT frequency spectrum analysis method, calculates rotor A, the relative fundamental vibration frequency vibration amplitude of B two side shafts A in real time synchronously
Ra, A
RbWith phase place P
Ra, P
Rb, and storage relative fundamental vibration frequency vibration amplitude of rotor two side shafts and phase value, wherein FFT is a fast fourier transform;
(3) the axle fundamental vibration amplitude real-time verification that shakes carries out current axle shake fundamental vibration amplitude and amplitude thresholds A in real time
T1Comparison, if the axle of any side of current rotor shakes the amplitude of fundamental vibration greater than A
T1, write down this T constantly so
1And carry out follow-up computational analysis;
(4) the axle mild property of the fundamental vibration amplitude trend checking of shaking is from T
1Constantly be truncated to T forward
0Fundamental vibration amplitude data are constantly calculated T
0Constantly to T
1The linear model best-fit slope of fundamental vibration amplitude data constantly, carry out the quadravalence polynomial expression best fit analysis of fundamental vibration amplitude data simultaneously, according to axle shake the slope a of fundamental vibration amplitude linear model best-fit and the square error e1 of quadravalence polynomial expression best-fit, judge whether the axle mild property of the fundamental vibration amplitude trend checking of shaking is passed through;
(5) the axle mild property of the fundamental vibration phase place trend checking of shaking is from T
1Constantly be truncated to T forward
0The fundamental vibration phase data of shaking of axle is constantly calculated T
0Constantly to T
1Fundamental vibration phase data is constantly carried out the best linear fit analysis; According to axle shake the slope m of linear model best-fit of fundamental vibration phase data and the square error e2 of linear model best-fit, judge whether the axle mild property of the fundamental vibration phase place trend checking of shaking is passed through;
(6) rotor thermal bending unbalance fault is judged, in conjunction with axle fundamental vibration amplitude real-time verification, axle the shake result of 3 real-time analysiss of fundamental vibration phase place trend mild property checking of the mild property checking of fundamental vibration amplitude trend and axle that shakes that shakes, judge whether drawing the large turbo-type generator group, rotor thermal bending unbalance fault takes place.
2. according to the described steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method of claim 1, it is characterized in that each passage technology parameter of described vibrating data collection card is 50ks/s, 24bit.
3. according to the described steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method of claim 1, it is characterized in that, shake data in real time computing and the storage of described axle is at the axle of the selected rotor both sides of unit vibration data relatively, utilize the FFT frequency spectrum analysis method, calculate the relative fundamental vibration frequency vibration amplitude of rotor two side shafts A in real time synchronously
Ra, A
RbWith phase place P
Ra, P
RbStorage relative fundamental vibration frequency vibration amplitude of rotor two side shafts and phase value were every 0.1 second storage one secondary data.
4. according to the described steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method of claim 1, it is characterized in that the described axle fundamental vibration amplitude real-time verification that shakes is to carry out current axle shake fundamental vibration amplitude and threshold value A in real time
T1Comparison, the setting shaft fundamental vibration amplitude thresholds A that shakes
T1, A
T1=62 μ m, the amplitude of fundamental vibration is greater than A if the axle of any side of current rotor shakes
T1, write down this T constantly so
1And carry out follow-up computational analysis.
5. according to the described steam-electric generating set rotor thermal bending unbalance fault real-time online of claim 1 diagnostic method, it is characterized in that the described mild property of the fundamental vibration amplitude trend checking of shaking is from T
1Constantly be truncated to T forward
0Fundamental vibration amplitude data constantly, | T
1-T
0|=P
T01, P
T01Be Preset Time segment length, P
T01=3600 seconds;
At T
0Constantly to T
1Fundamental vibration amplitude data are constantly carried out the best linear fit analysis, calculate T
0Constantly to T
1The best straight line slope a of representative input data on the least square method meaning of fundamental vibration amplitude data constantly;
At T
0Constantly to T
1Fundamental vibration amplitude data are constantly carried out quadravalence polynomial expression best fit analysis on the least square method meaning, according to formula
The square error e1 of computational data minimizes square error e1; Wherein, N is the shake data number of fundamental vibration amplitude data Y of axle, h
iBe i element of quadravalence polynomial expression best-fit, y
iBe shake i element of fundamental vibration amplitude data Y of axle;
If satisfy following two conditions simultaneously, judge that so the axle mild property of the fundamental vibration amplitude trend checking of shaking passes through, two conditions comprise: the shake slope a of fundamental vibration amplitude linear model of (1) axle falls into the slope range interval [a of setting
Min, a
Max] in, i.e. a
Min≤ a≤a
Max, a wherein
Min=0.005, a
Max=0.03; (2) the shake square error e1 of fundamental vibration amplitude quadravalence polynomial expression best-fit of axle falls into the scope interval [e1 of setting
Min, e1
Max] in, i.e. e1
Min≤ e1≤e1
Max, e1 wherein
Min=0, e1
Max=20.
6. according to the described steam-electric generating set rotor thermal bending unbalance fault real-time online of claim 1 diagnostic method, it is characterized in that the described mild property of the fundamental vibration phase place trend checking of shaking is from T
1Constantly be truncated to T forward
0Fundamental vibration phase data constantly, | T
1-T
0|=P
T01, P
T01Be Preset Time segment length, P
T01=3600 seconds.At T
0To T1 fundamental vibration phase data constantly, carry out the best linear fit analysis constantly; Calculate T
0Constantly to T
1The best straight line slope m of representative input data on the least square method meaning of fundamental vibration phase data constantly is according to formula
The square error e2 of computational data minimizes square error e2, obtains the slope m of linear model; Wherein, N is the shake data number of fundamental vibration phase data P of axle, g
iBe i element of best linear fit, p
iBe shake i the element of fundamental vibration phase data P of axle;
If satisfy following two conditions simultaneously, judge that so the axle mild property of the fundamental vibration phase place trend checking of shaking passes through, two conditions comprise: the shake slope m of fundamental vibration phase linearity model of (1) axle falls into the slope range interval [m of setting
Min, m
Max] in, i.e. m
Min≤ m≤m
Max, m wherein
Min=-0.003, m
Max=0.003; (2) the shake square error e2 of fundamental vibration phase linearity model best-fit of axle falls into the scope interval [e2 of setting
Min, e2
Max] in, i.e. e2
Min≤ e2≤e2
Max, e2 wherein
Min=0, e2
Max=20.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100766594A CN101451898B (en) | 2009-01-13 | 2009-01-13 | Steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100766594A CN101451898B (en) | 2009-01-13 | 2009-01-13 | Steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101451898A CN101451898A (en) | 2009-06-10 |
| CN101451898B true CN101451898B (en) | 2011-01-05 |
Family
ID=40734264
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009100766594A Expired - Fee Related CN101451898B (en) | 2009-01-13 | 2009-01-13 | Steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101451898B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102788662A (en) * | 2011-05-17 | 2012-11-21 | 大连光洋科技工程有限公司 | Vertical rotator dynamic balance detection system without additional sensor |
Families Citing this family (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102087139B (en) * | 2010-11-24 | 2012-10-10 | 华北电力大学 | Method for analyzing frequency components of low-frequency vibration of steam turbine generator unit in real time |
| CN102175440B (en) * | 2011-01-24 | 2012-08-29 | 华北电力大学 | Method for identifying hot bending failure of large-sized generator set in real time |
| CN102788663B (en) * | 2011-05-17 | 2016-01-20 | 大连科德数控有限公司 | Without the horizontal rotation body dynamic balance detection system of additional sensor |
| CN102790476B (en) * | 2011-05-17 | 2016-08-24 | 科德数控股份有限公司 | Dynamic balance analysis device |
| CN102288911B (en) * | 2011-06-03 | 2014-04-02 | 中国神华能源股份有限公司 | System and method for monitoring and diagnosing thermal generator set |
| CN102937460A (en) * | 2012-11-06 | 2013-02-20 | 昆山北极光电子科技有限公司 | Fault diagnosis system for rotating machinery |
| CN103698088B (en) * | 2013-11-14 | 2016-08-17 | 哈尔滨电机厂有限责任公司 | The method of testing of turbogenerator shaft asymmetric stiffness |
| CN103913272B (en) * | 2014-04-03 | 2016-05-18 | 华北电力大学 | A kind of quantitative Diagnosis method of rotor of turbogenerator set quality imbalance fault |
| CN104314617B (en) * | 2014-08-20 | 2016-02-17 | 中国神华能源股份有限公司 | Suppress or stop the gradual bending balance method of steam turbine High inter case |
| CN104215395B (en) * | 2014-09-09 | 2016-11-23 | 中国石油大学(北京) | The detection method of a kind of rotor unbalance fault and device |
| CN106199186B (en) * | 2015-05-05 | 2020-01-31 | 福建宁德核电有限公司 | phase testing device |
| CN105043667B (en) * | 2015-07-01 | 2017-10-20 | 广东石油化工学院 | Rotor unbalance vibration signal amplitude, phase real-time computing technique |
| US9978248B2 (en) * | 2015-10-02 | 2018-05-22 | Skf Canada Limited | Low cost online condition monitoring |
| CN105758460B (en) * | 2016-04-05 | 2017-11-03 | 西安交通大学 | A kind of disc type pull bar combined rotor thermal bending deformation and vibration-testing testing stand |
| CN106768260B (en) * | 2016-12-14 | 2019-05-31 | 广东石油化工学院 | It can inhibit the vibration signal maximum power frequency component real time detection algorithm of direct current disturbance |
| CN110261117A (en) * | 2018-03-12 | 2019-09-20 | 中国航发商用航空发动机有限责任公司 | Turbogenerator simulation experiment system |
| CN109855801B (en) * | 2019-03-29 | 2020-12-15 | 国网山东省电力公司电力科学研究院 | Quick dynamic balance method of pumped storage unit by combining accurate test weight with influence coefficient method |
| CN110907031A (en) * | 2019-12-04 | 2020-03-24 | 江苏方天电力技术有限公司 | Visual monitoring method for vibration variation of phase modulator based on data statistical analysis |
| CN111367260B (en) * | 2020-03-20 | 2022-06-17 | 上海电力大学 | Fault diagnosis device and method for rotor unbalance |
| CN112541155B (en) * | 2020-09-30 | 2022-06-14 | 上海电力大学 | Steam turbine rotor unbalance detection method based on confidence rule base |
| CN112287552B (en) * | 2020-10-30 | 2022-12-06 | 中国航空工业集团公司西安航空计算技术研究所 | Aero-engine vibration trend analysis method |
| CN112857808B (en) * | 2020-12-31 | 2024-02-23 | 沈阳德瓦特汽轮动力有限责任公司 | Verification method and device for steam turbine rotor modal analysis |
| CN113092116B (en) * | 2021-05-20 | 2022-10-28 | 江苏金沃伺服冲床有限公司 | Servo motor fault detection device for servo machinery and use method thereof |
| CN113435395B (en) * | 2021-07-13 | 2023-04-07 | 重庆大学 | Wind wheel fault self-adaptive diagnosis method based on self-adaptive Kalman filter |
| CN114089036B (en) * | 2021-11-09 | 2024-04-16 | 许继集团有限公司 | Key phase pulse intelligent acquisition device and acquisition method thereof |
| CN114322746A (en) * | 2021-12-27 | 2022-04-12 | 贵州电网有限责任公司 | Online testing device and method for flatness of end face of rotating machinery shaft vibration measurement |
| CN115376301A (en) * | 2022-07-25 | 2022-11-22 | 广东粤电博贺能源有限公司 | Blade fracture early warning method and system of steam turbine, electronic equipment and storage medium |
-
2009
- 2009-01-13 CN CN2009100766594A patent/CN101451898B/en not_active Expired - Fee Related
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102788662A (en) * | 2011-05-17 | 2012-11-21 | 大连光洋科技工程有限公司 | Vertical rotator dynamic balance detection system without additional sensor |
| CN102788662B (en) * | 2011-05-17 | 2016-04-06 | 科德数控股份有限公司 | Without the vertical revolving-body dynamic-balancing detection system of additional sensor |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101451898A (en) | 2009-06-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101451898B (en) | Steam-electric generating set rotor thermal bending unbalance fault real time diagnostic method | |
| CN101738293B (en) | Real-time diagnosis method of original quality imbalance fault of turbo generator set rotor | |
| CN101532911B (en) | Large steam turbine-generator set rotor crack fault real-time diagnosis method | |
| CN101430247B (en) | Real-time diagnosis method for random vibration fault of steam generator set | |
| CN101430239B (en) | Real-time diagnosis method for oil film whirl fault of large steam turbine-generator | |
| JP5897824B2 (en) | Turbomachine risk analysis system and machine-readable medium for storing machine-readable instructions for causing a computer to create an inspection recommendation for the turbomachine | |
| CN101429877B (en) | Real-time diagnosis method for angle misalignment fault of steam-electric generating set coupling | |
| CN102282450A (en) | Method and system for monitoring vibratory phenomena occurring in an aircraft's gas turbine engine in operation | |
| CN101403648A (en) | Steam flow excitation fault real-time diagnosis method for large steam turbine-generator | |
| CN102087139B (en) | Method for analyzing frequency components of low-frequency vibration of steam turbine generator unit in real time | |
| CN105841966A (en) | Turbo generator set vibration fault diagnosis method based on forward reasoning | |
| CN102175306B (en) | Method for identifying oil whipping fault of steam turbine generator unit in real-time | |
| CN102072764A (en) | Low-frequency vibration single-peak real-time analysis method for steam turbine generator set | |
| CN102096771B (en) | Real-time temperature variation tendency analysis method of steam turbine generator unit bearing bush | |
| CN102087140B (en) | Method for analyzing stability of low-frequency vibration main peak frequency of turbo generator set | |
| CN102175409B (en) | Real-time identification method for oil whirl fault of turbo generator set | |
| CN102012263B (en) | Method for identifying turbine unit rotor vibration in-phase component stability in real time | |
| CN102879084B (en) | On-line low-frequency vibration unstable state early warning method for steam turbine generator unit | |
| CN102175307A (en) | Real-time quantitative analysis method for low-frequency vibration spectrum array of steam turbine generator set | |
| CN102692303A (en) | High-efficiency identification method of steam excited vibration fault for steam turbine generator unit | |
| CN102022143B (en) | Method for analyzing low-frequency vibration stability of steam turbine generator unit in real time | |
| CN102879085B (en) | Real-time low-frequency vibration unstable state early warning method for steam turbine generator unit | |
| CN102052964B (en) | Real-time recognition method for vibration opposite-phase vector stability of turbogenerator unit rotor | |
| CN102680243B (en) | Online judgment method for steam flow shock excitation fault of steam turbine generator unit | |
| Ogbonnaya et al. | Optimizing gas turbine rotor shaft fault detection, identification and analysis for effective condition monitoring |
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 | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110105 Termination date: 20140113 |