Issue 
Mechanics & Industry
Volume 21, Number 1, 2020



Article Number  109  
Number of page(s)  15  
DOI  https://doi.org/10.1051/meca/2019076  
Published online  07 January 2020 
Regular Article
Computeraided design and loaded tooth contact analyses of bevel gear pair having straight teeth by different loaded torques
Department of Mechanical Engineering, University of Debrecen, Debrecen, Hungary
^{*} email: bodzassandor@eng.unideb.hu
Received:
8
May
2019
Accepted:
24
October
2019
The aim of this publication is to show how to integrate the designing process of straight bevel gears into a computer software so it can be further used for making the threedimensional (CAD) model of the gear pair. During meshing the driven and the driving gears engage gradually so contact points can be mathematically determined between the element pairs according to the developed coordinate system's arrangement. With mydeveloped computer software, I designed a certain bevel gear pair having straight teeth. Naturally, many different types of this pair could be designed. After that Loaded Tooth Contact Analysis (LTCA) have been done − by normal stress, normal deformation and normal elastic strain parameters. The aim of the LTCA is the analysis of the connection tooth zone in mechanical aspects by different loads. If the received parameters are not appropriate, you can return to the mechanical designing process where the starting parameters of the gear could be modified. Different load torques were used to determine the established mechanical parameters of the elements.
Key words: Bevel gear / module / LTCA / CAD / normal
© AFM, EDP Sciences 2019
1 Introduction
Straight bevel gears are applied widely in machinery (in vehicles, tool machines, for medical tools and machines etc.). They are used to connect shafts whose axes intersect in some angles, thus the meshing surfaces form a cone on which teeth are shaped (Fig. 1) [1–6].
One of the oldest process for producing gears applied even today is the Reinecker − Bilgram process. In this case, the vertical shaft of the machine intersects the centre of the pitch surface. The cutting tool makes an alternating movement while the operating pitch cone of the bevel gear unwinds the axoid of the deriving pitch cone (Fig. 2) [1,3–8].
During another production process of Heidenreich–Harbeck, the machine uses 2 cutting tools. This method is more advantageous because in this case the machine is moved in the axes perpendicular to the axes of the pitch surface. The machine is not set according to the λ dedendum angle, and the grip of the tool is stricter (Fig. 3) [1,3–8].
The paths of the cutting tools are not parallel during production, but they intersect in 2⋅θ angle, this is called planing angle, which can be calculated in the following way [5] (Fig. 4):(1)
Fig. 1 Angles and diameters of bevel gear. 
Fig. 4 Definition of planing angle. 
2 Defining the contact points
The profile curve of the bevel gear wheel is made as a result of the movement of the machining tool and the unwinding movement for the K_{1R} (x_{1R}, y_{1R}, z_{1R}) rotating coordinate system related to the driver gear wheel.
I am looking for the surface connected to K_{2R} coordinate system related to . I use the fact that the 2 surfaces during their movements coincide. Taking into consideration the following correlation [4,9]:(2)
I can state that the movement can be described with (ϕ_{1}) movement parameter [4,9,10].
The relative velocity between the 2 surfaces can be determined by the transformation between the rotating K_{1R} coordinate system of the driver wheel and the rotating K_{2R} coordinate system of the driven wheel [4,9]:(3)Transformation matrices between the K_{1R} and K_{2R} coordinate systems (Fig. 5):(4) (5)
Taking into consideration the correlation between the velocity vectors of the relative movement in K_{1R} and K_{2R} coordinate systems [4,9]:(6)in K_{1R} coordinate system, the relative velocity vector, based on (3) is:(7)where P_{1k} is the kinematic mapping matrix [4,9]:(8)
On the tooth surfaces of the meshing teeth, contact points mutually covering each other can be determined by solving the connection equation − which expresses the 1st Law of Contact − and the vectorscalar function simultaneously [4,9]:(9)
Fig. 5 Connection between the coordinate systems of the bevel gear drive pair. 
3 Designing the computer program
Input data for designing the drive pair are: m module, z_{1} number of teeth of the driver bevel gear wheel, z_{2} number of teeth of the driven gear wheel, c ^{*} root clearance factor and the α_{0} angle of contact [1,3–8].
For designing, the necessary geometrical correlations from technical literature was also used, which can be seen below based on Figures 1 and 5 [1,3–8]:

the largest pitch circle diameters:

half pitch angle of the driver gear wheel:

half pitch angle of the driven gear wheel:

effective pitch surface radius:

addendum on the largest diameter:

dedendum on the largest diameter:

the largest tip circle diameters:

the largest root circle diameters:

face width:

dedendum angle:

tip cone angles

root cone angles:
The output of my computer program is some calculated geometrical data of the bevel gear, and the profile curves (Fig. 6) in case of the smallest and largest diameters. Knowing the above mentioned manufacturing technologies [1,3–8] and the path of the cutting tool this program can determine the shape of the gear profiles.
The program saves the profiles of the bevel gear in .txt file format. Onto the given profile set of dots an interpolating Bspline can be fitted [11]. With the geometrical parameters of the drive pair and its profiles, the CAD model of the bevel gear was made by SolidWorks designing software (Fig. 7).
Fig. 6 Profile curves of the designed bevel gear pair (m = 15 mm, z_{1 }= 20, z_{2 }= 30, Σ = 90°). (a) Profile curves of the driver gear. (b) Profile curves of the driven gear. 
Fig. 7 CAD model of the bevel gear pair having straight teeth (m = 15 mm, z_{1} = 20, z_{2} = 30, Σ = 90°). 
4 LTCA analysis
The aims of the so called “Loaded Tooth Contact Analysis” (LTCA) are to model the connection of toothed drive pairs and to simulate the model by given loads and boundary conditions. During the LTCA in case of this given geometrical bevel gear (Figs. 6 and 7, Tab. 1) normal stress, normal elastic strain and normal deformation values were analysed on the surface of the driven and driver gear wheels [12–21].
Parameters of the designed bevel gear pair.
4.1 Material quality, FE mesh (finite element)
The property of the material used for our bevel gear is given in Table 2.
During defining the FE mesh in the contact zone, sphere volume (37 mm radius) including dense triangles (1 mm mesh dimension) has been applied on the toothed area (Fig. 8). Automatic meshing was used on the outside areas. The friction coefficient in the tooth contact zone is µ = 0.01, because I considered the liquid and rolling frictions on the contact zone.
Material properties.
Fig. 8 FE mesh. 
4.2 Loads and boundary conditions
Four coordinate systems have been defined for LTCAs: K_{s} − absolute static, K_{s1}–static related to the driver wheel, K_{s2}–static, related to the driven wheel and K_{c}– coordinate system in the tooth contact zone.
During the analyses, the gear having less number of teeth (driver gear) drives the gear having more number of teeth. Five degrees of freedom were fixed of the driver gear, only rotation along the rotational axis was allowed. The driver gear was loaded by M = 200–400 Nm torque by 50 Nm steps. Six degrees of freedom of the driven gear were fixed (Fig. 9).
Fig. 9 Setting loads and boundary conditions. 
Fig. 10 Defining normal stress. 
4.3 Normal stress analyses
In a fixed P point of the body, signals the stress vector awakening on the normal surface elements. It can be divided into a normal directional component and one perpendicular to it [22–24] (Fig. 10):(22)where the normal stress is [22–24]:(23)
Normal stress partitions and average normal stress values have been detected as an effect of the different loads, on the tooth surfaces of the driver and driven gear (Fig. 11).
Results are shown in the diagram of Figure 12. In absolute value, as an effect of the increasing load moment, normal stress values also increase on the tooth surfaces of the driver and driven gear.
The shapes of the two diagrams are the same and they are parallel. The normal stresses of the driven gear are a little higher than the normal stresses of the driver gear in absolute value. Because of the good comparison the same limits were set.
Fig. 11 Normal stress values. (a) M = 200 Nm: driver gear, MPa (top), driven gear, MPa (bottom). (b) M = 250 Nm: driver gear, MPa, driven gear, MPa. (c) M = 300 Nm: driver gear, MPa, driven gear, MPa. (d) M = 350 Nm: driver gear, MPa, driven gear, MPa. (e) M = 400 Nm: driver gear, MPa, driven gear, MPa. 
Fig. 11 Continued. 
Fig. 12 Load − normal stress diagram. 
4.4 Normal elastic strain analyses
The state of a body's deformation can be signalled by deformation tensor, which can be determined from the displacement vector as seen below [22–24]:(24) (25)where(26)
Based on the values of ε_{x}, ε_{y}, ε_{z} elastic strain per unit length and unit vectors related to them, the FE program calculates the normal direction resultant (ε_{n} ), which is perpendicular to the surface. As an effect of the loads on the surface of the driver and driven gear, normal elastic strain partitions and average normal elastic strain values occur (Fig. 13).
The results are shown in the diagram of Figure 14. You can see that, in absolute value, as an effect of increasing loads, normal elastic strain values also increase on the tooth surfaces of the driver and driven gears.
The shapes of the two diagrams are the same and they are parallel. The normal elastic strains of the driven gear are a little higher than the normal elastic strains of the driver gear in absolute value. Because of the good comparison the same limits were set.
Fig. 13 Normal elastic strain values as a result of increasing loads. (a) M = 200 Nm: driver gear, −0.0000114 (top), driven gear, −0.0000127 (bottom). (b) M = 250 Nm: driver gear, −0.0000141, driven gear, −0.0000154. (c) M = 300 Nm: driver gear, −0.0000168, driven gear, −0.0000183. (d) M = 350 Nm: driver gear, −0.0000197, driven gear, −0.0000214. (e) M = 400 Nm: driver gear, −0.0000227, driven gear, −0.0000246. 
Fig. 13 Continued. 
Fig. 14 Load − normal elastic strain diagram. 
4.5 Normal deformation analyses
As an effect of a load, the points of a solid body move, and the beginning state of the body differs from the state after full load application [22–24].
For instance, the position of an arbitrary P point of the body is given by the following placement vector, in the state before the load was applied [22–24]:(27)
Signals P′ the state of P point after loading, the related placement vector is the displacement vector is u_{P}. Based on Figure 15 you can state the following correlations [22–24]:(28)
The displacement vector is the function of the position of the points of the body before loading, i.e. the placement vector [22–24]:(29)
The sum of the displacement vectors related to the points of the body is called displacement field, where the displacement vector is [22–24]:(30)
The coordinates of the displacement:(31)
As an effect of the different loads, on the tooth surfaces of the driver and driven gear, the following normal deformation values (in this case x direction deformation) and average normal deformation values (the average of the x direction values) occur (Fig. 16).
Calculated results are shown in Figure 17. You can see that if the loads are increased the values of normal deformation also increase on the tooth surfaces of the driver and driven gear.
The shapes of the two diagrams are almost parallel. The normal deformations of the driver gear are higher than the normal deformations of the driven gear in absolute value. The teeth are in contact that is why the presages of this specific are different. Because of the good comparison the same limits were set.
Fig. 15 Definition of displacement vector. 
Fig. 16 Normal deformation values. (a) M = 200 Nm: driver gear, −0.0654 mm (top), driven gear, 0.00166 mm (bottom). (b) M = 250 Nm: driver gear, −0.0644 mm, driven gear, 0.00207 mm. (c) M = 300 Nm: driver gear, −0.0636 mm, driven gear, 0.00247 mm. (d) M = 350 Nm: driver gear, −0.0629 mm, driven gear, 0.00286 mm. (e) M = 400 Nm: driver gear, −0.0621 mm, driven gear, 0.00325 mm. 
Fig. 16 Continued. 
Fig. 17 Load–normal deformation diagram. 
5 Conclusion
Straight bevel gears are widely used in machinery where the axes of the two shafts intersect, and the faces of the gears themselves form a cone. They are used as load transmission equipment, to change the number of rotation etc. Their further development (from constructional and technological point of view) and their LTCA based on different mechanical properties are highly needed.
Based on the geometrical and production analyses of the straight bevel gear pair, I developed a computer program to be able to design the drive pair. This program calculates all the geometrical parameters of the drive pair and draws the profile curves (which occur during production), on the smallest and largest diameters.
The given profile curves can be analyzed from constructional point of view and they are important for making our CAD model.
After that a CAM (ComputerAided Manufacturing) software can be used to simulate the total production process of the drive pair.
It is important to have the CAD model for LTCA as well. These analyses are good for analyzing the mechanical properties on the driver and driven gear (such as stress, strains, elastic strain). LTCA are also important for having the appropriate constructional planning, to be able to choose the appropriate raw material, and to get important information about loads.
Normal stress, normal deformation and normal elastic strain values have been analyzed on the tooth surfaces of the driver and driven gear pair by different load torques. The results have been averaged, analyzed and shown in diagrams.
Nomenclature
: Relative velocity vector (mm/min^{−1})
: Normal vector of the surface in K_{1R} coordinate system
: Placement vector of the moving points of the profile curve
: Placement vector of moving point of conjugated surface
M_{2R,1R}, M_{1R,2R} : Translation matrices between K_{1R} and K_{2R} coordinate systems
P_{1k} : Matrix of kinematic mapping
c ^{*} : Root clearance factor
d_{0} : Largest pitch circle diameter (mm)
d_{f} : Largest tip circle diameter (mm)
d_{l} : Largest root circle diameter (mm)
j_{s} : Clearance at flank (mm)
K_{1R} (x _{1R}, y _{1R}, z _{1R}): Rotational coordinate system related to the driver gear wheel
K_{1S} (x_{1S}, y _{1S}, z _{1S}): Static coordinate system related to the driver gear wheel
K_{2R} (x _{2R}, y _{2R}, z _{2R}): Rotational coordinate system related to the driven gear wheel
K_{2S} (x _{2S}, y _{2S}, z _{2S}): Static coordinate system related to the driven gear wheel
R_{e} : Effective pitch surface radius (mm)
S_{ax} : Pitch circle face width on the largest diameter (mm)
t : Circular pitch on the largest pitch circle diameter (mm)
z_{1}, z_{2} : Number of teeth
ε_{n} : Normal elastic strain (mm)
ε_{x}, ε_{y}, ε_{z} : Elastic strain values per unit length (mm)
ϕ_{1R}, ϕ_{2R} : Angular displacement (°)
E : Distance between the vertex of cones (mm)
u, v, w : Coordinates of the displacement vector (mm)
: Parts of the differential operator
γ _{xy,} γ _{yx,} γ_{xz} : Shear specific elongation (mm)
Acknowledgments
This research was supported by the János Bolyai Research Scholarship of the Hungarian Academy of Sciences. The publication itself is partly supported by the EFOP3.6.116201600022 project, and it is cofinanced by the European Union and the European Social Fund.
References
 D.W. Dudley, Gear Handbook, MC Graw Hill Book Co., New York, 1962 [Google Scholar]
 R. Gołębski, A. Szarek, Diagnosis of the operational gear wheel wear, Technical Gazette 26, 658–661 (2019) [Google Scholar]
 J. Klingelnberg, Bevel Gear, Fundamentals and Applications, Springer, Berlin, 2016 [CrossRef] [Google Scholar]
 F.L. Litvin, A. Fuentes, Gear Geometry and Applied Theory, Cambridge University Press, Cambridge, 2004 [CrossRef] [Google Scholar]
 V. Rohonyi, Fogaskerékhajtások. Műszaki Könyvkiadó, Budapest, 1980 [Google Scholar]
 Z. Terplán, Gépelemek IV., Kézirat, Tankönyvkiadó, Budapest, 1975, p. 220 [Google Scholar]
 J. Argyris, A. Fuentes, F.L. Litvin, Computerized integrated approach for design and stress analysis of spiral bevel gears, Comput. Methods Appl. Mech. Eng. 191, 1057–1095 (2002) [Google Scholar]
 I. PopaMüller, Generation of eloid bevelgear on toothcutting machine with the generating facegear, 15th International Conference in Mechanical Engeneering, ClujNapoca, 2007 [Google Scholar]
 F.L. Litvin, Gear geometry and applied theory, Prentice Hall, Englewood Cliffs, NJ, 1994 [Google Scholar]
 A. Fuentes, J.L. Iserte, I. GonzalezPerez, F.T. SanchezMarin, Computerized design of advanced straight and skew bevel gears produced by precision forging, Comput. Methods Appl. Mech. Eng. 200, 2363–2377 (2011) [Google Scholar]
 I. Juhász, Görbékés felületek modellezése, Miskolci Egyetem, elektronikus jegyzet, p. 114 [Google Scholar]
 Z. Fangyan, Z. Mingde, Z. Weiqing, T. Rulong, G. Xiaodong, On the deformed tooth contact analysis for forged bevel gear modification, Mech. Mach. Theory 135, 192–207 (2019) [Google Scholar]
 G.D. Bibel, A. Kumar, S. Reddy, R. Handschuh, Contact stress analysis of spiral bevel gears using finite element analysis, J. Mech. Des 117, 235–240 (1995) [CrossRef] [Google Scholar]
 F.L. Litvin, A. Fuentes, I. GonzalesPerez, L. Carnevali, T.M. Sep, New version of NovikovWildhaber helical gears: computerized design, simulation of meshing and stress analysis, Comput. Methods Appl. Mech. Eng. 191, 5707–5740 (2002) [Google Scholar]
 F.L. Litvin, A. Fuentes, I. GonzalezPerez, L. Carvenali, K. Kawasaki, R.F. Handschuh, Modified involute helical gears: computerized design, simulation of meshing and stress analysis, Comput. Methods Appl. Mech. Eng. 192, 3619–3655 (2003) [Google Scholar]
 S. Peng, H. Ding, G. Zhang, J. Tang, Y. Tang, New determination to loaded transmission error of the spiral bevel gear considering multiple elastic deformation evaluations under different bearing supports, Mech. Mach. Theory 137, 37–52 (2019) [Google Scholar]
 I. PopaMüller, Casestudy and simulation by meshing of bevelgears, 14th International Conference in Mechanical Engeneering, TgMures, 2006 [Google Scholar]
 I. PopaMüller, Simulation of Generation by Meshing of Conical Gear Octoid II with Hobbing Machine HeidenreichHarbeck, 18th International Conference in Mechanical Engeneering, Baia Mare, 2010 [Google Scholar]
 J.K. Theodore, Tooth contact analysis of spiral bevel and hypoid gears under load, SAE Trans. 90, 2205–2216 (1981) [Google Scholar]
 H. Zehua, D. Han, P. Shandong, T. Yi, T. Jinyuan, Numerical determination to loaded tooth contact performances in consideration of misalignment for the spiral bevel gears, Mech. Mach. Theory 151, 343–355 (2019) [Google Scholar]
 W. Wagner, S. Schumann, B. Schlecht, Cosimulation of the Tooth Contact of Bevel Gears within a Multibody Simulation, Forsch Ingenieurwes 83, 425–433 (2019) [CrossRef] [Google Scholar]
 I. Kozák, Gy. Szeidl, Fejezetek a szilárdságtanból, Kézirat, 2008–2012, Miskolci Egyetem, p. 284. elektronikus jegyzet [Google Scholar]
 S. Moaveni, Finite Element Analysis, Theory and Application with ANSYS, Pearson Education Limited, London, 2015, p. 928 [Google Scholar]
 I. Páczelt, T. Szabó, A. Baksa, A végeselem módszer alapjai, Miskolci Egyetem, p. 243 [Google Scholar]
