CN106053034A - Rotary machine holographic diagnostics method based on frequency modulation information reconstruction - Google Patents

Abstract

The present invention relates to a rotary machine holographic diagnostics method based on frequency modulation information reconstruction. The method is characterized by firstly installing two displacement sensors of which the included angle is 90 degrees on a rotor vibration test cross-section, acquiring the displacement vibration signals in the X and Y directions, carrying out the ensemble empirical mode decomposition (EEMD) on the displacement vibration signals, obtaining a plurality of intrinsic mode functions (IMF), then carrying out the spectral analysis on the IMFs, and selecting the IMF<x>(t) and the IMF<y>(t) which contain a rotating frequency; adopting a direct orthogonal method to calculate the instantaneous frequencies IF<x>(t) and the IF<y>(t) of the IMF<x>(t) and the IMF<y>(t) respectively, extracting the frequency, amplitude and phase information of the instantaneous frequencies IF<x>(t) and the IF<y>(t) to synthesize a two-dimensional holospectrum, and effectively identifying the operation state of a rotor. The rotary machine holographic diagnostics method based on frequency modulation information reconstruction of the present invention comprehensively considers the vibration signals of the rotor in the horizontal and vertical directions, fully utilizes the phase information, fuses the frequency, amplitude and phase information of the signals organically, and is visual and complete in diagnosis result.

Description

A kind of rotating machinery holographic diagnostics method based on frequency modulated information reconstruct

Technical field

The present invention relates to rotary machinery fault diagnosis technical field, particularly to a kind of rotation based on frequency modulated information reconstruct Machinery holographic diagnostics method.

Background technology

Large rotating machinery is the key equipment in many industries such as oil, chemical industry, the energy, metallurgy, wherein rotor assembly It it is the core of rotating machinery.The common fault of rotor has imbalance, impact and rub, oil whirl etc., its rotor impact and rub The vibration that fault causes often shows nonlinear characteristic, is the failure mode that a class diagnosis difficulty is bigger.More existing letters Number analysis method tends not to accurately, differentiates intuitively the diagnosis of rotor operation state, particularly early-stage weak fault.

The vibration that rotor impact and rub fault causes is Flexural-Torsional Coupling Vibration, and vibration signal form is frequency-modulated form, abnormal Vibration information can be embodied in the instantaneous frequency of radial vibration signal, instantaneous frequency information is used, and is expected to improve rotor The accuracy of fault diagnosis, benefits to diagnosis rotor fault.Owing to the most on-the-spot slewing run all does not installs torsion Turn vibration measurement device, it is meant that cannot directly obtain the torsion vibration signal of slewing.By contrast, most large-scale rotation Equipment has installed rotor radial vibration monitor system the most, extracts torsion and vibration of rotor information, comprehensively examine from rotor radial vibrates Worry rotor, in vibration both horizontally and vertically, can provide new thinking, have important rotor impact and rub fault diagnosis Engineering and economic benefit.

Summary of the invention

In order to overcome the shortcoming of above-mentioned prior art, it is an object of the invention to provide a kind of based on frequency modulated information reconstruct Rotating machinery holographic diagnostics method, by radial vibration signal extraction torsion and vibration of rotor information, the two dimension of synthesis instantaneous frequency Holographic spectrum, thus realize the efficient diagnosis to rotor impact and rub fault.

For reaching above-mentioned purpose, the technical scheme that the present invention takes is:

A kind of rotating machinery holographic diagnostics method based on frequency modulated information reconstruct, comprises the following steps:

Step one, at the displacement transducer that rotor oscillation testing section two angles of installation are 90 °, it is stipulated that rotor rotates through In journey, the displacement transducer of first process is X-direction, and clockwise another displacement transducer after half-twist is Y-direction, two The mounting condition of individual displacement transducer, physical characteristic keep consistent, gather displacement vibration signal x (t) and the y (t) of rotor；

Step 2, carries out gathering empirical mode EEMD to X, displacement vibration signal x (t) of Y-direction and y (t) respectively and decomposes, Obtaining multiple Intrinsic mode functions IMF, Intrinsic mode functions IMF must is fulfilled for following two condition:

