CN102175409B - Real-time identification method for oil whirl fault of turbo generator set - Google Patents

Abstract

The invention discloses a real-time identification method for oil whirl fault of a turbo generator set, which belongs to the technical field of vibration state monitoring and fault diagnosing of rotating machinery. The method comprises the following steps: setting a time duration and a stepping length; acquiring the data of a generator set in the starting and accelerating process in real time; calculating and acquiring a vibration amplitude sequence and a vibration frequency sequence from low frequency to high frequency in the shaft relative vibration at the current moment by utilizing a rapid Fourier spectrum analysis method; when the preset duration is reached, obtaining the maximum sequence of the low-frequency vibration amplitude sequence, which is formed by arranged and stored data, and a frequency sequence corresponding to the maximum of the low-frequency vibration amplitude sequence; running a rotating speed frequency sequence; running a vibration amplitude sequence corresponding to the rotating speed frequency; and according to the value of each sequence, calculating and judging whether a support bearing at one side of the generator set rotor has oil whirl fault. By using the method provided by the invention, the real-time automatic online monitoring, analysis and judgment for the generator set rotor oil whirl fault are realized, thereby increasing the accuracy for identifying the generator set rotor oil whirl fault.

Description

Turbo-generator Set oil film whirl fault real-time identification method

Technical field

The invention belongs to rotating machinery vibrating condition monitoring and fault diagnosis technical field, relate in particular to a kind of Turbo-generator Set oil film whirl fault real-time identification method.

Background technology

Oil whirl is to have a strong impact on one of vibration fault of axle system stability, can make the rotor of turbogenerator set generation sudden vibration that is in operation.Axle be oil film in the radial journal bearing except producing the whole rotor of pressure support, also produce damping force and suppress rotor oscillation.Under certain condition, oil film whirl fault can take place in the little radial journal bearing of damping, reduces shaft system of unit stability.If there is enough big damping in the bearing, then rotor can be got back to its normal position and settles out; If the damping in the bearing is too small, then rotor will continue unstability, bigger instability vibration occur.Because the unit area load less (underloading) of bearing, the whirling motion meeting of rotor in bearing shell obviously increases, and low frequency component occurs in the frequency spectrum of vibration this moment, even low frequency component can become predominant frequency.

In the unit boosting velocity procedure, when the working rotor rotating speed when twice rotor first critical rotary speed is following, the low-frequency vibration frequency of vortex motion that oil whirl causes is along with the increase of rotating speed, also increases in proportion thereupon, is approximately the half the of working rotor frequency.To the research of oil film whirl fault, confirm that its vibration performance frequency omega is (0.42～0.48) ω according to both at home and abroad, ω is the running speed frequency of axle system.

At present; Turbo-generator Set oil film whirl fault identification need of work is accomplished through observing three-dimensional spectrogram by the expert with certain on-the-spot vibrating failure diagnosis experience; Objectivity is relatively poor; Subjectivity degree of dependence to the expert is higher, and can't accomplish the real-time automatic on-line monitoring of machine group rotor oil film whirl fault identification, analyzes and differentiate.Therefore, propose a kind of Turbo-generator Set oil film whirl fault real-time identification method and just seem very important.

Summary of the invention

The objective of the invention is to; To the problem that exists in the identification of present Turbo-generator Set oil film whirl fault; A kind of Turbo-generator Set oil film whirl fault real-time identification method is proposed; Real-time automatic on-line monitoring, analysis and the differentiation of realization machine group rotor oil film whirl fault, and then the safe operation of assurance Turbo-generator Set.

Technical scheme is that a kind of Turbo-generator Set oil film whirl fault real-time identification method is characterized in that said method comprises the following steps:

Step 1: the number n of setting data in duration T, stepping length t and the sequence;

Step 2: gather the data in the unit starting boosting velocity procedure in real time, said data comprise the tach signal and the key signal of the relative vibration data of the axle of machine group rotor one side radial journal bearing, rotor;

Step 3: utilize the fast Flourier frequency spectrum analysis method, the axle that calculates current time is vibration amplitude sequence and the vibration frequency sequence from the low frequency to the high frequency in the vibration relatively;

Step 4: from the vibration amplitude sequence, obtain unit operation speed-frequency f ^WfreqCorresponding vibration amplitude A ^Wfreq, storage unit running speed frequency f ^WfreqAnd corresponding vibration amplitude A ^Wfreq