Cite this article as: S. Bodzás, Computeraided design and loaded tooth contact analyses of bevel gear pair having straight teeth by different loaded torques, Mechanics & Industry 21, 109 (2020)
All Tables
All Figures
Fig. 1 Angles and diameters of bevel gear. 

In the text 
Fig. 2 Reinecker–Bilgram principle [6]. 

In the text 
Fig. 3 Heidenreich–Harbeck principle [6]. 

In the text 
Fig. 4 Definition of planing angle. 

In the text 
Fig. 5 Connection between the coordinate systems of the bevel gear drive pair. 

In the text 
Fig. 6 Profile curves of the designed bevel gear pair (m = 15 mm, z_{1 }= 20, z_{2 }= 30, Σ = 90°). (a) Profile curves of the driver gear. (b) Profile curves of the driven gear. 

In the text 
Fig. 7 CAD model of the bevel gear pair having straight teeth (m = 15 mm, z_{1} = 20, z_{2} = 30, Σ = 90°). 

In the text 
Fig. 8 FE mesh. 

In the text 
Fig. 9 Setting loads and boundary conditions. 

In the text 
Fig. 10 Defining normal stress. 

In the text 
Fig. 11 Normal stress values. (a) M = 200 Nm: driver gear, MPa (top), driven gear, MPa (bottom). (b) M = 250 Nm: driver gear, MPa, driven gear, MPa. (c) M = 300 Nm: driver gear, MPa, driven gear, MPa. (d) M = 350 Nm: driver gear, MPa, driven gear, MPa. (e) M = 400 Nm: driver gear, MPa, driven gear, MPa. 