Extreme point number on the most whole signal and zero crossing number is equal or at most difference one；

B. in whole time series, all Local modulus maximas the coenvelope that determines and by all local minizing points Determined by the average of lower envelope be zero；

Step 3, decomposes, to set empirical mode EEMD, each Intrinsic mode functions IMF obtained and carries out spectrum analysis, select Go out and comprise the Intrinsic mode functions IMF turning frequency_x(t)And IMF_y(t)；

Step 4, respectively to Intrinsic mode functions IMF_x(t)And IMF_y(t)Direct orthogonalization method is used to calculate instantaneous frequency IF_x(t)And IF_y(t), calculate process as follows:

Intrinsic mode functions IMF_x(t)Resolve into amplitude-modulated signal A (t) and FM signalThat is:

Order:Orthogonal function is defined from the experience frequency modulation component of signal, as follows:

Therefore, the instantaneous phase of signalCalculated by following formula and obtain:

By to instantaneous phaseDerivation, obtains instantaneous frequency IF_x(t)For:

It is calculated instantaneous frequency IF according still further to said process_y(t)；

Step 5, extracts instantaneous frequency IF_x(t)And IF_y(t)Frequency, amplitude, phase information, synthesize two-dimension holographic spectrum, so Effectively identifying rotor operation state afterwards, detailed process is as follows:

First, with interpolation respectively to X, instantaneous frequency IF of Y-direction_x(t)And IF_y(t)Carry out frequency, amplitude, phase place school Just；Then, amplitude A after correction is extracted_x、A_yWith phase place Merge in units of frequency order, obtain ellipse Equation is:

In formula, α n is the displacement transducer angle with horizontal direction of X-direction；

Change the A at different frequency order_x、A_y、Just obtaining corresponding ellipse, combination obtains two-dimension holographic spectrum；

Finally, calculate eccentricity e at each frequency order, major axis a and the oval information of short axle b, according to oval information Effectively identifying rotor operation state, calculation expression is as follows:

$e=\sqrt{a^2-b^2}/a,$

$a=(\sqrt{p+q}+\sqrt{p-q})/2,$

$b=(\sqrt{p+q}-\sqrt{p-q})/2,$

Wherein,

The invention have the benefit that

Being found by the holographic diagnostics method based on frequency modulated information reconstruct using the present invention to propose, rotor is properly functioning Time, each frequency order ellipse presents one group of the least ellipse of random amplitude；During rotor generation impact and rub fault, turn frequency ellipse Round and high frequency multiplication is oval substantially becomes big, and in the ellipse that one group of eccentricity is the biggest, some even straight line, this feature is permissible Effectively identify rotor impact and rub fault.The invention have the advantages that

A) present invention extracts torsion and vibration of rotor information by rotor radial vibration signal, and then realizes rotor fault and sentence Fixed, improve the utilization rate of rotor radial vibration signal；

B) comparing existing analysis on Torsional Vibration technology, the present invention just can complete without directly installing twisting vibration measurement apparatus additional Twisting vibration information retrieval, has saved hardware cost；

C) the holographic diagnostics method based on frequency modulated information reconstruct that the present invention proposes compensate for traditional analysis in rotation Deficiency in mechanical fault diagnosis, has considered rotor at vibration signal both horizontally and vertically, and has made full use of phase place Information, organically merges the frequency of signal, amplitude and phase information, and diagnostic result is the most complete.

Accompanying drawing explanation

Fig. 1 is embodiment 1 rotor-support-foundation system FEM (finite element) model.

Fig. 2 is the instantaneous frequency two-dimension holographic spectrogram of vibration signal under embodiment 1 rotor normal operating condition.

Fig. 3 is the instantaneous frequency two-dimension holographic spectrogram of embodiment 1 rotor impact and rub fault vibration signal.

Fig. 4 is embodiment 2 Bentley RK4 rotor experiment table structural representation.

Fig. 5 is the set empirical mode EEMD decomposition result of embodiment 2 rotor impact and rub fault vibration signal x (t).

Fig. 6 is the set empirical mode EEMD decomposition result of embodiment 2 rotor impact and rub fault vibration signal y (t).

Fig. 7 is embodiment 2 Intrinsic mode functions IMF_x(t)Time domain beamformer and amplitude spectrum.

Fig. 8 is embodiment 2 Intrinsic mode functions IMF_y(t)Time domain beamformer and amplitude spectrum.

Fig. 9 is embodiment 2 instantaneous frequency IF_x(t)Time domain beamformer and amplitude spectrum.

Figure 10 is embodiment 2 instantaneous frequency IF_y(t)Time domain beamformer and amplitude spectrum.

Figure 11 is the instantaneous frequency two-dimension holographic spectrogram of vibration signal under embodiment 2 rotor normal operating condition.

Figure 12 is the instantaneous frequency two-dimension holographic spectrogram of embodiment 2 rotor impact and rub fault vibration signal.

Detailed description of the invention

Describe the present invention below in conjunction with embodiment and accompanying drawing, for the effectiveness of the inventive method is described, point With rotor dynamics emulation signal and experimental signal, the present invention is not verified.

Embodiment 1, uses rotor dynamics emulation signal to verify, sets up rotor-support-foundation system FEM (finite element) model such as figure Shown in 1, rotor impact and rub fault is excited by rotor unbalance, carries out impact and rub fault simulation by changing sound spacing, Rotor-support-foundation system FEM (finite element) model basic parameter is shown in Table 1,

Table 1 rotor-support-foundation system FEM (finite element) model basic parameter table

The element stiffness matrix of rotor-support-foundation system, element mass matrix, unit gyroscopic couple battle array are as follows:

Element stiffness matrix is:

$\&lsqb;M_e\&rsqb;=\left[\begin{array}{cccccccccccc}\frac{EA}{l}& 0& 0& 0& 0& 0& -\frac{EA}{l}& 0& 0& 0& 0& 0\\ 0& \frac{12El}{l^3}& 0& 0& 0& \frac{6El}{l^2}& 0& -\frac{12El}{l^3}& 0& 0& 0& \frac{6El}{l^2}\\ 0& 0& \frac{12El}{l^3}& 0& -\frac{6El}{l^2}& 0& 0& 0& -\frac{12El}{l^3}& 0& -\frac{6El}{l^2}& 0\\ 0& 0& 0& \frac{GJ}{l}& 0& 0& 0& 0& 0& -\frac{GJ}{l}& 0& 0\\ 0& 0& -\frac{6El}{l^2}& 0& \frac{4El}{l}& 0& 0& 0& -\frac{6El}{l^2}& 0& \frac{2El}{l}& 0\\ 0& \frac{6El}{l^2}& 0& 0& 0& \frac{4El}{l}& 0& -\frac{6El}{l^2}& 0& 0& 0& \frac{2El}{l}\\ -\frac{EA}{l}& 0& 0& 0& 0& 0& \frac{EA}{l}& 0& 0& 0& 0& 0\\ 0& -\frac{12El}{l^3}& 0& 0& 0& -\frac{6El}{l^2}& 0& \frac{12El}{l^3}& 0& 0& 0& -\frac{6El}{l^2}\\ 0& 0& -\frac{12El}{l^3}& 0& \frac{6El}{l^2}& 0& 0& 0& \frac{12El}{l^3}& 0& \frac{6El}{l^2}& 0\\ 0& 0& 0& -\frac{GJ}{l}& 0& 0& 0& 0& 0& \frac{GJ}{l}& 0& 0\\ 0& 0& -\frac{6El}{l^2}& 0& \frac{2El}{l}& 0& 0& 0& \frac{6El}{l^2}& 0& \frac{4El}{l}& 0\\ 0& \frac{6El}{l^2}& 0& 0& 0& \frac{2El}{l}& 0& -\frac{6El}{l^2}& 0& 0& 0& \frac{4El}{l}\end{array}\right]$

Element mass matrix is:

$\&lsqb;M_e\&rsqb;=\left[\begin{array}{cccccccccccc}\frac{l\mathrm{\&mu;}}{3}& 0& 0& 0& 0& 0& \frac{l\mathrm{\&mu;}}{6}& 0& 0& 0& 0& 0\\ 0& \frac{13l\mathrm{\&mu;}}{35}& 0& 0& 0& \frac{11l^2\mathrm{\&mu;}}{210}& 0& \frac{9l\mathrm{\&mu;}}{70}& 0& 0& 0& -\frac{13l^2\mathrm{\&mu;}}{420}\\ 0& 0& \frac{13l\mathrm{\&mu;}}{35}& 0& -\frac{11l^2\mathrm{\&mu;}}{210}& 0& 0& 0& \frac{9l\mathrm{\&mu;}}{70}& 0& \frac{13l^2\mathrm{\&mu;}}{420}& 0\\ 0& 0& 0& \frac{{\mathrm{lr}}^2\mathrm{\&mu;}}{6}& 0& 0& 0& 0& 0& \frac{{\mathrm{lr}}^2\mathrm{\&mu;}}{12}& 0& 0\\ 0& 0& -\frac{11l^2\mathrm{\&mu;}}{210}& 0& \frac{l^3\mathrm{\&mu;}}{105}& 0& 0& 0& -\frac{13l^2\mathrm{\&mu;}}{420}& 0& -\frac{l^3\mathrm{\&mu;}}{420}& 0\\ 0& \frac{11l^2\mathrm{\&mu;}}{210}& 0& 0& 0& \frac{l^3\mathrm{\&mu;}}{105}& 0& \frac{13l^2\mathrm{\&mu;}}{420}& 0& 0& 0& -\frac{l^3\mathrm{\&mu;}}{140}\\ \frac{l\mathrm{\&mu;}}{6}& 0& 0& 0& 0& 0& \frac{l\mathrm{\&mu;}}{3}& 0& 0& 0& 0& 0\\ 0& \frac{9l\mathrm{\&mu;}}{70}& 0& 0& 0& \frac{13l^2\mathrm{\&mu;}}{420}& 0& \frac{13l\mathrm{\&mu;}}{35}& 0& 0& 0& -\frac{11l^2\mathrm{\&mu;}}{210}\\ 0& 0& \frac{9l\mathrm{\&mu;}}{70}& 0& -\frac{13l^2\mathrm{\&mu;}}{420}& 0& 0& 0& \frac{13l\mathrm{\&mu;}}{35}& 0& \frac{11l^2\mathrm{\&mu;}}{210}& 0\\ 0& 0& 0& \frac{{\mathrm{lr}}^2\mathrm{\&mu;}}{12}& 0& 0& 0& 0& 0& \frac{{\mathrm{lr}}^2\mathrm{\&mu;}}{6}& 0& 0\\ 0& 0& \frac{13l^2\mathrm{\&mu;}}{420}& 0& -\frac{l^3\mathrm{\&mu;}}{140}& 0& 0& 0& \frac{11l^2\mathrm{\&mu;}}{210}& 0& \frac{l^3\mathrm{\&mu;}}{105}& 0\\ 0& -\frac{13l^2\mathrm{\&mu;}}{420}& 0& 0& 0& -\frac{l^3\mathrm{\&mu;}}{140}& 0& -\frac{11l^2\mathrm{\&mu;}}{210}& 0& 0& 0& \frac{l^3\mathrm{\&mu;}}{105}\end{array}\right]$

Unit gyroscopic couple battle array is:

$\&lsqb;G_e\&rsqb;=-\frac{{\mathrm{\&omega;\&mu;r}}^2}{60l}\left[\begin{array}{cccccccccccc}0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& 36& 0& -3l& 0& 0& 0& -36& 0& 3l& 0\\ 0& -36& 0& 0& 0& -3l& 0& 36& 0& 0& 0& -3l\\ 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 3l& 0& 0& 0& 4l^2& 0& -3l& 0& 0& 0& l^2\\ 0& 0& 3l& 0& -4l^2& 0& 0& 0& -3l& 0& l^2& 0\\ 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& -36& 0& 3l& 0& 0& 0& 36& 0& 3l& 0\\ 0& 36& 0& 0& 0& 3l& 0& -36& 0& 0& 0& 3l\\ 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 3l& 0& 0& 0& -l^2& 0& -3l& 0& 0& 0& 4l^2\\ 0& 0& 3l& 0& l^2& 0& 0& 0& -3l& 0& -4l^2& 0\end{array}\right]$

