Damage Degree Detection of Cracks in a Locomotive Gear Transmission System

-e gear transmission system is a vital component of the locomotive bogie. During locomotive operation, the damage to the gear transmission system spreads rapidly and affects the locomotive’s operational safety. In this paper, a method is proposed to detect the degree of tooth root crack damage. First, a dynamic locomotive model with a gear transmission is built, and the vertical acceleration of the locomotive subsystem (car body, bogie frame, wheelset, and motor) vibrations is obtained under various degrees of tooth root crack damage on the gear transmission system. By comparing the characteristics of those signals, the subsystem that is more sensitive to the effect of the tooth root crack is found. -e characteristic parameters of the sensitive subsystem are calculated, and a multidimensional characteristic parameter matrix is established. -e multidimensional characteristic parameter matrix is optimized and reduced by principal component analysis (PCA). Using the Grey relational analysis method, the degree of tooth root crack damage is detected. -e proposed method demonstrates the ability to recognize the degree of tooth root crack damage.


Introduction
e gear transmission system transmits the torque of the motor and is a key element of the locomotive system.During locomotive operation, the gear transmission system is subjected to the mesh force of the gear pair and the wheelrail impact force.When the gear transmission system goes out of order, the damage spreads rapidly and affects the locomotive's operational safety [1].erefore, the detection method of the gear transmission system is the key to ensure locomotive safety.
In recent years, scholars worldwide have carried out much research on gear dynamics.Wu et al. [2] established a six-degree-of-freedom dynamic model and analyzed gear dynamic behavior with a change in the gear root crack.
rough the dynamic analysis of the gear transmission system, Endo et al. [3] noted that the change in the meshing stiffness and the geometric error of the tooth shape are important causes of tooth root crack failure.e dynamic behavior of gear wear failure is analyzed by Ding and Kahraman [4].Cui et al. [5] established the dynamic model of the gear rotor coupling system and studied the influence of the rotational speed on the dynamic response.In 2011, Ma and Chen [6] studied the nonlinear dynamics of a gear system with a tooth root crack.To research the spalling mechanism, Ma et al. established a dynamic model of a gear pair with local spalling [7].Bahk and Parker developed an analytical model used to study the influence of tooth profile modification on the planetary gear dynamics [8].In 2015, Li and Peng developed a multi-degreeof-freedom nonlinear gear model that considers the meshing force and tribological force of the gear pair [9].To analyze the influence of helical gear wear on the dynamic response of the gear transmission, Brethee et al. developed a comprehensive dynamic model that consists of 18 degrees of freedom [10].
At present, scholars have studied the fault diagnosis of the gear transmission system and put forward a series of effective diagnosis methods.To monitor the mechanical condition of a railway traction system, Kia et al. used timefrequency signal processing techniques to address the machine current signal [1].Shao et al. introduced the support vector machine (SVM) method for damage detection of the mechanical transmission system [11].In 2011 [12], Henao et al. put forward an electromagnetic-torque estimation method to detect tooth damage in a high-speed railway traction system.Cao et al. proposed a method for detecting gear damage based on empirical mode decomposition (EMD) and high-order cumulant (HOC) analysis [13].
is method is used to identify the damage to the gear transmission system under various operating conditions.In 2013, [14], Yang and Cheng proposed a damage classification method based on local mean decomposition and principle component analysis methods.In 2015, [15], Kia et al. presented a noninvasive technique using space vector analysis of the stator current to detect gear tooth surface damage faults.In 2016, [16], Heydarzadeh et al. used the discrete wavelet transform and deep neural networks to realize the automatic detection of gear faults.In 2017, [17], Chai et al. introduced a method using resonance-based sparse signal decomposition combined with envelope spectrum analysis to extract the fault features of motor current.Peng et al. realized the fault detection of wind turbine drivetrain gearboxes using the recurrent neural network [18].Ding et al. proposed a new method to detect the fault of planetary gearboxes based on local mean decomposition and permutation entropy [19].
However, the main object of those methods is to detect the existence of damage.
e detection of the degree of damage of rotating machinery is the focus of this field.In this paper, the locomotive subsystem which can be used to reflect the gear tooth root crack is found.e degree of tooth root crack damage is detected based on the PCA and Grey relational analysis methods.

Locomotive-Track Vertical Coupled
Dynamics Model with a Gear e locomotive is modeled as a multirigid body system running on the track.It consists of one car body, two bogie frames, four wheelsets, four motors, and a gear transmission system.ese components are connected by spring-damper elements.e car body connects with the bogie frame by a secondary suspension.e bogie frame is supported on the wheelset via the primary suspension.One end of the motor is fixed on the bogie frame while the other end of the motor is fixed on the wheelset.Each component of this mathematical model has two vibrational degrees of freedom, which are torsional and vertical.
e track subsystem contains the rail, the sleepers, the ballasts, and the subgrade.ese components are connected by spring-damper elements.e rail is regarded as a continuous Bernoulli-Euler beam.It is mounted on the sleepers by the rail pads.Zhai et al. [20] provide a method of modeling railway ballasts.In the new model, the shear stiffness and shear damping are introduced between the adjacent ballasts.Each component of this mathematical model has one vibrational degree of freedom, which is vertical.
e notation used in this study is shown in Tables 1-3.e schematic of the total dynamic system is depicted in Figure 1(a).e schematic of the gear transmission is depicted in Figure 1(b).Figure 1(c) is a detailed illustration of the gear transmission system.e gear transmission consists of motor, pinion, and gear.
e dashed circle represents the base circle of the gear.e rotor of the motor is linked with the pinion by the torsional spring-damper element.e pinion drives the gear to rotate through the engaging force.e gear is mounted on the wheelset.In this system, only the rotational motion of the components is considered.
As the locomotive is modeled as a multirigid body, the differential equations of the locomotive with the gear transmission can be established by D'Alembert's principle.
Vertical vibration of car body: Pitch motion of car body: Vertical vibration of bogie frame: Pitch motion of bogie frame: Vertical vibration of motor: Pitch motion of motor: Vertical vibration of wheelset: Rotational motion of motor rotor: Rotational motion of the pinion: 2 Shock and Vibration Rotational motion of the gear: Rotational motion of the wheelset: e forces and moments are given by the following formulas: Force between the car body and bogie: Force between the bogie and wheelset: Force between the bogie and motor: Force between the motor and wheelset: In equations ( 12)-( 15), the symbol "round" refers to the function for rounding.l c denotes the half distance between the centers of the two bogie frames.l t is the half distance between the centers of two wheelsets.l s1 represents the longitudinal distance from the motor center of gravity to the center of the frame.l s2 represents the longitudinal distance from the pinion shaft to the motor center of gravity.l s3 represents the longitudinal distance from the pinion shaft to the gear shaft.
e symbols P, F m , and F creep , respectively, represent the wheel-rail interaction force, the gear transmission mesh force, and the wheel-rail interaction creep force.e gear transmission mesh force, P and F creep are shown as follows: where the symbols δ represents the transmission error.e dots above the variables δ indicate a first-order differential operation.e transmission error δ can be calculated by the following formula: where R p and R g denote the radius of the pinion and gear base circles.Where θ p and θ g denote the angular displacements of the pinion and gear.e meaning of symbol Z moi and Z wi (i � 1 − 4) can be found in Table 1.In this paper, the manufacturing and the assembly error are ignored.e symbol G in equation ( 16) represents the wheel-rail contact constant.e symbol δZ(t) represents the wheel-rail elastic compressive deformation.Zhai and Cai [21] notes that the wheel-rail elastic compressive deformation includes the static compressive deformation, and it is calculated as follows: where the symbols Z w and Z r are the vertical distances of the wheel and the rail on the contact point, respectively.Torque that actuates the rotor Shock and Vibration 3

Bogie frame
Wheel Motor  Shock and Vibration e symbol Q in equation ( 17) represents the normal load of the wheel-rail contact.e symbol μ in equation (17) represents the adhesion coefficient of the wheel-rail interface.Ishikawa gives the formula as follows: where a, b, c, and d are calculated constants and v slip is the difference between the wheel rotational velocity and the locomotive velocity.

Gear Mesh Stiffness Model with Tooth Crack Damage.
During gear transmission system operation, the tooth root is subjected to repeated bending stress.After some time, cracks will appear at the tooth root site.When the teeth with tooth root cracks enter the engagement, the mesh stiffness changes.
In 1987, Yang and Lin [22] used the potential energy method to calculate the gear mesh stiffness.ey assumed that the gear meshing energy consists of three parts: U h (hertz energy), U b (bending potential energy), and U a (axial compressive energy).Tian [23] improved the mesh stiffness model in 2004.In this method, the U s (shear energy) was considered.In 2013, Chen and Shao [24] propose a much accurate method for calculating gear mesh stiffness.In this method, they consider the effect of the tooth root fillet-foundation deflection and the range of the cantilever beam is between the base circle to the dedendum circle.In order to calculate the mesh stiffness of tooth profile modifications and profile shifted gear, Ma et al. [25] revised the model put forward by Chen.In 2016, Ma and Zeng [26] put forward a method to calculate the mesh stiffness of gear pairs with tip relief.
During gear transmission system operation, the tooth root is subjected to repeated bending stress.After some time, cracks will appear at the tooth root site.When the teeth with tooth root cracks enter the engagement, the mesh stiffness changes.References [27,28] consider that a crack develops at the root of a single tooth.In most cases, the crack is a straight line.Based on the above study, the crack propagation path is assumed to be a smooth and continuous line.Reference [29] indicates that when the gear has a tooth root crack, parameters I x and A x will change.In 2014, Ma et al. [30,31] put forward a model to calculate the mesh stiffness of gear with tooth root crack considering the accurate transition curve of tooth root.In 2013, Mohammed et al. [32] assumed that the crack of the gear tooth is parabolic, and the accuracy of gear stiffness calculation is improved.In order to calculate the mesh stiffness of gear with tooth crack in planetary gearbox, Liang developed an improved beam model to analyze the effect of the crack evolution on gear mesh stiffness.
In this paper, the mesh stiffness of the healthy gear can be calculated by equation (22), where k h , k b , k s , k a represent the Hertzian, bending, shear, and axial compressive mesh stiffness.ose indicators calculated by the formula equation (23), where F represents the acting force at the contact point.F can be divided into F a and F b , where F a is the radial force and F b is the tangential force.e schematic of gear force is shown in Figure 2. G and E denote the shear modulus and Young's modulus, respectively.I x and A x are the area moment of inertia and section area, respectively, where the distance from the tooth root is x. e total mesh gear potential energy can be calculated by equation (24).Chaari et al. [33] noted that the deformation of the base circle has an effect on the gear mesh stiffness.e mesh stiffness can be calculated by equation (25).e coefficients L * , M * , P * and Q * can be calculated by equation (26), and X * i represents these coefficients L * , M * , P * , and Q * .e coefficients A i , B i , C i , and D i can be obtained by Table 4. θ f is an angle that refers to half the base circle thickness.r f and r int represent the basic circle radius and gear bore, respectively.h fi is calculated by equation (27).
Shock and Vibration References [27,28] consider that a crack develops at the root of a single tooth.In most cases, the crack is a straight line.Based on the above study, the crack propagation path is assumed to be a smooth and continuous line.Reference [29] indicates that when the gear has a tooth root crack, parameters I x and A x will change.ey can be calculated by equations (28) and (29), respectively.e parameter h x is the vertical distance from the point on the curve to the point of the tooth's centerline, as shown in Figure 2. e bending potential energy and shear potential energy can be calculated by equations ( 22), (28), and (29).We can then finally obtain the mesh stiffness of a gear with a tooth crack.

Detection of the Degree of Locomotive
Gearbox Gear Crack Damage

