A Diagnostic System for Speed-Varying Motor Rotary Faults

This study proposed an intelligent rotary fault diagnostic system for motors. A sensorless rotational speed detectionmethod and an improved dynamic structural neural network are used.Moreover, to increase the convergence speed of training, a terminal attractor method and a hybrid discriminant analysis are also adopted.The proposedmethod can be employed to detect the rotary frequencies of motors with varying speeds and can enhance the discrimination of motor faults. To conduct the experiments, this study used wireless sensor nodes to transmit vibration data and employedMATLAB towrite codes for functionalmodules, including the signal processing, sensorless rotational speed estimation, neural network, and stochastic process control chart. Additionally, Visual Basic software was used to create an integrated human-machine interface.The experimental results regarding the test of equipment faults indicated that the proposed novel diagnostic system can effectively estimate rotational speeds and provide superior ability of motor fault discrimination with fast training convergence.


Introduction
The abnormal vibrations observed during motor rotary faults are typically caused by incorrect shaft connection or friction between rotary components [1][2][3][4]. Traditionally, maintenance personnel had to shut down motors for examination and maintenance. However, the vibration detection method, which involves simply attaching an accelerometer to the motor [5,6] for analysis and measurement, fulfills the demand for real-time, online, and minimally invasive diagnosis.
For experiments regarding vibration faults, time-domain signals are typically converted into frequency-domain signals before analysis. For example, Iorgulescu et al. [7] employed a fast Fourier transform (FFT) to achieve this objective, and Garcia-Perez et al. [8] adopted power spectral density (PSD) to complete a conversion to frequency-domain signals for further analysis.
Previous studies of motor faults have typically focused on measurements of single rotational speeds; however, changes in rotational speed alter fault characteristic frequencies, resulting in incorrect diagnosis. In recent years, scholars have begun considering variations in rotational speeds when diagnosing faults. For instance, Teotrakool et al. [9] used wavelet packet decomposition to diagnose bearing faults under various rotational speeds. Arabaci and Bilgin [10] employed torque speed curves to diagnose rotor bar broken in motors. Also, Blodt et al. [11] conducted a time-frequency analysis to detect torque vibrations under various rotational speeds. These studies further adopted a motor current signature analysis to diagnose motor faults with various rotary frequencies.
However, these studies required the installation of tachometers and the a priori input of rotary frequencies or rotational speeds into diagnostic systems. Inspired by the research done by Ocak and Loparo [12], who employed the physical characteristics of frequency spectra, this study proposed a sensorless method that enables a system to automatically detect the motor rotational speed, thereby reducing the tachometer installation cost.
The development of neural networks in modern intelligent motor fault diagnostic systems is quite mature [13,14]. Such networks possess substantial parallel computing abilities, adaptability, fault tolerance, self-learning abilities, and memory abilities and are, hence, suitable for diagnosing motor faults. The disadvantages of traditional neural networks are that the structure such as the appropriate number of neurons is difficult to know and the training is time consuming. To solve the neuron choice problem, Chuang et al. [15] proposed a dynamic structural neural network (DSNN) that was applied to power quality classification. Still, the training of DSNN is time consuming.
In addition to the above mentioned considerations, this study also incorporated hybrid discriminant analysis (HDA). HDA is a combination of principal component analysis (PCA) and discriminant analysis (DA). These two analysis methods are commonly used in statistics, medicine, and pattern recognition to categorize data and reduce data dimensionality. Researchers [16,17], for example, used PCA in a motor fault analysis and Ayhan et al. [18] combined DA with a neural network for diagnosing the faults caused by rotor bar damage. HDA, which is based on PCA and DA, can be used to identify optimal data classification results. Converted data are then fed into a neural network to increase the speed of learning and discrimination during training.
In this study, motor vibration signals are firstly measured. Subsequently, the vibrational time-domain signals are converted into frequency-domain signals using the fast Fourier transform (FFT) algorithm. The sensorless rotational speed estimation method was used to obtain the current rotatory frequency and rotational speed from spectrograms as a reference for extracting motor fault characteristics. Also, a normalization procedure was used to preprocess the spectral characteristics. To speed up the convergence of training dynamic structural neural network, a hybrid discriminant analysis (HDA) was used to reduce the influence of insignificant characteristics. Different from the conventional motor diagnostic system, this study used the improved dynamic structural network with the terminal attractor method [19] to update the weights (hereafter, called HTAIDSNN) for fast convergence and proposed a novel diagnostic system which is capable of detecting rotary faults in motors at different motor speeds. Furthermore, in practice, to eliminate the influence of noises or outliers, this study also established a dual core diagnostic system that comprised the HTAIDSNN and statistical process control (SPC) techniques [20].