A kind of rotating machinery holographic diagnostics method based on frequency modulated information reconstruct, comprises the following steps:

Step one, at the displacement transducer that rotor oscillation testing section two angles of installation are 90 °, it is stipulated that rotor rotates through In journey, the displacement transducer of first process is X-direction, and clockwise another displacement transducer after half-twist is Y-direction, two The mounting condition of individual displacement transducer, physical characteristic keep consistent, gather displacement vibration signal x (t) and the y (t) of rotor；

Step 2, carries out gathering empirical mode EEMD to X, displacement vibration signal x (t) of Y-direction and y (t) respectively and decomposes, Obtaining multiple Intrinsic mode functions IMF, Intrinsic mode functions IMF must is fulfilled for following two condition:

Extreme point number on the most whole signal and zero crossing number is equal or at most difference one；

B. in whole time series, all Local modulus maximas the coenvelope that determines and by all local minizing points Determined by the average of lower envelope be zero；

Step 3, decomposes, to set empirical mode EEMD, each Intrinsic mode functions IMF obtained and carries out spectrum analysis, select Go out and comprise the Intrinsic mode functions IMF turning frequency_x(t)And IMF_y(t)；

Step 4, respectively to Intrinsic mode functions IMF_x(t)And IMF_y(t)Direct orthogonalization method is used to calculate instantaneous frequency IF_x(t)And IF_y(t), calculate process as follows:

Intrinsic mode functions IMF_x(t)Resolve into amplitude-modulated signal A (t) and FM signalThat is:

Order:Orthogonal function is defined from the experience frequency modulation component of signal, as follows:

Therefore, the instantaneous phase of signalCalculated by following formula and obtain:

By to instantaneous phaseDerivation, obtains instantaneous frequency IF_x(t)For:

It is calculated instantaneous frequency IF according still further to said process_y(t)；

Step 5, extracts instantaneous frequency IF_x(t)And IF_y(t)Frequency, amplitude, phase information, synthesize two-dimension holographic spectrum, so Effectively identifying rotor operation state afterwards, detailed process is as follows:

First, with interpolation respectively to X, instantaneous frequency IF of Y-direction_x(t)And IF_y(t)Carry out frequency, amplitude, phase place school Just；Then, amplitude A after correction is extracted_x、A_yWith phase place Merge in units of frequency order, obtain ellipse Equation is:

In formula, α n is the displacement transducer angle with horizontal direction of X-direction；

Change the A at different frequency order_x、A_y、Just obtaining corresponding ellipse, combination obtains two-dimension holographic spectrum, As shown in Figure 2 and Figure 3, Fig. 2 is analysis of vibration signal result under rotor normal operating condition, and Fig. 3 is that rotor impact and rub fault is shaken Dynamic signal analysis result；

Finally, calculate eccentricity e at each frequency order, major axis a and the oval information of short axle b, according to oval information Effectively identifying rotor operation state, calculation expression is as follows:

$e=\sqrt{a^2-b^2}/a,$

$a=(\sqrt{p+q}+\sqrt{p-q})/2,$

$b=(\sqrt{p+q}-\sqrt{p-q})/2,$

Wherein,

The effect of the present embodiment: the holographic diagnostics method based on frequency modulated information reconstruct using the present invention to propose is carried out point Analysis, comparison diagram 2 and Fig. 3, when rotor is properly functioning, each frequency order ellipse presents one group of the least ellipse of random amplitude, As shown in Figure 2；During rotor generation impact and rub fault, turn the oval substantially change of the most oval and high frequency multiplication big, in one group of eccentricity The biggest ellipse, as it is shown on figure 3, this feature can effectively identify rotor impact and rub fault.

Embodiment 2, carries out experimental verification by Bentley RK4 rotor experiment table, and Bentley RK4 rotor experiment table structure is shown Being intended to as shown in Figure 4, root diameter is 10mm, long 550mm, and front and back supporting span is 325mm, in terms of drive end, and this rotor system System rotates in the counterclockwise direction, at rotor across interior axially symmetric layout two for adding the balancing frame of balance mass, each Dish weight is 0.8kg, and on balancing frame, radius is that at 30mm, the most angularly 16 holes of distribution, for adding balance mass, are tested Touched the spacing model rotor impact and rub fault of frotton and rotor by adjustment, experiment rotating speed is set to 1700rpm.

A kind of rotating machinery holographic diagnostics method based on frequency modulated information reconstruct, comprises the following steps:

Step one, at the displacement transducer that rotor oscillation testing section two angles of installation are 90 °, it is stipulated that rotor rotates through In journey, the displacement transducer of first process is X-direction, and clockwise another displacement transducer after half-twist is Y-direction, two The mounting condition of individual displacement transducer, physical characteristic keep consistent, gather displacement vibration signal x (t) and the y (t) of rotor；

Step 2, carries out gathering empirical mode EEMD to X, displacement vibration signal x (t) of Y-direction and y (t) respectively and decomposes, Obtaining multiple Intrinsic mode functions IMF, as shown in Figure 5, Figure 6, Intrinsic mode functions IMF must is fulfilled for following two condition:

Extreme point number on the most whole signal and zero crossing number is equal or at most difference one；

B. in whole time series, all Local modulus maximas the coenvelope that determines and by all local minizing points Determined by the average of lower envelope be zero；

Step 3, decomposes, to set empirical mode EEMD, each Intrinsic mode functions IMF obtained and carries out spectrum analysis, select Go out and comprise the Intrinsic mode functions IMF turning frequency_x(t)And IMF_y(t), as shown in Figure 7, Figure 8, X-direction displacement vibration signal x (t) is passed through After set empirical mode EEMD is decomposed, the Intrinsic mode functions IMF3 obtained comprises a turn frequency, its time domain beamformer and amplitude spectrum As shown in Figure 7；Y-direction displacement vibration signal y (t) through set empirical mode EEMD decompose after, the Intrinsic mode functions obtained IMF2 comprises a turn frequency, and its time domain beamformer and amplitude spectrum are as shown in Figure 8；

Step 4, respectively to Intrinsic mode functions IMF_x(t)And IMF_y(t)Direct orthogonalization method is used to calculate instantaneous frequency IF_x(t)And IF_y(t), as shown in Figure 9, Figure 10, calculate process as follows:

Intrinsic mode functions IMF_x(t)Resolve into amplitude-modulated signal A (t) and FM signalThat is:

Order:Orthogonal function is defined from the experience frequency modulation component of signal, as follows:

Therefore, the instantaneous phase of signalCalculated by following formula and obtain:

By to instantaneous phaseDerivation, obtains instantaneous frequency IF_x(t)For:

It is calculated instantaneous frequency IF according still further to said process_y(t)；

Step 5, extracts instantaneous frequency IF_x(t)And IF_y(t)Frequency, amplitude, phase information, synthesize two-dimension holographic spectrum, so Effectively identifying rotor operation state afterwards, detailed process is as follows:

First, with interpolation respectively to X, instantaneous frequency IF of Y-direction_x(t)And IF_y(t)Carry out frequency, amplitude, phase place school Just；Then, amplitude A after correction is extracted_x、A_yWith phase place Merge in units of frequency order, obtain ellipse Equation is:

In formula, α n is the displacement transducer angle with horizontal direction of X-direction；

Change the A at different frequency order_x、A_y、Just obtaining corresponding ellipse, combination obtains two-dimension holographic spectrum, As shown in Figure 11, Figure 12, Figure 11 is analysis of vibration signal result under rotor normal operating condition, and Figure 12 is the event of rotor impact and rub Barrier analysis of vibration signal result；

Finally, calculate eccentricity e at each frequency order, major axis a and the oval information of short axle b, according to oval information Effectively identifying rotor operation state, calculation expression is as follows:

$e=\sqrt{a^2-b^2}/a,$

$a=(\sqrt{p+q}+\sqrt{p-q})/2,$

$b=(\sqrt{p+q}-\sqrt{p-q})/2,$

Wherein,

The effect of the present embodiment: the holographic diagnostics method based on frequency modulated information reconstruct using the present invention to propose is carried out point Analysis, contrasts Figure 11 and Figure 12, and when rotor is properly functioning, each frequency order ellipse presents the least ellipse of one group of random amplitude Circle, as shown in figure 11；During rotor generation impact and rub fault, turn that frequency is the most oval and high frequency multiplication is oval substantially becomes big, in one group from The ellipse that heart rate is the biggest, as shown in figure 12, this feature can effectively identify rotor impact and rub fault.

Above content is to combine concrete preferred implementation further description made for the present invention, it is impossible to assert The detailed description of the invention of the present invention is only limitted to this, for general technical staff of the technical field of the invention, is not taking off On the premise of present inventive concept, it is also possible to make some simple deduction or replace, all should be considered as belonging to the present invention by institute The claims submitted to determine scope of patent protection.

Claims (1)

1. a rotating machinery holographic diagnostics method based on frequency modulated information reconstruct, it is characterised in that comprise the following steps:

Step one, at the displacement transducer that rotor oscillation testing section two angles of installation are 90 °, it is stipulated that in rotor rotation process The displacement transducer of first process is X-direction, and clockwise another displacement transducer after half-twist is Y-direction, two positions The mounting condition of displacement sensor, physical characteristic keep consistent, gather displacement vibration signal x (t) and the y (t) of rotor；

Step 2, carries out gathering empirical mode EEMD to X, displacement vibration signal x (t) of Y-direction and y (t) respectively and decomposes, obtain Multiple Intrinsic mode functions IMF, Intrinsic mode functions IMF must are fulfilled for following two condition:

Extreme point number on the most whole signal and zero crossing number is equal or at most difference one；

B. in whole time series, all Local modulus maximas the coenvelope that determines and by institute of all local minizing points really The average of fixed lower envelope is zero；

Step 3, decomposes, to set empirical mode EEMD, each Intrinsic mode functions IMF obtained and carries out spectrum analysis, pick out bag Containing the Intrinsic mode functions IMF turning frequency_x(t)And IMF_y(t)；

Step 4, respectively to Intrinsic mode functions IMF_x(t)And IMF_y(t)Direct orthogonalization method is used to calculate instantaneous frequency IF_x(t) And IF_y(t), calculate process as follows:

Intrinsic mode functions IMF_x(t)Resolve into amplitude-modulated signal A (t) and FM signalThat is:

Order:Orthogonal function is defined from the experience frequency modulation component of signal, as follows:

Therefore, the instantaneous phase of signalCalculated by following formula and obtain:

By to instantaneous phaseDerivation, obtains instantaneous frequency IF_x(t)For:

It is calculated instantaneous frequency IF according still further to said process_y(t)；

Step 5, extracts instantaneous frequency IF_x(t)And IF_y(t)Frequency, amplitude, phase information, synthesize two-dimension holographic spectrum, the most right Rotor operation state effectively identifies, detailed process is as follows:

First, with interpolation respectively to X, instantaneous frequency IF of Y-direction_x(t)And IF_y(t)Carry out frequency, amplitude, phasing；So After, extract amplitude A after correction_x、A_yWith phase place Merge in units of frequency order, obtain elliptic equation For:

In formula, α n is the displacement transducer angle with horizontal direction of X-direction；

Change the A at different frequency order_x、A_y、Just obtaining corresponding ellipse, combination obtains two-dimension holographic spectrum；

Finally, calculate eccentricity e at each frequency order, major axis a and the oval information of short axle b, according to oval information to turning Sub-running status effectively identifies, calculation expression is as follows:

$e=\sqrt{a^2-b^2}/a,$

$a=(\sqrt{p+q}+\sqrt{p-q})/2,$

$b=(\sqrt{p+q}-\sqrt{p-q})/2,$

Wherein,