Feature Selection Based on PCA (Principal Component Analysis)
. Principal component analysis (PCA) is a statistical analysis method that maps multiple features into a few comprehensive features.To systematically analyze a practical problem, indicators that are extracted using a time-frequency domain analysis need to be considered.However, each indicator reflects some information about the problem, and there is a certain correlation between the indicators.Principal component analysis takes effective data dimensionality reduction technology to extract several comprehensive and irrelevant features that reflect the comprehensive test information as much as possible [34,35].e principal component analysis process is shown as follows: (1) e standardization of data X � X 1 X 2 • • • X P  T is a p-dimensional random vector, and T is the sample data, where i � 1, 2, . . ., n. Z is the standardization matrix that is calculated by equation (30): where (2) e calculation of the correlation coefficient matrix R of the standardization matrix Z is done in equation (32).
(3) Calculate the principal component U. e characteristic root λ of the correlation coefficient matrix R can be obtained by equation (33).e eigenvector of the characteristic root λ i is b i .I is the characteristic vector matrix.
e contributions ratio of the ith principal component can be calculated by equation (33).
e cumulative contribution rate of the former m components can be calculated by equation (35).e number of the principal components can be determined by the value of tt.Finally, we can obtain the principal component U by equation (36):

Grey Relational Analysis.
e Grey relational analysis is a method of quantitative description of the development and change of a system [36].It determines the state of the system by calculating the correlation degree between the standard mode and the unknown mode.e process of detecting the degree of damage using the Grey relational analysis is shown as follows: (1) Structure standard damage mode matrix.In this paper, the standard mode matrix is established using the simulation data of the locomotive car Table 4: Values of the coefficients of equation (25).
where k is the number of damage states and n is the number of indicators of each standard state.
In this paper, we choose a series of statistical indicators (such as RMS, kurtosis, S α , S r , and FM4) to analysis the degree of gearbox gear crack damage.RMS, kurtosis, S α , S r , and FM4 can be calculated by where x(n) is the amplitude of the signal and N is the number of data.x is the mean value of the data and d is the difference signal of x(n).d is the mean value of d. e detailed flowchart for correlation degree is shown in Figure 3.
(2) Calculation of the correlation coefficient.Assume X i (j) � (X 1 (j), X 2 (j), . . ., X n (j)) T is the characteristic vector calculated from the data to be checked.X i (j) is the ith characteristic parameter.e correlation coefficient can be calculated by where ρ is the resolution coefficient.In this paper, the value of ρ is 0.5.(3) Detection of damage degree.
e vector of the correlation degree c is shown as follows: e damage mode corresponding to the maximum element of c is the state of the data to be checked.

Detection of Degree of Locomotive Gearbox Gear Damage.
e process of detecting the degree of gearbox gear damage is shown as follows: (1) Decompose the fault range that needs to be detected into m parts.(2) Calculate the n indicators of the vertical acceleration of car body vibration under each part that was decomposed in Section 1. ese indicators make up the standard damage mode matrix.e matrix and the calculated indicators that are to be verified make up a matrix with the size of (m + 1) × n. (3) Use the method of PCA to reduce the dimension of matrix.(4) Solve the correlation degree vector between the signal to be verified and the standard mode using the Grey relational analysis method.(5) Determine the maximum value of the correlation degree vector corresponding to the degree of gear damage.e damage degree detection chart is shown in Figure 4.

Simulation and Result Discussion
With the development of the vertical dynamics model of the locomotive-track coupled with the gear transmission, the dynamic responses of the locomotive subsystem can be obtained.
We adopt the locomotive model reported earlier.e running speed of the locomotive is 100 km/h.e radius of the wheelset is 0.625 m.During locomotive operation, the mesh force of the gear transmission system is the primary excitation that causes locomotive vibrations.e notation for the main parameters of the gearbox is in Table 5.Using the above parameters, we can obtain the gear mesh frequency (f m ) and the rotation frequency of the pinion (f p ): 848.8 Hz and 36.9Hz, respectively.e total mesh stiffness of gear pair can be calculated by equations ( 22)- (29).When there is a crack in the gear tooth root, the mesh stiffness of the gear pair will change.e mesh stiffness of gear pair under different crack sizes is calculated and shown in Figure 5. e angle θ between the crack direction and the tooth center is 45 °.As shown, when there is a crack on the pinion, the mesh stiffness of the gear pair is decreased.

