Developing an Intelligent Fault Diagnosis of MF285 Tractor Gearbox Using Genetic Algorithm and Vibration Signals

This article investigates a fault detection system of MF285 Tractor gearbox empirically. After designing and constructing the laboratory set up, the vibration signals obtained using a Piezoelectric accelerometer which has been installed on the Bearing housings are related to rotary gear number 1 in two directions perpendicular to the shaft and in line with the shaft. The vector data were conducted in three different speeds of shaft 1500, 1000 and 2000 rpm and 130 repetitions were performed for each data vector state to increase the precision of neural network by using more data. Data captured were transformed to frequency domain for analyzing and input to the neural network by Fourier transform. To do neural network analysis, significant features were selected using a genetic algorithm and compatible neural network was designed with data captured. According to the results of the best output mode for each position of the sensor network in 1000, 1500 and 2000 rpm, totally for the six output models, all function parameters for MATLAB Software quality content calculated to evaluate network performance. These experiments showed that the overall mean correlation coefficient of the network to adapt to the mechanism of defect detection and classification system is equal to 99.9%.


Introduction
Rotating machines are the mainstay of several industrial applications worldwide.They can be found in automobile, water industry, wind turbines, agricultural machinery, etc. Failures associated with gears represent the cause of extended outages and are typically caused by gradual deterioration and wear [1].An important objective of condition-based maintenance (CBM) is to determine the optimal time for replacement or overhaul of a machine.The ability to accurately predict the remaining useful life of a machine system is crucial for its operation and can also be used to improve productivity [2].Vibration analysis is a powerful tool in machine fault diagnosis.It is a way to check the status of a machine using a special tool.Troubleshooting machines can be more identified in the future.Fault diagnosis of mechanical or electronic components or systems is a subject of expert system applications.Today, neuro-fuzzy devices are emerged by combining network based on adaptive neural networks and fuzzy models.These adaptive networks, based on fuzzy inference system, are called ANFIS.In such systems the desired output will be achieved by means of human knowledge, in the form of if-then rules for the desired input [2].Neuro Fuzzy system can be applied as a true model based on conventional mathematical tools and component identification.Many researchers [3][4][5][6][7][8][9][10][11][12] have employed various methods to detect weak periodic signals using the sensitivity of system parame-ters.Some of these methods are shock pulse monitoring, crest factor analysis, kurtosis, spectrum analysis, envelope spectrum analysis, wavelet analysis, filtration of signal, etc.Among the researches carried on machine condition monitoring with vibration analysis using fuzzy logic and neuro-fuzzy networks there are rotary machines fault diagnosis, troubleshooting automotive active suspension system and taking into account the non-linear behavior of hydraulic actuator and significantly the impaired controller, an engine control ventilation, fault diagnosis of bearing rotary, fault diagnosis of rolling bearing [13], tapered gearbox fault diagnosis [14].
In recent years, the development of mechanical equipments and maintenance has always been regarded as an essential aspect and different ways to improve safety and reduce vehicle maintenance have always been so important.On the other hand, technological advances in complicated, powerful and expensive machines caused the maintenance programs and systems to be improved [15][16][17].The aim of maintenance based on the status is to gain the symbols of the industrial machinery at a time of being in progress, to provide economic and safety maintenance or to make it stay active [18][19][20].Condition monitoring is performed using different techniques and methods.Each technique is particularly focused on some characteristic symptoms of the status and function of the machine and measures certain parameters that indicate the health status of the machine.Various condition monitoring techniques include vibration analysis, noise analysis, oil analysis and thermal analysis.At the technique of vibration analysis, vibration levels are measured using specific sensors [21,22].To this end, one of the parameters of displacement velocity or acceleration would be measured.The vibration amplitude of each machine monitors the health status of this machine.Furthermore, vibration has high stability and low influence of environmental factors and since it can be an appropriate parameter for condition monitoring methods.In addition, vibration analysis can be used to troubleshoot the machine.Most problems in rotary machines are displayed as vibrating symptoms or changes in the vibration characteristics of the system.Technical advances in mechanical equipment and strong machines, machines and industrial equipment complexity and the high cost of these equipments make it necessary to develop and consider new methods of maintenance including condition monitoring.What is certain is that developing and implementing the methods for machine condition monitoring have had slower rates as a comprehensive way to maintain and troubleshoot the complex mechanical systems during the past years in developing countries which mostly rely on the imports of industrial machinery, construction, and transportation.Due to increase of the accuracy and decrease of the mistakes resulted from human judgment, today, the modeling methods are used for defects classification and detection.Regarding the considerable and increasing development in science, the older modeling and relations, more or less, have lost their influence.Besides, considering the dynamical modeling which works in various parts, caused the use of typical simple models to be complicated.The un-linearity of the problems is also an important and unignorable factor.In this regard, several modeling methods have been proposed, including ant colony, fuzzy logic theory and effective methods like artificial neural networks and genetic algorithm.Each method has its own advantages and disadvantages [23][24][25].In developed agriculture, it is so important to do agricultural works just on time.Time is a vital factor in agricultural systems, which could affect total function [26,27].Therefore, it is necessary for the machines to be available at the time required.Whenever in Tractors of a certain farm a defect is occurred in the gear box, including the important issues that we will be encountered with, too much time and fee would be spent to transmit, open and close and repair the gear box, and since it is a vital part of the tractor and there is no more work without the gear box, so we will have to replace or repair it.However, sometimes a certain human error would cause the maintenance time and cost to be doubled.Since, for years, the experts have followed the discovery of defects in various parts of the tractor in order to reduce the time and cost of economy and raise the work efficiency.This study has been also conducted to make it easier to troubleshoot in gearbox of MF285 Tractor.

Materials and Methods
In this study a smart diagnostic system to make decisions faster than normal due to impair function machines and to increase efficiency and overall lifetime of the device is presented.After adjusting the environment settings, the test was performed in a healthy condition of gears and two defective conditions of the gear worn and the gear broken, and in three iterations of input to the axle along with installing the accelerometer sensor in both directions of Perpendicular to the axis and aligned with the axis on the bearing.For each case, the sample test was repeated 130 times and the data related to each sample were recorded and registered.In total, 2340 experimental data set were recorded as the files on Excel software and it took 20 days to collect data, as a whole.After designing and building profiles needed to position the pieces, all the necessary settings on the device and on the desktop were conducted to start the data analyzing process (Figure 1).
Once the engine rotation been measured using the optical rpm gauge, we provide different rotations in 3 various conditions and in two sensor directions in motor, implementing the dimmer, and record the relevant data.
The gear vibrating signals are measured by the accelerometer installed on bearings on the shaft.The engine has one central axis which is coupled with the gear input and the resulted rotation speed of motor is transformed according to the proportion of the base gear.In order to make failures on gear, we used the gear with 44 gears in condition of gear 1 on the input axis (Figure 2).
In order to collect the vibrating signals, containing the healthy condition and fraction and breakage defects in both perpendicular sensor and axial sensor directions in three motor inputs 1500, 1000 and 2000 and under the loading for 4 seconds, the data were collected in all conditions.In all conditions, every 4 seconds 10,240 samples of voltage were measured in time range of 0.0004.The software used, ARMA, recorded 10,240 voltage data in time range of 4 seconds based on the adjustments in an Excel file.After gaining the data related to healthy condition, the defects and failures considered were applied and the vibrating data were gained, separately.

Signal Processing
Signal processing and the feature extracting part are the most important parts in the condition monitoring procedure, based on which the system feature including being healthy or fractured and type of the system are de-  termined.The current method is to provide the sufficient features, signal transferring from time to frequency domain.For this purpose necessary codes have been transformed in the form of MATLAB software.This program selects Excel files containing the vibration data (based on time/voltage) and filters the data, eliminates the noises and finally calculated the statistical parameters of time and frequency extent, implementing the filter used.Therefore, a fast Fourier function (FFF) of 10,240 points was calculated for each time signal, separately.Also, the power spectral density (PSD) and phase angle (FFT) of vibration signals were calculated using the software MATLAB.

Feature Extraction
The main aim in feature extraction methods is to provide a way or arrangement that could monitor the components hidden in a signal.In this study, 10,240 data sample from time extent formed the data collection.Also, the frequency extent data points, let's say PSD and phase angel FFT, were each 10,240.These signals were used to extract 33 features of healthy and failure conditions.11 features (T1 -T11) were from statistical parameters in time extent and the remaining (P1 -P11) and (A1 -A11) were from statistical parameters in frequency extent (Table 1).
Then, using MATLAB software, the diagrams related to vibration spectrum in frequency extent in healthy and failure conditions were gained due to get analyzed and compared to each other.In order to make the differences in the frequency of normal and defective devices more clear, either axial rotation 1000 rpm and 2000 rpm on the two healthy gear broken at the same diagrams were classified (Figure 3).
Figure 3 shows that in safety condition in both 1000 and 2000 rpm, the amplitude for each diagram remains unchanged during the work and only one or two cases have a maximum that could be justifiable due to the gear involvement.But in unhealthy conditions while the gear with broken gear is replaced with the safe one, the frequency-voltage diagram will immediately have many maximum points which are occurred almost in a regular order.It could be observed that in increasing the rotary axis rotation there are some fixed amplitude increase, while in the broken gear the increase in amplitude is iterated in an unusual way in certain coefficients.The difference resulted from the broken gear has been illustrated in frequency-voltage diagrams against the time-voltage diagram.

Classification of Defects
In MATLAB software using the neural network, the chromosomal genetic with 20 gens was designed and the In which i(j) is a signal for j = 1, 2, …, J. J indicates the number of data points.proper neural network structure was created by determining the number of layers, neurons, epoxies, biases and gens implementing the numbered gens capacity, which studied using various statistical parameters of effi-E.EBRAHIMI ET AL. 156 ciency and error percent of neural network structure.Applying the genetic-neural network hybrid, the necessary classification and calculations were performed to create a network compatible with the test.In this study, one chromosome with 20 gens was designed.The first gen determines the number of inputs.Because, it tried to investigate each input parameter on the model output and if any of these inputs have trivial impact on the model output it would be deleted from the model and in other words the unvalued inputs genetic algorithm will be eliminated (Figure 4).
Gens 2 and 9 illustrate two considerable issues from the neural network structure, one shows the number of layers hidden in network and the other illustrates the number of neurons in each hidden layer, and the amount could be variable from 0 to 12.The transfer functions include: competitive functions, two-valued stimulation threshold, two-valued symmetrical stimulation threshold, sigmoid, converse, positive linear, linear, radial, saturated linear, symmetrical saturated linear, triangular and sigmoid hyperbolic tangent.The 18th gen determined the type of training function, and the amount of this gen could be variable from 0 to 18  algorithm, Levenberg-Marquardt algorithm, the scaled combined gradient algorithm, quasi-Newton back propagation algorithm, monomentum algorithm, combined gradient algorithm Flathcher-Reeves and etc. the 19th gen determines the learning function that have the range of 0 to 14.The learning functions include: bias reduction, weigh or bias of descent gradient, the descent gradient with Monomentum weigh or bias, weigh given, the rule of learning hob with descent weigh, the Instar weigh, the Kohnen weigh, the vector quantization learning weigh LVQ1, LVQ2, outstar weigh, percepteron weigh and bias, percepteron normalized weigh and bias, self-organizer mapping weigh, collective self-organizer mapping weigh and the rule of Vidro-Hauf weigh and bias learning.The 20st gen determines the epoxies to train each network and the range is variable from 10 to 500.Regarding the genetic algorithm adjustments, the best population number for each generation was selected as 50.So, the scale of terminating the genetic algorithm was selected as 100 and the competitive method and two-pointed integration were implemented.All data were included 30 percent of test and 70 percent of training and all the calculations were performed by hybrid genetic-neural network.The optimizing program was used in software MATLAB version Release 2010a, as a whole.In order to classify and analyze the failures and relevant parameters, using the neural network and gaining effective and necessary inputs to be applied in neural-genetic network, the nerve solution was implemented.This software is a sort of developed environment for neural networks, which includes the models based on network design, implementation using advanced learning methods like conjugate gradient, Lonberg-Marquadet and the error propagation during the time.This designing, training and establishing the neural networks (learning under supervision and training without supervision), is used to perform an extent range of tasks like: data analysis, classification, multi-variable regression and time series prediction.Implementing the nerve solution software sensitivity analysis the considered inputs from 55 inputs we had calculated, were gained and the analysis diagram of sensitivity to the average for all data related to the condition of perpendicular to shaft at 1000 rpm was drawn .According to previous analysis using MATLAB for which the results were recorded as the same names in different columns in Excel files, the name of each column is significant for the next calculation.Implementing the name of these variables, the relevant addresses were used as the neural network inputs (Table 2), and then using genetic algorithm in software MATLAB, the best designing structure and it's features were registered (Table 3).

Result
To evaluate the accuracy of the neural network model,    5 and 6), (Table 4).The best network outputs for each sensor condition are divided into three 1000, 1500 and 2000 rpm.There are, totally, 6 output models for which all the system function statistical parameters have been calculated (Table 4).Due to compare two diagnostic positions of each model, general description of the correlation coefficient in each case is given in diagram (Figure 7).The highest correlation coefficients in Model 4 and 5 are about 99% and the lowest one in model ( 6) is about 97%.

Conclusion
In this study a network-compatible device has been designed using artificial neural networks, to detect and classify defects in a Massey Ferguson 285 tractor gear.According to the results, the best network outputs for each sensor mode are divided into three categories in 1000, 1500 and 2000 rpm.Totally, outputs for six cases and the function of all the quality parameters have been computed in order to troubleshoot the gear 1 for gear box in tractor Massey Ferguson 285 in the worst case (aligned with the axis at 2000 rpm) with a correlation coefficient of 97% and the mean square error of 0.000047%.The numbers indicate the high accuracy and low error in trial, respectively.These experim s showed that the overall ent   mean correlation coefficient of the network is equal to 99%, to adapt to the mechanism of defect detection and classification system, and according to the definition of correlation coefficient, this number is near 1 that is the desired mode and the mean error occurred in defect detection is nearly 0. These statistical results show the efficiency and applicability of this method to detect the defects in gearbox.

Figure 4 .
Figure 4.A sample of genetic algorithm structure and procedure to form the best neural network.

5 6
psd-value-y-FFT phase angle-y 2 Perpendicular-1500 rpm fft-values-FFT phase angle-psd-value 3 Perpendicular to the axis at 2000 rpm psd-value-fft-values-fft-values-y-psd-value 4 Aligned with the axis at 1000 rpm fft-values-psd-value-FFT phase angle-y Aligned with the axis at 1500 rpm phase angle-y Aligned with the axis at 2000 rpm FFT phase phase-FFT phase angle-fft-values

Figure 5 .Figure 6 .
Figure 5. Mean-square error network output results at 1000 rpm in the status of perpendicular on the axis.

Figure 7 .
Figure 7.The overall correlation diagram based on tested models.

Table 1 . Features of time and frequency domain.
. The training functions include: the combined gradient algorithm Powell Beale, the combined gradient algorithm Polak-Ribirere, Resilient Back propagation, BFGS algorithm, One Step Secant