Next Article in Journal
On the Modelling of Tethered Mobile Robots as Redundant Manipulators
Previous Article in Journal
The Operation of UAV Propulsion Motors in the Presence of High External Magnetic Fields
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Method for Health Indicator Evaluation for Condition Monitoring of Industrial Robot Gears

by
Corbinian Nentwich
* and
Gunther Reinhart
Institute for Machine Tools and Industrial Management, Boltzmannstraße 15, 85747 Garching, Germany
*
Author to whom correspondence should be addressed.
Submission received: 29 April 2021 / Revised: 27 May 2021 / Accepted: 7 June 2021 / Published: 9 June 2021

Abstract

:
Condition monitoring of industrial robots has the potential to decrease downtimes in highly automated production systems. In this context, we propose a new method to evaluate health indicators for this application and suggest a new health indicator (HI) based on vibration data measurements, Short-time Fourier transform and Z-scores. By executing the method, we find that the proposed health indicator can detect varying faults better, has lower temperature sensitivity and works better in instationary velocity regimes compared to several state-of-the-art HIs. A discussion of the validity of the results concludes our contribution.

1. Introduction

Industrial robots are a fundamental part of highly automated production systems, which can be found in the automotive or electronics industry [1]. Since they operate in complex production cells and as a part of linear production lines, robot malfunctions lead to long downtimes for repair or replacement and, hence, to increased costs. In particular, robot gear faults are responsible for the longest downtimes because they often require the replacement of the whole robot [2]. The condition monitoring (CM) of these gears offers the potential to resolve this issue. CM is the monitoring of an asset’s health using sensor data. The health state represents a wear reserve before a failure occurs. This health state is quantified with a health indicator (HI). A significant monitored change in this health indicator can be used as a decision-making aid in the planning of maintenance actions [3].

1.1. State-of-the-Art

In recent years, different HIs based on vibration data for several industrial robot components, such as bearings, gears and motors, and their specific faults have been investigated. Furthermore, several approaches to cope with instationary signals in CM have been presented. The next two sections give a short overview of these topics followed by a section stating the contribution of our publication.

1.1.1. Vibration-Based Robot Condition Monitoring

A fault detection method was developed in [4], which first uses a novel phase-based, time-domain averaging method to remove the deterministic part of the vibration signal. Subsequently, the root mean square (RMS) and power spectrum entropy of the remaining residual signal are calculated as health indicators. A vibration signal based CM system for SCARA robots was implemented in [5], which in the first step uses statistical HIs of the time-domain signal to detect the occurrence of a defect and in the second step uses an artificial neural network to diagnose the fault type. A three-layer architecture for remote fault diagnosis of industrial robot gearboxes was proposed using vibration signals in [6]. In the diagnosis layer, the authors present a performance evaluation approach using a support vector machine (SVM), a remaining useful life prediction by a Markov model and a fault-type diagnosis based on a Bayesian network. The degenerative behavior of an industrial robot gear was observed with vibration sensors by [7] as well as [8] in accelerated wear tests. After pre-processing the signals using order tracking and spectral auto-correlation, the characteristic fault frequencies were calculated and monitored by root mean square analysis, which revealed a trend correlating with increasing wear. In addition to the installation of accelerometers, other additional data sources were investigated in this context. The acoustic emission technology was used to detect robot gearbox faults based on the ball spinning and ball passing frequency of the bearings in [9]. The changes of the RMS-HI and characteristic frequencies for functional and broken strain gears of industrial robots were investigated in [10]. The classification and regression performance of different data-driven models based on frequency-domain data and principal component analysis for dimensionality reduction was evaluated in [11].

1.1.2. Time–Frequency-Based Health Indicators

In addition to vibration data based CM approaches for industrial robots, there also exist several publications considering HIs from the time–frequency-domain. Here, approaches based on the Short-time Fourier transform (STFT), Wavelet transform (WT) or Hilbert Huang transform (HHT) can be divided. STFT is used to derive two HIs named Prominence and Compliance in [12] to detect bearing faults based on their characteristic fault frequencies. A similarity measure between the STFT spectrograms based on standard deviation and correlation is combined with a simple classifier in [13] to detect bearing faults. The same objective was pursued in [14] by means of the marginal time integration of STFTs. Bearing fault classification by means of non-negative matrix factorization or convolutional neural networks and STFT was evaluated in [15,16].
In the field of WT, several approaches exist for different assets. CM of brushless DC motors is investigated based on energies for characteristic frequencies based on both STFT and WT in [17]. A decomposition rate is used in [18] for CM of electric drives based on WT. RMS and Kurtosis are calculated for the WT coefficients for broken bar fault detection in electric drives and combined with a neural network for fault classification in [19]. Bearing fault classification was performed with an SVM based on WT in combination with singular value decomposition for dimensionality reduction in [20]. The spectra of WT coefficients were the basis for the calculation of statistical HIs and frequency specific energy values for the CM of bearing faults in [21]. The similarities of continuous WT spectra are used as an HI for bearing fault detection [22]. The permutation entropy derived from flexible analytical wavelet transform was used as a feature for an SVM for bearing fault classification [23]. Impulse factor, Kurtosis and RMS based on WT coefficients were used for bearing fault detection of helicopters [24]. Statistical features and Hoelder’s exponent were derived from WT coefficients for milling tool health state monitoring. Here, the HIs were the input for an SVM and Decision Tree classifier [25]. HIs were derived by a convolutional neural network for milling tool condition monitoring based on the wavelet decomposition in [26]. Energies of WT coefficients were also used for detecting generator and gear faults in wind turbines [27]. Different entropy-based and statistical features were used in [28] for gearbox health monitoring in combination with an SVM. Energy and entropy values derived from WT for characteristic frequencies are applied for gearbox condition monitoring in [29].
In [30], the Shannon entropy based on HHT was used for the CM of gears. HHT was also used in [31] to derive HIs by an autoencoder based on the Marginal Hilbert Spectrum. A component dependent frequency energy based on HHT was used as a label in [32] for a CNN-based regression model trained on raw vibration time series data for bearing fault detection. Different statistical and entropy-based HIs were calculated from the Intrinsic mode functions (IMFs) derived by HHT in [33].

1.2. Contribution to the State-of-the-Art

However, none of these publications assess vibration data-based HIs’ ability to detect faults in an industry-like industrial robot application setting. It is characterized by changing robot axes’ velocities, changing temperatures of the gears due to unbalanced robot utilization and unknown robot gear fault types. This is why we present a new HI for robot gear condition monitoring, which potentially copes with these characteristics. Furthermore, we propose a method to evaluate the suitability of HIs for the task of robot gear condition monitoring. We apply this method on the newly formulated HI and several HIs from the state-of-the-art.

2. Materials and Methods

