CN114528525B - Mechanical fault diagnosis method based on maximum weighted kurtosis blind deconvolution - Google Patents
Mechanical fault diagnosis method based on maximum weighted kurtosis blind deconvolution Download PDFInfo
- Publication number
- CN114528525B CN114528525B CN202210028718.6A CN202210028718A CN114528525B CN 114528525 B CN114528525 B CN 114528525B CN 202210028718 A CN202210028718 A CN 202210028718A CN 114528525 B CN114528525 B CN 114528525B
- Authority
- CN
- China
- Prior art keywords
- signal
- kurtosis
- filter
- blind deconvolution
- weighted kurtosis
- 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
- 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/15—Correlation function computation including computation of convolution operations
-
- 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/021—Gearings
-
- 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
-
- 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/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Computing Systems (AREA)
- Acoustics & Sound (AREA)
- Wind Motors (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention relates to a mechanical fault diagnosis method based on maximum weighted kurtosis blind deconvolution, belongs to the technical field of wind turbine generator fault diagnosis, and provides a new blind deconvolution method, namely maximum weighted kurtosis blind deconvolution. The heavily weighted kurtosis has good robustness to single or a small amount of strong impact interference in a fault signal, and prior knowledge of a fault impact sequence to be recovered is not needed. Based on the method, the maximum weighted kurtosis blind deconvolution method can effectively solve the problem that a classical kurtosis maximization-based method tends to recover a single leading impact instead of a gear fault impact sequence, and meanwhile, compared with a common non-full 'blind' (relying on fault characteristic frequency prior), the method has stronger applicability in the aspect of industrial equipment gear fault diagnosis. And the application case in wind turbine fault diagnosis proves the effectiveness of the method for gear fault diagnosis.
Description
Technical Field
The invention belongs to the technical field of wind turbine generator system fault diagnosis, and particularly relates to a mechanical fault diagnosis method based on maximum weighted kurtosis blind deconvolution.
Background
The wind turbine is the main equipment of the wind power plant, the price accounts for 74-82% of the total investment of the wind power plant, and the maintenance cost of the wind turbine becomes the main operation cost of the wind power plant due to high equipment failure rate and high maintenance cost caused by severe operation environment. Reducing maintenance cost of the wind turbine is an important way for improving economic benefits of wind power plant operation, and in order to effectively reduce maintenance cost of the wind turbine, wind power enterprises introduce technologies such as state monitoring, fault diagnosis and state maintenance in a dispute.
In a fault diagnosis link, a gear fault characteristic signal is usually weak under the influence of noise, a transmission path and the like, and the current fault diagnosis cannot effectively diagnose the gear fault of the wind turbine generator.
Therefore, at the present stage, a mechanical fault diagnosis method based on blind deconvolution of the maximum weighted kurtosis needs to be designed to solve the above problems.
Disclosure of Invention
The invention aims to provide a mechanical fault diagnosis method based on maximum weighted kurtosis blind deconvolution, which is used for solving the technical problems in the prior art, wherein in the fault diagnosis link, a gear fault characteristic signal is usually weak under the influence of noise, a transmission path and the like, and the current fault diagnosis cannot effectively diagnose the gear fault of a wind turbine generator.
In order to achieve the purpose, the technical scheme of the invention is as follows:
the mechanical fault diagnosis method based on the maximum weighted kurtosis blind deconvolution comprises the following steps:
s1: inputting a measuring signal x, and randomly initializing a filter h coefficient;
s2: filtering the measurement signal x through a filter h to obtain a filtering signal s, and calculating the weighted kurtosis RK(s) of the filtering signal s;
s3: continuously updating the h coefficient of the filter by iterative solution of blind deconvolution of maximum weighted kurtosis;
s4: repeating the steps S2 and S3 to enable the weighted kurtosis RK (S) to reach the maximum and meet the iteration times;
s5: and the filtering signal s corresponding to the maximum weighted kurtosis RK(s) is the target filtering signal.
Further, in step S2, the step of calculating the weighted kurtosis RK (S) of the filtered signal S is as follows:
(1) M equal division is carried out on the filtering signal s to obtain each sub-section signal s m (M =1, …, M), M being any positive integer not exceeding the signal length;
(2) Calculating Kurt of each sub-section signal m ;
(3) For Kurt m The ascending order is arranged and expressed as a vector, namely: kurt is a very powerful tool asc ;
(4) Calculate Kurt m The weight of the sum, i.e.:
(5) To W m The descending order is performed, which is also expressed in vector form, i.e.: w is a group of desc ;
(6) Using rearranged weight directionQuantity W desc For rearranged Kurt vector asc Weighting to obtain the weighted kurtosis RK of the signal:
RK=Kurt asc ·(W desc ) T (2)。
further, M is 4.
Further, step S3 specifically includes the following steps:
the maximum heavily weighted kurtosis blind deconvolution solved filter coefficients are represented as:
to obtain the maximum value of RK, the partial derivative of RK on the filter coefficients is made equal to 0, i.e.:
first of all, the first step is to,
wherein Kurt is θ For the kurtosis of theta stage filtering signal, the effective approximate solution of equation is obtained by iterative solution and mode of continuously updating filter coefficient, and the equation is utilizedObtaining M filters, and weighting and summing the M filters to obtain an updated filter, wherein the process is as follows:
where N is the filter signal length, L is the filter length, x θ ,s θ Respectively corresponding to the theta-th section of the signal subsection and the filtered signal; thus, an updated filter is obtained:
a storage medium having stored thereon a computer program which, when executed, performs a method of mechanical fault diagnosis based on blind deconvolution of maximum heavily weighted kurtosis as described above.
An electronic device comprises a processor and a memory, wherein the memory is used for storing executable commands of the processor, and the processor executes the executable commands to realize the mechanical fault diagnosis method based on the maximum weighted kurtosis blind deconvolution.
Compared with the prior art, the invention has the beneficial effects that:
one of the benefits of the present solution is to provide a new blind deconvolution method, i.e., blind deconvolution with maximum weighted kurtosis. The heavily weighted kurtosis has good robustness to single or a small amount of strong impact interference in a fault signal, and prior knowledge of a fault impact sequence to be recovered is not needed. Based on the method, the maximum weighted kurtosis blind deconvolution method can effectively solve the problem that a classical kurtosis maximization-based method tends to recover a single leading impact instead of a gear fault impact sequence, and meanwhile, compared with a common non-full 'blind' (relying on fault characteristic frequency prior), the method has stronger applicability in the aspect of industrial equipment gear fault diagnosis. And the application case in wind turbine fault diagnosis proves the effectiveness of the method for gear fault diagnosis.
Drawings
FIG. 1 is a schematic diagram of a measurement signal and a filtering signal of a gearbox of a wind power transmission system according to an embodiment of the application.
Fig. 2 is a schematic diagram illustrating a process of calculating a weighted kurtosis of a filtered signal according to an embodiment of the present disclosure.
Fig. 3 is a schematic diagram of a method according to the present embodiment.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to fig. 1 to 3 of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example (b):
as shown in fig. 2, the process of calculating the weighted kurtosis of the filtered signal is as follows:
(1) M equal division is carried out on the filtering signal s to obtain each sub-section signal s m (m=1,…,M)。
(2) Calculating Kurt of each sub-section signal m 。
(3) For Kurt m The ascending order is arranged and expressed as a vector, namely: kurt is a very powerful tool asc 。
(4) Calculate Kurt m The weight of the sum, i.e.:
(5) To W m The descending order is performed, which is also expressed in vector form, i.e.: w is a group of desc 。
(6) Using the rearranged weight vector W desc For rearranged Kurt vector asc Weighting to obtain the weighted kurtosis RK of the signal:
RK=Kurt asc ·(W desc ) T (2)
the value of M does not depend on a specific formula, theoretically, any positive integer not exceeding the signal length can enable the algorithm to normally run, but each sub-integerThe length of the segment signal cannot be too small, otherwise the Kurt of the sub-segment signal m The physical significance will be lost. Analysis of multiple sets of simulation signals and actual gear vibration signals shows that most of data analysis can achieve ideal results when M =4, so that the value of M is recommended to be 4 when the method is used.
Based on this, the blind deconvolution of the maximum heavily weighted kurtosis to solve FIR filter coefficients can be represented as a maximization problem as follows:
to find the maximum value of RK, let RK equal 0 for the partial derivative of the filter coefficients, i.e.:
first of all, the first step is to,
wherein Kurt is θ For the kurtosis of the filtered signal in section θ, it can be seen that equation (4) is highly nonlinear and difficult to solve directly. Therefore, an effective approximate solution to the equation is solved by iteratively solving and continuously updating the filter coefficients. In addition, to simplify the operation, use is made ofObtaining M filters, and weighting and summing the M filters to obtain an updated filter, wherein the process comprises the following steps:
where N is the filter signal length, L is the filter length, x θ ,s θ Respectively corresponding to the theta-th section of the signal subsection and the filtered signal. Thus, an updated filter is obtained:
as shown in fig. 3, in summary, the overall flow of the proposed method is as follows:
(1) Input measurement signal x, specify parameters M =4, randomly initialize filter h coefficients, such as:
(2) And filtering the measurement signal x through a filter h to obtain a filtering signal s, and calculating the weighted kurtosis RK(s) of the filtering signal according to a formula (2).
(3) The filter h coefficients are updated according to equation (9).
(4) Repeating steps 2 and 3 maximizes RK(s) and satisfies the number of iterations.
(5) And 4, the filtering signal corresponding to the maximum RK(s) in the step 4 is the target filtering signal.
The method is researched and verified through the wind power transmission system fault diagnosis application.
This data comes from vibration monitoring of a wind turbine drive train. In the signal acquisition process, the sampling frequency is 97656Hz, and the rotating speed of the output pinion shaft of the gearbox is 1800rpm measured by the tachometer. Intercepting the signal length 51200 points, as shown in fig. 1 (a), the gear fault information is completely masked.
The filtered signal of the method is shown in fig. 1 (b). The frequency of the periodic impact sequence recovered by the method is 29.6Hz (1/0.0338 Hz), and is very close to the theoretical fault characteristic frequency (30 Hz) of the pinion. Therefore, it can be preliminarily considered that the pinion gear has a failure. And a certain gear tooth of the pinion is found to be damaged in the subsequent shutdown and unpacking detection, so that the effectiveness of the method is verified.
In summary, the present solution provides a new blind deconvolution method, i.e., a blind deconvolution with the maximum weighted kurtosis. The heavily weighted kurtosis has good robustness to single or a small amount of strong impact interference in a fault signal, and prior knowledge of a fault impact sequence to be recovered is not needed. Based on the method, the maximum weighted kurtosis blind deconvolution method can effectively solve the problem that a classical kurtosis maximization-based method tends to recover a single leading impact instead of a gear fault impact sequence, and meanwhile, compared with a common non-full 'blind' (relying on fault characteristic frequency prior), the method has stronger applicability in the aspect of industrial equipment gear fault diagnosis. And the application case in wind turbine fault diagnosis proves the effectiveness of the method for gear fault diagnosis.
The above are preferred embodiments of the present invention, and all changes made according to the technical scheme of the present invention that produce functional effects do not exceed the scope of the technical scheme of the present invention belong to the protection scope of the present invention.
Claims (4)
1. The mechanical fault diagnosis method based on the maximum weighted kurtosis blind deconvolution is characterized by comprising the following steps of:
s1: inputting a measuring signal x, and randomly initializing a filter h coefficient;
s2: filtering the measurement signal x through a filter h to obtain a filtering signal s, and calculating the weighted kurtosis RK(s) of the filtering signal s;
s3: continuously updating the h coefficient of the filter by iterative solution of blind deconvolution of the maximum weighted kurtosis;
s4: repeating the steps S2 and S3 to enable the weighted kurtosis RK (S) to reach the maximum and meet the iteration number;
s5: the filtering signal s corresponding to the maximum weighted kurtosis RK(s) is the target filtering signal;
in step S2, the step of calculating the weighted kurtosis RK (S) of the filtered signal S is as follows:
(1) M is equally divided to the filtering signal s to obtain eachSub-section signal s m Wherein M =1, …, M is any positive integer not exceeding the signal length;
(2) Calculating Kurt of each sub-section signal m ;
(3) For Kurt m The ascending order is arranged and expressed as a vector, namely: kurt is a very powerful tool asc ;
(4) Calculate Kurt m The weight of the sum, i.e.:
(5) To W m The descending order is performed, which is also expressed in vector form, i.e.: w desc ;
(6) Using rearranged weight vector W desc For rearranged kurtosis vector Kurt asc And weighting to obtain the weighted kurtosis RK of the signal:
RK=Kurt asc ·(W desc ) T (2);
the step S3 is specifically as follows:
the maximum weighted kurtosis blind deconvolution solved filter coefficient is expressed as:
to find the maximum value of RK, let RK equal 0 for the partial derivative of the filter coefficients, i.e.:
first of all, when a user wants to use the apparatus,
wherein Kurt is θ For the kurtosis, pass, of the theta-th filtered signalSolving the effective approximate solution of the equation by over-iterative solution and continuously updating the filter coefficient, and utilizingObtaining M filters, and weighting and summing the M filters to obtain an updated filter, wherein the process is as follows:
where N is the filter signal length, L is the filter length, x θ ,s θ Respectively corresponding to the theta-th section of the signal subsection and the filtered signal; thus, an updated filter is obtained:
2. the method of claim 1, wherein M is 4.
3. A storage medium having stored thereon a computer program which, when executed, performs a method of diagnosing a mechanical failure based on blind deconvolution of maximum weighted kurtosis as claimed in claim 1 or 2.
4. An electronic device, comprising a processor and a memory, wherein the memory is used for storing executable commands of the processor, and the processor executes the executable commands to realize the method for diagnosing the mechanical fault based on the blind deconvolution of the maximum weighted kurtosis as claimed in claim 1 or 2.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210028718.6A CN114528525B (en) | 2022-01-11 | 2022-01-11 | Mechanical fault diagnosis method based on maximum weighted kurtosis blind deconvolution |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210028718.6A CN114528525B (en) | 2022-01-11 | 2022-01-11 | Mechanical fault diagnosis method based on maximum weighted kurtosis blind deconvolution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114528525A CN114528525A (en) | 2022-05-24 |
CN114528525B true CN114528525B (en) | 2023-03-28 |
Family
ID=81621607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210028718.6A Active CN114528525B (en) | 2022-01-11 | 2022-01-11 | Mechanical fault diagnosis method based on maximum weighted kurtosis blind deconvolution |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114528525B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115859091B (en) * | 2022-11-01 | 2023-05-26 | 哈尔滨工业大学 | Bearing fault feature extraction method, electronic equipment and storage medium |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108536877A (en) * | 2017-03-01 | 2018-09-14 | 新疆金风科技股份有限公司 | The appraisal procedure and system of vibration of wind generating set degree |
CN111623982A (en) * | 2020-06-15 | 2020-09-04 | 大连理工大学 | Planetary gearbox early fault diagnosis method based on APEWT and IMOMEDA |
CN113536226A (en) * | 2021-07-14 | 2021-10-22 | 东南大学 | Blind deconvolution algorithm for enhancing rotating mechanical fault signal characteristics |
CN113702042A (en) * | 2021-07-14 | 2021-11-26 | 西南交通大学 | Mechanical fault diagnosis method and system based on enhanced minimum entropy deconvolution |
CN113702030A (en) * | 2021-07-14 | 2021-11-26 | 西南交通大学 | Mechanical fault diagnosis method and system based on maximum average window kurtosis blind deconvolution |
CN113702037A (en) * | 2021-07-14 | 2021-11-26 | 西南交通大学 | Re-weighted spectrum kurtosis method based on subband rearrangement and set dual-tree complex wavelet packet transformation |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8175291B2 (en) * | 2007-12-19 | 2012-05-08 | Qualcomm Incorporated | Systems, methods, and apparatus for multi-microphone based speech enhancement |
-
2022
- 2022-01-11 CN CN202210028718.6A patent/CN114528525B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108536877A (en) * | 2017-03-01 | 2018-09-14 | 新疆金风科技股份有限公司 | The appraisal procedure and system of vibration of wind generating set degree |
CN111623982A (en) * | 2020-06-15 | 2020-09-04 | 大连理工大学 | Planetary gearbox early fault diagnosis method based on APEWT and IMOMEDA |
CN113536226A (en) * | 2021-07-14 | 2021-10-22 | 东南大学 | Blind deconvolution algorithm for enhancing rotating mechanical fault signal characteristics |
CN113702042A (en) * | 2021-07-14 | 2021-11-26 | 西南交通大学 | Mechanical fault diagnosis method and system based on enhanced minimum entropy deconvolution |
CN113702030A (en) * | 2021-07-14 | 2021-11-26 | 西南交通大学 | Mechanical fault diagnosis method and system based on maximum average window kurtosis blind deconvolution |
CN113702037A (en) * | 2021-07-14 | 2021-11-26 | 西南交通大学 | Re-weighted spectrum kurtosis method based on subband rearrangement and set dual-tree complex wavelet packet transformation |
Non-Patent Citations (5)
Title |
---|
Kaixuan Liang 等.Maximum average kurtosis deconvolution and its application for the impulsive fault feature enhancement of rotating machinery.《Mechanical Systems and Signal Processing》.2021,第149卷第1-16页. * |
Lei Wu 等.Maximum reweighted-kurtosis deconvolution:a fully blind and adaptive method for restoration of gear fault impulse trains.《Measurement Science and Technology》.2022,第33卷(第33期),第1-12页. * |
刘文朋 ; 廖英英 ; 杨绍普 ; 刘永强 ; 顾晓辉 ; .一种基于多点峭度谱和最大相关峭度解卷积的滚动轴承故障诊断方法.振动与冲击.2019,(第02期),第146-151页. * |
吴磊 等.基于最大重加权峭度盲解卷积的风电故障诊断.《中国机械工程》.2022,第33卷(第33期),第2356-2363页. * |
吴磊 等.基于最大重加权峭度盲解卷积的风电机组齿轮故障诊断 .《中国机械工程》.2022,第1-8页. * |
Also Published As
Publication number | Publication date |
---|---|
CN114528525A (en) | 2022-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Research and application of improved adaptive MOMEDA fault diagnosis method | |
Liu et al. | A hybrid fault diagnosis method based on second generation wavelet de-noising and local mean decomposition for rotating machinery | |
Zheng et al. | Sparse elitist group lasso denoising in frequency domain for bearing fault diagnosis | |
Cai et al. | Sparsity-enhanced signal decomposition via generalized minimax-concave penalty for gearbox fault diagnosis | |
CN103645052B (en) | A kind of gearbox of wind turbine remote online status monitoring and lifetime estimation method | |
Feng et al. | A novel order spectrum-based Vold-Kalman filter bandwidth selection scheme for fault diagnosis of gearbox in offshore wind turbines | |
CN114528525B (en) | Mechanical fault diagnosis method based on maximum weighted kurtosis blind deconvolution | |
Liu et al. | Naturally damaged wind turbine blade bearing fault detection using novel iterative nonlinear filter and morphological analysis | |
CN111507221B (en) | Gear signal denoising method based on VMD and maximum overlapping discrete wavelet packet transformation | |
He et al. | Feature extraction of gear and bearing compound faults based on vibration signal sparse decomposition | |
Wang et al. | Multiwavelet construction via an adaptive symmetric lifting scheme and its applications for rotating machinery fault diagnosis | |
CN105841792A (en) | Gear pressure angle direction local vibration signal acquisition method based on microsensor | |
CN113390631A (en) | Fault diagnosis method for gearbox of diesel engine | |
CN113250911A (en) | Fan blade fault diagnosis method based on VMD decomposition algorithm | |
CN113702030A (en) | Mechanical fault diagnosis method and system based on maximum average window kurtosis blind deconvolution | |
CN111582248A (en) | SVD-based gearbox signal noise reduction method | |
CN107121493A (en) | A kind of equipment damage detection method and system | |
CN114166503B (en) | Gear health monitoring index construction method based on order sideband product spectrum | |
CN112484998B (en) | Wind turbine generator bearing fault diagnosis method based on synchronous modal spectrum | |
Yu | Early fault detection for gear shaft and planetary gear based on wavelet and hidden Markov modeling | |
CN112161807A (en) | Fault diagnosis method and device for speed change gear box and storage medium | |
CN116610939B (en) | Asymmetric penalty sparse regularization pulse extraction method under Coiflet discrete wavelet | |
CN113465916B (en) | Gear tooth state evaluation method, device, equipment and medium of planetary gear train | |
Zhou et al. | Application of morphological filter in purification of rotor center's orbit | |
CN113358212B (en) | Electromechanical fault diagnosis method and system based on relative harmonic order and modeling method |
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 |