Dynamic Assessment of Vibration of Tooth Modification Gearbox Using Grey Bootstrap Method

The correlation analysis between gear modification and vibration characteristics of transmission system was difficult to quantify; a novel small sample vibration of gearbox prediction method based on grey system theory and bootstrap theory was presented. Themethod characterized vibration base feature of tooth modification gearbox by developing dynamic uncertainty, estimated true value, and systematic error measure, and these parameters could indirectly dynamically evaluate the effect of tooth modification. Themethod can evaluate the vibration signal of gearboxwith installation of no toothmodification gear and topologicalmodification gear, respectively, considering that 100% reliability is the constraints condition andminimumaverage uncertainty is the target value. Computer simulation and experiment results showed that vibration amplitude of gearbox was decreased partly due to topological toothmodification, and each value of average dynamic uncertainty, mean true value, and systematic errormeasure was smaller than the no toothmodification value.The study provided an important guide for toothmodification, dynamic performance optimization.


Introduction
With the development of high speed and heavy load gear transmission, more and more people pay attention to the problems of vibration and noise.It demands harshly the vibration in aviation, marine, and other areas of transmission system; how to assess and reduce transmission system vibration is a direction of current research.Modification technology as the effective means of improving the performance of gear is paid attention early; it can better improve meshing condition by tooth surface modification, effectively prevent edge contact, and reduce vibration and noise [1][2][3].
In the transmission process, the gearbox vibration has become the main incentive source of transmission system because of the gear mesh stiffness and error incentive, so the prediction of the dynamic characteristics of the gearbox vibration is extremely important [4][5][6].Kato [7] introduced the FEM/BEM method of predicting vibration and noise radiation, which was consistent with the experimental results.Zhou et al. [8] proposed formula-fitting method for predicting noise radiation, which can directly predict noise by the basic system parameters of the gearbox.Tuma [9] analyzed the incentive component and prediction method of the gearbox vibration noise and proposed the method of vibration isolation from the tooth geometric shapes, stiffness of the gearbox, and so on.However, the solution of gearbox vibration noise quantitative assessment needs to consume a lot of computing resources, time, and effort.Xia et al. [10] established a dynamic assessment model GBM(1, 1) for manufacturing errors using grey bootstrap method, in order to solve the problem of predicting manufacturing errors with poor information.Wang et al. [11] gained the predictions of calibration data of all error sources by a grey bootstrap fusion model, which can accurately predict measurements of the material Brinell hardness.Grey bootstrap method has been used in the aspects of dynamic assessment of bearing vibration [12], uncertain assessment of virtual instrument in small sample measuring [13], and dynamic predicting of the shake of radial sloshing error [14], which coincided with the actual measurements.But gearbox vibration dynamic predicting by grey theory and grey bootstrap has not been studied [15,16].

Shock and Vibration
This paper established grey bootstrap dynamic assessment model GBM(1, 1) [17] by grey predicting model GM(1, 1) [18,19] and bootstrap theory and dynamically assessed gearbox vibration in a small sample with poor information.Compared to no modification gear, this paper indirectly analyzed the impact of modification gear on the dynamic characteristics of the transmission system and demonstrated the accuracy of the dynamic assessment by conducted experimental investigation.

Dynamic Predicting Model GBM(1, 1)
2.1.Principle of GBM (1,1).Due to machine errors, heat treatment distortions, variation of cutting forces, and other unpredictable factors reduce gear quality and cause unfavorable displacement of tooth contact and increased transmission errors, resulting in edge contact, highly concentrated stresses, and vibration.Tooth modification is used to prevent increased levels of noise and vibration due to the deformation of the teeth.So it is important to evaluate the uncertainty in measurement of vibration of gearbox.The principle of grey bootstrap model GBM(1, 1) is shown in Figure 1, which is a method combining the bootstrap theory with grey predicting model GM(1, 1) to predict measurement uncertainties.

Building Bootstrap Resampling Samples of Vibration.
In the gearbox, gear and bearing are the factors causing vibrations.By replacing different modification gears, gearbox vibration will be much different.In the large-volume gear grinding process, the processed products in production process can be seen as an infinite population; at the time  replacing the gear of the gearbox, data series vector of gearbox vibration can be obtained and given by Y = { ()} ;  = 1, 2, . . ., where () is the original data of the vibration of the gearbox at the time .To estimate the uncertainty at the time  by means of  vibration data, which are close and before the time  (including the time ), the data series vector at the time  is constituted and is expressed by where  is the time  and  is the data number of subseries Y  of dynamic predicting.In bootstrap statistics, Y  is the initial sample and  is the capacity of sample.
According to the bootstrap theory, randomly extracting  times vibration data of gearbox with equal probability from Y  on sampling with replacement, a bootstrap sample with  vibration data is formed.Repeating  times in a row, the  bootstrap resampling samples can be obtained and numerically represented as where B  is the th bootstrap sample, and

Calculating Parameters of GBM.
According to the grey system theory, the first-order accumulated generating operation of B  is defined as The grey generated model can be described by the differential equation which is given by where  is time variable,  1 and  2 are the coefficients to be estimated, and in addition  1 is not equal to zero.Substitute the increment with the differential, so Suppose average generated series vector is In the initial condition   ( −  + 1) =   ( −  + 1), suppose grey differential equation least-squares solution is where the coefficients,  1 and  2 , are given by So at the th time the predicted value is as follows: b () = ŝ () − ŝ ( − 1) ;  =  + 1.
Series vector of the  data at the th time is as follows: The  data of Ŝ+1 are established in the histogram; the probability density distribution curve can be obtained.At the th time, let  data be divided into  groups by a certain interval;   is the class midvalue of th group; the frequency in the group is   .Then, in the time  + 1 the true predicted value can be calculated as a weighted average: 2.4.Output Assessment Indicators Value.The model characterizes vibration of gearbox by developing six parameters, such as estimated true value, interval dynamic uncertainty, average uncertainty, average true value, and system trend error measure.Assume the significance level is ; the confidence level can be given as follows: At the time  + 1 an estimated interval of true value of vibration at the confidence degree  is as follows: (15) where  /2 is the value of gearbox vibration   corresponding to a probability /2;  1−/2 is the value of gearbox vibration   corresponding to a probability 1−/2;   is the low boundary of the estimated interval;   is the upper boundary of the estimated interval.
When predicting, the system inputs the first data  = ; if there are ℎ vibration data outside the estimated interval [  ,   ], the predicted accuracy   is expressed as follows: The expanded uncertainty of gearbox vibration at time  is defined as where  is the estimated uncertainty at the time  under the  confidence level, and it is, namely, dynamic uncertainty.
If  = 100%, the maximum  will be received, but the larger  is, the farther [  ,   ] get away from the true value.Therefore, the average uncertainty is defined as In order to assess the size of the vibration value of the gearbox,  0 is defined as Here  0 is the mean true value of the gearbox vibration.In order to assess the impact of trends in terms of vibration,  0 is defined as where  0 is the range value of the estimated value of vibration and is, namely, the system error measure.

Gearbox Vibration Signals
Gearbox contains gears, bearings, shafts, and other components and during operation gear will produce meshing impact and cause vibration corresponding to gear mesh frequency.Vibration signals of gearbox are measured from the bearing housing, containing the meshing vibration response of the gear and also including vibration produced by the rolling bearing and the response of other incentives in gearbox which is expressed as where  gear () is the meshing vibration response of the gear under the normal condition;  bearing () is the meshing vibration response of the bearing; and   () is the vibration response of the random incentives.Namely, in the certain condition, the speed and load are determined, gear meshing vibration response can be expressed as gear mesh frequency  mesh , and its harmonic frequency in the form of series is expressed as where   is the meshing frequency component amplitude of order ; Φ  is phase;   () and   () are additional amplitude and phase modulation generated errors incentives.The mesh frequency  mesh can be expressed as where   is gear rotation frequency,  is gear tooth number, and  is gear speed (r/min).The torque on gearbox output shaft is controlled by magnetic loader according to its current size, the torque and speed sensor measures the torque and speed of the output shaft, and sound pressure meter records the noise.Voltage signal measured by acceleration sensor is processed by the data acquisition instrument.Internal structure of the gearbox is shown in Figure 3, consisting of a two-stage parallel shaft helical gear transmission; gearbox basic parameters are shown in Table 1.During the transmission process, because of the effect of the gear mesh stiffness and error excitation, the fluctuation of gear meshing force can cause the vibration of the system.The module of vibration acceleration vector sum in different directions on the box is defined as the vibration intensity  rms which is expressed as

Experimental Study
where ∑   , ∑   , and ∑   are the sum of acceleration vibration of each measuring point in different directions under the same work condition and   ,   , and   represent measuring points in different directions, respectively.Gearbox vibration is affected by the working condition; for example, different speed and torque have their own vibration effects, respectively.Suppose the speed of the motor is 300 r/min, the current of magnetic loader is 0.46 A, and the output torque is 125 NM; the test apparatus is shown in  2) and no modification driven gears which are, respectively, installed in the gearbox are detected.The amplitude of vibration, respectively, is 0.02037 g and 0.03056 g, and gearbox vibration amplitude of tooth modification compared to no modification has decreased.Considering, from a single piece of modification and no modification gear, tooth modification has effect on gearbox vibration, for a detailed assessment of its vibration reduction effects, the grey dynamic self-assessment method has been used.

Gearbox Vibration Dynamic Assessment.
Taking 30 sets of no modification gear and 30 sets of modification gear installed in the gearbox, respectively, the original data  of gearbox vibration is measured.The results of dynamic assessment are shown in Figures 7 and 8, where  is the index of the groups of gears to be tested.
By calculating, the minimum parameter of average uncertainty  and evaluation results are obtained, which is shown in Table 3.The gearbox vibration mean value  0mean of no tooth modification is large; the gearbox vibration mean value  0mean of modification gear is decreased.
Figures 7(a  20th set of data; this process is reflected by  0 .For example, the system error is large, where  0 equals 0.00388 g in Table 3.In Figure 8(b), the volatility has a slight increase from the 20th set of data, but, compared to Figure 7(b), the volatility is somewhat flat and the system error is smaller, where  0 equals 0.00253 g.All the above results show that GBM(1, 1) can dynamically assess the instantaneous impact of system error on the gearbox vibration.Figures 7(c) and 8(c) are the dynamic description of the uncertainty  of .The larger the random volatility of  is, the larger the value of  is.The uncertainty  reflects the degree of data difference between time  and time  + 1, which is independent of the change trend of  0 .It shows that GBM(1, 1) can dynamically separate systematic errors, but it does not regard the impact of systematic errors in the dynamic assessment.
Table 3 is an overall assessment of the gearbox vibration before and after gear modification, compared to the reduction amplitude of three parameters such as dynamic uncertainty, the average true value, and system error measure.The reduction amplitude of the system error measure  0 is more obvious, and it shows that the system error measure  0 is more sensitive to gearbox vibration, so it can be used as an effective parameter to judge the gear modification effect and help with the dynamic assessment of gear modification effect.

Conclusions
Based on grey prediction model GM(1, 1) and bootstrap sampling method, to establish grey bootstrap dynamic assessment model GBM(1, 1) of gearbox vibration with modification gear, it is a better solution to the problem of dynamic assessment of gear modification.Evaluation index system of six parameters has been established, such as estimated true value, interval dynamic uncertainty, average uncertainty, average true value, and system trend error measure.It is dynamically assessed with 100% reliability as constraints condition and minimum average uncertainty as the target value.According to the parameter of GBM(1, 1) assessment shows       that a reasonable tooth modification has significant damping effect on gearbox and the amplitude is decreased.However, the effect of modification is different in other conditions; it is difficult to determine the optimal tooth modification in a variety of conditions.Modification parameters must be based on the actual working conditions in order to verify and correct the gearbox vibration to a minimum.The analysis of this paper provides the reference to determining the parameters of gear modification and dynamic assessment for the gearbox vibration.

Figure 2 :
Figure 2: Schematic illustration of gearbox vibration test platform.

Figure 4 .
Acceleration sensor in the gearbox bearing cap outer radial direction extracts vibration acceleration signals; the signals are dealt with through the signal conditioning and sent to the data acquisition system; the sampling frequency of the vibration signal is 5000 Hz, without considering the impact of other errors in the measurement.The vibration signals of gearbox with 30 modification active gears (topographic modification parameters are shown in Table

Figures 5 (
Figures 5(a) and 6(a) are gearbox vibration time domain curves of no modification gear and modification gear, respectively.In Figure 5(a), the magnification factor is 10 in the measurement.Figures 5(b) and 6(b) are gearbox vibration frequency domain curves of no modification and modification gear, respectively.The main vibration frequencies of the gearbox are 115, 230, 345, and 460 Hz.The amplitude of vibration, respectively, is 0.02037 g and 0.03056 g, and gearbox vibration amplitude of tooth modification compared to no modification has decreased.Considering, from a single piece of modification and no modification gear, tooth modification has effect on gearbox vibration, for a detailed assessment of its vibration reduction effects, the grey dynamic self-assessment method has been used.
) and 8(a) are the dynamic description of random volatility of the original data  for the estimated interval [  ,   ] which reflected the trend of .[  ,   ] can completely envelop the volatility of , indicating that GBM(1, 1) can better describe the transient trend of gearbox vibration acceleration.Figures 7(b) and 8(b) are the dynamic description of estimated true value  0 for the trend change of .In Figure 7(b), the volatility is increasing significantly from the

Figure 5 :
Figure 5: Vibration signal detection of gearbox with no modification gear.

Figure 6 :
Figure 6: Vibration signal detection of gearbox with modification gear.
Original data  and estimated interval [  ,   ]

Figure 7 :
Figure 7: Dynamic evaluation for vibration acceleration of gearbox with no modification gear.
Original data  and estimated interval [  ,   ]

Figure 8 :
Figure 8: Dynamic evaluation for vibration acceleration of gearbox with modification gear.

Table 1 :
Parameters of the gearbox.
[20]Gearbox Vibration Testing.Load platform consists of the gearbox, magnetic loader, adjustable motor, torque sensor, acceleration sensor, and sound level meter, whose load platform structure diagram has been shown in Figure2[20].

Table 2 :
Parameters of driven gear 2.