Vibration Characteristics of Motor Rotary Faults and the Normalization
where represents the synchronous speed of the motors, represents the supply frequency, P represents the number of poles, S represents the slip rate, represents the actual speed, and represents the rotary frequency.
After obtaining the rotary frequency of a motor, motor fault can be diagnosed using the rotary frequency and amplitudes related to frequency multipliers. Three most common types of motor rotary faults are bent shafts, imbalanced rotors, and loose screw components. The spectral characteristics of these types of faults are depicted in Figure 1. As indicated in Figure 1, different types of faults correspond to different characteristic frequency spectra; thus, the motor fault type can be determined by examining the specific pattern of a frequency spectrum.
The corresponding spectral characteristics examined during motor fault diagnoses are called fault characteristic parameters. Equation (2) is the mathematical notion of the types of faults, fault characteristics, and fault characteristic parameters where ( ) denotes the fault type, represents the fault characteristics, ,1 , ,2 , . . . , , , . . ., and , represent the fault characteristic parameters of the th-segment signal, represents the th characteristic frequency, and denotes the number of characteristic frequencies.
When faults occur, the amplitudes of specific frequency multipliers change accordingly. The fault characteristic parameters used in this study are defined as 0.5 , 1 , 1.5 , 2 , 2.5 , 3 , and 3.5 . As shown in (3), a vibration eigen-matrix of a healthy motor is compiled using fault characteristics based on the -segment vibrational data collected from healthy motors collected in this study: ] . ( The means and standard deviations of amplitudes of these characteristic frequencies are then calculated by using The amplitude of each characteristic frequency of thesegment signal is normalized by ] . When normalizing the vibration eigen-matrices of different fault types, the amplitude means and standard deviations of the characteristic frequencies for healthy motors are employed as the basis of normalization. Thus, all the normalization processes are performed using the same basis. As shown in (8), (9), and (10), the normalized fault eigenmatrices , , and are obtained by normalizing the fault eigen-matrices of bent shafts, imbalanced rotors, and loose components, respectively. The compiled normalized fault eigen-matrices are then used as training data for the neural networks. Fault eigen-matrices with a greater number of segments provide more training data, thereby increasing the discrimination ability of the networks. Consider the following: ] .

Sensorless Rotational Speed Estimation
Traditionally, an additional tachometer or an encoder is employed for measuring motor rotational speed. To reduce the measurement cost, this study proposes a sensorless rotational speed estimation method that determined motor rotational speed by analyzing the vibrational frequencydomain signals of motors. Previous literature [12] has indicated that the harmonics of rotary frequency exhibit greater amplitudes in the vibration frequency spectra of motors. Considering this finding, this study first identifies the local maximums in relevant spectrograms and their corresponding frequencies.
After identifying the frequencies that corresponded to the local maximums, an indicator with an initial value of 0 is assigned to each frequency shown in Figure 2. If the corresponding frequency for one of the local maximums has a multiplicative relationship with a larger frequency, 1 is added to . If two relationships are determined, 2 is added to the value, and so on. Because the actual frequency values do not exhibit exact multiplicative relationships, errors less than 3 Hz are allowed in this study. Finally, the frequency corresponding to the greatest ( = 2) is selected as the estimated rotatory frequency. For example, a frequency of 30 Hz (the corresponding rotational speed is 1800 rpm) is lower than the frequencies of 43.8, 60.7, and 120.2 Hz, with the frequencies of 60.7 and 120.2 Hz being two and four times (second and fourth harmonic generation) higher than that of 30 Hz; thus, is set as 2. If the frequency is 60.7 Hz, 120.2 Hz is a frequency that doubles the original value (second harmonic); therefore, is set as 1. The frequency with the greatest is selected as the desired rotatory frequency as indicated in Figure 3.
To verify the speed estimation method using the above mentioned concept, a four-pole, one-hp, three-phase induction motor manufactured by the TECO Electric & Machinery Company (Taiwan) is adopted for this experiment. The experiment performed ten times of measurement and estimation for rotational speeds from 900 rpm to 1800 rpm. Table 1 summarizes the statistics of the estimated motor rotational speeds. The standard deviations are smaller than 4 rpm for each rotational speed. Moreover, the estimation errors of rotational speeds are smaller than 2 rpm and thus the accuracy of the proposed sensorless rotational speed estimation method is confirmed.     [15]. Figure 4 shows a schematic diagram of the DSNN framework on a pseudo-3D plane. The network is capable of adjusting the number of neurons in the hidden layer, which can be calculated using the following equation: where ℎ represents the number of neurons in the hidden layer, represents the length of the hidden layer, represents the height of the hidden layer, represents the width of the hidden layer, and represents the number to which the neurons can be expanded, which can prevent neurons from increasing without limitation. Assume that there are neurons in the input layer, neurons in the output layer, and neurons in the th hidden layer, for = 1, 2, . . . , . The total connection in Figure 4 indicates that all neurons in hidden layer are interconnected. Equation (12) shows the output vectors of the neurons in the first hidden layer. The output vectors for neurons in the other hidden layers are presented in (13), and the output vectors of the output layer are expressed in (14) as follows: where ℎ 1 , ℎ , and are the outputs of the neurons in the first and th hidden layer and the output layer, respectively.