Acceleration Comparison of the Car Body, Frame, and Wheelset.
e dynamic model of the locomotive is Shock and Vibration a large-scale system that contains the car body, bogie, and wheelset.To determine the component with an acceleration that shows more sensitivity to a tooth root crack, we compare the acceleration data of the car body, bogie frame, and wheelset.Figure 6(a) is the time domain diagram of vertical vibration acceleration of car body when there is no damage on the gear.Figure 6(b) is the time domain diagram of vertical vibration acceleration of car body when the crack length of the pinion is 4 mm.By comparing these diagrams, we find that when there is a crack on the root of the pinion, the amplitude of the acceleration of the car body is bigger.
Figure 7(a) is the frequency spectrum diagram (in logarithm form) of vertical vibration acceleration of car body when there is no damage on the gear.Figure 7(b) is the frequency spectrum diagram of vertical vibration acceleration of car body when the crack length of the pinion is 4 mm.And the frequency range of these two diagrams is 0-10000 Hz.As we can see, the mesh frequency (f m � 848.8 Hz) and the harmonics of the gear pair is apparent.In order to better observe the low frequency part of Figure 7, we restrict the abscissa range of Figure 7 between 0 and 1200 Hz and obtain Figure 8.As we can see, when there is a crack on the pinion, the rotation frequency of the pinion (f p � 36.9Hz) and the harmonics of the rotation frequency of the pinion appear.
Figure 9(a) is the time domain diagram of vertical vibration acceleration of car body when there is no damage on the gear.Figure 9(b) is the time domain diagram of vertical vibration acceleration of car body when the crack length of the pinion is 4 mm.By comparing these diagrams, we find that when there is a crack on the root of the pinion, the amplitude of the acceleration of the car body is bigger.
Figure 10(a) is the frequency spectrum diagram (in logarithm form) of vertical vibration acceleration of car body when there is no damage on the gear.Figure 10(b) is the frequency spectrum diagram of vertical vibration acceleration of car body when the crack length of the pinion is 4 mm.And the frequency range of these two diagrams is 0-10000 Hz.As we can see, the mesh frequency (f m � 848.8 Hz) and the harmonics of the gear pair are apparent.In order to better observe the low frequency part of Figure 10, we restrict the abscissa range of Figure 10 between 0 and 1200 Hz and obtain Figure 11.As we can see, when there is a crack on the pinion the rotation frequency of the pinion (f p � 36.9Hz) and the harmonics of the rotation frequency of the pinion appear.Shock and Vibration frequency spectrum diagram of vertical vibration acceleration of car body when the crack length of the pinion is 4 mm.And the frequency range of these two diagrams is 0-10000 Hz.As we can see, the mesh frequency (f m 848.8 Hz) and the harmonics of the gear pair is apparent.In order to better observe the low frequency part of Figure 13, we restrict the abscissa range of Figure 13 between 0 and 1200 Hz and obtain Figure 14.As we can see, when there is a crack on the pinion, the rotation frequency of the pinion (f p 36.9 Hz) and the harmonics of the rotation frequency of the pinion appear.In addition, the di erence in the gear pair mesh frequency and the pinion rotation frequency is appear.By comparing the vertical acceleration of the car body, bogie, and wheelset vibrations, we nd that when there is a crack on the pinion, the amplitude of the acceleration of these vibrations are larger in the time history and the frequency spectrum.e rotation frequency of the pinion and the harmonics of the rotation frequency of the pinion appear.Due to the in uence of the spring-damper system, the vertical acceleration of the car body and frame vibrations attenuates much more than in the wheelset.
e vertical acceleration of the car body is more sensitive to the e ect of the tooth root crack.

