A Novel Personalized Diagnosis Methodology Using Numerical Simulation and an Intelligent Method to Detect Faults in a Shaft

Personalized medicine is a hot topic to develop a medical procedure for healthcare. Motivated by molecular dynamics simulation-based personalized medicine, we propose a novel numerical simulation-based personalized diagnosis methodology and explain the fundamental procedures. As an example, a personalized fault diagnosis method is developed using the finite element method (FEM), wavelet packet transform (WPT) and support vector machine (SVM) to detect faults in a shaft. The shaft unbalance, misalignment, rub-impact and the combination of rub-impact and unbalance are investigated using the present method. The method includes three steps. In the first step, Theil’s inequality coefficient (TIC)-based FE model updating technique is employed to determine the boundary conditions, and the fault-induced FE model of the faulty shaft is constructed. Further, the vibration signals of the faulty shaft are obtained using numerical simulation. In the second step, WPT is employed to decompose the vibration signal into several signal components. Specific time-domain feature parameters of all of the signal components are calculated to generate the training samples to train the SVM. Finally, the measured vibration signal and its components decomposed by WPT serve as a test sample to the trained SVM. The fault types are finally determined. In the simulation of a simple shaft, the classification accuracy rates of unbalance, misalignment, rub-impact and the combination of rub-impact and unbalance are 93%, 95%, 89% and 91%, respectively, whereas in the experimental investigations, these decreased to 82%, 87%, 73% and 79%. In order to increase the fault diagnosis precision and general applicability, further works are continuously improving the personalized diagnosis methodology and the corresponding specific methods.


Introduction
Rotating machinery has been widely used in transmission machinery, and its running state directly affects the entire machine, including reliability and stability.Thus, the fault diagnosis of the rotating machinery is increasingly having attention paid.Over the past several decades, new fault detection methods [1][2][3][4][5][6][7] and new criteria [8][9][10] were proposed, which greatly enriched the field of fault diagnosis.He et al. [11] proposed a data mining-based diagnostic system using acoustic emission (AE) signals to detect faults in full ceramic bearings.They further developed a series of methods to identify faults in plastic bearings [12], gearbox tooth cut faults [13], planetary gearbox faults [14], etc.
For the fault detection of the shaft, various researchers have reported feasible methods.Mohammed et al. [15] investigated a novel route using artificial neural networks (ANN) and power spectral density (PSD) to detect the cracks in a rotating shaft.As is well known, intelligent methods play an important role in detecting faults in mechanical components, such as genetic algorithms (GAs), neural networks (NNs), support vector machine (SVM), etc., and those approaches have been applied to detect faults using experimental samples of all kinds of faults [7].However, for the lack of suitable training samples that represent every possible fault in all kinds of practical mechanical systems, the above intelligent techniques have not been agreeable applied to detect faults in mechanical systems.
In order to understand the fault effects in depth, many researchers proposed numerical simulation for fault diagnosis.Torkaman et al. [16] adopted the three-dimensional time-step finite element method (3D-TSFEM) to simulate the fault of a switched reluctance motor (SRM), and the new diagnosis index associated with power losses was obtained in their results.To detect faults online in nonlinear continuous systems, Bregon et al. [17] used the simulation and state observer models to obtain the final fault diagnosis results.Gong et al. [18] investigated rotor dynamic stability analysis on the hydraulic turbine by numerical simulation and also demonstrated how each factor affects the dynamic character of the labyrinth system.Xiang et al. [19,20] carried out two kinds of wavelet-based numerical simulation models to calculate the dynamic responses of a shaft.Tannous et al. [21] performed rotor-stator contact simulations using 3D numerical modeling.Baccarini et al. [22] proposed a dynamic numerical model to simulate mechanical faults in induction machines, and the simulation results confirmed the validity of the model.For numerical simulation, on the one hand, less time and equipment are used to get a large amount of experimental data, especially for those that are very difficult to conduct in a real machine for carrying out the experimental investigation.On the other hand, this allows the investigator to determine the fault samples for all types of faults under the complex running conditions.
In order to perform the intelligent methods, the selection of feature vectors is one of the key problems.Recently, time domain or time frequency domain feature indexes were adopted to create feature vectors.Chen et al. [23] proposed an approach that has six indexes (e.g., mean value, root mean square, standard deviation, skewness, kurtosis and shape indicator) in the time domain and two indexes (e.g., mean frequency and standard deviation frequency) in the frequency domain to construct the feature sets to train the SVM.Liu et al. [24] used a hybrid time frequency analysis method to get the feature information of gear faults.
Wavelet packet transform (WPT) is a useful tool in handling non-stationary signals to detect faults in mechanical systems [25].As a signal processing tool, WPT decomposes the original signal into multi-layers, and the same frequency bandwidth is offered in each layer.Xiang et al. [26] propose a new method based on the wavelet transform for the detection of damages in plate-like and shell-like structures.Diego and Barros [27] employed WPT to analyze the magnitude of the harmonic distortion and the corresponding sub-harmonics in power systems.Pan et al. [28] proposed a novel approach based on improved WPT and support vector data description (SVDD) to assess bearing performance degradation.
To obtain the features of the ignition pattern, Vong and Wong [29] used WPT to decompose the engine ignition signal and then employed multi-class least squares SVM to classify the fault types.It is worth pointing out the coefficients of WPT and the corresponding single branch reconstruction signals on all of the decomposing levels, which play a key role in singularity detection and feature extraction.
As a novel pattern recognition approach, support vector machine (SVM), was proposed in the 1990s by Cortes and Vapnik, and now, it has been widely used in fault classification [30].The basic theory for SVM is structural risk minimization (SRM) with small training samples, which are suitable for the classification and regression of the nonlinear and high dimensional problems.However, the lack of faulty training samples has greatly limited the fault classification using SVM in real-world applications.
Recently, personalized medicine has become a hot topic to develop a medical procedure that separates patients into different groups, with medical decisions, practices, interventions and/or products being tailored to the individual patient based on his/her predicted response or risk of disease [31].It is also a new direction for placing the patient at the center of healthcare, and the theoretical basis is molecular dynamics simulation from person to person [32].Motivated by molecular dynamics simulation-based personalized medicine, the fault detection strategy for mechanical systems might be developed for the personalized diagnosis era.The most probable tools might be numerical simulation, the big data technique or a combination of the two.
From the above descriptions, it might be concluded that if the limitation of faulty training samples can be inexpensively overcome by numerical simulation using the dynamic model of mechanical systems, the intelligent methods, such as SVM, GAs and NNs, will be directly applied to detect faults for machines under all kinds of working conditions.In this paper, we present the basic idea of a new personalized fault diagnosis methodology focused on the fault diagnosis of shafts.Based on the present methodology, this paper is arranged as follows: The basic idea and the procedures of the personalized diagnosis methodology are presented in Section 2. In Section 3, finite element (FE) models of shafts are constructed; an example is presented; and the shaft unbalance, misalignment, rub-impact and the combination of rub-impact and unbalance are considered.To further apply the present method to real-world usage, the experimental investigations and further works are discussed in Section 4. Finally, the conclusion of this study is given in Section 5.

The Basic Idea for the Personalized Diagnosis Methodology
Personalized medicine is a kind of accurate diagnosis technology based on individual differences with the development of molecular biology in recent years [33], which provides a new idea for condition monitoring and fault diagnosis in mechanical systems.Recently, the intelligent diagnosis methods have become one of the hot and promising research areas [34].The valuable application of intelligent methods has reflected and provided the realization of personalized fault diagnosis in mechanical systems.The early development of mechanical diagnosis and condition monitoring was based on the faulty samples extracted from actual faults in mechanical systems.However, for the complex mechanical system and complicated fault mechanism, there are too many complexity vibration signals to be analyzed.Furthermore, in many faulty cases, it is very difficult to obtain all of the faulty samples when the machine is running.Therefore, the lack of faulty samples will lead to the failure to detect faults in real-world mechanical systems.
Recently, it has been possible to conduct a simulation model of a mechanical system with faults, such as shafts in a rotor-bearing system, bearings and gears in a gearbox, etc., to obtain the faulty samples with various faults of different types using numerical simulation.
Based on the above technique development, the personalized diagnosis methodology is proposed, and the fundamental procedure is shown in Figure 1.In order to realize personalized diagnosis, the present methodology combines the numerical simulation model and intelligent diagnosis methods to detect faults in mechanical systems.Numerical simulation plays a key role in the present methodology to obtain the simulated faulty samples to replace the measured ones.Therefore, the damage location prediction procedure is as follows: (1) Construct the FE model of mechanical systems with faults: In the simulation, the faulty model (e.g., unbalance, misalignment, rub-impact, damages, etc.) and the dynamic model of the mechanical system are constructed using finite element method (FEM).To decrease the difference between the measured vibration signals and the computed (by the FE model) ones of the same mechanical systems, the FE model updating technique using the experimental investigation of the intact mechanical systems [35] should be applied to reduce the difference.After the updating process, the resulting FE model of the mechanical systems with faults can be used to calculate the simulation vibration signals.
(2) Obtain the faulty samples: The FE model of the mechanical systems with faults is calculated for the vibration signal in the time domain.The WPT is then applied to decompose the original vibration signals into multi-layers, and the corresponding sub-signals (components) are finally obtained.To generate the faulty training samples, the time domain indexes, such as the standard deviation (SD), peak, kurtosis, clearance factor (CF), impulse factor (IF), etc. (summarized in Table 1), are calculated from the sub-signals.Finally, the faulty training samples formed by the time domain indexes severe as inputs to train the intelligent diagnosis models using intelligent diagnosis methods.

Feature Equation
Standard deviation x std 1 x is the data; N is the number of data points; x m is the mean value of x; σ s is the standard deviation.
(3) The personalized fault diagnosis methodology: The measured vibration signals of the real-world mechanical systems and the corresponding time domain indexes are calculated using the same procedures described in Step (2).Then, the time domain indexes of the real-world mechanical systems are input into the trained intelligent diagnosis models.The outputs indicate the healthy condition or the types of faults.It is worth pointing out that if the precision of the FE model of the mechanical systems can be guaranteed, the present fault diagnosis methodology could be applied to any type of real-world mechanical system for any working conditions.Therefore, the present methodology could be called the personalized fault diagnosis methodology.(2) Obtain the faulty samples: The FE model of the mechanical systems with faults is calculated for the vibration signal in the time domain.The WPT is then applied to decompose the original vibration signals into multi-layers, and the corresponding sub-signals (components) are finally obtained.To generate the faulty training samples, the time domain indexes, such as the standard deviation (SD), peak, kurtosis, clearance factor (CF), impulse factor (IF), etc. (summarized in Table 1), are calculated from the sub-signals.Finally, the faulty training samples formed by the time domain indexes severe as inputs to train the intelligent diagnosis models using intelligent diagnosis methods.

Feature Equation
Standard deviation std x n 1 x is the data; N is the number of data points; xm is the mean value of x; s σ is the standard deviation.
(3) The personalized fault diagnosis methodology: The measured vibration signals of the real-world mechanical systems and the corresponding time domain indexes are calculated using the same procedures described in Step (2).Then, the time domain indexes of the real-world mechanical systems are input into the trained intelligent diagnosis models.The outputs indicate the healthy condition or the types of faults.It is worth pointing out that if the precision of the FE model of the mechanical systems can be guaranteed, the present fault diagnosis methodology could be applied to any type of real-world mechanical system for any working conditions.Therefore, the present methodology could be called the personalized fault diagnosis methodology.

The Basic Principle of WPT
Wavelet transform (WT) is the process of the decomposition and reconstruction of the original signal both in the time and frequency spaces, which is widely used to extract the fault feature.However, at the high frequency bands, WT fails to give a good performance.Wavelet packet transforms (WPT) are the better choice to decompose the high-frequency bands.The typical binary tree structure is shown in Figure 2. Index (i, j) represents the decomposition signal, where i is the layer and j (j = 2i − 1) is the node in the i-th layer, respectively.For instance, (0,0) represents the primary signal; (1,0) is the low frequency component in the first layer; and the corresponding high frequency component is (1,1).

The Basic Principle of WPT
Wavelet transform (WT) is the process of the decomposition and reconstruction of the original signal both in the time and frequency spaces, which is widely used to extract the fault feature.However, at the high frequency bands, WT fails to give a good performance.Wavelet packet transforms (WPT) are the better choice to decompose the high-frequency bands.The typical binary tree structure is shown in Figure 2. Index (i, j) represents the decomposition signal, where i is the layer and j (j = 2i − 1) is the node in the i-th layer, respectively.For instance, (0,0) represents the primary signal; (1,0) is the low frequency component in the first layer; and the corresponding high frequency component is (1,1).

A Brief Review of SVM
Consider a training set S: where xiR l , yi{−1,+1} and l is the number of samples.The aim of SVM is to determine an optimal hyper plane for separating one from the others by using the training dataset.To get the ideal hyper plane, the dual optimization problem is often mentioned in SVM as: in which α i is the Lagrange multiplier coefficient obtained by dealing with the dual optimization in the process of the SVM training; K(xi, yj) is referred to as the kernel function; C > 0 is the error penalty parameter.There are many forms of the kernel function.The radial basis function (RBF) kernel is employed in this paper because of the highest accuracy rate of classification.How to choose the tradeoff parameters p1 (the width of RBF) and C is a difficult task in the application of SVM.A possible way is suggested by Xiang et al. [36] in which the simulation investigation should proceed first to obtain the relative best parameters p1 and C, and hence, these parameters were employed for the realworld structures.

A Brief Review of SVM
Consider a training set S: where x i ∈R l , y i ∈{−1,+1} and l is the number of samples.The aim of SVM is to determine an optimal hyper plane for separating one from the others by using the training dataset.To get the ideal hyper plane, the dual optimization problem is often mentioned in SVM as: in which α i is the Lagrange multiplier coefficient obtained by dealing with the dual optimization in the process of the SVM training; K(x i , y j ) is referred to as the kernel function; C > 0 is the error penalty parameter.There are many forms of the kernel function.The radial basis function (RBF) kernel is employed in this paper because of the highest accuracy rate of classification.How to choose the tradeoff parameters p 1 (the width of RBF) and C is a difficult task in the application of SVM.A possible way is suggested by Xiang et al. [36] in which the simulation investigation should proceed first to obtain the relative best parameters p 1 and C, and hence, these parameters were employed for the real-world structures.

The Proposed Method to Detect Faults in Shafts
In this section, the basic idea for the personalized diagnosis methodology is applied to detect faults in shafts.Both WPT and SVM are employed to generate the personalized diagnosis method, and the flowchart of the proposed fault diagnosis method is shown in Figure 3. Three main steps are given in detail as follow: 1 Step (1) Build the FE model and obtain the simulation signals: The numerical model of the shaft can be constructed by commercial Finite element analysis (FEA) software ANSYS (version 10.0, ANSYS Inc., Pittsburgh, PA, USA, 2011), and the simulation signals will be obtained by applying the specific loads.The FE model updating technique (e.g., Theil's inequality coefficient, TIC [37]) is used to determine the stiffness and damping coefficients of the boundary condition (as shown in Figure 4) using the measured signals of the intact mechanical systems to reduce the difference between the simulations and real-world systems.If the difference between the two signals is within the agreeable range (in engineering applications, lower than five percent is commonly accepted), the FE model of shafts will well represent the real-world ones.Then, the faulty models will be inserted into the FE model to simulate the dynamic response of the shaft with faults.(2) Decompose and reconstruct the vibration signal to train SVM: WPT is employed to decompose and reconstruct the vibration signal.In the present work, the Daubechies wavelet is employed as the basis function to decompose the raw signal at a certain level to obtain a reasonable number of frequency bands.For each frequency band in the maximum layer, the single branch reconstruction technique is performed by wavelet packet reconstruction to obtain the signals at each frequency band.In order to detect the shaft faults, five time domain feature parameters, such as SD, peak, kurtosis, CF and IF, as shown in Table 1, are calculated from the signals at each frequency band to generate the feature vector as training samples to obtain the trained SVM.(2) Decompose and reconstruct the vibration signal to train SVM: WPT is employed to decompose and reconstruct the vibration signal.In the present work, the Daubechies wavelet is employed as the basis function to decompose the raw signal at a certain level to obtain a reasonable number of frequency bands.For each frequency band in the maximum layer, the single branch reconstruction technique is performed by wavelet packet reconstruction to obtain the signals at each frequency band.In order to detect the shaft faults, five time domain feature parameters, such as SD, peak, kurtosis, CF and IF, as shown in Table 1, are calculated from the signals at each frequency band to generate the feature vector as training samples to obtain the trained SVM.
(3) Detect the faults using the measured vibration signals and the trained SVM: The signal of the faulty shaft is measured from the real-world rotor-bearing system.The signal will be decomposed and reconstructed by WPT at the last layer to form the sub-signals.The five time domain feature parameters of each sub-signal will further be calculated, and then, the feature vectors serve as a testing sample for the trained SVM.The fault pattern will finally be obtained.

The Simulation Model for a Shaft with Different Faults
In this section, the one-dimensional finite element model is constructed using commercial FEA software ANSYS.The geometry of the shaft is shown in Figure 4.In the FE model, the beam element (BEAM188), combine element (COMBINE14) and mass element (MASS21) are employed, and the whole constraints of the outer side of the bearing are conducted in the model.The material properties are: Young's modulus E = 2.1 × 10 11 Pa, Poisson's ratio µ = 0.3, material density ρ = 7.8 × 10 3 kg/m 3 , respectively.It is noted that this paper introduces the dynamic response simulation for the no-load operation of a shaft with rotary speed provided by the drive system.However, for a running shaft with full load, both the rotary speed and torsion will be considered.

Case Investigations Using Numerical Simulation
Four different faults including shaft unbalance, misalignment, rub-impact and the combination of rub-impact and unbalance are considered in this section.

Unbalance and Misalignment
According to the theory of the fault mechanism, the centrifugal forces F c (x) and F c (y) of shaft unbalance are obtained from the following equations in the x-direction and y-direction as [38]: where m, e and ω refer to the eccentric mass, eccentric distance and angular velocity, respectively.For the misalignment fault of the shaft, the exciting forces F m (x) and F m (y) occurring on the coupling along the x-direction and y-direction are represented by [38]: where m c is the coupling mass and δ is the parallel misalignment between the shafts.∆L is the installation distance between the two half couplings.α is the angular misalignment between the shafts.Suppose the shaft speed Ω = 1800 rpm (randomly selected); the typical finite element simulation signal with unbalance and misalignment faults is shown in Figure 5a,b, respectively.We point out that the length of the signal is 8000 data points, and the first sub-data points (400 data points) are shown in Figure 5.To simulate the unbalance fault, let the eccentric mass m = 0.002 kg, the eccentric distance e = 0.033 m and the corresponding centrifugal forces F c (x) = 2.34cos(188t)N and F c (y) = 2.34sin(188t)N, calculated using Equation (3), be applied to the disc mass element.4), be applied to the end of shaft shown in Figure 4.

Rub-Impact and the Combination of Rub-Impact and Unbalance
Rub-impact occurs in a rotating rotor-bearing system with radial clearance and the process of contacting and separation between the rotor and stator.During the rotating process, the transient impact forces will excite when the rub-impact fault occurs on the rotor-bearing system, and the shape of the force impulse is modeled as a triangular form according to [39].Therefore, in this section, the exciting force of the rub-impact in the x-direction and the y-direction is applied with an amplitude of 30 and 60 N, respectively, shown in Figure 6 near the end of shaft, and the corresponding simulation signal with the rub-impact fault is shown in Figure 5c.To simulate the misalignment fault, let the coupling mass m = 0.05 kg, the parallel misalignment distance δ = 0.001 m, the installation distance between the two half couplings ∆L = 0.002 m, the angular misalignment α = 10 • and the corresponding exciting forces F m (x) = 4.79sin(376t)N and F m (y) = 4.79cos(376t)N, calculated using Equation ( 4), be applied to the end of shaft shown in Figure 4.

Rub-Impact and the Combination of Rub-Impact and Unbalance
Rub-impact occurs in a rotating rotor-bearing system with radial clearance and the process of contacting and separation between the rotor and stator.During the rotating process, the transient impact forces will excite when the rub-impact fault occurs on the rotor-bearing system, and the shape of the force impulse is modeled as a triangular form according to [39].Therefore, in this section, the exciting force of the rub-impact in the x-direction and the y-direction is applied with an amplitude of 30 and 60 N, respectively, shown in Figure 6 near the end of shaft, and the corresponding simulation signal with the rub-impact fault is shown in Figure 5c.
For the fault, the combination of rub-impact and unbalance, the two kinds of forces are applied at the same time using the same finite element model of unbalance and rub-impact simulations.The simulation signal with the combination of rub-impact and unbalance fault is shown in Figure 5d.

Fault Detection Using the Present Method
In the real-world practice, certain additive noises, which might possibly be interferences from the environment and the equipment of mechanical system, are inevitable with the running of rotating machinery.To verify the present personalized diagnosis methodology in the real-world applications, the independent realization of artificial white Gaussian noise (AWGN) is added in the simulation using the command awgn in MATLAB (version 2010, MathWorks, Natick, MA, USA, 2010) as: where Data Noise−free and Data Noise are the noise-free data and noisy data and SNR is the signal-to-noise ratio (dB).
Figure 7 shows the results of adding artificial white Gaussian noise using Equation ( 5) under SNR = 97 dB.
In the present, the signals with the four fault types, such as unbalance, misalignment, rub-impact and the combination of rub-impact and unbalance faults, shown in Figure 7a-d, respectively, are decomposed with the WPT using Daubechies wavelet (Db18) into three layers.Therefore, eight frequency bands are obtained, and the corresponding five time domain feature parameters are calculated using the formulas in Table 1; the values are listed in Table 2.As shown in Table 2, for each fault type of the shaft, 5 × 8 = 40 features are forming a sub-feature vector.In order to obtain satisfactory classification results, we use 20 sub-data points (8000 data points in total) of the simulation signal, and hence, the length of the feature vector is 800.Therefore, the data used for training samples (including four types of faults) are an 800 × 4 matrix.In the simulation, the 'testing sample' is numerically represented by an 800 × 1 vector, which is the same simulation signal adding another SNR = 87 from the training samples shown in Figure 8a-d.To distinguish the four faults numerically, unbalance, misalignment, rub-impact and rub-impact based on unbalance faults are labeled from 1 to 4, respectively.In the process of classification, the adjustable parameters of SVM are selected as p1 = 1.05 and C = 10 [36].We use the SVM toolkit trained by Franc et al. [40], and four loops are necessary to obtain the four damage parameters.The final classification result is listed in Table 3; from the table, we can see that the accuracy rates of unbalance, misalignment and rub-impact are 93%, 95% and 89%.For the combination of rub-impact and unbalance fault, the result attains 91%, which verifies the proposed personalized diagnosis methodology.It is noted that the accuracy rates are given automatically by the SVM toolkit.More details about the multi-classification using the SVM toolkit for fault detection can be seen in [41].As shown in Table 2, for each fault type of the shaft, 5 × 8 = 40 features are forming a sub-feature vector.In order to obtain satisfactory classification results, we use 20 sub-data points (8000 data points in total) of the simulation signal, and hence, the length of the feature vector is 800.Therefore, the data used for training samples (including four types of faults) are an 800 × 4 matrix.In the simulation, the 'testing sample' is numerically represented by an 800 × 1 vector, which is the same simulation signal adding another SNR = 87 from the training samples shown in Figure 8a-d.To distinguish the four faults numerically, unbalance, misalignment, rub-impact and rub-impact based on unbalance faults are labeled from 1 to 4, respectively.In the process of classification, the adjustable parameters of SVM are selected as p 1 = 1.05 and C = 10 [36].We use the SVM toolkit trained by Franc et al. [40], and four loops are necessary to obtain the four damage parameters.The final classification result is listed in Table 3; from the table, we can see that the accuracy rates of unbalance, misalignment and rub-impact are 93%, 95% and 89%.For the combination of rub-impact and unbalance fault, the result attains 91%, which verifies the proposed personalized diagnosis methodology.It is noted that the accuracy rates are given automatically by the SVM toolkit.More details about the multi-classification using the SVM toolkit for fault detection can be seen in [41].

Experimental Investigations
To investigate the effectiveness of the proposed personalized fault diagnosis methodology of the shaft, an experimental investigation is conducted in this section.The experimental setup of a rotor-bearing system under normal load conditions is shown in Figure 9.It includes displacement sensor, spindle, rotor controller, data acquisition instrument, data processing software, etc.For the experimental rotor-bearing system, the length of the shaft and the disk mass are the same as the geometry model as shown in Figure 4.In the experimental investigation, the rotation speed is kept constant at 1800 rpm, as well as in numerical simulation; the sampling frequency fs = 6000 Hz, and the signal length is 8000 points.

Experimental Investigations
To investigate the effectiveness of the proposed personalized fault diagnosis methodology of the shaft, an experimental investigation is conducted in this section.The experimental setup of a rotor-bearing system under normal load conditions is shown in Figure 9.It includes displacement sensor, spindle, rotor controller, data acquisition instrument, data processing software, etc.For the experimental rotor-bearing system, the length of the shaft and the disk mass are the same as the geometry model as shown in Figure 4.In the experimental investigation, the rotation speed is kept constant at 1800 rpm, as well as in numerical simulation; the sampling frequency f s = 6000 Hz, and the signal length is 8000 points.

Confirm the Stiffness and Damping Coefficients
In real-world practice, a large difference between the measured vibration signal and the computed one (by the finite element method) is inevitable.Therefore, in order to check on the closeness between simulation signals and experimental vibration signals and to ensure the final stiffness and damping coefficients, the finite element model updating technique called TIC [37] is used.The value of TIC serves as the objective function as: x n y n x y x n y n (6) where ρ( , ) x y is the coefficient between the measured signal and the average of the simulated signal, x(n) is the measured signal and ( ) m y n is the mean of the simulation signal.As a figure-of-merit, the value of TIC ρ( , ) x y has a numerical value between zero and one.If x y is near zero, this corresponds to accurate predictions by the model.Generally, in the engineering field, ρ( , ) x y < 0.4 will lead to a satisfactory result [37].
The classic linear model with four spring and four damping coefficients (shown in Figure 4b) is employed to model bearings in the present work.In the present work, the experimental setup is the geometry symmetry, and the two bearing can be treated as the same model.Nine combinations are carefully selected and listed in Table 4.For the present experimental setup, the test signal is simple, and ρ( , ) x y < 0.4 can be easily satisfied.Therefore, we just present nine combinations to fulfill the alternative optimization.

