^{1}

^{1}

^{1}

The paper describes the dynamic characteristic and fatigue life of the gear meshing. Use the ANSYS to deal with the problem of motion errors which are the result of the wear and the curve caused by the rough surface of the gear generated by the low machining accuracy in the process. And through the analysis, the value and the changing trend of the equivalent stress, fatigue life, and fatigue sensitivity with different friction coefficients can be known. Then it provides an important theoretical basis to strengthen the strength and stability, reduce the deformation of the gear, and improve fatigue life.

Gear transmission is the most widely applied transmissions among the mechanical transmissions. It has a strict transmission ratio and transfers smoothly. Because the transmission ratio can be kept constant, the gear transmission is widely applied in many kinds of gearboxes, heavy equipment, and other areas [

The common methods of gear processing can be divided into copying and generating by the theory. The gear processed by copying has low accuracy and the accuracy is always below 11 degree. Because the generating bases on the theory of gear meshing, the accuracy is higher. During the transmission process of gears with varieties accuracy, the different friction coefficient between the faces of gears has a direct effect on the gears’ life. So, during the transmission, the face of gear often becomes invalid. Lightly there are gear wearing and gear pitting, and worse there are gear bonding even tooth fracture. These failures all affect the transmission ratio and the stability of the gear transmission system so that the life of the mechanical system decreases.

Scholars at home and abroad have conducted in-depth research on this issue and have achieved fruitful results. Chaari et al. used the Weber method material mechanics to analyze the effect of tooth surface flaking on the meshing stiffness by the gear deformation and hertz contact deformation [

To solve these problems, in this paper, a pair of ordinary straight tooth cylindrical involute gears have been taken as research object. The 3-D geometry model was established by UG, and the ANSYS was used to analyze the transient contact dynamic characteristics. The deformation, stress change, and fatigue change have been obtained from the analysis. During the process of the gear meshing transmission, it provides a common conclusion and solution for the great force of the gear tooth, serves deformation and fatigue failure.

Transient dynamics analysis is used to determine the dynamic response of the time-varying load structure. The transient dynamics analysis can determine the time-varying displacement, strain, stress, and force response of the structure under the random combination of the steady load, transient load, and simple harmonic load.

The motion equation of transient dynamics analysis:

[ M ] { u ¨ } + [ C ] { u ˙ } + [ K ] { u } = { F ( t ) } (1)

The load { F ( t ) } is the function of time. [ M ] is the mass matrix. [ C ] is the damping matrix. [ K ] is the stiffness matrix. { u } is the displacement matrix.

Based on the linear static structure analysis, the driving wheel and the driven wheel are set as elastic bodies A and B respectively. When the A and B contacts, there will be innumerable contact points and they can make respective displacement { U A } and { U B } under the external load { P A } and { P B } . According to the basic theory of Finite Element, the balanced equation can be obtained:

[ K A ] { U A } = { P A } + { R A } (2)

[ K B ] { U B } = { P B } + { R B } (3)

[ K A ] and [ K B ] are the stiffness matrixes of A and B; { U A } and { U B } are the nodes displacement vectors; { P A } and { P B } are the external load of A and B; { R A } and { R B } are the contact force vectors of A and B.

According to the analysis of the balanced equations, once the characteristics of materials and external load have been determined, the [ K A ] , [ K B ] , { P A } and { P B } also have been determined. But there are still four unknown parameters that cannot be solved by two equations. So the elastic body contact nature should be added, and different natures correspond to different displacement and supplementary equations of contact force relationship. Last the balanced equations are solved.

For the process of solving equations with ANSYS, the stiffness matrixes [ K ] and external load { P } are taken into the balanced equations, and according to the hypothesis of contact nature, the contact displacement { U } and contact force { R } are solved. Then the obtained { U } and { R } are compared with the hypothesis contact nature. If they are not the same, other hypotheses will be needed and the results will be calculated again until the results are the same with the hypothesis [

In this paper, all the gears are a pair of standard mounting external meshing and standard spur gears. The main parameters are Z 1 = 22 , Z 2 = 33 , α = 20 ∘ , h a * = a , m = 2.5 mm , C * = 0.25 , b = 20 mm . Among these gears, Z 1 was the driving wheel and Z 2 was driven wheel.

The CG toolbox in UG was used to set gear meshing. In order to reduce the difficulty of mesh dividing and the requirements for computer hardware and increase the speed of calculation, under the condition of guaranteeing the accuracy, the keyways and chamfers on gears were omitted. The three-D model is shown in

The established gear engagement model was imported into ANSYS. The material was defined as structure steel. Considering the real working condition of gears, the type of contact of the contact gear surface was set as Frictional, and the friction coefficient was initialized as 0.01. Among this, the Contac was set as a driving wheel and Target was set as a driven wheel [

the 3D tetrahedral mesh intelligent division was used. So there were 23,519 notes and 4065 elements, shown as

According to the calculation of the finite element model, the calculation results were imported into the postprocessor and the equivalent stress distribution cloud diagram during gear meshing with the friction coefficient 0.01 was obtained, shown as

Changing the friction coefficient in the contact set to simulate the gear meshing in different processing accuracy, the equivalent stress distribution cloud diagram during gear meshing with the friction coefficient 0.05, 0.1, 0.15, and 0.2 could be obtained, shown in

From

With the increasing of the friction coefficient, it means that the processing accuracy of the gear surface gets decreased, the surface gets roughness and the

wear of gear gets increased. So the equivalent stress in the area of gear meshing increases obviously and the minimum equivalent stress in the external area of the gear engagement decreases.

The reason is the increase of the friction coefficient, and it causes the friction drag increased during the relative motion of the gear faces. It makes the gear faces in the gear meshing area sustains more load.

The calculation results were imported into the Fatigue Tool. In addition, the type of fatigue load was set as Zero-Based and the fatigue stress factor was also set as 0.8, shown as

After post-processing, the fatigue life and sensitivity with friction coefficient 0.01, 0.05, 0.1, 0.15, and 0.2 were obtained, shown as

When the friction coefficient keeps static, from the fatigue life cloud diagram, the tooth root got fatigue failure first. With the increase of friction coefficient, the minimum fatigue life was getting shorter. So it meant that the processing accuracy gets lower, the gear gets fatigue invalid more easily.

In the curve of fatigue sensitivity, the abscissa represents the number of cycles. The whole graph represents, under a certain number of cycles, the maximum fatigue stress which can be sustained by the components with the alternating loads. When the friction coefficient gets higher and the cycle number of alternating loads keeps static, the maximum fatigue stress sustained by the gear decreases.

During the process of the gear meshing, the relative sliding speed at the nodal line during the driving and driven wheels are opposite, which means that the sliding speed at the nodal line is zero. At this time, the friction coefficient at the nodal line is the static friction coefficient. But the static friction coefficient is higher than the dynamic friction coefficient. According to the surface friction equation, when the friction coefficient is higher, and it causes the friction is higher. So the possibility of fatigue failure is higher [

F = f ⋅ N (4)

Within the increasing of the friction coefficient, the change limitation of the load is getting lower. It means that, under the condition of the critical fatigue life, a higher friction coefficient causes lower processing accuracy. Then it causes the load range sustained by the components to get smaller, so, under the condition of a long time running, fatigue failure happens more easily.

According to the transient contact dynamic analysis of the gear meshing, the stress distribution, change the trend, and the range of fatigue failure has been obtained. The conclusions are shown as follow:

According to the transient contact dynamic analysis, the position of the maximum stress during the gear meshing has been obtained; changing the friction coefficient of the contact area to simulate the different processing accuracy, the trend of equivalent stress change of the gear meshing has also been obtained.

After the analysis of the gear meshing fatigue life, the area, under the different friction coefficient during the gear meshing, where the fatigue failure happens more easily has been obtained. And the limit value of the load change under the critical fatigue condition has been obtained.

This can provide a theoretical basis for the analysis of the strength analysis, fatigue life prediction, and the adjustment of the processing plans.

The authors express their sincere thanks to the workmates giving a lot of help to this study.

The authors declare no conflicts of interest regarding the publication of this paper.

Zou, X., Yang, C.W. and Liu, B.L. (2020) The Contact and Transient Dynamic Analysis of Gear Meshing with ANSYS. Open Access Library Journal, 7: e6771. https://doi.org/10.4236/oalib.1106771