This section is divided in two parts. First, the newly developed HI is presented. Afterwards, the methodologies to evaluate the HI’s performance and data sets used in this context are explained.

2.1. Time–Frequency-Domain-Based Z-Score

The concept of the newly designed HI is based on two cornerstones. To deal with instationary velocity regimes, which are found in robot applications due to the typical movement patterns of a robot, the HI is based on time–frequency-domain data. Simultaneously, the HI must take into account a certain variance of this data due to environmental changes such as temperature fluctuations. This is realized by the concept of Z-scores, a common similarity measure from statistics [34]. The process to calculate the new HI is depicted in Figure 1.
In detail, the new HI is based on high-frequency sampled acceleration sensor data. Data from one measurement are transformed to a time–frequency-spectrogram by usage of the STFT, which is calculated according to Equation (1). Here, τ and ω are time and frequency indices, x ( n ) is the time series signal of the vibration signal at timestep n and w is a windowing function with the length R.
s p e c ( τ , ω ) = | n = x ( n ) w ( n τ R ) e j ω n |
To set up the HI, a certain number of vibration signal spectrograms must be collected for the robot to capture its signal signature in a healthy state with its stochastic variations. This takes place in an initialization phase. For this, initially, two measurements must be collected. In this context, a measurement is defined as the collection of vibration data over one single movement. Based on this data, the two spectrograms are calculated. To determine whether this reference quantity of two spectrograms captures the stochastic variation of the signal, the overall mean (Equation (2)) and standard deviation (Equation (3)) of the spectrograms are calculated.
s p e c ( τ , ω ) a v g = 1 k i = 0 k s p e c ( τ , ω ) i
s t d s p e c , o v e r a l l = 1 0.5 F T τ = 0 T ω = 0 0.5 F i = 0 k ( s p e c ( τ , ω ) i s p e c ( τ , ω ) a v g ) 2 k
In these formulas, k describes the number of measurements in the reference quantity. T is the time length of each measurement, F is the sampling frequency and s p e c ( τ , ω ) a v g is the average value of s p e c ( τ , ω ) over measurements 0 to k. Afterwards, one measurement is added to the reference quantity at a time, and again a v g s p e c , o v e r a l l and s t d s p e c , o v e r a l l are calculated. Plotting these standard deviations over the number of measurements in the reference quantity usually first shows an increase in s t d s p e c , o v e r a l l and then a saturation as can be seen in Figure 2. If this saturation is reached, the reference quantity can sufficiently represent the stochastic behavior of the signal signature. In the shown example, this saturation is reached after 5 measurements.
After the initialization, an HI can be determined based on a newly collected measurement. For this, the measurement’s spectrogram overall Z-score is determined according to Equation (4).
H I m e a s = 1 0.5 F T τ = 0 T ω = 0 0.5 F | s p e c ( τ , ω ) m e a s s p e c ( τ , ω ) a v g , r e f s p e c ( τ , ω ) s t d , r e f |
In this context, s p e c ( τ , ω ) a v g , r e f and s p e c ( τ , ω ) s t d , r e f are the mean value and the standard deviation of s p e c ( τ , ω ) for all measurements in the reference quantity. In Figure 3, the STFT and Z-score spectrograms of exemplary vibration measurements from a healthy and a faulty robot gear are depicted. The Z-score-based spectrogram of the faulty measurement shows more prominent changes compared to the STFT-based spectrogram.

2.2. Hi Evaluation Method

To compare the ability of the newly designed HI to cope with industrial robot application characteristics, we followed a three step approach. First of all, we investigated how well the designed HI can detect different kinds of faults in comparison to HIs from the state-of-the-art. Second, we investigated the temperature sensitivity of HIs from the state-of-the-art meeting this criterion and our HI. Third, we investigated the trend behavior of HIs showing a low temperature sensitivity on data from two accelerated wear tests. These three steps are now described more precisely. The overall process of our investigations is also described in Figure 4.

2.2.1. Varying Fault Detection Analysis

We used the FEMTO data set, which is described in detail in [35], to select HIs capable of detecting different faults. The data set is available in [36]. This data set provides run to failure vibration data from 16 identical bearings and for different faults and working conditions defined by the applied load and the rotational speed. The acceleration sensor sampled data with 25.6 kHz, one measurement has a length of 0.1 s and measurements were taken in equidistant timesteps of 10 s for all bearings. The test run for one bearing ended when the signal from the acceleration sensor exceeded 20 g. Therefore, different numbers of measurements are available per bearing ranging from 230 to 2803. We calculated the HIs summarized in Table 1 for all measurements of one sensor. These HIs were derived from several review papers regarding gearbox and bearing CM [37,38,39,40] and the publications mentioned in Section 1. Therefore, the HI calculation was based either on the raw acceleration signal, an enveloped signal as described in [41] or the residual signal as suggested by [4]. Additionally, the newly designed HI presented in Section 2 was calculated for the measurements based on the raw signals.
To detect whether these HIs are sensitive to multiple faults, different techniques can be applied. In addition to filter techniques, ensemble, wrapper and embedded methods exist [42]. However, the latter three techniques combine classification or regression models with HIs for their evaluation. Hence, this evaluation is always dependent on the used models. Thus, we chose to use filter methods for the evaluation. Here, different figures of merit for regression and classification tasks can be applied, such as trendability, robustness, monotony or discriminance [42]. To combine these different performance indicators, we fitted different basic functions on the HIs calculated for the last 20 percent of measurements per bearing. These functions were first and second degree polynoms, exponential and sigmoid functions. For each of the fits, we calculated the R2 value. This means that we received four R2 values per HI and bearing. High R2 values of these fits correlate with a high trendability, monotony, robustness and discriminance, which are desirable for HIs. To evaluate whether an HI can detect several damages, we considered only the best R2 value per HI and bearing. We plotted the statistics of these 16 remaining R2 values per HI as a boxplot. Suitable HIs should show high R2 values with low variance.

2.2.2. Temperature Sensitivity Analysis

HIs showing this behavior were analyzed regarding their temperature sensitivity. For this purpose, we acquired vibration data from an industrial robot test rig. This test rig consists of a KUKA KR510 industrial robot with an attached load of 365 kg. We attached acceleration sensors close to the gearboxes as shown in Figure 5 on the right side. These sensors have a sampling rate of 26 kHz. The acceleration direction of the sensors was orthogonal to their contact area. For data acquisition, the robot performed a trajectory where each joint was moved individually at different speeds in an angle area of 10°, as described in Figure 6, and for different gear temperatures in the range of 25 °C and 60 °C and 5 °C steps. One measurement per axis lasted 8 s. The gear temperature was measured at the gearbox cap with an infrared thermometer. For each temperature step, four measurements were made. For each measurement at each temperature step, the remaining HIs were calculated. To determine the temperature sensitivity, we divided the average HI values calculated from measurements at the highest gear temperatures by the values calculated from measurements at the lowest temperature. HIs with a high sensitivity were eliminated for the last step.

