CN113447264B - Online acoustic monitoring and diagnosing method for tooth breakage fault of planetary gear box - Google Patents

Online acoustic monitoring and diagnosing method for tooth breakage fault of planetary gear box Download PDF

Info

Publication number
CN113447264B
CN113447264B CN202110799918.7A CN202110799918A CN113447264B CN 113447264 B CN113447264 B CN 113447264B CN 202110799918 A CN202110799918 A CN 202110799918A CN 113447264 B CN113447264 B CN 113447264B
Authority
CN
China
Prior art keywords
fault
sound pressure
planetary gearbox
gear box
planetary gear
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.)
Active
Application number
CN202110799918.7A
Other languages
Chinese (zh)
Other versions
CN113447264A (en
Inventor
刘超
姚家驰
蒋东翔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN202110799918.7A priority Critical patent/CN113447264B/en
Publication of CN113447264A publication Critical patent/CN113447264A/en
Application granted granted Critical
Publication of CN113447264B publication Critical patent/CN113447264B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/021Gearings
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/028Acoustic or vibration analysis

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

An online acoustic monitoring and diagnosing method for a broken tooth fault of a planetary gear box comprises the steps of intercepting sound pressure signals of a plurality of working cycles of the planetary gear box through a shearing, filtering and normalization preprocessing method; decomposing the preprocessed sound pressure signal by using a multi-resolution singular value decomposition method, setting the optimal decomposition layer number, and obtaining a plurality of components through decomposition; calculating fault characteristic frequency of a broken tooth theory according to the structure of the planetary gearbox, constructing a reference signal based on the fault characteristic frequency, and selecting a fault component from a plurality of components by using an autoregressive model designed by the reference signal; and carrying out envelope spectrum analysis on the fault component, calculating the amplitude ratio of the fault characteristic frequency spectral line to the whole frequency spectrum spectral line, setting an amplitude ratio threshold value, and indicating that the tooth breaking fault occurs when the calculated amplitude ratio is greater than the amplitude ratio threshold value. According to the method, the state of the planetary gear box can be accurately monitored and diagnosed on line only by one sound pressure sensor, other signals do not need to be additionally measured or other monitoring devices do not need to be adopted, and the reliability of fault diagnosis is high.

Description

Online acoustic monitoring and diagnosing method for tooth breakage fault of planetary gear box
Technical Field
The invention relates to the technical field of acoustic monitoring and fault diagnosis of a planetary gear box, in particular to an online acoustic monitoring and diagnosis method for a broken tooth fault of the planetary gear box, which is applied to an actual engineering mechanical system.
Background
As a core device of a mechanical system, the working state of the planetary gearbox is related to the safe and stable operation of the whole mechanical system. The online monitoring and diagnosing system for the planetary gearbox can help users to find and eliminate faults in time, and plays an important role in ensuring stable and safe operation of a mechanical system. The traditional planetary gearbox monitoring technology is mainly a vibration monitoring technology, is widely applied to mechanical equipment systems such as motors, rails, automobiles, machine tools and the like, and obtains a good fault diagnosis effect, for example, in a published patent document, "a gearbox fault diagnosis method based on resonance baseband wide fourier decomposition" (CN 112539933A), a resonance baseband wide fourier decomposition method is used for monitoring and diagnosing gearbox faults by using vibration signals.
The vibration monitoring technology adopts a contact type measuring mode, however, in some practical engineering applications, the vibration monitoring technology is not suitable for being adopted when mechanical equipment is running or high corrosion, high-temperature and humid environment exist and the installation position is limited so that the accelerometer cannot be conveniently installed. At this time, a non-contact acoustic monitoring and diagnosis technology needs to be researched, for example, in a published patent document "fault diagnosis method for central transmission bevel gear of aircraft engine based on acoustic test" (CN 112284720 a), sound pressure signals are used to perform fault diagnosis for the central transmission bevel gear of aircraft engine by using order analysis, hilbert transform, envelope spectrum and cepstrum technology, however, the method needs to additionally measure rotating speed signals and introduce a sound conduit system, and the problem of endpoint effect and modal aliasing exists in an empirical mode decomposition algorithm used in hilbert transform, which affects the accuracy of fault diagnosis.
The structure of the planetary gear box is complex, the rotating speed and the load can be changed in practical engineering application, sound signals generated by gear meshing and vibration in the planetary gear box can be transmitted outwards through a complex transmission path, the sound pressure signals received by a sound pressure sensor outside the planetary gear box are weak, interference noise is high, visual and reliable sensitive characteristic information is difficult to extract from the existing acoustic monitoring and diagnosing method, and effective on-line acoustic monitoring and diagnosing can not be carried out on the running state of the planetary gear box.
Disclosure of Invention
The invention aims to provide an online acoustic monitoring and diagnosing method for a broken tooth fault of a planetary gear box, aiming at solving the technical problem of timely and accurately judging the operation state and the maintenance requirement of the planetary gear box, aiming at solving the technical problem of the existing acoustic monitoring and diagnosing technology in the process of monitoring and diagnosing the operation state of a complex planetary gear box, and providing accurate reference for the maintenance and the replacement of the planetary gear box.
The technical scheme of the invention is as follows:
an online acoustic monitoring and diagnosing method for a broken tooth fault of a planetary gearbox is characterized by comprising the following steps:
s1: collecting sound pressure signals:
measuring sound pressure signal P of planetary gear box at close distance position for a period of time by using sound pressure sensor 1
S2: from sound pressure signals P by window functions 1 Sound pressure signal P of middle-shear planetary gear box in t working cycle time 2 Removing the high-frequency interference signal by using a low-pass filter to obtain a filtered sound pressure signal P 3 The sound pressure signal is normalized to [ -1, +1 ] by the following equation]Obtaining a normalized sound pressure signal P 4
Figure BDA0003164280480000021
In the formula (1), P 4 For normalized sound pressure signal, P 3max Is a sound pressure signal P 3 Maximum value of, P 3min Is a sound pressure signal P 3 Minimum value of (d);
s3: for sound pressure signal P 4 And (3) decomposing:
determining the optimal decomposition layer number based on a noise energy ratio method, wherein a calculation formula of a noise energy ratio xi is as follows:
Figure BDA0003164280480000022
in the formula (2), xi is the noise energy ratio, and the xi value range is [0, 0.01 ]]J is the number of decomposition layers, E J Energy of the component decomposed for the J-th layer, E t Is the total energy of the signal to be decomposed;
calculating to obtain the optimal decomposition layer number J and J +1 decomposition components;
s4: constructing a reference signal:
calculating theoretical fault characteristic frequency f when tooth breaking fault occurs according to structure of planetary gear box fd The reference signal x (t) is constructed by:
Figure BDA0003164280480000023
in the formula (3), f fd Phi is a theoretical fault characteristic frequency, phi is a phase, and t is a working cycle time;
s5: selecting a fault component:
i) converting the reference signal x (t) into an autoregressive model:
Figure BDA0003164280480000024
in the formula (4), β i Is the coefficient of the autoregressive model, k is the order of the model, ε t Is a residual error;
ii) Z-transforming equation (4) to obtain a system function h (Z) as:
Figure BDA0003164280480000031
in formula (5), z is a complex variable;
iii) passing J +1 decomposition components through a system function H (z), and outputting a fault component most relevant to the reference signal and recording the fault component as C 1
S6: calculating the amplitude ratio:
i) calculating the envelope spectrum ES (C) 1 ):
ES(C 1 )=FT(|C 1 +j×h|) (6)
In equation (6), FT is Fourier transform,
Figure BDA0003164280480000032
Figure BDA0003164280480000033
ii) selecting a frequency spectrum of a frequency range, and calculating the amplitude ratio by using the following formula:
Figure BDA0003164280480000034
in the formula (7), A is the amplitude, A all Is the sum of the amplitudes of the entire spectral lines;
q is a value between 0 and 1, when Q is close to 0, the running state of the planetary gearbox is good, and when Q is close to 1, the running state of the planetary gearbox is poor;
s7: outputting a state monitoring result:
setting a fault alarm threshold Q 0 Comparing the amplitude ratio Q and Q obtained in step S6 0 If Q is greater than or equal to Q 0 An alarm command is sent to prompt the planetary gearbox to be repaired or replaced, otherwise, the steps S1 to S6 are repeatedly executed to monitor the running state of the planetary gearbox until the planetary gearbox is repaired or replaced.
In the technical scheme, the sound pressure signal P of the planetary gear box at the close-distance position for a period of time is measured 1 The close-range position refers to a position of the sound pressure sensor 1-3cm away from the planetary gearbox; the window function is preferablyBy usingA rectangular window, hanning window or hamming window function; the t work cycle time is preferably 8-12 work cycle times of the planetary gearbox; the low-pass filter for removing the high-frequency interference signals refers to the high-frequency interference signals which are larger than or equal to 400 Hz.
Further, the optimal number of decomposition layers J obtained by calculation in step S3 is preferably 7.
Further, the frequency spectrum of the frequency range in step S5 is preferably 20-25 times the theoretical fault characteristic frequency f fd The spectrum of the range.
Further, the failure alarm threshold Q in step S7 0 The value range of (A) is preferably 0.8-0.9 times of the amplitude ratio Q when the tooth breaking fault occurs.
Compared with the prior art, the invention has the following advantages and prominent technical effects: the online acoustic monitoring and diagnosing method provided by the invention can effectively monitor the running state of the planetary gear box in real time only by one sound pressure sensor, can early warn possible faults, does not need to additionally measure other signals or adopt other monitoring devices, and can timely and accurately judge the technical problems of the running state and the maintenance requirement of the planetary gear box, thereby providing accurate reference for the maintenance and the replacement of the planetary gear box. The acoustic monitoring and diagnosing method provided by the invention has no end effect and modal aliasing problems in an empirical mode decomposition algorithm, can accurately extract fault components and has high fault diagnosis reliability; the sound pressure sensor adopted by the invention can be arranged under the working state of the planetary gear box, is convenient and reliable to install, can greatly improve the convenience of online monitoring and diagnosis of the planetary gear box, and is suitable for practical engineering application.
Drawings
FIG. 1 is a flow chart of an online acoustic monitoring and diagnosing method for a planetary gearbox in an embodiment.
FIG. 2 shows the fault component C in the example 1 The time domain waveform of (a).
FIG. 3 shows the fault component C in the embodiment 1 The envelope spectrum of (a).
Detailed Description
The principles and operation of the present invention will be further described with reference to the drawings and examples, which are provided for illustration only and are not intended to limit the scope of the present invention.
The invention provides an online acoustic monitoring method for the operating state of a planetary gearbox, which specifically comprises the following steps:
s1: collecting sound pressure signals:
measuring sound pressure signal P of planetary gear box at close distance position for a period of time by using sound pressure sensor 1 (ii) a The close-up position is generally at the position of 1-3cm of the planetary gearbox;
s2: from sound pressure signals P by window functions 1 Sound pressure signal P of middle shear output planetary gear box in t working cycle time 2 Removing the high-frequency interference signal by using a low-pass filter to obtain a filtered sound pressure signal P 3 The sound pressure signal is normalized to [ -1, +1 ] by the following equation]Obtaining a normalized sound pressure signal P 4
Figure BDA0003164280480000041
In the formula (1), P 4 For normalized sound pressure signal, P 3max Is a sound pressure signal P 3 Maximum value of, P 3min Is a sound pressure signal P 3 Minimum value of (d);
the window function is used for intercepting a sound pressure signal P of the planetary gear box in t working cycle time 2 Typically a rectangular window, a hanning window or a hamming window function. The t duty cycle times are preferably 8-12 duty cycle times for the planetary gearbox. And the low-pass filter is used for removing the high-frequency interference signals with the frequency of more than or equal to 400 Hz.
S3: sound pressure signal P by using multi-resolution singular value decomposition method 4 And (3) decomposing:
sound pressure signal P by multiresolution singular value decomposition method 4 The number of the components obtained after decomposition is determined by the number of decomposition layers, if the number of the decomposition layers is large, the calculated amount is large, and if the number of the decomposition layers is small, the decomposition components contain more interference noise components to influence the monitoring and diagnosis effect; determining the optimal decomposition layer number based on a noise energy ratio method, wherein a calculation formula of a noise energy ratio xi is as follows:
Figure BDA0003164280480000051
in the formula (2), xi is the noise energy ratio, J is the number of decomposition layers, and E J Energy of the component decomposed for the J-th layer, E t Is the total energy of the signal to be decomposed;
when the noise energy ratio xi is less than or equal to 1%, the interference noise component contained in the decomposed component does not affect the monitoring and diagnosis effect, and the value range of the noise energy ratio xi is set to be [0, 0.01 ]]Calculating to obtain the optimal decomposition layer number J of 7, and performing multi-resolution singular value decomposition on the normalized sound pressure signal P 4 After decomposition, 8 decomposition components can be obtained;
s4: constructing a reference signal:
calculating theoretical fault characteristic frequency f when tooth breaking fault occurs according to structure of planetary gear box fd Then, thenThe reference signal x (t) is constructed by:
Figure BDA0003164280480000052
in the formula (3), f fd Phi is the theoretical fault characteristic frequency and phi is the phase;
s5: selecting a fault component:
i) converting the reference signal x (t) into an autoregressive model:
Figure BDA0003164280480000053
in the formula (4), β i Is the coefficient of the autoregressive model, k is the order of the model, ε t Is a residual error;
ii) Z-transforming equation (4) to obtain a system function h (Z) as:
Figure BDA0003164280480000061
in formula (5), z is a complex variable;
iii) passing J +1 decomposition components through a system function H (z), and outputting a fault component most relevant to the reference signal and recording the fault component as C 1
S6: calculating an amplitude ratio:
i) calculating the envelope spectrum ES (C) 1 ):
ES(C 1 )=FT(|C 1 +j×h|) (6)
In equation (6), FT is Fourier transform,
Figure BDA0003164280480000062
Figure BDA0003164280480000063
ii) selecting a frequency spectrum of a frequency range, namely extracting 20-25 times of theoretical fault characteristic frequency f fd The following frequency spectrum, benefitThe amplitude ratio is calculated using the following equation:
Figure BDA0003164280480000064
in the formula (7), A is the amplitude, A all Is the sum of the amplitudes of the entire spectral lines;
q is a value between 0 and 1, when Q is close to 0, the running state of the planetary gearbox is good, and when Q is close to 1, the running state of the planetary gearbox is poor;
s7: outputting a state monitoring result:
setting a fault alarm threshold Q 0 Threshold value Q of fault alarm 0 The value range of (A) is generally 0.8-0.9 times of the amplitude ratio Q when the tooth breaking fault occurs. Comparing the amplitude ratio Q and Q obtained in step S6 0 If Q is greater than or equal to Q 0 If so, sending an alarm command to prompt the planetary gearbox to be maintained or replaced; otherwise, the steps S1 to S6 are repeatedly performed to monitor the planetary gearbox operation state until the planetary gearbox is repaired or replaced.
Example (b):
the embodiment of the invention takes a three-stage planetary gearbox in practical engineering application as an example, and the method provided by the invention is adopted to carry out online acoustic monitoring diagnosis on the three-stage planetary gearbox, and the method specifically comprises the following steps:
s1: collecting sound pressure signals:
measuring sound pressure signal P at a position 1cm away from the three-stage planetary gear box by using 1 sound pressure sensor 1 The sampling frequency is set to be 25.6kHz, the sensitivity of the sound pressure sensor is-25.8 dB (0dB ═ 1V/Pa, at 1kHz), and because the amount of the collected test data is large, only part of the data is shown in table 1;
TABLE 1 Sound pressure Signal P 1
Data points 1 2 3 4 5 6 7 8 9
P 1 (mV) 0.0397 0.0306 0.0200 0.0167 0.0186 0.0203 0.0208 0.0181 0.0116
Data points 10 11 12 13 14 15 16 17 18
P 1 (mV) 0.0007 -0.0068 -0.0093 -0.0104 -0.0084 -0.0054 -0.0013 0.0042 0.0026
S2: using rectangular windows to derive the sound pressure signal P 1 Acoustic pressure signal P of 10 working cycles of three-stage planetary gear box is cut out in middle 2 Removing high-frequency interference signals larger than or equal to 400Hz by using a 6-order Butterworth low-pass filter to obtain a filtered sound pressure signal P 3 The sound pressure signal is normalized to [ -1, +1 ] by the following equation]Obtaining a normalized sound pressure signal P 4
Figure BDA0003164280480000071
In the formula (1), P 4 For normalized sound pressure signal, P 3max Is a sound pressure signal P 3 Maximum value of, P 3min Is a sound pressure signal P 3 Minimum value of (d);
due to the sound pressure signal P 4 The data volume is large, and only part of the data is shown in table 2;
TABLE 2 Sound pressure Signal P 4
Data points 1 2 3 4 5 6 7 8 9
P 4 (mV) -0.0580 -0.0547 -0.0513 -0.0478 -0.0443 -0.0407 -0.0371 -0.0334 -0.0298
Data points 10 11 12 13 14 15 16 17 18
P 4 (mV) -0.0261 -0.0224 -0.0187 -0.0151 -0.0114 -0.0078 -0.0042 -0.0006 0.0030
S3: sound pressure signal P by using multi-resolution singular value decomposition method 4 And (3) decomposing:
determining the optimal decomposition layer number based on a noise energy ratio method, wherein a calculation formula of a noise energy ratio xi is as follows:
Figure BDA0003164280480000072
in the formula (2), xi is the noise energy ratio, J is the number of decomposition layers, and E J Energy of the component decomposed for the J-th layer, E t Is the total energy of the signal to be decomposed;
when the noise energy ratio xi is less than or equal to 1%, the interference noise component contained in the decomposed component does not affect the monitoring and diagnosis effect, and the value range of the noise energy ratio xi is set to be [0, 0.01 ]]Calculating to obtain the optimal decomposition layer number J of 7, and decomposing the sound pressure signal P by using a multi-resolution singular value decomposition method 4 8 decomposition components are obtained after decomposition;
s4: constructing a reference signal:
by calculatingTheoretical fault characteristic frequency f when tooth breaking fault occurs to the outer gear ring of the three-stage planetary gear box fd At 0.13Hz, and phase phi is set to 0, the constructed reference signal is:
Figure BDA0003164280480000081
s5: selecting a fault component:
i) converting the reference signal x (t) into an autoregressive model:
Figure BDA0003164280480000082
in the formula (4), β i Is the coefficient of the autoregressive model, k is the order of the model, ε t Is a residual error;
ii) Z-transforming equation (4) to obtain a system function h (Z) as:
Figure BDA0003164280480000083
in formula (5), z is a complex variable;
iii) passing 8 decomposition components through a system function H (z), and outputting a fault component most relevant to the reference signal, and recording the fault component as C 1
In this embodiment, the order k of the autoregressive model is 5, and the coefficients of the 5-order autoregressive model obtained by calculation are [ -2.02567, 1.0843, -0.0999, 0.0494, -0.0082]The mean square error values obtained by passing 8 decomposition components through a system function H (z) are respectively [4.5880e-12, 7.5273e-13, 4.2112e-13, 3.1127e-13, 2.5798e-13, 2.2672e-13, 2.0618e-13, 5.3462e-08]Wherein the minimum mean square error value of the 7 th decomposition component indicates that the 7 th decomposition component is a fault component, which is denoted as C 1 The time domain waveform is shown in fig. 2;
s6: calculating an amplitude ratio:
i) calculating the envelope spectrum ES (C) using the following equation 1 ):
ES(C 1 )=FT(|C 1 +j×h|) (6)
In equation (6), FT is Fourier transform,
Figure BDA0003164280480000091
Figure BDA0003164280480000092
ii) in this embodiment, the spectrum below 2Hz is extracted and analyzed, and the amplitude ratio is calculated by the following formula:
Figure BDA0003164280480000093
in the formula (7), A is the amplitude, A all Is the sum of the amplitudes of the entire spectral lines;
q is a value between 0 and 1, when Q is close to 0, the operation state of the three-stage planetary gearbox is good, when Q is close to 1, the operation state of the three-stage planetary gearbox is poor, and a fault component C is obtained through calculation 1 The envelope spectrum of (a) is as shown in fig. 3, and the amplitude ratio Q when a tooth breakage fault occurs is 0.1089;
s7: sending alarm information:
setting a fault alarm threshold Q 0 Comparing the amplitude ratio Q and Q obtained at S6 as 0.0926 0 After the size of the (B) is more than or equal to Q 0 And sending alarm information to prompt that the three-stage planetary gear box has a tooth breaking fault and needs to be maintained or replaced.

