Next Article in Journal
NG-GAN: A Robust Noise-Generation Generative Adversarial Network for Generating Old-Image Noise
Previous Article in Journal
Weigh-in-Motion System Based on an Improved Kalman and LSTM-Attention Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Root Strain Response Characteristics of Inner Ring of Planetary Gear Transmission System with Crack Fault

College of Mechanical Engineering, Xinjiang University, Urumqi 830047, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(1), 253; https://doi.org/10.3390/s23010253
Submission received: 22 November 2022 / Revised: 16 December 2022 / Accepted: 23 December 2022 / Published: 26 December 2022
(This article belongs to the Section Intelligent Sensors)

Abstract

:
This paper established the system dynamics model for two kinds of tooth cracks of different depths of the sun gear and inner gear ring to study the influence mechanism of crack failure on the tooth root strain of the planetary gear transmission system. Combined with the finite element model of the inner gear ring, the tooth root strain of the ring was solved. Experiments verified the correctness of the solution method. The root strain under the crack fault of the sun gear and the tooth crack fault of the inner gear ring is analyzed, and the following conclusions are drawn: Periodic fault impact occurs in the strain signal of the tooth root of the inner gear ring during the crack fault of the sun gear root. The fault can be extracted by the fast spectral kurtosis method (FSK), and the fault components are used to determine whether the sun gear cracks. The Lempel–Ziv index showed a tendency to increase gradually during the process of solar wheel crack deepening, which could be used as the damage index of crack depth. The results can provide a basis and reference for fault diagnosis.

1. Introduction

Planetary transmissions are widely used in aerospace, automotive, and heavy industrial applications such as helicopters, wind turbines, and heavy trucks due to their large transmission ratio and strong load-bearing capacity [1,2]. However, the planetary gear transmission system, complicated structure, and poor working conditions led to an increased equipment failure rate. The gear crack fault is an early failure, occupying a great proportion of the fault, and the internal tooth circle is an essential part of the planetary gear transmission system. Therefore, it is vital for fault diagnosis to study the influence of cracks on the root strain of the inner ring and to judge the degree of crack damage.
There have been many studies on the crack failure of the planetary gear system. In the fault simulation analysis, Chaari et al. [3] simulated the tooth pitting and cracks of the solar gear of the planetary gearbox. They analyzed the influence of the fault on the gear mesh stiffness. In addition, they compared the dynamic response of a healthy planetary gearbox with that of a planetary gearbox with eccentricity and contour errors by modeling [4]. Park et al. [5] used the finite element model to study the influence of defects on the carrier of the planetary gear set from the aspect of stress distribution. Yuksel and Kahraman [6] established a computational model to study the effect of surface wear on the dynamic behavior of planetary gearboxes. Zhang Jun et al. Han [7] analyzed the influence of ring gear crack fault on dynamic load distribution between external meshing gear pairs. Yang [8] considered the influence of tooth crack opening state on meshing stiffness and dynamic response of spur gear pair and proved that ring crack faults rarely affect the dynamic load distribution between external meshing gear pairs. Li Zhanwei et al. [9] Li Zhanwei et al. discussed the effects of the depth, length, and height of crack. The equivalent stress, contact pressure, and displacement of the tooth are also analyzed under different crack types using the FE method. Barszcz et al. [10] present the application of the spectral kurtosis technique for detecting a tooth crack in the planetary gear of a wind turbine. The authors propose a method based on spectral kurtosis, which yields good results. Their method was able to detect the existence of the tooth crack several weeks before the gear failure. Vicuna [11] proposed a phenomenological model to simulate the vibration that could be measured by a hypothetical sensor mounted outside the ring gear. Feng and Zuo [12,13] studied the spectral structure of planetary gearboxes’ vibration signals, provided a gear damage signal model, and demonstrated the signal model using experimental and industrial signals. Wu et al. [14] studied a planetary gearbox using a multi-body dynamics model. They observed that the dynamic response of the model depends on the interaction of many components inside the gearbox. Liang et al. [15] used a modified cantilever beam model to represent the external gear tooth and derive the analytical equations of the bending, shear, and axial compressive stiffness. A crack propagation model is developed, and the mesh stiffness reduction is quantified when a crack occurs in the sun or planet gear. Lei Yaguo et al. [16] studied the solution of time-varying meshing stiffness under the fault condition of a planetary wheel, established the fault dynamics model, and analyzed the spectrum characteristics under the fault condition. Xiao Zhengming et al. [17] adopted the improved energy method to establish a system dynamics model considering time-varying parameters, verified that the higher the model’s accuracy, the more pronounced the fault characteristics, and revealed the fault mechanism of early tooth root cracks. The modal sensitivity analysis is conducted using a three-dimensional dynamic model of a planetary geared rotor system for the number of planet gears, planet mistuning, mass of planet gears, gear mesh stiffness and planet gear speed [18]. Sanchez-Espiga et al. [19] proposes a numerical approach to the problem of the calculation of the load sharing in planetary transmissions by measuring the strains in the root of the sun gear teeth.
To sum up, the existing research mainly focused on various errors to reduce the influence of the whole model and the characteristic of the system of fault signal transmission. However, the study of the inner gear ring is insufficient. The planetary gear transmission system’s failure can often collect relevant fault signals on the internal gear ring. The ring gear is a critical component in planetary transmission. An easy-to-crack fault occurs. However, when analyzing fault response characteristics of vibration signals, it is easy to be affected by the transmission path. The fault components are weakened or disappear by modulation, which seriously affects the crack fault analysis. Therefore, different fault carriers are selected to reflect the fault state. As the signal carrier of fault, the strain signal is more direct than the vibration signal. The reason is that the change of time-varying meshing stiffness under crack failure will indirectly lead to the change of inner ring root strain. In contrast, the time-varying meshing stiffness caused by crack failure will directly change the strain signal of the inner ring root, which is more evident than the fault information transmitted by the vibration signal. Therefore, studying the root strain analysis of planetary transmission gear ring with the crack fault in fault diagnosis is essential.
In this paper, the dynamic model of the planetary gear transmission system was established, the failure excitation of the inner gear ring under the crack failure of the sun wheel was analyzed, the finite element model of the inner gear ring was established, the root strain of inner gear ring was solved, and the influence mechanism of root crack failure of a sun wheel and inner gear ring on the root strain of inner gear ring of planetary gear transmission was studied.