2.2.3. Accelerated Wear Test Analysis

Here, we calculated the remaining HIs for measurements from two data sets from accelerated robot wear tests to see how these HIs perform in a more industry like setting and how they cope with instationary velocity behavior. The first data set was collected during a time range of approximately one year with an ABB robot of type RB 6600-255/2.55. During the data acquisition, the robot performed an isolated movement of the second axis in an angle area of 150° for each measurement. Vibration data were only acquired with a sensor attached axially at the robot axis 2 gearbox. At the end of the experiment, the gearbox was dismantled and faults on the bearings and the shafts of the gear were found. A total of 2290 measurements, equally distributed over time, were taken for our analysis from this data set. One measurement lasted 1.6 s and the sampling rate was 10 kHz. More detailed information about this experiment can be found in [7,8]. The second data set was derived from another experiment. Here, the second axis of an ABB IRB 7600-340/2.8 was moved in an angle area of 80° continuously over the time frame of three months. The vibration sensor attached to the gearbox cap of axis 2 sampled with 20 kHz and one measurement lasted 2.15 s. The measurement setup is presented on the left side in Figure 5. The experiment ended after a roller element of a bearing had cracked and had blocked the gear. In this time range, 920 vibration measurements were taken in total in equidistant time steps. The faults, which occurred in both experiments, can be seen in Figure 7. In both experiments, environmental conditions such as load and trajectory were kept constant. Fluctuations of the temperature were kept at a minimum due to the constant movements of the robots. In this way, signal changes are likely to be correlated to increasing wear.

3. Results

This section is divided in three parts. First, the results from the varying fault detection experiments are shown. Secondly, the results from the temperature sensitivity analysis are presented. Finally, the application of the HIs on the two accelerated wear tests is described.

3.1. Varying Fault Detection Analysis

From the 16 bearing experiments, the HIs presented in Table 1 were calculated. We used the first 100 measurements per bearing as the reference quantity for the Z-score-HI and set R to 128. Figure 8 shows the R2 values for a selection of different HIs as a box plot. The R2 statistics for all HIs can be found in Appendix A. The abbreviations of the HIs are explained in Table 1. The PtP-, Peak-, RMS-, Std- and Z-score-HI show the highest R2 values on average. They also show the lowest variance between the different bearings. This means that these HIs detect different faults most reliably. Other HIs show also high trend values but only for some of the bearings. HIs derived from the frequency-domain (DomF, SpC, SpE, SpF, SpRO) perform worse compared to HIs from the time-domain. The preprocessing steps of enveloping the signal or calculating the residual signal do not affect the HI trend behavior significantly, which can be seen in Table A1, Table A2 and Table A3. The TDI-, and DWTRMS-HI for specific frequency bands also show high average values with changing variance (see Table A4). If these HIs would be used for robot gear condition monitoring, the progress of all frequency band specific HIs would have to be tracked as different faults stimulate changes in different frequency bands.

3.2. Temperature Sensitivity Analysis

Based on this result, we conducted the temperature sensitivity analysis for the PtP-, Peak-, RMS-, Std-, TDI-, DWTRMS- and Z-score-HI. Here, we used one measurement per temperature step as the reference quantity for the Z-score-HI and set R to 128. Figure 9 shows the change of the HIs per axis in percent for the PtP-, Peak-, RMS-, Std- and Z-score-HI. The RMS- and Z-score-HI show the lowest temperature sensitivity overall. Figure 10 shows the results for the DWTRMS-HIs. Here, high sensitivities for different detail coefficient DWTRMS-HIs exist. Figure A1 shows the temperature sensitivity of the TDI-HIs of different frequency bands. Here, a similar result can be seen compared to the DWTRMS-HIs. The data of Figure 10 and Figure A1 can also be found in Table A5 and Table A6. In general, the data from axis 4 show the highest temperature sensitivity for all HIs. The comparably higher sensitivity of the HI values derived from data at axis 4 can be related to the robot trajectory. During the trajectory, the robot arm was stretched out, which leads to greater elasticity at the position of the sensor at axis 4. This can cause increased vibrations, which are magnified under changing temperature influences. Given the results of the temperature sensitivity analysis, we analyzed the data sets from the accelerated wear tests with only the RMS- and the Z-score-HI. The other HIs were excluded due to their high temperature sensitivity. Even though some frequency band specific DWTRMS-HIs and TDI-HIs show low sensitivity, they were excluded as robot gear faults do not have to stimulate these frequency bands with low sensitivity.

3.3. Accelerated Wear Tests Analysis

In this analysis, we used the first 100 measurements as the reference quantity for the Z-score-HI and set R to 256. For smoothing, we applied a rolling average with a window length of 15 on both HI series. The progress of the HIs in the accelerated wear test of the ABB IRB 7600 is shown in Figure 11. Both HIs show a plateau with increased values at the end of the experiment. It can be assumed that, at this point in time, faults have already been present. Here, the increased HI values over a longer time period could have been used as a decision criterion for maintenance actions.
The measurements at the very end show decreased values again. We assume that this decrease is correlated to a part of the bearing roller. In the end of the experiment, one of the roller elements showed a large pit. During the measurements showing the higher HI values this detached part of the roller element could have been still slightly fixed at the roller element and thus could have caused high vibration. After full detachment, this noise level decreased again. For the measurements before the plateau, the RMS-HI shows higher fluctuations compared to the Z-score-HI. For instance, the RMS-HI shows a first high peak around measurement 100. Such peaks could lead to false alarms in a condition monitoring scenario and should be avoided.
The progress of the HIs in the other accelerated wear test performed with the ABB IRB 6600 is shown in Figure 12. Here, the Z-score-HI shows a trending behavior and the RMS shows a stationary progress. Both HIs show a high increase during the last measurements. In this experiment, the trending behavior of the Z-score could have been a criterion to execute maintenance actions. This information is not present in the RMS-progress. Based on the fact that the Z-score showed a better trend behavior in the ABB IRB 6600 experiment and less noisy behavior in the ABB IRB 7600 experiment, we suggest the use of the Z-score-HI for the condition monitoring of robot gears.

4. Discussion