Step 5: from the vibration frequency sequence, intercepting all less than unit operation speed-frequency f ^WfreqFrequency, form the low-frequency vibration frequency sequence

Simultaneously, from the vibration amplitude sequence, intercepting all less than unit operation speed-frequency f ^WfreqThe corresponding amplitude of frequency, form low-frequency vibration amplitude sequence

Storage low-frequency vibration frequency sequence

With low-frequency vibration amplitude sequence

Wherein, i=1,2 ..., n;

Step 6: calculate said low-frequency vibration amplitude sequence In maximal value A ^Mlf, and obtain this maximal value A ^MlfCorresponding frequency f ^MlfStorage A ^MlfAnd f ^Mlf

Step 7: judge whether to reach and set duration T, if then execution in step 8; Otherwise, increase by 1 stepping length t, return step 2;

Step 8: according to the sequencing of storage time, with the low-frequency vibration amplitude sequence of each stepping length storage

In maximal value A ^Mlf, line up low-frequency vibration amplitude sequence maximal value sequence

Low-frequency vibration amplitude sequence with each stepping length storage

In maximal value A ^MlfCorresponding frequency f ^Mlf, line up the corresponding frequency sequence of low-frequency vibration amplitude sequence maximal value

Running speed frequency f with each stepping length storage ^Wfreq, line up the running speed frequency sequence

Running speed frequency f with each stepping length storage ^WfreqCorresponding vibration amplitude A ^Wfreq, line up the corresponding vibration amplitude sequence of running speed frequency

Wherein,

Step 9: calculate low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^MlfAnd low-frequency vibration amplitude sequence maximal value sequence

Unstable state parameter ε ^Amlf

Step 10: calculate the corresponding vibration amplitude sequence of running speed frequency

Entropy E ^Wfreq

Step 11: calculate the frequency corresponding to the maximum amplitude of vibration frequency sequence sequence the frequency of each frequency and operating speed serial

corresponding operating speed frequency ratio

and the results are arranged in the frequency ratio sequence

Step 12: from frequency ratio sequence

, obtain frequency ratio maximal value and frequency ratio minimum value

Step 13: whether deteminate machine group rotor one side radial journal bearing oil film whirl fault takes place.

Said calculating low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^MlfUtilize formula

$E^{\mathrm{mlf}}=\underset{j=1}{\overset{m}{\mathrm{\Σ}}}\left[{\left(A_j^{\mathrm{mlf}}\right)}^2\ln \left({\left(A_j^{\mathrm{mlf}}\right)}^2\right)\right]$

Wherein, During as ;

j=1; 2; 3 ..., m;

Said low-frequency vibration amplitude sequence maximal value sequence

Unstable state parameter ε ^AmlfUtilize formula

ε ^amlf＝(S ^amlf+0.5-μ _A)/σ _A

Wherein, S ^AmlfIt is low-frequency vibration amplitude sequence maximal value sequence

The backward number; μ _AIt is low-frequency vibration amplitude sequence maximal value sequence The theoretical average of backward number, μ _A=m (m-1)/4; σ _AIt is low-frequency vibration amplitude sequence maximal value sequence

Backward count the theoretical standard deviation,

J=1,2,3 ..., m;

The vibration amplitude that said calculating running speed frequency is corresponding

Entropy E ^WfreqUtilize formula

$E^{\mathrm{wfreq}}=\underset{j=1}{\overset{m}{\mathrm{\Σ}}}\left[{\left(A_j^{\mathrm{wfreq}}\right)}^2\ln \left({\left(A_j^{\mathrm{wfreq}}\right)}^2\right)\right]$

Wherein, During as

;

j=1; 2; 3; ..., m,

Whether said deteminate machine group rotor one side radial journal bearing oil film whirl fault takes place specifically is, if low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^Mlf, low-frequency vibration amplitude sequence maximal value sequence

Unstable state parameter ε ^Amlf, the corresponding vibration amplitude sequence of running speed frequency

Entropy E ^Wfreq, the frequency ratio sequence In the frequency ratio maximal value

With the frequency ratio minimum value

Satisfy following condition simultaneously:

1) low-frequency vibration amplitude sequence maximal value sequence

Unstable state parameter ε ^AmlfAbsolute value greater than first setting value;

