Prerpints.org logo
Preprint
Article

A Computational Model for Analysing the Wear Behaviour of Polymer Gears

Altmetrics

Downloads

179

Views

131

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

23 January 2024

Posted:

24 January 2024

You are already at the latest version

Alerts
Abstract
This study presents the computational model to determine the wear behaviour of polymer gears. Using PrePoMax finite element numerical calculation software, a proposed computational model was built to predict wear behaviour based on Archard's wear model. This allows the calculation of the wear depth in each loading cycle with constant mesh updating. The developed computational model has been evaluated on the spur gear pair, where the pinion made of POM has meshed with a support gear made of steel. The computational results that were obtained were compared with the analytical results according to the VDI 2736 guidelines. Based on this comparison, it was concluded that the proposed computational model could be used to simulate the wear behaviour of contacting mechanical elements like gears, bearings, etc. The main advantage of the model, if compared to the standardised procedure according to the VDI 2736 guidelines, is the geometry updating after a chosen number of loading cycles, which enables a more accurate prediction of wear behaviour under rolling/sliding loading conditions.
Keywords: 
Subject: Engineering  -   Mechanical Engineering

1. Introduction

Polymer gears are used widely in many engineering applications, such as office appliances, mechatronic devices, household facilities, computer and laboratory equipment, medical instruments, etc. [1,2,3,4,5,6,7]. In these applications, the contact of two mechanical elements (i.e. gear flanks) is concentrated on a small contact area, resulting in higher contact stresses, which often exceed the yield stress of the material. Due to repeated contact loading during the exploitation, the local cyclic plastic deformations may lead to surface damage, usually termed surface wear [8,9,10,11].
In recent years, many wear models relating to rolling/sliding contacting mechanical elements were proposed. However, Archard’s model [12] is widely used in engineering applications. The model requires information about the coefficient of wear, relative sliding distance and load and has been applied to gear wear prediction by many researchers, who have derived equations that connect the wear volume or wear depth with the operating conditions and material properties of the gear tooth. Flodin et al. [13,14] derived a load distribution equation to determine the transferred load between the spur and helical gear pairs, using Herztzian contact theory to calculate the contact stress and the single point tracking method for relative sliding distance. Later, this model was extended to analyse the wear behaviour of double helical gears [15], worm gears [16], internal gears [17], hypoid bevel gears [18,19,20] and planetary gears [21,22].
Discretising the Archard model allows for a more detailed analysis of the wear behaviour of materials, especially where wear varies across the material surface. Põdra et al. [23] used the finite element method for simulating sliding wear between the rotation disk and pin. For sliding-rolling applications, such as gears, Lundvall et al. [24,25] derived time increment equations to predict wear using the nonsmooth Newton method. Employing the finite element method opens more possibilities to improve accuracy, which is why Bajpai et al. [26] combined the finite element method and the Archard model, using the iterative method to consider contact changes due to wear. Because mechanical elements like gears, bearings etc., are dynamically loaded components, it is also useful to consider time-varying sliding distance and time-varying contact load in the wear analysis to make results more accurate [27,28,29,30,31,32].
In recent years, polymer gears have become more widely used due to their good properties such as lubrication-free operation, corrosion resistance, lightweight, arbitrary geometry and low cost of manufacture. They are mostly found in gears within small appliances, where they can be combined with a variety of other materials [33]. As they are under cyclic load, they are dynamically loaded components. As the gears roll and slide between the tooth flanks of the gears, erosion and wear occur on the surfaces over time, which can be determined with experiments, analytically or numerically. The analytical calculation of polymer gears can be found in the VDI 2736 standard (Association of German Engineers) [34]. This covers the basic calculation and the presentation of the properties of polymer gears. Among those things, the wear depth of the gears can be calculated using this standard. Although this standard can be used to obtain certain values, it is limited to constant operating conditions, which, unlike gears from metallic materials, can change rapidly due to their poorer mechanical properties. As mentioned, the finite element method provides more accurate results and allows rapid change of material properties and operating conditions. Lin et al. [35] investigated the interaction between contact loads and tooth wear in POM (poly-oxy-methylene) and nylon tooth pairs, using the finite element method to analyse tooth stiffness and a quasi-linear iterative method to calculate wear depth. The result shows that peak load increases as the number of cycles reaches a threshold limit.
In the framework of the presented study, a computational model was built using open-source finite element numerical calculation software PrePoMax to predict wear behaviour based on Archard’s wear model. This allows the calculation of the wear depth in each loading cycle with constant mesh updating. The computational results that are obtained are compared with the analytical results according to the VDI 2736 guidelines.