The discussion is divided in four parts. First, some remarks regarding our designed HI are given. Afterwards, three parts make up the Results subsections.
To derive the spectrograms required for the Z-score-HI, the length of the window function must be defined. High values for R result in a high frequency resolution and low values in a high time resolution. For the individual experiments, we chose window lengths that lead to a good compromise between time and frequency resolution by inspecting spectrograms created with different window lengths. We chose window lengths that lead to spectrograms appearing the least noisy in a visual inspection. In an industrial setting, an automated approach should be developed for this dependent on the robot’s trajectory and the used sensor.
The motivation to use the FEMTO data set to investigate HI performance was to assess HIs’ capability to detect multiple faults. Within a robot gearbox, which are mostly RV reducers, not only bearings but also the gear teeth can have faults. Such faults are not taken into account by our analysis explicitly. However, the bearing faults present in the FEMTO data set, e.g., pitting, are similar to typical gear teeth or shaft damage from a signal analysis point of view. Damage from all components modulate the acceleration signals at a specific frequency and its sidebands. Exactly this capability to track such changes in the signal was investigated in our analysis. There also exist HIs that track energy changes at the specific component fault frequencies. Such HIs were excluded from our analysis because expert knowledge about the geometric characteristics of the gears, e.g., the bearing diameters or the number of roller elements, is required to calculate these HIs. This expert knowledge is usually not available to industrial robot users. We also excluded HIs that could be derived automatically from machine learning models, such as autoencoders, as the physical interpretation of these HIs is difficult and hence a transferability between different robot systems is questionable from our point of view.
Regarding the results of the temperature sensitivity analysis, it must be pointed out that the results are valid only for the chosen robot trajectory. As the dynamic behavior of the robot changes within its working space, this analysis should be performed individually for trajectories and robot systems. However, from a theoretical point of view, the Z-score-HI possesses the ability to cope with these temperature fluctuations independently of the trajectory. Temperature variations lead to variance in the STFT spectrograms. This variance is taken into account in the s p e c ( τ , ω ) a v g , r e f and s p e c ( τ , ω ) s t d , r e f during the initialization phase. Hence, Z-score-HIs derived from measurements from functional robot gears and different temperatures will show only little differences in the Z-score-HI value. This becomes more clear considering Figure 13. Here, the STFT and Z-score spectrograms from two vibration measurements of the temperature sensitivity experiment are shown. On the left side, the spectrograms from a cold gear measurement are depicted. On the right side, the spectrograms from a warm gear measurement are shown. Differences are visible in the STFT spectrograms around seconds 1 and 2. No differences are visible in the Z-score spectrograms. The scales of the STFT spectrograms reach from −5 to 0 and the scales of the Z-score spectrograms from 0 to 1.5. Hence, the relative changes of the STFT spectrograms are bigger compared to the Z-score spectrograms. In this example, the total relative change in energy in the STFT spectrogram is 9.15 percent, whereas the total relative change in the Z-score spectrogram is just 1.63 percent.
Finally, the results from the accelerated wear tests show noisy progress over time. This hinders a simple or automated detection of faults in a condition monitoring behavior. To establish an automated CM system, a suitable trend detection in combination with an outlier detection system must be set up. A trend detection system could identify HI progress shown as in Figure 12, whereas an outlier detection system could detect progress as depicted in Figure 11. The development of such a system also marks the outlook of our future work.

5. Conclusions

Condition monitoring of robot gears has the potential to decrease production system downtimes. The state-of-the-art provides many health indicators to track the health state of gears. We analyzed these health indicators regarding specific requirements rising from typical industrial robot applications. These requirements are the ability to detect different faults, low temperature sensitivity and the capability to deal with instationary velocity behavior. Additionally, we suggested a new health indicator based on STFT spectrograms and Z-scores that can cope with these requirements. Our analysis showed that the RMS health indicator and our suggested health indicator meet the defined requirements the best. Data from accelerated wear tests show that for an automatic condition monitoring system a combination of a trend detection and an outlier detection system that can deal with a noisy signal is required.

Author Contributions

Conceptualization, C.N. and G.R.; methodology, C.N.; software, C.N.; validation, C.N.; formal analysis, C.N.; investigation, C.N.; resources, G.R.; data curation, C.N.; writing—original draft preparation, C.N.; writing—review and editing, G.R.; visualization, C.N.; supervision, G.R.; project administration, C.N.; funding acquisition, G.R. All authors have read and agreed to the published version of the manuscript.

Funding

We express our gratitude to the Bavarian Ministry of Economic Affairs, Regional Development, and Energy for the funding of our research. The formulated outlook will be investigated as part of the research project “KIVI” (grant number IUK-1809-0008 IUK597/003) and will be further developed and implemented.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to confidentiality reasons.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table A1. R2 statistics for HIs derived from the normal signal.
Table A1. R2 statistics for HIs derived from the normal signal.
CrFDomFImpFKurtMarFMeanMedMedFPeakPtPRMSSkewSpCSpESpFSpROStdZ-Score
Mean0.2310.3540.2870.3360.2960.0340.0630.4670.8220.8440.8870.2420.5990.4880.5140.4910.8870.934
Std0.2460.3510.2890.3340.2980.0760.1470.3110.1730.1450.2150.2890.2980.2810.2250.3220.2150.080
Min0.0050.0130.0070.0040.0060.0030.0070.0240.2590.3970.0820.0040.0250.0670.1350.0350.0820.644
Max0.7060.9980.7640.8690.7730.3230.6180.9410.9870.9850.9900.8880.9440.9710.9000.9840.9900.983
Table A2. R2 statistics for HIs derived from the enveloped signal.
Table A2. R2 statistics for HIs derived from the enveloped signal.
CrFDomFImpFKurtMarFMeanMedMedFPeakPtPRMSSkewSpCSpESpFSpROStd
Mean0.2150.1390.2840.2990.2960.7760.8190.4960.8160.8160.8720.2750.6050.4640.4620.5140.898
Std0.2290.2480.2830.3280.2910.2960.2480.3180.1780.1780.2270.3030.2860.2980.2760.3090.131
Min0.0050.0010.0060.0010.0050.0100.0140.0230.2460.2460.0740.0080.0460.0130.0120.0070.427
Max0.6350.9970.7410.9190.7660.9770.9830.9810.9880.9880.9890.9050.9390.9780.9010.9870.987
Table A3. R2 statistics for HIs derived from the residual signal as suggestet by [4].
Table A3. R2 statistics for HIs derived from the residual signal as suggestet by [4].
CrFDomFImpFKurtMarFMeanMedMedFPeakPtPRMSSkewSpCSpFSpROStdSpE
Mean0.3290.4230.3480.3650.3550.0320.0880.6080.8470.8590.8840.2370.7010.6090.6050.8840.534
Std0.2700.3480.3030.3370.3120.0890.1870.2920.1780.1530.2210.3190.1960.2020.2130.2210.316
Min0.0100.0210.0030.0070.0180.0010.0010.0160.2110.3280.0820.0030.4250.1470.1320.0830.023
Max0.8020.9920.8280.8580.8300.3740.6710.9880.9870.9870.9900.9430.9590.9300.9770.9900.969
Table A4. R2 statistics for HIs derived from the time–frequency-domain. DWT-, IMF-, and TDI-based HIs were calculated for different frequency bands. The frequency bands are encoded in the abbreviation of the HI name. Large numbers correspond to high frequency bands for TDI-HIs and low frequency bands for DWT- and IMF-HIs.
Table A4. R2 statistics for HIs derived from the time–frequency-domain. DWT-, IMF-, and TDI-based HIs were calculated for different frequency bands. The frequency bands are encoded in the abbreviation of the HI name. Large numbers correspond to high frequency bands for TDI-HIs and low frequency bands for DWT- and IMF-HIs.
MeanStdMinMax
DWTRMS50.8910.2070.0990.990
DWTRMS40.8890.2230.0320.988
TDI330.8820.0870.5830.965
DWTRMS60.8820.2160.1170.991
TDI430.8710.1560.2960.983
IMFRMS20.8690.2270.0140.991
TDI360.8670.2030.0950.976
TDI140.8650.1540.3270.970
TDI350.8640.2020.0960.981
TDI220.8620.1960.1350.963
TDI340.8610.1970.1140.974
TDI230.8600.1810.2000.961
TDI210.8590.1970.1470.971
DWTRMS30.8590.2450.0190.998
TDI370.8570.2080.0720.969
TDI390.8530.2210.0170.975
TDI150.8470.2180.0350.982
IMFRMS10.8450.2100.0910.978
TDI440.8440.2230.0040.984
TDI130.8420.2240.0150.975
TDI120.8410.2180.0220.969
TDI400.8400.2150.0720.969
TDI240.8370.2230.0270.954
TDI410.8360.2210.0610.970
TDI160.8330.2280.0220.969
TDI420.8300.2210.0720.964
TDI180.8300.2370.0440.976
TDI100.8300.2140.0250.966
TDI50.8270.2730.0320.995
TDI60.8250.2630.0150.983
TDI80.8250.2060.0510.948
TDI90.8240.2090.0470.954
TDI40.8230.2730.1510.995
TDI450.8200.2200.0130.986
TDI200.8180.2770.0660.963
TDI380.8160.2590.0380.972
TDI320.8150.2720.0900.968
TDI460.8130.2180.0230.987
TDI70.8130.2640.0080.973
TDI470.8120.2170.0130.988
TDI110.8120.2510.0170.966
TDI170.8080.2930.0050.967
DWTRMS20.8040.2970.1010.981
IMFRMS30.8030.2430.0250.961
TDI30.8030.2900.1480.987
TDI490.7920.2150.0180.990
TDI190.7820.3040.0260.959
TDI280.7820.2700.0140.960
TDI500.7800.2160.0260.989
TDI520.7700.2160.0250.989
TDI300.7690.2930.0490.964
TDI250.7690.2900.0250.966
TDI530.7670.2150.0250.989
TDI540.7660.2150.0250.990
TDI480.7650.2520.0240.988
TDI550.7650.2150.0280.989
TDI580.7620.2150.0260.990
TDI570.7610.2160.0280.989
TDI600.7560.2220.0280.989
TDI310.7470.3120.0730.971
TDI510.7470.2530.0220.989
TDI290.7440.2970.0320.961
TDI20.7440.3190.0690.977
DWTRMS10.7400.3060.0260.978
TDI270.7390.2920.0270.954
TDI630.7360.2530.0280.989
TDI620.7360.2540.0260.990
TDI590.7350.2540.0240.990
TDI560.7350.2540.0230.989
TDI610.7340.2540.0240.989
TDI640.7320.2540.0240.989
TDI260.7120.3420.0430.964
IMFRMS40.6650.3350.0350.985
HHTentr0.6460.3170.0210.982
TDI10.5880.3130.0110.937
IMFRMS50.5670.3130.0160.934
IMFRMS60.4600.2230.1100.829
DecompRate0.4580.2920.0750.974
DWTKurt50.3830.2960.0140.821
DWTKurt20.3810.3180.0070.815
DWTEntr60.3800.3190.0170.991
DWTImpF20.3770.3080.0050.806
DWTKurt40.3610.3150.0070.851
DWTKurt60.3530.3300.0030.860
DWTImpF60.3270.2800.0010.796
DWTImpF50.3260.2740.0010.766
IMFRMS70.3260.2710.0260.863
TDI00.3250.2460.0370.906
DWTImpF30.3250.2600.0170.701
IMFKurt10.3200.3140.0090.816
DWTKurt30.3190.2960.0070.827
DWTImpF40.3020.2730.0040.748
IMFImpF10.2880.2800.0090.838
IMFKurt30.2870.2880.0170.843
IMFKurt40.2860.2220.0240.677
DWTKurt10.2820.2170.0120.722
IMFImpF40.2770.2020.0140.601
IMFKurt20.2630.2740.0100.689
IMFEntr20.2540.2650.0120.867
DWTEntr50.2540.2910.0010.989
DWTImpF10.2510.2130.0060.660
IMFImpF30.2420.2530.0160.761
IMFEntr10.2320.2800.0090.944
IMFImpF20.2240.1850.0110.602
IMFEntr30.2070.2690.0130.888
IMFKurt50.1840.2150.0070.822
IMFRMS80.1810.2380.0000.682
IMFImpF50.1640.1870.0070.693
IMFEntr40.1530.1380.0050.524
DWTEntr40.1360.2350.0050.969
IMFImpF60.1210.1180.0100.377
IMFKurt60.1210.1250.0100.428
IMFEntr60.1050.1020.0090.393
DWTEntr30.1010.1870.0030.783
IMFEntr50.0950.0950.0020.367
IMFEntr70.0890.0920.0030.276
DWTEntr20.0620.1110.0010.417
IMFImpF70.0480.0390.0020.134
IMFKurt70.0420.0330.0020.106
IMFEntr80.0320.1110.0000.461
DWTEntr10.0310.0340.0030.141
IMFImpF80.0260.0280.0000.078
IMFKurt80.0210.0310.0000.114
IMFRMS90.0000.0000.0000.000
IMFKurt90.0000.0000.0000.000
IMFImpF90.0000.0000.0000.000
IMFEntr90.0000.0000.0000.000
IMFRMS100.0000.0000.0000.000
IMFImpF100.0000.0000.0000.000
IMFKurt100.0000.0000.0000.000
IMFEntr100.0000.0000.0000.000
IMFRMS110.000NaN0.0000.000
IMFImpF110.000NaN0.0000.000
IMFKurt110.000NaN0.0000.000
IMFEntr110.000NaN0.0000.000
Table A5. Temperature sensitivity of the different DWTRMS-HIs.
Table A5. Temperature sensitivity of the different DWTRMS-HIs.
1234
DWTRMS14.94473441.83325013.22581727.640727
DWTRMS29.43177961.78438620.96640552.444752
DWTRMS311.10087067.90634117.37229064.555873
DWTRMS416.17632288.06401519.29471582.614996
DWTRMS513.63168687.43976317.45168190.271799
DWTRMS610.01579770.88787019.14173393.391435
Table A6. Temperature sensitivity of the different TDI-HIs.
Table A6. Temperature sensitivity of the different TDI-HIs.
Axis 1Axis 2Axis 3Axis 4
TDI07.8876803.6236264.0324289.856386
TDI17.34523013.7801921.41932915.104011
TDI25.70221746.38768723.87637224.639909
TDI36.72135547.59046022.69175642.390072
TDI410.87382959.62150324.15864138.671919
TDI513.12775274.78321917.39807554.988992
TDI614.67051659.57625512.95342475.739203
TDI718.16177957.59683117.54513265.254182
TDI815.09043370.50183011.00450868.946535
TDI919.13081994.75656710.26787363.083458
TDI1019.817337140.43464822.21261148.494003
TDI1131.037879152.47370528.13432175.366125
TDI1233.67587393.05755918.95318498.385403
TDI1326.25089279.86600217.15270692.669107
TDI1411.63589658.28693324.51006986.809206
TDI1523.37737946.90790316.49655372.502066
TDI1630.06421355.24385511.11979959.119520
TDI1721.26874364.3325887.93880379.776374
TDI189.80863952.9331434.19785687.106639
TDI193.42040857.1023013.50714094.175140
TDI200.06262958.6977285.538016107.584154
TDI2117.43381560.09360512.872000112.916761
TDI2223.98043865.0381029.486835164.461250
TDI2323.74931891.05538911.379480152.353906
TDI2411.483126109.84274621.619409121.442989
TDI256.59670685.65375323.250270111.067700
TDI2616.05735068.79413220.72456688.922105
TDI2716.55767068.37353918.973086114.047417
TDI2819.37207565.93743013.256130140.639806
TDI2922.61342768.08111827.002980129.175849
TDI3015.44124168.25534931.219217111.829705
TDI314.72954668.35345828.824115116.071704
TDI321.11263069.44068318.666493104.393062
TDI333.01001866.10735413.594603104.544706
TDI346.22313649.77121214.597765135.411471
TDI3510.11904233.5365904.257545167.094574
TDI360.77707933.9634670.393118158.518025
TDI370.89358447.21186718.367647110.996278
TDI388.78940639.55085333.93813655.254421
TDI3918.2865514.90381723.61578135.301715
TDI4028.76752534.2240596.48552447.061077
TDI4121.67049539.7284169.69609333.227046
TDI424.43781759.76590818.54928879.919517
TDI430.14752352.93938919.84825283.294446
TDI440.44605753.53944619.35132679.595920
TDI450.65708252.35651319.03209277.824663
TDI460.19994755.87380019.52865675.827413
TDI471.38877557.01016819.01422576.394986
TDI482.09711757.77428618.88982276.077127
TDI493.16492658.95906518.69679675.814731
TDI503.85151459.59262118.41040475.719379
TDI514.57758760.33944618.30818075.543086
TDI525.16025860.87709518.21134075.458546
TDI535.65931561.36514418.06912675.348426
TDI546.07712861.76121817.97663275.218693
TDI556.43642662.07761217.94554975.121278
TDI566.76325762.36123017.85147375.050770
TDI577.08073062.59144717.78136774.975471
TDI587.36966762.78164617.78367374.880725
TDI597.62292462.94683017.71025974.853578
TDI607.82853263.07838417.67063774.815067
TDI617.97454463.17774617.69392574.733074
TDI628.07217363.25969817.63693374.758740
TDI638.12885863.30895917.62371574.750804
TDI648.14715863.32221017.66445674.689287
Figure A1. Temperature sensitivity for the different TDI-HIs and robot axes.
Figure A1. Temperature sensitivity for the different TDI-HIs and robot axes.
Robotics 10 00080 g0a1