Claims (8)

1. An online acoustic monitoring and diagnosing method for a broken tooth fault of a planetary gearbox is characterized by comprising the following steps:
s1: collecting sound pressure signals:
measuring sound pressure signal P of planetary gear box at close distance position for a period of time by using sound pressure sensor 1
S2: from sound pressure signals P by window functions 1 Sound pressure signal P of middle-shear planetary gear box in t working cycle time 2 Using low-pass filtersRemoving the high-frequency interference signal to obtain a filtered sound pressure signal P 3 The sound pressure signal is normalized to [ -1, +1 ] by the following equation]Obtaining a normalized sound pressure signal P 4
Figure FDA0003589909540000011
In the formula (1), P 4 For normalized sound pressure signal, P 3max Is a sound pressure signal P 3 Maximum value of, P 3min Is a sound pressure signal P 3 Minimum value of (d);
s3: sound pressure signal P by using multi-resolution singular value decomposition method 4 And (3) decomposing:
determining the optimal decomposition layer number based on a noise energy ratio method, wherein a calculation formula of a noise energy ratio xi is as follows:
Figure FDA0003589909540000012
in the formula (2), xi is the noise energy ratio, and the xi value range is [0, 0.01 ]]J is the number of decomposition layers, E i Energy of the decomposed component for the ith layer, E t Is the total energy of the signal to be decomposed;
calculating to obtain the optimal decomposition layer number J and J +1 decomposition components;
s4: constructing a reference signal:
calculating theoretical fault characteristic frequency f when tooth breaking fault occurs according to structure of planetary gear box fd The reference signal x (t) is constructed by:
Figure FDA0003589909540000013
in the formula (3), f fd Phi is a theoretical fault characteristic frequency, phi is a phase, and t is a working cycle time;
s5: selecting a fault component:
converting the reference signal x (t) into an autoregressive model:
Figure FDA0003589909540000021
in the formula (4), β i Is the coefficient of the autoregressive model, k is the order of the model, ε t Is a residual error;
ii) Z-transforming equation (4) to obtain a system function h (Z) as:
Figure FDA0003589909540000022
in formula (5), z is a complex variable;
iii) passing J +1 decomposition components through a system function H (z), and outputting a fault component most relevant to the reference signal and recording the fault component as C 1
S6: calculating an amplitude ratio:
i) calculating the envelope spectrum ES (C) 1 ):
ES(C 1 )=FT(|C 1 +j×h|) (6)
In equation (6), FT is Fourier transform,
Figure FDA0003589909540000023
ii) selecting a frequency spectrum of a frequency range, and calculating the amplitude ratio Q by using the following formula:
Figure FDA0003589909540000024
in the formula (7), A is the amplitude, A all Is the sum of the amplitudes of the entire spectral lines;
q is a value between 0 and 1, when Q is close to 0, the running state of the planetary gearbox is good, and when Q is close to 1, the running state of the planetary gearbox is poor;
s7: outputting a state monitoring result:
is provided withFailure alarm threshold Q 0 Comparing the amplitude ratio Q and Q obtained in step S6 0 If Q is greater than or equal to Q 0 And sending an alarm command to prompt the planetary gearbox to be repaired or replaced, otherwise, repeatedly executing the steps S1 to S6 until the planetary gearbox is repaired or replaced.
2. The on-line acoustic monitoring and diagnosing method for the broken tooth fault of the planetary gearbox as claimed in claim 1, wherein: the close-range position refers to a position where the sound pressure sensor is 1-3cm away from the planetary gear box.
3. The on-line acoustic monitoring and diagnosing method for the broken tooth fault of the planetary gearbox as claimed in claim 1, wherein the method comprises the following steps: the window function is a rectangular window, a hanning window or a hamming window.
4. The on-line acoustic monitoring and diagnosing method for the broken tooth fault of the planetary gearbox as claimed in claim 1, wherein the method comprises the following steps: the t work cycle time is 8-12 work cycle times of the planetary gearbox.
5. The on-line acoustic monitoring and diagnosing method for the broken tooth fault of the planetary gearbox as claimed in claim 1, wherein the method comprises the following steps: the high-frequency interference signal is a signal greater than or equal to 400 Hz.
6. The on-line acoustic monitoring and diagnosing method for the broken tooth fault of the planetary gearbox as claimed in claim 1, wherein the method comprises the following steps: the optimal number of decomposition layers J is 7.
7. The on-line acoustic monitoring and diagnosing method for the broken tooth fault of the planetary gearbox as claimed in claim 1, wherein the method comprises the following steps: the frequency spectrum of the frequency range is 20-25 times of theoretical fault characteristic frequency f fd The spectrum of the range.
8. The on-line acoustic monitoring and diagnosing method for the broken tooth fault of the planetary gearbox as claimed in claim 1, wherein the method comprises the following steps: the failure alarm thresholdValue Q 0 The value range of (A) is 0.8-0.9 times of the amplitude ratio Q when the tooth breaking fault occurs.
CN202110799918.7A 2021-07-15 2021-07-15 Online acoustic monitoring and diagnosing method for tooth breakage fault of planetary gear box Active CN113447264B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110799918.7A CN113447264B (en) 2021-07-15 2021-07-15 Online acoustic monitoring and diagnosing method for tooth breakage fault of planetary gear box

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110799918.7A CN113447264B (en) 2021-07-15 2021-07-15 Online acoustic monitoring and diagnosing method for tooth breakage fault of planetary gear box