2) low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^MlfGreater than the corresponding vibration amplitude sequence of running speed frequency

Entropy E ^Wfreq

3) Frequency ratio sequence

the frequency ratio of the maximum value

is less than the second set value;

4) frequency ratio sequence

The frequency ratio of the minimum value

greater than the third set value; it is determined that the unit supports the rotor side bearing oil whirl failure occurs; otherwise, determining the unit supports the rotor side bearing oil whirl failure did not occur.

Said first setting value is N _1-α/2(0,1), promptly probability is the standardized normal distribution variate-value of (1-α/2), wherein, sets α/2=0.5%.

Said second setting value is 0.55.

Said the 3rd setting value is 0.4.

The present invention has realized real-time automatic on-line monitoring, the analysis of machine group rotor oil film whirl fault and has differentiated, improved the accuracy of machine group rotor oil film whirl fault identification.

Description of drawings

Fig. 1 is a Turbo-generator Set oil film whirl fault real-time identification method flow diagram;

Fig. 2 is a Turbo-generator Set oil film whirl fault real-time identification synoptic diagram;

Fig. 3 is low-frequency vibration amplitude sequence maximal value sequence data figure;

Fig. 4 is the corresponding vibration amplitude sequence data figure of running speed frequency.

Embodiment

Below in conjunction with accompanying drawing, preferred embodiment is elaborated.Should be emphasized that following explanation only is exemplary, rather than in order to limit scope of the present invention and application thereof.

Fig. 1 is a Turbo-generator Set oil film whirl fault real-time identification method flow diagram.Among Fig. 1, Turbo-generator Set oil film whirl fault real-time identification method provided by the invention comprises the following steps:

Step 1: set duration T=200 second, stepping length t=1 second, the number n of data=50 in the sequence; In addition, first setting value is N _0.995(0,1)=2.5762, the second setting value is that 0.55, the three setting value is 0.4.

Step 2: gather the data in the unit starting boosting velocity procedure in real time, said data comprise the tach signal and the key signal of the relative vibration data of the axle of machine group rotor one side radial journal bearing, rotor.

Fig. 2 is a Turbo-generator Set oil film whirl fault real-time identification synoptic diagram.Among Fig. 2, the key signal that relative vibration signal of turbine generator unit shaft that method provided by the invention needs and analysis of vibration signal processing need can obtain from the supervisory instrument (TSI) of configuration Turbo-generator Set or can obtain from professional vibrating data collection conditioning device.In the present embodiment, the key signal that relative vibration signal of turbine generator unit shaft and analysis of vibration signal processing need obtains from the professional vibrating data collection conditioning device that links to each other with vibration transducer.Among Fig. 2, in the slot that high-speed data acquisition card insertion industrial microcomputer (IPC) provides.Requirement according to high-speed data acquisition card; The key signal that specialty vibrating data collection conditioning device processing relative vibration signal of turbine generator unit shaft and analysis of vibration signal processing need, the high-speed data acquisition card in the key signal input IPC that relative vibration signal of turbine generator unit shaft after treatment and analysis of vibration signal processing need.

Design concrete Turbo-generator Set oil whirl real-time identification program according to the present invention, the real-time analysis program is installed in the industrial microcomputer (IPC).Once diagnosis cyclic process in the Turbo-generator Set oil whirl real-time identification program comprises that the real-time collection boosting velocity procedure data that relate in the diagnostic method, data in real time calculate that storage, boosting velocity procedure low-frequency vibration correlation parameter are calculated in real time, running frequency vibration amplitude entropy calculates in real time, frequency ratio is calculated in real time and series of computation analysis verification link such as oil film whirl fault real time discriminating.Utilize Turbo-generator Set oil whirl real-time identification sequential monitoring to analyze high pressure rotor A side radial journal bearing whether oil film whirl fault takes place.Industrial microcomputer (IPC) adopts the vibration at high speed data collecting card; Gather the data in the unit starting boosting velocity procedure in real time, data comprise the tach signal and the key signal of near the relative vibration data of the axle that records the Turbo-generator Set high pressure rotor A side radial journal bearing, rotor.Each passage technology parameter of vibrating data collection card is 50ks/s, 24bit.

Step 3: utilize the fast Flourier frequency spectrum analysis method, the axle that calculates current time is vibration amplitude sequence and the vibration frequency sequence from the low frequency to the high frequency in the vibration relatively.

The relative vibration data of axle to machine group rotor one side; Utilize fast Flourier (FFT) frequency spectrum analysis method, calculate vibration amplitude sequence (amplitude unit is μ m) and vibration frequency sequence (cps is Hz) from the low frequency to the high frequency in the relative vibration of current time axle.Data in above-mentioned vibration frequency sequence and the vibration amplitude sequence are one to one on order, and promptly the corresponding relation of frequency and amplitude is clear and definite.

Step 4: from the vibration amplitude sequence, obtain unit operation speed-frequency f ^WfreqCorresponding vibration amplitude A ^Wfreq, storage unit running speed frequency f ^WfreqAnd corresponding vibration amplitude A ^Wfreq

In the unit starting process, running speed respective frequencies f ^WfreqBe and increase progressively variation.From the vibration amplitude sequence, obtain unit operation speed-frequency f ^WfreqCorresponding vibration amplitude A ^Wfreq, real-time storage unit operation speed-frequency f ^WfreqAnd corresponding vibration amplitude A ^Wfreq, data are that every separated t=1 stores once second.

Step 5: from the vibration frequency sequence, intercepting all less than unit operation speed-frequency f ^WfreqFrequency, form the low-frequency vibration frequency sequence

Simultaneously, from the vibration amplitude sequence, intercepting all less than unit operation speed-frequency f ^WfreqThe corresponding amplitude of frequency, form low-frequency vibration amplitude sequence

Storage low-frequency vibration frequency sequence

With low-frequency vibration amplitude sequence

Wherein, i=1,2 ..., 50.

Usually, homemade Turbo-generator Set working speed frequency f ^WfreqBe 50Hz.Therefore, from the vibration frequency sequence, intercepting all less than unit working speed frequency f ^WfreqFrequency, promptly from the vibration frequency sequence, all frequencies less than 50Hz of intercepting form the low-frequency vibration frequency sequence

And storage.

From the vibration amplitude sequence, intercepting all less than unit operation speed-frequency f ^WfreqThe corresponding amplitude of frequency, promptly from the vibration amplitude sequence, all frequencies of intercepting form low-frequency vibration amplitude sequence less than the corresponding amplitude of the frequency of 50Hz

And storage.

In the implementation process, vibrating data collection frequency and image data amount be can set, low-frequency vibration frequency sequence, low-frequency vibration amplitude sequence data number n=50 made.

Step 6: calculate said low-frequency vibration amplitude sequence

In maximal value A ^Mlf, and obtain this maximal value A ^MlfCorresponding frequency f ^MlfStorage A ^MlfAnd f ^Mlf

Owing to set n=50; Therefore in low-frequency vibration amplitude sequence

, have 50 data.Calculate the maximal value in 50 data, be low-frequency vibration amplitude sequence

In maximal value A ^MlfIn addition, obtain this maximal value A ^MlfCorresponding frequency f ^MlfAt last, storage low-frequency vibration amplitude sequence

In maximal value A ^MlfAnd corresponding frequency f ^Mlf

Step 7: judge whether to reach and set duration T=200 second, if then execution in step 8; Otherwise, increase by 1 stepping length t=1 second, return step 2.

Step 8:, the data of each stepping length storage are lined up sequence according to the sequencing of storage time.The data that wherein need line up sequence comprise:

A, according to the sequencing of storage time, with the low-frequency vibration amplitude sequence of each stepping length storage

In maximal value A ^Mlf, line up low-frequency vibration amplitude sequence maximal value sequence

Owing to set duration T=200 second, stepping length t=1 second, the low-frequency vibration amplitude sequence that each stepping length t=1 stores second

In maximal value A ^MlfHave 200.With these 200 low-frequency vibration amplitude sequences

In maximal value A ^Mlf, line up sequence according to time storage sequencing, be designated as low-frequency vibration amplitude sequence maximal value sequence

B, with the low-frequency vibration amplitude sequence of each stepping length storage

In maximal value A ^MlfCorresponding frequency f ^Mlf, line up the corresponding frequency sequence of low-frequency vibration amplitude sequence maximal value

Equally, each stepping length t=1 low-frequency vibration amplitude sequence of storing second

In maximal value A ^MlfCorresponding frequency f ^MlfAlso there are 200.With these 200 frequency f ^Mlf, line up sequence according to time storage sequencing, be designated as the corresponding frequency sequence of low-frequency vibration amplitude sequence maximal value

C, with the running speed frequency f of each stepping length storage ^Wfreq, line up the running speed frequency sequence

The running speed frequency f that each stepping length t=1 stores second ^WfreqThere are 200, with these 200 running speed frequency f ^Wfreq, line up sequence according to time storage sequencing, be designated as the running speed frequency sequence

D, with the running speed frequency f of each stepping length storage ^WfreqCorresponding vibration amplitude A ^Wfreq, line up the corresponding vibration amplitude sequence of running speed frequency

The running speed frequency f that each stepping length t=1 stores second ^WfreqCorresponding vibration amplitude A ^WfreqThere are 200, with these 200 vibration amplitude A ^WfreqLine up sequence according to time storage sequencing, be designated as the corresponding vibration amplitude sequence of running speed frequency

In above-mentioned each sequence, j=1,2 ..., 200.

Step 9: calculate low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^MlfAnd low-frequency vibration amplitude sequence maximal value sequence Unstable state parameter ε ^Amlf

Calculate low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^MlfUtilize formula

Wherein, regulation is worked as

The time,

J=1,2,3 ..., m;

Calculate low-frequency vibration amplitude sequence maximal value sequence

Unstable state parameter ε ^AmlfUtilize formula ε ^Amlf=(S ^Amlf+ 0.5-μ _A)/σ _AWherein, S ^AmlfIt is low-frequency vibration amplitude sequence maximal value sequence

The backward number; Backward is to being meant that in a data sequence, the front and back position and the size order of a logarithm are opposite, and promptly the number of front is greater than the number of back; The backward number is meant the right sum of backward in the data sequence.μ _AIt is low-frequency vibration amplitude sequence maximal value sequence

The theoretical average of backward number, μ _A=m (m-1)/4; σ _AIt is low-frequency vibration amplitude sequence maximal value sequence

Backward count the theoretical standard deviation,

J=1,2,3 ..., m;

Step 10: calculate the corresponding vibration amplitude sequence of running speed frequency

Entropy E ^Wfreq

Calculate the corresponding vibration amplitude sequence of running speed frequency Entropy E ^WfreqUtilize formula

Wherein, regulation is worked as

The time, J=1,2,3 ..., m,

Step 11: calculate the frequency corresponding to the maximum amplitude of vibration frequency sequence sequence the frequency of each frequency and operating speed serial

corresponding operating speed frequency ratio

and the result is stored in the frequency ratio sequence.

Using the formula

Calculate the sequence corresponding to the maximum amplitude of low frequency vibration frequency sequence

the frequency of each frequency and operating speed serial

corresponding operating speed frequency ratio.And result of calculation

lined up frequency ratio sequence

wherein; J=1; 2; 3; ..., 200.

Step 12: from frequency ratio sequence , obtain frequency ratio maximal value and frequency ratio minimum value

Frequency ratio maximal value

frequency ratio minimum value

wherein; J=1; 2; 3; ..., 200.

Step 13: whether deteminate machine group rotor one side radial journal bearing oil film whirl fault takes place.

The process whether deteminate machine group rotor one side radial journal bearing oil film whirl fault takes place is that whether deteminate machine group rotor one side radial journal bearing oil film whirl fault takes place specifically is, if low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^Mlf, low-frequency vibration amplitude sequence maximal value sequence

Unstable state parameter ε ^Amlf, the corresponding vibration amplitude sequence of running speed frequency

Entropy E ^Wfreq, the frequency ratio sequence

In the frequency ratio maximal value

With the frequency ratio minimum value

Satisfy following condition simultaneously:

1) low-frequency vibration amplitude sequence maximal value sequence Unstable state parameter ε ^AmlfAbsolute value greater than first setting value;

2) low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^MlfGreater than the corresponding vibration amplitude sequence of running speed frequency

Entropy E ^Wfreq

3) Frequency ratio sequence

the frequency ratio of the maximum value

is less than the second set value;

4) frequency ratio sequence

the frequency ratio of the minimum value greater than the third set value;

Deteminate machine group rotor one side radial journal bearing generation oil film whirl fault then; Otherwise oil film whirl fault does not take place in deteminate machine group rotor one side radial journal bearing.

