CN111975784A - Joint robot fault diagnosis method based on current and vibration signals - Google Patents
Joint robot fault diagnosis method based on current and vibration signals Download PDFInfo
- Publication number
- CN111975784A CN111975784A CN202010918331.9A CN202010918331A CN111975784A CN 111975784 A CN111975784 A CN 111975784A CN 202010918331 A CN202010918331 A CN 202010918331A CN 111975784 A CN111975784 A CN 111975784A
- Authority
- CN
- China
- Prior art keywords
- joint
- vibration
- robot
- sequence
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000003745 diagnosis Methods 0.000 title claims abstract description 15
- 238000001228 spectrum Methods 0.000 claims abstract description 18
- 238000005070 sampling Methods 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 238000001514 detection method Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000006641 stabilisation Effects 0.000 abstract description 2
- 238000011105 stabilization Methods 0.000 abstract description 2
- 230000001360 synchronised effect Effects 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000000605 extraction Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 239000003638 chemical reducing agent Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000010183 spectrum analysis Methods 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J19/00—Accessories fitted to manipulators, e.g. for monitoring, for viewing; Safety devices combined with or specially adapted for use in connection with manipulators
- B25J19/0095—Means or methods for testing manipulators
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1674—Programme controls characterised by safety, monitoring, diagnostic
Landscapes
- Engineering & Computer Science (AREA)
- Robotics (AREA)
- Mechanical Engineering (AREA)
- Manipulator (AREA)
Abstract
The invention discloses a joint robot fault diagnosis method based on current and vibration signals, which comprises the steps of firstly picking up a motor current signal and a joint vibration signal of a joint to be researched of a robot; filtering the motor current signal to extract the instantaneous frequency of the joint motor current signal; calculating a real-time rotation angle of a joint motor shaft by using the instantaneous frequency extracted from the joint motor current signal; carrying out low-pass filtering processing on the vibration signal; determining the maximum order ratio range and equally dividing the rotation angle to obtain an equal-angle sampling time sequence and carrying out equal-angle sampling on the filtered vibration signal; finally, Fourier transformation is carried out on the equal-angle sampling sequence of the vibration signals to obtain a robot joint vibration signal order spectrum and the robot joint vibration signal order spectrum is analyzed to realize robot joint fault diagnosis; the first-order frequency required by vibration order ratio analysis is extracted very conveniently and quickly, and vibration signal stabilization processing can be completed only by ensuring synchronous acquisition of vibration and current, so that robot joint fault diagnosis and state detection are completed.
Description
Technical Field
The invention relates to a joint robot fault diagnosis method based on current and vibration signals, and belongs to the technical field of state monitoring and fault diagnosis of industrial robots.
Background
The six-degree-of-freedom serial industrial robot has very wide application in industrial automation production, and parts abrasion caused by long-time operation and the occurrence of some emergencies can cause the robot to stop operating suddenly, so that the operation of a production line is damaged, and loss is caused to enterprises. The main devices of the robot joint are a servo motor and a reducer, and the monitoring of the robot and the traditional motor and reducer state monitoring method can be used for reference. The traditional motor state monitoring adopts the vibration signal spectrum analysis; for a six-degree-of-freedom series robot, a vibration signal is a typical non-stationary signal, and is generally firstly subjected to stationary processing, and the stationary processing mode generally adopts order ratio tracking. However, for the robot, the key phase pulse is difficult to obtain, the joint reducer has a complex structure, and the time-frequency instantaneous frequency extraction difficulty is high, so that the order ratio analysis from the two aspects has a plurality of problems, and the instantaneous frequency extracted by the current signals acquired synchronously replaces the time-frequency analysis to obtain the instantaneous conversion frequency relatively easily.
When the robot operates, the rotating speed of the robot is always in a changing state, and joint vibration signals of the robot are typical non-stationary signals; the robot joint structure is complex, so it is not practical to realize order tracking through hardware. Meanwhile, the joint servo motor encoder feedback pulse acquisition difficulty is high, and the cost is high. In addition, because the robot joint comprises a motor, a speed reducer and other parts, and each joint is in a connection state, the frequency conversion of the motor shaft in a frequency spectrum is not obvious when the robot joint vibrates a signal, the extraction difficulty is high, an equiangular sampling sequence is difficult to obtain, and the order ratio analysis of the vibration signal is difficult. The robot joint motor current signal acquisition cost is low, the acquisition is convenient, and the order analysis by using the current signal to extract the instantaneous frequency has great significance.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a joint robot fault diagnosis method based on current and vibration signals, which does not need key phase signals and extraction of instantaneous frequency of the vibration signals, and can complete the order analysis of the vibration signals only by synchronously acquiring the current and the vibration signals and extracting the instantaneous frequency of the current signals.
The vibration signal order ratio analysis method based on the current comprises the following specific steps:
(1) synchronously acquiring a motor current signal I (t) and a vibration signal S (t) of an ith joint of the robot in the process of rotating any angle, wherein t is sampling time;
acquiring a motor current signal I (t) of an ith joint of the robot through a current sensor, and acquiring a vibration signal S (t) of the ith joint of the robot through an acceleration sensor;
the arbitrary angle refers to an arbitrary value within an angle range allowed by a joint, and the angle range is described on a robot specification or a robot demonstrator;
(2) filtering and deburring the current signal I (t) by adopting a method of combining zero-phase filtering with singular value denoising to obtain a filtered current signal I1(t);
(3) By the formulaCalculating the current signal I1(t) phase phi (t) by the formulaCalculating the current signal I1(t) instantaneous frequency f (t);
(4) by the formulaWherein p is the number of pole pairs of the motor, the instantaneous rotating speed n (t) of the motor shaft is obtained through calculation, and the instantaneous rotating speed n (t) is obtained through a formulaCalculating to obtain an instantaneous rotation angle theta (t) of the motor shaft;
(5) by the formulaCalculating an equal-angle sampling interval delta theta; maximum order ratio vmaxIs an integer greater than or equal to 1;
(6) sampling the instantaneous rotation angle theta (t) according to the equal-angle sampling interval delta theta obtained by calculation in the step (5) to obtain an equal-angle instantaneous rotation angle sequence theta1(t);
(7) Cubic spline interpolation is carried out on the instantaneous corner theta (t) to obtain an equal-angle instantaneous corner sequence theta1(T) time series T;
(8) carrying out low-pass filtering on the vibration signal S (t) acquired in the step (1) to obtain a filtered signal S1(t) for S1(T) carrying out cubic spline interpolation to obtain a vibration sequence S of the time sequence T2(t);
(9) By the formulaMaximum order ratio vmaxIs an integer greater than or equal to 1; calculating a sequence of order ratios L, where n is the vibration sequence S2(t) length;
(10) for vibration sequence S2(t) performing Fast Fourier Transform (FFT) to obtain a ratio spectrum sequence S3(t);
Taking the order ratio sequence L as an abscissa, and taking the order ratio spectrum sequence S3(t) drawing a scale spectrum with the first m values as vertical coordinates, wherein the m value is half of the length of the scale sequence L, and if L is an odd number, rounding the half of the length of the scale sequence L downwards;
and observing whether fault characteristic order ratios of all parts of the robot joint exist in order ratio components in the order ratio spectrum, if so, indicating that the parts have faults, and if not, indicating that the joints have no faults, thereby completing vibration signal order ratio analysis and realizing robot joint fault diagnosis and state detection.
The invention has the beneficial effects that:
1. the current signal instantaneous frequency extraction approach is superior to that of key phase pulse calculation and extraction from vibration signal time frequency spectrum;
2. the method can conveniently and quickly realize the order ratio analysis of the vibration signals;
3. the first-order frequency required by vibration order ratio analysis is extracted very conveniently and quickly, and vibration signal stabilization processing can be completed only by ensuring synchronous acquisition of vibration and current, so that joint state detection and fault diagnosis of the robot are realized.
Drawings
FIG. 1 is a schematic view of a robot configuration;
FIG. 2 is a schematic diagram of the original waveform of the current signal I (t) in the state of 60 degrees of articulation 2;
fig. 3 is a schematic diagram of an original waveform of a vibration signal s (t) in a state of 60 ° of joint movement of fig. 2;
FIG. 4 is a current signal I (t) after filtering for the current signal I (t) at 60 degrees of articulation 21(t) a waveform diagram;
FIG. 5 is a current signal I (t) after filtering for a current signal I (t) at 60 degrees of articulation 21(t) a waveform of the instantaneous frequency f (t);
fig. 6 is a waveform diagram of an instantaneous rotation angle θ (t) calculated by an instantaneous frequency f (t);
FIG. 7 is a sequence S after resampling of the vibration signal2(t) a waveform diagram;
FIG. 8 is a schematic diagram of a vibration signal scale spectrum;
in fig. 1: the system comprises a 1-1 st joint, a 2-2 nd joint, a 3-connecting arm I, a 4-3 rd joint, a 5-4 th joint, a 6-connecting arm II, a 7-5 th joint, an 8-6 th joint, a 9-electric cabinet, a 10-current sensor and an 11-acceleration sensor.
Detailed Description
The present invention is further illustrated by the following examples, but the scope of the invention is not limited to the above-described examples.
Example 1: the vibration signal order ratio analysis method based on the current comprises the following steps:
the Qianjiang QJR6-1 welding robot is adopted, and fig. 1 is a schematic structural diagram of the robot, which comprises a joint I1, a joint II 2, a connecting arm I3, a joint III 4, a joint IV 5, a connecting arm II 6, a joint V7, a joint VI 8, an electric cabinet 9, a current sensor 10 and an acceleration sensor 11; the implementation object is a joint II, the motion angle is 60 degrees, and the specific operation flow is as follows:
1. in the state that a single joint of a second joint moves by 60 degrees, a current sensor is adopted to obtain a motor current signal I (t) of the second joint of the robot, as shown in figure 2, an acceleration sensor 11 is adopted to obtain a vibration signal S (t) of the second joint, as shown in figure 3, t is sampling time, an acquisition card is NI9234, acquisition software is SignalExpress, and sampling frequency is 25600 Hz;
2. the method for realizing zero-phase filtering and singular value denoising by utilizing matlab software is used for carrying out filtering and deburring processing on the current signal I (t) to obtain a filtered current signal I1(t); wherein the zero phase filter parameters are as follows: chebyshev type 1, 5-order, 0.2 cutoff low pass, 0.1 equal ripple; a time-frequency matrix used in the singular value denoising process is obtained through short-time Fourier transform, the length of a short-time Fourier transform window is 1024, and the step length is 128; the results are shown in FIG. 4;
the following calculation processes and output results are all completed in matlab software;
3. by the formulaCalculating the current signal I1(t) phase phi (t) by the formulaCalculating the current signal I1(t) instantaneous frequency f (t), the results are shown in FIG. 5;
4. by the formulaWherein the number p of pole pairs of the motor is equal to 5, the instantaneous rotating speed n (t) of the motor shaft is obtained through calculation, and the instantaneous rotating speed n (t) is calculated through a formulaThe instantaneous rotation angle θ (t) is calculated, and the result is shown in fig. 6;
5. taking the maximum order ratio vmaxEqual to 30, by formulaCalculating equal-angle sampling interval delta theta to be equal to 6; sampling the instantaneous rotation angle theta (t) according to the equal-angle sampling interval of 6 degrees to obtain an equal-angle instantaneous rotation angle sequence theta1(t); cubic spline interpolation is carried out on the instantaneous turning angle theta (t) by utilizing matlab software to obtain the equal-angle instantaneous turning angle theta1A time series T at (T);
6. calculating the frequency f corresponding to the highest rotation speed of the motor1Namely: converting the unit r/min of the rotating speed into degree/s, and realizing the cutoff frequency f by utilizing matlab software1The low-pass filter filters the vibration signal S (t) to obtain a filtered vibration signal S1(t); realization of S pair by utilizing matlab software1(T) carrying out cubic spline interpolation to obtain a vibration sequence S at the time sequence T2(t), as shown in FIG. 7;
7. by the formulaCalculating a sequence of order ratios L, where n is the vibration sequence S2(t) has a length equal to 817;
8. vibration sequence S by utilizing matlab software2(t) performing Fast Fourier Transform (FFT) to obtain a ratio spectrum sequence S3(t); taking the order ratio sequence L as an abscissa, and taking the order ratio spectrum sequence S3(t) (m is half of the length of the sequence L equal to 408.5, and 408)) as the ordinate, and a scale spectrum is plotted, and as shown in FIG. 8, the extracted scale spectrum is used for spectrum analysis, and the scale spectrum of 1, 2, 3, 4, and 16 orders are observedAnd the specific components respectively correspond to the order corresponding to the conversion frequency of the motor shaft of the second joint, the order corresponding to the 2, 3 and 4 frequency doubling of the conversion frequency and the order corresponding to the meshing frequency of the sun wheel and the planet wheel, and other obvious orders do not appear in addition, so that the second joint is judged to have no fault. The method is used for carrying out order ratio analysis on different joints to realize state monitoring and early fault diagnosis on each joint arm of the robot, and the condition that the robot is suddenly stopped due to the expiration of the service life of parts or some emergencies to damage a production line and influence production is prevented.
Claims (3)
1. A joint robot fault diagnosis method based on current and vibration signals is characterized by comprising the following specific steps:
(1) synchronously acquiring a motor current signal I (t) and a vibration signal S (t) of an ith joint of the robot in the process of rotating any angle, wherein t is sampling time;
(2) filtering and deburring the current signal I (t) by adopting a method of combining zero-phase filtering with singular value denoising to obtain a filtered current signal I1(t);
(3) By the formulaCalculating the current signal I1(t) phase phi (t) by the formulaCalculating the current signal I1(t) instantaneous frequency f (t);
(4) by the formulaWherein p is the number of pole pairs of the motor, the instantaneous rotating speed n (t) of the motor shaft is obtained through calculation, and the instantaneous rotating speed n (t) is obtained through a formulaCalculating to obtain an instantaneous rotation angle theta (t) of the motor shaft;
(5) by passingFormula (II)Calculating an equal-angle sampling interval delta theta; maximum order ratio vmaxIs an integer greater than or equal to 1;
(6) sampling the instantaneous rotation angle theta (t) according to the equal-angle sampling interval delta theta obtained by calculation in the step (5) to obtain an equal-angle instantaneous rotation angle sequence theta1(t);
(7) Cubic spline interpolation is carried out on the instantaneous corner theta (t) to obtain an equal-angle instantaneous corner sequence theta1(T) time series T;
(8) carrying out low-pass filtering on the vibration signal S (t) acquired in the step (1) to obtain a filtered signal S1(t) for S1(T) carrying out cubic spline interpolation to obtain a vibration sequence S of the time sequence T2(t);
(9) By the formulaWherein N ═ 0, 1, 2, 3,. and N-1, maximum order ratio vmaxIs an integer greater than or equal to 1, calculating a sequence of order ratios L, where n is a vibration sequence S2(t) length;
(10) for vibration sequence S2(t) performing fast Fourier transform to obtain a ratio spectrum sequence S3(t);
Taking the order ratio sequence L as an abscissa, and taking the order ratio spectrum sequence S3(t) drawing a scale spectrum with the first m values as vertical coordinates, wherein the m value is half of the length of the scale sequence L, and if L is an odd number, rounding the half of the length of the scale sequence L downwards;
and observing whether fault characteristic order ratios of all parts of the robot joint exist in order ratio components in the order ratio spectrum, if so, indicating that the parts have faults, and if not, indicating that the joints have no faults, thereby completing vibration signal order ratio analysis and realizing robot joint fault diagnosis and state detection.
2. The joint robot fault diagnosis method based on current and vibration signals according to claim 1, characterized in that: and acquiring a motor current signal I (t) of the ith joint of the robot through a current sensor.
3. The joint robot fault diagnosis method based on current and vibration signals according to claim 1, characterized in that: and acquiring a vibration signal S (t) of the ith joint of the robot through an acceleration sensor.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010918331.9A CN111975784B (en) | 2020-09-03 | 2020-09-03 | Joint robot fault diagnosis method based on current and vibration signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010918331.9A CN111975784B (en) | 2020-09-03 | 2020-09-03 | Joint robot fault diagnosis method based on current and vibration signals |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111975784A true CN111975784A (en) | 2020-11-24 |
CN111975784B CN111975784B (en) | 2022-12-13 |
Family
ID=73447529
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010918331.9A Active CN111975784B (en) | 2020-09-03 | 2020-09-03 | Joint robot fault diagnosis method based on current and vibration signals |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111975784B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113405795A (en) * | 2021-06-22 | 2021-09-17 | 昆明理工大学 | Weak fault identification method for joint RV reducer |
CN113985276A (en) * | 2021-10-18 | 2022-01-28 | 上海电气风电集团股份有限公司 | Fault diagnosis method and device of wind generating set |
CN114441125A (en) * | 2021-12-24 | 2022-05-06 | 北京化工大学 | Non-contact type turbine machinery high-rotating-speed blade vibration multi-mode decomposition method |
CN114454213A (en) * | 2022-01-26 | 2022-05-10 | 东华大学 | Industrial robot joint current signal abnormity detection method based on EEMD-HT-kurtosis analysis |
CN114608821A (en) * | 2022-03-01 | 2022-06-10 | 北京辰宇航康科技有限公司 | Method and device for determining order spectrum of vibration signal of rotating machine |
CN114700988A (en) * | 2022-05-09 | 2022-07-05 | 西安交通大学 | Joint action division method for health monitoring of mechanical transmission part of industrial robot |
CN114970227A (en) * | 2021-11-10 | 2022-08-30 | 昆明理工大学 | Industrial robot joint vibration solving method based on joint current and modal analysis |
CN116533253A (en) * | 2023-07-03 | 2023-08-04 | 佛山智能装备技术研究院 | Industrial robot fault diagnosis method based on feedback current spectrum analysis |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103018043A (en) * | 2012-11-16 | 2013-04-03 | 东南大学 | Fault diagnosis method of variable-speed bearing |
CN103868694A (en) * | 2014-03-26 | 2014-06-18 | 东南大学 | Embedded variable-rotation-speed bearing fault diagnosis device |
US20140324375A1 (en) * | 2013-04-24 | 2014-10-30 | Seiko Epson Corporation | Physical quantity detection circuit, physical quantity detection device, electronic apparatus and moving object |
CN106124197A (en) * | 2016-06-14 | 2016-11-16 | 昆明理工大学 | A kind of epicyclic gearbox sun gear partial fault detection method and system |
CN108120598A (en) * | 2017-12-19 | 2018-06-05 | 胡文扬 | Square phase-couple and the bearing incipient fault detection method for improving bispectrum algorithm |
CN108638128A (en) * | 2018-05-24 | 2018-10-12 | 哈工大机器人(合肥)国际创新研究院 | A kind of real-time method for monitoring abnormality and its system of industrial robot |
CN110617964A (en) * | 2019-07-29 | 2019-12-27 | 中国铁道科学研究院集团有限公司城市轨道交通中心 | Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing |
CN111469128A (en) * | 2020-04-21 | 2020-07-31 | 昆明理工大学 | Current coupling signal separation and extraction method for articulated robot |
-
2020
- 2020-09-03 CN CN202010918331.9A patent/CN111975784B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103018043A (en) * | 2012-11-16 | 2013-04-03 | 东南大学 | Fault diagnosis method of variable-speed bearing |
US20140324375A1 (en) * | 2013-04-24 | 2014-10-30 | Seiko Epson Corporation | Physical quantity detection circuit, physical quantity detection device, electronic apparatus and moving object |
CN103868694A (en) * | 2014-03-26 | 2014-06-18 | 东南大学 | Embedded variable-rotation-speed bearing fault diagnosis device |
CN106124197A (en) * | 2016-06-14 | 2016-11-16 | 昆明理工大学 | A kind of epicyclic gearbox sun gear partial fault detection method and system |
CN108120598A (en) * | 2017-12-19 | 2018-06-05 | 胡文扬 | Square phase-couple and the bearing incipient fault detection method for improving bispectrum algorithm |
CN108638128A (en) * | 2018-05-24 | 2018-10-12 | 哈工大机器人(合肥)国际创新研究院 | A kind of real-time method for monitoring abnormality and its system of industrial robot |
CN110617964A (en) * | 2019-07-29 | 2019-12-27 | 中国铁道科学研究院集团有限公司城市轨道交通中心 | Synchronous compression transformation order ratio analysis method for fault diagnosis of rolling bearing |
CN111469128A (en) * | 2020-04-21 | 2020-07-31 | 昆明理工大学 | Current coupling signal separation and extraction method for articulated robot |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113405795B (en) * | 2021-06-22 | 2023-12-12 | 昆明理工大学 | Method for identifying weak faults of joint RV reducer |
CN113405795A (en) * | 2021-06-22 | 2021-09-17 | 昆明理工大学 | Weak fault identification method for joint RV reducer |
CN113985276A (en) * | 2021-10-18 | 2022-01-28 | 上海电气风电集团股份有限公司 | Fault diagnosis method and device of wind generating set |
CN113985276B (en) * | 2021-10-18 | 2024-02-27 | 上海电气风电集团股份有限公司 | Fault diagnosis method and device for wind generating set |
CN114970227B (en) * | 2021-11-10 | 2024-01-05 | 昆明理工大学 | Joint vibration solving method of industrial robot based on joint current and modal analysis |
CN114970227A (en) * | 2021-11-10 | 2022-08-30 | 昆明理工大学 | Industrial robot joint vibration solving method based on joint current and modal analysis |
CN114441125A (en) * | 2021-12-24 | 2022-05-06 | 北京化工大学 | Non-contact type turbine machinery high-rotating-speed blade vibration multi-mode decomposition method |
CN114441125B (en) * | 2021-12-24 | 2022-10-28 | 北京化工大学 | Non-contact type turbine machinery high-rotating-speed blade vibration multi-mode decomposition method |
CN114454213A (en) * | 2022-01-26 | 2022-05-10 | 东华大学 | Industrial robot joint current signal abnormity detection method based on EEMD-HT-kurtosis analysis |
CN114454213B (en) * | 2022-01-26 | 2023-12-12 | 东华大学 | Industrial robot joint current signal abnormality detection method based on EEMD-HT-kurtosis analysis |
CN114608821A (en) * | 2022-03-01 | 2022-06-10 | 北京辰宇航康科技有限公司 | Method and device for determining order spectrum of vibration signal of rotating machine |
CN114700988B (en) * | 2022-05-09 | 2023-10-24 | 西安交通大学 | Joint motion dividing method for health monitoring of industrial robot transmission part |
CN114700988A (en) * | 2022-05-09 | 2022-07-05 | 西安交通大学 | Joint action division method for health monitoring of mechanical transmission part of industrial robot |
CN116533253A (en) * | 2023-07-03 | 2023-08-04 | 佛山智能装备技术研究院 | Industrial robot fault diagnosis method based on feedback current spectrum analysis |
Also Published As
Publication number | Publication date |
---|---|
CN111975784B (en) | 2022-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111975784B (en) | Joint robot fault diagnosis method based on current and vibration signals | |
Li et al. | Detection of gear cracks in a complex gearbox of wind turbines using supervised bounded component analysis of vibration signals collected from multi-channel sensors | |
CN111469128B (en) | Current coupling signal separation and extraction method for articulated robot | |
CN113405795B (en) | Method for identifying weak faults of joint RV reducer | |
CN104634571B (en) | A kind of Fault Diagnosis of Roller Bearings based on LCD MF | |
CN104819841B (en) | Built-in-coding-information-based single sensing flexible angle-domain averaging method | |
CN108844733B (en) | Gear state monitoring index extraction method based on KL divergence and root mean square value | |
CN110514437B (en) | Planetary gearbox vibration signal separation method for fault diagnosis | |
CN106124197A (en) | A kind of epicyclic gearbox sun gear partial fault detection method and system | |
CN101972170B (en) | Self-adapting filter for least square support vector machine and filtering method thereof | |
CN114509159A (en) | Order ratio tracking analysis method, system and computer readable storage medium | |
CN111510026A (en) | Output torque estimation method and system for permanent magnet synchronous motor | |
CN115683617A (en) | Improved variable working condition fault diagnosis method and system | |
CN114298110A (en) | Rolling bearing fault diagnosis method and system based on interpretable 1DCNN model | |
CN107224284A (en) | The noise detecting method and system of a kind of digital electrocardiosignal | |
CN111582248A (en) | SVD-based gearbox signal noise reduction method | |
CN116734980A (en) | Non-contact type mechanical vibration monitoring and fault diagnosis method based on event camera | |
CN114112396B (en) | Bearing fault diagnosis method under rotating speed fluctuation quasi-stationary working condition by adopting current signal analysis | |
Wolf et al. | Using the traction drive as the sensor to evaluate and track deterioration in electrified vehicle gearboxes | |
US20230008153A1 (en) | Vibration compensation controller with neural network band-pass filters for bearingless permanent magnet synchronous motor | |
CN112526339B (en) | Polynomial-phase transformation-based squirrel-cage motor rotor broken bar fault diagnosis method | |
CN114705422A (en) | Harmonic reducer torsional mode testing method and system | |
CN115683580A (en) | Rotating machinery fault diagnosis method based on phase domain adaptive filtering | |
CN102798462B (en) | No-time-mark order tracking method based on self-demodulation transform | |
Li et al. | Gear fault diagnosis based on order tracking and Hilbert-Huang transform |
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 |