Publications (2)

Publication Number Publication Date
CN113447264A CN113447264A (en) 2021-09-28
CN113447264B true CN113447264B (en) 2022-08-09

Family

ID=77816381

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110799918.7A Active CN113447264B (en) 2021-07-15 2021-07-15 Online acoustic monitoring and diagnosing method for tooth breakage fault of planetary gear box

Country Status (1)

Country Link
CN (1) CN113447264B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005300517A (en) * 2004-03-17 2005-10-27 Nippon Steel Corp Diagnostic method and device for abnormality of roller bearing
CN107246967A (en) * 2017-07-07 2017-10-13 武汉钢铁有限公司 Signal processing method and device for gear box arrangement fault diagnosis
CN109029977A (en) * 2018-07-12 2018-12-18 福州大学 A kind of epicyclic gearbox Incipient Fault Diagnosis method based on VMD-AMCKD
CN112284720A (en) * 2020-10-16 2021-01-29 中国航发四川燃气涡轮研究院 Acoustic test-based fault diagnosis method for central transmission bevel gear of aircraft engine
CN113074935A (en) * 2021-04-01 2021-07-06 西华大学 Acoustic separation and diagnosis method for impact fault characteristics of gearbox

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005300517A (en) * 2004-03-17 2005-10-27 Nippon Steel Corp Diagnostic method and device for abnormality of roller bearing
CN107246967A (en) * 2017-07-07 2017-10-13 武汉钢铁有限公司 Signal processing method and device for gear box arrangement fault diagnosis
CN109029977A (en) * 2018-07-12 2018-12-18 福州大学 A kind of epicyclic gearbox Incipient Fault Diagnosis method based on VMD-AMCKD
CN112284720A (en) * 2020-10-16 2021-01-29 中国航发四川燃气涡轮研究院 Acoustic test-based fault diagnosis method for central transmission bevel gear of aircraft engine
CN113074935A (en) * 2021-04-01 2021-07-06 西华大学 Acoustic separation and diagnosis method for impact fault characteristics of gearbox

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Yao Jiachi, et al..Fault detection of complex planetary gearbox using acoustic signals.《Meaurement》.2021,(第178期),第3节. *
Yao Jiachi, et al..Fault diagnosis of planetary gearbox based on acoustic signals.《Applied Acoustics》.2021,(第181期),第2、3节. *
Zhang Gang.Research on a Noise Reduction Method Based on Multi-Resolution Singular Value Decomposition.《Applied Sciences》.2020,第2、3节. *

Also Published As

Publication number Publication date
CN113447264A (en) 2021-09-28

Similar Documents

Publication Publication Date Title
CN108507787B (en) Wind power gear speed increasing box fault diagnosis test platform and method based on multi-feature fusion
Yang et al. Vibration feature extraction techniques for fault diagnosis of rotating machinery: a literature survey
Li et al. A new noise-controlled second-order enhanced stochastic resonance method with its application in wind turbine drivetrain fault diagnosis
Guo et al. Faulty bearing signal recovery from large noise using a hybrid method based on spectral kurtosis and ensemble empirical mode decomposition
CN104048825B (en) A kind of gearbox of wind turbine Fault Locating Method of Multi-sensor Fusion
Yoon et al. On the use of a single piezoelectric strain sensor for wind turbine planetary gearbox fault diagnosis
CN201311360Y (en) Vibrating online monitoring failure diagnostic device
Farokhzad Vibration based fault detection of centrifugal pump by fast fourier transform and adaptive neuro-fuzzy inference system
CN104977047A (en) Wind turbine online condition monitoring and health assessment system and method thereof based on vibration and oil
CN105003453A (en) Online monitoring and fault diagnosis system of mine fan
CN204113701U (en) A kind of mine fan on-line monitoring and fault diagnosis system
Yao et al. Fault detection of complex planetary gearbox using acoustic signals
CN109596349A (en) A kind of decelerator trouble diagnostic method based on VMD and PCT
CN114371001B (en) Gear box fault defect detection system
CN107229269A (en) A kind of wind-driven generator wheel-box method for diagnosing faults of depth belief network
CN106769053A (en) A kind of hydraulic turbine fault diagnosis system and method based on acoustic emission signal
CN107588802A (en) Hoisting equipment reductor vibration monitor system and method based on fibre optical sensor
CN114964769A (en) Wind power gear box vibration signal fault diagnosis method
CN113447264B (en) Online acoustic monitoring and diagnosing method for tooth breakage fault of planetary gear box
CN113280910A (en) Real-time monitoring method and system for long product production line equipment
CN116738314A (en) Air compressor fault detection method based on neural network
CN113688791B (en) Method for identifying CMS abnormal data packet of wind turbine generator
Yao et al. A low-frequency fault detection method for low-speed planetary gearbox based on acoustic signals
Yang et al. Intelligent diagnosis technology of wind turbine drive system based on neural network
CN116499742A (en) Rotary machine fault monitoring method based on ensemble learning

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