2. Construction of Planetary Transmission System Model

2.1. The Parameters of the Planetary Transmission Gearbox

This paper selects the planetary gear transmission system in the wind power generation testbed produced by SQI company for research. The straight-tooth planetary gear train is 2K-H type. The structural diagram is shown in the figure, including the sun gear S, the planet gear P, the planet carrier C, and the inner ring gear R. During the operation, the sun gear is used as the input and the planet carrier as the output components, and the four planet gears are uniformly distributed around the sun gear. The parameters of the above components are shown in Table 1. The ring gear rim thickness is 20 mm. In Figure 1, the planetary wheel is the input, and the sun wheel is the output. The gear ratio is 1.24. The planetary gear system plays an increasing role in the wind turbine gearbox.

2.2. Construction of System Dynamics Model

The complex coupling relationship between the planetary gear transmission system components makes modeling more challenging. Therefore, the gear transmission system is reasonably simplified. First of all, the appliance sometimes has variable stiffness of the spring links the parts, through which to reflect the meshing relationship between the parts. Secondly, the parts are regarded as rigid, and the support of the planet wheel and sun wheel and the fixing of the inner gear ring are simulated by springs with appropriate stiffness. Finally, a translation-coupled dynamic model was established without considering assembly errors, tooth clearance, and tooth machining errors.
The center of the sun wheel in the planetary gear system is taken as the center of the coordinate system, and the 100 teeth on the inner gear ring are numbered from 1 on the positive Y axis counterclockwise. Here, as shown in the figure, 1/4 of the whole inner gear ring is numbered 25 as an indication. The three degrees of freedom in transverse, longitudinal, and torsional directions are taken into account for each part of the transmission system: fs, ds, fr, ds, fc, and dc are the transverse and longitudinal micro displacements of the solar wheel, inner gear ring, and planetary frame respectively. fpi and dpi are the radials and tangential micro displacements of the planetary wheel. qs, qpi, qc, and qr are the torsional micro displacements of the sun wheel, planetary wheel, planetary frame, and inner gear ring. The 21 DOF displacement array of the system is as follows
X = f s , d s , θ s , f pi , d pi , θ pi , f c , d c , θ c , f r , d r , θ r
kr, kc, kpi and ks are the supporting stiffness of the inner gear ring, planet frame, planet wheel and sun wheel, respectively. krq, kcq, and ksq are the torsional stiffness of the inner gear ring, planetary frame, and solar wheel, respectively, where i = 1, 2, 3, 4.
As shown in Figure 2, when mesh stiffness of an internal spur gear pair with ring gear rim deformation is included, the relative displacement of each component in planetary transmission can be expressed as
δ spi = ( f p i f s ) sin φ spi + ( d s d pi ) cos φ spi + u s u pi δ rpi = ( f p i f r ) sin φ rpi + ( d r d p i ) cos φ rpi + d r θ r u pi δ cpx = f c f p i ( u c u p i ) sin φ rpi δ cpy = d c d p i ( u c u p i ) cos φ rpi δ spu = ( f c f p i ) sin φ i + ( d p i d c ) cos φ i + u c u pi δ x p i = f p i f c + u c sin φ i δ y p i = d p i d c u c cos φ i
where δspi is the projection of the relative displacement of the solar wheel and the planetary wheel along the meshing line; spi = αsi, where αs is the engagement angle of the external gear pair; δrpi: projection of relative displacement of planetary wheel and gear ring wheel along the direction of their meshing line; rpi = αr + i, αr is the engagement angle of the internal meshing gear pair; δcpx, δcpy, δcpu are the projection of relative displacement of the planetary shelf and planetary wheel along fc, dc, and θc directions, respectively; δxpi, δypi are the projection of the relative displacement of the planetary wheel and the planetary shelf along the fi and di directions.
The input/output torques Tin and Tout (input is positive, the output is negative) of the sun wheel and planetary frame, the mass of the sun wheel, planetary frame, planetary wheel, and inner gear ring wheel are ms, mc, mp and mr respectively. The moment of inertia is Ic, Ip and Ir, respectively. The differential equation of dynamics can be obtained by Newton’s Equation of motion using the lumped-mass method as follows:
m s f ¨ s i = 1 4 k spi δ spi sin φ spi + k sx f s = 0 m s d ¨ s + i = 1 4 k s p i δ s p i cos φ s p i + k sy d s = 0 ( I s / R 2 s ) θ ¨ s R s + i = 1 4 k sp i η sp i + k s θ θ s R s = T out / R s m c f ¨ c i = 1 4 k p i δ c p i + k cx f c = 0 m c d ¨ c + i = 1 4 k p i δ c p i + k cy d c = 0 ( I c / R 2 c ) θ ¨ c R c + i = 1 4 k p i δ cpu + k c θ θ c R c = T in / R c m p f ¨ c k sp i δ spi sin φ s p i k rp i δ rpi sin φ r p i k p i δ x p i = 0 m p d ¨ c k sp i δ spi cos φ s p i k rp i δ rpi cos φ r p i k p i δ y p i = 0 ( I p / R 2 c ) θ ¨ c R c k sp i δ spi k rp i δ rpi = 0 m r f ¨ r i = 1 4 k rp i η rp i sin ( α rp i + φ p i ) + k r f r = 0 m r d ¨ r i = 1 4 k rp i η rp i cos ( α rp i + φ p i ) + k r d r = 0 ( I r / R 2 r ) θ ¨ r R r + i = 1 4 k rp i η rp i + k r θ θ r R r = 0

