Investigation of the Effects of Profile Shift in Helical Gear Mechanisms with Analytical and Numerical Methods

In this paper, the effects of profile shift in cylindrical helical gear mechanisms have been investigated with numerical and analytical calculations. The mathematical model for computer simulation of gears has been designed and the numerical calculations have been carried out. Analytical calculations have been made with an excel program which was designed at different profile shift coefficients for a selected mechanism. Analytical calculations of the same mechanism have been verified by using ANSYS 14.5. The results of analytical and numerical solutions have been compared to profile shift coefficients.


Introduction
Cylindrical gears are essential machine elements and widely used in many fields, such as wind power, coastal ships, automotive products, and general gear boxes.
The helical gears offer considerable advantages for high speed and heavy-duty gear drives.For helical gears the total contact ratio is higher owing to the axial tooth overlap.Helical gears, therefore, tend to have greater load-carrying capacity compared with spur gears which material, quality of manufacture, lubrication, type of loading, and other factors are the same.They also possess the ability to a much smoother and quieter running [1] [2]. the separating and tangential (driving) loads of spur gears, where suitable means can be provided to take this thrust, such as thrust collars or ball or tapered-roller bearings; it is no great disadvantage.The axial load increases as the helix angle increases.Helix angles typically range from 15˚ to 45˚.At high values of helix angle (for 45˚), the improvement in noise tends to peak [3].

Helical gears produce an end thrust along the axis of the shafts in addition to
There are two basic strength calculations for two primary fatigue related failure modes in the design of gear.The first is to calculate the maximum bending stress which leads to fracture at the root of the gear tooth.The most common type of surface failure is pitting which is caused by the repeating high contact stress (Hertz surface pressure).The second job is to calculate the maximum surface pressure [3].
Profile shift is used by designers for many reasons including the elimination of undercut when hobbing or shaping gears with few teeth, increasing the load carrying capacity or the lifetime of the gears, decreasing the noise level, increasing or decreasing the center distance and changing tooth thickness.Profile shift-sometimes known as "addendum modification"-is the displacement of the basic rack (or cutting tool) datum line from the reference diameter of the gear.The size of the profile shift is usually made to be non-dimensional by dividing it by the normal module, and it is then defined by the profile shift coefficient "x".A positive profile shift increases the tooth thickness while a negative profile shift reduces tooth thickness.Some of the studies done recently on profile shift in cylindrical gear mechanisms have been mentioned below.
Antal et al. presented a new method for the determination of the specific addendum modifications at helical gears.The method is based on the equalization of the relative velocities at the points where the meshing of the teeth begins and ends [4].
Senthil Kumar et al. worked on optimization of asymmetric spur gear drives to improve the bending load capacity.They concluded that the optimum value of profile shift coefficient for pinion increased with the increase in the values of speed ratio and teeth number in pinion [5].Mallesh et al. investigated the effect of profile shift in asymmetric spur gears on tooth root bending stress by finite element analysis [6].

Magalhaes et al. investigated the influence of tooth profile and oil formulation
on gear power loss.In the study mentioned the power losses reduction has been obtained using two different approaches: using lower modulus helical gears and significant positive profiles shifts and using gear oil formulations with different base oils [7].
Baglioni et al. investigated the influence of the addendum modification on spur gear efficiency [8].
Li et al. worked on optimal selection of addendum modification coefficients of involute cylindrical gears [9].
In this paper, the effects of profile shift in cylindrical involute helical gear mechanisms have been investigated with numerical and analytical calculations.
The mathematical model for computer simulation of helical gears has been designed and the numerical calculations have been carried out.Analytical calculations were made with an excel program which was designed at different profile shift coefficients for a selected gear mechanism.Analytical calculations of the same mechanism were verified by using ANSYS 14.5.In this paper, the effects of profile shift coefficient on parameters such as center distance, normal pressure angle, transverse contact ratio, pinion tooth root stress, gear tooth root stress, Hertzian contact stress (surface pressure at tooth flank) in helical gear mechanisms have been investigated with numerical and analytical methods.Analytical solutions have been carried out according to DIN 3990.Numerical solutions have been performed using the finite element method.The most realistic modeling of the operation of the all gear mechanism has been done with the aid of three-dimensional finite-element analysis.It has also not been applied any external load.The results of analytical and numerical solutions have been compared to profile shift coefficients [10].

Methodology
In analytical solutions, initially, module and other size calculations are done for helical gear mechanism given in Table 1.In gear systems, there has not been any change except profile shift coefficients.Input parameters of gear systems, the main geometric dimensions and tooth forces remained constant.Center distances, pressure angles, transverse contact ratios, pinion tooth root stresses, gear tooth root stresses, Hertzian contact stresses have been calculated according to DIN 3990.
According to the data in Table 1, the module and the number of teeth of gear were chosen as m = 5 and z g = 59.Pinion profile shift coefficient (x p ) is selected according to the following equation proposed by DIN 3992 where x g represents gear profile shift coefficient [11].
Thus, tooth root stresses for the pinion and the gear have been expected in the near values.The most frequently used numerical method is the finite element method, as can be seen from many papers [12].Numerical solutions were made with the program ANSYS 14.5.
Tooth geometries were made according to dimensional solutions obtained during analytical calculations.Accordingly, tooth geometry change of the pinion depending on profile shift coefficient is shown in Figure 1.
In finite element analysis, for the purpose of the closest modeling to the actual situation, the whole gear mechanism is modeled as three-dimensional to prevent application of an external force.In finite element models of geometries, ANSYS SOLID 186 elements with 20 nodes were used.
In all models examined, edges of the elements have generally a length of 2 mm.In tooth regions where are in contact with each other, edges of the elements have a length of 0.5 mm.The 0.5 mm of edge length has been defined based on a grid independence study and the details of this study can be found in Table 2.
The total number of nodes for each model is over two million.The total number of elements is over five hundred thousand.The finite element model created for non-shifted gear mechanism is shown in Figure 2.   At analytical solution, input torque is calculated as In the central part of the pinion (colored red), torque is applied clockwise.All degrees of freedom of blue surface in the middle of the gear are restricted in order to maintain balance and to read stresses in analysis (Figure 3).
In meshing surfaces of teeth, frictionless contacts are defined as target (colored blue) and contact (colored red) (Figure 4).Thus, the teeth surfaces can slide over one another and leave each other.In models, due to contacts, the nonlinear analysis was performed for solutions.

Results
Profile shift coefficient for each mechanism was implemented by increasing or decreasing x = 0.1 for models used in the analytical solution.Results of all the mechanisms calculated in these values are shown in Table 3 where σ p , σ g , and σ H represents the tooth root stress of pinion, tooth root stress of gear and Hertzian contact stress respectively.x t is the total profile shift coefficient.
In the analytical calculation, numerical solutions are made for four different gear mechanisms.These are negative shifted (x t = −0.5),non-shifted (x t = 0), and positive shifted (x t = +0.5 and x t = +1) gear mechanisms.Maximum tooth root stress and Hertzian contact stress for each solution is calculated.The results which are obtained by solving the general equation of motion in ANSYS software are given in Table 4.
The distribution of stress on the teeth for couple and pinion in case of non-shifted gear mechanism are shown in Figure 5. Hertzian pressure distribution on the tooth surface of pinion for different profile shift values are shown in Figure 6.The following conclusions have been reached in the analytical solution.According to non-shifted gear mechanism, the center distance has decreased approximately 0.27% for negative profile shift of each x = −0.1.According to non-shifted gear mechanism, the center distance has increased approximately 0.25% for positive profile shift of each x = +0.1.
According to non-shifted gear mechanism, the pressure angle has decreased approximately 2.24% for negative profile shift of each x = −0.1.According to non-shifted gear mechanism, the pressure angle has increased approximately 1.88% for positive profile shift of each x = +0.1.
Due to small number of teeth of the pinion, only positive profile shift is applied to prevent undercut in pinion.According to non-shifted gear mechanism, the pinion tooth root stress has decreased approximately 5.44% for positive profile shift of each x = +0.1.
According to non-shifted gear mechanism, the gear tooth root stress has increased approximately 3.63% for negative profile shift of each x = −0.1.According to non-shifted gear mechanism, the gear tooth root stress has decreased approximately 1.52% for positive profile shift of each x = +0.1.
According to non-shifted gear mechanism, Hertzian contact stress has increased approximately 1.14% for negative profile shift of each x = −0.1.According to non-shifted gear mechanism, Hertzian contact stress has decreased approximately 0.83% for positive profile shift of each x = +0.1.
Profile shift coefficient for each mechanism was implemented by increasing or decreasing x = 0.5 for models used in the numerical solution.
Tooth root stresses and Hertzian contact stresses are increased with negative profile shift, while they are decreased with positive profile shift.
Compared to the numerical simulation results, the pinion tooth root stress has decreased approximately 3.42% for each x = +0.1 positive profile shift coefficient, according to non-shifted gear mechanism.
If the analytical and the numerical results are compared, the analytical results are found to be at higher values.This is because the analytical calculation is done for the most critical case.As seen from Figure 7 and Figure 8, the analytical-numerical result differences from negative profile shifting towards positive profile shifting have been decreased.
While the difference for tooth root stress of non-shifted pinion is 20.72%, the difference is 6.58% for positive shifted pinion of x p = 0.5 (Figure 7).While the difference for tooth root stress of negative shifted gear of x g = −0.69 is 18.85%, the difference is 2% for positive shifted gear of x g = 0.5 (Figure 8).As seen from Figure 9, Hertzian contact stresses in the analytical solution have been decreased approximately 28% compared to numerical results.

Conclusion
Consequently, in this study, the results of the analytical solution, have been found to be larger according to the results of the numerical solution.The differences of tooth root stresses found through analytical and numerical methods is How to cite this paper: Karadere, G. and Yilmaz, I. (2018) Investigation of the Effects of Profile Shift in Helical Gear Mechanisms with Analytical and Numerical Methods.World Journal of Mechanics, 8, 200-209.https://doi.org/10.4236/wjm.2018.85015G. Karadere, I. Yilmaz DOI: 10.4236/wjm.2018.85015201 World Journal of Mechanics

Figure 1 .
Figure 1.Influence of pinion profile shift coefficient on the pinion tooth geometry (for numerical solution).

Figure 2 .
Figure 2. Finite element model for non-shifted gear mechanism (discretization of the volume by finite elements).

Figure 5 .
Figure 5.The distribution of von-Mises stress on the teeth for couple and pinion.

Figure 6 .
Figure 6.Hertzian pressure distribution on the tooth surface of pinion for different profile shift values.

Figure 7 .
Figure 7. Influence of pinion profile shift coefficient on pinion tooth root stress (Comparison of analytical and numerical solutions).

Figure 8 .
Figure 8. Influence of gear profile shift coefficient on gear tooth root stress (Comparison of analytical and numerical solutions).

Figure 9 .
Figure 9. Influence of total profile shift coefficient on Hertzian contact stress (Comparison of analytical and numerical solutions).

Table 1 .
Properties of selected in volute helical gear mechanism.

Table 2 .
Grid independece study of the non-shifted case.

Table 3 .
Analytical solution results for profile shifted gear mechanisms.

Table 4 .
Numerical solution results for profile shifted gear mechanisms.