Detection of Degree of the Gear Box Gear Crack Damage by the Vertical Acceleration of the Car Body Vibration.
rough the above analysis, we nd that the vertical acceleration of the car body vibration is more sensitive to the e ect of the tooth root crack.erefore, we choose the vertical acceleration of the car body vibration to detect the degree of gear box gear crack damage.In this paper, the crack length of the pinion from 0 to 5 mm is simulated and analyzed.e length of the crack grows in 0.5 mm increments.e gear pair mesh sti ness under di erent crack sizes is calculated and shown in Figure 3.As shown, with an increase in tooth root crack length, the gear meshing sti ness decreases.
e process of detecting the degree of gearbox gear damage is shown as follows: (1) Decompose the fault range that needs to be detected into m parts.In this paper, the crack length of the pinion is increased from 0 to 5 mm in 0.5 mm increments.erefore, the damage states that it is 11.(2) Calculate the 5 indicators (RMS, kurtosis, S α , S r , and FM4) of the vertical acceleration of car body vibration under each damage states.And those indicators of matrix are shown in Table 6.ese indicators make up the standard damage mode matrix.
e acceleration of the car body vibration with a crack length of 3.25 mm was veri ed. e 5 indicators of the signal to be veri ed are shown in Table 7. e standard damage mode matrix and the calculated indicators that are to be veri ed make up a matrix with the size of 12 × 5.
(3) Use the method of PCA to reduce the dimension of matrix.e principal component analysis method proposed in the Section 3.1 is applied to deal with the standard damage mode matrix.
(4) Use the Grey relational analysis method proposed in the Section 3.2 to solve the correlation degree vector between the signal to be veri ed and the standard   The normalization value of the correlation degree vector The damage degree of gear box gear mode.So, we can obtain the normalization curve of the correlation degree vector, as shown in Figure 12. e maximum correlation degree vector is 0.9765.e crack length corresponding to the maximum correlation degree is 3 mm.erefore, the crack length of the detected data is between 2.5 and 3.5 mm.It can be seen from Figure 12 that the correlation curve is consistent.
e value of the correlation degree becomes larger when the value of the standard damage mode approaches the data to be verified.

Conclusions
In this paper, a vertical dynamic model of the locomotivetrack coupled with the gear transmission is established.Using a comprehensive model, the dynamic response of the locomotive subsystem is obtained under various degrees of tooth root crack damage in the gear transmission system.By comparing the vertical acceleration for vibrations in an undamaged locomotive and in one with a gear tooth root crack, we find that when there is a gear tooth root crack on the pinion of the gear transmission system, the rotation frequency and the harmonics of the rotation frequency of the pinion appear.
e vertical acceleration of the car body vibration is more sensitive to the effect of the tooth root crack.e degree of tooth root crack damage is detected using PCA and the Grey relational analysis methods.
e track excitation is inevitable when a locomotive operating at track.However, in this paper, track excitation is not considered.
e dynamic response of a locomotive subsystem to track excitation will be investigated in the future.

Figure 1 :
Figure 1: Vertical dynamics model of locomotive-track coupled with gear transmission: (a) total dynamic system, (b) schematic of the gear transmission, and (c) detail illustration of the gear transmission system.

Figure 6 :Figure 7 :Figure 8 :
Figure 6: Time domain diagram of vertical vibration acceleration of car body: (a) health gear (b) the crack length of the pinion is 4 mm.

Figure 13 (Figure 9 :Figure 10 :Figure 11 :
Figure 13(a) is the frequency spectrum diagram (in logarithm form) of vertical vibration acceleration of car body when there is no damage on the gear.Figure 13(b) is the

Figure 12 :Figure 13 :
Figure 12: Time domain diagram of vertical vibration acceleration of wheelset: (a) health gear (b) the crack length of the pinion is 4 mm.

Figure 14 :
Figure 14: Frequency spectrum diagram of vertical vibration acceleration of wheelset (frequency range: 0-1200): (a) health gear (b) the crack length of the pinion is 4 mm.

Figure 15 :
Figure 15: Normalization curve of the correlation degree vector.

Table 1 :
Notation used in the locomotive subsystem.

Table 2 :
Notation used in the track subsystem.

Table 3 :
Notation used in the gear transmission system.

Table 5 :
Main parameters of the gearbox.

Table 7 :
Indicators of the vertical acceleration of car body vibration when crack length is 3.25.