According to the first setting value N that sets in the step 1 _0.995(0,1)=2.5762, the second setting value 0.55 and the 3rd setting value 0.4 are supposed low-frequency vibration amplitude sequence maximal value sequence Unstable state parameter ε ^Amlf=4.67606; Low-frequency vibration amplitude sequence maximal value sequence

As shown in Figure 3, low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^Mlf=1.06283E+7; The vibration amplitude sequence that the running speed frequency is corresponding

Vibration amplitude sequence as shown in Figure 4, that the running speed frequency is corresponding Entropy E ^Wfreq=3.12365E+6; The frequency ratio sequence In the frequency ratio maximal value

The frequency ratio sequence

In the frequency ratio minimum value

The above results satisfies condition simultaneously | ε ^Amlf|＞N _0.995(0,1)=2.5762, E ^Mlf＞E ^Wfreq,

And

Then can judge: machine group rotor one side radial journal bearing generation oil film whirl fault from result of calculation.

The present invention utilizes the data in the unit starting boosting velocity procedure of real-time collection; Real-time automatic on-line monitoring, analysis and differentiation to machine group rotor oil film whirl fault; Realized the automatic on-line identification of machine group rotor oil film whirl fault; Improve the accuracy of machine group rotor oil film whirl fault identification, guaranteed the safe operation of Turbo-generator Set.

The above; Be merely the preferable embodiment of the present invention, but protection scope of the present invention is not limited thereto, any technician who is familiar with the present technique field is in the technical scope that the present invention discloses; The variation that can expect easily or replacement all should be encompassed within protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with the protection domain of claim.

Claims (7)

1. a Turbo-generator Set oil film whirl fault real-time identification method is characterized in that said method comprises the following steps:

Step 1: the number n of setting data in duration T, stepping length t and the sequence;

Step 2: gather the data in the unit starting boosting velocity procedure in real time, said data comprise the tach signal and the key signal of the relative vibration data of the axle of machine group rotor one side radial journal bearing, rotor;

Step 3: utilize the fast Flourier frequency spectrum analysis method, the axle that calculates current time is vibration amplitude sequence and the vibration frequency sequence from the low frequency to the high frequency in the vibration relatively;

Step 4: from the vibration amplitude sequence, obtain unit operation speed-frequency f ^WfreqCorresponding vibration amplitude A ^Wfreq, storage unit running speed frequency f ^WfreqAnd corresponding vibration amplitude A ^Wfreq

Step 5: from the vibration frequency sequence, intercepting all less than unit operation speed-frequency f ^WfreqFrequency, form the low-frequency vibration frequency sequence

Simultaneously, from the vibration amplitude sequence, intercepting all less than unit operation speed-frequency f ^WfreqThe corresponding amplitude of frequency, form low-frequency vibration amplitude sequence

Storage low-frequency vibration frequency sequence

With low-frequency vibration amplitude sequence

Wherein, i=1,2 ..., n;

Step 6: calculate said low-frequency vibration amplitude sequence

In maximal value A ^Mlf, and obtain this maximal value A ^MlfCorresponding frequency f ^MlfStorage A ^MlfAnd f ^Mlf

Step 7: judge whether to reach and set duration T, if then execution in step 8; Otherwise, increase by 1 stepping length t, return step 2;

Step 8: according to the sequencing of storage time, with the low-frequency vibration amplitude sequence of each stepping length storage

In maximal value A ^Mlf, line up low-frequency vibration amplitude sequence maximal value sequence

Low-frequency vibration amplitude sequence with each stepping length storage In maximal value A ^MlfCorresponding frequency f ^Mlf, line up the corresponding frequency sequence of low-frequency vibration amplitude sequence maximal value

Running speed frequency f with each stepping length storage ^Wfreq, line up the running speed frequency sequence

Running speed frequency f with each stepping length storage ^WfreqCorresponding vibration amplitude A ^Wfreq, line up the corresponding vibration amplitude sequence of running speed frequency

Wherein, $j=\mathrm{1,2},...,\frac{T}{t};$

Step 9: calculate low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^MlfAnd low-frequency vibration amplitude sequence maximal value sequence

Unstable state parameter ε ^Amlf

Step 10: calculate the corresponding vibration amplitude sequence of running speed frequency