The Terminal Attractor Method.
The error function is defined as (15), where is the expected output and is the actual output vector. The number of neurons in the output layer is represented by = 1, 2, . . . , . Therefore, During the learning process for a neural network, the steepest descent method is the option most commonly employed for reducing the error function in (15)   in a local optimal solution. This leads to the fact that the error cannot converge to zero within a limited time. To resolve this issue, Wang and Hsu incorporated heuristic terminal attractors into back-propagation neural networks in 1991 [19]. Assume that is the aggregated weight vector and let the derivative of satisfy where 0 < < 1. Consequently, the derivative of error function becomes which determines a terminal attractor. Using the concept, the update law of the weight vector in this paper is given as where is the learning rate. Hereafter, the proposed improved dynamic structural neural network is abbreviated as HTAID-SNN (heuristic terminal attractor improved dynamic system neural networks) and the learning rate is set to 0.9.

Applying Hybrid Determinant Analysis to HTAIDSNN.
To identify the eigen-matrix of principal components, ( ) , which possesses various classes of data, is considered as where denotes the class number in sequence, denotes the times of collection in the th class, denotes the number of factors and characteristics considered, and The covariance matrix for the data is presented using (22), where represents the transpose or transposition symbol.

Mathematical Problems in Engineering
Finally, through (23), the transformation matrix PCA is obtained following PCA: Using (19), (20), and (21), this study independently establishes matrices for within-class variance (24) and betweenclass variance (25). The ratio of to is larger, and the discrimination is better. Finally, using the Fisher DA method, the transformation matrix FDA is obtained by using (26) Although between-class variances can increase during performing PCA, within-class variances also increase. Therefore, DA appears to be a more appropriate method for relevant computations. However, for data classes with a limited sample size or uneven sample distribution, PCA is superior to DA. To enhance the discrimination of data analysis, HDA [20] is used to obtain the optimal transformation matrix for given data. Equation (27) shows the transformation matrix of HDA: where represents the identity matrix and and represent parameters that are adjustable between 0 and 1. Beginning at = 0 and = 0, the parameters can increase by 0.1 to obtain the optimal HDA transformation matrix HDA . Equation (28) shows the new matrix HDA ( ) obtained after the transformation: To perform the HDA, firstly the seven frequencymultiplier amplitudes considered in Section 2 are measured 20 times for various motor faults. The normalized fault eigenmatrices are then compiled using (7), (8), (9), and (10). These matrices are substituted into (19), where the total number of classes is 4 and the times or frequency for every class is equal to 20. The covariance matrix , within-class variance matrix , and between-class variance matrix are obtained by substituting them into HDA. The adjustable parameters and are set to 0.4 and 0.7, respectively. Finally, a new fault eigen-matrix is obtained with (28). Generally, when using hybrid determinants or principal components,  the distribution of the first two components is analyzed [20], because the components are the most suitable for representing data traits. Figure 5 shows the distribution of the first two hybrid discriminant components. Figure 5 shows that HDA provides comprehensive data classification and small within-class distance or variance and enhanced data discrimination. When clustering or grouping for fault classification is used, the relationship between only two component values is typically considered. This study combined the HTAIDSNN with the fault eigen-matrices converted by HDA. Based on the compilation of the fault eigenmatrices, this study established a neural network training matrix, as shown in Figure 6, where hd denotes the th hybrid discriminant or discrimination and = 1, 2, . . . , 7. hd is used as the input neuron value. Regarding the parameters of neural networks, the HTAIDSNN possessed seven neurons in the input layer and four neurons in the output layer. Equation (11) shows the framework parameters of the hidden layer, where = 4, = 5, and = 2. Furthermore, the number to which the neurons could be expanded is = 2; therefore, the initial number of neurons in each hidden layer is 10. The initial weight and bias or threshold values are both random numbers ranging between −1 and 1. The terminal criterion for concluding training is 5000 iterations or a minimum error smaller than 10 −6 . These criteria prevented excessive convergence times when the system is constrained in local optimal solutions. Figure 7 represents a comparison of the error convergence curves for the data processed by the HTAIDSNN and  HTAIDSNN  · · · · · · · · · · · · · · · b 1,1 b   HTAIDSNN with HDA has superior convergence performance. The improvement of the convergence rate is mainly due to the specific classification information provided by HDA. Table 2 shows the statistics of fault classification using the HTAIDSNN-based methods. From the table, the performance of the HTAIDSNN with HDA is superior to the other two methods.