In the text 
Fig. 11 Continued. 

In the text 
Fig. 12 Load − normal stress diagram. 

In the text 
Fig. 13 Normal elastic strain values as a result of increasing loads. (a) M = 200 Nm: driver gear, −0.0000114 (top), driven gear, −0.0000127 (bottom). (b) M = 250 Nm: driver gear, −0.0000141, driven gear, −0.0000154. (c) M = 300 Nm: driver gear, −0.0000168, driven gear, −0.0000183. (d) M = 350 Nm: driver gear, −0.0000197, driven gear, −0.0000214. (e) M = 400 Nm: driver gear, −0.0000227, driven gear, −0.0000246. 

In the text 
Fig. 13 Continued. 

In the text 
Fig. 14 Load − normal elastic strain diagram. 

In the text 
Fig. 15 Definition of displacement vector. 

In the text 
Fig. 16 Normal deformation values. (a) M = 200 Nm: driver gear, −0.0654 mm (top), driven gear, 0.00166 mm (bottom). (b) M = 250 Nm: driver gear, −0.0644 mm, driven gear, 0.00207 mm. (c) M = 300 Nm: driver gear, −0.0636 mm, driven gear, 0.00247 mm. (d) M = 350 Nm: driver gear, −0.0629 mm, driven gear, 0.00286 mm. (e) M = 400 Nm: driver gear, −0.0621 mm, driven gear, 0.00325 mm. 

In the text 
Fig. 16 Continued. 

In the text 
Fig. 17 Load–normal deformation diagram. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.