Next Article in Journal
Orientation Modeling Using Quaternions and Rational Trigonometry
Next Article in Special Issue
Optimal Design for Vibration Mitigation of a Planar Parallel Mechanism for a Fast Automatic Machine
Previous Article in Journal
Assessment of Industry 4.0 for Modern Manufacturing Ecosystem: A Systematic Survey of Surveys
Previous Article in Special Issue
Embedded Payload Solutions in UAVs for Medium and Small Package Delivery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Friction-Induced Efficiency Losses and Wear Evolution in Hypoid Gears

Dipartimento di Ingegneria Civile e Industriale, Università di Pisa, Largo Lucio Lazzarino 1, 56122 Pisa, Italy
*
Author to whom correspondence should be addressed.
Machines 2022, 10(9), 748; https://doi.org/10.3390/machines10090748
Submission received: 30 July 2022 / Revised: 24 August 2022 / Accepted: 25 August 2022 / Published: 29 August 2022
(This article belongs to the Special Issue Advances of Machine Design in Italy 2022)

Abstract

:
A correct methodology to evaluate the friction coefficient in lubricated gear pairs is paramount for both the estimation of energy losses and the prediction of wear. In the first part of the paper, a methodology for estimating the coefficient of friction with a semi-empirical formulation is presented, and its results are also employed to analyze mechanical efficiency losses in a hypoid gearset. Hypoid gears have complex tooth surface geometries, and the entraining kinematics of the lubricant is quite involved. The second part of the paper showcases a simulated wear investigation based on the Archard model. The main focus is on the impact of the frequency adopted for updating the worn geometry of the gear and pinion teeth on the fidelity and consistency of the tribological outcomes. These are measured in terms of overall quantity of material removed and characteristics of the loaded contact pattern. More in detail, a sensitivity analysis is presented that compares the total wear of a hypoid gearset after 30 million cycles estimated using different geometry update steps. Contact pressures, which are necessary to perform the aforementioned analyses, are calculated through an accurate, state-of-the-art loaded tooth contact analysis solver.

1. Introduction

Hypoid gears are associated with high sliding speeds during meshing due to their screwing relative motion. This brings about peculiar and, in general, more severe wear patterns than other gear types. Such wear can result in a contact pattern that is drastically different from the designed one. In fact, since an optimal contact pattern is the outcome of a micro-geometry optimization process [1,2,3], a variation of just a few tens of microns can significantly change the contact properties. Hence, wear must be cautiously considered for gears operating in long-life applications.
Park et al. [4] introduced a wear model based on Archard’s law interfaced with a finite element (FEM)-based contact model. The very same loaded tooth contact analysis (LTCA) tool, Ansol Transmission3D [5], is employed in this work. However, the software’s capabilities have been greatly improved over the last few years, thus allowing a more accurate prediction of the contact properties and, therefore, more accurate analyses of both lubrication and wear.
The study in [6] proposed a wear model interfaced with a semi-analytical contact model developed by Kahraman and Kolivand [7]. This tool greatly reduced the computation cost to perform the LTCA simulations while trading some accuracy for efficiency. Park [8] developed a “patching” surface interpolation technique to predict wear, which allows the employment of fewer time-step discretizations of the meshing cycle. More recently, Ref. [9] investigated wear in hypoid gears and its experimental correlation with the loaded transmission error.
In the first part of this paper, a semi-empirical formulation is introduced and employed to analyze lubrication and friction-induced efficiency losses during gear meshing. Despite the simplicity of the proposed model, the estimates it provides are effective and have been validated with experimental data in [10].
Efficiency losses are an important aspect of geared transmissions. They can be classified into load-independent and load-dependent losses. The first category includes the so-called churning and windage losses (i.e., pumping of the lubricant between the mating members and splashing caused by inertial effects). The second category is related to rolling and sliding frictional losses. The rolling losses are usually negligible compared to the sliding ones, especially in hypoid gears. Our work proposes an estimation method for the friction coefficient under different lubrication regimes, which is the fundamental parameter required to accurately predict friction-induced efficiency losses.
Many friction prediction models for hypoid gears are available in the literature. Contribution [11] evaluated the lubrication performance under different possible contact paths on bevel gears, [12] carried out an elasto-hydrodynamic lubrication analysis on hypoid gears under relatively high loads, while [13] analyzed the lubrication of hypoid gears taking into account the three-dimensional surface roughness. However, despite their accurate tribological analyses, many of the cited contributions fall short in providing reliable contact analysis results. In fact, an accurate LTCA tool has a paramount importance for the subsequent estimation of the friction coefficient.
In the second part of the paper, a wear investigation is performed. The main goal is to establish how often the worn tooth flank geometry should be updated (properly accounting for the removed material) to achieve a reliable and consistent prediction of the final wear and contact pattern. To this end, we computed wear predictions after 30 million meshing cycles (counted on the pinion) using different geometry update steps. This allowed the performance of a sensitivity analysis and provide quantitative indications on the minimum number of geometry updates required to provide a reliable prediction with limited computational burden.

2. Estimation of the Lubricated Friction Coefficient

There are several factors that influence friction, such as the lubrication regime (full film, mixed, boundary), the behavior of the lubricant with varying operating conditions (temperature, pressure, shear rate) and the surrounding environment (the boundary conditions of the lubricated contact).
The friction coefficient f is evaluated from boundary to full-film lubrication conditions using a load sharing function g ( Λ ) , related to the portion of load supported by the full film. According to [14], f can be calculated as:
f = f h g ( Λ ) 1.2 + f b ( 1 g ( Λ ) )
where f h and f b are the friction coefficients related to the hydrodynamic (full fluid) and the boundary lubrication conditions, respectively. Λ is the ratio between the film thickness h and the equivalent surface roughness of the contacting bodies:
Λ = h R q 1 2 + R q 2 2
Usually, the film thickness h is assumed to be the central film thickness h c , and R q 1 and R q 2 are the root mean square roughnesses of the two surfaces.
Several expressions can be found for g ( Λ ) in the literature. The formulation in [14] is employed in this work, which reads:
g ( Λ ) = 0.84 Λ 0.23 if Λ 2 ; g ( Λ ) = 1 if Λ > 2
The coefficient of friction f b is considered constant (with f b = 0.08 being the typical value [14]), while the mean value of f h can be evaluated as the ratio between the shear stress τ and the mean contact pressure p m :
f h = τ p m
For elliptical nonconformal contacts, p m = F π a b , where a and b are the axis semi-widths of the contact ellipse and F is the normal load.
The estimation of τ is performed according to the Eyring constitutive model:
τ = τ E sinh η γ ˙ τ E 1 ,
where τ E is the Eyring shear stress, η is the dynamic viscosity, and the shear strain rate is approximated as γ ˙ = Δ u h c , with Δ u being the sliding velocity between the two mating teeth at the nominal contact point.
The dynamic viscosity variation with temperature T and pressure p is described by the empirical formulas introduced by Vogel and Gold [15]:
η = η 0 e B T C + α P p ,
where the viscosity-pressure coefficient is calculated as follows:
α P = 10 5 a 1 + a 2 T + ( b 1 + b 2 T ) ( p 10 5 )
The remaining coefficients in Equations (6) and (7) need to be estimated empirically, based on the properties of the specific lubricant.
Different expressions for the evaluation of the central film thickness can be found in the literature. Four different lubrication regimes can be observed for nonconformal contacts, depending on the elastic deformation of the bodies and the variation of viscosity with pressure. The four regimes are usually indicated as isoviscous-rigid IR, piezoviscous-rigid PR, isoviscous-elastic IE and piezoviscous-elastic PE. The formulas reported in [16] were elaborated in [17] for the more general case in which the entraining velocity is not collinear with any principal direction. The final expressions employed in this work read as follows:
h c I R = 128 η 0 2 u 2 R e 3 F 2 R s R e 0.131 arctan R s 2 R e + 1.683 2 1 + 2 R e 3 R s 2
h c P R = 141 η 0 1.25 α 0.375 u 1.25 R e 1.5 F 0.875 1 e 0.0387 R s R e
h c I E = 11.15 η 0 0.66 u 0.66 R e 0.766 F 0.213 E 0.447 1 0.72 e 0.28 ( R s R e ) 2 π
h c P E = 3.61 η 0 0.68 α 0.53 u 0.68 R e 0.446 E 0.087 F 0.063 1 0.61 e 0.73 ( R s R e ) 2 π
Here, E is the equivalent elastic modulus, u is the entraining velocity, and R e and R s are the equivalent radii of curvature parallel and perpendicular to u, respectively. An accurate description of u, R e and R s is provided in Section 3. The proper lubrication regime is determined in a practical way by taking the highest of the values given by Equations (8)–(11).
Thermal effects are included using a reduced value for the central film thickness obtained by multiplying it by a dimensionless reduction factor Φ . Several models for Φ are available. The one described in [18] is used here:
Φ = 1 1 + 0.1 ( 1 + 8.33 S 0.83 ) L 0.64
where S = Δ u / u and L = β η 0 u 2 k is the dimensionless thermal loading parameter.
The coefficients that characterize the 75W90 oil, which is the lubricant used in this study, are shown in Table 1. Those experimental values have been obtained from [19].
It is worth remarking that the formulas introduced in this section provide a first approximation of the lubricant film thickness and coefficient of friction. In particular, very high sliding speeds are associated with large increments of temperature that may cause a significant decrease in viscosity and hence in film thickness. The problem in its full generality would be also time-dependent, but some preliminary film thickness evaluation can be made considering a stationary situation with constant speed and geometry per contact configuration of the mating teeth. It may be worth noting in passing that the film thickness would also be influenced by the presence of spin, which is neglected here for simplicity.