Description of Experimental Equipment.
The equipment employed in this study to conduct the motor vibration experiment included a four-pole, one-hp, three-phase induction motor manufactured by the TECO Electric & Machinery Company (Taiwan). An A7200-MA VFD (TECO Electric & Machinery Company) was also installed to change the rotational speed, and the source frequency was 60 Hz. Figure 8 shows the wireless sensor nodes that are attached to the motor in the motor vibration experiment. An MMA7260 three-axis accelerometer was installed on a wireless sensor node (employed for transmitting vibration data) to measure vibration displacement. During the experiment, screw nuts were installed on the motor flying wheel as unbalanced weights to generate the vibration fault signals that had frequency characteristics similar to that of an imbalanced rotor or bent shaft. The screws on the motor equipment stand were adjusted to emulate   vibration fault signals resulting from loose components as shown in Figure 8.

Description of System
Procedures. As depicted in Figure 9, this study initially retrieves the motor vibration signals using wireless sensor nodes and obtains the current rotary frequency through the sensorless rotational speed estimation method. Subsequently, the fault characteristics are extracted from the frequency-multiplier amplitudes of the rotary frequency, and the characteristic data are normalized. Finally, after HDA conversion, the neural network feeds the compiled hybrid discriminant fault eigen-matrices into the HTAIDSNN for training of learning and memorizing for the relationships between the fault characteristics and fault classes.
In addition to the HTAIDSNN mentioned previously, control charts for stochastic process control (SPC) are included for fault diagnosis. The result of this combination is a dual core diagnostic system that employed parallel processing [21].
Mathematical Problems in Engineering  The SPC technique distinguishes between random outliers and abnormal statuses or states regarding the quality of monitored products to provide warnings concerning abnormal trends in the production process. Since the workplace is in general a noisy environment, the SPC charting used in this work is to monitor the shift and drift of the characteristic frequencies and rule out the abruptly environmental noise influence. The rule adopted in this paper is as follows: if there is only one point out of control limit, a warning is assigned. An abnormality in this study is defined as the situations in which two of three continuous or consecutive points are larger than control limit.
In the dual core diagnostic system, if the neural network and control chart are free from abnormalities, the system generates an LV1 signal to indicate the normal state. If either the neural network or the control chart generates an abnormal signal, an LV2 signal is generated to indicate that the motor may be faulty and require reconfirmation or stay alerted. If both are abnormal, an LV3 signal is generated to indicate that the motor is faulty and should be immediately repaired by technicians.

Description of the Human-Machine
Interface. The human-machine interface designed in this study primarily includes a motor diagnostic system. Figure 10 shows the motor monitoring interface designed for the imbalanced rotor experiment conducted at a rotary frequency of 30 Hz. This interface comprises a time-domain graph, a spectrogram, and a rotational speed estimation graph for the vibration signal. The frequency spectrum displays the peaks of the recent vibration signals at 0.5 , 1 , 1.5 , 2 , 2.5 , 3 , and 3.5 . Finally, the speed estimation graph displays the five most recent rotational speeds. The tolerable rotational speed change is 50 rpm. If the change of speed is larger than 50 rpm, the rotational frequency and speed are recalculated accordingly. To exclude the influences of noise or disturbance, the changes on the rotational speed are only displayed after the occurrence of the changes had been verified. Figure 11 shows the interface designed for the diagnostic system used in the imbalanced rotor experiment conducted at a rotary frequency of 30 Hz. This interface includes a neural network detection and Shewhart control chart. Both the neural network and the control chart in the figure detect an abnormality, and the overall diagnosis is marked in red (LV3). Figure 12 shows vibration signal, frequency spectrum, and the estimated rotational speed of a bent-shaft motor. Figure 13 indicates the diagnosis information of the motor. In the figure, the neural network detects an abnormality, and the Shewhart chart points out that there is a point out of control limit. By the rule defined in this work, this situation is acceptable unless the consecutive point is also out of control limit. Consequently, the overall diagnosis is classified as LV2 marked in yellow.

Conclusion
This study proposes a novel intelligent diagnostic system for speed-varying motors, which improves the convergence rate of traditional dynamic structural neural network training by introducing the terminal attractor and the hybrid discriminant analysis. The implemented system integrates the required functions: sensorless speed estimation, neural network-based diagnosis, and SPC charting, with a humanmachine interface to perform rotary fault diagnosis. The simulation results indicate that the proposed system is equipped with quick convergence rate for neural network training. Also, motor rotary faults are detected effectively by the developed diagnostic function and statistics process control chart.