EDP Sciences logo
Submit your paper
Search
Menu
All issues Volume 21 / No 1 (2020) Mechanics & Industry, 21 1 (2020) 109 Full HTML
Open Access
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

© 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)

thumbnail Fig. 1

Angles and diameters of bevel gear.
thumbnail Fig. 2

Reinecker–Bilgram principle [6].
thumbnail Fig. 3

Heidenreich–Harbeck principle [6].
thumbnail 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 K1R (x1R, y1R, z1R) rotating coordinate system related to the driver gear wheel.

I am looking for the surface connected to K2R 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 K1R coordinate system of the driver wheel and the rotating K2R coordinate system of the driven wheel [4,9]:(3)Transformation matrices between the K1R and K2R coordinate systems (Fig. 5):(4) (5)

Taking into consideration the correlation between the velocity vectors of the relative movement in K1R and K2R coordinate systems [4,9]:(6)in K1R coordinate system, the relative velocity vector, based on (3) is:(7)where P1k 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 vector-scalar function simultaneously [4,9]:(9)

thumbnail 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, z1 number of teeth of the driver bevel gear wheel, z2 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:

(10)

  • half pitch angle of the driver gear wheel:

(11)

  • half pitch angle of the driven gear wheel:

(12)

  • effective pitch surface radius:

(13)

  • addendum on the largest diameter:

(14)

  • dedendum on the largest diameter:

(15)

  • the largest tip circle diameters:

(16)

  • the largest root circle diameters:

(17)

  • face width:

(18)

  • dedendum angle:

(19)

  • tip cone angles

(20)

  • root cone angles:

(21)

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 B-spline 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).

thumbnail Fig. 6

Profile curves of the designed bevel gear pair (m = 15 mm, z1 = 20, z2 = 30, Σ = 90°). (a) Profile curves of the driver gear. (b) Profile curves of the driven gear.
thumbnail Fig. 7

CAD model of the bevel gear pair having straight teeth (m  = 15 mm, z1 = 20, z2 = 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].

Table 1

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.

Table 2

Material properties.

thumbnail Fig. 8

FE mesh.

4.2 Loads and boundary conditions

Four coordinate systems have been defined for LTCAs: Ks − absolute static, Ks1–static related to the driver wheel, Ks2–static, related to the driven wheel and Kc– 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).

thumbnail Fig. 9

Setting loads and boundary conditions.
thumbnail 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.

thumbnail 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.
thumbnail Fig. 11

Continued.
thumbnail 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.

thumbnail 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.
thumbnail Fig. 13

Continued.
thumbnail 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 uP. 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.

thumbnail Fig. 15

Definition of displacement vector.
thumbnail 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.
thumbnail Fig. 16

Continued.
thumbnail 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 (Computer-Aided 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 K1R coordinate system

: Placement vector of the moving points of the profile curve

: Placement vector of moving point of conjugated surface

: Placement vector

: Normal stress vector (MPa)

: Shear stress vector (MPa)

: Differential operator

M2R,1R, M1R,2R : Translation matrices between K1R and K2R coordinate systems

P1k : Matrix of kinematic mapping

c * : Root clearance factor

: Displacement vector (mm)

σn : Normal stress (MPa)

2 ⋅ θ : Angle of planning (°)

A : Deformation tensor

b : Face width (mm)

d0 : Largest pitch circle diameter (mm)

df : Largest tip circle diameter (mm)

dl : Largest root circle diameter (mm)

f0 : Largest addendum (mm)

i21, i : Transmission

js : Clearance at flank (mm)

K1R (x 1R, y 1R, z 1R): Rotational coordinate system related to the driver gear wheel

K1S (x1S, y 1S, z 1S): Static coordinate system related to the driver gear wheel

K2R (x 2R, y 2R, z 2R): Rotational coordinate system related to the driven gear wheel

K2S (x 2S, y 2S, z 2S): Static coordinate system related to the driven gear wheel

l0 : Largest dedendum (mm)

m : Module (mm)

M: Load moment (Nm)

Re : Effective pitch surface radius (mm)

Sax : Pitch circle face width on the largest diameter (mm)

t : Circular pitch on the largest pitch circle diameter (mm)

z1, z2 : Number of teeth

α0 : Angle of contact (°)

δ0 : Pitch angle (°)

δf : Tip cone angle (°)

δl : Root cone angle (°)

εn : Normal elastic strain (mm)

εx, εy, εz : Elastic strain values per unit length (mm)

λ : Dedendum angle (°)

μ : Friction coefficient

Σ : Shaft angle (°)

ϕ : Addendum angle (°)

ϕ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)

x, y, z : Coordinates (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 EFOP-3.6.1-16-2016-00022 project, and it is co-financed by the European Union and the European Social Fund.

References

  1. D.W. Dudley, Gear Handbook, MC Graw Hill Book Co., New York, 1962 [Google Scholar]
  2. R. Gołębski, A. Szarek, Diagnosis of the operational gear wheel wear, Technical Gazette 26, 658–661 (2019) [Google Scholar]
  3. J. Klingelnberg, Bevel Gear, Fundamentals and Applications, Springer, Berlin, 2016 [CrossRef] [Google Scholar]
  4. F.L. Litvin, A. Fuentes, Gear Geometry and Applied Theory, Cambridge University Press, Cambridge, 2004 [Google Scholar]
  5. V. Rohonyi, Fogaskerékhajtások. Műszaki Könyvkiadó, Budapest, 1980 [Google Scholar]
  6. Z. Terplán, Gépelemek IV., Kézirat, Tankönyvkiadó, Budapest, 1975, p. 220 [Google Scholar]
  7. 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]
  8. I. Popa-Müller, Generation of eloid bevel-gear on tooth-cutting machine with the generating face-gear, 15th International Conference in Mechanical Engeneering, Cluj-Napoca, 2007 [Google Scholar]
  9. F.L. Litvin, Gear geometry and applied theory, Prentice Hall, Englewood Cliffs, NJ, 1994 [Google Scholar]
  10. A. Fuentes, J.L. Iserte, I. Gonzalez-Perez, F.T. Sanchez-Marin, Computerized design of advanced straight and skew bevel gears produced by precision forging, Comput. Methods Appl. Mech. Eng. 200, 2363–2377 (2011) [Google Scholar]
  11. I. Juhász, Görbékés felületek modellezése, Miskolci Egyetem, elektronikus jegyzet, p. 114 [Google Scholar]
  12. 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]
  13. 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]
  14. F.L. Litvin, A. Fuentes, I. Gonzales-Perez, L. Carnevali, T.M. Sep, New version of Novikov-Wildhaber helical gears: computerized design, simulation of meshing and stress analysis, Comput. Methods Appl. Mech. Eng. 191, 5707–5740 (2002) [Google Scholar]
  15. F.L. Litvin, A. Fuentes, I. Gonzalez-Perez, 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]
  16. 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]
  17. I. Popa-Müller, Case-study and simulation by meshing of bevel-gears, 14th International Conference in Mechanical Engeneering, Tg-Mures, 2006 [Google Scholar]
  18. I. Popa-Müller, Simulation of Generation by Meshing of Conical Gear Octoid II with Hobbing Machine Heidenreich-Harbeck, 18th International Conference in Mechanical Engeneering, Baia Mare, 2010 [Google Scholar]
  19. J.K. Theodore, Tooth contact analysis of spiral bevel and hypoid gears under load, SAE Trans. 90, 2205–2216 (1981) [Google Scholar]
  20. 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]
  21. W. Wagner, S. Schumann, B. Schlecht, Co-simulation of the Tooth Contact of Bevel Gears within a Multibody Simulation, Forsch Ingenieurwes 83, 425–433 (2019) [CrossRef] [Google Scholar]
  22. I. Kozák, Gy. Szeidl, Fejezetek a szilárdságtanból, Kézirat, 2008–2012, Miskolci Egyetem, p. 284. elektronikus jegyzet [Google Scholar]
  23. S. Moaveni, Finite Element Analysis, Theory and Application with ANSYS, Pearson Education Limited, London, 2015, p. 928 [Google Scholar]
  24. 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, Computer-aided design and loaded tooth contact analyses of bevel gear pair having straight teeth by different loaded torques, Mechanics & Industry 21, 109 (2020)

All Tables

Table 1

Parameters of the designed bevel gear pair.

In the text
Table 2

Material properties.

In the text

All Figures

thumbnail Fig. 1

Angles and diameters of bevel gear.
In the text
thumbnail Fig. 2

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

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

Definition of planing angle.
In the text
thumbnail Fig. 5

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

Profile curves of the designed bevel gear pair (m = 15 mm, z1 = 20, z2 = 30, Σ = 90°). (a) Profile curves of the driver gear. (b) Profile curves of the driven gear.
In the text
thumbnail Fig. 7

CAD model of the bevel gear pair having straight teeth (m  = 15 mm, z1 = 20, z2 = 30, Σ = 90°).
In the text
thumbnail Fig. 8

FE mesh.
In the text
thumbnail Fig. 9

Setting loads and boundary conditions.
In the text
thumbnail Fig. 10

Defining normal stress.
In the text
thumbnail 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
thumbnail Fig. 11

Continued.
In the text
thumbnail Fig. 12

Load − normal stress diagram.
In the text
thumbnail 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
thumbnail Fig. 13

Continued.
In the text
thumbnail Fig. 14

Load − normal elastic strain diagram.
In the text
thumbnail Fig. 15

Definition of displacement vector.
In the text
thumbnail 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
thumbnail Fig. 16

Continued.
In the text
thumbnail Fig. 17

Load–normal deformation diagram.
In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.