2. Materials and Methods

2.1. Computational modelling

In this study, the PrePoMax open-source pre- and post-processor for the finite element method (FEM) [36] was used to analyse the wear behaviour of rolling/sliding contacting mechanical elements (i.e. mating gear flanks). The solver used by PrePoMax is the open-source CalculiX FEM solver [37] that does not support the wear analysis by itself, so the wear analysis was implemented into the PrePoMax. The linear Archard wear model [38] for predicting the volume of worn material was chosen for this purpose. The model was originally developed based on experimental results of cyclic wear of metal specimens under dry contact conditions. The model is represented by the following equation:
V = k · F H · s
where V is the volume of worn material, k is the dimensionless wear coefficient, F is the normal contact force, H is the surface hardness in contact, and s represents the sliding distance. To account for changes in the geometry of the contact surfaces during wear cycles, the aforementioned model, which neglects the geometry change, can be divided into multiple steps when used in numerical analyses using the FEM. In a multi-step numerical procedure, where the depth of wear needs to be calculated for each wear cycle [39], Eq. (1) can be expressed as follows:
h i = ∫ k H · p i · d s
where h i represents the depth of wear in the i-th wear cycle, p i denotes the normal contact pressure, and d s represents the differential sliding distance. The total wear depth h over n wear cycles is then computed by the Eq. (3), summing up the wear depths of individual cycles. Such an approach enables the wear analysis to be independent of the model geometry, boundary conditions, and loading.
h = ∑ i = 1 n h i
The multi-step method implemented in PrePoMax is based on the following assumptions: (i) the entire wear process can be divided into multiple repetitive wear cycles, (ii) wear occurs only on the softer surface in the contact pair (slave surface), (iii) the change in contact surface geometry during one wear cycle is negligible and (iv) the dimensionless wear coefficient is constant.
For each wear cycle, a finite element model specific to that cycle needs to be prepared. The model must accurately capture the evolving contact conditions during the cycle, requiring a nonlinear contact analysis solved in multiple increments. For an accurate result, the number of increments must be large enough to result in a smooth movement of the contact pressure over the finite element mesh of the contact surfaces.
The wear calculation for the cycle begins after the nonlinear contact analysis is completed. All wear parameters are computed only at those slave surface finite element nodes that were in active contact with the master surface during the wear cycle. The normal contact pressure p i , j at the node j and the nodal relative sliding distance ∆ s i , j are determined from the resulting contact fields. The wear depth at node h i , j is then computed using the expression:
h i , j = k j H j · p i , j · ∆ s i , j
where the material properties kj and Hj are taken at node j. Processing of each wear cycle concludes with the calculation of the change in the contact surface geometry due to wear, determined by the wear displacements at the nodes calculated from the wear depth. The node wear displacements due to wear in each cycle h ⃑ i , j are calculated by a scalar product:
h ⃑ i , j = − h i , j · n ⃑ j
where where n ⃑ j represents the node’s outward normal vector. The node’s outward normal vector is determined as an average value of the normal vectors of the FEs sharing the node j (Figure 1). At the start of the next wear cycle, the wear displacements computed in the previous cycle are used to update the mesh. Thus, the simulation of the new wear cycle is carried out on a mesh where the wear from the previous cycles is considered
Applying the wear displacements to the mesh changes only the position of the surface nodes and thus deforms the shape of the finite elements, making all surface elements on the wear surface a little thinner. If the accumulated wear displacement is large enough, this might lead to unusable finite element shapes. This is especially problematic when a dense mesh is used on the contact surfaces to obtain an accurate contact pressure field since even a small wear displacement will deform a small finite element beyond being usable. To overcome this problem, local or global remeshing of the domain can be applied, or a method for mesh redistribution can be used. In PrePoMax software, a boundary displacement method (BDM) [40] for mesh redistribution was introduced. This is a two-step procedure where the surface nodes are first moved, and then new positions of the internal nodes are computed. From many different approaches to the BDM, the positions of the internal nodes were determined based on the solution of the linear elastic solid body problem. This was the most convenient approach since it can be solved with the FEM method on the existing mesh. The method was integrated into PrePoMax by an additional BDM simulation step that is added after each wear cycle. In this simulation step, the mesh from the previous wear step is taken, uniform elastic properties are assigned to all elements, and a prescribed boundary condition is applied to all nodes on the surface. In the nodes where the wear occurs, the value of the prescribed displacement vector is equal to the wear displacement h ⃑ i , j and equal to 0 ⃑ for all other nodes. After the BDM simulation step completes, the computed nodal displacements for all nodes from the BDM step are used to deform the mesh before the next wear step is started.
Running multiple wear steps one after the other and deforming the mesh after each step can lead to the occurrence of high-frequency oscillations (roughness) of the mesh (see Figure 1), which results in an irregular contact pressure field. This causes an irregular wear displacement field, which additionally enhances the roughness of the contact surfaces and renders the results of the wear analysis unusable. A mesh smoothing algorithm on the wear displacement field was introduced to remove the high-frequency oscillations from the mesh. The smoothing algorithm is based on the Laplacian mesh smoothing procedure [41], where the new position of the node x ⃑ j is computed from the positions of its neighbouring nodes x ⃑ k as:
x ⃑ j = 1 m ∑ k = 1 m x ⃑ k
where m equals the number of node neighbours. Mesh smoothing is applied after the wear displacement field is computed and before the BDM step is prepared. The effect of the smoothing step can be increased by running the smoothing algorithm multiple times, which can be defined by the user. All these features make wear analyses applicable and robust for a large variety of problems, and the user can prepare the wear analysis without additional coding using only the built-in user interface. This makes the implementation usable for many engineers since all these features are already available in the latest release version of the PrePoMax v1.4.0, available online [42].

