CN114813117B - Fault diagnosis method and device for RV reducer - Google Patents
Fault diagnosis method and device for RV reducer Download PDFInfo
- Publication number
- CN114813117B CN114813117B CN202210385732.1A CN202210385732A CN114813117B CN 114813117 B CN114813117 B CN 114813117B CN 202210385732 A CN202210385732 A CN 202210385732A CN 114813117 B CN114813117 B CN 114813117B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- reducer
- fault diagnosis
- fault
- 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
Links
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
The invention discloses a fault diagnosis method and device for an RV reducer. The invention firstly extracts a signal of the whole cycle movement, then carries out time-frequency analysis on the whole cycle signal by continuous wavelet transform, and intercepts the signal of the stable stage in the ascending movement. By intercepting the step, the non-stationary vibration component in the analyzed signal is reduced, so that the deviation between each theoretical calculation value and each actual calculation value of the signal is reduced as much as possible, and the non-stationary signal is approximately converted into a stationary signal. Then, considering the NeighCoeff method, setting the threshold value by taking the adjacent coefficients as a whole, not only can more characteristic information be kept, but also the noise reduction effect is better than that of the wavelet threshold value. Therefore, the invention uses the NeighCoeff method to filter and reduce noise of the intercepted signal while keeping more impact information. And finally, carrying out FFT on the noise reduction signal, and carrying out fault diagnosis on the transformed frequency domain signal.
Description
Technical Field
The invention relates to a fault diagnosis method and device for an RV reducer, and belongs to the technical field of state monitoring and fault diagnosis of mechanical equipment.
Background
With the progress of modern manufacturing industry, high precision, high rotation speed and stable control gradually become the development direction of industrial robots, and the cooperation of high-precision speed reducers cannot be avoided. The speed reducer is used as a core component of the joint of the industrial robot, and the health condition of the speed reducer determines whether the industrial robot can normally operate. The RV reducer (Rotate Vector reducer) is the most commonly used reducer for robots due to its small size, large transmission ratio, high efficiency, etc. Since the RV reducer is sealed inside the joint as a high-precision core component, when the reducer fails, how to quickly complete the failure diagnosis becomes an inevitable problem.
At present, the fault diagnosis research on rolling bearings, or gears, has never been interrupted. Most researches on RV reducers relate to stress analysis, dynamic modeling, testing devices and the like, and fault diagnosis is rarely related. The reason is that the working condition of the RV reducer is usually intermittent reciprocating motion and variable-speed motion, so that the fault signal of the RV reducer is a non-stationary signal, and the fault characteristics cannot be accurately identified by the traditional signal processing method.
Disclosure of Invention
The invention provides a fault diagnosis method and device for an RV reducer, which are used for realizing fault diagnosis of the RV reducer under the condition of reciprocating periodic intermittent motion.
The technical scheme of the invention is as follows: a fault diagnosis method for an RV reducer, comprising:
s1, an acceleration sensor 2 is used for picking up an observation signal S (n) of mechanical vibration in the reciprocating intermittent motion process of a joint arm 1 and calculating the relevant theoretical characteristic frequency of an RV reducer; wherein n is the sampling time of the vibration signal s (n);
s2, intercepting the whole-cycle signal of the observation signal S (n) obtained in the step S1 to obtain a periodic motion signal S simultaneously comprising an ascending stage and a descending stage 1 (n)；
S3, obtaining the periodic motion signal S obtained in the step S2 1 (n) performing continuous wavelet transform, and synchronously intercepting periodic motion signal s according to time-frequency diagram 1 (n) rising phase, approximately stationary vibration signal of sample length m, as preprocessing signal s 2 (n)；
S4, preprocessing signal S obtained in the step S3 2 (n) filtering to obtain a filtered signal x (n);
and S5, carrying out FFT on the filtering signal x (n) obtained in the step S4, and carrying out fault diagnosis on the frequency domain signal after transformation.
Of the RV reducerThe theoretical characteristic frequency of interest including the sun gear revolution f 1 Planet wheel frequency conversion f 2 Support plate frequency conversion f 3 Planet wheel failure frequency
Engagement frequency f of sun gear and planet gear 1c 。
The specific steps for calculating the relevant theoretical characteristic frequency of the RV reducer are as follows:
s1.1, knowing the preset speed n of the articulated arm 1 3 By calculating the rotational speed n 3 Conversion to a rotary frequency, and the knuckle arm rotary frequency is equivalent to a supporting disc rotary frequency f 3 (ii) a According to the formula
Calculating to obtain the main shaft rotating speed n of the motor 5 1 (ii) a Wherein n is 1 The spindle speed of the motor 5; z is a linear or branched member 1 Is the number of sun gear teeth, Z 2 Number of teeth of planet gear, Z 4 The number of teeth of the pin gear; />
S1.2, by following the formula
The revolution frequency f of the sun gear is obtained through calculation 1 ；
S1.3 according to the formula
And &>
Calculating the planet wheel rotation frequency f of the RV reducer 2 And the engagement frequency f of the sun wheel and the planet wheel 1c (ii) a Wherein Z is 3 The number of teeth of the cycloid gear is shown;
s1.4, failure frequency of planetary gear
For the rotational frequency of the planet wheel relative to the planet carrier, the formula is
In the continuous wavelet transformation, morlet wavelets are selected as mother wavelets.
The sample length m is more than or equal to 10F max (ii) a Wherein, F max Representing 2 times the highest frequency value among the relevant theoretical characteristic frequencies of the RV reducer.
Preprocessing the signal s using the NeighCoeff algorithm 2 (n) filtering.
The fault diagnosis of the transformed frequency domain signal specifically comprises: carrying out multi-tooth wear fault diagnosis on the planet gear on the frequency domain signal after FFT, and when the frequency domain spectral line and the planet gear fault frequency
When the absolute value of the error is more than 0.2Hz, the multi-tooth wear fault of the planet wheel is determined not to occur; when the error between the frequency domain spectral line and the fault frequency of the planet wheel is between minus 0.2Hz and 0Hz or the value is 0, the planet wheel is determined to be a multi-tooth wear fault; when the error between the frequency domain spectral line and the planet wheel fault frequency is between 0Hz and 0.2Hz, cepstrum analysis is carried out on the filtering signal x (n), and fault diagnosis is realized.
A fault diagnosis device for an RV reducer, comprising:
the signal pickup module is used for picking up an observation signal s (n) of mechanical vibration in the reciprocating intermittent motion process of the articulated arm 1 through the acceleration sensor 2 and calculating the relevant theoretical characteristic frequency of the RV reducer; wherein n is the sampling time of the vibration signal s (n);
a first obtaining module for intercepting the observation signal s (n) obtained by the signal pickup module to obtain a periodic motion signal s including a rising stage and a falling stage 1 (n)；
An intercepting module for intercepting the periodic motion signal s obtained by the first obtaining module 1 (n) performing continuous wavelet transform, and synchronously intercepting periodic motion signal s according to time-frequency diagram 1 (n) rising phase, approximately stationary vibration signal of sample length m, as preprocessing signal s 2 (n)；
A second obtaining module for obtaining the preprocessed signal s 2 (n) filtering to obtain a filtered signal x (n);
and the fault diagnosis module is used for carrying out FFT on the filtering signal x (n) obtained by the third obtaining module and carrying out fault diagnosis on the frequency domain signal after transformation.
A processor for executing a program, wherein the program when executed performs any of the above-described fault diagnosis methods for a RV reducer.
A computer-readable storage medium including a stored program, wherein when the program runs, the apparatus where the computer-readable storage medium is located is controlled to execute the fault diagnosis method for the RV reducer.
The invention has the beneficial effects that: the method for intercepting the steady-phase signals through time-frequency analysis provided by the invention has the advantages that the error between the actual calculation result and the theoretical calculation value is smaller; the method can complete the fault identification of the RV reducer under the condition of unstable vibration.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is an actual layout diagram of an RV reducer fault simulation experiment table in the invention;
FIG. 3 is a 3D schematic diagram of an RV reducer fault simulation experiment table according to the invention;
FIG. 4 is a diagram showing a planet wheel failure in the present invention;
FIG. 5 is a time domain waveform of a vibration observation signal picked up by a sensor when the RV reducer fault simulation experiment table runs, and the operation of the vibration observation signal after the whole-cycle signal is intercepted is carried out;
FIG. 6 is a time-frequency spectrum of vibration observation signals picked up by a sensor when the bearing fault simulation experiment table runs; and intercepts the stable part in the rising phase of the observed signal (intercepts about 0.5s of signal);
FIG. 7 is a time domain waveform of a vibration observation signal after a stabilization phase is intercepted by time-frequency analysis according to the present invention;
FIG. 8 is a time domain waveform of a steady-state signal filtered by a NeighCoeff algorithm according to the present invention;
FIG. 9 is a frequency domain waveform of a filtered signal after FFT in accordance with the present invention;
FIG. 10 is a cepstral diagram of the filtered signal after cepstral analysis in accordance with the present invention;
in FIGS. 2 and 3, the reference numbers are as follows: the test bench comprises a joint arm 1, an acceleration sensor 2, a speed reducer retainer 3, an RV speed reducer 4, a motor 5 and a test bench base 6.
Detailed Description
The invention will be further described with reference to the following figures and examples, but the scope of the invention is not limited thereto.
Example 1: as shown in fig. 1-10, a fault diagnosis method for an RV reducer includes:
s1, firstly, mounting an acceleration sensor 2 on the upper surface of a reducer retainer 3, picking up an observation signal S (n) of mechanical vibration in the reciprocating intermittent motion process of a joint arm 1 through the acceleration sensor 2, and calculating the relevant theoretical characteristic frequency of an RV reducer; wherein n is the sampling time of the vibration signal s (n);
s2, intercepting the whole-cycle signal of the observation signal S (n) obtained in the step S1 to obtain a periodic motion signal S simultaneously comprising an ascending stage and a descending stage 1 (n)；
S3, obtaining the periodic motion signal S obtained in the step S2 1 (n) performing continuous wavelet transform, and synchronously intercepting periodic motion signal s according to time-frequency diagram 1 (n) rising phase, approximately stationary vibration signal of sample length m, as preprocessing signal s 2 (n)；
S4, preprocessing signal S obtained in the step S3 2 (n) filtering to obtain a filtered signal x (n);
and S5, carrying out FFT on the filtering signal x (n) obtained in the step S4, and carrying out fault diagnosis on the frequency domain signal after transformation.
Further, it may be provided that the relevant theoretical characteristic frequency of the RV reducer comprises a sun gear revolution frequency f 1 Planetary gear frequency conversion f 2 Supporting disk frequency conversion f 3 Planet wheel failure frequency
Engagement frequency f of sun gear and planet gear 1c 。
Since the RV reducer fault is mostly occurred in the primary speed reducing mechanism (i.e. the planetary gear train), further, the specific steps of calculating the relevant theoretical characteristic frequency of the RV reducer can be set as follows:
s1.1, knowing the preset speed n of the articulated arm 1 3 By calculating the rotational speed n 3 Conversion to a rotary frequency, and the knuckle arm rotary frequency is equivalent to a supporting disc rotary frequency f 3 (ii) a According to the formula
Calculating to obtain the main shaft rotating speed n of the motor 5 1 (ii) a Wherein n is 1 The spindle speed of the motor 5; z 1 Is the number of sun gear teeth, Z 2 Number of teeth of planet gear, Z 4 The number of teeth of the pin gear; />
S1.2, by following the formula
The revolution frequency f of the sun gear is obtained through calculation 1 ；
S1.3, according to the formula
And &>
Calculating the rotation frequency f of the planet gear of the RV reducer 2 And the engagement frequency (i.e. primary engagement frequency) f of the sun and planet gears 1c (ii) a Wherein Z is 3 The number of teeth of the cycloid gear is shown;
s1.4, failure frequency of planetary gear
For the rotational frequency of the planet wheel relative to the planet carrier, the calculation formula is
Further, in the continuous wavelet transform, a Morlet wavelet may be used as a mother wavelet.
Further, the sample length m is set to be more than or equal to 10F max (ii) a Wherein, F max Representing 2 times the highest frequency value among the relevant theoretical characteristic frequencies of the RV reducer. Because the signal collected by the RV reducer is a transient impact signal, when the intercepted sample length does not meet the requirement, the amplitude at the moment of impact cannot be captured, obviously, the sample length mode intercepted by the method also considers the signal frequency which is more than 10 times on the basis of considering 2 times of the highest frequency value in the related theoretical characteristic frequencies of the RV reducer, thereby not only ensuring the frequency of the signal not to be distorted, but also further ensuring the amplitude of the signal not to be distorted, and ensuring that the collected signal can be fit to the reality. Specifically, in the embodiment of the present invention, the maximum frequency of the analysis is 465.87Hz, and considering that 2-fold frequency may be needed, the number of sample points should be not less than 9300, so in summary, the length m of the sample obtained by the analysis is 12.5k (i.e. 0.5 s).
Further, it may be arranged to use the NeighCoeff algorithm on the preprocessed signal s 2 (n) filtering.
Further, the fault diagnosis of the transformed frequency domain signal may be set to specifically: carrying out multi-tooth wear fault diagnosis on the planet gear on the frequency domain signal after FFT, and when the frequency domain spectral line and the planet gear fault frequency
When the absolute value of the error is more than 0.2Hz (namely more than or equal to 0.2 Hz), determining that the multi-tooth wear fault of the planet wheel does not occur; when the error between the frequency domain spectral line and the fault frequency of the planet wheel is between-0.2 Hz and 0Hz or the value is 0 (namely the value between-0.2 Hz and 0Hz or the value is 0), the multi-tooth wear fault of the planet wheel is determined; when the error between the frequency domain spectral line and the planet wheel fault frequency is between 0Hz and 0.2Hz (namely, a numerical value between 0Hz and 0.2 Hz), cepstrum analysis is carried out on the filtering signal x (n), and fault diagnosis is realized.
Still further, in step S3, the specific steps of performing continuous wavelet transform on the periodic signal are as follows:
s3.1, setting function psi ∈ L 2 (R)∩L 1 (R) and
from the psi, a cluster of functions can be derived by scaling and translation:
in the formula, a, b belongs to R, a is not equal to 0 and is called psi a,b Is a continuous wavelet, a is a scale factor, b is a translation factor, and psi is a mother wavelet. a is used for changing the shape of the continuous wavelet, and b is used for changing the displacement of the continuous wavelet;
s3.2, for any function, the continuous wavelet transform is defined as:
in the formula, < f, ψ a,b The inner product of the two functions is represented by > x,
represents the complex conjugate of ψ (t);
s3.3, each wavelet transformation coefficient W can be obtained through the sub-wavelet with the scale a and the translation b and the signal inner product a,b (a, b) when the signal is more similar to the wavelet, the larger the coefficient value, the more the feature will converge somewhere on the time scale phase plane into a high amplitude block of energy, otherwise the energy will diverge. Finally obtaining a coefficient matrix W through different scale transformation f (a, b) are capable of characterizing a two-dimensional time scale of the signal. The Morlet wavelet is selected as the mother wavelet considering that it is similar to the shape of the impact signal generated when the rotary machine fails.
In S4, the preprocessing signal S obtained in the step S3 is subjected to NeighCoeff algorithm 2 (n) is carried outFiltering, specifically:
s4.1, performing discrete wavelet transform on a signal containing noise;
s4.2, grouping the wavelet coefficients into blocks with length L for each scale j
S4.3, for each block
Estimating new coefficients thereof by using a contraction rule;
in the formula: λ is a parameter used to adjust the threshold, and the length L is L = ln n, where n is the signal length;
is determined by the following formula:
and S4.4, performing inverse wavelet transform on the obtained wavelet coefficient to obtain a new signal, and completing signal noise reduction.
The working condition of the existing industrial robot is usually intermittent reciprocating motion, so that the acquired signal is usually characterized by unstable vibration. The invention firstly extracts a signal of the whole cycle movement, then carries out time-frequency analysis on the whole cycle signal by continuous wavelet transform, and intercepts the signal of the stable stage in the ascending movement. By intercepting the step, the non-stationary vibration component in the analyzed signal is reduced, so that the deviation between each theoretical calculation value and each actual calculation value of the signal is reduced as much as possible, and the non-stationary signal is approximately converted into a stationary signal. Then, considering the NeighCoeff method, setting the threshold value by taking the adjacent coefficients as a whole, not only can more characteristic information be kept, but also the noise reduction effect is better than that of the wavelet threshold value. Therefore, the invention uses the NeighCoeff method to carry out filtering and noise reduction on the intercepted signal while reserving more impact information. And finally, carrying out FFT on the noise reduction signal, and carrying out fault diagnosis on the transformed frequency domain signal.
Example 2: as shown in fig. 1-10, and further, in conjunction with experimental data, the present invention presents alternative embodiments as follows:
a mechanical fault diagnosis method for an RV reducer is characterized in that a test bed used in the example is an RV reducer test bed for simulating joint motion of a robot, and FIG. 3 shows collection positions of the test bed and sensors, wherein the test bed comprises a joint arm 1, an acceleration sensor 2, a reducer retainer 3, an RV reducer 4, a motor 5 and a test bed base 6; wherein, a speed reducer retainer 3 is arranged on the test bed base 6, an RV speed reducer 4 is arranged through the speed reducer retainer 3, and an output shaft of the motor 5 is connected with the knuckle arm 1 through the RV speed reducer 4. RV40E model reduction gear is selected for use in the experiment and is fixed in the test bench with the pin wheel, and its main parameter includes that the reduction ratio is 121, planetary gear figure 2, sun gear number of teeth 12, planet wheel number of teeth 42, cycloid wheel number of teeth 39 and pin wheel number of teeth 40. FIG. 4 is a physical diagram of a planet wheel containing a multi-tooth wear failure. The acquisition system is composed of an NI-USB9234 acquisition card and a one-way acceleration sensor. The vibration signal is acquired by an accelerator sensor placed on the reducer holder, and the sampling frequency is 25.6kHz. Presetting an experiment: the articulated arm is fixed on a supporting disk in the RV reducer 4, and the movement range is 0-90 degrees (single lifting or descending is 90 degrees). The swing arm running speed is 100 DEG/s. The result of each characteristic frequency can be directly calculated according to the parameters as follows: sun gear frequency conversion f 1 39.1Hz, and the planet wheel frequency conversion is f 2 11.09Hz, supporting disk frequency conversion of 0.28Hz and primary transmission meshing frequency f 1c 465.87Hz, planet wheel fault frequency
Is 10.81Hz. The specific diagnosis method comprises the following steps:
s1, firstly, an acceleration sensor 2 is arranged on the upper surface of a reducer retainer 3, an observation signal S (n) of mechanical vibration in the reciprocating intermittent motion process of a robot joint arm 1 is picked up through the acceleration sensor 2, and relevant theoretical characteristic frequency of an RV reducer is calculated; wherein n is the sampling time of the vibration signal s (n);
s2, intercepting the whole-cycle signal of the observation signal S (n) sampled in the step S1 to obtain a periodic motion signal S simultaneously comprising an ascending stage and a descending stage 1 (n) as shown in FIG. 5;
s3, passing formula
The periodic motion signal S obtained in the step S2 1 (n) performing continuous wavelet transform to obtain coefficient matrix W f (a, b) and plotting the time-frequency diagram of the signal, as shown in fig. 6. And intercepting the signal of the stable stage in the ascending motion according to the time-frequency diagram (approximately intercepting the signal between 0.5s and 1 s), thereby obtaining a preprocessed signal s 2 (n) as shown in FIG. 7;
s4, the NeighCoeff algorithm filtering operation is as follows: first on the preprocessed signal s 2 (n) obtaining wavelet coefficients using a discrete wavelet transform. Then according to the formula
And a formula->
The wavelet coefficient is corrected, and then the reconstructed signal is filtered to obtain a filtered signal x (n), as shown in fig. 8;
s5, FFT is performed on the filtered signal x (n) in step S4, to obtain fig. 9. The fault frequency in fig. 9 is 10.94Hz and frequency multiplication, and the actual fault frequency (10.94 Hz) is between the planetary wheel rotation frequency f 2 (11.09 Hz) and theoretical planetary gear fault frequency
(10.81 Hz), the error is between 0Hz and 0.2Hz, and the fault diagnosis can not be accurately finished;
s6, considering that the frequency spectrum can not judge the fault, performing cepstrum analysis on the filtering signal x (n) in the step S5, wherein a cepstrum is shown as 10; two more obvious spectral lines exist simultaneously, the frequency corresponding to the left spectral line is 11.09Hz and represents the frequency conversion of the planet wheel, the frequency corresponding to the right spectral line is 10.94Hz, and the spectral lines with the frequencies of 11.09Hz and 10.94Hz simultaneously appear in the cepstrum, so that the condition that the multi-tooth abrasion fault of the planet wheel occurs can be determined.
Example 3: a fault diagnosis device for an RV reducer, comprising:
the signal pickup module is used for picking up an observation signal s (n) of mechanical vibration in the reciprocating intermittent motion process of the articulated arm 1 through the acceleration sensor 2 and calculating the relevant theoretical characteristic frequency of the RV reducer; wherein n is the sampling time of the vibration signal s (n);
a first obtaining module for intercepting the observation signal s (n) obtained by the signal pickup module to obtain a periodic motion signal s including a rising stage and a falling stage 1 (n)；
An intercepting module for intercepting the periodic motion signal s obtained by the first obtaining module 1 (n) performing continuous wavelet transform, and synchronously intercepting periodic motion signal s according to time-frequency diagram 1 (n) rising phase, approximately stationary vibration signal of sample length m, as preprocessing signal s 2 (n)；
A second obtaining module for obtaining the preprocessed signal s 2 (n) filtering to obtain a filtered signal x (n);
and the fault diagnosis module is used for carrying out FFT on the filtering signal x (n) obtained by the third obtaining module and carrying out fault diagnosis on the frequency domain signal after transformation.
Example 4: a processor for executing a program, wherein the program when executed performs any of the above-described fault diagnosis methods for a RV reducer.
Example 5: a computer-readable storage medium including a stored program, wherein when the program is run, an apparatus in which the computer-readable storage medium is located is controlled to execute any one of the above-described fault diagnosis methods for a RV reducer.
The above-mentioned serial numbers of the embodiments of the present invention are only for description, and do not represent the advantages and disadvantages of the embodiments.
In the above embodiments of the present invention, the descriptions of the respective embodiments have respective emphasis, and for parts that are not described in detail in a certain embodiment, reference may be made to related descriptions of other embodiments.
While the present invention has been described in detail with reference to the embodiments shown in the drawings, the present invention is not limited to the embodiments, and various changes can be made without departing from the spirit of the present invention within the knowledge of those skilled in the art.
Claims (9)
1. A fault diagnosis method for an RV reducer, characterized by: the method comprises the following steps:
s1, an acceleration sensor (2) is used for picking up an observation signal S (n) of mechanical vibration in the reciprocating intermittent motion process of a joint arm (1), and calculating the relevant theoretical characteristic frequency of an RV reducer; wherein n is the sampling time of the vibration signal s (n);
s2, intercepting the whole-cycle signal of the observation signal S (n) obtained in the step S1 to obtain a periodic motion signal S simultaneously comprising an ascending stage and a descending stage 1 (n)；
S3, obtaining the periodic motion signal S obtained in the step S2 1 (n) performing continuous wavelet transform, and synchronously intercepting periodic motion signal s according to time-frequency diagram 1 (n) rising phase, approximately stationary vibration signal of sample length m, as preprocessing signal s 2 (n)；
S4, preprocessing signal S obtained in the step S3 2 (n) filtering to obtain a filtered signal x (n);
s5, performing FFT on the filtering signal x (n) obtained in the step S4, and performing fault diagnosis on the transformed frequency domain signal;
the related theoretical characteristic frequency of the RV reducer comprises a sun gear rotating frequency f 1 Planet wheel frequency conversion f 2 And a supportDisk rotation frequency f 3 Planet gear fault frequency f s r Sun gear and planet gear engagement frequency f 1c 。
2. The fault diagnosis method for the RV reducer according to claim 1, characterized in that: the specific steps for calculating the relevant theoretical characteristic frequency of the RV reducer are as follows:
s1.1, knowing the preset rotation speed n of the articulated arm (1) 3 By calculating the rotational speed n 3 Conversion to a rotary frequency, and the knuckle arm rotary frequency is equivalent to a supporting disc rotary frequency f 3 (ii) a According to the formula
The rotating speed n of the main shaft of the motor (5) is obtained by calculation 1 (ii) a Wherein n is 1 The rotating speed of a main shaft of the motor (5); z 1 Is the number of sun gear teeth, Z 2 Number of teeth of planet gear, Z 4 The number of teeth of the pin gear;
s1.2, by following the formula
The rotation frequency f of the sun wheel is obtained through calculation 1 ；
S1.3, according to the formula
And &>
Calculating the planet wheel rotation frequency f of the RV reducer 2 And the engagement frequency f of the sun wheel and the planet wheel 1c (ii) a Wherein, Z 3 Is the number of cycloid gear teeth;
s1.4, failure frequency f of planet wheel s r The calculation formula is f for the rotation frequency of the planet wheel relative to the planet carrier s r ＝f 2 -f 3 。
3. The fault diagnosis method for the RV reducer according to claim 1, characterized in that: in the continuous wavelet transformation, morlet wavelets are selected as mother wavelets.
4. The fault diagnosis method for an RV reducer according to claim 1, characterized in that: the sample length m is more than or equal to 10F max (ii) a Wherein, F max Representing 2 times the highest frequency value among the relevant theoretical characteristic frequencies of the RV reducer.
5. The fault diagnosis method for an RV reducer according to claim 1, characterized in that: preprocessing the signal s using the NeighCoeff algorithm 2 And (n) filtering.
6. The fault diagnosis method for an RV reducer according to claim 1, characterized in that: the fault diagnosis of the transformed frequency domain signal specifically comprises: carrying out multi-tooth wear fault diagnosis on the planet gear on the frequency domain signals after FFT, and carrying out multi-tooth wear fault diagnosis on the frequency domain signals after FFT when the frequency of 1 doubling of the frequency domain spectral line and the fault frequency f of the planet gear s r When the absolute value of the error is more than 0.2Hz, the multi-tooth abrasion fault of the planet wheel is determined not to occur; when the error between the frequency multiplication of 1 of the frequency domain spectral line and the fault frequency of the planet wheel is between-0.2 Hz and 0Hz or the value is 0, the multi-tooth wear fault of the planet wheel is determined; when the error between the frequency multiplication of the frequency domain spectral line 1 and the planet wheel fault frequency is between 0Hz and 0.2Hz, cepstrum analysis is carried out on the filtering signal x (n), and fault diagnosis is realized.
7. A failure diagnosis device for an RV reducer, characterized in that: the method comprises the following steps:
the signal pickup module is used for picking up an observation signal s (n) of mechanical vibration in the reciprocating intermittent motion process of the articulated arm (1) through the acceleration sensor (2) and calculating the relevant theoretical characteristic frequency of the RV reducer; wherein n is the sampling time of the vibration signal s (n);
a first obtaining module for intercepting the whole cycle signal of the observation signal s (n) obtained by the signal pickup module to obtain a cycle motion signal s including a rising stage and a falling stage 1 (n)；
An intercepting module for intercepting the periodic motion signal s obtained by the first obtaining module 1 (n) performing continuous wavelet transform, and synchronously intercepting periodic motion signal s according to time-frequency diagram 1 (n) a rising phase, a nearly stationary vibration signal of sample length m as a preprocessed signal s 2 (n)；
A second obtaining module for obtaining the preprocessed signal s 2 (n) filtering to obtain a filtered signal x (n);
the fault diagnosis module is used for carrying out FFT on the filtering signal x (n) obtained by the third obtaining module and carrying out fault diagnosis on the frequency domain signal after transformation;
the relevant theoretical characteristic frequency of the RV reducer comprises the revolution frequency f of the sun gear 1 Planetary gear frequency conversion f 2 Support plate frequency conversion f 3 Planet wheel failure frequency f s r Sun gear and planet gear engagement frequency f 1c 。
8. A processor, characterized in that: the processor is configured to execute a program, wherein the program executes the fault diagnosis method for the RV reducer according to any one of claims 1 to 6.
9. A computer-readable storage medium, characterized in that: the computer-readable storage medium includes a stored program, wherein the program, when executed, controls an apparatus in which the computer-readable storage medium is located to perform the fault diagnosis method for the RV reducer according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210385732.1A CN114813117B (en) | 2022-04-13 | 2022-04-13 | Fault diagnosis method and device for RV reducer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210385732.1A CN114813117B (en) | 2022-04-13 | 2022-04-13 | Fault diagnosis method and device for RV reducer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114813117A CN114813117A (en) | 2022-07-29 |
| CN114813117B true CN114813117B (en) | 2023-04-07 |
Family
ID=82537028
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210385732.1A Active CN114813117B (en) | 2022-04-13 | 2022-04-13 | Fault diagnosis method and device for RV reducer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114813117B (en) |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5365595B2 (en) * | 2010-09-15 | 2013-12-11 | 株式会社安川電機 | Speed reducer abnormality determination method, abnormality determination apparatus, robot, and robot system |
| JP5482699B2 (en) * | 2011-03-16 | 2014-05-07 | トヨタ自動車株式会社 | Reduction gear abnormality diagnosis device and reduction gear abnormality diagnosis method |
| CN103018043A (en) * | 2012-11-16 | 2013-04-03 | 东南大学 | Fault diagnosis method of variable-speed bearing |
| CN113405795B (en) * | 2021-06-22 | 2023-12-12 | 昆明理工大学 | Method for identifying weak faults of joint RV reducer |
| CN113639999A (en) * | 2021-08-16 | 2021-11-12 | 北京信息科技大学 | Fault diagnosis method, system, medium and equipment for rolling bearing under variable rotating speed |
-
2022
- 2022-04-13 CN CN202210385732.1A patent/CN114813117B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN114813117A (en) | 2022-07-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113405795B (en) | Method for identifying weak faults of joint RV reducer | |
| CN111975784B (en) | Joint robot fault diagnosis method based on current and vibration signals | |
| CN111665051A (en) | Bearing fault diagnosis method under strong noise variable-speed condition based on energy weight method | |
| Wang et al. | Multi-scale enveloping order spectrogram for rotating machine health diagnosis | |
| Cheng et al. | Envelope deformation in computed order tracking and error in order analysis | |
| CN109520738B (en) | Rotating machinery rolling bearing fault diagnosis method based on order spectrum and envelope spectrum | |
| CN106959211A (en) | The fatigue test device and its method of testing of a kind of gear engagement | |
| CN110487547B (en) | Rolling bearing fault diagnosis method under variable working conditions based on vibration diagram and transfer learning | |
| Jaber et al. | Industrial robot backlash fault diagnosis based on discrete wavelet transform and artificial neural network | |
| CN108015797B (en) | A kind of RV speed reducer drive error on-line monitoring method | |
| CN111879508B (en) | Method and device for estimating instantaneous rotating speed of rotating machine based on time-frequency transformation and storage medium | |
| CN106124197A (en) | A kind of epicyclic gearbox sun gear partial fault detection method and system | |
| CN108036938A (en) | A kind of industrial robot gear reducer vibration mechine and test method | |
| CN113702043A (en) | Planetary gear box fault diagnosis method based on POVMD and FDTW under time-varying rotation speed | |
| Shang et al. | An intelligent fault diagnosis system for newly assembled transmission | |
| CN114813117B (en) | Fault diagnosis method and device for RV reducer | |
| WO2020255446A1 (en) | Operation sound diagnosis system, operation sound diagnosis method, and machine learning device for operation sound diagnosis system | |
| CN111413095A (en) | Instantaneous angular velocity-based distributed fault diagnosis and analysis method for planetary bearing | |
| CN112327957B (en) | Method and device for controlling low-frequency vibration multi-order line spectrum of rotor | |
| CN112287492B (en) | Dynamic characteristic identification method for five-axis linkage machining center swinging turntable | |
| CN218885427U (en) | Gear fault simulation experiment table | |
| CN113340631A (en) | Torsional vibration testing device and signal analysis method | |
| CN108106840A (en) | A kind of Novel industrial robot RV speed reducer method for testing performances | |
| CN115421400A (en) | Digital twin virtual-real consistency discrimination method for service health state intelligent recognition decision | |
| Zhou et al. | Current-based feed axis condition monitoring and fault diagnosis |
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 |