2.3. Solution of Time-Varying Meshing Stiffness

ANSYS software was used to solve the meshing stiffness of the planetary gear transmission system by using the finite element method, as shown in the Figure 3. Taking the three-dimensional contact model of the planetary wheel and the solar wheel as an example, torque was applied to the driving wheel to fix the slave wheel. The meshes of the existing model are divided, and the meshes of the non-contact area are sparse. The contact surface is CONTA174, the target surface is discretized by the TARG170 element, and the main gear element is the SOLID185 element.
When the boundary conditions are applied, the total constraint is applied to the inner ring of the driving wheel, the lateral freedom of the inner ring of the driving wheel is constrained, and the torque is applied to the driving wheel. It should be noted that when torque is applied to the driving wheel, it should be applied to the position of the inner ring node of the driving wheel. The coordinate system needs to be converted into a cylindrical coordinate system, and the torque is decomposed into each node for application. In the load application process, the torque Ta to the inner ring node of the driving wheel is adopted.
T a = F n r b n r n
where Ta is the torque applied to a single node; n is the number of nodes. rn inner ring radius of the main driving wheel; rb is the radius of the base circle of the main driving wheel.
The inner ring deformation μ1 of the driving wheel was calculated by the deformation program, the driving wheel deformation Angle δ1 was obtained by the formula, and the equivalent deformation μ of the meshing line was finally obtained.
σ 1 = μ 1 r n μ = δ 1 r b
where δ1 is the main wheel deformation angle, rn is the main wheel inner ring radius, μ1 is the main wheel inner ring deformation, μ is the equivalent deformation at the meshing line.
Meshing stiffness is calculated as follows
k n = F n μ = T θ r b 2
The meshing period n of a tooth on the sun wheel is divided equally, and the meshing angle of the planet wheel is adjusted according to the meshing period after the equal division to obtain the meshing state of n sun wheels and the planet wheel. In addition, the meshing period n of each tooth on the planetary wheel was equally divided, the meshing angle of the planetary wheel was adjusted according to the meshing period, and the meshing state between the planetary wheel and the inner gear ring was obtained.
The solution result through the meshing stiffness calculation formula is shown in Figure 4 and Figure 5. As can be seen from the figure, the meshing process presents an alternation of two-tooth meshing to single-tooth meshing to two-tooth meshing. In the alternation process, meshing stiffness decreases sharply when meshing to a single tooth, and stiffness excitation is extremely significant.

2.4. Calculation of Gear Fault Stiffness

During the planetary gear transmission system operation, the gears are constantly engaged, and the gear teeth receive the continuous impact of the meshing excitation. Structural fatigue and cracks easily occur at the gear’s root, as shown in Figure 6, e is the crack depth, and gc is the crack angle. The meshing stiffness calculation method in Section 2.3 is used to establish a gear model with a tooth root crack here to solve the time-varying meshing stiffness under the fault of the tooth root crack. In the past, Refs. [9,10] used different gear crack modeling methods to calculate the stiffness of cracked gear. In this paper, the crack modeling method was simplified to facilitate the modeling and solving, and the crack model with uniform propagation depth along the tooth width was established. When cracks appear, the gear’s moment of inertia and effective cross-sectional area change, and the gear tooth model has greater flexibility. Therefore, the time-varying mesh stiffness decreases with the increase in crack depth. The meshing stiffness of the gear with the crack fault is shown in Figure 7.

3. Dynamic Load Analysis of Crack Fault

The solar wheel was used as the input, the speed was set at 600 r/min, and the load was 500 N·m. The meshing excitation was solved and analyzed. The Figure 8 shows that the dynamic excitation under normal working conditions shows periodic fluctuations. The dynamic excitation fluctuates up and down around the mean value of 2145 N, with a fluctuation period T of 0.00457 s and a fluctuation amplitude of 25 N. When there is a solar wheel crack fault, periodic fault impact appears in the dynamic excitation time domain, and the impact characteristics are instantaneous decrease and instantaneous increase. The instantaneous reduction of impact characteristics is caused by the instantaneous reduction of meshing stiffness when the fault tooth engages with a planetary wheel. Meanwhile, the dynamic excitation decreases with the decrease in stiffness, as shown in Figure 9b. At the same time, the dynamic excitation of the other three planetary wheels will increase instantaneously. This increase is relatively small compared to reducing the planetary wheels’ dynamic excitation with fault teeth engagement. According to the time domain analysis, the dynamic excitation still fluctuates around the mean value of 2145.61 N in the fault condition, and the time interval between the fault tooth engaging the same planetary wheel again is a complete fault cycle Tg = 0.128 s. At the same time, with the increase of crack depth, the mean value of dynamic excitation remains unchanged, and the fault impact amplitude increases.
Solar wheel single tooth crack fault, dynamic excitation fault frequency is solved as follows
f gs = ( n s n c ) / 60
where ns is the sun wheel speed, nc is the rotational speed of the planetary frame. In the time domain history of dynamic excitation, every time the solar wheel meshed with the same planetary wheel, the meshing excitation suddenly decreased once and increased three times. In the corresponding frequency domain analysis, three kinds of fault signals, fgs, 2fgs and 3fgs, appeared. In the operating condition of 600 r/min speed, fgs was 7.81 HZ, and frequency coupling appeared at 4fgs with a relatively small amplitude. Crack faults lead to fault components in the low-frequency band of the signal spectrum, and the amplitude of fault components increases with the increase of crack depth.

4. Solution of Response of Finite Element Model of Inner Ring