2.2. Geometrical and material parameters of the analysed gear pair

The developed computational model has been evaluated on the spur gear pair, where the pinion made of POM has meshed with a support gear made of steel. The basic geometrical and material parameters of the analysed gear pair are presented in Table 1.

2.3. FE model of gear pair

The geometry of the gear pair is shown in Figure 2, where the pinion is shown in grey, and the support gear is shown in green. Here, only one tooth of the gear and three teeth of the pinion are presented. It was due to the smaller stiffness of the pinion material (POM) if compared to the gear material (steel). The analysis of one load cycle was split into three simulation steps, as shown in Figure 3: static step (contact initialisation), wear step (computation of wear) and BDM step (mesh redistribution). In the simplified model presented in Figure 2, the reference points RP1 and RP2 have been defined in both axes of the gear pair. Here, RP1 was connected rigidly to the gear shaft, with the surfaces of the gear marked in red. Additionally, at this reference point, a torque of T = 1.6 Nm was applied acting to the pinion tooth, and a boundary condition was imposed, which allowed the pinion to move only about its axis. The point RP2 was connected rigidly to the pinion surfaces marked in blue. At this reference point, a boundary condition was prescribed first in the static step, which prevented displacement and rotation in all directions of the coordinate system. As sliding was required for the wear calculation procedure, the boundary condition in RP1 in the second step (Wear step) required a rotation around the pinion axis of ρ = 0.285 rad. The angle of rotation was prescribed in such a way that the gear only reached the position where the gear started to return to the already worn surface.
The computational model implemented into the PrePoMax software (see Section 2.1) was then used to analyse the wear behaviour of mating gear flanks. The numerical simulation of gear meshing consisted of three separate simulation steps, which followed one after the other. This simulation then repeats to account for the selected number of load cycles. Simulating each load cycle in multi-thousand load cycle analysis would require too much computational effort. At the same time, during one load cycle, the geometry of the model does not change significantly; thus, multiple load cycles can be merged into one simulation cycle that updates the mesh.
Archard’s model implemented in the PrePoMax software is limited to a 3D wear calculation, as it calculates the volume of wear. However, as the pinion and gear analysed in this study do not vary in thickness, it would be reasonable to simplify the model to 2D, but this was not possible. Therefore, a reduced thickness of 0.25 mm was assumed to reduce the computational time, and consequently, the load (torque) has also been reduced proportionally. Furthermore, an additional boundary condition has been added to both gears to prevent the gears from moving in the thickness direction. To best approximate the 2D model, a finite element mesh with prismatic finite elements was created using the boundary layer function, as shown in Figure 4. The overall model was meshed with a global mesh size of 1 mm. However, for better accuracy, the model around the contact was meshed with a local mesh size of 0.05 mm, which resulted in 3300 parabolic prismatic finite elements. The model used a linear overclosure-pressure contact relation with stiffness K = 105 N/mm3 and coefficient of friction μ = 0.25 mm.
Two computational models have been considered in the numerical simulation. In the first model (Model 1), a wear calculation after 10,000 load cycles was determined. First, for the purpose of comparison with VDI 2736, the 10,000 load cycles were calculated in one simulation cycle, which means that the geometry was not updated during the calculation. The purpose of this computational model was to compare the results with the VDI 2736 calculation, which does not consider the geometry change. In the second model (Model 2), 10,000 load cycles were used again. However, if compared to the first model, the number of cycles was divided into ten simulation cycles, meaning that after every 1,000 load cycles, the geometry was updated. The wear results were thus transferred directly to the finite element mesh, and the wear was recalculated on the updated geometry.

2.4. Analytical approach according to VDI 2736

Following the calculation according to VDI 2736, a wear behaviour can be obtained for the tooth flank profile of the polymer gear. The calculation is position-dependent, as the wear is different at each point due to the involute shape of the tooth. Here, the wear starts to decrease from the tip of the tooth up to the pitch point. Due to the involute shape of the tooth, sliding dominates at the starting engagement point, then the tooth rolls at the pitch circle, and later, towards the root of the tooth, sliding increases gradually again, and rolling decreases. Therefore, in the pitch point region, it can be observed that the wear is theoretically zero. Considering this fact, the local wear of the gear flank can be obtained as follows:
W l o c a l = F n b w ⋅ N L ⋅ ζ ⋅ k w
where bw is a common face width, NL is the number of load cycles, ζ is the specific local sliding of gear flanks, kw is the wear coefficient, and Fn is the normal force acting on the tooth flank
F n = 2 ⋅ T d d w ⋅ cos ⁡ α t
where Td is the nominal torque, dw is the pitch circle diameter, and αt is the pressure angle in the transverse section. It follows from Eq. (8) that the normal force Fn is a constant throughout the engagement line of gear flanks. However, its tangential and radial components change, resulting in a combination of sliding and rolling of gear flanks. When sliding is involved, a specific local sliding ζ appears, which should be considered when determining the local wear using Eg. (7). Figure 5 shows the sliding conditions in an arbitrary mating point Y before pitch point C. The circular velocities in an arbitrary mating point Y (vy1= ry1⋅ω1 and vy2= ry2⋅ω2), which are acting perpendicular to the lines O1Y and O2Y, can be distributed in the directions of common normal (vny1 and vny2) and common tangent (vty1 and vty2). Here, the normal velocity components must be the same, i.e. vny1 = vny2. Furthermore, the tangential components vty1 and vty2 can be obtained from the triangles O1T1Y and O2T2Y and the hatched triangles.
The sliding velocities of gear flanks in an arbitrary mating point Y (vgy1 for a pinion and vgy2 for a gear) are defined as a difference between the tangential velocity components vty1 and vty2 (vgy1 = vty1 − vty2; vgy2 = vty2 − vty1). From triangles O1T1C and O2T2C and with consideration O1C = dw1/2 and O2C = dw2/2, it follows:
v gy 1 = −   ω 1 + ω 2 ⋅ YC ¯   ;   v gy 2 = +   ω 1 + ω 2 ⋅ YC ¯
Because the sum (ω1 + ω2) is a constant value, it is clear from Eq. (9) that the sliding velocity vgY only depends on the distance of point Y from the pitch point C. When an arbitrary mating point Y coincides with pitch point C, the sliding velocity equals zero.
Based on the sliding velocity vgy1 in an arbitrary mating point Y, the specific local sliding of gear flanks ζ is defined as [44]:
ζ = v gy 1 v ty 1 = 1 + z 1 z 2 Y C ¯ r w 1 · s i n α t + Y C ¯
Following the calculation procedure according to the VDI 2736, the wear behaviour of gear flanks is position-dependent due to different sliding velocities along the engagement line of mating gears. It is clear that the sliding velocities and, consequently, the wear of gear flanks are the highest at the starting and ending engagement points and then decrease toward pitch point C.

3. Results and discussion

Figure 6 shows the wear of the analysed gear made of POM after 10,000 loading cycles considering Model 1, where the geometry change due to wear of the gear flank was not considered in the numerical simulation. It is clear that the depth of the wear is the highest at the starting and ending engagement points, which corresponds to the theoretical background as described in Section 2.4. Similar conclusions can also be made for Model 2, which considered the geometry change during the numerical simulations (see Figure 7).
Figure 8 shows the comparison between the initial and worn geometry of the analysed POM gear after 10,000 loading cycles (a scaling factor of 3 was considered). It is clear that the main differences appear at the starting and ending mating points, where the wear of the gear flank is the highest.
Figure 9 shows the comparison between computational results (for Models 1 and 2) and the standardised approach according to VDI 2736. It should be noted that the VDI 2736 is based on the wear properties of the materials obtained from a rotating pin-bearing disc test, whereas the final equation for the calculation of the wear depth is based solely on the involute profile of the analysed gear. Based on the obtained computational results using the proposed computational model and analytical results according to VDI 2736, it can be concluded that both approaches give a similar wear behaviour of the analysed POM gear (the maximum wear appears at the starting and ending mating points). It is also evident that the computational results according to Model 1 (geometry changes due to wear of the gear flank were not considered) better describe the wear behaviour if compared to VDI 2736. However, it should be pointed out that, in real applications of the gear drive, the geometry of gear flanks constantly changes due to wear, which may affect the subsequent wear phenomenon of gear flanks significantly. By adding sensible updates to the geometry, we can get closer to more realistic wear values. Furthermore, it is possible to observe the change in wear compared to the model without a geometry update. In Figure 7, the green line shows the path of the wear curve of the Model 2. In the computational model, oscillations of the wear depth values between adjacent nodes have occurred. Wear is dependent directly on the contact pressure field, which means that if the contact pressure field is smooth, then the calculated wear field is smooth, and vice versa. Therefore, a smoothing function was applied to the contact pressure field, which acts as an average filter, averaging the results between adjacent nodes. For that reason, the obtained wear curves are corrected with dotted lines in Figure 9. The values at pitch point C are slightly higher, but this is not only due to smoothing but also due to the new tooth geometry, which is no longer ideal because of the tooth deformation. This results in sliding instead of rolling. It is possible to observe a marked change at the root of the tooth. The wear starts to decrease with the number of cycles also due to the change in geometry. Relief apparently occurs at the root of the tooth as the maximum pressure is reduced or redistributed.

4. Conclusions

In the framework of this study, the comprehensive computational model to analyse the wear behaviour of polymer gears was developed and evaluated on the spur gear pair, where the pinion made of POM has meshed with a support gear made of steel. A proposed computational model was built using PrePoMax finite element numerical calculation software considering Archard’s wear model to predict the wear behaviour of analysed contacting surfaces. Based on the obtained computational results and their comparison to the analytical results according to the VDI 2736 guideline, the following conclusions could be made:
  • The proposed computational model could be used to simulate the wear behaviour of contacting mechanical elements like gears, bearings, etc.
  • The use of the boundary displacement method (BDM) in the framework of PrePoMax software provides information on the geometry change that alters the operating conditions in subsequent load cycles.
  • Due to the contact problem and the formation of contact pressure peaks in the finite element mesh, it is necessary to consider mesh smoothing in the model to allow a smooth distribution of wear over the surface. This avoids additional convergence problems in the use of BDM, but it does result in the averaging of values in locations where the differences in wear between adjacent element nodes should actually be larger (such as pitch point) and is not entirely correct. The validity of the results should, therefore, be checked with a model without smoothing.
  • With the computational model using a multi-step geometry update, more accurate results were obtained, which show a reduced depth of wear at the root of the tooth. In the pitch point region, a non-zero value of the wear depth appears as the number of cycles increases. The tooth deflection and the new tooth flank geometry have a major impact on this.
  • The main advantage of the model, if compared to the standardised procedure according to the VDI 2736 guidelines, is the geometry updating after a certain number of loading cycles, which enables a more accurate prediction of wear behaviour under rolling/sliding loading conditions.
  • In the future, a comprehensive 3D numerical model will be developed to analyse the interaction of meshing gears. This model aims to compare the results obtained from experimental testing with those derived from numerical simulations.
  • Potential future advancements of the existing Archard’s wear model could include incorporating a contact stress-dependent wear coefficient, as demonstrated in prior research. Currently, the wear coefficient remains constant throughout the analysis. Additionally, the model could be enhanced by establishing a relationship between the wear parameter and wear depth in order to account for diverse surface improvement techniques.

Author Contributions

Conceptualisation, S.G.; methodology, S.G. and A.I.; software, A.I. and M.B.; validation, S.G.; formal analysis, A.I. and M.B.; investigation, A.I.; writing—original draft preparation, A.I and S.G.; writing—review and editing, all authors; visualisation, A.I.; supervision, S.G.; project administration, S.G.; funding acquisition, S.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovenian Research Agency (ARIS), Research Core Funding No. P2-0063.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors acknowledge the financial support of the Research Core Funding (No. P2-0063) from the Slovenian Research Agency.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pogačnik, A.; Tavčar, J. An Accelerated Multilevel Test and Design Procedure for Polymer Gears. Mater. Des. 2015, 65, 961–973. [Google Scholar] [CrossRef]
  2. Matkovič, S.; Pogačnik, A.; Kalin, M. Wear-Coefficient Analyses for Polymer-Gear Life-Time Predictions: A Critical Appraisal of Methodologies. Wear 2021, 480–481, 203944. [Google Scholar] [CrossRef]
  3. Muratović, E.; Muminović, A.; Delić, M.; Pervan, N.; Muminović, A.J.; Šarić, I. Potential and Design Parameters of Polyvinylidene Fluoride in Gear Applications. Polymers (Basel). 2023, 15, 4275. [Google Scholar] [CrossRef] [PubMed]
  4. Tosangthum, N.; Krataitong, R.; Wila, P.; Koiprasert, H.; Buncham, K.; Kansuwan, P.; Manonukul, A.; Sheppard, P. Dry Rolling-Sliding Wear Behavior of ER9 Wheel and R260 Rail Couple under Different Operating Conditions. Wear 2023, 518–519, 204636. [Google Scholar] [CrossRef]
  5. Hu, Y.; Guo, L.C.; Maiorino, M.; Liu, J.P.; Ding, H.H.; Lewis, R.; Meli, E.; Rindi, A.; Liu, Q.Y.; Wang, W.J. Comparison of Wear and Rolling Contact Fatigue Behaviours of Bainitic and Pearlitic Rails under Various Rolling-Sliding Conditions. Wear 2020, 460–461, 203455. [Google Scholar] [CrossRef]
  6. Zhang, B.; Liu, H.; Zhu, C.; Ge, Y. Simulation of the Fatigue-Wear Coupling Mechanism of an Aviation Gear. Friction 2021, 9, 1616–1634. [Google Scholar] [CrossRef]
  7. Zhu, W.T.; Guo, L.C.; Shi, L.B.; Cai, Z.B.; Li, Q.L.; Liu, Q.Y.; Wang, W.J. Wear and Damage Transitions of Two Kinds of Wheel Materials in the Rolling-Sliding Contact. Wear 2018, 398–399, 79–89. [Google Scholar] [CrossRef]
  8. Pei, X.; Pu, W.; Zhang, Y.; Huang, L. Surface Topography and Friction Coefficient Evolution during Sliding Wear in a Mixed Lubricated Rolling-Sliding Contact. Tribol. Int. 2019, 137, 303–312. [Google Scholar] [CrossRef]
  9. Seo, J.-W.; Jun, H.-K.; Kwon, S.-J.; Lee, D.-H. Rolling Contact Fatigue and Wear of Two Different Rail Steels under Rolling–Sliding Contact. Int. J. Fatigue 2016, 83, 184–194. [Google Scholar] [CrossRef]
  10. Xu, L.; Fan, X.; Wei, S.; Liu, D.; Zhou, H.; Zhang, G.; Zhou, Y. Microstructure and Wear Properties of High-Speed Steel with High Molybdenum Content under Rolling-Sliding Wear. Tribol. Int. 2017, 116, 39–46. [Google Scholar] [CrossRef]
  11. Ramalho, A.; Esteves, M.; Marta, P. Friction and Wear Behaviour of Rolling–Sliding Steel Contacts. Wear 2013, 302, 1468–1480. [Google Scholar] [CrossRef]
  12. Archard, J.F. Contact and Rubbing of Flat Surfaces. J. Appl. Phys. 1953, 24, 981–988. [Google Scholar] [CrossRef]
  13. Flodin, A.; Andersson, S. Simulation of Mild Wear in Spur Gears. Wear 1997, 207, 16–23. [Google Scholar] [CrossRef]
  14. Flodin, A.; Andersson, S. Simulation of Mild Wear in Helical Gears. Wear 2000, 241, 123–128. [Google Scholar] [CrossRef]
  15. Zhou, C.; Wang, H. An Adhesive Wear Prediction Method for Double Helical Gears Based on Enhanced Coordinate Transformation and Generalized Sliding Distance Model. Mech. Mach. Theory 2018, 128, 58–83. [Google Scholar] [CrossRef]
  16. Jbily, D.; Guingand, M.; de Vaujany, J.-P. A Wear Model for Worm Gear. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2016, 230, 1290–1302. [Google Scholar] [CrossRef]
  17. Tunalioğlu, M.Ş.; Tuç, B. Theoretical and Experimental Investigation of Wear in Internal Gears. Wear 2014, 309, 208–215. [Google Scholar] [CrossRef]
  18. 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]
  19. 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]
  20. 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]
  21. Shen, Z.; Qiao, B.; Yang, L.; Luo, W.; Yang, Z.; Chen, X. Fault Mechanism and Dynamic Modeling of Planetary Gear with Gear Wear. Mech. Mach. Theory 2021, 155, 104098. [Google Scholar] [CrossRef]
  22. Shen, Z.; Qiao, B.; Yang, L.; Luo, W.; Chen, X. Evaluating the Influence of Tooth Surface Wear on TVMS of Planetary Gear Set. Mech. Mach. Theory 2019, 136, 206–223. [Google Scholar] [CrossRef]
  23. Põdra, P.; Andersson, S. Simulating Sliding Wear with Finite Element Method. Tribol. Int. 1999, 32, 71–81. [Google Scholar] [CrossRef]
  24. Lundvall, O.; Klarbring, A. PREDICTION OF TRANSMISSION ERROR IN SPUR GEARS AS A CONSEQUENCE OF WEAR *. Mech. Struct. Mach. 2001, 29, 431–449. [Google Scholar] [CrossRef]
  25. Lundvall, O.; Klarbring, A. SIMULATION OF WEAR BY USE OF A NONSMOOTH NEWTON METHOD—A SPUR GEAR APPLICATION1-2 * †. Mech. Struct. Mach. 2001, 29, 223–238. [Google Scholar] [CrossRef]
  26. 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]
  27. Ding, H.; Kahraman, A. Interactions between Nonlinear Spur Gear Dynamics and Surface Wear. J. Sound Vib. 2007, 307, 662–679. [Google Scholar] [CrossRef]
  28. Zhou, C.; Dong, X.; Wang, H.; Liu, Z. Time-Varying Mesh Stiffness Model of a Modified Gear–Rack Drive with Tooth Friction and Wear. J. Brazilian Soc. Mech. Sci. Eng. 2022, 44, 213. [Google Scholar] [CrossRef]
  29. Huangfu, Y.; Zhao, Z.; Ma, H.; Han, H.; Chen, K. Effects of Tooth Modifications on the Dynamic Characteristics of Thin-Rimmed Gears under Surface Wear. Mech. Mach. Theory 2020, 150, 103870. [Google Scholar] [CrossRef]
  30. Liu, X.; Yang, Y.; Zhang, J. Investigation on Coupling Effects between Surface Wear and Dynamics in a Spur Gear System. Tribol. Int. 2016, 101, 383–394. [Google Scholar] [CrossRef]
  31. He, Z.; Hu, Y.; Zheng, X.; Yu, Y. A Calculation Method for Tooth Wear Depth Based on the Finite Element Method That Considers the Dynamic Mesh Force. Machines 2022, 10, 69. [Google Scholar] [CrossRef]
  32. Osman, T.; Velex, P. Static and Dynamic Simulations of Mild Abrasive Wear in Wide-Faced Solid Spur and Helical Gears. Mech. Mach. Theory 2010, 45, 911–924. [Google Scholar] [CrossRef]
  33. Trobentar, B.; Glodež, S.; Flašker, J.; Zafošnik, B. The Influence of Surface Coatings on the Tooth Tip Deflection of Polymer Gears. Mater. Tehnol. 2016, 50, 517–522. [Google Scholar] [CrossRef]
  34. VDI 2736, Thermoplastic Gear Wheels, German Standard; 2019.
  35. Lin, A.-D.; Kuang, J.-H. Dynamic Interaction between Contact Loads and Tooth Wear of Engaged Plastic Gear Pairs. Int. J. Mech. Sci. 2008, 50, 205–213. [Google Scholar] [CrossRef]
  36. Borovinšek, M. Exploit the Power of the CalculiX FEM Solver Using a Modern Graphical Interface 2023.
  37. CALCULIX, V. 2. 2. A Free Software Three-Dimensional Structural Finite Element Program 2023.
  38. The Wear of Metals under Unlubricated Conditions. Proc. R. Soc. London. Ser. A. Math. Phys. Sci. 1956, 236, 397–410. [CrossRef]
  39. Shen, X.; Cao, L.; Li, R. Numerical Simulation of Sliding Wear Based on Archard Model. In Proceedings of the 2010 International Conference on Mechanic Automation and Control Engineering; IEEE, June 2010; pp. 325–329. [Google Scholar]
  40. Perron, N.; Dellinger, N.; Prévereaud, Y.; Balat-Pichelin, M. 3D Mesh Displacement Strategy to Simulate the Thermal Degradation of Materials under Atmospheric Reentry Conditions. Acta Astronaut. 2022, 199, 293–312. [Google Scholar] [CrossRef]
  41. Hansen, G.A.; Douglass, R.W.; Zardecki, A. Mesh Enhancement; PUBLISHED BY IMPERIAL COLLEGE PRESS AND DISTRIBUTED BY WORLD SCIENTIFIC PUBLISHING CO., 2005; ISBN 978-1-86094-487-1.
  42. Borovinšek, M. Downloads - PrePoMax 2023.
  43. Bončina, T.; Polanec, B.; Zupanič, F.; Glodež, S. Wear Behaviour of Multilayer Al-PVD-Coated Polymer Gears. Polymers (Basel). 2022, 14, 4751. [Google Scholar] [CrossRef] [PubMed]
  44. Vincenzo Vullo Gears: Geometric and Kinematic Design; Springer Verlag, 2020.
Preprints 97110 g001
Figure 1. Outward normal vectors of the finite elements and the common node.
Figure 1. Outward normal vectors of the finite elements and the common node.
Preprints 97110 g001
Preprints 97110 g002
Figure 2. Computational model of gear pair.
Figure 2. Computational model of gear pair.
Preprints 97110 g002
Preprints 97110 g003
Figure 3. Flowchart of the computational analysis.
Figure 3. Flowchart of the computational analysis.
Preprints 97110 g003
Preprints 97110 g004
Figure 4. Finite element mesh of gear pair.
Figure 4. Finite element mesh of gear pair.
Preprints 97110 g004
Preprints 97110 g005
Figure 5. Sliding conditions in an arbitrary mating point Y: a) before pitch point, b) after pitch point.
Figure 5. Sliding conditions in an arbitrary mating point Y: a) before pitch point, b) after pitch point.
Preprints 97110 g005
Preprints 97110 g006
Figure 6. Wear of POM-gear after 10,000 cycles (Model 1).
Figure 6. Wear of POM-gear after 10,000 cycles (Model 1).
Preprints 97110 g006
Preprints 97110 g007
Figure 7. Wear of POM gear after 10,000 cycles (Model 2).
Figure 7. Wear of POM gear after 10,000 cycles (Model 2).
Preprints 97110 g007
Preprints 97110 g008
Figure 8. Worn geometry of POM-gear after 10,000 cycles (Model 2).
Figure 8. Worn geometry of POM-gear after 10,000 cycles (Model 2).
Preprints 97110 g008
Preprints 97110 g009
Figure 9. Wear of POM-gear after 10,000 cycles (Comparison between computational results and standardised approach according to VDI 2736).
Figure 9. Wear of POM-gear after 10,000 cycles (Comparison between computational results and standardised approach according to VDI 2736).
Preprints 97110 g009
Table 1. Basic parameters of the analysed gear pair [43].
Table 1. Basic parameters of the analysed gear pair [43].
Parameter Tested gear Supported gear
Material POM Steel (16MnCr5)
Normal module m 2.5 mm 2.5 mm
Pressure angle αn 20°
Helix angle β 0°
Number of teeth z 36 36
Tooth width b 14 mm 14 mm
Profile shift coefficient x 0
Centre distance a 90 mm
Basic rack profile ISO 53
Young’s modulus E 2,600 MPa 210,000 MPa
Poisson’s ratio ν 0.386 0.280
Lubrication Dry (no lubricated)
Wear coefficient kw (VDI 2736) 3.4⋅10−6 mm3/Nm
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated