Sage Journals: Discover world-class research

Abstract

The lubrication characteristics of a cylindrical gear transmission with a variable hyperbolic circular-arc-tooth-trace (VH-CATT) are important theoretical bases for the tribology design of the gear and its fatigue life prediction. Existing lubrication characteristics studies of it are based on simplified linear contact models, which are limited to one specific meshing position only. Therefore, the mathematical models for tooth contact analysis (TCA) and loaded tooth contact analysis (LTCA) are established to calculate most time-varying parameters in a meshing period, such as the kinematic parameters, spatial geometric parameters, and contact parameters, are calculated. Accordingly, a thermo-elastohydrodynamic lubrication (TEHL) model of an elliptical contact of this gear transmission is presented. Using all the proposed mathematical models, the TEHL results with smooth and rough tooth surfaces in a meshing period are analyzed. The obtained results reveal that these models can be used to predict the distribution of the film thickness, pressure, and temperature rise on the smooth and rough tooth surfaces of VH-CATT cylindrical gear transmission in each contact area of a meshing period. And it shows that the time-varying parameters have great influences on lubrication performance. Furthermore, the lubrication mechanism in a meshing period is investigated systematically in VH-CATT cylindrical gear transmission.

Keywords

VH-CATT cylindrical gear transmission TCA elliptical contacts time-varying parameters TEHL

Introduction

VH-CATT cylindrical gear is a new type of gear developed based on the Gleason spiral bevel gear, as shown in Figure 1. Compared with traditional spur and helical gears, VH-CATT cylindrical gears have a strong bearing capacity, higher transmission efficiency, lower installation precision requirements, and no axial force.^1,2 Because of its excellent performance, VH-CATT cylindrical gears have significant application prospects. In recent years, the exploration of VH-CATT cylindrical gear transmission progressed considerably. Tseng and Tsay³ and Wu et al.⁴ deduced the full tooth surface equation for VH-CATT cylindrical gears with transition tooth surface, and Zhao et al.⁵ established an accurate three-dimensional mathematical model. Then, Wu et al.⁶ studied the effect of machine tool errors on the tooth surface errors of this gear. Generally, poor lubrication is the main factor leading to frictional wear of gears. In a meshing period, most time-varying parameters, such as the gear geometry, speed, contact geometry, and load at each meshing moment, constantly vary. Concurrently, with the influences of pressure, temperature field, and tooth surface roughness in the lubricated contact area, the surface elastic deformation, lubricating oil viscosity and density will significantly change. Therefore, it is necessary to develop a TEHL model of VH-CATT cylindrical gears considering most time-varying parameters to study the lubricant behavior in the meshing process.

Figure 1.

Model of VH-CATT cylindrical gear.

In 1916, for the first time, Martin⁷ used a simplified model to apply the Reynolds equation for analyzing gear lubrication. With the maturity of gear technology and the elastohydrodynamic lubrication (EHL) theory as well as significant improvement in numerical algorithms and computer performance, the application to EHL theory has spanned a century and is still providing innovations. For instances, Xu and Kahraman⁸ and Liu et al.⁹ developed the EHL-line and TEHL-line contact model of a hypoid bevel gear and helical gear transmission, respectively. It should be noted that these studies mainly focused on theoretical conjugate tooth surfaces. In practical industrial applications, owing to the deviation caused by gear assembly and machining errors, elliptical contacts are formed on gear contact traces instead of theoretical line contacts. Therefore, a highly realistic elliptical contact was introduced into the EHL model of hypoid gears¹⁰ and the applicability of this model was subsequently expanded to non-Newtonian fluid.¹¹ To achieve further closeness to a real scenario, the present study examines the lubrication mechanism of VH-CATT cylindrical gears with elliptical contact.

Following the determination of the effects of tooth surface roughness and thermal effect factor on gear lubrication performance, many scholars have progressively conducted EHL model studies with these on various gears. After Bobach et al.¹² introduced the roughness into the TEHL of an involute spur gear, Xue et al.¹³ further considered load changes in the TEHL analysis of this gear. The lubrication performance of helical gears was studied by using a shear-thinning TEHL model incorporating ultimate shear stress.¹⁴ For planetary gear, the model of EHL was established, which consider the effects of roughness.¹⁵ After Pu et al.¹⁶ established the EHL model considering roughness of a spiral bevel gear, Zhang et al.¹⁷ predicted the coefficient of friction and tooth surface wear of this gear using TEHL technology. After establishing a TEHL model covering the full film, mixed, and boundary lubrication, Simon presented the methods to improve the lubrication performance of hypoid gear based on multi-objective optimization,¹⁸ optimization of manufacture parameters,¹⁹ and machine tool setting parameters.²⁰ Therefore, when analyzing the lubrication performances of VH-CATT cylindrical gears, the effects of the tooth surface roughness and the temperature rise cannot be ignored.

Most recently, numerous scholars have considered the time-varying parameters in the meshing process when analyzing the lubrication performance of gear transmission. For instance, the effect of load and speed were introduced into the TEHL model of a helical gear.²¹ In addition, Wang et al.²² presented a TEHL model of a spiral bevel gear, considering of the velocity vector and the surface roughness. Considering certain time-varying parameters along the meshing trace of a spiral bevel gear, the progressive mesh densification was introduced in the numerical calculation of a EHL model²³ and the behavior of dynamical wear,²⁴ friction, temperature, and contact fatigue²⁵ was subsequently discussed. Then, an EHL analysis of spiral bevel gears was conducted after considering the contact point migration due to the deformation under a load.²⁶ These studies serve as important references for the establishment of an accurate TEHL model of VH-CATT gear transmission.

The coupling of kinematic, spatial geometric, and contact parameters increases the complexity of lubrication analyses of VH-CATT cylindrical gears. However, existing studies on the EHL²⁷ and TEHL²⁸ of VH-CATT cylindrical gears are based on simplified linear contact models, which are limited to a specific meshing position. These studies ignore not only the time-varying parameters along the contact trace during gear meshing but also the tooth surface roughness. Therefore, this paper proposes a method for comprehensively calculating the time-varying parameters on a contact trace in a meshing cycle, which is developed using TCA and LTCA models. Subsequently, a TEHL model of VH-CATT cylindrical gears is established, and these time-varying parameters are extracted as the input parameters of the TEHL model. After the numerical calculation, the film thickness, pressure, and temperature rise distributions of a VH-CATT cylindrical gear pair under the effects of the time-varying parameters in a meshing period along the contact trace are analyzed in detail.

Contact analysis of VH-CATT cylindrical gears

To simulate the meshing of a gear pair under real working conditions, TCA and LTCA mathematical models based on the tooth surface equation of VH-CATT cylindrical gear transmission are established. The contact trace as well as time-varying parameters such as the main curvature, entrainment speed, slide-roll ratio, contact ellipse, load distribution, and maximum Hertzian contact pressure in the meshing process, obtained via TCA and LTCA in a meshing period are used as the input parameters for the TEHL model.

Tooth surface equation

VH-CATT cylindrical gears are machined using a double-edged milling cutter fixed on a rotating cutter head. The machining principle is shown in Figure 2.

Figure 2.

Principle of rotating cutter-head milling.

The tooth surface equation⁵ of the VH-CATT cylindrical gears is written as equation (1).

Where the upper (lower) symbols of ∓ and ± represent the convex (concave) tooth surface of the gears machined by the inner (outer) cutting edge; α is the pressure angle; θ is the cutter rotation angle, φ is the workpiece rotation angle; m is the gear module; R is the radius of the rotating cutter head (radius of the tooth trace), R₁ is the pitch radius.

{\begin{matrix} x_{0} = [- (R \mp \frac{π}{4} m \pm u_{0} \sin α) \cos θ + R + R_{1} φ] \\ \cos φ + (u_{0} \cos α - R_{1}) \sin φ, \\ y_{0} = [(R \mp \frac{π}{4} m \pm u_{0} \sin α) \cos θ - R - R_{1} φ] \\ \sin φ + (u_{0} \cos α - R_{1}) \cos φ, \\ z_{0} = (R \mp \frac{π}{4} m \pm u_{0} \sin α) \sin θ, \\ u_{0} = \mp \sin α \frac{\cos θ (R \mp \frac{π}{4} m) - (R_{1} φ + R)}{\cos θ} \end{matrix}

(1)

TCA model

As shown in Figure 3, the spatial coordinate system of the VH-CATT cylindrical gear pair is established based on the spatial configuration during meshing. o_fx_fy_fz_f is the fixed coordinate system; o_px_py_pz_p and o_gx_gy_gz_g are the workpiece dynamic coordinate systems fixed with driving and driven gears, respectively; o_mx_my_mz_m and o_nx_ny_nz_n are the intermediate coordinate systems in the coordinate transformation; l is the center distance; ϕ_p and ϕ_g are the meshing angles of the driving and driven gear, respectively.

Figure 3.

Coordinate systems used for movement of driving and driven tooth surfaces.

It is assumed that Σ_p (the concave tooth surface) in o_px_py_pz_p represents the family of surfaces r _p and that its unit normal vector is n _p. Similarly, Σ_g (the convex tooth surface) in o_gx_gy_gz_g denotes the family of surfaces r _g, and its unit normal vector is n _g. In o_fx_fy_fz_f, the vector and unit normal vector of Σ_p and Σ_g are denoted as ( r _fp, n _fp) and ( r _fg, n _fg), respectively.

Thus, the TCA mathematical model of the VH-CATT cylindrical gear pair in the meshing process is expressed as:

{\begin{matrix} r_{fp} (φ_{p}, θ_{p}, ϕ_{p}) - r_{fg} (φ_{g}, θ_{g}, ϕ_{g}) = 0 \\ n_{fp} (φ_{p}, θ_{p}, ϕ_{p}) - n_{fg} (φ_{g}, θ_{g}, ϕ_{g}) = 0 \end{matrix}

(2)

The transformation of Σ_p into o_fx_fy_fz_f is completed by:

{Σ_{p}} : {\begin{matrix} r_{fp} = M_{fp} r_{p} \\ n_{fp} = L_{fp} n_{p} \end{matrix}

(3)

Where M (4 × 4) is the transformation matrix of r _p from o_px_py_pz_p to o_fx_fy_fz_f, and L _fp (3 × 3) is the transformation matrix corresponding to n _p of Σ_p.

The transformation of Σ_g into o_fx_fy_fz_f is completed by

{Σ_{g}} : {\begin{matrix} r_{fg} = M_{fn} M_{nm} M_{mg} r_{g} \\ n_{fg} = L_{fn} L_{nm} L_{mg} n_{g} \end{matrix}

(4)

Where M _mg (4 × 4), M _nm (4 × 4), and M _fn (4 × 4) are the transformation matrices of r _g from o_gx_gy_gz_g to o_fx_fy_fz_f, and L _mg (3 × 3), L _nm (3 × 3), and L _fn (3 × 3) are the transformation matrices corresponding to n _g of Σ_g.

Based on the existence theorem of implicit function,²⁹ five independent equations with five unknown parameters are derived from equation (2), as follows:

{\begin{matrix} f_{i} [φ_{p} (ϕ_{p}), θ_{p} (ϕ_{p}), φ_{g} (ϕ_{p}), θ_{g} (ϕ_{p}), ϕ_{g} (ϕ_{p})] = 0 \\ f_{i} \in C \end{matrix}, (i = 1, 2, 3, 4, 5)

(5)

Taking the design parameters of the VH-CATT cylindrical gear pair listed in Table 1 as a case, all data reflecting the gear meshing status are solved through numerical iterations, the results as shown in Figure 4.

Table 1.

Design parameters.

Parameters	Numerical value
Module of gear pair (mm)	8
Number of driving gear teeth	21
Number of driven gear teeth	29
Pressure angle (°)	20
Tooth width (mm)	80
Radius of the tooth trace (mm)	500

Figure 4.

Contact trace in (a) coordinate system o_px_py_pz_p and (b) coordinate system o_gx_gy_gz_g.

The curvature of a VH-CATT cylindrical gear is an important spatial geometric parameter affecting its lubrication behavior. As shown in Figure 5, the principal orientations of Σ_p and Σ_g are t ^a_p, t ^b_p, t ^a_g, and t ^b_g, respectively. Correspondingly, the principal curvatures are K^a_p, K^b_p, K^a_g, and K^b_g.

Figure 5.

Curvature relationship of tooth surface.

Thus, the mathematical expression of the main curvature of the VH-CATT cylindrical gears can be derived as:

{\begin{matrix} K_{i}^{a} = K_{iA} + \sqrt{{K_{iA}}^{2} - K_{iG}} \\ K_{i}^{b} = K_{iA} - \sqrt{{K_{iA}}^{2} - K_{iG}} \end{matrix}, (i = p, g)

(6)

Where K_iG is the Gaussian curvature and K_iA is the average curvature of the tooth surface, which can be deduced by the surface theory.³⁰

Generally, velocity parameters have a major influence on the TEHL analysis of gear transmission. Therefore, the entrainment velocity, u _e, the absolute velocities, v _fp and v _fg, and the slide-roll ratio, S, of the current meshing point are introduced into the TEHL model as velocity parameters. The v _fp and v _fg, of the driving and driven gears, respectively, are obtained using the transmission characteristics of the VH-CATT cylindrical gear pair:

{\begin{matrix} v_{f p} = ω_{f p} \times r_{f p} = - ω_{p} (y_{f p} + R_{p}) i + ω_{p} x_{f p} j \\ v_{f g} = ω_{f g} \times r_{f g} = ω_{g} (y_{f g} - R_{g}) i - ω_{g} x_{f g} j \end{matrix}

(7)

Where ω _fp and ω _fg are the angular velocities of the meshing point, and R_p and R_g are the radii of the pitch circles.

The sliding velocity, v _r, at the meshing point can be calculated as follows:

\begin{matrix} v_{r} = v_{fp} - v_{fg} = - [ω_{p} (y_{fp} + R_{p}) + ω_{g} (y_{fg} - R_{g})] i \\ + (ω_{p} x_{fp} + ω_{g} x_{fg}) j \end{matrix}

(8)

The method for calculation u _e at any meshing point is as follows:

u_{e} = \frac{v_{fp} \cdot l_{r} \cdot l_{r} + v_{fg} \cdot l_{r} \cdot l_{r}}{2}

(9)

Where, l _r is the unit vector of the v _r.

Thus, the calculation method of S is as follows:

S = \frac{2 | v_{r} |}{| u_{e} |}

(10)

LTCA model

Generally, the real working conditions of gear pair can be simulated using LTCA model. The LTCA model of the VH-CATT cylindrical gear transmission is established based on the method proposed in Klima et al.³¹ and Jiang and Fang,³² which is illustrated in Figure 6.

Figure 6.

LTCA model of VH-CATT cylindrical gears.

In a meshing period, the VH-CATT cylindrical gear pair is in a contact state at a specific time, assuming that the contact load between Σ_p and Σ_g is distributed along the main direction of the Hertz contact ellipse. In Figure 6, the tooth curve is the intersection of the tooth surface and the normal plane along the relative main direction, ik (k = I, II) is a contact point, and jk (k = I, II) is a discrete point along the relative main direction. The gaps and contact loads of the contact points on each tooth surface in the meshing satisfy the following conditions:

{\begin{matrix} F_{k} f_{k} + w_{k} = Z + d_{k}, (k = I, II) \\ \sum_{j = 1}^{n} f_{jI} + \sum_{j = 1}^{n} f_{jII} = F \\ d_{jk} = 0, (if f_{jk} > 0) \\ d_{jk} > 0, (if f_{jk} = 0) \end{matrix}

(11)

Where F _k (n × n) is the assembled flexibility matrix determined using the finite element method, f _k (1 × n) and d _k (1 × n) are column matrixes composed of elements [f_1k, f_2k, …, f_nk] and [d_1k, d_2k, …, d_nk], respectively, w _k (1 × n) is a column matrix that contains the initial tooth gaps, and Z (1 × n) represents the approach of the tooth surface along the action line direction. f_jk (j = 1, 2, …, n) is the contact load at point j on tooth k, F is the total load, and d_jk (j = 1, 2, …, n) is the final tooth gap at point j on tooth k.

Thus, the maximum Hertzian contact pressure, p_h, in the contact area can be deduced as follows:

p_{h} = \frac{3 F}{2 π ab}

(12)

Where a (or b) is the minor (or major) semi-axis of the Hertz contact ellipse.

The a and b are expressed as follows:

{\begin{matrix} a = α_{ab} \sqrt[3]{\frac{3 F}{4 A} (\frac{1 - μ_{p}^{2}}{E_{p}} + \frac{1 - μ_{g}^{2}}{E_{g}})} \\ b = β_{ab} \sqrt[3]{\frac{3 F}{4 A} (\frac{1 - μ_{p}^{2}}{E_{p}} + \frac{1 - μ_{g}^{2}}{E_{g}})} \end{matrix}

(13)

Where E_p and E_g are the elastic moduli of the gear material, μ_p and μ_g are Poisson’s ratios of the gear material, and α_ab, A, and β_ab are coefficients.

Basic equations of the TEHL model

Referring to the model proposed in Zhu,³³ the Reynolds equation of the VH-CATT cylindrical gear pair with an elliptical contact is as follows:

\frac{\partial}{\partial x} (\frac{ρ h^{3}}{η} \frac{\partial p}{\partial x}) + \frac{\partial}{\partial y} (\frac{ρ h^{3}}{η} \frac{\partial p}{\partial y}) = 12 | u_{e} | \frac{\partial (ρ h)}{\partial x} + 12 \frac{\partial (ρ h)}{\partial t}

(14)

Where η is the viscosity of the lubricant, ρ is the density of the lubricant, h is the film thickness, p is the film pressure, and t is the meshing time.

Accordingly, the nondimensional Reynolds equation is

\frac{\partial}{\partial X} (ε_{X} \frac{\partial P}{\partial X}) + \frac{\partial}{\partial Y} (ε_{Y} \frac{\partial P}{\partial Y}) = \frac{\partial (\bar{ρ} H)}{\partial X} + \frac{\partial (\bar{ρ} H)}{\partial \bar{t}}

(15)

Where

ε_{X} = \frac{a}{12 | u_{e} |} \frac{p_{h}}{η_{0}} (\frac{\bar{ρ} H^{3}}{\bar{η}}), ε_{Y} = \frac{a^{3}}{12 b^{2} | u_{e} |} \frac{p_{h}}{η_{0}} (\frac{\bar{ρ} H^{3}}{\bar{η}})

(16)

Where η₀ is the initial viscosity of the lubricant.

The boundary conditions of equation (15) are

\begin{matrix} P (X_{in}, Y) = P (X_{out}, Y) = 0, P (X, Y_{in}) = P (X, Y_{out}) = 0, \\ \frac{\partial P}{\partial X} = \frac{\partial P}{\partial Y} = 0, X_{in} \leq X \leq X_{out}, Y_{in} \leq Y \leq Y_{out} \end{matrix}

(17)

To ensure that the TEHL model replicates a real scenario, the dynamic surface topography of the VH-CATT cylindrical gears is introduced into the film thickness equation in the form of the three-dimensional random roughness distribution. The film thickness equation is as follows:

h = h_{0} (t) + \frac{x^{2}}{2 R_{x}} + \frac{y^{2}}{2 R_{y}} + δ_{p} (x, y, t) + δ_{g} (x, y, t) + V (x, y, t)

(18)

Where h₀ is the initial central film thickness; R_x and R_y are the equivalent radii of the curvature at the current meshing point in the directions of the tooth profile and tooth trace, respectively, which are derived from equation (6); x²/2R_x and y²/2R_y are the transient geometries before elastic deformation; δ_p and δ_g are the three-dimensional random roughness of Σ_p and Σ_g, respectively. δ_p and δ_g are not introduced into equation (18) for smooth tooth surfaces.

The local contact deformation, V, is calculated using the Boussinesq integral equation:

V (x, y, t) = \frac{2}{π E'} \int \int_{Ω} \frac{p (ξ, ψ)}{\sqrt{{(x - ξ)}^{2} + {(y - ψ)}^{2}}} d ξ d ψ

(19)

Where $E'$ is the integrated elastic modulus.

Therefore, the nondimensional film thickness equation is:

\begin{matrix} H (X, Y, \bar{t}) = H_{0} (\bar{t}) + \frac{X^{2}}{2} + e_{k} \frac{Y^{2}}{2} + {\bar{δ}}_{p} (X, Y, \bar{t}) \\ + {\bar{δ}}_{g} (X, Y, \bar{t}) + L \int_{X_{in}}^{X_{out}} \int_{Y_{in}}^{Y_{out}} \\ \frac{P (\bar{ξ}, \bar{ψ})}{\sqrt{(X - \bar{ξ}) + (Y - \bar{ψ})}} d \bar{ξ} d \bar{ψ} \end{matrix}

(20)

Where e_k is the ellipticity, and e_k = R_x/R_y. L is the contact deformation coefficient of the contact ellipse, and L = 3WR_x/a²bE′.

Considering the influence of temperature change on the lubrication of VH-CATT cylindrical gears, the three-dimensional energy³⁴ is introduced into the TEHL model based on the theory of moving heat over a semi-infinite solid. The three-dimensional energy equation in the x, y, and z directions (z direction represents the film thickness direction) is as follows:

\begin{matrix} c_{a} ρ (u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} - w \frac{\partial T}{\partial z}) = κ (\frac{\partial^{2} T}{\partial z^{2}} + \frac{\partial^{2} T}{\partial y^{2}}) \\ - \frac{T}{ρ} \frac{\partial ρ}{\partial T} (u \frac{\partial p}{\partial x} + v \frac{\partial p}{\partial y}) + η [{(\frac{\partial u}{\partial z})}^{2} + {(\frac{\partial v}{\partial z})}^{2}] \end{matrix}

(21)

Where C_a and κ are the isobaric heat capacity and thermal conductivity of the lubricant, respectively, and u, v, and w are the flow rates along the x, y, and z directions, respectively.

To obtain the solution of the energy equation, the velocity field must be determined first. When the pressure distribution and the film thickness are known and the lubricant viscosity changes along the z direction, the flow velocities, u, v, and w, in the x, y, and z directions, respectively, are expressed as:

{\begin{matrix} \begin{matrix} u = | v_{fg} | + \frac{\partial p}{\partial x} (\int_{0}^{z} \frac{z}{η} dz - \frac{\int_{0}^{h} \frac{z}{η} dz}{\int_{0}^{h} \frac{1}{η} dz} \int_{0}^{z} \frac{1}{η} dz) + \frac{| v_{fp} | - | v_{fg} |}{\int_{0}^{h} \frac{1}{η} dz} \int_{0}^{z} \frac{1}{η} dz \\ v = \frac{\partial p}{\partial y} (\int_{0}^{z} \frac{z}{η} dz - \frac{\int_{0}^{h} \frac{z}{η} dz}{\int_{0}^{h} \frac{1}{η} dz} \int_{0}^{z} \frac{1}{η} dz) \\ w = - \frac{1}{ρ} \int_{0}^{z} [\frac{\partial (ρ u)}{\partial x} + \frac{\partial (ρ v)}{\partial z}] dz \end{matrix} \end{matrix}

(22)

The fundamental equation for calculating the temperatures of the upper and lower tooth surfaces in the meshing is as follows:

{\begin{matrix} T (x, y, 0) = T_{0} + κ {(π ρ_{p} c_{p} | v_{f p} | κ_{p})}^{- 0.5} \int_{- \infty}^{x} \frac{\partial T}{\partial z} \frac{d s}{{(x - s)}^{0.5}} \\ T (x, y, h) = T_{0} + κ {(π ρ_{g} c_{g} | v_{f g} | κ_{g})}^{- 0.5} \int_{- \infty}^{x} \frac{\partial T}{\partial z} \frac{d s}{{(x - s)}^{0.5}} \end{matrix}

(23)

Where T₀ is the initial temperature, (ρ_p, ρ_g), (c_p, c_g), and (κ_p, κ_g) are the density, specific heat capacity, and thermal conductivity of the tooth surface material of Σ_p and Σ_g, respectively.

Viscosity is used as the dependent variable of pressure and temperature to reflect their relationship. The Roelands viscosity–pressure–temperature equation is used to calculate the viscosity:

η = η_{0} \exp {(\ln η_{0} + 9.67) [{(1 + \frac{p}{p_{0}})}^{Z_{η}} {(\frac{T - 138}{T_{0} - 138})}^{- 1.1} - 1]}

(24)

Where Z_η is the coefficient of the equation, Z_η= α₀/5.1(lnη₀+ 9.67), and α₀ is the viscosity-pressure coefficient.

The Dowson–Higginson density–pressure–temperature equation is selected to reflect the specific effects of pressure and temperature on density, and the equation is as follows:

ρ = ρ_{0} [1 + \frac{0.6 p}{1 + 1.7 p} + D (T - T_{0})]

(25)

Where ρ₀ and D are the thermal expansion rate and initial density of the lubricant, respectively.

The load is balanced by the integrating the pressure over the entire solution region. The load-balancing equation is

\int_{x_{in}}^{x_{out}} \int_{y_{in}}^{y_{out}} p (x, y, t) dxdy = F (t)

(26)

Numerical computation method

Figure 7 shows the program flow of above models. The Reynolds equation and the energy equation are discretized iteratively by combining the central and forward difference methods. The Jacobi bipolar iteration and the Gauss–Seidel iteration methods is used in high-pressure and low-pressure regions, respectively. To prevent divergence in the iteration process, the relaxation factor, c₁, of the correcting pressure iteration is limited to 0.005–0.01, and the correction factor, c₂, of initial central film thickness is limited to 0.001–0.01. The convergence precisions of the pressure and the temperature are $| \frac{p_{i, j}^{New} - p_{i, j}^{Old}}{p_{i, j}^{New}} | 10^{- 5}$ and $| \frac{T_{i, j}^{New} - T_{i, j}^{Old}}{T_{i, j}^{New}} | 10^{- 5}$ for smooth surface solutions, respectively, and ≤10⁻⁴ for rough surface solutions. The solving domain is defined as −2.0a ≤ x ≤ 2.0a, −2.0b ≤ y ≤ 2.0b. The computing grid covering this domain consists of 299 ×299 equidistant nodes. The number of layers in the direction of the film thickness is 5. To get numerical solution of approximate actual condition, the roughness of both friction surfaces, is moved with the meshing time.^23,25 Since VH-CATT cylindrical gear has similar precision grade and working condition with the spiral bevel gear, the roughness is set as 0.3 μm according to the value in Pu et al.¹⁶ and Wang et al.²³ One of the three-dimensional rough surfaces is shown in Figure 8.

Figure 7.

Calculation flowchart.

Figure 8.

Three-dimensional rough surface profiles.

Model validation

In order to verify the present TEHL model of VH-CATT cylindrical gear transmission, the film thickness distribution solved by the present model is compared with the data solved by Wang et al.²² and measured by Omasta et al.³⁵ To ensure that the conditions for validation are identical, the same input parameters, solving domain, and the computing grid are obtained via from Wang et al.²² Table 2 and Figure 9 show that the central film thickness and film thickness distribution calculated by present model are substantially consistent with numerical simulation²² and measured result³⁵ at same working condition. Thus, the TEHL model has a certain computational accuracy, and the error is small compared with the results in the literature.

Table 2.

Central film thickness comparison.

h-cen (nm)	Data sources
388	Calculated by present TEHL model
392	Calculated by Wang et al.²²
395	Measured by Omasta et al.³⁵

Figure 9.

Comparison of film thickness distributions calculated by (a) present model, (b) Wang et al.,²² and (c) measured by Omasta et al.³⁵

Results and analysis

Results of TCA and LTCA

The design parameters of the VH-CATT cylindrical gear pair used for the TCA and the LTCA are listed in Table 1. In addition, the input rotational speed, n_in, of the driving gear is constant at 500 rpm, and the input torque of the driven gear is set as 232 Nm. The analysis results for the time-varying parameters along the contact trace in the VH-CATT cylindrical gears are generalized in Figure 10.

Figure 10.

The (a) composite curvature radii along directions of tooth profile and tooth trace, (b) sliding velocities of driving gear (red line) and driven gear (blue line), (c) entrainment velocities, (d) slide-roll ratio, (e) gaps between contacting teeth, (f) load distributions and maximum Hertz contact pressures, and (g) contact ellipses on tooth surface.

As shown in Figure 10(a), R_x reaches a maximum near the pitch point and exhibits a parabola-like variation with the meshing process. In comparison, R_y gradually increases linearly with the meshing time. As shown in Figure 10(b), the | v _fp| increases gradually while | v _fg| decrease, and reaches a same velocity near the pitch point. According to the calculation results in Figure 10(c) and (d), | u _e| and S gradually increase with the meshing time, and the value of S is 0 near the pitch point. Figure 10(e) gives the gaps between contacting teeth along the contact trace in a meshing period. As shown in Figure 10(f), both F and p_h have a sharp rise and fall in the middle of meshing period. Because of the transformation of the meshing process between single-tooth or double-teeth meshing, the area size of contact ellipse abruptly increases and decreases when the meshing process first enters the single-tooth and double-teeth meshing areas, respectively. A certain proportion of the transient contact ellipses is uniformly extracted from the path of contact and is projected on the tooth surface of driving gear, as shown in Figure 10(g).

Note that the rotation angle increment, Δϕ, of the gear is converted to the corresponding time increment, Δt, based on the input rotational speed, n_in, that is, Δt = Δϕ/n_in. In general, the meshing process of the studied gear pair alternately involves single or double teeth. Moreover, the values of some time-varying parameters between the single-tooth and double-teeth meshing areas rise and fall sharply. Therefore, a total of nine time nodes, t_i (i = 1, 2, 3, …, 9), are selected as typical examples for the analysis of the lubrication characteristics in a meshing period as shown in Figure 10(f). Note that the entire gear meshing process is disaggregated into three stages, which are the first meshing stage (t₁ ₋ t₃), the second meshing stage (t₄ ₋ t₆), and the third meshing stage (t₇ ₋ t₉), respectively.

Results of TEHL with smooth tooth surfaces

The time-varying parameters in Figure 10 and the lubrication parameters in Table 3 are taken as the input parameters of the TEHL model. Figures 11 to 13 give the numerical calculation results of the film thickness h, pressure p, and temperature rise, T, along the X- and Y-axes in a complete meshing period.

Table 3.

Material parameters of gear pair and lubricant.

Lubricant	Numerical value	Gear pair	Numerical value
Initial viscosity (Pa s)	0.09	Initial temperature (K)	303
Initial density (kg/m³)	846	Thermal conductivity	46
Specific heat (J/kg/K)	2000	Specific heat (J/kg/K)	470
Thermal conductivity (W/m K)	0.14	Density (kg/m³)	7850
Thermal expansivity (K⁻¹)	−0.00065	Elastic modulus (Pa)	2.07 × 10¹¹
Viscosity–pressure coefficient (Pa⁻¹)	2.3 × 10⁻⁸	Poisson’s ratio	0.3

Figure 11.

Film thickness distributions along X- and Y-axis directions of contact ellipse in meshing period.

Figure 12.

Pressure distributions along X- and Y-axis directions of contact ellipse in meshing period.

Figure 13.

Temperature rise distributions along X- and Y-axis direction of contact ellipse in meshing period.

h, p, and T, exhibit different variation trends in different meshing stages, as observed from Figures 11 to 13. The R_x, R_y, | u _e| gradually increase, and p_h decreases sharply in t₁ − t₃, so h increases significantly, whereas p and T decrease rapidly with t. In addition, S increase gradually while T decreases rapidly, enabling the viscosity of the lubrication to be recovered. Thus, the lubrication performance of the VH-CATT cylindrical gears in t₁ − t₃ is significantly improved with the meshing process. The sensitivity of h and p of numerical calculation results to the time-varying parameters in t₄ − t₆ is similar to that in t₁ − t₃. However, S reaches 0 near the pitch point (it’s near t₅), so T manifest as changing trend of “V” shape in this stage. T in Δ(t₄ − t₅) increases more than the decrement in Δ(t₅ − t₆), and this trend can be explained by the variation in p_h in this stage, which is shown in Figure 10. Notably, in t₇ − t₉, the variation laws of h, p, and T are different from those in t₁ − t₃ and t₄ − t₆. h and p continue to increase (or decrease) slightly during Δ(t₇ − t₈). Correspondingly, because the contact ellipse and R_x continue to decrease and p_h exhibits a slight increasing trend in Δ(t₈ − t₉), the variation trends of h and p are opposite to those in the previous period. In addition, with the increment of S (in Figure 10) in this stage, effect of lubricant viscosity on T gradually less than thermal expansion, so T increases continuously.

To evaluate the lubrication performance in the different meshing stages more intuitively, the contour lines of h, p, and T under smooth tooth surfaces are plotted in Figure 14.

Figure 14.

The contours of (a) film thickness h, (b) pressure p, and (c) temperature rise T in meshing period.

Figure 14 shows the significant differences in h, p, and T at different meshing moments, which are caused by the variations in the time-varying parameters at each moment. Therefore, the positions of the minimum film thickness, maximum pressure, and maximum temperature rise deviate with the meshing time, and at each meshing instant, a prominent squeezing effect occurs along the flow direction (X-axis) of the lubricant. During different meshing stages, that is, when the VH-CATT cylindrical gear pair operates under the transmission condition of a single-tooth and double-teeth alternately carrying the load, the h, p, and T change abruptly in the calculation area. This is because of the rapid variations in the load and the contact ellipse in a short time.

Based on the TEHL results for the VH-CATT cylindrical gears considering smooth tooth surfaces, Figure 15 gives the variations in the minimum film thickness h-min, center film thickness h-cen, maximum pressure p-max, and maximum temperature rise T-max in a meshing period. Here, define λ as the increment in the TEHL results at the current moment compared with that in the previous one, that is, λ = 100(C_r −P_r)/P_r%, where P_r and C_r represent the TEHL values at the previous and current meshing moments, respectively.

Figure 15.

The (a) minimum film thickness, (b) central film thickness (c) maximum pressure, and (d) maximum temperature rise.

First, Figure 15 shows that the lubrication condition is the least favorable at the instant when the gears initially mesh, when the occurrence of friction and wear failure are the most probable. The lubrication condition gradually improves as the meshing progresses. Second, because the time-varying parameters change significantly in Δ(t₁ − t₂), λ increases or decreases rapidly in this period. The first peak values (0.63 ms) of the four λ curves are 12.49%, 8.02%, −7.61%, and −7.10%, respectively. Concurrently, it is noted that the λ values clearly fluctuate at 3.74 and 5.71 ms caused by the switch between single-tooth and double-teeth load. Finally, when t is 4.62 ms, the gear pair is located near the pitch point, and S is 0 at this time. Thus, h-min,p-max, and T-max have good extreme values of obvious lubrication near the pitch point, which are 0.8036 μm, 0.6948 GPa, and 38.6125°C.

Results of TEHL with rough tooth surfaces

To predict the lubrication characteristics more accurately in a meshing period of a real machined surface, this section introduces the three-dimensional surface roughness based on the input parameters of smooth TEHL. The TEHL results considering rough tooth surfaces for the VH-CATT cylindrical gears are summarized in Figures 16 to 18.

Figure 16.

Film thickness distributions along X- and Y-axis directions of contact ellipse in meshing period.

Figure 17.

Pressure distributions along X- and Y-axis directions of contact ellipse in meshing period.

Figure 18.

Temperature rise distributions along X- and Y-axis directions of contact ellipse in meshing period.

Figures 16 to 18 show that the TEHL numerical simulation results for rough tooth surfaces are significantly different from those of smooth tooth surfaces, which are shown in Figures 11 to 13. Under the effect of roughness, although the variation laws of the TEHL results for smooth tooth surfaces still apply in each meshing stage for rough tooth surfaces, the latter exhibit clear oscillations. The temperature rise between the different meshing stages for rough tooth surfaces exhibits numerous divergent peaks, as seen in Figure 18. In a meshing period, the smaller film thickness and larger pressure are, the higher temperature is. Therefore, the first meshing stage is most likely subjected to the influence of direct contact of the rough peaks by comparing the peaks of temperature at different meshing stages.

To facilitate direct observation of the lubrication performance of the VH-CATT gear pair under thermal and roughness effects, Figure 19 presents contour maps.

Figure 19.

The contours of (a) film thickness h, (b) pressure p, and (c) temperature rise T in meshing period.

In Figure 19(a), the squeezing effect caused by the extrusion of lubrication at each meshing instant is still evident; however, the shape of the film fluctuates owing to the roughness of the gear surface. The pressure in Figure 19(b) exhibits abrupt changes between different meshing stages. A large ratio of the contact load to the contact ellipse implies a high viscosity of lubrication and a strong influence of the roughness on the pressure. In Figure 19(c), the influence of roughness on the lubrication performance of the VH-CATT cylindrical gears can be manifested by the number of dark red spots in the contact ellipse. Note that the there are some dark red spots randomly distributed on the right side of the contact ellipse at t₈ and t₉ in red rectangular box of Figure 19(c). The reasons for this phenomenon are the random distribution of roughness and the poor lubrication performance in t₇ − t₉, so it is still possible to damage the oil film by the roughness of the tooth surface.

Conclusions

In this study, the TEHL model of VH-CATT cylindrical gear transmission in elliptical contact has been established by considering the effects of time-varying parameters. The numerical method for calculating the time-varying parameters of this gear transmission in a meshing period has been proposed. After introducing the time-varying parameters in a meshing period, the lubrication characteristics considering smooth and rough tooth surfaces have been investigated. Based on numerical simulation, the following conclusions can be drawn:

The kinematic parameters, spatial geometric parameters, and contact parameters of VH-CATT cylindrical gear transmission constantly vary in a meshing period. To obtain these time-varying parameters, the TCA and LTCA models are established. Considering the effects of temperature and time-varying parameters for the analysis of the lubrication characteristics, a TEHL model of the VH-CATT cylindrical gears in elliptical contact is presented.

The TEHL results for smooth tooth surfaces show that the lubrication characteristics of the VH-CATT cylindrical gears are directly affected by time-varying parameters. The lubrication characteristics in the entire meshing process fluctuate because the load distribution and Hertz contact ellipse vary significantly. The lubrication performance at the beginning of the first meshing stage—when the film thickness is the smallest, the pressure is the largest, and the temperature rise is highest—is the worst. However, the performance improves rapidly as the meshing progresses. The second meshing stage presents the best lubrication performance. The lubrication performance of the third meshing stage is better than that of the first meshing stage and inferior to that of the second meshing stage.

The TEHL results for rough tooth surfaces show that the lubrication characteristics obtained after considering roughness are closer to the actual lubrication characteristics. The overall trends of the variations in the lubrication performance in a meshing period for rough tooth surfaces are similar to those observed for smooth tooth surfaces. However, the positions of the minimum film thickness, maximum pressure, and maximum temperature rise shift are more stochastic than the corresponding positions for smooth tooth surfaces. The stage most likely subjected to the influence of direct contact due to rough peaks is the first meshing stage, followed by the third meshing stage.

Footnotes

Appendix

Handling Editor: Chenhui Liang

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The present study was founded by the National Natural Science Foundation of China Project Nos. 51375320 and 51875370, the Foundation of the State Key Laboratory of Mechanical Transmission of Chongqing University Project No. SKLMT-KFKT-201807, and the Natural Science Foundation of Gansu Province Project No. 18JR3RA140.

ORCID iD

Pei Luo

References

Tseng

Tsay

. Contact characteristics of cylindrical gears with curvilinear shaped teeth. Mech Mach Theory 2004; 39: 905–919.

Andrei

Walton

, et al. Analysis of gear quality criteria and performance of curved face width spur gears. In: The 10th International Conference on Tribology “ROTRIB’07”, Bucharest, Romania, 2008, vol.14, pp.28–32.

Tseng

Tsay

. Mathematical model and surface deviation of cylindrical gears with curvilinear shaped teeth cut by a hob cutter. J Mech Des 2004; 127: 271–277.

Chen

Tsay

, et al. Contact characteristics of circular-arc curvilinear tooth gear drives. J Mech Des 2009; 131: 1–8.

Zhao

Hou

Duan

, et al. Research on the forming theory analysis and digital model of circular arc gear shaped by rotary cutter. J Sichuan Univ 2016; 48: 119–125.

Hou

, et al. Milling machine error modelling and analysis in the machining of circular-arc-tooth-trace cylindrical gears. Trans FAMENA 2021; 44: 13–29.

Martin

. Lubrication of gear teeth. Engineering 1916; 102: 119–121.

Kahraman

. Prediction of friction-related power losses of hypoid gear pairs. Proc IMechE, Part K: J Multi-body Dynamics 2007; 221: 387–400.

Liu

Zhu

Liu

. A micro-TEHL finite line contact model for a helical Gear pair. Adv Mech Eng 2015; 7: 04790.

10.

Mohammadpour

Theodossiades

Rahnejat

. Elastohydrodynamic lubrication of hypoid gear pairs at high loads. Proc IMechE, Part J: J Engineering Tribology 2012; 226: 183–198.

11.

Mohammadpour

Theodossiades

Rahnejat

, et al. Non-Newtonian mixed elastohydrodynamics of differential hypoid gears at high loads. Meccanica 2014; 49: 1115–1138.

12.

Bobach

Beilicke

Bartel

, et al. Thermal elastohydrodynamic simulation of involute spur gears incorporating mixed friction. Tribol Int 2012; 48: 191–206.

13.

Xue

Qin

. The scuffing load capacity of involute spur gear systems based on dynamic loads and transient thermal elastohydrodynamic lubrication. Tribol Int 2014; 79: 74–83.

14.

Beilicke

Bobach

Bartel

. Transient thermal elastohydrodynamic simulation of a DLC coated helical gear pair considering limiting shear stress behavior of the lubricant. Tribol Int 2016; 97: 136–150.

15.

Wang

Xiang

, et al. Effect of roughness on meshing power loss of planetary gear set considering elasto-hydrodynamic lubrication. Adv Mech Eng 2020; 12: 1–1.

16.

Wang

Yang

, et al. Mixed elastohydrodynamic lubrication with three-dimensional machined roughness in spiral bevel and hypoid gears. J Tribol 2015; 137: 1–11.

17.

Zhang

Liu

Fang

. On the prediction of friction coefficient and wear in spiral bevel gears with mixed TEHL. Tribol Int 2017; 115: 535–545.

18.

Simon

. Multi-objective optimization of hypoid gears to improve operating characteristics. Mech Mach Theory 2020; 146: 103727.

19.

Simon

. Improved mixed elastohydrodynamic lubrication of hypoid gears by the optimization of manufacture parameters. Wear 2019; 438-439: 102722.

20.

Simon

. Influence of machine tool setting parameters on EHD lubrication in hypoid gears. Mech Mach Theory 2009; 44: 923–937.

21.

Wang

Shao

Wang

, et al. Effect of load and speed on warship power rear drive system helical gear anti-scuffing properties in thermal elasto-hydrodynamic lubrication. Adv Mech Eng 2017; 9: 1–10.

22.

Wang

Ren

Zhang

, et al. A mixed TEHL model for the prediction of thermal effect on lubrication performance in spiral bevel gears. Tribol Trans 2020; 63: 314–324.

23.

Wang

, et al. Numerical simulation of transient mixed elastohydrodynamic lubrication for spiral bevel gears. Tribol Int 2019; 139: 67–77.

24.

Pei

Huang

, et al. Dynamical wear prediction along meshing path in mixed lubrication of spiral bevel gears. Adv Mech Eng 2020; 12: 1–16.

25.

Wang

Zhang

, et al. Transient behaviors of friction, temperature and fatigue in different contact trajectories for spiral bevel gears. Tribol Int 2020; 141: 105965.

26.

Sun

Liu

Zhao

, et al. EHL analysis of spiral bevel gear pairs considering the contact point migration due to deformation under load. Math Probl Eng 2020; 2020: 1–19.

27.

Yongqiao

Shuhong

, et al. Numerical analysis of isothermal elastohydrodynamic lubrication of cylindrical gears with variable hyperbolic circular arc and tooth trace. Trans FAMENA 2018; 42: 61–72.

28.

Wei

Liu

, et al. Research on thermal elastohydrodynamic lubrication of cylindrical gears with curvilinear shaped teeth. J Chin Soc Mech Eng 2018; 39: 451–458.

29.

Litvin

Krylov

Erikhov

. Generation of tooth surfaces by two-parameter enveloping. Mech Mach Theory 1975; 10: 365–373.

30.

Litvin

Fuentes

. Gear geometry and applied theory. Cambridge: Cambridge University, 2004. 153–194.

31.

Klima

Zhang

Fang

. Analysis of tooth contact and load distribution of helical gears with crossed axes. Mech Mach Theory 1999; 34: 41–57.

32.

Jiang

Fang

. Design and analysis of modified cylindrical gears with a higher-order transmission error. Mech Mach Theory 2015; 88: 141–152.

33.

Zhu

. On some aspects of numerical solutions of thin-film and mixed elastohydrodynamic lubrication. Proc IMechE, Part J: J Engineering Tribology 2007; 221: 561–579.

34.

Zhu

Wen

. A full numerical solution for the thermoelastohydrodynamic problem in elliptical contacts. J Tribol 1984; 106: 246–254.

35.

Omasta

Křupka

Hartl

. Effect of surface velocity directions on elastohydrodynamic film shape. Tribol Trans 2013; 56: 301–309.

TEHL analysis of VH-CATT cylindrical gear transmission in elliptical contact considering time-varying parameters

Abstract

Keywords

Introduction

Contact analysis of VH-CATT cylindrical gears

Tooth surface equation

TCA model

LTCA model

Basic equations of the TEHL model

Numerical computation method

Model validation

Results and analysis

Results of TCA and LTCA

Results of TEHL with smooth tooth surfaces

Results of TEHL with rough tooth surfaces

Conclusions

Footnotes

Appendix

Declaration of conflicting interests

Funding

ORCID iD

References