The tooth root crack of the sun gear and the tooth root crack of the inner ring was selected to analyze, aiming at the crack failure of the planetary gear transmission system. Analysis to solve the sun round the inner ring gear vibration response under the crack fault, the finite element model of trouble-free with the internal ring gear, the sun wheel crack fault conditions of extracting dynamic incentives, according to the tooth-bearing contact analysis (LTCA) uniform, applied to tooth surface contact, by ordinal position in tooth surface applied load equivalent gear meshing process. The dynamic equation of the inner ring under the action of meshing excitation is:
M v ¨ + c v ˙ + K v = G
where M is the mass matrix, C is the damping matrix, K is the stiffness matrix, V is the displacement vector of the node, G is the external load matrix.
For modal extraction of the finite element model, the first 50 modal matrices were selected for decoupling, and then the Newmark-β time integral method was adopted to solve Equation (8).
The main location of the solar wheel crack fault is the tooth of the solar wheel, so the dynamic meshing excitation under the fault condition is collected and applied to the finite element model of the inner tooth ring to solve the fault vibration response, as shown in Figure 10. The fault location of the inner gear ring crack is on the inner gear ring, and the dynamic meshing excitation under normal working conditions is collected and applied to the finite element model of the inner gear ring crack to solve the problem.

5. Experimental Verification

The wind turbine test bench produced by SQI company was used for the experiment, which was composed of a magnetic powder brake, two-stage fixed-axis reducer, wind power planetary accelerator, and drive motor. Strain signals were collected by a grating dynamic monitoring system, as shown in Figure 11.
In order to analyze the accuracy of the simulation signal, collect the time-domain history of the experimental strain signal and compare it with the simulation signal, observe the period of the time-domain history of the strain and the change of the strain value, compare the error between the simulation signal and the experimental signal, set the speed of the input shaft of the sun gear at 52.2 r/min and the load torque at 59 N·m for simulation and experiment, measure the strain data in the experiment respectively and simulate and solve the root strain signal. The time domain history comparison results of tooth root strain are shown in the Figure 12.
Comparing the simulation signal with the experimental signal, the time domain history of the signal can be divided into the meshing region and the non-meshing region. Meshes into the recess action of the meshing area of the planet wheel sensor location, rim tooth root, first squash and then stretch the process. As the chart shows, the simulation signal and experimental signal strain signals are in compression and tension of the trend, but as a result of the experiment are certain errors in the test method, the simulation results of the maximum tensile strain for 29 με. In the experiment, the maximum tensile strain is 26 με, and the error is 10.3%. In the non-meshing region, the experimental signal is compared with the simulation signal, and the simulation signal decreases to near zero with a slight fluctuation. In the experimental signal, the strain attenuation is relatively slow because the gear ring gradually recovers to the state before deformation when the gear teeth are gnawing out. The error is within tolerable limits. Therefore, the experimental results show that the simulation results are reliable.

6. Strain Analysis

6.1. Analysis of Failure Strain Results of Sun Gear Tooth

The dynamic strain method was used to solve the crack fault of the planetary gear transmission system. The crack fault of the solar wheel affects the strain signal results of the root of the inner ring of the planetary transmission system, as shown in the Figure 13.
Solving process, as shown by the sun around a crack of the gear fault, the fault gear and planetary gear, the meshing stiffness of volatility and fault gear meshing planetary wheel at the same time with the inner ring gear meshing, affected by the fault gear, planetary wheel and the ring gear meshing to stimulate a cyclical fluctuation, fluctuation cycle as shown in Figure 8. The input speed was set at 600 r/min and the load at 500 N·m. According to the finite element model solution process, the fault meshing excitation was applied to the finite element model of the inner gear ring. The solution results of strain extraction were as following Figure 14.
Under normal working conditions, because the planetary gear transmission system is a structure of four planetary wheels, and the four planetary wheels keep rotating with the planetary frame, the planetary wheels engage the tooth root strain collection position one by one. The strain shows periodic fluctuation, with the fluctuation cycle Tr being 0.1142 s, namely a quarter of the mesh cycle. The inner gear ring is supported by eight evenly distributed (Kr) fixed. As shown above, the tooth root stress extraction location selection between the adjacent two support and the time domain history can be divided into planetary gear extraction into the two supporting points to lead mesh deformation. The mesh over the strain extracts the location where the two supports are caused by the deformation of two kinds of vibration signal fluctuation.
By analyzing the time domain history of normal and fault working conditions, the dynamic excitation decreases when the sun gear tooth cracks occur, which further leads to the change of the tooth root strain. As shown in Figure 15, when the planetary wheel was about to engage two teeth near the extraction position of tooth root strain, it was not affected by the faulty teeth. Under the influence of compression deformation of the tooth ring, the compressive strain was not significantly different from that under normal working conditions. As the planet gear meshed through the tooth root strain extraction position, the cracked tooth meshed with the planet gear simultaneously, and the tensile strain weakened compared with the normal working condition. In dynamic excitation analysis, increased crack depth leads to different weakening degrees of dynamic excitation. Dynamic excitation weakens compared with normal working conditions, and the failure tensile strain also decreases, and the strain reduction degree is inversely proportional to the crack depth.
The root strain of the fault conditions a significant change in the mesh deformation area. This sets the planet wheel failure tooth contact with the sun at the same time and the inner gear tooth root circle strain mesh patch for the initial position until the planet wheel and inner gear ring patch and fault signal tooth meshing Tg for complete failure cycle at the same time. Again the sun wheel crack failure strain cycle for the Tg solution is as follows
T g = Z s + Z r , Z s n s N Z s
where Zs and Zr are the number of teeth of the solar wheel and planetary wheel, Ns is the sun wheel speed, and N is the number of planetary wheels.
As shown in Figure 15 (enlarging the deformation by 10,000 times), in the meshing deformation zone, when the planetary wheel was about to engage the sun wheel fault tooth, the tooth teeth on both sides of the strain collection point were affected by the compression deformation of the tooth ring. A large compressive strain appeared at the extraction position of the tooth root strain. As the planetary wheel meshed through the fault tooth, the meshing excitation caused the tensile deformation of the strain collection point, and the strain signal increased sharply. After the planetary wheel gradually gnaws the fault tooth, the strain signal fluctuates around 0. Analysis of tooth root stress wave curves under different meshing states because the internal gear ring structure has symmetry. The number 1~25 teeth are the research object, the tooth root position near the alveolar strain signal acquisition and analysis of the failure of meshes into gear. The meshing tooth and tooth broke down after the tooth of the three states of strain signal, setting the fault tooth in meshing, The planetary wheel meshes with 7~8 teeth on the inner ring. As shown in Figure 16, combined with trend analysis, the variations of dynamic incentive at will, while the link meshes into the tooth state failure within the planet wheel and near six teeth on the ring gear meshing, 5 and 6~8 tooth root position before subjected to tensile deformation, the rest of the tooth root position affected by compression deformation and dynamic incentive to keep normal at this time, so the working condition of two kinds of strain signals keep highly consistent. In fault gear meshing condition, 8 tooth root stress signal under the effect of fault appeared weak. It is caused by fault under dynamic excitation attenuation compression deformation is reduced, compared to the compressive strain of collected signals under normal working conditions of a less. The rest of the tooth’s position by compressive deformation is reduced, compressive strain decreases, and the tensile deformation decreases at the same time, as well as the decrease of strain. After the failure teeth were gnawing, the dynamic excitation remained normal, and the groove position of 9~10 teeth was subjected to compression deformation. The root strain had little difference from the normal condition, and the tension strain signals of other root positions did the same. The two adjacent supports (14 teeth~25 teeth) in the three states are in the vibration deformation zone. Due to the support’s fixed action, the tooth root’s strain signal always fluctuates weakly around 0, and the deformation is small.

6.2. Fault Strain Analysis of Inner Ring Teeth Crack

The specific flow of strain solution for cracks at the root of the inner ring is as following Figure 17.
Compared with the crack fault of the sun wheel, the meshing excitation under the crack fault of the tooth root of the inner ring is normal. Therefore, the system dynamics model is set at the normal working condition of 600 r/min speed and 500 N·m load. The meshing excitation is extracted, and the finite element model with cracks at a tooth root of the inner ring is established. Finally, the corresponding dynamic strain method calculates and analyses the results.
According to the symmetry of the gear transmission system, the 25 teeth marked on the inner ring were taken as the research object. The crack diagram of the inner ring tooth is shown in Figure 18, and the time domain history of the tooth groove strain was collected for comparative analysis, as shown in Figure 19.
Mesh deformation zone and local amplification, as shown in Figure 20. This can be divided into four a, b, c, d, and e (five) parts. The planetary wheel meshes into fault support but does not mesh into the tooth in the strain extraction location (Figure 20a). Figure 21 state 2, namely the planet wheel mesh to signal extraction location, signal acquisition position affected by gear ring tensile deformation, and tooth root part of the strain increase. The planetary wheel is about to engage in the strain acquisition position (Figure 20b), and the gear ring is convex and deformed at this position. The strain signal shows a sharp decrease under the action of compression deformation at the tooth groove position. The overstrain acquisition position of the planet wheel meshing (Figure 20c) and the tooth groove position was subjected to tension and compression deformation simultaneously. The strain increased suddenly, corresponding to state 3 in Figure 21. The planetary wheel is about to snap out the location of signal extraction (Figure 20d), which is subjected to tensile deformation, and the strain signal decreases. After the planetary wheel gnaws out the signal collection point (Figure 20e), the strain signal fluctuates near 0 with the deformation of the inner gear ring.
In the meshing deformation area, the crack fault on the inner ring changes the original inner ring structure. Dynamic excitation is applied to the tooth surface in the planetary gear meshing failure process, and larger deformation is concentrated at the crack position. As a result, the amplitude of compression deformation and tensile deformation at the position of the crack tooth groove decreases to a certain extent. When the fault tooth is located in the vibration deformation area, corresponding to state 1 in Figure 21, the strain phase at the tooth groove also appears to have certain weaknesses in normal working conditions.
In total, 25 teeth, including fault teeth on the tooth ring, were selected to be numbered, and a strain fluctuation curve was drawn. Typical states are shown in Figure 22. When the planetary gear meshes to state 2, the fault tooth does not engage, the position of the fault tooth is affected by the compression deformation of the tooth ring, and the tensile strain decreases. When meshing to state 3, dynamic excitation is applied to the fault tooth. The meshing stiffness decreases, and the tension strain and compression strain in the meshing deformation zone decreases. When the planetary gear runs to state 4, the corresponding position of the fault tooth is affected by the tensile action of the tooth ring, and the tensile strain decreases.
The planetary gear tooth crack fault is similar to the sun gear tooth crack fault, both changing the meshing excitation between the planetary gear and the inner gear ring. When the inner gear tooth crack fault occurs, the meshing excitation is not affected, so it is classified as the planetary gear and sun gear crack fault.

6.3. Crack Damage Discrimination

According to the strain signal of the inner ring under crack fault, crack damage discrimination was carried out on the planetary gear transmission system, and the specific process is as following Figure 23.
According to the fault mechanism analyzed above, when the sun wheel crack fault occurs, it causes a periodic fault impact on the strain signal of the tooth root on the tooth ring. The expression of the inner tooth root crack fault is that the amplitude of the strain signal caused by the fault tooth changes. Therefore, the root strain signal is first collected, and the fault type is distinguished according to the fault frequency component in the time domain signal. If there is no obvious fault frequency in the signal component, then the crack fault on the inner ring is judged by comparing the amplitude fluctuation of the normal signal.
After the strain signal collection, the signal components are analyzed. Kurtosis is sensitive to impact signals and can be used to evaluate the strength of vibration impact components such as bearings and gears. However, the kurtosis index is easily disturbed by noise and cannot reflect the change characteristics of the impact. In order to solve this problem, the spectral kurtosis (SK) index was proposed in recent years, which can reflect the intensity of the transient impact and indicate its frequency [20]. It is a method that combines the kurtosis index with spectrum analysis, which can detect fault impact signals very effectively. This method was first introduced into the field of fault diagnosis by Dwyer [21] as a theoretical statistical method. After that, Antoni. J [22] provides the formal theoretical basis for applying SK to fault diagnosis, thus filling the theoretical basis blank of SK in the field of fault diagnosis. A spectral kurtosis method based on the short-time Fourier transform (STFT) is proposed and applied to diagnose actual mechanical faults, as shown in Figure 24.
The fast spectral kurtosis diagram is analyzed. The maximum kurtosis value is 12.6, the corresponding center frequency is 815 Hz, and the bandwidth is 328 Hz. The processed fault signals are shown in the Figure 25 below.
According to the sun gear tooth crack fault cycle, Tg solved above is the corresponding fault frequency fg. In the signal envelope diagram of the sun gear tooth crack, the fault frequency, and its frequency doubling are obvious.
After processing the fault signal of the inner gear ring crack, the main components are four times the planetary frame rotation frequency and frequency doubling, and meshing frequency and no obvious fault components appear.
The fast spectrum kurtosis processing results of the fault strain signals of the sun wheel crack and the inner gear ring crack were compared, which was consistent with the characteristics of the above analysis. When the crack fault was determined, the fault component could be identified as the sun wheel or the inner gear ring according to the characteristics.
According to the signal characteristics of different parts, different methods are used to judge the degree of crack damage. If there is a solar wheel fault impact component in the crack fault signal, the complexity index sensitive to the periodic fault impact is selected to determine the damage degree.
As shown in Figure 26, 1, 2 and 3, respectively, refer to the three states of the normal, shallow crack fault, and deep crack fault. With the increase of crack depth, the index increases accordingly. Therefore, this index is used to measure the damage degree of the sun gear tooth under crack conditions.
If it is determined that there is no fault impact component in the crack fault signal, then the signal amplitude is further compared. Under the known crack fault state of the inner ring, the tension-compression deformation trend in the strain signal decreases. With the increase of crack depth, the tension-compression strain value decreases more. According to the characteristic of the fault strain of the inner ring crack, the root mean square value, which is sensitive to signal amplitude change, is selected to judge the degree of crack damage.
Figure 27 shows states 1, 2, and 3 represent the normal, shallow, and deep inner ring cracks, respectively. As the amplitude of the strain signals in the crack state changes, the changing trend of the signal amplitude is described by the root mean square value. As the crack depth increases, the root-mean-square value of strain signals in the inner ring crack root decreases continuously. This verifies the above analysis that the strain signals in the inner ring crack depth increase continuously.

7. Conclusions

The system dynamics model was established by taking the crack fault in the gear ring of planetary transmission as the research object. The strain signal was solved by combining the finite element model of the inner gear ring. The correctness of the model was verified through experiments, and the crack fault was added to analyze the strain signal of the solved tooth root, and the conclusion was reached:
(1)
When the sun wheel tooth cracks occur, the fault signal of tooth root strain will occur only if the tooth has meshed with the fault tooth. In contrast, when the tooth root cracks occur in the inner ring tooth, the strain at the tooth root position will fluctuate due to the influence of the crack, regardless of whether there is a planetary wheel meshing with the cracked tooth.
(2)
The strain signals under the sun wheel and inner ring crack fault can be distinguished by the fast spectrum kurtosis method (FSK).
(3)
According to the strain signal characteristics of the sun wheel’s crack fault and the inner ring’s crack fault, the damage degree of the crack fault of the sun wheel can be identified by the Lempel–Ziv algorithm. The damage degree of the crack fault of the inner ring can be identified by the root mean square index.
In this paper, the finite element analysis model considering the flexibility of the gear ring is established to study the crack damage response characteristics of planetary transmission systems from the vibration signal and strain signal of two fault carriers. However, there are still other aspects to be studied. (1) This paper analyzes the planetary gear transmission crack fault. However, the effects of planetary gear transmission crack propagation in different directions on system performance and the life changes of gear components after propagation can be further studied. (2) This paper analyzes the fault response characteristics of the crack and does not conduct quantitative research and discrimination on the fault damage degree. Subsequent research can be carried out for the crack damage degree discrimination.

Author Contributions

Writing—original draft preparation, methodology, L.C.; validation and editing, X.Z.; data curation and editing, L.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under Grant No. 51665004, in part by the Natural Science Foundation of Xinjiang Uygur Autonomous Region under Grant 2021B01003-1, in part by the Natural Science Foundation of Xinjiang Uygur Autonomous Region under Grant 2021D01C050.

Institutional Review Board Statement

Not application.

Informed Consent Statement

Not application.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McNames, J. Fourier series analysis of epicyclic gearbox vibration. J. Vib. Acoust. 2001, 124, 150–152. [Google Scholar] [CrossRef]
  2. Cammalleri, M.; Castellano, A. Analysis of hybrid vehicle transmissions with any number of modes and planetary gearing: Kinematics, power flows, mechanical power losses. Mech. Mach. Theory 2021, 162, 104350. [Google Scholar] [CrossRef]
  3. Chaari, F.; Fakhfakh, T.; Haddar, M. Analytical investigation on the effect of gear teeth faults on the dynamic response of a planetary gear set. Noise Vib. Worldw. 2009, 37, 9–15. [Google Scholar] [CrossRef]
  4. Chaari, F.; Fakhfakh, T.; Hbaieb, R.; Louati, J.; Haddar, M. Inflfluence of manufacturing errors on the dynamic behavior of planetary gears. Int. J. Adv. Manuf. Technol. 2006, 27, 738–746. [Google Scholar] [CrossRef]
  5. Park, S.; Lee, J.; Moon, U.; Kim, D. Failure analysis of a planetary gear carrier of 1200HP transmission. Eng. Fail. Anal. 2010, 17, 521–529. [Google Scholar] [CrossRef]
  6. Yuksel, C.; Kahraman, A. Dynamic tooth loads of planetary gear sets having tooth profile wear. Mech. Mach. Theory 2004, 39, 695–715. [Google Scholar] [CrossRef]
  7. Han, H.Z.; Zhao, Z.F.; Tian, H.X. Fault feature analysis of planetary gear set influenced by cracked gear tooth and pass effect of the planet gears. Eng. Fail. Anal. 2021, 2021, 105162. [Google Scholar] [CrossRef]
  8. Yang, L.T.; Wang, L.M.; Yu, W.N.; Shao, Y.M. Investigation of tooth crack opening state on time varying meshing stiffness and dynamic response of spur gear pair. Eng. Fail. Anal. 2021, 121, 105181. [Google Scholar] [CrossRef]
  9. Li, Z.; Ma, H.; Feng, M.; Zhu, Y.; Wen, B. Meshing characteristics of spur gear pair under different crack types. Eng. Fail. Anal. 2017, 80, 123–140. [Google Scholar] [CrossRef]
  10. Barszcz, T.; Randall, R.B. Application of spectral kurtosis for detection of a tooth crack in the planetary gear of a wind turbine. Mech. Syst. Signal Process. 2009, 23, 1352–1365. [Google Scholar] [CrossRef]
  11. Vicuña, C.M. Theoretical frequency analysis of vibrations from planetary gearboxes. Forsch Ing. 2012, 76, 15–31. [Google Scholar]
  12. Feng, Z.P.; Zuo, M.J. Vibration signal models for fault diagnosis of planetary gearboxes. J. Vib. Acoust. 2012, 331, 4919–4939. [Google Scholar] [CrossRef]
  13. Feng, Z.P.; Zuo, M.J. Fault diagnosis of planetary gearboxes via torsional vibration signal analysis. Mech. Syst. Signal Process. 2013, 36, 401–421. [Google Scholar] [CrossRef]
  14. Wu, X.; Meagher, J.; Sommer, A. A differential planetary gear model with backlash and teeth damage. Shock Vib. 2011, 5, 203–215. [Google Scholar]
  15. Liang, X.; Zuo, M.J.; Pandey, M. Analytically evaluating the influence of crack on the mesh stiffness of a planetary gear set. Mech. Mach. Theory 2014, 76, 20–38. [Google Scholar] [CrossRef]
  16. Lei, Y.G.; Lu, X.; Liu, Z.Y. Vibration Signal Simulation and Fault Diagnosis of Planetary Gearboxes Based on Transmission Mechanism Analysis. J. Mech. Eng. 2016, 52, 111–122. [Google Scholar] [CrossRef]
  17. Xiao, Z.M.; Huan, L.R.; Cao, J.X. Influence of crack fault on the dynamical characteristics of a planetary gear transmission system. J. Vib. Shock 2020, 39, 188–194. [Google Scholar]
  18. Tatar, A.; Schwingshackl, C.W.; Friswell, M.I. Modal sensitivity of three-dimensional planetary geared rotor systems to planet gear parameters. Appl. Math. Model. 2023, 113, 309–332. [Google Scholar] [CrossRef]
  19. Sanchez-Espiga, J.; Fernandez-del-Rincon, A.; Iglesias, M.; Viadero, F. Planetary gear transmissions load sharing measurement from tooth root strains: Numerical evaluation of mesh phasing influence. Mech. Mach. Theory 2021, 163, 104370. [Google Scholar] [CrossRef]
  20. Antoni, J. The spectral kurtosis: A useful tool for characterising non-stationary signals. Mech. Syst. Signal Process. 2006, 20, 282–307. [Google Scholar] [CrossRef]
  21. Dwyer, R. Detection of non-Gaussian signals by frequency domain kurtosis estimation. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP’83), Boston, MA, USA, 14–16 April 1983; Volume 8, pp. 607–610. [Google Scholar]
  22. Antoni, J. Fast computation of the kurtogram for the detection of transient faults. Mech. Syst. Signal Process. 2007, 21, 108–124. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the planetary gear train structure.
Figure 1. Schematic diagram of the planetary gear train structure.
Sensors 23 00253 g001
Figure 2. System dynamics model.
Figure 2. System dynamics model.
Sensors 23 00253 g002
Figure 3. Meshing stiffness solution model.
Figure 3. Meshing stiffness solution model.
Sensors 23 00253 g003
Figure 4. Meshing stiffness under normal conditions.
Figure 4. Meshing stiffness under normal conditions.
Sensors 23 00253 g004
Figure 5. The meshing stiffness of inner gear ring and planetary gear.
Figure 5. The meshing stiffness of inner gear ring and planetary gear.
Sensors 23 00253 g005
Figure 6. Meshing root crack model.
Figure 6. Meshing root crack model.
Sensors 23 00253 g006
Figure 7. Meshing stiffness under different fault states.
Figure 7. Meshing stiffness under different fault states.
Sensors 23 00253 g007
Figure 8. (a) Time domain history of the normal dynamic loads. (b) Time domain history of the shallow crack dynamic load. (c) Time domain history of the deep crack dynamic excitation.
Figure 8. (a) Time domain history of the normal dynamic loads. (b) Time domain history of the shallow crack dynamic load. (c) Time domain history of the deep crack dynamic excitation.
Sensors 23 00253 g008
Figure 9. (a) Frequency domain history of the normal dynamic loads. (b) Comparison diagram of the dynamic excitation frequency domain under crack fault.
Figure 9. (a) Frequency domain history of the normal dynamic loads. (b) Comparison diagram of the dynamic excitation frequency domain under crack fault.
Sensors 23 00253 g009
Figure 10. Finite element model and schematic diagram of dynamic load application.
Figure 10. Finite element model and schematic diagram of dynamic load application.
Sensors 23 00253 g010
Figure 11. Wind turbine test bench and signal acquisition system.
Figure 11. Wind turbine test bench and signal acquisition system.
Sensors 23 00253 g011
Figure 12. (a) Simulation and experiment signal time domain history. (b) Enlarged diagram of the meshing area.
Figure 12. (a) Simulation and experiment signal time domain history. (b) Enlarged diagram of the meshing area.
Sensors 23 00253 g012aSensors 23 00253 g012b
Figure 13. Flow chart of the dynamic strain solution of the inner ring under the crack fault of the solar wheel.
Figure 13. Flow chart of the dynamic strain solution of the inner ring under the crack fault of the solar wheel.
Sensors 23 00253 g013
Figure 14. (a) Time domain history of the gear ring root strain under normal working conditions. (b) Time domain history of tooth root strain of sun gear under crack condition.
Figure 14. (a) Time domain history of the gear ring root strain under normal working conditions. (b) Time domain history of tooth root strain of sun gear under crack condition.
Sensors 23 00253 g014
Figure 15. (a) The planetary gear is about to snap into the faulty tooth inner ring deformation. (b) Planetary gear meshing fault inner ring deformation. (c) The inner ring is deformed before the planetary gear gnawing the faulty tooth.
Figure 15. (a) The planetary gear is about to snap into the faulty tooth inner ring deformation. (b) Planetary gear meshing fault inner ring deformation. (c) The inner ring is deformed before the planetary gear gnawing the faulty tooth.
Sensors 23 00253 g015aSensors 23 00253 g015b
Figure 16. (a) The stress fluctuation curve of the fault tooth is about to be inserted. (b) Strain fluctuation curve of the snagging fault tooth. (c) Strain fluctuation curve of the gnawed fault tooth.
Figure 16. (a) The stress fluctuation curve of the fault tooth is about to be inserted. (b) Strain fluctuation curve of the snagging fault tooth. (c) Strain fluctuation curve of the gnawed fault tooth.
Sensors 23 00253 g016aSensors 23 00253 g016b
Figure 17. Flow chart of dynamic strain solution for inner ring root crack fault.
Figure 17. Flow chart of dynamic strain solution for inner ring root crack fault.
Sensors 23 00253 g017
Figure 18. Diagram of tooth root crack model of the inner ring.
Figure 18. Diagram of tooth root crack model of the inner ring.
Sensors 23 00253 g018
Figure 19. Time-domain history diagram of the strain signal of the inner ring tooth root crack.
Figure 19. Time-domain history diagram of the strain signal of the inner ring tooth root crack.
Sensors 23 00253 g019
Figure 20. Enlarged view of the partial meshing deformation zone.
Figure 20. Enlarged view of the partial meshing deformation zone.
Sensors 23 00253 g020
Figure 21. (a) State 1. (b) State 2. (c) State 3. (d) State 4.
Figure 21. (a) State 1. (b) State 2. (c) State 3. (d) State 4.
Sensors 23 00253 g021
Figure 22. (a) State 2 strain wave curve. (b) State 3 strain wave curve. (c) State 4 strain wave curve.
Figure 22. (a) State 2 strain wave curve. (b) State 3 strain wave curve. (c) State 4 strain wave curve.
Sensors 23 00253 g022
Figure 23. Damage discrimination flow chart.
Figure 23. Damage discrimination flow chart.
Sensors 23 00253 g023
Figure 24. Fast spectral kurtograph.
Figure 24. Fast spectral kurtograph.
Sensors 23 00253 g024
Figure 25. (a) Sun gear tooth crack fault signal. (b) Inner ring crack fault signal.
Figure 25. (a) Sun gear tooth crack fault signal. (b) Inner ring crack fault signal.
Sensors 23 00253 g025
Figure 26. Complexity in different states.
Figure 26. Complexity in different states.
Sensors 23 00253 g026
Figure 27. Root-mean-square values in different states.
Figure 27. Root-mean-square values in different states.
Sensors 23 00253 g027
Table 1. Planetary gear transmission system parameters.
Table 1. Planetary gear transmission system parameters.
NameSunPlanetRing
Tooth number2836100
Tooth width/mm101010
Modulus/mm111
Pressure angle/°202020
Mass/kg0.03560.05440.432
Moment of inertia/(kg·m2)0.000470.001210.143
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, L.; Zhang, X.; Wang, L. Research on Root Strain Response Characteristics of Inner Ring of Planetary Gear Transmission System with Crack Fault. Sensors 2023, 23, 253. https://doi.org/10.3390/s23010253

AMA Style

Chen L, Zhang X, Wang L. Research on Root Strain Response Characteristics of Inner Ring of Planetary Gear Transmission System with Crack Fault. Sensors. 2023; 23(1):253. https://doi.org/10.3390/s23010253

Chicago/Turabian Style

Chen, Lan, Xiangfeng Zhang, and Lizhong Wang. 2023. "Research on Root Strain Response Characteristics of Inner Ring of Planetary Gear Transmission System with Crack Fault" Sensors 23, no. 1: 253. https://doi.org/10.3390/s23010253

APA Style

Chen, L., Zhang, X., & Wang, L. (2023). Research on Root Strain Response Characteristics of Inner Ring of Planetary Gear Transmission System with Crack Fault. Sensors, 23(1), 253. https://doi.org/10.3390/s23010253

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

Article Metrics

Back to TopTop