References

  1. Krockenberger, O. Industrial Robots for the Automotive Industry. Sae Tech. Pap. Ser. 1996. [Google Scholar] [CrossRef]
  2. Lee, J.; Wu, F.; Zhao, W.; Ghaffari, M.; Liao, L.; Siegel, D. Prognostics and health management design for rotary machinery systems—Reviews, methodology and applications. Mech. Syst. Signal Process. 2014, 42, 314–334. [Google Scholar] [CrossRef]
  3. ISO. DIN ISO 17359:2018-05, Zustandsüberwachung und -Diagnostik von Maschinen—Allgemeine Anleitungen (ISO_17359:2018); Beuth Verlag GmbH: Berlin, Germany, 2018. [Google Scholar] [CrossRef]
  4. Kim, Y.; Park, J.; Na, K.; Yuan, H.; Youn, B.D.; Kang, C.S. Phase-based time domain averaging (PTDA) for fault detection of a gearbox in an industrial robot using vibration signals. Mech. Syst. Signal Process. 2020, 138, 106544. [Google Scholar] [CrossRef]
  5. Jaber, A.A. Design of an Intelligent Embedded System for Condition Monitoring of an Industrial Robot; Springer Theses, Recognizing Outstanding Ph.D. Research; Springer International Publishing: Cham, Switzerland, 2017. [Google Scholar] [CrossRef]
  6. Zhi, H.; Shang, Y. Remote performance evaluation, life prediction and fault diagnosis of RV reducer for industrial robot. J. Phys. Conf. Ser. 2020, 1676, 012212. [Google Scholar] [CrossRef]
  7. Hugo, D.; Benjamin, S. Robot Condition Monitoring: A first Step in Condition Monitoring for Robotic Applications. Master’s Thesis, Lulea University of Technology, Lulea, Sweden, 2017. [Google Scholar]
  8. Karlsson, M.; Hörnqvist, F. Robot Condition Monitoring and Production Simulation. Master’s Thesis, Lulea University of Technology, Lulea, Sweden, 2018. [Google Scholar]
  9. Liu, X.; Wu, X.; Liu, C.; Liu, T. Research on condition monitoring of speed reducer of industrial robot with acoustic emission. Trans. Can. Soc. Mech. Eng. 2016, 40, 1041–1049. [Google Scholar] [CrossRef]
  10. Sun, H.; Zhang, J. Health Monitoring of Strain Wave Gear on Industrial Robots. In Proceedings of the 2019 IEEE 8th Data Driven Control and Learning Systems Conference (DDCLS), Dali, China, 24–27 May 2019. [Google Scholar]
  11. Nentwich, C.; Junker, S.; Reinhart, G. Data-driven Models for Fault Classification and Prediction of Industrial Robots. Procedia CIRP 2020, 93, 1055–1060. [Google Scholar] [CrossRef]
  12. Jahagirdar, A.C.; Gupta, K.K. Cumulative Distribution Sharpness Profiling Based Bearing Fault Diagnosis Framework Under Variable Speed Conditions. IEEE Sens. J. 2021. [Google Scholar] [CrossRef]
  13. Attoui, I.; Boutasseta, N.; Fergani, N. Novel Machinery Monitoring Strategy Based on Time–Frequency Domain Similarity Measurement With Limited Labeled Data. IEEE Trans. Instrum. Meas. 2021, 70, 1–8. [Google Scholar] [CrossRef]
  14. Cocconcelli, M.; Zimroz, R.; Rubini, R.; Bartelmus, W. STFT Based Approach for Ball Bearing Fault Detection in a Varying Speed Motor. In Condition Monitoring of Machinery in Non-Stationary Operations; Fakhfakh, T., Bartelmus, W., Chaari, F., Zimroz, R., Haddar, M., Eds.; Springer: Berlin/Heidelberg, Germany, 2012; Volume 34, pp. 41–50. [Google Scholar] [CrossRef]
  15. Gao, H.; Liang, L.; Chen, X.; Xu, G. Feature extraction and recognition for rolling element bearing fault utilizing short-time Fourier transform and non-negative matrix factorization. Chin. J. Mech. Eng. 2015, 28, 96–105. [Google Scholar] [CrossRef]
  16. Jian, B.L.; Su, X.Y.; Yau, H.T. Bearing Fault Diagnosis Based on Chaotic Dynamic Errors in Key Components. IEEE Access 2021, 9, 53509–53517. [Google Scholar] [CrossRef]
  17. Vippala, S.R.; Bhat, S.; Reddy, A.A. Condition Monitoring of BLDC Motor Using Short Time Fourier Transform. In Proceedings of the 2021 IEEE Second International Conference on Control, Measurement and Instrumentation (CMI), Kolkata, India, 8–10 January 2021; pp. 110–115. [Google Scholar] [CrossRef]
  18. Veerendra, A.S.; Mohamed, M.R.; Punya Sekhar, C. A novel fault—Detection methodology of proposed reduced switch MLI fed induction motor drive using discrete wavelet transforms. Int. Trans. Electr. Energy Syst. 2021, 31. [Google Scholar] [CrossRef]
  19. Defdaf, M.; Berrabah, F.; Chebabhi, A.; Cherif, B.D.E. A new transform discrete wavelet technique based on artificial neural network for induction motor broken rotor bar faults diagnosis. Int. Trans. Electr. Energy Syst. 2021, 31. [Google Scholar] [CrossRef]
  20. Zhu, H.; He, Z.; Wei, J.; Wang, J.; Zhou, H. Bearing Fault Feature Extraction and Fault Diagnosis Method Based on Feature Fusion. Sensors 2021, 21, 2524. [Google Scholar] [CrossRef]
  21. Kotsanidis, K.; Benardos, P. Rolling element bearings fault classification based on feature extraction from acceleration data and artificial neural networks. Iop Conf. Ser. Mater. Sci. Eng. 2021, 1037, 012008. [Google Scholar] [CrossRef]
  22. Skariah, A.; Pradeep, R.; Rejith, R.; Bijudas, C.R. Health monitoring of rolling element bearings using improved wavelet cross spectrum technique and support vector machines. Tribol. Int. 2021, 154, 106650. [Google Scholar] [CrossRef]
  23. Sharma, S.; Tiwari, S.K.; Singh, S. Integrated approach based on flexible analytical wavelet transform and permutation entropy for fault detection in rotary machines. Measurement 2021, 169, 108389. [Google Scholar] [CrossRef]
  24. Elasha, F.; Li, X.; Mba, D.; Ogundare, A.; Ojolo, S. A Novel Condition Indicator for Bearing Fault Detection Within Helicopter Transmission. J. Vib. Eng. Technol. 2021, 9, 215–224. [Google Scholar] [CrossRef]
  25. Mohanraj, T.; Yerchuru, J.; Krishnan, H.; Nithin Aravind, R.S.; Yameni, R. Development of tool condition monitoring system in end milling process using wavelet features and Hoelder’s exponent with machine learning algorithms. Measurement 2021, 173, 108671. [Google Scholar] [CrossRef]
  26. Duan, J.; Duan, J.; Zhou, H.; Zhan, X.; Li, T.; Shi, T. Multi-frequency-band deep CNN model for tool wear prediction. Meas. Sci. Technol. 2021, 32, 065009. [Google Scholar] [CrossRef]
  27. Gómez, M.J.; Marklund, P.; Strombergsson, D.; Castejón, C.; García-Prada, J.C. Analysis of Vibration Signals of Drivetrain Failures in Wind Turbines for Condition Monitoring. Exp. Tech. 2021, 45, 1–12. [Google Scholar] [CrossRef]
  28. Suresh, S.; Naidu, V.P.S. Gearbox Health Condition Monitoring Using DWT Features. In Proceedings of the 6th National Symposium on Rotor Dynamics, Bangalore, India, 2–3 July 2019; Rao, J.S., Arun Kumar, V., Jana, S., Eds.; Springer: Singapore, 2021; Volume 329, pp. 361–374. [Google Scholar]
  29. de Sena, A.P.C.; de Freitas, I.S.; Filho, A.C.L.; Sobrinho, C.A.N. Fuzzy diagnostics for gearbox failures based on induction motor current and wavelet entropy. J. Braz. Soc. Mech. Sci. Eng. 2021, 43. [Google Scholar] [CrossRef]
  30. Yu, D.; Yang, Y.; Cheng, J. Application of time—Frequency entropy method based on Hilbert-Huang transform to gear fault diagnosis. Measurement 2007, 40, 823–830. [Google Scholar] [CrossRef]
  31. Mao, W.; He, J.; Zuo, M.J. Predicting Remaining Useful Life of Rolling Bearings Based on Deep Feature Representation and Transfer Learning. IEEE Trans. Instrum. Meas. 2020, 69, 1594–1608. [Google Scholar] [CrossRef]
  32. Cheng, C.; Ma, G.; Zhang, Y.; Sun, M.; Teng, F.; Ding, H.; Yuan, Y. A Deep Learning-Based Remaining Useful Life Prediction Approach for Bearings. IEEE/ASME Trans. Mechatron. 2020, 25, 1243–1254. [Google Scholar] [CrossRef] [Green Version]
  33. Thakker, H.T.; Dave, V.; Vakharia, V.; Singh, S. Fault Diagnosis of Ball Bearing Using Hilbert Huang Transform and LASSO Feature Ranking Technique. Iop Conf. Ser. Mater. Sci. Eng. 2020, 841, 012006. [Google Scholar] [CrossRef]
  34. Larsen, R.J.; Marx, M.L. An Introduction to Mathematical Statistics and Its Applications, 3rd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2001. [Google Scholar]
  35. Nectoux, P.; Gouriveau, R.; Medjaher, K.; Ramasso, E.; Chebel-Morello, B.; Zerhouni, N.; Varnier, C. PRONOSTIA: An experimental platform for bearings accelerated degradation tests. In Proceedings of the IEEE International Conference on Prognostics and Health Management, PHM’12, Denver, CO, USA, 23–27 September 2012. [Google Scholar]
  36. Nectoux, P.; Gouriveau, R.; Medjaher, K.; Ramasso, E.; Chebel-Morello, B.; Zerhouni, N.; Varnier, C. PHM IEEE 2012 Data Challenge Data Set. Available online: https://github.com/wkzs111/phm-ieee-2012-data-challenge-dataset (accessed on 19 May 2021).
  37. Arun, P.; Lincon, S.A.; Prabhakaran, N. Detection and Characterization of Bearing Faults from the Frequency Domain Features of Vibration. IETE J. Res. 2018, 64, 634–647. [Google Scholar] [CrossRef]
  38. Caesarendra, W.; Tjahjowidodo, T. A Review of Feature Extraction Methods in Vibration-Based Condition Monitoring and Its Application for Degradation Trend Estimation of Low-Speed Slew Bearing. Machines 2017, 5, 21. [Google Scholar] [CrossRef]
  39. Vecer, P.; Kreidl, M.; Smid, R. Condition Indicators for Gearbox Condition Monitoring Systems. Acta Polytech. 2005, 45, 35–42. [Google Scholar] [CrossRef]
  40. Zhu, J.; Nostrand, T.; Spiegel, C.; Morton, B. Survey of condition indicators for condition monitoring systems. In Proceedings of the Annual Conference of the Prognostics and Health Management Society 2014, Fort Worth, TX, USA, 27 September–3 October 2014. [Google Scholar]
  41. Geropp, B. Envelope Analysis—A Signal Analysis Technique for Early Detection and Isolation of Machine Faults. IFAC Proc. Vol. 1997, 30, 977–981. [Google Scholar] [CrossRef]
  42. Lei, Y.; Li, N.; Guo, L.; Li, N.; Yan, T.; Lin, J. Machinery health prognostics: A systematic review from data acquisition to RUL prediction. Mech. Syst. Signal Process. 2018, 104, 799–834. [Google Scholar] [CrossRef]
