CN103412057A - Weak oscillation attenuation signal parameter identification based on stochastic resonance and moving least square - Google Patents
Weak oscillation attenuation signal parameter identification based on stochastic resonance and moving least square Download PDFInfo
- Publication number
- CN103412057A CN103412057A CN2013103064369A CN201310306436A CN103412057A CN 103412057 A CN103412057 A CN 103412057A CN 2013103064369 A CN2013103064369 A CN 2013103064369A CN 201310306436 A CN201310306436 A CN 201310306436A CN 103412057 A CN103412057 A CN 103412057A
- Authority
- CN
- China
- Prior art keywords
- signal
- accidental resonance
- oscillatory extinction
- moving
- frequency
- 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
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention relates to weak oscillation attenuation signal parameter identification based on stochastic resonance and moving least square. Firstly, directed at the characteristics of an oscillation attenuation shock signal, information entropy is employed as an immune optimization measure to perform algorithm improvement on parameter adjustment stochastic resonance, thus solving the problem of difficult calculation of the signal to noise ratio quantification measure of an oscillation attenuation signal. Then, two-step stochastic resonance and moving least squares fitting methods are employed to acquire the failure frequency, the oscillation frequency and damping of the weak oscillation attenuation signal and identify the primary fault of a rolling bearing. The technical scheme involved in the invention is of important significance to operating state monitoring and maintenance of rolling bearings and other rotating machineries, and also can be used for mode recognition and other applications.
Description
Technical field
The invention belongs to the mechanical fault diagnosis field, be specifically related to the weak oscillatory extinction signal parameter identification based on accidental resonance and Moving Least Squares.
Background technology
The oscillatory extinction signal is common rotating machinery such as rolling bearing, gear distress characterization signal.The detection of faint oscillatory extinction signal during for initial failure, accidental resonance is domestic and international main solution about the weak signal feature extraction at present, yet the research of accidental resonance is confined to the cosine signal that time domain is unlimited basically to be detected, and the oscillatory extinction signal has the advantages that time domain is tightly propped up, traditional accidental resonance is applied to the problem that faint oscillatory extinction signal demand solves quantitative evaluation parameter adjustment accidental resonance effect.The signal to noise ratio (S/N ratio) of estimating as cosine signal etc. is because the characteristics that the frequency domain of oscillatory extinction signal does not tightly prop up are difficult to calculate.On the other hand, what accidental resonance obtained is generally the failure-frequency of oscillatory extinction signal, thus judgement expression power fault.In order to obtain the quantification information of primary fault and fault degree, need further to obtain damping, oscillation frequency and the amplitude of oscillatory extinction signal.Yet after accidental resonance, signal has become rectangular signal, non-linear amplification has in various degree also been arranged simultaneously.Therefore the inversion method that needs research oscillatory extinction signal.
Summary of the invention
In order to overcome the shortcoming of above-mentioned prior art, the object of the present invention is to provide the weak oscillatory extinction signal parameter identification based on accidental resonance and Moving Least Squares, inverting obtains failure-frequency, oscillatory extinction frequency and the damping of oscillatory extinction signal, the primary fault of identification rolling bearing, significant to monitoring running state and the maintenance of the rotating machineries such as rolling bearing.
In order to achieve the above object, the technical solution used in the present invention is:
Weak oscillatory extinction signal parameter identification based on accidental resonance and Moving Least Squares comprises the following steps:
The first step, carry out the accidental resonance processing to the vibration acceleration signal gathered, and obtains the failure-frequency of vibration damping signal,
The failure-frequency of first pre-estimation vibration signal, in estimation range, select the parameter adjustment accidental resonance that is applicable to large parameter signal to process the vibration acceleration signal gathered,
Wherein, in parameter adjustment accidental resonance process, for the characteristics of oscillatory extinction impact signal, adopt and estimate and immune optimization algorithm based on information entropy, obtain optimum accidental resonance result;
Second step, according to the result of first step accidental resonance, intercept original vibration acceleration signal, and carry out secondary accidental resonance acquisition oscillatory extinction frequency,
According to the result of first step accidental resonance, obtain positional information and cycle information that each oscillatory extinction signal occurs, the vibration acceleration signal of acquired original is carried out to the windowing intercept operation,
Signal to intercepting carries out estimating with the second parameter of immune optimization algorithm and adjusting the accidental resonance operation based on information entropy, obtains the oscillatory extinction frequency of oscillatory extinction signal;
The 3rd step, utilize oscillatory extinction signal oscillating frequency of fadings that second step obtains as parameter, adopts the Moving Least Squares method be with ginseng signal retrieve, acquisition oscillatory extinction signal damping,
According to the result of second step accidental resonance, obtain positional information and cycle information that each oscillatory extinction signal occurs,
Long and the initial position of the window that Moving Least Squares is set according to the positional information obtained and cycle information, then carry out the match of Moving Least Squares method, obtains damping and the amplitude of oscillatory extinction signal.
Advantage of the present invention is: select information entropy to estimate as the immune optimization of parameter adjustment accidental resonance, solved the oscillatory extinction Signal-to-Noise and estimated the problem that is difficult to calculate, obtain accidental resonance effect accurately.By twice accidental resonance, can extract failure-frequency, oscillatory extinction frequency, damping and the amplitude information of signal, can provide diagnosing primary fault and quantification to pass judgment on the foundation of fault degree to the tester.
The accompanying drawing explanation
Fig. 1 is process flow diagram of the present invention.
The vibration acceleration signal figure of Fig. 2 for gathering.
Fig. 3 is the figure as a result of parameter adjustment accidental resonance for the first time.
Fig. 4 is the vibration acceleration signal figure of intercepting.
Fig. 5 is the figure as a result of parameter adjustment accidental resonance for the first time.
Fig. 6 is rolling bearing peak index design sketch.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in detail.
With reference to accompanying drawing 1, the weak oscillatory extinction signal parameter identification based on accidental resonance and Moving Least Squares comprises the following steps:
The first step, carry out the accidental resonance processing to the vibration acceleration signal gathered as shown in Figure 2, obtains failure-frequency, positional information and the oscillatory extinction signal number of vibration damping signal,
Start, reading out data, the failure-frequency of first pre-estimation vibration signal, carry out parameter initialization, in estimation range, selects the parameter adjustment accidental resonance that is applicable to large parameter signal to process the vibration acceleration signal gathered,
Wherein, in parameter adjustment accidental resonance process, characteristics for the oscillatory extinction impact signal, employing is estimated and immune optimization algorithm based on information entropy, obtain optimum accidental resonance result, thereby obtain failure-frequency, positional information and oscillatory extinction signal number information, the T as shown in Figure 3 of vibration damping signal
0For the inaction interval of vibration damping signal, τ
0Positional information for the vibration damping signal;
Second step, according to the result of first step accidental resonance, intercept original vibration acceleration signal, and carry out secondary accidental resonance acquisition oscillatory extinction frequency,
According to the result of first step accidental resonance, obtain positional information and cycle information that each oscillatory extinction signal occurs, the vibration acceleration signal of acquired original is carried out to intercept signal, intercept signal as shown in Figure 4,
Estimate the scope of oscillation frequency, and in this scope, intercept signal is carried out estimating with the second parameter of immune optimization algorithm and adjusting the accidental resonance operation based on information entropy, obtain the oscillatory extinction frequency of oscillatory extinction signal, result as shown in Figure 5;
The 3rd step, utilize oscillatory extinction signal oscillating frequency of fadings that second step obtains as parameter, adopts the Moving Least Squares method be with ginseng signal retrieve, acquisition oscillatory extinction signal damping, fitting result as shown in Figure 6,
Window length and the initial position of Moving Least Squares are set according to failure-frequency, positional information and the oscillatory extinction signal number information of first step accidental resonance result acquisition, utilize the oscillatory extinction signal oscillating frequency of fadings that second step obtains, as Moving Least Squares method known parameters, intercept signal is carried out to match, obtain damping and the amplitude of oscillatory extinction signal.
Claims (1)
1. based on the weak oscillatory extinction signal parameter identification of accidental resonance and Moving Least Squares, it is characterized in that, comprise the following steps:
The first step, carry out the accidental resonance processing to the vibration acceleration signal gathered, and obtains the failure-frequency of vibration damping signal,
The failure-frequency of first pre-estimation vibration signal, in estimation range, select the parameter adjustment accidental resonance that is applicable to large parameter signal to process the vibration acceleration signal gathered,
Wherein, in parameter adjustment accidental resonance process, for the characteristics of oscillatory extinction impact signal, adopt and estimate and immune optimization algorithm based on information entropy, obtain optimum accidental resonance result;
Second step, according to the result of first step accidental resonance, intercept original vibration acceleration signal, and carry out secondary accidental resonance acquisition oscillatory extinction frequency,
According to the result of first step accidental resonance, obtain positional information and cycle information that each oscillatory extinction signal occurs, the vibration acceleration signal of acquired original is carried out to the windowing intercept operation,
Signal to intercepting carries out estimating with the second parameter of immune optimization algorithm and adjusting the accidental resonance operation based on information entropy, obtains the oscillatory extinction frequency of oscillatory extinction signal;
The 3rd step, utilize oscillatory extinction signal oscillating frequency of fadings that second step obtains as parameter, adopts the Moving Least Squares method be with ginseng signal retrieve, acquisition oscillatory extinction signal damping,
According to the result of second step accidental resonance, obtain positional information and cycle information that each oscillatory extinction signal occurs,
Long and the initial position of the window that Moving Least Squares is set according to the positional information obtained and cycle information, then carry out the match of Moving Least Squares method, obtains damping and the amplitude of oscillatory extinction signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310306436.9A CN103412057B (en) | 2013-07-19 | 2013-07-19 | Based on the weak oscillation damping signal parameter identification of accidental resonance and Moving Least Squares |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310306436.9A CN103412057B (en) | 2013-07-19 | 2013-07-19 | Based on the weak oscillation damping signal parameter identification of accidental resonance and Moving Least Squares |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103412057A true CN103412057A (en) | 2013-11-27 |
CN103412057B CN103412057B (en) | 2015-11-25 |
Family
ID=49605083
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310306436.9A Active CN103412057B (en) | 2013-07-19 | 2013-07-19 | Based on the weak oscillation damping signal parameter identification of accidental resonance and Moving Least Squares |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103412057B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104298878A (en) * | 2014-10-14 | 2015-01-21 | 天津工业大学 | Self-adaptive searching method for obtaining optimal parameters of stochastic resonance system |
CN105067105A (en) * | 2015-05-04 | 2015-11-18 | 西安交通大学 | Kinetic parameter identification method utilizing rotary machine start and stop transient signal feature |
CN107666328A (en) * | 2017-10-09 | 2018-02-06 | 中国电子科技集团公司第二十研究所 | Low signal-to-noise ratio satellite communication signals method of reseptance |
CN111220383A (en) * | 2019-11-14 | 2020-06-02 | 东华大学 | Transmission equipment weak fault detection method based on oscillation resonance theory |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101762746A (en) * | 2009-12-29 | 2010-06-30 | 西安交通大学 | Cyclostationary demodulation method of complex gear transmission chain current signal based on entropy optimization |
CN102175915A (en) * | 2010-12-17 | 2011-09-07 | 西安交通大学 | Current characteristic frequency extraction method of machine tool transmission system with immune random resonance |
CN102254223A (en) * | 2011-06-29 | 2011-11-23 | 西安交通大学 | Ant colony algorithm-based self-adapting stochastic resonance system parameter selecting method |
CN102735330A (en) * | 2012-06-15 | 2012-10-17 | 天津大学 | Self-adaption stochastic resonance weak signal detecting method based on particle swarm optimization algorithm |
CN102778356A (en) * | 2012-08-06 | 2012-11-14 | 西安交通大学 | Enhanced stochastic resonance system and diagnostic method for machinery fault based on system |
-
2013
- 2013-07-19 CN CN201310306436.9A patent/CN103412057B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101762746A (en) * | 2009-12-29 | 2010-06-30 | 西安交通大学 | Cyclostationary demodulation method of complex gear transmission chain current signal based on entropy optimization |
CN102175915A (en) * | 2010-12-17 | 2011-09-07 | 西安交通大学 | Current characteristic frequency extraction method of machine tool transmission system with immune random resonance |
CN102254223A (en) * | 2011-06-29 | 2011-11-23 | 西安交通大学 | Ant colony algorithm-based self-adapting stochastic resonance system parameter selecting method |
CN102735330A (en) * | 2012-06-15 | 2012-10-17 | 天津大学 | Self-adaption stochastic resonance weak signal detecting method based on particle swarm optimization algorithm |
CN102778356A (en) * | 2012-08-06 | 2012-11-14 | 西安交通大学 | Enhanced stochastic resonance system and diagnostic method for machinery fault based on system |
Non-Patent Citations (5)
Title |
---|
LI QIANG ET AL.: "Engineering signal processing based on adaptive step-changed stochastic resonance", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》 * |
MARK D. MCDONNELL ET AL.: "Optimal information transmission in nonlinear arrays through suprathreshold stochastic resonance", 《PHYSICS LETTERS A》 * |
***等: "基于近似熵测度的自适应随机共振研究", 《物理学报》 * |
杨宁等: "自适应随机共振在微弱信号检测中的应用", 《机械强度》 * |
王晶等: "采用遗传算法的自适应随机共振***弱信号检测方法研究", 《西安交通大学学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104298878A (en) * | 2014-10-14 | 2015-01-21 | 天津工业大学 | Self-adaptive searching method for obtaining optimal parameters of stochastic resonance system |
CN105067105A (en) * | 2015-05-04 | 2015-11-18 | 西安交通大学 | Kinetic parameter identification method utilizing rotary machine start and stop transient signal feature |
CN105067105B (en) * | 2015-05-04 | 2018-04-17 | 西安交通大学 | Utilize the kinetic parameter identification method of rotating machinery start and stop car transient signal feature |
CN107666328A (en) * | 2017-10-09 | 2018-02-06 | 中国电子科技集团公司第二十研究所 | Low signal-to-noise ratio satellite communication signals method of reseptance |
CN111220383A (en) * | 2019-11-14 | 2020-06-02 | 东华大学 | Transmission equipment weak fault detection method based on oscillation resonance theory |
Also Published As
Publication number | Publication date |
---|---|
CN103412057B (en) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Miao et al. | Application of an improved maximum correlated kurtosis deconvolution method for fault diagnosis of rolling element bearings | |
CN101561314B (en) | Method for testing stochastic resonance-chaotic weak signal | |
Ding et al. | Multi-objective iterative optimization algorithm based optimal wavelet filter selection for multi-fault diagnosis of rolling element bearings | |
CN105784366A (en) | Wind turbine generator bearing fault diagnosis method under variable speed | |
CN103412057A (en) | Weak oscillation attenuation signal parameter identification based on stochastic resonance and moving least square | |
Gao et al. | Impulsive gear fault diagnosis using adaptive Morlet wavelet filter based on alpha-stable distribution and kurtogram | |
CN107941510A (en) | Extracting method based on the angularly Rolling Bearing Fault Character of dual sampling | |
CN105388012A (en) | Order tracking method based on nonlinear frequency modulation wavelet transformation | |
Xiao et al. | A fault frequency bands location method based on improved fast spectral correlation to extract fault features in axial piston pump bearings | |
CN102721462B (en) | Method for quickly computing Bode plot and Nyquist plot of rotary mechanical vehicle starting and parking processes | |
CN102928225A (en) | Self-adaptive resonance demodulation method for fault diagnosis of rolling bearing | |
CN103217248A (en) | Method for detecting tensile force of bridge steel cable | |
Zhao et al. | Vibration health monitoring of rolling bearings under variable speed conditions by novel demodulation technique | |
CN104217112A (en) | Multi-type signal-based power system low-frequency oscillation analysis method | |
CN103267652B (en) | Intelligent online diagnosis method for early failures of equipment | |
Guo et al. | Ball screw fault detection and location based on outlier and instantaneous rotational frequency estimation | |
Yao et al. | A Fault Feature Extraction Method for Rolling Bearing Based on Pulse Adaptive Time‐Frequency Transform | |
CN103743969A (en) | Measurement noise suppression method for near-field electromagnetic radiation interference measurement signals | |
Zheng et al. | Zero-phase filter-based adaptive Fourier decomposition and its application to fault diagnosis of rolling bearing | |
Meng et al. | Maximum cyclic Gini index deconvolution for rolling bearing fault diagnosis | |
Dhillon | Weighted multi-sensor data level fusion method of vibration signal based on correlation function | |
Zhu et al. | Fault diagnosis method for rolling element bearings under variable speed based on TKEO and Fast-SC | |
CN114297579A (en) | Blind source separation method based on cyclostationary measure | |
Li et al. | Optimal maximum cyclostationary blind deconvolution for bearing fault detection | |
CN102798462B (en) | No-time-mark order tracking method based on self-demodulation transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |