Research Article  Open Access
Research on Fault Evolution Feature Extraction and Identification of Sun Gear Cracks in Planetary Gearbox Based on Volterra HighOrder Kernel Generalized Frequency Response Graphic Analysis
Abstract
In this paper, relying on the Volterra series nonlinear system model and the highorder kernel Hilbert’s reconstructed kernel fast solved algorithm, a fault feature frequency domain identification method based on Volterra highorder kernel generalized frequency response graph analysis is proposed. Firstly, the method uses the system input and output vibration signals to determine the Volterra model. Then, the Volterra highorder kernel function is solved quickly by reproducing kernel Hilbert space method, and the generalized frequency response function is used to identify the model. Finally, multidimensional highorder spectral pattern analysis is used to separate and extract the fault and degree characteristic information implied by frequency and phase coupling in the thirdorder kernel function. Following the theoretical approach, in the experimental part, this paper uses the planetary gearbox fault loading test rig to complete the data collection and establishes the Volterra experimental model through the measured data. The generalized frequency responses of each order kernel function are compared and analyzed and the capability of distinguishing and the adaptability of different order kernel functions for the degree of crack failure are discussed. The effects of changing the memory length of the Volterra model and the order of the kernel function on the recognition result are verified. The final experimental results show that the use of reproducing kernel Hilbert space can effectively avoid the dimension disaster problem that occurs in the highorder kernel solution process. Moreover, the thirdorder kernel can describe more intuitively the nonlinear system model under multifactor coupling than the secondorder kernel. Finally, Volterra series model the thirdorder kernel’s generalized frequency response can effectively distinguish between nondefective and faulty gears, and its resolution is enough to distinguish the degree of failure of gear cracks.
1. Introduction
Among the advantages of planetary gear transmission system, high transmission ratio, strong load capacity, small size, and light weight contributed to its success in many industrial applications such as wind power equipment, engineering machinery, and aerospace [1, 2]. However, planetary gear transmission systems are subjected to hostile longlasting working conditions, e.g., variable load, frequent start and stop, and high speed, which inevitably leads to different types of failures of the critical parts of planetary gearboxes [3].
Statistical analysis highlighted that the gear crack plays an important role among the most important planetary gearbox failure factors [4]. Moreover, the gear crack easily triggers other types of failures, eventually leading to a chain reaction from small to major faults. Therefore, timely diagnosis and prediction of early failures of gear cracks, advanced control and resolution of crack propagation, and timely planned maintenance will effectively avoid or reduce the probability of a complete machine accident caused by a planetary gearbox crack failure. Therefore, the study and establishment of a set of feature extraction and prediction techniques for early and weak failures of planetary gearbox gears are significant [5].
In general, when the internal gear of the planetary gearbox fails, most of the vibration signals of the box body show obvious nonstationary [6–8], nonlinear, and complex features [9]. Therefore, the key to this paper is to choose a nonlinear analysis method that is suitable for the system characteristics to pick up useful information from the vibration signal and complete the crack fault feature identification. At present, among the commonly used analysis technique for fault diagnosis of planetary gearboxes there are new cepstrum method [10], spectral kurtosis (SK) method [11], wavelet transform (WT) [12, 13], singular value decomposition (SVD) [14], empirical mode decomposition (EMD) [15], and adaptive stochastic resonance (ASR) [16, 17]. Some important results have been achieved through research on fault feature identification using the above methods. Barzcz et al. proved the probability of applying spectral kurtosis to detect toothlike cracks on planetary gear rings of wind turbines [11]. Similarly, Sayena et al. used intricate Morlet wavelets to extract features that can distinguish between healthy carriers and faults in planetary gearboxes [18]. Jiang et al. propose a signal denoising technology based on adaptive Morlet wavelet and singular value decomposition (SVD) and applied it to extract pulse characteristics of second stage reducer planetary gearbox for wind turbine [19]. Based on a set of EMD and an energy segmentation algorithm, Feng et al. proposed a combined amplitude and frequency demodulation technology to diagnose sun gear failures using a planetary gearbox test bench [20]. Further to this, Lei et al. discussed an adaptive stochastic resonance technology to diagnose sun gear faults [15]. From the direction of the above research results, many studies have focused on the qualitative analysis of typical fault types of planetary gearboxes. Of course, it does not involve the light or severity of the fault, nor does it provide a more indepth study of the mechanism and crack evolution of the fault crack.
More recently, many researches focused on investigating the degrees of failure. Ren et al. proposed a method combining the threedimensional waterfall spectrum and a variable scale wavelet to analyze timefrequency domain characteristics of three different crack depths in a rotor [21]. Moreover, to calculate the timevarying mesh stiffness, Jiateng et al. proposed an analyticalfinite element method called the assiststress intensity factor gear contact model [22]. Haitao Wang et al. attempted to establish a nonlinear system model for rolling bearings using the Volterra series theory and used a multipulse excitation algorithm to quickly estimate the Volterra secondorder kernel and their generalized frequency responses [23]. Finally, the highresolution feature extraction and fault degree identification of the inner ring, outer ring, and rolling element pitting failure of the rolling bearing are finally achieved by analyzing the numerical characteristics of the broadspectrum frequency response and bispectrum diagonal slice. Although the above studies have achieved specific results in the identification of faults, they still focus on the failure of a single component and do not fully meet the actual engineering needs, that asks for the extraction and prediction of the vibration signal characteristics of the cracks and the degree of evolution of planetary gearbox gears in various mechanical equipment under actual conditions.
In response to the aboveunsolved problems and considering our previous research results, this paper proposes a novel identification method for crack fault evolutionary characteristics of planetary gearboxes based on the Volterra highorder kernel function theory [24, 25]. The main academic ideas of this method can be summarized as follows: First, a nonlinear Volterra system model of the system is established by providing vibration signals of planetary gearboxes as the input and output; second, through the rational construction of the Volterra series Hilbert reproducing kernels [26–28], the highorder kernels of the Volterra model are calculated. Avoiding the “dimensional disaster” problem that may occur [29], third, through the generalized frequency response characteristics of highorder kernel of the system Volterra model, the crack faults and degree characteristics of planetary gearbox gears are identified efficiently. Finally, several fault evolution feature extraction experiments based on planetary gearbox fault loading experimental platform are simulated and carried out through numerical examples. Further, the correctness, universality, and superiority of the proposed method are verified proving that it further enriches and promotes fundamental research in the field of mechanical and electrical equipment fault diagnosis.
2. Theoretical Background
2.1. Volterra Series Theory
The Volterra series theory was first proposed from the perspective of mathematics and used to study some integrodifferential equations. However, starting from the basic theory of signals and systems, according to the Volterra series, the output of a system can be expressed in terms of an infinitelength algebraic sum of multidimensional convolution integrals of its input values, thereby defining the system model in a specific and unique way.
For a linear system, the system input and output relations can be expressed in the time domain as the following convolution operations:
where , are system input and output, respectively, and is impulse response transfer function (IRF).
In the same way, for any continuous nonlinear system, if the input energy of the system is limited under the zero initial condition, the Volterra series represents the response of the system.
where is kernel of the order Volterra series and it is the generalization of impulse response transfer function (GIRF) of the linear system in the highdimensional space. By comparing the linear system with the nonlinear system representation, with regard to linear system, represents the system property of the linear system, that is, the system transfer function; similarly, for a nonlinear system, the set of , , represents the system properties of nonlinear system. Therefore, one can infer that the system properties of the nonlinear system can be determined and represented by the kernel function of the Volterra nonlinear model of the system.
If formula (refer to formula (2)) is discretized in the time domain, a discrete time domain Volterra expression can be obtained:where , is the system input and output, and ) is the th order kernel function of the discrete Volterra series.
To study the kernel function of Volterra series in more depth, the analysis in the time domain is not enough. Therefore, the concept of generalized frequency response function (GFRF) is introduced to realize the research of kernel function in the frequency domain. It is defined as the multidimensional Fourier transform of the kernel function of Volterra series, as follows:
Therefore, the higherorder kernel functions of Volterra nonlinear system model represent the basic characteristics of the nonlinear system. At this stage of the analysis, how to quickly and efficiently determine the order Volterra kernel functions of the system to be identified becomes another important and difficult point. Searching in the literature for the solution of the Volterra kernel function shows that most of the calculation methods mainly adopt traditional approximation algorithms such as the leastsquares method, multitone excitation, and neural network. This will inevitably expose the deficiencies of this type of method, such as large computational complexity, high complexity, and highorder kernels that are prone to dimension disasters. Therefore, for solving the highorder Volterra kernel function of a nonlinear coupled system under complex conditions, such as the complex planetary gearboxes, this paper proposes a novel Volterra series kernel method based on reproducing kernel Hilbert space (RKHS). This calculation method can more easily determine the Volterra kernel function of nonlinear systems, especially highorder kernel.
2.2. Improved Calculation Method for Solving HighOrder Kernel
2.2.1. Volterra Series in Hilbert Space
The discrete form of the bounded Volterra series is defined as where m is length of the system memory and t is time. Let the nonlinear system have an input vector U and ask to compute the Volterra series kernel to the pth order.
Define the spatial projection as If we let , different nonlinear systems can only correspond to a unique set of Volterra kernel vectors, such as . If is given, the Volterra series can be represented as a compact linear form .
Let all the nonlinear timeinvariant dynamic systems (LTIs) that can be modeled by the Volterra series be the space H, then the Volterra series can be expressed aswhere , is spatial projection vector input by the system, , i.e., is an infinite series.Since a linear space can be called a Hilbert space after establishing an inner product, the space at this time is Hilbert space.
If the input function in Hilbert space is , then there isThe input function is constructed from the spatial projection η(u) of the original system input value . It can be shown that function also has the following two important properties.
Fix v, then is an element in the Hilbert space H, where (·) represents any one of input vector.
Proof: Let , , then (11) can be transformed to The control can be obtained .
For any element in the Hilbert space H, , and its inner product are values where is at .
Proof:
According to the kernel Hilbert space theory, the linear functional of the Hilbert space with both properties mentioned above is called a “reproducing kernel Hilbert space” (RKHS). According to the RKHS theory, the nonlinear system can be described as an element in Hilbert space; assuming some input values and corresponding output values of the nonlinear system were known, then the system can be reconstructed. It is a coincidence that this is the process of modeling nonlinear systems using the Volterra series theory, proving that there is interchangeability between the two. Moreover, this important theoretical discovery lays on solid mathematical basis for solving the kernel function of the nonlinear system Volterra model.
Assume that the system input value is ) and the observation value after passing through the nonlinear system is . Then, according to the second property of the RKHS, it is known that . The solution for in (13) requires the value of is known for any input, as in
According to linear space theory, a nonlinear model to be identified corresponds to the output value of the input , which is the projection of the response of any input signal of the system in , , the linear subspace into which the system responds. According to the Hilbert space theory of reproducing kernels, the better approximation of the Volterra series representation of a nonlinear model can be represented by a linear combination of , , i.e.,
2.2.2. Explicit Expressions of Each Intermediate Parameter and Volterra Series Kernel
This paper replaces the (·) in (15) with for a more clear and visual representation of the Volterra series kernel .
ifAmong them there is . By rewriting (16) into a matrix, we get
From the solutionwhere represents the generalized inverse of the matrix . According to (16),
After comparing the coefficients of (20) and (16), we can get the Volterra series kernel as
This solution highlights that improved method for solving the kernel of the Volterra series model based on RKHS avoids the exponentially explosive dimensionality of the iterations of conventional calculation methods and also solves Volterra highorder kernel problems, providing a new way of thinking.
2.2.3. Examples and Simulations
Let the excitation signal input by the system be a random signal uniformly distributed over the range . Observe 500 sets of output data. Then use the Volterra series expansion formula identified by the proposed method as
(1) Comparison of Time Domain Output. For nonlinear systems, this paper applies excitation signals
The number of sampling points was 120, and the time domain response of the system output based on the simulation is shown in Figure 1; except for some points in the middle and the tail of the signal, the system well approximates the output at most sampling points.
(2) Comparison of Generalized Frequency Response of FirstOrder Kernel Functions. Because the effect of the DC component on the amplitudefrequency response is the same as that of the phasefrequency response, only the firstorder kernel exists in the Volterra series of the nonlinear system when the zeroorder or highorder Volterra kernels are ignored. Therefore, in the case of neglecting these Volterra kernels, only the firstorder kernel exists in the Volterra series of nonlinear systems. The comparison results of the frequency response of the firstorder kernel functions of an ideal system with those of the identified model are shown in Figure 2 and Figure 3, respectively. Figures show that the amplitude error at any point never exceeds 1.2 dB, while the phase angle error never exceeds 4 degrees.
(3) Comparison of Generalized Frequency Response of SecondOrder Kernel Function. The comparison between the firstorder kernel function amplitudefrequency response and the firstorder kernel function phasefrequency response of the above system can be intuitively obtained. The system identification result has a good identification accuracy from the perspective of frequency response analysis. However, because the identified system is a nonlinear system, the frequency response of the firstorder kernel function can only reflect the linear characteristics of the system. Therefore, the similarity of the frequency response of the firstorder kernel function does not mean that the identified model can completely represent the original system. For this reason, the input signal and the output signal in this paper are used to solve the secondorder kernel function amplitudefrequency response of the system (Figure 4(a)) and the secondorder kernel function phasefrequency response (Figure 4(b)).
(a) 3D diagram of amplitudefrequency response curve of Volterra secondorder
(b) 3D diagram of phasefrequency response curve of Volterra secondorder
The frequency response of the secondorder kernel function of the system obtained by reproducing kernel Hilbert space method is shown in Figure 5(a) and its phasefrequency response in Figure 5(b). The amplitude frequency and phasefrequency errors between the identification model and the real system model are shown in Figures 6(a) and 6(b).
(a) 3D diagram of amplitudefrequency response curve of Volterra secondorder
(b) 3D diagram of phasefrequency response curve of Volterra secondorder
(a) Magnitude response error
(b) Phasefrequency response error
To sum up, the correctness, generalization, and universality of the Volterra kernel function algorithm based on RKHS are proved by comparison and analysis of the firstorder kernel, secondorder kernel amplitude frequency, and phasefrequency response of the nonlinear simulation example model. By analogy, according to the calculation method described in this paper, the thirdorder kernel, fourthorder kernel, and even higher orders of the nonlinear system Volterra model can be solved, thus providing a reliable and complete mathematical theoretical basis for the subsequent acquisition, extraction, and tracking algorithm for the evolution characteristics of the early weak crack faults in the planetary gearbox gear.
Section 3 presents two basic studies on the mechanism of crack failure of the planetary gearbox and the theoretical model of planetary gear system and also discusses the influence of the crack evolution on the frequency response of the dynamic model.
3. Failure Mechanism Analysis
3.1. Planetary Gear System Dynamic Modeling
According to the dynamic model of planetary gear transmission system shown in Figure 7, the system consists of a sun gear, a planet carrier, a ring gear, and three planetary gears. , , and are rotation angle of the planet carrier, the sun gear, and the planet wheels, respectively; , , , are the meshing stiffness and meshing damping coefficient of the sun gear, the ring gear, and the planet gears, respectively; , are the meshing gear pair side clearance; and are input torque and load torque, respectively.
From the Lagrange equation, the differential equation for the dynamic model shown in Figure 7 can be obtained:where is dynamic moment of inertia of the planet gear ; is moment of inertia of the planet carrier; is base circle radius of the sun wheel; is radius of the planet’s base circle; is planet carrier radius; is quality of the planet’s wheels; is pressure angle; , are integrated meshing error of each gear pair; represents a gap nonlinear function and can be expressed asTo facilitate the solution, introduce relative displacement coordinates:Defining the time nominal scale and the displacement nominal scale , the dimensionaldynamic differential equation is
where ; is planet carrier equivalent radius; is meshing damping ratio of planet i and sun gear; is meshing damping ratio of the planetary gear i and the ring gear; and are equivalent masses of the sun gear and the planet gear, respectively, on their respective base circle radii.
Research on the meshing stiffness of the sun gear failure shows that when the gear is slightly broken, it will cause a change in the meshing stiffness. Therefore, the gear failure factor is introduced. In the normal state , the severe fault condition of the broken tooth , the slight fault state of the gear crack , and the timevarying meshing stiffness of the meshing gear pair is uniformly expressed as
Because the crack fault of the gear is mainly presented on the equation of the meshing stiffness, the relevant research on the meshing stiffness of the fault of a gear broken tooth by Shi Lichen et al. is used for reference [30]. In this paper, take and , respectively, to establish the sun gear wheel crack failure meshing stiffness equation, then its expression is as follows:
3.2. Dynamic Model Response Analysis
Based on the established planetary gearbox system dynamics model described in Section 3.1 and on the system parameters listed in Tables 1 and 2, by solving and simulating the dynamic equations, the model response spectrum of the sun gear broken tooth can be obtained as shown in Figures 8 and 9.


When a 0.2 mm gear crack fault occurs in the sun gear (Figure 8), due to the amplitude modulation effect of the planetary gear revolution, the frequency spectrum in the spectrum diagram appears as the side band of the planet carrier transponder ; i.e., the amplitude at the site is more obvious (being , an integer). For example, consider the normalized frequency . The peak at ≈ 0.008916, is large, and the side frequency with the sun gear failure frequency as the modulation frequency also appears in the spectrum. If the normalized frequency , ≈ 0.009541, , and at the peak, side band is obvious. That is, the amplitude in the spectrogram appears at the normalized frequency (being , , an integer). Comparing the results with those of Figure 9, when a 0.4 mm gear crack fault occurs in the sun gear, the crack fault frequency calculated by the theoretical model is basically very close to the characteristic frequency at 0.2 mm fault. The maximum normalized frequency discrimination does not exceed 0.0003, so a highresolution and reasonable discrimination of the difference in the degree of minor failure of the crack cannot be made.
At this stage, by summarizing the mechanism analysis of the cracks and the results of the frequency characteristics of the dynamic response of different degrees of failure, one can note that different degrees of early microcrack fault eigen frequency have very high similarity and low frequency feature resolution, which is not enough to effectively distinguish faults.
Therefore, in Section 4, this paper carries out an experimental study of the fault feature extraction and identification method for planetary gearbox gears based on Volterra highorder kernel generalized frequency response graphical analysis. The experimental results are compared and analyzed and they prove that the proposed fault diagnosis method has certain innovation and advantages in the identification of crack fault feature resolution and degree.
4. Experimental Study
The experimental study was carried out using the planetary gearbox fault diagnosis test bench shown in Figure 10, and the detailed parameters of the planetary gearbox are listed in Table 1. In order to simulate the gear crack fault, two cracks of different depths (0.2 and 0.4 mm) were machined on one of the gears of the first stage sun gear. In the experiment, two accelerometers installed on the motor shell by using a magnetic base acquired the vibration signal. The acceleration sensors are installed at the driving end of the motor shell and the 12 o’clock position of the fan end, respectively. The load of the driving end is 0, the drive end speed is 1979 r/min, and the vibration signal enters the 32bit data acquisition card through the signal conditioning unit to collect the highspeed data. The experimental observation data is sampled and 4000 sets of data are collected to carry on the experimental analysis and processing.
The experimental technology program continues the path developed in the previous research [23]. That is, Hilbert’s reproducing kernel space method is first used to calculate and solve the input and output data. The firstorder and secondorder kernel functions of the Volterra model of planetary gearbox nondefective gear and that of the two faulted gears are obtained, respectively. Moreover, use the kernel function to draw the bispectrum threedimensional map of the two gears shown in Figure 11.
(a) Nondefective gear
(b) 0.2 mm crack depth
(c) 0.4 mm crack depth
The comparison of the results in Figure 11 shows that there is a small difference in the shape of the peaks and the energy distribution among the bispectrum 3D maps of the nondefective gear and that of defective gears. Moreover, because of the too low resolution, the curves do not give enough information about fault feature information and details of its fault evolution. Therefore, to observe the features shown in Figure 11, diagonal slices are reduced in dimension obtaining a twodimensional map (see Figure 12). The map highlights that the frequencies of fault features corresponding to nondefective gears and defective gears with the 0.2 mm and 0.4 mm cracks are all around 70 Hz, and effective distinctions among defects cannot be made. It can be seen that the spectra of secondorder kernel and diagonal slice diagrams are difficult to provide a certain resolution for distinguishing crack failures. The relevant conclusions obtained in the failure mechanism analysis of the third part of this paper are verified. The above experiments can initially prove that in the process of modeling and analyzing Volterra of a more complex nonlinear system such as a planetary gearbox, because of the uncertainty of the coupling factors between the components and the transmission path in the gearbox, the Volterra model’s secondorder kernel function and its corresponding generalized frequency response eventually lose the ability to identify faults.
(a) Nondefective gear
(b) 0.2 mm crack depth
(c) 0.4 mm crack depth
The comparison of results in Figure 12 shows that although the characteristic frequency points can be captured, the correlation calculation result graph cannot extract fault features efficiently. Combining the characteristics of the Volterra series itself, there are two variable parameters for the nonlinear coupling resolution of the Volterra model built by the system to adjust its characteristics: one is the model memory length and the other is the model order . In this research, the change of the length of the memory, from to and , was tried first and the calculation results are given in Figures 13–16.
(a) Nondefective gear
(b) 0.2 mm crack depth
(c) 0.4 mm crack depth
(a) Nondefective gear
(b) 0.2 mm crack depth
(c) 0.4 mm crack depth
(a) Nondefective gear
(b) 0.2 mm crack depth
(c) 0.4 mm crack depth
(a) Nondefective gear
(b) 0.2 mm crack depth
(c) 0.4 mm crack depth
The plots show that changing the memory length improves the gear fault identification of planetary gearboxes slightly because the change highlights the fault frequency point, which corresponds to the highest peak. However, the amplitude of other noise frequency points also increases; because of its irregular and random behavior, the equipment failure still cannot be identified. Therefore, one can conclude that changing the system memory length of the Volterra model is also not ideal for feature extraction and identification of gear cracks in planetary gearboxes.
In summary, the study readjusted the experimental scheme and tried to increase the Volterra order from 2 to 3. However, the problems faced are very serious; that is, the computational complexity of highorder Volterra series kernels becomes more difficult with the increase of order. The dimensions of its parameters and computational data increase exponentially and fall into the vortex of dimension disasters. In order to solve the problem of dimensional disasters, this paper used the reproducing kernel Hilbert space method to calculate the Volterra series higherorder kernel. Data acquisition for input and output of the system is carried out before the experiment is extended. The thirdorder kernel of the Volterra nonlinear model of the planetary gearbox is obtained by commercial software MATLAB, and the generalized frequency response graph of the thirdorder kernel is drawn in a 4D plot, using the dimensionality reduction technique to visualize the results (see Figure 17). This technique allows to extract and compare the characteristics of gear crack faults through the generalized frequency response graph of the thirdorder kernel so that the failure coupling characteristics of Volterra highorder kernel functions can be expressed more clearly and intuitively.
(a) Nondefective gear
(b) 0.2 mm crack depth
(c) 0.4 mm crack depth
By observing the evolution of the Volterra thirdorder kernel generalized frequency response of the epicyclic gearbox as shown in Figure 17, it can be seen that the highorder spectra of nondefective gear and the highorder spectra with gear cracks of 0.2 mm and 0.4 mm are very different. Comparing the plots in Figures 17(a), 17(b), and 17(c), some fault features emerge. For example, as can be seen from the whole, the shape of the nondefective gear is similar to a wavy, reversed cone, but that of the defective gear is like a column. Moreover, the colors highlight that the highorder spectral set in nondefective gear tends to be light red and yellow, whereas those of gears with crack failure are mainly blue and there are a few dark red areas.
5. Discussion of Experimental Results
The experimental part is divided into three stages.
In the first stage, the main research results and ideas of the research group were applied. The experimental analysis of the cracks in the planetary gearbox is performed through the Volterra secondorder kernel generalized frequency response characteristics. Although the algorithm can obtain the frequency response graph of the system, the characteristics of the crack failure and evolution cannot be extracted clearly from the feature analysis of diagonal slice graphics. The secondorder kernel fault characteristics of the cracks and evolution of gears in planetary gearboxes do not fully comply with the frequency response characteristics of pitting failures in the inner and outer rings of rolling bearings because of the ambiguity of the frequency of the fault. Moreover, there is a broad, random interference which causes significant difficulties in fault feature extraction. These issues show that the secondorder kernel generalized frequency response of the Volterra model is not suitable for fault diagnosis analysis of nonlinear coupled system like a planetary gearbox.
The second stage mainly analyzes the reasons why the first phase of the experiment did not meet the expectations of the study and its solutions. Through a new analysis of the Volterra series theory, this paper found an improved method for improving the identification resolution of the Volterra model, that is, increase the memory length of Volterra model and the order of the kernel function. In subsequent experiments, the memory length was first scaled up. The experimental results show that the increase of memory length contributes to the nonlinear coupling characteristics of Volterra secondorder kernel function and its generalized frequency response, but some shortcomings such as large frequencycoupled crosstalk randomness still exist.
In the third stage, mainly based on the summarization of relevant experiences in the first and second stages, a tentative study was conducted by increasing the order of kernel function n. The experimental scheme has been improved from the original secondorder kernel to the thirdorder kernel, but subject to the constraints of traditional algorithms, the thirdorder kernel calculation is more difficult. Therefore, in this experimental scheme, the thirdorder kernel solution of the planetary gearbox crack failure Volterra model was solved using the related results based on the improved calculation method of the Hilbert regenerative kernel studied in Section 2. By drawing the dimensionality reduction curve of higher order kernel generalized frequency response graph, the fault characteristics are more clearly defined, and it is found that the graph has higher resolution.
Combining all the above experimental contents, one can conclude that the Volterra highorder kernel generalized frequency response method can effectively extract fault features of planetary gearbox cracks and it is also possible to perform highresolution graphical analysis on the evolution of early weak faults, thus introducing a new method of fault diagnosis.
6. Conclusion
This paper addresses a series of nonlinear problems related to the diagnosis of gear cracks in planetary gearboxes. By introducing Volterra nonlinear system modeling theory and proposing a highorder kernel improved algorithm based on reproducing kernel Hilbert space, the problem of dimension disasters is easily avoided in the process of solving highorder kernel functions of the Volterra model. Based on the failure mechanism and experimental comparison of the measured data of the planetary gearbox experimental loading platform, the research highlights that the Volterra secondorder kernel function does not have enough capacity to represent a nonlinear frequency coupling system. Through repeated attempts in the experiment, the experimental scheme is adjusted for increasing the memory length of the Volterra model and increasing the order of the kernel function. The experiments show that the method of improving the model order is more suitable for the acquisition of crack fault characteristics of planetary gearbox gears under actual conditions than increasing the length of Volterra memory. A comparative test between a nondefective gear and two defective gears with crack depths of 0.2 mm and 0.4 mm, respectively, proved the conclusions. Finally, a new method for nonlinear frequency domain coupled feature identification of planetary gearboxes based on the generalized frequency response of higherorder kernel functions of Volterra series theory is proposed. This new method has very important theoretical and practical research significance because of the very high frequency response resolution of fault signatures in nonlinear systems.
Data Availability
The data [the data supporting the conclusions] used to support the findings of this study may be released upon application to the Xi’an University of Architecture and Technology, Electrical and Mechanical College, who can be contacted at 1639093204@qq.com.
Conflicts of Interest
The authors declare that are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China under Grant 51105292 and in part by the Shaanxi Provincial Natural Science Foundation of China under Grant 2017JQ5044.
References
 H. Niu, X. Zhang, X. Zhao, and C. Hou, “Action mechanism of planetary gearbox's typical faults on tooth root strain of ring gear,” Journal of Vibration and Shock, vol. 36, no. 9, pp. 253–260, 2017. View at: Google Scholar
 Y. Li, Y. Yang, G. Li, M. Xu, and W. Huang, “A fault diagnosis scheme for planetary gearboxes using modified multiscale symbolic dynamic entropy and mRMR feature selection,” Mechanical Systems and Signal Processing, vol. 91, pp. 295–312, 2017. View at: Publisher Site  Google Scholar
 C. Zhao and Z. Feng, “Application of multidomain sparse features for fault identification of planetary gearbox,” Measurement, vol. 104, pp. 169–179, 2017. View at: Publisher Site  Google Scholar
 X. Liang, M. J. Zuo, and Z. Feng, “Dynamic modeling of gearbox faults: A review,” Mechanical Systems and Signal Processing, vol. 98, pp. 852–876, 2018. View at: Publisher Site  Google Scholar
 Z. Feng and M. Liang, “Complex signal analysis for planetary gearbox fault diagnosis via shift invariant dictionary learning,” Measurement, vol. 90, pp. 382–395, 2016. View at: Publisher Site  Google Scholar
 A. Singh and A. Parey, “Gearbox fault diagnosis under nonstationary conditions with independent angular resampling technique applied to vibration and sound emission signals,” Applied Acoustics, 2017. View at: Google Scholar
 R. Zimroz and A. Bartkowiak, “Two simple multivariate procedures for monitoring planetary gearboxes in nonstationary operating conditions,” Mechanical Systems and Signal Processing, vol. 38, no. 1, pp. 237–247, 2013. View at: Publisher Site  Google Scholar
 Z. Feng, X. Chen, and M. Liang, “Joint envelope and frequency order spectrum analysis based on iterative generalized demodulation for planetary gearbox fault diagnosis under nonstationary conditions,” Mechanical Systems and Signal Processing, vol. 7677, pp. 242–264, 2016. View at: Publisher Site  Google Scholar
 K. Feng, K. Wang, Q. Ni, M. J. Zuo, and D. Wei, “A phase angle based diagnostic scheme to planetary gear faults diagnostics under nonstationary operational conditions,” Journal of Sound and Vibration, vol. 408, pp. 190–209, 2017. View at: Publisher Site  Google Scholar
 R. B. Randall, “A history of cepstrum analysis and its application to mechanical problems,” Mechanical Systems and Signal Processing, vol. 97, pp. 3–19, 2017. View at: Publisher Site  Google Scholar
 T. Barszcz and R. B. Randall, “Application of spectral kurtosis for detection of a tooth crack in the planetary gear of a wind turbine,” Mechanical Systems and Signal Processing, vol. 23, no. 4, pp. 1352–1365, 2009. View at: Publisher Site  Google Scholar
 P. D. Samuel and D. J. Pines, “Constrained adaptive lifting and the CAL4 metric for helicopter transmission diagnostics,” Journal of Sound and Vibration, vol. 319, no. 12, pp. 698–718, 2009. View at: Publisher Site  Google Scholar
 W. Teng, X. Ding, X. Zhang, Y. Liu, and Z. Ma, “Multifault detection and failure analysis of wind turbine gearbox using complex wavelet transform,” Journal of Renewable Energy, vol. 93, pp. 591–598, 2016. View at: Publisher Site  Google Scholar
 W.X. Yang and P. W. Tse, “Development of an advanced noise reduction method for vibration analysis based on singular value decomposition,” NDT & E International, vol. 36, no. 6, pp. 419–432, 2003. View at: Publisher Site  Google Scholar
 Y. G. Lei, D. Han, J. Lin, and Z. J. He, “Planetary gearbox fault diagnosis using an adaptive stochastic resonance method,” Mechanical Systems and Signal Processing, vol. 38, no. 1, pp. 113–124, 2013. View at: Publisher Site  Google Scholar
 Z. Wang, Z. Han, F. Gu, J. X. Gu, and S. Ning, “A novel procedure for diagnosing multiple faults in rotating machinery,” ISA Transactions, vol. 55, pp. 208–218, 2015. View at: Publisher Site  Google Scholar
 B. Liu, S. Riemenschneider, and Y. Xu, “Gearbox fault diagnosis using empirical mode decomposition and Hilbert spectrum,” Mechanical Systems and Signal Processing, vol. 20, no. 3, pp. 718–734, 2006. View at: Publisher Site  Google Scholar
 A. Saxena, B. Wu, and G. Vachtsevanos, “A methodology for analyzing vibration data from planetary gear systems using complex morlet wavelets,” in Proceedings of the American Control Conference (ACC '05), pp. 4730–4735, Portland, Ore, USA, June 2005. View at: Google Scholar
 Y. Jiang, B. Tang, Y. Qin, and W. Liu, “Feature extraction method of wind turbine based on adaptive Morlet wavelet and SVD,” Journal of Renewable Energy, vol. 36, no. 8, pp. 2146–2153, 2011. View at: Publisher Site  Google Scholar
 Z. P. Feng, M. Liang, Y. Zhang, and S. M. Hou, “Fault diagnosis for wind turbine planetary gearboxes via demodulation analysis based on ensemble empirical mode decomposition and energy separation,” Journal of Renewable Energy, vol. 47, pp. 112–126, 2012. View at: Publisher Site  Google Scholar
 Z. Ren, S. Zhou, E. Chunhui, M. Gong, B. Li, and B. Wen, “Crack fault diagnosis of rotor systems using wavelet transforms,” Computers and Electrical Engineering, vol. 45, pp. 33–41, 2015. View at: Publisher Site  Google Scholar
 J. Wu, Y. Yang, X. Yang, and J. Cheng, “Fault feature analysis of cracked gear based on LOD and analyticalFE method,” Mechanical Systems and Signal Processing, vol. 98, pp. 951–967, 2018. View at: Publisher Site  Google Scholar
 H. Wang, Z. Kang, L. Shi, K. Wang, and X. Zhang, “Research and application of volterra series theory in rolling bearing fault state feature extraction,” Journal of Vibroengineering, vol. 20, no. 1, pp. 189–201, 2018. View at: Publisher Site  Google Scholar
 Y. Wei, Y. Chen, and X. Shi, “A spectral collocation method for multidimensional nonlinear weakly singular Volterra integral equation,” Journal of Computational and Applied Mathematics, vol. 331, pp. 52–63, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 F. Mirzaee and S. F. Hoseini, “A new collocation approach for solving systems of highorder linear Volterra integrodifferential equations with variable coefficients,” Applied Mathematics and Computation, vol. 311, pp. 272–282, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 M. Mohammadi and R. Mokhtari, “Solving the generalized regularized long wave equation on the basis of a reproducing kernel space,” Journal of Computational and Applied Mathematics, vol. 235, no. 14, pp. 4003–4014, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 H. Beyrami, T. Lotfi, and K. Mahdiani, “Stability and error analysis of the reproducing kernel Hilbert space method for the solution of weakly singular Volterra integral equation on graded mesh,” Applied Numerical Mathematics, vol. 120, pp. 197–214, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 E. Babolian, S. Javadi, and E. Moradi, “Error analysis of reproducing kernel Hilbert space method for solving functional integral equations,” Journal of Computational and Applied Mathematics, vol. 300, pp. 300–311, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 J. Cao, L. Chen, J. Zhang, and W. Cao, “Fault diagnosis of complex system based on nonlinear frequency spectrum fusion,” Measurement, vol. 46, no. 1, pp. 125–131, 2013. View at: Publisher Site  Google Scholar
 S. Lichen, L. Kun, W. Haitao, and L. Yang, “Fault response characteristics of tractor planetary gearbox based on dynamical simulation and its validation,” Transactions of the Chinese Society of Agricultural Engineering, vol. 34, no. 7, pp. 66–74, 2018. View at: Google Scholar
Copyright
Copyright © 2018 Haitao Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.