Figure 1. Process to derive the Z-score-based HI.
Figure 1. Process to derive the Z-score-based HI.
Robotics 10 00080 g001
Figure 2. Saturation of the standard deviation in the time–frequency spectrograms.
Figure 2. Saturation of the standard deviation in the time–frequency spectrograms.
Robotics 10 00080 g002
Figure 3. Comparison of STFT and Z-score spectrograms from healthy and faulty robot gear measurements.
Figure 3. Comparison of STFT and Z-score spectrograms from healthy and faulty robot gear measurements.
Robotics 10 00080 g003
Figure 4. Overall process of the evaluation method.
Figure 4. Overall process of the evaluation method.
Robotics 10 00080 g004
Figure 5. Robot test beds.
Figure 5. Robot test beds.
Robotics 10 00080 g005
Figure 6. Measurement trajectory for the temperature sensitivity analysis.
Figure 6. Measurement trajectory for the temperature sensitivity analysis.
Robotics 10 00080 g006
Figure 7. Faults in the accelerated wear tests, lower image following [7].
Figure 7. Faults in the accelerated wear tests, lower image following [7].
Robotics 10 00080 g007
Figure 8. R2 values for different HIs and bearings from the FEMTO data set based on raw signals.
Figure 8. R2 values for different HIs and bearings from the FEMTO data set based on raw signals.
Robotics 10 00080 g008
Figure 9. Temperature sensitivity for different HIs and robot axes.
Figure 9. Temperature sensitivity for different HIs and robot axes.
Robotics 10 00080 g009
Figure 10. Temperature sensitivity for different DWTRMS-HIs and robot axes.
Figure 10. Temperature sensitivity for different DWTRMS-HIs and robot axes.
Robotics 10 00080 g010
Figure 11. Z-score-HI and RMS-HI for the IRB 7600 experiment.
Figure 11. Z-score-HI and RMS-HI for the IRB 7600 experiment.
Robotics 10 00080 g011
Figure 12. Z-score-HI and RMS-HI for the IRB 6600 experiment.
Figure 12. Z-score-HI and RMS-HI for the IRB 6600 experiment.
Robotics 10 00080 g012
Figure 13. Spectrograms from cold and warm gear measurements.
Figure 13. Spectrograms from cold and warm gear measurements.
Robotics 10 00080 g013
Table 1. Calculated HIs.
Table 1. Calculated HIs.
HI NameHI AbbreviationHI Source
Crest FactorCrF[40]
Dominant FrequencyDomF[37]
Impulse FactorImpF[38]
KurtosisKurt[39]
Margin FactorMarF[38]
MeanMean[40]
MedianMed[40]
Median FrequencyMedF[37]
PeakPeak[39]
Peak to PeakPtP[39]
Root Mean SquareRMS[39]
SkewnessSkew[40]
Spectral CentroidSpC[37]
Spectral FluxSpF[37]
Spectral RolloverSpRO[37]
Spectral EntropySpE[4]
Standard DeviationStd[39]
Discrete Wavelet RMSDWTRMS[19]
Discrete Wavelet Impulse FactorDWTImpF[19]
Discrete Wavelet KurtosisDWTKurt[19,38]
Discrete Wavelet EntropyDWTEntr[4,19]
Discrete Wavelet Decomposition RateDecompRate[18]
Hilbert Huang EntropyHHTEntr[30]
Intrinsic Mode Function RMSIMFRMS[33]
Intrinsice Mode Function Impulse FactorIMFImpF[33]
Intrinsic Mode Function KurtosisIMFKurt[33]
Intrinsic Mode Function EntropyIMFEntr[33]
Time Domain IntegralTDI[14]
Z-scoreZ-score-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nentwich, C.; Reinhart, G. A Method for Health Indicator Evaluation for Condition Monitoring of Industrial Robot Gears. Robotics 2021, 10, 80. https://doi.org/10.3390/robotics10020080

AMA Style

Nentwich C, Reinhart G. A Method for Health Indicator Evaluation for Condition Monitoring of Industrial Robot Gears. Robotics. 2021; 10(2):80. https://doi.org/10.3390/robotics10020080

Chicago/Turabian Style

Nentwich, Corbinian, and Gunther Reinhart. 2021. "A Method for Health Indicator Evaluation for Condition Monitoring of Industrial Robot Gears" Robotics 10, no. 2: 80. https://doi.org/10.3390/robotics10020080

APA Style

Nentwich, C., & Reinhart, G. (2021). A Method for Health Indicator Evaluation for Condition Monitoring of Industrial Robot Gears. Robotics, 10(2), 80. https://doi.org/10.3390/robotics10020080

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop