Next Article in Journal
Advanced Autonomous Machines and Design Developments
Previous Article in Journal
Tracking and Counting of Tomato at Different Growth Period Using an Improving YOLO-Deepsort Network for Inspection Robot
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Kriging Model to Gear Wear Calculation under Mixed Elastohydrodynamic Lubrication

School of Mechanical Engineering and Automation, Northeastern University, Shenyang 110819, China
*
Author to whom correspondence should be addressed.
Machines 2022, 10(6), 490; https://doi.org/10.3390/machines10060490
Submission received: 17 May 2022 / Revised: 14 June 2022 / Accepted: 17 June 2022 / Published: 19 June 2022
(This article belongs to the Section Friction and Tribology)

Abstract

:
The calculation of tooth wear under mixed elastohydrodynamic lubrication is very complex and requires consideration of many conditions such as load distribution in the tooth meshing zone, micro-convex elastoplastic deformation and tooth surface temperature. The accurate calculation of tooth wear requires a lot of time and effort. In order to calculate tooth face wear under mixed elastomeric flow lubrication quickly and accurately, a new proxy model of tooth face wear is developed using the Kriging method. The pressure distribution required for the wear calculation was obtained utilizing the modified Reynolds equation and ZMC elasto-plastic model. The numerical calculation model of gear wear was derived using the modified Archard wear model. The Kriging model was used to construct a proxy model between gear parameters and tooth wear, and the degree of approximation and goodness of fit of the Kriging model were investigated. The results are as follows. The wear depth at each position is different, the smallest at the pitch, the largest near the tooth root, and the pinion has a larger wear depth than the gear. The Kriging model is highly efficient and accurate in its computation and overcomes the shortage of excessive time spent on the calculation of numerical calculation models.

1. Introduction

Gears are widely employed in aircraft, rail transportation, marine equipment, energy machinery, and other areas as an important basic mechanical component. Tooth surface wear is a typical failure mode in heavy-duty gear systems. Tooth surface wear increases transmission error and diminishes transmission accuracy [1,2], produces vibration and noise [3,4,5], affects transmission efficiency and smoothness, and speeds up gear failure [6]. The frictional wear mechanism of heavy-duty gears and the wear calculation problem under mixed elastohydrodynamic lubrication are difficulties that must be addressed in the design, manufacture, and application of novel gears. Holm, Burwell, and Archard et al. produced consecutive adhesive wear calculation formulae, the most extensively utilized of which being Archard’s generalized wear calculation equation [7]. Flodin [8] used the single-point observation approach to compute the contact pressure and sliding velocity of gears. Lundvall et al. [9,10] applied the finite element approach to compute tooth contact pressure and sliding distance for spur gears, but the sliding distance was calculated inefficiently. Park introduced a semi-analytic contact algorithm [11] and a surface interpolation approach [12], both of which enhanced the efficiency of gear wear estimations greatly. Heavy-duty transmission systems usually operate under mixed or boundary lubrication. Johnson et al. [13] initially suggested the notion of load sharing, where the working load is thought to be carried by the lubricant film between the surface and the micro-convex under mixed elastohydrodynamic lubrication. Bearings, cams and gears often operate in a mixed lubrication condition and wear can damage the surface and affect the clearance. Sander and Allmaier et al. [14] studied the friction of a shaft diameter bearing when moving from fluid to mixed lubrication by a combination of experiments and simulations, and analyzed the effect of the limit of constant boundary coefficient and the flow factor. Sharma et al. [15] investigated the effect of wear on the mixed lubrication performance of tapered journal bearings and obtained the static and dynamic performance characteristics of the bearings under different external loads, wear depths and different taper angles. Zhang et al. [16] proposed a method to correct the rounding of the roller ends of cylindrical roller bearings using wear. The mixed lubrication model combines the average flow model with the Kogut and Etsion model for contact pressure on rough surfaces. The results show that the optimized rolls have a uniform pressure distribution and the sharp increase in contact pressure of the micro-convex disappears. Winkler et al. [17] developed a numerical simulation method for studying the wear evolution of thrust roller bearings under mixed elastohydrodynamic lubrication. During the simulation, a finite element-based three-dimensional EHL model couples the elastohydrodynamic lubrication and the microscopic surface topography. The Archard wear model, the Greenwood–Williamson load-sharing model, and the Sugimura surface adaption model are all used in the numerical simulation technique. The aforementioned models explore the wear behavior of bearings under four load situations and two distinct mineral oil types working settings. Zhu [18] analyzed the dynamics and wear characteristics of the cam and guide barrel friction pair of a gasoline engine by computer programming, and predicted the minimum oil film thickness and friction coefficient, taking into account the transient effect and mixed lubrication. It is also pointed out that the extrusion effect of the lubricant and the surface roughness of the cam and guide barrel have a significant influence on the frictional wear performance. Wu et al. [19] developed a calculation equation for gear wear that takes into account gear dynamics and elastohydrodynamic lubrication. Considering the influence of surface roughness on elastohydrodynamic lubrication, Patir and Cheng [20,21] modified the Reynolds equation proposed by Dowson [22]. Majumdar and Hamrock [23] explored the fluctuation law of the load and the oil film thickness carried by the micro-convex contact with surface roughness using the modified Reynolds equation and the rough surface contact model extended by Greenwood and Tripp [24]. Zhao et al. [25] introduced a complete rough surface contact model that accommodates elastic, plastic, and elastoplastic deformation. The model is adaptable to a large variety of loads and surface roughness values. Masjedi and Khonsari derived a numerical model for oil film thickness and micro-convex contact ratio based on Zhao’s elastic–plastic rough surface contact model [26]. They then combined the aforementioned numerical model and the theory of thermal desorption of lubricating layer to develop a computer model for the interfacial friction coefficient and wear rate under mixed lubrication [27], which they used to calculate steady-state adhesive wear of spur gear teeth [28].
In mixed elastohydrodynamic lubrication, gear wear is a progressive and cumulative process with a long time horizon to reach the limit value, and the numerical simulation step size offered by the method must not be too big to maintain acceptable accuracy. Both of these considerations make calculating the wear volume time-consuming. When a large sample size is necessary to develop the relationship between gear parameters and wear, the numerical simulation algorithm is no longer appropriate. Approximate models can be applied instead of numerical wear models to tackle the aforesaid problem. Response surface approach [29], adaptive neural network method [30], Kriging method [31], and other approximation methods are extensively utilized. The response surface approach generates a fixed form of function, which is ineffective for dealing with complicated situations. The adaptive neural network’s resilience is poor. The Kriging model, on the other hand, performs better when dealing with situations involving significant nonlinearity and abrupt local response, and it does not require a huge number of design sample points. Laurenceau et al. [32] compared the Kriging model to the response surface approach in aviation dynamics design. The Kriging model was established by Feng et al. [33] to calculate the reliability of the main shaft vibration, and the study’s findings revealed that the Kriging operation is more efficient and that the method is suitable for the reliable calculation of substantially nonlinear systems. Yang et al. [34] established the Kriging model by Latin Hypercube Sampling to solve the problem of multiple responses and implicit functions for gear trim optimization. The calculation of gear wear under mixed elastohydrodynamic lubrication involves a large number of parameters with complex relationships between them. The calculation of gear wear requires the solution of a nonlinear system of equations, which is very complex and time-consuming. Especially when a large number of wear calculations are required, such as the study of gear wear mechanism and gear transmission reliability, traditional numerical calculations cannot meet the requirements. In order to solve the above problems, this paper establishes a Kriging model for gear wear calculation under mixed elastohydrodynamic lubrication by taking advantage of the Kriging model to obtain a proxy model with high prediction accuracy and fast solution speed based on a small amount of numerical simulation results. The model is also validated.
The remainder of the paper is laid out as follows. A numerical model of tooth wear under mixed elastohydrodynamic lubrication is built in Section 2. Section 3 describes some basic theories of the Kriging method. In Section 4, the Kriging method is used to create a proxy model for tooth surface wear under mixed elastomeric lubrication, and the computational results are discussed and analyzed. Finally, Section 5 covers the main conclusions.

2. Gear Wear Model under Mixed Elastohydrodynamic Lubrication

2.1. Mixed Elastohydrodynamic Lubrication Model

In mixed elastohydrodynamic lubrication, the contact surface is a rough surface, and the Reynolds equation should be adjusted to consider surface roughness effects. The surface load is borne by the micro-convex and the oil film in mixed elastohydrodynamic lubrication issue. The surface micro-convex may be deformed elastically, elasto-plastically, or plastically when the contact surface is deformed elastically under the load. As a result, the supplementary equation utilized to express the connection between the micro-convex contact pressure and the interface gap must be solved simultaneously.
Since the tooth surface load is supplied by the oil film and the micro-convex, the sum of the oil film pressure p h and the micro-convex contact pressure p a equals the total pressure p , i.e.,
p = p h + p a

2.1.1. Governing Equations

In line-contact elastohydrodynamic lubrication, Patir and Cheng amended the Newtonian fluid steady-state Reynolds equation to include roughness effects. The equation is as follows [22]:
x ( x ρ h 3 12 μ p h x ) = u ( ρ h T ) x
The variables included in the above equation are the oil film thickness h, the lubricant density   ρ , the lubricant viscosity μ, the rolling speed u, the pressure-flow coefficient x , and the average clearance between the contact surfaces h T . Equation (2) considers the influence of fluid compressibility. The first-order ordinary differential equation of the Equation (2) can be represented as follows:
d p h d x = 12 μ u ρ h T k r x ρ h 3
where k r is a constant that needs to be obtained. Substituting x corresponding the isothermal surface in reference [20] to Equation (2), we obtain:
d p h d x = 12 μ u ( h T ρ 1 k r ) h 3 ( 1 0.9 e 0.56 h / σ ) 1
where the surface roughness σ is defined as the standard deviation of the height of the surface micro-convex. The heights of the micro-convex on the surface obey a Gaussian distribution and the expression for h T is as follows [26]:
h T = 0.5 h [ 1 + e r f ( h 2 σ ) ] + σ 2 π e x p ( h 2 2 σ 2 )
The density–pressure function and the viscosity–pressure function can be described as follows [26]:
ρ ρ 0 = 1 + 0.6 p h 1 + 1.7 p h
μ μ 0 = exp { ( ln μ 0 + 9.67 ) [ 1 + ( 1 + 5.1 × 10 9 p h ) Z ] }
where ρ 0 is the lubricant density at standard atmospheric pressure, μ 0 is the lubricant viscosity at standard atmospheric pressure and Z is the viscosity–pressure index.
Define the following dimensionless covariates:
X = x b ,   H = h R ,   μ ¯ = μ μ 0 ,   ρ ¯ = ρ ρ 0 , P h = 4 R p h E e f f b = 1 E e f f 2 π W p h ,   U = μ 0 μ E e f f R ,   W = w E e f f R ,   G = α E e f f ,   σ ¯ = σ R
The parameters above are the Hertzian contact half-width b, the effective modulus of elasticity E e f f ( E e f f = 2 E e q ), the equivalent radius R, and the load w. In dimensionless form, the Reynolds equation and load balance equation may be expressed as:
d P h d X = 48 μ ¯ U ( H T ρ ¯ 1 K r ) H 3 ( 1 0.9 e 0.56 H / σ ¯ ) 1
π 2 = X m i n X e n d P h ( X ) d X + X m i n X e n d P a ( X ) d X
where X m i n and X e n d are the entrance and exit locations. Considering the elastic deformation and contact deformation of the surface, the expression for the lubricant film is [26]:
h ( x ) = h 0 + x 2 2 R 2 π E e f f x m i n x e n d p ln ( x s ) 2 d s
where h 0 is a constant which is related to the total pressure p. The dimensionless form is given by:
H ( X ) = H 00 + 4 W π [ X 2 1 π X m i n X e n d P ln ( X S ) 2 d S ]

2.1.2. Rough Surface Pressure Distribution Function

The rough surface contact pressure is the contact pressure between micro-convexes, which must account for both elastic, elastic–plastic, and completely plastic deformation of the micro-convex contact. Zhao’s elastic–plastic rough surface contact model [25] is applied to derive a pressure distribution function on the rough surface:
β ¯ = β R ,   n ¯ = n R 2 ,   V = h d E e f f
P a = 4 R E e f f b p a = 2 3 n ¯ β ¯ 0.5 σ ¯ 1.5 W 0.5 ( σ ¯ σ ¯ s ) I 1 I 2 ( z * I 1 ) 1.5 e 0.5 ( σ ¯ σ ¯ s z * ) 2 d z * + 2 π V n ¯ β ¯ σ ¯ W 0.5 ( σ ¯ σ ¯ s ) I 3 ( z * I 1 ) e 0.5 ( σ ¯ σ ¯ s z * ) 2 d z * + π V n ¯ β ¯ σ ¯ W 0.5 ( σ ¯ σ ¯ s ) I 2 I 3 ( z * I 1 ) e 0.5 ( σ ¯ σ ¯ s z * ) 2 × [ 1 0.6 ln w ¯ 2 ln ( z * I 1 ) ln w ¯ 2 ln w ¯ 1 ] × { 1 2 [ ( z * I 1 ) w ¯ 1 w ¯ 2 w ¯ 1 ] 3 + 3 [ ( z * I 1 ) w ¯ 1 w ¯ 2 w ¯ 1 ] 2 } d z *
where
I 1 = H y ¯ s σ ¯ ,   I 2 = H y ¯ s + w ¯ 1 σ ¯ ,   I 3 = H y ¯ s + w ¯ 2 σ ¯ ,   σ ¯ s = σ s R ,   y ¯ s = y s R ,   w ¯ 1 = ( 0.6 π V ) 2 β ¯
Furthermore, w ¯ 2 = 54 w ¯ 1 . Since n ¯ β ¯ σ ¯ = n β σ , either of these input parameters can be represented by the other parameters. The micro-convex density can be described by the following function:
n ¯ = 0.05 / ( β ¯ σ ¯ )

2.2. Geometry and Kinematics

The tooth profile of a spur gear is involute. The gear meshing relationship is shown in Figure 1.
The curvature radius of the pinion is R 1 and the curvature radius of the gear is R 2 . They can be obtained by the radii of the pinion and gear pitch circle R b 1 , R b 2 :
{ R 1 = R b 1 sin α + y R 2 = R b 2 sin α y
where y is the distance from the current meshing point to the pitch on the meshing line. The linear velocities at the meshing point are:
{ u 1 = π n 1 R 1 / 30 u 2 = π n 2 R 2 / 30
The linear velocities of the pinion and gear are u 1 , u 2 . The rotational speeds are n 1 , n 2 . The sliding distances of the pinion and gear are s 1 , s 2 which can be estimated according to the following formula:
{ s 1 = 2 a H | u 1 u 2 | u 1 s 2 = 2 a H | u 1 u 2 | u 2
where a H is the contact half-width and can be obtained by the following function:
a H = 4 F · R π · E e q

2.3. Wear Model

A model for calculating the wear rate under mixed lubrication is derived on the basis of the Archard wear model [27]:
Ω l u b = K Ψ ( L a 100 ) F u s v
where Ω l u b is the volume wear rate under the lubricated condition, L a is the surface micro-convex contact ratio, F is the contact load, u s is the relative sliding velocity, Ψ is the oil film deficit coefficient, and the expression is as follows:
Ψ = 1 e x p [ X u s t 0 e x p ( E a R g T s ) ]
where X is the diameter of the lubricant molecule, t 0 is the basic time of molecular vibration in the adsorption state, E a is the heat of adsorption of the lubricant at the interface, R g is the molar gas constant, and T s is the interface flash temperature. Based on the flash temperature theory of Tian and Kennedy [35], the formula of the tooth surface flash temperature is:
T S = T 0 + Δ T
Δ T = 2 a H q π ( K 1 1 + P e 1 + K 2 1 + P e 2 )
q = f c u s p ( L a / 100 ) + u s Λ l i m p ( 1 L a / 100 )
where K 1 , K 2 are the thermal conductivities of the two contact surfaces, q is the average heat generated at the contact point, P e 1 , P e 2 are the Peclet coefficients of the two contact surfaces, Λ l i m is the ultimate shear stress coefficient, and f c is the micro-convex friction coefficient.
Noting that p = F / A and Ω l u b = A h , A is the contact area, we can obtain the rate of wear depth to sliding distance as shown in the following equation:
h u s = K Ψ ( L a 100 ) p v
Therefore, the wear depth of the gear after n 0 revolutions can be expressed as:
h = n = 1 n 0 K Ψ ( L a 100 ) p v s
The cumulative wear of gears is minor when the number of engagements is low, and the geometry and pressure distribution of the tooth surface can be estimated to remain constant. The contact parameters must be updated when the number of meshes is large enough and the cumulative wear reaches a certain level.

2.4. Numerical Procedure

The Newton–Raphson algorithm is used to solve for the pressure and oil film thickness after the Equations (9), (10), (12) and (14) are discretized using the finite difference approach. The boundary conditions are X m i n = 4 at the inlet and the outlet location at several nodes after X = 1 . At each node, the film thickness and micro-convex contact pressure are determined. After forming the Jacobi matrix, the film pressure, K r , and H 00 at each node are obtained by solving the system of equations. The new value of the total pressure is used to update the film thickness, and iterations are repeated until the results converge. The contact pressure, micro-convex contact ratio and gear sliding distance, as well as other parameters, are imported into the modified Archard wear model for the calculation of the wear depth. Each wear cycle includes 1000 revolutions of the pinion.
It should be noted that the steady-state infinite line contact elastohydrodynamic lubrication model is established in this paper. In involute spur gear elastohydrodynamic lubrication, the calculation of the finite line contact model is very close to that of the infinite line contact model, thus the infinite line contact model is chosen to simplify the calculation in this paper. In the 2D case, we may expand the minimum lubricant gap, so the wear will be larger outside the y = 0 plane. Since the load and contact radius of curvature varies with time, leading to a damping effect of the lubricant during the film thickness change, the calculated value of the film thickness may be slightly lower in this model because the time-varying effect is not considered, resulting in a higher wear depth of the gear. In addition, the model ignores the effects of non-Newtonian rheological properties of the fluid and thermal effects on the film thickness.
The specific calculation process is shown in Figure 2.

3. Kriging Model

The Kriging model is a semi-parametric difference model that does not need the state function to be specified in any particular manner. Compared to numerical calculations, the model of response and model parameters built by the Kriging method can quickly and accurately predict the response values when given a sample [36]. The Kriging model is expressed as:
G ( x ) = h = 1 p β h g h ( x ) + z ( x ) = g T ( x ) β + z ( x )
where g T ( x ) β is the regression model, β is the regression coefficient and g(x) is the polynomial function of the random variable x. z(x) is the stochastic process function, reflecting the approximation of the local deviation, which has a mean of 0 and variance σ 2 , and its covariance matrix is
Cov [ z ( x i ) , z ( x j ) ] = σ 2 R ( x i , x j ; θ )
where R ( x i , x j ; θ ) is the correlation function which is usually expressed by the Gaussian correlation equation:
R ( x i , x j ; θ ) = m = 1 M e x p [ θ m ( x i m x j m ) 2 ]
where x i m is the mth component of the sample point x i .
Given a set of training samples S = [ x 1 , x 2 , · · · , x N ] , the response values are calculated according to the corresponding real functional functions and expressed in vector form as Y = [ g ( x 1 ) , g ( x 2 ) , · · · , g ( x N ) ] T , the unknown parameters β , σ 2 can be estimated by maximum likelihood estimation:
β ^ = ( G T R 1 G ) 1 G T R 1 Y
σ ^ 2 = 1 N ( Y G T β ) T R 1 ( Y G T β )
G = [ g ( x 1 ) , g ( x 2 ) , · · · , g ( x N ) ] T
Applying the Kriging model, the predicted response value of the point x to be measured is obtained as:
G ^ ( x ) = g ( x ) β ^ + r ( x ) T γ
σ G ( x ) 2 = σ ^ 2 [ 1 + u T ( x ) ( G T R 1 G ) 1 u ( x ) r T ( x ) R 1 r ( x ) ]
R = ( R ( x i , x j ; θ ) ) N × N
γ = R 1 ( Y G β ^ )
r ( x ) = [ R ( x 1 , x ; θ ) , · · · , R ( x N , x ; θ ) ] T
u ( x ) = G T R 1 r ( x ) g ( x )
The training random sample points, denoted by S, are generated by the Latin Hypercube Sampling method, and their actual functional response values, denoted by Y, need to be calculated by calculation. Using the training sample points S and the actual functional response value Y, the Kriging model is built. The Kriging model’s regression model is the main polynomial, with a Gaussian correlation function as the correlation function. The flowchart for constructing the gear wear Kriging model under mixed elastohydrodynamic lubrication is shown in Figure 3. The input values of the Kriging model are W, U, G, σ ¯ , β ¯ , V and wear time t. By taking the time parameter into account in the Kriging model, the wear trend can be predicted quickly and accurately. The output value of the model is the wear depth of gears.

4. Results and Discussion

The parameters of gear dimensions and material properties required to calculate the wear of the gear are shown in Table 1. The lubricant is SAE 30.
The torque, speed, surface hardness, surface roughness and wear coefficient are important factors affecting gear wear. Therefore they are used as random variables, and each random variable obeys normal distribution. The mean and standard deviation of the random variables are shown in Table 2.
The variation of the central film pressure along the line of action is shown in Figure 4. The central film pressure in the single tooth-meshing zone is greater than that in the double tooth-meshing zone. The central film pressure is also higher from the beginning of engagement because of the small roll suction speed and equivalent radius of curvature, although the load is smaller. The central film pressure in the single tooth-meshing zone decreases along the line of action, which is mainly caused by the change of the roll suction speed and equivalent radius of curvature.
The curve of the central film thickness along the line of action is shown in Figure 5. The film thickness is thinner at the beginning of engagement although the load is smaller due to the equivalent radius of curvature and the roll suction speed. The central film thickness has an increasing trend along the line of action, which is attributed to the increase in the roll suction speed. However, it has an abrupt drop in the middle part, which is due to a sudden change in load during the transition between the double- and single-tooth meshing of gears.
The gears’ wear depth after 1.5 × 10 5 revolutions of gear operation is shown in Figure 6. The pitch wear depth is the smallest, the root wear depth exceeds the top wear depth, and the pinion wear depth exceeds the gear wear depth. Compared with Ref. [8], the wear behavior of gears under mixed elastohydrodynamic lubrication is similar to that of dry friction. However, tooth wear under mixed lubrication is very much smaller than under dry friction for the same tooth profile parameters and operating parameters. This indicates that the rational use of lubricants can greatly reduce the reduction of tooth wear.
A Kriging model was developed using 100 samples. The variation curves of the maximum wear of gears are shown in Figure 7a,b. The accuracy of the model is an analysis of the progressive wear process of both the precise model and the proxy model where the random input values are deterministic and equal to the mean. As can be seen from Figure 7a,b, the wear depth increases cumulatively with time and the gear wear depth simulated by the Kriging model is very close to the actual wear depth. The accuracy of the Kriging model in approximating the actual wear depth is very high.
The goodness of fit refers to taking a certain number of samples outside the training samples to form a test set when the time is deterministic, and then the actual values of gear wear are compared with the predicted values of the Kriging model on the test set. Thirty sets of data were tested in the aforementioned Kriging model. The actual values of the samples were compared with the predicted values of the Kriging model, and the goodness of fit of the Kriging model was 0.9974. The actual wear depth and predicted wear depth of wear using the Kriging model are shown in Figure 8.
The wear depth from the Kriging model is compared to the results of another pair of gears derived from the numerical calculation to validate the correctness and generality of the Kriging model. The parameters are shown in Table 3.
It can be found that the results after 1 × 10 5 revolutions from the Kriging model correlate quite well with those from the calculation of the pinion wear depth for the mixed elastohydrodynamic lubrication, particularly near the pitch in Figure 9. In comparison to numerical calculation, it takes only a few minutes to obtain the wear depth of gears by the Kriging model, and the model’s accuracy in approximating real wear is quite high.
An important application of the Kriging model is to solve the reliability problem. In reliability studies, the time it takes to set up the original sample accounts for a large percentage of the total, usually about 99 percent or more. For the reliability problem of gear wear under mixed elastic flow lubrication, the Monte Carlo method is a widely accepted method, but it requires a large number of samples. If this method is used for reliability analysis of the implicit response function of gear wear, a large number of calculations are difficult to complete in a short time due to the low efficiency of a single calculation. Assuming that it takes 90 min to complete a gear wear calculation under mixed elastohydrodynamic lubrication and the study of gear wear reliability requires a sample space with an original sample size of 1 × 10 5 , the number of numerical calculations and the time required for the Monte Carlo method and the Kriging model is shown in Table 4. The time of the direct MCS method is unacceptable. It can be seen that the established Kriging model for gear wear with mixed elastohydrodynamic lubrication can accurately predict the approximate values of sample points, reduce the number of numerical calculations, and save a lot of time cost for the further study of gear wear.

5. Conclusions

In this paper, a model for computing gear wear under mixed elastohydrodynamic lubrication is constructed using the modified Reynolds equation, surface deformation model, elastoplastic micro-contact model, and modified Archard wear model. An approximate proxy model is created using the Kriging method to establish the direct relationship between gear parameters and wear, providing a novel way to predict gear wear. We have come to the following conclusions:
  • Under mixed elastohydrodynamic lubrication, the tooth surface load is supplied by the oil film and the micro-convex of the tooth surface, and tooth surface wear occurs in the micro-convex contact;
  • The smallest wear occurs at the pitch, the wear at the root is more than the wear at the top, and the pinion wear is greater than the gear wear;
  • The Kriging model may replace the numerical wear model with a limited number of samples, and the fit is excellent. This gear wear model may be used to increase the accuracy and efficiency of calculations in mixed elastohydrodynamic lubrication and to forecast unknown wear.

Author Contributions

Conceptualization, Z.S. and K.D.; methodology, K.D.; validation, K.D.; data curation, K.D.; writing—original draft preparation, K.D.; writing—review and editing, K.D. and Z.S.; supervision, Z.S.; project administration, Z.S.; funding acquisition, Z.S. 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, (grant no. 51775097).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ding, H. Dynamic Wear Models for Gear Systems. Ph.D. Thesis, The Ohio State University, Columbus, OH, USA, 2007. [Google Scholar]
  2. Atanasiu, V.; Oprisan, C.; Leohchi, D. The effect of tooth wear on the dynamic transmission error of helical gears with smaller number of pinion teeth. Appl. Mech. Mater. 2014, 657, 649–653. [Google Scholar] [CrossRef]
  3. Kuang, J.H.; Lin, A.D. The effect of tooth wear on the vibration spectrum of a spur gear pair. J. Vib. Acoust. 2001, 123, 311–317. [Google Scholar] [CrossRef]
  4. Choy, F.K.; Polyshchuk, V.; Zakrajsek, J.J.; Handschuh, R.F.; Townsend, D.P. Analysis of the effects of surface pitting and wear on the vibration of a gear transmission system. Tribol. Int. 1996, 29, 77–83. [Google Scholar] [CrossRef]
  5. Mackaldener, M.; Flodin, A.; Andersson, S. GDN-5 Robust noise characteristics of gears due to their application, manufacturing errors and wear(Gear dynamics and noise). In Proceedings of the JSME International Conference on Motion and Power Transmissions I(01-202), Fukuoka, Japan, 15 November 2001; pp. 21–26. [Google Scholar]
  6. Chen, Y.; Matubara, M. GSD-2 Effect of automatic transmission fluid on pitting fatigue strength of carburized gears (Gear Strength and durability). In Proceedings of the JSME International Conference on Motion and Power Transmissions I(01-202), Fukuoka, Japan, 15 November 2001; pp. 151–155. [Google Scholar]
  7. Archard, J.F. Contact and Rubbing of Flat Surfaces. J. Appl. Phys. 1953, 24, 981–988. [Google Scholar] [CrossRef]
  8. Flodin, A.; Anderson, S. Simulation of mild wear in spur gears. Wear 1997, 207, 16–23. [Google Scholar] [CrossRef]
  9. Lundvall, O. Simulation of Contact, Friction, and Wear in Gears-a Finite Element Approach. Ph.D. Thesis, Linköping University, Linköping, Sweden, 2000. [Google Scholar]
  10. Lundvall, O.; Klarbring, A. Simulation of wear by use of a nonsmooth newton method-a spur gear application. Wear 2003, 254, 1216–1232. [Google Scholar]
  11. Park, D.; Kolivand, M.; Kahraman, A. Prediction of surface wear of hypoid gears using a semi-analytical contact model. Mech. Mach. Theory 2012, 52, 180–194. [Google Scholar] [CrossRef]
  12. Park, D.; Kolivand, M.; Kahraman, A. An approximate method to predict surface wear of hypoid gears using surface interpolation. Mech. Mach. Theory 2014, 71, 64–78. [Google Scholar] [CrossRef]
  13. Johnson, K.L.; Greenwood, J.A.; Poon, S.Y. A simple theory of asperity contact in elastohydro-dynamic lubrication. Wear 1972, 19, 91–108. [Google Scholar] [CrossRef]
  14. Sander, D.E.; Allmaier, H.; Priebsch, H.H.; Witt, M.; Skiadas, A. Simulation of journal bearing friction in severe mixed lubrication—Validation and effect of surface smoothing due to running-in. Tribol. Int. 2016, 96, 173–183. [Google Scholar] [CrossRef] [Green Version]
  15. Sharma, S.C.; Phalle, V.M.; Jain, S.C. Influence of wear on the performance of a multirecess conical hybrid journal bearing compensated with orifice restrictor. Tribol. Int. 2011, 44, 1754–1764. [Google Scholar] [CrossRef]
  16. Zhang, Y.; Cao, H.; Kovalev, A.; Meng, Y. Numerical running-in method for modifying cylindrical roller profile under mixed lubrication of finite line contacts. J. Tribol. 2019, 141, 041401. [Google Scholar] [CrossRef]
  17. Winkler, A.; Marian, M.; Tremmel, S.; Wartzack, S. Numerical Modeling of Wear in a Thrust Roller Bearing under Mixed Elastohydrodynamic Lubrication. Lubricants 2020, 8, 58. [Google Scholar] [CrossRef]
  18. Zhu, G. Valve Trains—Design Studies, Wider Aspects and Future Developments. Eng. Tribol. 1993, 26, 183–211. [Google Scholar]
  19. Wu, S.; Cheng, H.S. Sliding wear calculation in spur gears. J. Tribol. 1993, 115, 493–500. [Google Scholar] [CrossRef]
  20. Patir, N.; Cheng, H.S. An average flow model for determining effects of three-dimensional roughness on partial hydrodynamic lubrication. J. Lubr. Technol. 1978, 100, 12–17. [Google Scholar] [CrossRef]
  21. Patir, N.; Cheng, H.S. Application of average flow model to lubrication between rough sliding surfaces. J. Tribol. 1979, 101, 220–230. [Google Scholar] [CrossRef]
  22. Dowson, D. A generalized Reynolds equation for fluid-film lubrication. Int. J. Mech. Sci. 1962, 4, 159–170. [Google Scholar] [CrossRef]
  23. Majumdar, B.C.; Hamrock, B.J. Effect of surface roughness on elastohydrodynamic line contact. J. Lubr. Technol. 1982, 104, 401–409. [Google Scholar] [CrossRef]
  24. Greenwood, J.A.; Tripp, J.H. The contact of two nominally flat rough surfaces. P I Mech. Eng. 1970, 185, 625–663. [Google Scholar] [CrossRef]
  25. Zhao, Y.; Maietta, D.M.; Chang, L. An Asperity microcontact model incorporating the transition from elastic deformation to fully plastic flow. J. Tribol. 2000, 122, 86–93. [Google Scholar] [CrossRef]
  26. Masjedi, M.; Khonsari, M.M. Film thickness and asperity load formulas for line-contact elastohydrodynamic lubrication with provision for surface roughness. J. Tribol. 2012, 134, 011503. [Google Scholar] [CrossRef]
  27. Masjedi, M.; Khonsari, M.M. An engineering approach for rapid evaluation of traction coefficient and wear in mixed EHL. Tribol. Int. 2015, 92, 184–190. [Google Scholar] [CrossRef]
  28. Masjedi, M.; Khonsari, M.M. On the prediction of steady-state wear rate in spur gears. Wear 2015, 342–343, 234–243. [Google Scholar] [CrossRef]
  29. Cheng, J.; Li, Q.S. Application of the response surface methods to solve inverse reliability problems with implicit response functions. Comput. Mech. 2009, 43, 451–459. [Google Scholar] [CrossRef]
  30. Papadrakakis, M.; Lagaros, N.D. Reliability-based structural optimization using neural networks and Monte Carlo simulation. Comput. Method Appl. M 2002, 191, 3491–3507. [Google Scholar] [CrossRef]
  31. Zhu, Z.F.; Du, X.P. Reliability analysis with Monte Carlo simulation and dependent kriging predictions. J. Mech. Des. 2016, 138, 121403. [Google Scholar] [CrossRef]
  32. Laurenceau, J.; Sagaut, P. Building efficient response surfaces of aerodynamic functions with kriging and cokriging. AIAA J. 2008, 46, 498–507. [Google Scholar] [CrossRef] [Green Version]
  33. Feng, J.L.; Sun, Z.L.; Zhao, J.; Zhang, J.; Liang, C.F. Vibration reliability analysis on a spindle system based on AK-MCS method. Chin. J. Vib. Shock 2019, 38, 135–140. [Google Scholar]
  34. Yang, L.; Tong, C.; Chen, C.; Guo, Q.P. Vibration reduction optimization of gear modification based on Kriging model and genetic algorithm. Chin. J. Aerosp. Power 2017, 32, 1412–1418. [Google Scholar]
  35. Tian, X. Maximum and Average Flash Temperatures in Sliding Contacts. J. Tribol. 1994, 116, 167–174. [Google Scholar] [CrossRef]
  36. Matheron, G. The intrinsic random functions and their applications. Adv. Appl. Probab. 1973, 5, 439–468. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Gear meshing relationship.
Figure 1. Gear meshing relationship.
Machines 10 00490 g001
Figure 2. Flowchart of the calculation of the gear wear under mixed elastohydrodynamic lubrication.
Figure 2. Flowchart of the calculation of the gear wear under mixed elastohydrodynamic lubrication.
Machines 10 00490 g002
Figure 3. Flowchart for building the Kriging model.
Figure 3. Flowchart for building the Kriging model.
Machines 10 00490 g003
Figure 4. The variation of the central film pressure.
Figure 4. The variation of the central film pressure.
Machines 10 00490 g004
Figure 5. The curve of the central film thickness.
Figure 5. The curve of the central film thickness.
Machines 10 00490 g005
Figure 6. Wear depth of the pinion and gear.
Figure 6. Wear depth of the pinion and gear.
Machines 10 00490 g006
Figure 7. Calculation accuracy of the Kriging model: (a) Pinion; (b) Gear.
Figure 7. Calculation accuracy of the Kriging model: (a) Pinion; (b) Gear.
Machines 10 00490 g007
Figure 8. Actual wear depth and prediction based on Kriging model.
Figure 8. Actual wear depth and prediction based on Kriging model.
Machines 10 00490 g008
Figure 9. Wear depth of the Kriging model compared to the results of the numerical calculation: (a) Pinion; (b) Gear.
Figure 9. Wear depth of the Kriging model compared to the results of the numerical calculation: (a) Pinion; (b) Gear.
Machines 10 00490 g009
Table 1. Parameters required for calculation.
Table 1. Parameters required for calculation.
ParametersPinionGear
Number   of   teeth   n 1 , 2 1624
Module m (mm)4.54.5
Pressure   angle   α (deg) 20 ° 20 °
Tooth width b (mm)1515
Center   distance   a w (mm)91.591.5
Modification coefficient0.86350.5
Modulus of elasticity E (Gpa)210210
Poisson s   ratio   ν 0.30.3
Table 2. Distribution information of variables.
Table 2. Distribution information of variables.
VariableDistributionMeanStandard Deviation
Torque T (N m)Normal30210
Speed v (rpm)Normal15030
Hardness hd (GPa)Normal2.350.2
Roughness   σ   ( μ m )Normal0.250.02
Wear coefficient KNormal 5 × 10 4 1 × 10 4
Table 3. Main parameters required for the Kriging model.
Table 3. Main parameters required for the Kriging model.
ParametersPinionGear
Number   of   teeth   n 1 , 2 3296
Module m (mm)55
Pressure   angle   α (deg) 20 ° 20 °
Tooth width b (mm)5050
Contact load F (KN)2525
Speed v (rpm)100100
Hardness hd (GPa)2.32.3
Roughness   σ   ( μ m )11
Wear coefficient K 5 × 10 4 5 × 10 4
Table 4. Consuming time of the MCS method and Kriging model.
Table 4. Consuming time of the MCS method and Kriging model.
MethodNumber of CalculationsTime
Kriging MCS100150 h
Direct MCS 1 × 10 5 17 years
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dong, K.; Sun, Z. Application of Kriging Model to Gear Wear Calculation under Mixed Elastohydrodynamic Lubrication. Machines 2022, 10, 490. https://doi.org/10.3390/machines10060490

AMA Style

Dong K, Sun Z. Application of Kriging Model to Gear Wear Calculation under Mixed Elastohydrodynamic Lubrication. Machines. 2022; 10(6):490. https://doi.org/10.3390/machines10060490

Chicago/Turabian Style

Dong, Kunpeng, and Zhili Sun. 2022. "Application of Kriging Model to Gear Wear Calculation under Mixed Elastohydrodynamic Lubrication" Machines 10, no. 6: 490. https://doi.org/10.3390/machines10060490

APA Style

Dong, K., & Sun, Z. (2022). Application of Kriging Model to Gear Wear Calculation under Mixed Elastohydrodynamic Lubrication. Machines, 10(6), 490. https://doi.org/10.3390/machines10060490

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