Confirm the Stiffness and Damping Coefficients
In real-world practice, a large difference between the measured vibration signal and the computed one (by the finite element method) is inevitable.Therefore, in order to check on the closeness between simulation signals and experimental vibration signals and to ensure the final stiffness and damping coefficients, the finite element model updating technique called TIC [37] is used.The value of TIC serves as the objective function as: where ρ(x, y) is the coefficient between the measured signal and the average of the simulated signal, x(n) is the measured signal and y m (n) is the mean of the simulation signal.
As a figure-of-merit, the value of TIC ρ(x, y) has a numerical value between zero and one.If ρ(x, y) is near zero, this corresponds to accurate predictions by the model.Generally, in the engineering field, ρ(x, y) < 0.4 will lead to a satisfactory result [37].
The classic linear model with four spring and four damping coefficients (shown in Figure 4b) is employed to model bearings in the present work.In the present work, the experimental setup is the geometry symmetry, and the two bearing can be treated as the same model.Nine combinations are carefully selected and listed in Table 4.For the present experimental setup, the test signal is simple, and ρ(x, y) < 0.4 can be easily satisfied.Therefore, we just present nine combinations to fulfill the alternative optimization.

Obtain the Measured Signals and Classification Results
The four different faulty signals are collected, including shaft unbalance, misalignment, rub-impact and the combination of rub-impact and unbalance.Figure 10 shows the four faults of the shaft in the experimental operation.In order to conduct the unbalance fault, we add a bob-weight (2 g) mounted on the disk, as shown in Figure 10a.Figure 10b shows the metal washer (which change the height and angle of shaft) added to the bottle of the motor to simulate the misalignment fault.We operate the rub-impact fault of the shaft using a screw contacting the shaft surface, as shown in Figure 10c.Figure 10d shows the combination of rub-impact and unbalance.The vibration signals of shaft unbalance, misalignment, rub-impact and the combination of rub-impact and unbalance are shown in Figure 11a-d   Table 6 presents the final classification results, i.e., the classification accuracy ratios for the misalignment and the unbalance are 82% and 87%, respectively, whereas, for the rub-impact and the combination of rub-impact and unbalance, the classification accuracy ratios are 73% and 79%.The final results show that the proposed personalized fault diagnosis methodology might be used in practical applications.However, further works, such as the high performance finite element model updating techniques, the contact surface stiffness and damping parameters, the two-and three-dimensional geometry models of complex mechanical systems, etc., should be investigated in depth to improve the classification accuracy ratios and extend the application of the present methodology.The 8000 data points collected for the different faults are divide into 20 segments, respectively, and each corresponding segment contains 400 data points.WPT using Daubechies wavelet (Db18) is employed to decompose each segment into three layers, and the data processing procedures are similar to those used in Section 3.2.Table 5 shows the typical five time domain feature parameters for the eight frequency bands in the first segment (400 data points) gained from four faults (shown in Figure 11).Table 6 presents the final classification results, i.e., the classification accuracy ratios for the misalignment and the unbalance are 82% and 87%, respectively, whereas, for the rub-impact and the combination of rub-impact and unbalance, the classification accuracy ratios are 73% and 79%.The final results show that the proposed personalized fault diagnosis methodology might be used in practical applications.However, further works, such as the high performance finite element model updating techniques, the contact surface stiffness and damping parameters, the two-and three-dimensional geometry models of complex mechanical systems, etc., should be investigated in depth to improve the classification accuracy ratios and extend the application of the present methodology.

Conclusions
Intelligent diagnosis methods play a key role to prevent and detect faults in mechanical systems, except for the lack of faulty training samples.Motivated by molecular dynamics simulation-based personalized medicine, we develop the concept of the numerical simulation-based personalized diagnosis methodology and, more specifically, a personalized fault diagnosis method using FEM, WPT and SVM to detect faults in a shaft.Numerical and experimental investigations are performed, and the fault diagnosis of the shafts is given.In the simulation, the classification accuracy rates of unbalance, misalignment, rub-impact and the combination of rub-impact and unbalance are 93%, 95%, 89% and 91%, respectively.However, in the experimental investigations, this decreased to 82%, 87%, 73% and 79%.The results verified that the proposed method has the ability to distinguish the fault types.The numerical simulation-based personalized diagnosis methodology has many advantages, such as the activation of intelligent diagnosis methods by providing the complete faulty samples; the resolution of diagnosis reliability for a single machine stems from assembly variation with dynamic behaviors, influences from the operation environment, etc.
However, through the presented methodology, which is more likely to be expanded into more complex mechanical systems, such as gears, bearings, transmission systems, etc., the development of a standard FE model updating technique is the most important to obtain the agreeable numerical simulation faulty samples under different work conditions for the same type of mechanical systems.Furthermore, the computational cost is small for the present one-dimensional (1D) FE model employed.For complex mechanical systems, the calculation process is time consuming, but the proposed methodology only requires pre-computing the training samples, which serve as the inputs into the intelligent diagnostic models, such as SVM, GAs, NNs, etc.

Figure 1 .
Figure 1.The fundamental procedure of the personalized diagnosis methodology.Figure 1.The fundamental procedure of the personalized diagnosis methodology.

Figure 1 .
Figure 1.The fundamental procedure of the personalized diagnosis methodology.Figure 1.The fundamental procedure of the personalized diagnosis methodology.

2 ClassifiedFigure 3 .
Figure 3.The flowchart of the proposed fault diagnosis method.Solid and dotted arrows mean the training and testing paths to realize the method.
Appl.Sci.2016, 6, 414 7 of 19percent is commonly accepted), the FE model of shafts will well represent the real-world ones.Then, the faulty models will be inserted into the FE model to simulate the dynamic response of the shaft with faults.

Figure 6 .
Figure 6.The exciting force in each direction of the rub-impact: (a) the exciting force in x-direction; (b) the exciting force in y-direction.

Figure 9 .
Figure 9.The experimental setup of a rotor-bearing system under normal load conditions.

Figure 9 .
Figure 9.The experimental setup of a rotor-bearing system under normal load conditions.

Table 4 . 1 K
The nine different combinations of stiffness and damping coefficients and the corresponding ρ(x, y) 1 .(N/m), C (N•m/s).The first 400 data points from the experimental and the simulation signals are collected to validate the nine combinations of stiffness and damping coefficients.The values of TIC ρ(x, y) are shown in the last column of Table 4, and the least value ρ(x, y) = 0.1821 indicates the eight coefficients in the combination number; D is the relatively best choice to update the finite element model.The four stiffness coefficients and the four damping coefficients of the rotor-bearing system in the x+, x− direction and the y+, y− direction are: K x+ = 16.3 × 10 7 N/m, C x+ = 1500 N•S/m, K x− = 26.5 × 10 7 N/m, C x− = 1100 N•S/m and K y+ = 15.3 × 10 7 N/m, C y+ = 1600 N•S/m, K y− = 25.5 × 10 7 N/m and C y− = 1000 N•S/m, respectively.

Figure 10 .
Figure 10.Four faults of a shaft in the experimental operation.(a) Unbalance; (b) misalignment; (c) rub-impact; (d) the combination of rub-impact and unbalance.

Figure 10 .
Figure 10.Four faults of a shaft in the experimental operation.(a) Unbalance; (b) misalignment; (c) rub-impact; (d) the combination of rub-impact and unbalance.

Figure 11 .
Figure 11.The measured signals with different kinds of faults.(a) The measured signals with unbalance fault; (b) the measured signals with misalignment fault; (c) the measured signals with rubimpact fault; (d) the measured signals with the combination of rub-impact and unbalance fault.

Figure 11 .
Figure 11.The measured signals with different kinds of faults.(a) The measured signals with unbalance fault; (b) the measured signals with misalignment fault; (c) the measured signals with rub-impact fault; (d) the measured signals with the combination of rub-impact and unbalance fault.

Table 1 .
Five time domain feature parameters 1 .

Table 1 .
Five time domain feature parameters 1 .

Table 2 .
The calculation results in each frequency band of each fault.

Table 3 .
Fault recognition results with simulation.

Table 3 .
Fault recognition results with simulation.

Combination of Rub-Impact and Unbalance The Calculated Results in 8 Frequency Bands of Compound Faults
, respectively.

Table 6 .
The different types of faults' recognition results.

Table 6 .
The different types of faults' recognition results.