Entropy E ^Wfreq

Step 11: calculate the frequency corresponding to the maximum amplitude of vibration frequency sequence sequence

The frequency of each frequency and operating speed serial

corresponding operating speed frequency ratio

and The results are arranged in the frequency ratio of the sequence

Step 12: from frequency ratio sequence

, obtain frequency ratio maximal value

and frequency ratio minimum value

Step 13: whether deteminate machine group rotor one side radial journal bearing oil film whirl fault takes place; Specifically be, if low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^Mlf, low-frequency vibration amplitude sequence maximal value sequence Unstable state parameter ε ^Amlf, the corresponding vibration amplitude sequence of running speed frequency

Entropy E ^Wfreq, the frequency ratio sequence

In the frequency ratio maximal value

With the frequency ratio minimum value

Satisfy following condition simultaneously:

1) low-frequency vibration amplitude sequence maximal value sequence

Unstable state parameter ε ^AmlfAbsolute value greater than first setting value;

2) low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^MlfGreater than the corresponding vibration amplitude sequence of running speed frequency

Entropy E ^Wfreq

3) Frequency ratio sequence

the frequency ratio of the maximum value

is less than the second set value;

4) frequency ratio sequence the frequency ratio of the minimum value

greater than the third set value;

Deteminate machine group rotor one side radial journal bearing generation oil film whirl fault then; Otherwise oil film whirl fault does not take place in deteminate machine group rotor one side radial journal bearing.

2. a kind of Turbo-generator Set oil film whirl fault real-time identification method according to claim 1 is characterized in that said calculating low-frequency vibration amplitude sequence maximal value sequence

Entropy E ^MlfUtilize formula

$E^{\mathrm{mlf}}=\underset{j=1}{\overset{m}{\mathrm{\Σ}}}\left[{\left(A_j^{\mathrm{mlf}}\right)}^2\ln \left({\left(A_j^{\mathrm{mlf}}\right)}^2\right)\right],$

Wherein, when ${\left(A_j^{\mathrm{Mlf}}\right)}^2=0$ The time, $\mathrm{Ln}\left({\left(A_j^{\mathrm{Mlf}}\right)}^2\right)=0;$ J=1,2,3 ..., m; $m=\frac{T}{t}.$

3. a kind of Turbo-generator Set oil film whirl fault real-time identification method according to claim 1 is characterized in that said low-frequency vibration amplitude sequence maximal value sequence Unstable state parameter ε ^AmlfUtilize formula

ε ^amlf=(S ^amlf+0.5-μ _A)σ _A，

Wherein, S ^AmlfIt is low-frequency vibration amplitude sequence maximal value sequence

The backward number; μ _AIt is low-frequency vibration amplitude sequence maximal value sequence

The theoretical average of backward number, μ _A=m (m-1)/4; σ _AIt is low-frequency vibration amplitude sequence maximal value sequence

Backward count the theoretical standard deviation,

J=1,2,3 ..., m; $m=\frac{T}{t}.$

4. a kind of Turbo-generator Set oil film whirl fault real-time identification method according to claim 1 is characterized in that the vibration amplitude that said calculating running speed frequency is corresponding

Entropy E ^WfreqUtilize formula

$E^{\mathrm{wfreq}}=\underset{j=1}{\overset{m}{\mathrm{\Σ}}}\left[{\left(A_j^{\mathrm{wfreq}}\right)}^2\ln \left({\left(A_j^{\mathrm{wfreq}}\right)}^2\right)\right],$

Wherein, when ${\left(A_j^{\mathrm{Wfreq}}\right)}^2=0$ The time, $\mathrm{Ln}\left({\left(A_j^{\mathrm{Wfreq}}\right)}^2\right)=0;$ J=1,2,3 ..., m, $m=\frac{T}{t}.$

5. a kind of Turbo-generator Set oil film whirl fault real-time identification method according to claim 1 is characterized in that said first setting value is N _1-α/2(0,1), promptly probability is the standardized normal distribution variate-value of (1-α/2), wherein, sets α/2=0.5%.

6. a kind of Turbo-generator Set oil film whirl fault real-time identification method according to claim 1 is characterized in that said second setting value is 0.55.

7. a kind of Turbo-generator Set oil film whirl fault real-time identification method according to claim 1 is characterized in that said the 3rd setting value is 0.4.

Images