3. Generation and Kinematic Analysis of Hypoid Gears

The face-milled tooth flank geometries are digitally synthesized through simulation of the generation process (by envelope) during finishing (grinding). The enveloping motions are defined by the kinematic parameters (machine settings) of the 9-axis Gleason hypoid generator and the geometry of the grinding tool (tool settings). A special program developed by the authors [20] allows the obtaining of those settings from the basic macro-geometry data. A plot of the hypoid gearset under investigation is shown in Figure 1.
Through geometric modelling and simulated loaded contact analysis (LTCA), it is also possible to compute fundamental kinematic and geometric information at each contact point, such as (subscript j = 1 / 2 refers to the gear/pinion, respectively; superscript i = 1 / 2 refers to the first/second direction on each surface):
  • velocities u 1 and u 2 of the two bodies (teeth) at the contact point;
  • angular velocities ω 1 and ω 2 of the bodies;
  • principal curvatures of the bodies K 1 ( 1 ) , K 1 ( 2 ) , K 2 ( 1 ) and K 2 ( 2 ) and their directions τ 1 ( 1 ) , τ 1 ( 2 ) , τ 2 ( 1 ) and τ 2 ( 2 ) .
  • principal relative curvatures K x and K y and their associated directions τ x and τ y . The direction along the largest semi-width of the contact ellipse is denoted as τ y .
The principal normal curvature K j ( i ) of body j along direction τ i is considered positive/negative if the normal plane through τ i cuts body j producing a section with a convex/concave boundary curve.
To employ the formulas described in Section 2, the entraining velocity u of the lubricant needs to be computed. In addition, the direction of such velocity does not coincide with any of the principal directions at the contact point. However, given the angular velocities of the bodies, their linear velocities at the contact point and their curvatures K 1 x , K 1 y for the gear wheel and K 2 x , K 2 y for the pinion, respectively, along the principal relative directions τ x and τ y of the equivalent contact, the components of the entraining velocity vector can be obtained as [17,21]:
u x = ( ω 2 y ω 1 y ) 1 K 1 x + K 2 x + 1 2 ( u 1 x u 2 x ) K 2 x K 1 x K 1 x + K 2 x
u y = ( ω 1 x ω 2 x ) 1 K 1 y + K 2 y + 1 2 ( u 1 y u 2 y ) K 2 y K 1 y K 1 y + K 2 y
These are the components of the entraining velocity vector in the local frame with its axes aligned with the principal directions of the equivalent contact. The norm of the entraining velocity vector u = u x 2 + u y 2 is employed for the estimation of the friction coefficient formulas introduced in the previous section. It is important to remark that the validity of Equations (13) and (14) is subject to the following hypotheses:
(1)
body 1 (tooth of the gear wheel) is always convex along both its principal directions (which is true under the so-called Drive operating conditions, i.e., the gear convex tooth sides make contact with the pinion concave tooth sides);
(2)
the cross product τ x × τ y yields a unit vector always pointing inside body 1.
As a final step, the radii of curvature along the entraining and the side-leakage directions ( R e and R s ) can be obtained by [22]:
R e = ( K x cos 2 θ + K y sin 2 θ ) 1 R s = ( K x sin 2 θ + K y cos 2 θ ) 1 ,
where θ = arctan ( u y / u x ) is the angle between the entraining velocity and τ x .

4. Numerical Application on a Hypoid Gearset

An application of the aforementioned friction coefficient model is presented in this section. The basic data of the hypoid gearset employed in our study are shown in Table 2.
The machine-tool settings computed by a special program developed by the authors, and described in detail in [20], are given as input to Transmission3D to perform an accurate LTCA analysis. More in detail, the contact simulation is performed by discretizing the overall meshing cycle into N steps = 15 contact configurations (time steps). The analysis results are then extracted and post-processed to carry out the calculations for the estimation of the friction coefficient. We assumed a pinion torque of 250 Nm and a pinion speed of 2000 rpm.
Figure 2 shows the load shared by a mating tooth pair over a mesh cycle (as a function of the pinion rotation angle φ 2 ). Figure 3 shows the entraining velocity of the lubricant. Figure 4 shows the evolution of the entraining and the side-leakage radii. Those values are strongly dependent on the micron-level flank deviations from the conjugate surfaces. As a matter of fact, micro-geometry can drastically change the evolution of the contact path, which may result in different contact zones and thus different local curvatures.
In Figure 5 we can observe the different central film thicknesses calculated for the different lubrication regimes. At each contact configuration, the lubrication regime is the piezo-viscous-elastic one, probably due to the high loading of the gearset. This results in an analogous trend of the Λ factor, which can be observed in Figure 6. It is clear that the meshing tooth surfaces always operate under mixed lubrication conditions. To increase Λ , a lower surface roughness of the tooth flanks should be aimed for. Finally, the coefficient of friction is shown in Figure 7.
The efficiency loss can be computed by numerical integration of the frictional power losses over the meshing cycle:
E l = i = 1 N steps 1 1 2 ( P i + P i + 1 ) Δ t N steps 1
where P i = k = 1 n T f i k F i k u i k ( s ) is the power loss at the i-th time step and f i k , F i k , u i k ( s ) are, respectively, the friction coefficient, normal contact force and sliding velocity of the k-th tooth pair in engagement at the i-th time step.
The efficiency estimation has been performed under different torques and speeds, namely in a torque range between 50 and 300 Nm and in a speed range between 500 and 3000 rpm. The results are shown in Figure 8. According to our model, the efficiency decreases at large torque values and increases at higher speed. Those results are consistent with a mixed lubrication regime: a larger torque decreases the film thickness, resulting in a larger probability of asperity contacts. On the other hand, a higher speed increases the film thickness due to hydrodynamic effects, which mitigates the adverse effect of the boundary friction coefficient f b (cfr. Equation (1)).

5. Wear Model

The wear pattern is extracted directly by Transmission3D. Its wear implementation is based on the well-known Archard model:
d w d s = k P
where w is the wear depth, s is the sliding distance, k is the wear coefficient, and P is the local pressure. The numerical value of the coefficient k needs to be determined experimentally. Kahraman et al. [23] suggested a value of 9.65 × 10 19 m 2 / N , which was experimentally derived for helical gears under mixed lubrication conditions. Due to the complexity and cost of an experimental test rig, this value has been used also for the hypoid gears investigated in this work. Under comparable operating conditions in terms of loading and peripheral speeds, hypoid gears generally have more sliding than helical gears, thus we can hypothesize that the value of k is higher for the former. This might be caused by the fact that the larger sliding speeds entail an increased heating of the lubricant, which decreases the film thickness and thus induces more asperity contacts (more wear) under a mixed lubrication regime. However, obtaining an accurate k value for hypoid gears is out of the scope of the present work.
After a certain number of meshing cycles, wear can noticeably modify the tooth flanks, thus the tooth geometry needs to be updated (before simulating further wear) to re-evaluate the contact pressure distribution, which is very sensitive to micro-geometry variations. The geometric update is performed by superimposing the post-processed wear pattern (evaluated at points arranged in a grid on the tooth flanks) onto the previous tooth surfaces. At the j-th grid point:
p j ( k ) = p j ( k 1 ) + n j ( k 1 ) w j ( k ) ,
where:
  • p j ( k ) contains the coordinates of the grid point after the k-th geometry update;
  • p j ( k 1 ) and n j ( k 1 ) are the grid point coordinates and unit normal components after the ( k 1 ) -th (previous) geometry update;
  • w j ( k ) is the k-th local wear depth.
The updated points p ( k ) are then best-fit by a NURBS surface [24], which allows us to also compute the updated unit normals n ( k ) . The updated points and normals are then fed back to a tooth mesher for Transmission3D, and the next wear simulation step is executed.

6. Wear Analysis Application

In this section, we present the results of an investigation aimed at assessing the impact of the frequency in updating the geometry of the worn tooth surfaces on the prediction of the final accumulated wear and contact patterns. The contact pattern represents the envelope of all the instantaneous contact zones over the meshing cycle; the wear pattern is the corresponding wear distribution on the tooth surface (after a certain number of wear cycles). Each geometry update necessarily calls for a new LTCA simulation, hence frequent updates can easily lead to a significant cost in terms of CPU time for a complete wear simulation. Our goal is to evaluate the best trade-off between the number of geometry updates and the computational burden. We assume here a service life of 30 Mc (where 1 Mc = 10 6 cycles), counted as revolutions of the hypoid pinion. The initial contact patterns of pinion and gear are shown in Figure 9.
As a first step, serving as a baseline result, a very coarse analysis is carried out, where the predicted wear is evaluated assuming that the initial tooth geometry does not change during the whole service life. After calculating wear distribution, a geometry update is performed only once and the resulting contact pressures are computed. The corresponding wear and pressure patterns are shown in Figure 10 for the pinion and in Figure 11 for the gear. The contact patterns exhibit remarkable differences with respect to the ones in Figure 9: a significant shift of the pressure peaks away from the tooth flank center is evident (Figure 10a and Figure 11a), and it is due to an unrealistic wear indentation effect.
In the subsequent analysis, a similar wear simulation is performed. This time, however, the tooth flank modifications due to the worn-out material, and the corresponding variation in contact pressures, are gradually updated according to a stepwise approach, where each step consists of 5 Mc. In other words, the tooth flank geometry is iteratively modified every 5 Mc as per Equation (18). This update strategy, in practice, is equivalent to updating after a given wear volume, i.e., the wear volume at the first wear cycle, multiplied by the number of wear cycles [25]. The evolution of the pinion contact and wear patterns are shown, respectively, in Figure 12 and Figure 13. The associated gear patterns are shown in Figure 14 and Figure 15. As expected, the overall pattern grows larger along both the profile and face directions of the active flank. However, a localization of the contact pressures near the fillet portion of the pinion (close to, but away from, the gear tip) is evident. In the same area, also the wear pattern features a localization of its peaks. Figure 16 shows the pinion instantaneous contact zones at a specific meshing configuration before and after wear. The corresponding instantaneous contact zones for the gear are shown in Figure 17. It can be noticed how wear induced a greater overlap ratio, increasing the number of mating tooth pairs (from one to three). Pressure localization is evident here as well, which appears to be caused by an abrupt change to the local curvature due to wear.
To gather more data for a sensitivity analysis, wear simulations have been performed also with geometry updates every 10 Mc, 1 Mc and 0.5 Mc. A side-by-side comparison of the final contact patterns is shown in Figure 18. For brevity, only the pinion contact patterns are shown. It can be observed that the patterns obtained with 6, 30 and 60 update steps, i.e., every 5 Mc, 1 Mc, and 0.5 Mc, seem almost identical. Even with just three updates, i.e., every 10 Mc, the final pattern is already reasonably close to our most accurate analysis.
Figure 19 depicts a more quantitative result of this analysis by showing the maximum contact pressure ( P m a x ) registered during meshing as a function of the wear cycles. An analogous representation of the maximum wear depth w m a x registered on the pinion is shown in Figure 20. The graphs show that a geometry update performed every 5 Mc may be a good choice for balancing an accurate prediction, both in terms of contact pressure and wear depth, with affordable computations in terms of CPU time.

7. Conclusions

In the first part of the paper, a semi-empirical formulation of the friction coefficient has been proposed to estimate friction-induced losses in lubricated hypoid gears, whose tooth surfaces have been digitally synthesized through simulation of the grinding process. The calculations leverage both special procedures developed by the authors to compute basic geometric and kinematic quantities and an ad hoc LTCA simulation tool for an accurate prediction of normal loads and contact pressures. Our streamlined methodology is suitable for integration into an automatic optimization pipeline where the complex interactions between the flank geometry and the coefficient of friction are profitably accounted for to develop more efficient geared transmissions.
In the second part of the paper, a wear study has been presented. The aim was to assess the impact of the frequency adopted for updating the worn geometry of the gears on the accuracy and consistency of the tribological results. Interestingly, in this preliminary study, updating the geometry every five million cycles (corresponding to a wear variation of about 5 micrometers) turned out to be sufficient to suitably capture the evolution of both the maximum contact pressure and the maximum wear depth.

Author Contributions

Conceptualization, A.A., M.G. (Marco Gabiccini), L.M., F.D.P. and E.C.; Data curation, E.G.; Methodology, E.C.; Software, E.G.; Supervision, A.A., M.G. (Marco Gabiccini), M.G. (Massimo Guiggiani), L.M., F.D.P. and E.C.; Writing—original draft, E.G.; Writing—review and editing, A.A. and M.G. (Marco Gabiccini). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Artoni, A.; Bracci, A.; Gabiccini, M.; Guiggiani, M. Optimization of the loaded contact pattern in hypoid gears by automatic topography modification. J. Mech. Des. 2009, 131, 0110081–0110089. [Google Scholar] [CrossRef]
  2. Vivet, M.; Mundo, D.; Tamarozzi, T.; DEsmet, W. An analytical model for accurate and numerically efficient tooth contact analysis under load, applied to face-milled spiral bevel gears. Mech. Mach. Theory 2018, 130, 137–156. [Google Scholar] [CrossRef]
  3. Wang, Q.; Zhou, C.; Gui, L.; Fan, Z. Optimization of the loaded contact pattern of spiral bevel and hypoid gears based on a kriging model. Mech. Mach. Theory 2018, 122, 432–449. [Google Scholar] [CrossRef]
  4. Park, D.; Kahraman, A. A surface wear model for hypoid gear pairs. Wear 2009, 267, 1595–1604. [Google Scholar] [CrossRef]
  5. Advanced Numerical Solutions. Transmission3D User’s Manual v. 2.54. 2020. Available online: http://ansol.us/Products/TX3 (accessed on 1 January 2021).
  6. 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]
  7. Kolivand, M.; Kahraman, A. A load distribution model for hypoid gears using ease-off topography and shell theory. Mech. Mach. Theory 2009, 44, 1848–1865. [Google Scholar] [CrossRef]
  8. 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]
  9. Huang, D.; Wang, Z.; Kubo, A. Hypoid gear integrated wear model and experimental verification design and test. Int. J. Mech. Sci. 2020, 166, 105228. [Google Scholar] [CrossRef]
  10. Grabovic, E.; Ciulli, E.; Artoni, A.; Gabiccini, M. A model for the prediction of frictional power losses in hypoid gears. In Proceedings of the 4th International Conference of IFToMM, Italy, Napoli, 7–9 September 2022. in press. [Google Scholar]
  11. Cao, W.; Pu, W.; Wang, J.; Xiao, K. Effect of contact path on the mixed lubrication performance, friction and contact fatigue in spiral bevel gears. Tribol. Int. 2018, 123, 359–371. [Google Scholar] [CrossRef]
  12. Mohammadpour, M.; Theodossiades, S.; Rahnejat, H. Elastohydrodynamic lubrication of hypoid gear pairs at high loads. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2012, 226, 183–198. [Google Scholar] [CrossRef] [Green Version]
  13. Pu, W.; Wang, J.; Yang, R.; Zhu, D. Mixed Elastohydrodynamic Lubrication with Three-Dimensional Machined Roughness in Spiral Bevel and Hypoid Gears. J. Tribol. 2015, 137. [Google Scholar] [CrossRef]
  14. Castro, J.; Seabra, J. Coefficient of friction in mixed film lubrication: Gears versus twin-discs. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2007, 221, 399–411. [Google Scholar] [CrossRef]
  15. Gold, P.W.; Schmidt, A.; Dicke, H.; Loos, J.; Assmann, C. Viscosity–pressure–temperature behaviour of mineral and synthetic oils. J. Synth. Lubr. 2001, 18, 51–79. [Google Scholar] [CrossRef]
  16. Hamrock, B. Fundamentals of Fluid Film Lubrication; McGraw-Hill: New York, NY, USA, 1994; pp. 499–505. [Google Scholar]
  17. Bassani, R.; Ciulli, E.; Guiggiani, M.; Piccigallo, B.; Andrei, G. Elastohydrodynamic lubrication in face gears drives. In Proceedings of the International Tribology Conference, Nagasaki, Japan, 29 October–2 November 2000. [Google Scholar]
  18. Castro, J.; Seabra, J. Scuffing and lubricant film breakdown in FZG gears Part I. Analytical and experimental approach. Wear 1998, 215, 104–113. [Google Scholar] [CrossRef]
  19. Hoppert, M. Analytische und Experimentelle Untersuchungen zum Wirkungsgradverhalten von Achsgetrieben. Ph.D. Thesis, Fakultät für Maschinenbau der Technischen Universität Ilmenau, Ilmenau, Germany, 2015. [Google Scholar]
  20. Grabovic, E.; Artoni, A.; Gabiccini, M. Holistic Optimal Design Of Face-Milled Hypoid Gearsets. In Proceedings of the IDETC-CIE 2022, St. Louis, MO, USA, 14–17 August 2022. number IDETC2022-89598. [Google Scholar]
  21. Ciulli, E. Formulas for Entraining Velocity in Lubricated Line Contacts. J. Tribol. 2002, 124, 856–858. [Google Scholar] [CrossRef]
  22. Gohar, R. Elastohydrodynamics; E. Horwood Halsted Press: Chichester, UK; New York, NY, USA, 1988.
  23. Bajpai, P.; Kahraman, A.; Anderson, N.E. A Surface Wear Prediction Methodology for Parallel-Axis Gear Pairs. J. Tribol. 2004, 126, 597–605. [Google Scholar] [CrossRef]
  24. Piegl, L.; Tiller, W. The NURBS Book, 2nd ed.; Springer: New York, NY, USA, 1996. [Google Scholar]
  25. Mattei, L.; Di Puccio, F. Influence of the wear partition factor on wear evolution modelling of sliding surfaces. Int. J. Mech. Sci. 2015, 99, 72–88. [Google Scholar] [CrossRef]
Figure 1. Hypoid gearset. The dashed line represents the instantaneous screw axis of the relative motion.
Figure 1. Hypoid gearset. The dashed line represents the instantaneous screw axis of the relative motion.
Machines 10 00748 g001
Figure 2. Load shared by the contacting teeth over a mesh cycle.
Figure 2. Load shared by the contacting teeth over a mesh cycle.
Machines 10 00748 g002
Figure 3. Entraining velocity u of the lubricant over a mesh cycle.
Figure 3. Entraining velocity u of the lubricant over a mesh cycle.
Machines 10 00748 g003
Figure 4. Entraining ( R e ) and side leakage ( R s ) radii over a mesh cycle.
Figure 4. Entraining ( R e ) and side leakage ( R s ) radii over a mesh cycle.
Machines 10 00748 g004
Figure 5. Central film thickness ( h c ) values over a mesh cycle.
Figure 5. Central film thickness ( h c ) values over a mesh cycle.
Machines 10 00748 g005
Figure 6. Lambda factor ( Λ ) over a mesh cycle.
Figure 6. Lambda factor ( Λ ) over a mesh cycle.
Machines 10 00748 g006
Figure 7. Coefficient of friction (f) over a mesh cycle.
Figure 7. Coefficient of friction (f) over a mesh cycle.
Machines 10 00748 g007
Figure 8. Efficiency map as a function of pinion torque and pinion speed.
Figure 8. Efficiency map as a function of pinion torque and pinion speed.
Machines 10 00748 g008
Figure 9. Initial contact patterns. (a) Pinion. (b) Gear.
Figure 9. Initial contact patterns. (a) Pinion. (b) Gear.
Machines 10 00748 g009
Figure 10. Pinion patterns after 30 M cycles. Only one geometry update is performed. (a) Contact pattern. (b) Wear pattern.
Figure 10. Pinion patterns after 30 M cycles. Only one geometry update is performed. (a) Contact pattern. (b) Wear pattern.
Machines 10 00748 g010
Figure 11. Gear patterns after 30 Mc (pinion revolutions) cycles. Only one geometry update is performed. (a) Contact pattern. (b) Wear pattern.
Figure 11. Gear patterns after 30 Mc (pinion revolutions) cycles. Only one geometry update is performed. (a) Contact pattern. (b) Wear pattern.
Machines 10 00748 g011
Figure 12. Evolution of the pinion contact pattern with a geometry update every 5 Mc. (a) 5 Mc. (b) 10 Mc. (c) 15 Mc. (d) 20 Mc. (e) 25 Mc. (f) 30 Mc.
Figure 12. Evolution of the pinion contact pattern with a geometry update every 5 Mc. (a) 5 Mc. (b) 10 Mc. (c) 15 Mc. (d) 20 Mc. (e) 25 Mc. (f) 30 Mc.
Machines 10 00748 g012
Figure 13. Evolution of the pinion (cumulative) wear pattern with a geometry update every 5 Mc. (a) 5 Mc. (b) 10 Mc. (c) 15 Mc. (d) 20 Mc. (e) 25 Mc. (f) 30 Mc.
Figure 13. Evolution of the pinion (cumulative) wear pattern with a geometry update every 5 Mc. (a) 5 Mc. (b) 10 Mc. (c) 15 Mc. (d) 20 Mc. (e) 25 Mc. (f) 30 Mc.
Machines 10 00748 g013
Figure 14. Evolution of the gear contact pattern with a geometry update every 5 Mc. (a) 5 Mc. (b) 10 Mc. (c) 15 Mc. (d) 20 Mc. (e) 25 Mc. (f) 30 Mc.
Figure 14. Evolution of the gear contact pattern with a geometry update every 5 Mc. (a) 5 Mc. (b) 10 Mc. (c) 15 Mc. (d) 20 Mc. (e) 25 Mc. (f) 30 Mc.
Machines 10 00748 g014
Figure 15. Evolution of the gear (cumulative) wear pattern with a geometry update every 5 Mc. (a) 5 Mc. (b) 10 Mc. (c) 15 Mc. (d) 20 Mc. (e) 25 Mc. (f) 30 Mc.
Figure 15. Evolution of the gear (cumulative) wear pattern with a geometry update every 5 Mc. (a) 5 Mc. (b) 10 Mc. (c) 15 Mc. (d) 20 Mc. (e) 25 Mc. (f) 30 Mc.
Machines 10 00748 g015
Figure 16. Instantaneous contact area(s) at the same meshing configuration on the pinion before and after wear simulation (performed updating the geometry every 5 Mc). (a) Before wear. (b) After wear.
Figure 16. Instantaneous contact area(s) at the same meshing configuration on the pinion before and after wear simulation (performed updating the geometry every 5 Mc). (a) Before wear. (b) After wear.
Machines 10 00748 g016
Figure 17. Instantaneous contact area(s) at the same meshing configuration on the gear before and after wear simulation (performed updating the geometry every 5 Mc). (a) Before wear. (b) After wear.
Figure 17. Instantaneous contact area(s) at the same meshing configuration on the gear before and after wear simulation (performed updating the geometry every 5 Mc). (a) Before wear. (b) After wear.
Machines 10 00748 g017
Figure 18. Comparison of the final pinion contact patterns resulting from wear simulations with different numbers of geometry updates. (a) Update after 30 Mc. (b) Update every 10 Mc. (c) Update every 5 Mc. (d) Update every 1 Mc. (e) Update every 0.5 Mc.
Figure 18. Comparison of the final pinion contact patterns resulting from wear simulations with different numbers of geometry updates. (a) Update after 30 Mc. (b) Update every 10 Mc. (c) Update every 5 Mc. (d) Update every 1 Mc. (e) Update every 0.5 Mc.
Machines 10 00748 g018
Figure 19. Maximum contact pressure from wear simulations with different geometry update frequencies.
Figure 19. Maximum contact pressure from wear simulations with different geometry update frequencies.
Machines 10 00748 g019
Figure 20. Maximum pinion wear depth from wear simulations with different geometry update frequencies.
Figure 20. Maximum pinion wear depth from wear simulations with different geometry update frequencies.
Machines 10 00748 g020
Table 1. Coefficients of the 75W90 oil.
Table 1. Coefficients of the 75W90 oil.
SymbolValue
β 2.86 · 10 2 [ K 1 ]
Thermal conductivityk 0.14 [ W mK ]
Eyring shear stress τ E 8 [ MPa ]
Vogel coefficients
η 0 5.90 · 10 5 [ Pa s ]
B1205.5 [ C ]
C 125.3 [ C ]
Pressure-viscosity coefficients
a 1 4.929 · 10 2 [bar]
a 2 3.901 [ bar C ]
b 1 2.479 · 10 2
b 2 2.354 · 10 4 [ 1 C ]
Table 2. Hypoid basic parameters.
Table 2. Hypoid basic parameters.
ParameterValue
Hypoid offset30 mm
Shaft angle 90 deg
Hand (pinion)Left
Gear typeFormate
Pinion teeth13
Gear teeth44
Spiral angle (pinion)45.5 deg
Face width (gear)31.87 mm
Pitch diameter (gear)192.36 mm
Cutter mean radius63.5 mm
Nominal pressure angle20 deg
Profile shift coefficient0.524
Profile crowning (pinion)60 μ m
Lengthwise crowning (pinion)100 μ m
Surface roughness0.45 μ m
Operating oil temperature100 C
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Grabovic, E.; Artoni, A.; Gabiccini, M.; Guiggiani, M.; Mattei, L.; Di Puccio, F.; Ciulli, E. Friction-Induced Efficiency Losses and Wear Evolution in Hypoid Gears. Machines 2022, 10, 748. https://doi.org/10.3390/machines10090748

AMA Style

Grabovic E, Artoni A, Gabiccini M, Guiggiani M, Mattei L, Di Puccio F, Ciulli E. Friction-Induced Efficiency Losses and Wear Evolution in Hypoid Gears. Machines. 2022; 10(9):748. https://doi.org/10.3390/machines10090748

Chicago/Turabian Style

Grabovic, Eugeniu, Alessio Artoni, Marco Gabiccini, Massimo Guiggiani, Lorenza Mattei, Francesca Di Puccio, and Enrico Ciulli. 2022. "Friction-Induced Efficiency Losses and Wear Evolution in Hypoid Gears" Machines 10, no. 9: 748. https://doi.org/10.3390/machines10090748

APA Style

Grabovic, E., Artoni, A., Gabiccini, M., Guiggiani, M., Mattei, L., Di Puccio, F., & Ciulli, E. (2022). Friction-Induced Efficiency Losses and Wear Evolution in Hypoid Gears. Machines, 10(9), 748. https://doi.org/10.3390/machines10090748

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