CN107525671A - A kind of wind-powered electricity generation driving-chain combined failure character separation and discrimination method - Google Patents
A kind of wind-powered electricity generation driving-chain combined failure character separation and discrimination method Download PDFInfo
- Publication number
- CN107525671A CN107525671A CN201710628322.4A CN201710628322A CN107525671A CN 107525671 A CN107525671 A CN 107525671A CN 201710628322 A CN201710628322 A CN 201710628322A CN 107525671 A CN107525671 A CN 107525671A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- signal
- mfrac
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Mobile Radio Communication Systems (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
A kind of wind-powered electricity generation driving-chain combined failure character separation and discrimination method, step are as follows:Vibration signal is obtained using acceleration transducer;Using multi-wavelet packets transform, will be decomposed in the vibration signal whole frequency range of collection;To arrange entropy as evaluation index, select satisfactory single supported signal to be reconstructed respectively, complete noise reduction and the combined failure separation of signal;Using energy operator demodulation method to reconstructing signal transacting, the identification of fault message is completed.
Description
Technical field
The present invention relates to a kind of wind-powered electricity generation driving-chain combined failure character separation method.
Background technology
With the fast development of wind generating technology, the function of wind power equipment becomes strong, volume becomes big, degree of intelligence uprises,
Corresponding operation and maintenance cost is also continuously increased, and is increased with Wind turbines active time, and fault rate dramatically increases, data statistics
Proving, failure is concentrated mainly on the parts such as wind-powered electricity generation electric control system, transmission system, blade, once breaking down can cause to stop
Machine, huge economic loss is brought, wherein, driving unit fault causes fault power time most long, it is therefore necessary to wind-powered electricity generation
The running status of transmission chain system is monitored and diagnosed.
In actual motion, wind power equipment is often due to a variety of causes causes various faults while occurred, it is also possible to a kind of
The operation conditions of the wind power equipment that changes of failure, and then cause associated components to break down.Gear-box is as wind-power electricity generation
The important transmission parts of equipment, nonserviceable lower operation when, the signal that collects is mostly the fault-signal modulated, this be by
In when a failure occurs it, fault-signal is mainly shown as:Frequency and its multiple-frequency modulation are turned by axle around the meshing frequency of gear;Axle
Hold when failing, the vibration signal collected can be by periodic temporary impact pulse signal modulation.Therefore, how accurately,
Quickly fault signature is extracted, is the key to the diagnosis of wind power equipment running status, energy operator demodulation is as a kind of letter
Single quick demodulation method, can effectively solve the problem that the modulation phenomenon of signal.But due to energy operator demodulation method, it is vulnerable to noise
Influence, and the actual vibration signal obtained of engineering is mostly complicated amplitude-modulation frequency-modulation signal, and contain strong ambient noise, directly
Energy operator demodulation is carried out, it is difficult to extract whole fault signatures.
Chinese invention patent CN103900816A " a kind of wind power generating set Method for Bearing Fault Diagnosis " uses wavelet packet
Method is decomposed to vibration signal, is carried out noise reduction to the high frequency coefficient after decomposition using soft-threshold, is then carried out signal reconstruction,
Carry out three layers of decomposition of wavelet packet again to the signal after reconstruct, obtain the energy of each frequency band of third layer, and feature is used as using the energy
Vector input BP neural network, realizes Wind turbines bearing failure diagnosis.The patent is realized to signal using method of wavelet packet
Noise reduction and failure whether accurate judgement, but can not really realize the separation and identification of the combined failure of modulation.Chinese invention
Patent CN102937522A " a kind of gear-box combined failure diagnosis and system " uses dual-tree complex wavelet transform method decomposed signal,
Multiple subbands are obtained, each subband signal is analyzed using energy operator method respectively, extracts fault signature, though the patent can identify again
Failure is closed, but because the subband signal after decomposition is without relevant evaluation index, each subband signal need to be handled successively, cause to count
The result that evaluation time is long and accuracy is not high.Document " is based on LMD self-adapting multi-dimensions morphology and Teager energy operator methods
Application in bearing failure diagnosis " proposes to be decomposed into multiple subband signals using LMD (part mean decomposition method), and utilizes
Multi-scale morphology filtering is handled each subband, is demodulated after processing using energy operator, obtains bearing fault characteristics frequency,
This method can effectively extract fault message, but equally exist to the subband of decomposition without relevant evaluation On Index, need to be right successively
Each subband signal is handled, and causes to calculate the problem of time is long and precision is not high.
The content of the invention
The purpose of the present invention is to overcome the disadvantage mentioned above of prior art, proposes a kind of wind-powered electricity generation driving-chain combined failure feature point
From with discrimination method.
Present invention introduces multi-wavelet packets algorithm, multi-wavelet packets decomposition method can by signal in whole frequency ranges it is fine
Change and decompose.The characteristics of present invention contains multiple time-frequency characteristics basic functions using multi-wavelet packets, pass through of basic function and fault signature
Match somebody with somebody, obtain comprehensive fault message, realize the identification and diagnosis of failure.The present invention is using multi-wavelet packets transform to containing compound event
The modulated signal of barrier is decomposed, and feature reconstruction is carried out as evaluation index, the qualified single supported signal of selection to arrange entropy,
The fault-signal reconstructed with energy operator demodulation method, demodulation analysis, demodulation frequency spectrum is obtained, extracts fault signature, realizes event
Hinder the separation and diagnosis of feature.
The present invention includes the feature using multi-wavelet packets method separation signal, the single supported signal of selection, reconstruction signal, energy and calculated
Subsolution is adjusted with identifying the steps such as fault signature, specific as follows:
1. using vibration acceleration sensor collection wind-powered electricity generation driving-chain vibration signal x (t), wherein t represents collection signal institute
The corresponding time;
2. the vibration signal x (t) that step 1 collects is decomposed using multi-wavelet packets method, it is specific as follows:
(1) multi-wavelet pretreatment vibration signal x (t);
Because the vibration signal x (t) is one-dimensional signal, and m ultiwavelet generally comprises multiple scaling functions and wavelet function, because
This carries out multi-wavelet pretreatment using repeated sampling method to vibration signal x (t), obtains two-dimentional multi-wavelet pretreatment signal x1
(t), t represents the time corresponding to collection signal.
(2) to two-dimentional multi-wavelet pretreatment signal x1(t) multi-wavelet packets decomposition is carried out;
M ultiwavelet decomposition formula is as follows:
In formula, sj,kRepresent the scale coefficient of jth layer decomposition, dj,kRepresent the wavelet coefficient of jth layer decomposition, Hn-2kFor signal
Carry out m ultiwavelet and decompose low-pass filter coefficients, Gn-2kM ultiwavelet is carried out for signal and decomposes high-pass filter coefficient, and j is small wavelength-division
Solve the number of plies, sj-1,nRepresent the scale coefficient of -1 layer of decomposition of jth, dj-1,nThe wavelet coefficient of -1 layer of decomposition of jth is represented, n is sampled point
Number, k are the sampling number k=0,1 ..., n-1 that m ultiwavelet decomposes jth layer coefficients;
Scaling function and wavelet function are decomposed simultaneously, wherein third layer decomposes to obtain 4 groups of scaling functions and 4 groups of small echo letters
Count, totally 16 single supported signals.
3. to arrange entropy as evaluation index, single supported signal of the acquisition of analytical procedure 2, chooses and is enriched containing characteristic information respectively
Single supported signal carry out the reconstruct of multi-wavelet packets list branch, obtain multi-wavelet packets list branch reconstruction signal;
(1) arrangement entropy is solved first, is comprised the following steps that:
1) phase space of reconstruction signal
It is assumed that discrete-time series x (t), t=1,2 ... and .., N }, phase space reconfiguration is carried out to it, obtained:
In formula:Y is the phase space of reconstruct
N is discrete time data length
M is Embedded dimensions
τ is time delay
K be reconstruct component number, k=N-m+1
X (j) is restructuring matrix jth row component
J be reconstruct phase space any row component, j=1,2 ... k
2) primary signal is reconfigured
In matrix row component Y (j,:) it is considered as a reconstruct component, share k reconstruct component, k=N-m+1, to Y
(j,:) in each element [x (j) x (j+ τ) ... x (j+ (m-1) τ)] according to rearranging from small to large, i1 i2...inTable
Show the index position for rearranging rear element column, there is unique group code sequence after rearrangement:S (l)={ i1
i2...id, in formula, l=1,2 ..., k, and k≤m!, the probability by the way that every kind of sequence is calculated is p1 p2…pk;
3) value of arrangement entropy is solved
The arrangement entropy H of Discrete Stochastic time series { x (t), t=1,2 ... .., N }pDefined by formula:J=1,2 ..., k, pjFor the probability of j-th of sequence;
4) normalization arrangement entropy
Data are normalized:HP=HP(m)/ln(m!).
(2) the single supported signal enriched containing characteristic information is chosen:
It is 0.79 to set the threshold value of arrangement entropy according to Chebyshev inequality, and arrangement entropy is more than 0.79 single supported signal conduct
Reconstruction signal, data of the present invention select single supported signal that two arrangement entropys are more than 0.79.
4. reconstruction signal, single supported signal is reconstructed using m ultiwavelet, obtains reconstruction signal ya(t), wherein, a represents reconstruct letter
Number number, t represents the time corresponding to reconstruction signal;
As the inverse transformation process decomposed, how small the abundant single supported signal of the characteristic information to selection carry out for m ultiwavelet reconstruct
Ripple bag reconstructs, and specific formula is as follows:
K=0,1 ..., N-1
In formula, sj,nRepresent scale coefficient, dj,nRepresent wavelet coefficient, Hk-2nChecked the number for letter and carry out m ultiwavelet reconstruction low pass
Ripple device coefficient, Gk-2nHigh-pass filter coefficient is reconstructed to carry out m ultiwavelet to signal, j is the wavelet reconstruction number of plies, and N is sampled point
Number, sj-1,nRepresent the scale coefficient of -1 layer of decomposition of jth, dj-1,nRepresent the wavelet coefficient of -1 layer of jth decomposition, n is sampling number, k
The sampling number k=0,1 of jth layer coefficients is decomposed for m ultiwavelet ..., n-1.
5. demodulation analysis, the reconstruction signal y obtained using energy operator method demodulation step 4a(t) solved corresponding to, obtaining
Adjust spectrum signature;
Energy operator demodulation method step is as follows:
(1) for reconstruction signal ya(t) its energy operator ψ, is definedCFor:
Wherein, t represents the time corresponding to reconstruction signal,For reconstruction signal ya(t) time t is asked single order,
Second-order differential obtains.
(2) instantaneous amplitude and instantaneous frequency of amplitude modulationfrequency modulation time signal are solved using energy operator:
Wherein, t represents the time corresponding to reconstruction signal, i.e., the time corresponding to wind-powered electricity generation driving-chain, and a (t) is instantaneous amplitude, wi
For instantaneous frequency.
6. demodulation spectra characteristic frequency and the fault characteristic frequency in wind-powered electricity generation driving-chain running that step 5 obtains are contrasted,
The running status of wind-powered electricity generation driving-chain is judged, wherein bearing fault frequency, gear distress frequency calculation formula are as follows:
(1) bearing fault frequency:
1) outer ring failure-frequency fW
2) inner ring failure-frequency fN
3) rolling element failure-frequency fG
4) retainer failure-frequency fB
In formula:d0For rolling element diameter;D is pitch diameter;α is the nominal contact angle of bearing;frFor the rotating speed frequency of bearing
Rate;Z is the rolling element number of bearing.
(2) gear distress frequency:
Meshing frequency fM
fM=f1Z1=f2Z2
In formula:, frFor the speed-frequency of gear;Z is number of gear teeth;N gear rotational speeds, unit r/min;f1For driving wheel
Speed-frequency;f2For secondary speed frequency;Z1For the number of teeth of driving wheel;Z2For the number of teeth of driven pulley.
Compared with prior art, the present invention is taking into full account wind-powered electricity generation drive mechanism complexity, multiple faults frequency modulation(PFM) easily occurs
On the premise of problem, wind-powered electricity generation driving-chain combined failure is diagnosed with reference to multi-wavelet packets algorithm and energy operator demodulating algorithm.Wherein,
Energy operator demodulation method can effectively solve the problem that the modulation phenomenon of signal as a kind of simple and quick demodulation method, but due to
Energy operator demodulation method, is vulnerable to influence of noise.Multi-wavelet packets transform is used as a kind of efficient, accurate decomposition method, relatively
It can meet the characteristics such as orthogonality, compact sup-port, symmetry, high-order vanishing moment simultaneously in wavelet packet analysis, can effectively solve the problem that energy
The problem of amount operator demodulation method is easily influenceed by ambient noise, and the concept of arrangement entropy is introduced, evaluation of the arrangement entropy as signal
Index quantification, the quick single supported signal for determining to include fault message, the selection for single branch reconstruction signal provide scientific basis, therefore
The present invention is applied in wind-powered electricity generation driving-chain combined failure character separation and identification process, you can is precisely separating combined failure, again may be used
Accurate identification fault signature.Also, the present invention reconstructs single supported signal using multi-wavelet packets, and the reconstruct of multi-wavelet packets list branch can either be real
Existing signal de-noising and can enough completes the character separation of combined failure.
Brief description of the drawings
Fig. 1 is the overview flow chart of the present invention;
Fig. 2 is multi-wavelet packets character separation flow chart;
Fig. 3 is that multi-wavelet packets list branch reconstructs flow chart;
Fig. 4 is the result figure individually identified to combined failure using energy operator demodulating algorithm, and Fig. 4 a are time-domain diagram, Fig. 4 b
To demodulate spectrogram;
Fig. 5 is the result figure to combined failure identification using the present invention, and Fig. 5 a are single time-domain diagram of branch reconstruction signal 1, Fig. 5 b
Spectrogram is demodulated to be corresponding, Fig. 5 c are single time-domain diagram of branch reconstruction signal 2, and Fig. 5 d are corresponding demodulation spectrogram.
Embodiment
The present invention is further illustrated below in conjunction with the drawings and specific embodiments.
As shown in figure 1, combined failure character separation of the present invention and discrimination method step are as follows:
1. using vibration acceleration sensor collection wind-powered electricity generation driving-chain vibration signal x (t), wherein t represents collection signal institute
The corresponding time;
2. being decomposed using multi-wavelet packets method to described vibration signal x (t), flow is as shown in Fig. 2 specifically such as
Under:
(1) multi-wavelet pretreatment is carried out to vibration signal x (t);
Because the signal x (t) is one-dimensional signal, and m ultiwavelet generally comprises multiple scaling functions and wavelet function, therefore adopts
Multi-wavelet pretreatment is carried out to vibration signal x (t) with repeated sampling method, obtains two-dimentional multi-wavelet pretreatment signal x1(t), t
Represent the time corresponding to collection signal;
(2) to x1(t) multi-wavelet packets decomposition is carried out;
M ultiwavelet decomposition formula is as follows:
K=0,1 ..., n-1
In formula, sj,kRepresent scale coefficient, dj,kRepresent wavelet coefficient, Hn-2kM ultiwavelet, which is carried out, for signal decomposes LPF
Device coefficient, Gn-2kM ultiwavelet is carried out for signal and decomposes high-pass filter coefficient, and j is the wavelet decomposition number of plies, and n is sampling number t generations
Time corresponding to table wind-powered electricity generation driving-chain vibration signal;In (implication that K please be illustrate) formula, sj,kRepresent the yardstick system of jth layer decomposition
Number, dj,kRepresent the wavelet coefficient of jth layer decomposition, Hn-2kM ultiwavelet, which is carried out, for signal decomposes low-pass filter coefficients, Gn-2kFor letter
Number carrying out m ultiwavelet decomposes high-pass filter coefficient, and j is the wavelet decomposition number of plies, sj-1,nThe scale coefficient of -1 layer of decomposition of jth is represented,
dj-1,nThe wavelet coefficient of -1 layer of decomposition of jth is represented, n is sampling number, and k is the sampling number k=that m ultiwavelet decomposes jth layer coefficients
0,1,…,n-1;
Scaling function and wavelet function are decomposed simultaneously, wherein third layer decomposes to obtain 4 groups of scaling functions and 4 groups of small echo letters
Count, totally 16 single supported signals.
3. to arrange entropy as evaluation index, single supported signal of the acquisition of analytical procedure 2, chooses qualified Dan Zhixin respectively
Number the reconstruct of multi-wavelet packets list branch is carried out, obtain single branch reconstruction signal;
(1) arrangement entropy is solved first, is comprised the following steps that:
1) phase space of reconstruction signal
It is assumed that discrete-time series x (t), t=1,2 ... and .., N }, phase space reconfiguration is carried out to it, obtained:
In formula:Y is the phase space of reconstruct
N is discrete time data length
M is Embedded dimensions
τ is time delay
K be reconstruct component number, k=N-m+1
X (j) is restructuring matrix jth row component
J be reconstruct phase space any row component, j=1,2 ... k
2) primary signal is reconfigured
In matrix row component Y (j,:) it is considered as a reconstruct component, share k reconstruct component, k=N-m+1, to Y
(j,:) in each element [x (j) x (j+ τ) ... x (j+ (m-1) τ)] according to rearranging from small to large, i1 i2...inTable
Show the index position for rearranging rear element column, there is unique group code sequence after rearrangement:S (l)={ i1
i2...id, in formula, l=1,2 ..., k, and k≤m!, the probability by the way that every kind of sequence is calculated is p1 p2…pk。
3) value of arrangement entropy is solved
The arrangement entropy H of Discrete Stochastic time series { x (t), t=1,2 ... .., N }pDefined by formula:J=1,2 ..., k, pjFor the probability of j-th of sequence;
4) normalization arrangement entropy
Data are normalized:HP=HP(m)/ln(m!).
(2) the single supported signal enriched containing characteristic information is chosen:
It is 0.79 to set the threshold value of arrangement entropy according to Chebyshev inequality, and arrangement entropy is more than 0.79 single supported signal conduct
Reconstruction signal, data of the present invention select single supported signal that two arrangement entropys are more than 0.79;
It is reconstruction signal to determine LLH1, LHH2.
4. the single supported signal of reconstruct.LLH1, LHH2 are determined as reconstruction signal as evaluation index to arrange entropy, with multi-wavelet packets
Reconfiguration principle for instruct, single supported signal is reconstructed, complete reconstruct after, carry out m ultiwavelet post processing, flow as shown in figure 3,
Specific formula for calculation is as follows.
K=0,1 ..., N-1
In formula, sj,nRepresent scale coefficient, dj,nRepresent wavelet coefficient, Hk-2nChecked the number for letter and carry out m ultiwavelet reconstruction low pass
Ripple device coefficient, Gk-2nHigh-pass filter coefficient is reconstructed to carry out m ultiwavelet to signal, j is the wavelet reconstruction number of plies, and N is sampled point
Number, sj-1,nRepresent the scale coefficient of -1 layer of decomposition of jth, dj-1,nRepresent the wavelet coefficient of -1 layer of jth decomposition, n is sampling number, k
The sampling number k=0,1 of jth layer coefficients is decomposed for m ultiwavelet ..., n-1.
5. demodulation analysis, the reconstruction signal y obtained using energy operator demodulation method analytical procedure 4a(t), corresponded to
Demodulation spectrum signature;
Energy operator demodulation method step is as follows:
(1) for reconstruction signal ya(t) its energy operator ψ, is definedCFor:
Wherein, t represents the time corresponding to wind-powered electricity generation driving-chain vibration signal,For time signal yaPair (t) when
Between t ask single order, second-order differential to obtain.
(2) instantaneous amplitude and instantaneous frequency of amplitude modulationfrequency modulation time signal are solved using energy operator:
Wherein, t represents the time corresponding to wind-powered electricity generation driving-chain vibration signal, and a (t) is instantaneous amplitude, wiFor instantaneous frequency.
6. judge failure
Contrast step 5 obtain demodulation spectra frequecy characteristic with wind-powered electricity generation driving-chain running fault characteristic frequency therefore
Hinder frequency, judge the running status of Wind turbines, the present invention is recognized with wind-powered electricity generation driving-chain bearing fault to its running status:
Bearing fault frequency:
1) outer ring failure-frequency fW
2) inner ring failure-frequency fN
3) rolling element failure-frequency fG
4) retainer failure-frequency fB
In formula:d0For rolling element diameter;D is pitch diameter;α is the nominal contact angle of bearing;frFor the rotating speed frequency of bearing
Rate;Z is the rolling element number of bearing.
The Rolling Bearing Fault Character frequency of table 1
Fig. 4 b are that Fig. 4 a are obtained using the processing of energy operator demodulation method, the failure-frequency run according to known device:It is interior
Circle failure-frequency is 121.9Hz, is determined as bearing inner race failure.Because energy operator demodulates self-defect, using present invention side
Method handles wind-powered electricity generation driving-chain vibration signal, as a result as shown in Fig. 5 a, Fig. 5 b, Fig. 5 c, Fig. 5 d:Outer ring failure-frequency 76.88Hz, says
There is inside and outside circle combined failure in the bright bearing, it was demonstrated that the present invention can realize the separation and identification of combined failure simultaneously.
Claims (6)
1. a kind of wind-powered electricity generation driving-chain combined failure character separation and discrimination method, it is characterised in that:How small described method use is
Ripple bag method separates the feature of signal, selects single supported signal, reconstruction signal, and energy operator demodulation and identification fault signature,
Comprise the following steps that:
(1) using vibration acceleration sensor collection wind-powered electricity generation driving-chain vibration signal x (t), wherein t is represented corresponding to collection signal
Time;
(2) decomposed in vibration signal x (t) whole frequency ranges collected using multi-wavelet packets method to step (1);
(3) to arrange entropy as evaluation index, single supported signal of analytical procedure (2) acquisition, the list enriched containing characteristic information is chosen
Supported signal carries out the reconstruct of multi-wavelet packets list branch, obtains multi-wavelet packets list branch reconstruction signal;
(4) the selected single supported signal of reconstruct, single supported signal is reconstructed using multi-wavelet packets method respectively, obtains reconstruction signal ya(t),
Wherein a=1,2 ..., a represent the number of reconstruction signal, and t represents the time corresponding to reconstruction signal;
(5) demodulation analysis, the reconstruction signal y obtained using energy operator method demodulation step (4)a(t) demodulated corresponding to, obtaining
Spectrum signature;
(6) demodulation spectra characteristic frequency and the fault characteristic frequency in wind-powered electricity generation driving-chain running that step (5) obtains are contrasted, it is right
The running status of wind-powered electricity generation driving-chain judges.
2. wind-powered electricity generation driving-chain combined failure character separation according to claim 1 and discrimination method, it is characterised in that:It is described
The method that the vibration signal x (t) that step (2) is collected using multi-wavelet packets method to step (1) is decomposed is as follows:
(1) multi-wavelet pretreatment vibration signal x (t);
Multi-wavelet pretreatment is carried out to vibration signal x (t) using repeated sampling method, obtains two-dimentional multi-wavelet pretreatment signal x1
(t), t represents the time corresponding to collection signal;
(2) to two-dimentional multi-wavelet pretreatment signal x1(t) multi-wavelet packets decomposition is carried out;
M ultiwavelet decomposition formula is as follows:
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>n</mi>
</munder>
<msub>
<mi>H</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>s</mi>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>k</mi>
</munder>
<msub>
<mi>G</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
</mrow>
In formula, sj,kRepresent the scale coefficient of jth layer decomposition, dj,kRepresent the wavelet coefficient of jth layer decomposition, Hn-2kCarried out for signal
M ultiwavelet decomposes low-pass filter coefficients, Gn-2kM ultiwavelet is carried out for signal and decomposes high-pass filter coefficient, and j is wavelet decomposition layer
Number, sj-1,nRepresent the scale coefficient of -1 layer of decomposition of jth, dj-1,nRepresent the wavelet coefficient of -1 layer of jth decomposition, n is sampling number, k
The sampling number k=0,1 of jth layer coefficients is decomposed for m ultiwavelet ..., n-1;
Scaling function and wavelet function are decomposed simultaneously, wherein third layer decomposes to obtain 4 groups of scaling functions and 4 groups of wavelet functions,
Totally 16 single supported signals.
3. wind-powered electricity generation driving-chain combined failure character separation according to claim 1 and discrimination method, it is characterised in that:It is described
The step of (3) to arrange entropy as evaluation index, analytical procedure (2) obtain single supported signal, choose containing characteristic information enrich
Single supported signal carries out the reconstruct of multi-wavelet packets list branch, and the method for obtaining multi-wavelet packets list branch reconstruction signal is as follows:
(1) arrangement entropy is solved first:
1) phase space of reconstruction signal
It is assumed that discrete-time series x (t), t=1,2 ... and .., N }, phase space reconfiguration is carried out to it, obtained:
<mrow>
<mi>Y</mi>
<mo>=</mo>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>+</mo>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>k</mi>
</mrow>
In formula:Y is the phase space of reconstruct, and N is discrete time data length, and m is Embedded dimensions, and τ is time delay, and k is reconstruct
The number of component, k=N-m+1, x (j) are restructuring matrix jth row component, and j is the phase space any row component of reconstruct, j=1,
2,...k;
2) primary signal is reconfigured
In matrix row component Y (j,:) it is considered as a reconstruct component, share k reconstruct component, k=N-m+1, to Y
(j,:) in each element [x (j) x (j+ τ) ... x (j+ (m-1) τ)] according to rearranging from small to large, i1 i2 ...
inExpression rearranges the index position of rear element column, there is unique group code sequence after rearrangement:S (l)={ i1 i2
... id, in formula, l=1,2 ..., k, and k≤m!, the probability by the way that every kind of sequence is calculated is p1 p2 … pk;
3) value of arrangement entropy is solved
The arrangement entropy H of Discrete Stochastic time series { x (t), t=1,2 ... .., N }pDefined by formula:
J=1,2 ..., k, pjFor the probability of j-th of sequence;
4) normalization arrangement entropy
Data are normalized:HP=HP(m)/ln(m!);
(2) the single supported signal enriched containing characteristic information is chosen:
It is 0.79 to set the threshold value of arrangement entropy according to Chebyshev inequality, and single supported signal of the arrangement entropy more than 0.79 is as reconstruct
Signal, data of the present invention select single supported signal that two arrangement entropys are more than 0.79.
4. wind-powered electricity generation driving-chain combined failure character separation according to claim 1 and discrimination method, it is characterised in that:It is described
The step of the selected signal of (4) reconstruct, method is as follows:
The single supported signal enriched to the characteristic information that step (3) is chosen carries out multi-wavelet packets reconstruct, and specific formula is as follows:
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>n</mi>
</munder>
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>s</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>n</mi>
</munder>
<msub>
<mi>G</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>d</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
K=0,1 ..., N-1;
In formula, sj,nRepresent scale coefficient, dj,nRepresent wavelet coefficient, Hk-2nChecked the number for letter and carry out m ultiwavelet reconstruction low pass filter
Coefficient, Gk-2nHigh-pass filter coefficient is reconstructed to carry out m ultiwavelet to signal, j is the wavelet reconstruction number of plies, and N is sampling number,
sj-1,nRepresent the scale coefficient of -1 layer of decomposition of jth, dj-1,nThe wavelet coefficient of -1 layer of decomposition of jth is represented, n is sampling number, and k is
M ultiwavelet decomposes the sampling number k=0,1 of jth layer coefficients ..., n-1.
5. wind-powered electricity generation driving-chain combined failure character separation according to claim 1 and discrimination method, it is characterised in that:It is described
The step of (5) use the obtained reconstruction signal y of energy operator method demodulation step (4)a(t) demodulation spectrum signature corresponding to, obtaining
Method it is as follows:
Energy operator demodulation method step is as follows:
(1) for reconstruction signal ya(t) its energy operator ψ, is definedCFor:
<mrow>
<msub>
<mi>&psi;</mi>
<mi>C</mi>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>y</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>dy</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<msub>
<mi>y</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>dt</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Wherein, t represents the time corresponding to reconstruction signal,For reconstruction signal ya(t) single order, second order are asked to time t
Differential obtains;
(2) instantaneous amplitude and instantaneous frequency of amplitude modulationfrequency modulation time signal are solved using energy operator:
<mrow>
<mo>|</mo>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>a</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>a</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</msqrt>
</mfrac>
<mo>,</mo>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>a</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>a</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mfrac>
</msqrt>
</mrow>
Wherein, t represents the time corresponding to reconstruction signal, i.e., the time corresponding to wind-powered electricity generation driving-chain, and a (t) is instantaneous amplitude, wiFor wink
When frequency.
6. wind-powered electricity generation driving-chain combined failure character separation according to claim 1 and discrimination method, it is characterised in that:It is described
The step of (6) contrast step (5) demodulation spectra characteristic frequency and wind-powered electricity generation driving-chain running for obtaining in fault characteristic frequency,
Judge the running status of Wind turbines, wherein bearing fault frequency, gear distress frequency calculation formula is as follows:
(1) bearing fault frequency:
1) outer ring failure-frequency fW
<mrow>
<msub>
<mi>f</mi>
<mi>W</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mrow>
<mi>D</mi>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mi>Z</mi>
</mrow>
2) inner ring failure-frequency fN
<mrow>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mrow>
<mi>D</mi>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mi>Z</mi>
</mrow>
3) rolling element failure-frequency fG
<mrow>
<msub>
<mi>f</mi>
<mi>G</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>D</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mrow>
<mi>D</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
</mrow>
4) retainer failure-frequency fB
<mrow>
<msub>
<mi>f</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
</mrow>
<mi>D</mi>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
</mrow>
In formula:d0For rolling element diameter;D is pitch diameter;α is the nominal contact angle of bearing;frFor the speed-frequency of bearing;Z is
The rolling element number of bearing;
(2) gear distress frequency:
Meshing frequency fM
<mrow>
<msub>
<mi>f</mi>
<mi>M</mi>
</msub>
<mo>=</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mi>Z</mi>
<mo>=</mo>
<mfrac>
<mi>n</mi>
<mn>60</mn>
</mfrac>
<mi>Z</mi>
</mrow>
fM=f1Z1=f2Z2
In formula:, frFor the speed-frequency of gear;Z is number of gear teeth;N gear rotational speeds, unit r/min;f1For driving wheel rotating speed
Frequency;f2For secondary speed frequency;Z1For the number of teeth of driving wheel;Z2For the number of teeth of driven pulley.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710628322.4A CN107525671B (en) | 2017-07-28 | 2017-07-28 | Method for separating and identifying compound fault characteristics of transmission chain of wind turbine generator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710628322.4A CN107525671B (en) | 2017-07-28 | 2017-07-28 | Method for separating and identifying compound fault characteristics of transmission chain of wind turbine generator |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107525671A true CN107525671A (en) | 2017-12-29 |
CN107525671B CN107525671B (en) | 2020-12-18 |
Family
ID=60680279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710628322.4A Active CN107525671B (en) | 2017-07-28 | 2017-07-28 | Method for separating and identifying compound fault characteristics of transmission chain of wind turbine generator |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107525671B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110160778A (en) * | 2019-05-24 | 2019-08-23 | 武汉工程大学 | Gearbox fault state identification method based on sequential hypothesis testing |
CN110470462A (en) * | 2019-08-22 | 2019-11-19 | 苏州旋械感知信息科技有限公司 | One kind being based on C0The reconstructing method of the mechanical system fault features of complexity |
CN110595751A (en) * | 2019-09-19 | 2019-12-20 | 华东理工大学 | Early fault characteristic wavelet reconstruction method guided by Gini index and application thereof |
CN112146142A (en) * | 2019-06-27 | 2020-12-29 | 宁波方太厨具有限公司 | Method for recognizing sound fault of range hood |
CN112629850A (en) * | 2020-12-06 | 2021-04-09 | 北京工业大学 | Fault diagnosis method for variable-speed planetary gearbox |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6993377B2 (en) * | 2002-02-22 | 2006-01-31 | The Board Of Trustees Of The University Of Arkansas | Method for diagnosing heart disease, predicting sudden death, and analyzing treatment response using multifractal analysis |
WO2010006265A2 (en) * | 2008-07-10 | 2010-01-14 | Texas Heart Institute | Method and system for temperature analysis to provide an early marker of congestive heart failure progress that precedes a patient's symptoms |
CN102082754A (en) * | 2009-11-26 | 2011-06-01 | 中兴通讯股份有限公司 | OFDM (Orthogonal Frequency Division Multiplexing) channel estimation method and device |
CN102570979A (en) * | 2011-12-20 | 2012-07-11 | 重庆大学 | Iterative Teager energy operator demodulation method and system |
CN102937522A (en) * | 2012-08-30 | 2013-02-20 | 桂林电子科技大学 | Composite fault diagnosis method and system of gear case |
CN103076177A (en) * | 2013-01-16 | 2013-05-01 | 昆明理工大学 | Rolling bearing fault detection method based on vibration detection |
CN103499437A (en) * | 2013-09-11 | 2014-01-08 | 西安交通大学 | Rotary machine fault detection method of dual-tree complex wavelet transformation with adjustable quality factors |
WO2014096702A1 (en) * | 2012-12-19 | 2014-06-26 | Electricite De France | Locating of one or more defects in an electrochemical assembly |
CN103900816A (en) * | 2014-04-14 | 2014-07-02 | 上海电机学院 | Method for diagnosing bearing breakdown of wind generating set |
CN104374575A (en) * | 2014-11-25 | 2015-02-25 | 沈阳化工大学 | Wind turbine main bearing fault diagnosis method based on blind source separation |
CN104459388A (en) * | 2014-11-26 | 2015-03-25 | 国家电网公司 | Permanent magnetic direct-drive wind power generation system integrated fault diagnosis method |
CN105021277A (en) * | 2015-08-12 | 2015-11-04 | 黑龙江大学 | Wavelet-packet-correlation-dimension-combination-based vibration signal feature extraction method of high-voltage circuit breaker |
CN105758644A (en) * | 2016-05-16 | 2016-07-13 | 上海电力学院 | Rolling bearing fault diagnosis method based on variation mode decomposition and permutation entropy |
CN105784353A (en) * | 2016-03-25 | 2016-07-20 | 上海电机学院 | Fault diagnosis method for gear case of aerogenerator |
CN106228979A (en) * | 2016-08-16 | 2016-12-14 | 重庆大学 | A kind of abnormal sound in public places feature extraction and recognition methods |
CN105258940B (en) * | 2015-11-30 | 2018-02-09 | 上海无线电设备研究所 | The standardization m ultiwavelet that mechanical breakdown is quantitatively extracted and multi-wavelet packets transform method |
-
2017
- 2017-07-28 CN CN201710628322.4A patent/CN107525671B/en active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6993377B2 (en) * | 2002-02-22 | 2006-01-31 | The Board Of Trustees Of The University Of Arkansas | Method for diagnosing heart disease, predicting sudden death, and analyzing treatment response using multifractal analysis |
WO2010006265A2 (en) * | 2008-07-10 | 2010-01-14 | Texas Heart Institute | Method and system for temperature analysis to provide an early marker of congestive heart failure progress that precedes a patient's symptoms |
CN102082754A (en) * | 2009-11-26 | 2011-06-01 | 中兴通讯股份有限公司 | OFDM (Orthogonal Frequency Division Multiplexing) channel estimation method and device |
CN102570979A (en) * | 2011-12-20 | 2012-07-11 | 重庆大学 | Iterative Teager energy operator demodulation method and system |
CN102937522B (en) * | 2012-08-30 | 2015-01-21 | 桂林电子科技大学 | Composite fault diagnosis method and system of gear case |
CN102937522A (en) * | 2012-08-30 | 2013-02-20 | 桂林电子科技大学 | Composite fault diagnosis method and system of gear case |
WO2014096702A1 (en) * | 2012-12-19 | 2014-06-26 | Electricite De France | Locating of one or more defects in an electrochemical assembly |
CN103076177A (en) * | 2013-01-16 | 2013-05-01 | 昆明理工大学 | Rolling bearing fault detection method based on vibration detection |
CN103499437A (en) * | 2013-09-11 | 2014-01-08 | 西安交通大学 | Rotary machine fault detection method of dual-tree complex wavelet transformation with adjustable quality factors |
CN103900816A (en) * | 2014-04-14 | 2014-07-02 | 上海电机学院 | Method for diagnosing bearing breakdown of wind generating set |
CN104374575A (en) * | 2014-11-25 | 2015-02-25 | 沈阳化工大学 | Wind turbine main bearing fault diagnosis method based on blind source separation |
CN104459388A (en) * | 2014-11-26 | 2015-03-25 | 国家电网公司 | Permanent magnetic direct-drive wind power generation system integrated fault diagnosis method |
CN105021277A (en) * | 2015-08-12 | 2015-11-04 | 黑龙江大学 | Wavelet-packet-correlation-dimension-combination-based vibration signal feature extraction method of high-voltage circuit breaker |
CN105258940B (en) * | 2015-11-30 | 2018-02-09 | 上海无线电设备研究所 | The standardization m ultiwavelet that mechanical breakdown is quantitatively extracted and multi-wavelet packets transform method |
CN105784353A (en) * | 2016-03-25 | 2016-07-20 | 上海电机学院 | Fault diagnosis method for gear case of aerogenerator |
CN105758644A (en) * | 2016-05-16 | 2016-07-13 | 上海电力学院 | Rolling bearing fault diagnosis method based on variation mode decomposition and permutation entropy |
CN106228979A (en) * | 2016-08-16 | 2016-12-14 | 重庆大学 | A kind of abnormal sound in public places feature extraction and recognition methods |
Non-Patent Citations (2)
Title |
---|
刘朦月: ""基于振动信号的电机轴承故障诊断方法研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
杨志伊等: "《设备状态监测与故障诊断》", 30 June 2006 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110160778A (en) * | 2019-05-24 | 2019-08-23 | 武汉工程大学 | Gearbox fault state identification method based on sequential hypothesis testing |
CN112146142A (en) * | 2019-06-27 | 2020-12-29 | 宁波方太厨具有限公司 | Method for recognizing sound fault of range hood |
CN112146142B (en) * | 2019-06-27 | 2022-01-25 | 宁波方太厨具有限公司 | Method for recognizing sound fault of range hood |
CN110470462A (en) * | 2019-08-22 | 2019-11-19 | 苏州旋械感知信息科技有限公司 | One kind being based on C0The reconstructing method of the mechanical system fault features of complexity |
CN110470462B (en) * | 2019-08-22 | 2021-09-28 | 苏州旋械感知信息科技有限公司 | Based on C0Reconstruction method of early fault characteristics of mechanical system with complexity |
CN110595751A (en) * | 2019-09-19 | 2019-12-20 | 华东理工大学 | Early fault characteristic wavelet reconstruction method guided by Gini index and application thereof |
CN112629850A (en) * | 2020-12-06 | 2021-04-09 | 北京工业大学 | Fault diagnosis method for variable-speed planetary gearbox |
CN112629850B (en) * | 2020-12-06 | 2022-04-22 | 北京工业大学 | Fault diagnosis method for variable-speed planetary gearbox |
Also Published As
Publication number | Publication date |
---|---|
CN107525671B (en) | 2020-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107525671A (en) | A kind of wind-powered electricity generation driving-chain combined failure character separation and discrimination method | |
CN102721545B (en) | Rolling bearing failure diagnostic method based on multi-characteristic parameter | |
CN103499437B (en) | The rotating machinery fault detection method of adjustable quality factor dual-tree complex wavelet transform | |
CN102937522B (en) | Composite fault diagnosis method and system of gear case | |
CN104655423B (en) | A kind of Fault Diagnosis of Roller Bearings based on time-frequency domain multi-dimensional vibration Fusion Features | |
CN103810394B (en) | Rotating equipment failures signal singular values decomposes the method for designing of noise reduction | |
CN111160167A (en) | Spindle fault classification and identification method based on S-transform deep convolutional neural network | |
CN105758644A (en) | Rolling bearing fault diagnosis method based on variation mode decomposition and permutation entropy | |
CN105760839A (en) | Bearing fault diagnosis method based on multi-feature manifold learning and support vector machine | |
CN104316323B (en) | Method for confirming optimal resonance frequency band based on period target | |
CN104483127A (en) | Method for extracting weak fault characteristic information of planetary gear | |
CN107702922A (en) | Based on LCD and the Fault Diagnosis of Roller Bearings for stacking autocoder | |
CN104330257B (en) | A kind of planetary transmission system method for diagnosing faults | |
CN107608936A (en) | A kind of epicyclic gearbox combined failure feature extracting method | |
CN110046476A (en) | The ternary two of rolling bearing fault is into the sparse diagnostic method of Fractal Wavelet | |
CN109187023A (en) | A kind of automobile current generator bearing method for diagnosing faults | |
CN109765052B (en) | Planetary gearbox early fault diagnosis method based on GOA-ASR | |
CN108388908A (en) | Rolling bearing impulse fault diagnostic method based on K-SVD and sliding window noise reduction | |
CN110320040A (en) | A kind of Rolling Bearing Fault Character extracting method based on IITD and AMCKD | |
CN104330258A (en) | Method for identifying grey relational degree of rolling bearing fault based on characteristic parameters | |
CN105928702A (en) | Variable working condition gear case bearing fault diagnosis method based on form component analysis | |
CN109100144B (en) | Automobile hub bearing fault feature extraction method based on optimal quality factor selection | |
CN110146282A (en) | Epicyclic gearbox method for diagnosing faults based on multiple domain storehouse denoising autocoding network | |
CN106525426A (en) | Weak signal enhancement detection method based on complementary stochastic resonance filter | |
CN106092564A (en) | The gear failure diagnosing method demodulated based on ESMD and energy operator |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |