Extended calculation method for determining the pitting load carrying capacity of bevel and hypoid gears

The flank load carrying capacity of bevel and hypoid gears is mainly limited by the failure modes pitting, scuffing, tooth flank fracture and the phenomenon micropitting. By application of a standardized calculation method, e.g. according to the international standard ISO 10300:2014, a first estimation of the flank load carrying capacity can be made based on the macro geometry of the bevel or hypoid gear set. According to method B of ISO 10300:2014 the complexity of the real geometry of bevel and hypoid gears is reduced to a virtual cylindrical gear geometry. The load carrying capacity regarding scuffing, micropitting and tooth flank fracture can be determined by using the virtual cylindrical gear geometry along the path of contact. However, the determination of the pitting load carrying capacity is carried out on a single representative point on the path of contact of the virtual cylindrical gear. This paper shows an extended calculation method for the determination of the pitting load carrying capacity of bevel and hypoid gears along the path of contact of the virtual cylindrical gear geometry. Due to the calculation along the path of contact the extended method allows a more precise estimation of the pitting load carrying capacity than the current standard calculation method ISO 10300-2:2014 using the same input data. Within this paper all relevant factors of the extended calculation method are explained in detail. Furthermore, the verification of the extended calculation method with calculation results of an intense validated loaded tooth contact analysis, corresponding to method A of ISO 10300-2:2014, is presented.


Introduction
If rotary movements and torque have to be transmitted angularly in a drive train, bevel gears are of particular importance. A cone represents the basic body of bevel gears (see Fig. 1a). A bevel gear stage consists of an input gear (pinion) and an output gear (wheel). The input gear usually has a fewer number of teeth than the output gear [2]. Due to their tapered basic body, various shaft angles can be realized. The most commonly used shaft angle is 90°. An axis offset (also called hypoid offset) is also possible in order to optimally adapt the gear set to the respective application (see Fig. 1b). Due to the geometric flexibility regarding shaft angle and hypoid offset, bevel gears have a very wide range of applications. A typical large field of application are axle drive systems for passenger cars or commercial vehicles. But also, bevel and hypoid gears are widely used in industrial, marine and aeronautical applications. In the following, within the term bevel gear hypoid gears are included. The maximum transferable torque is mainly limited by the load carrying capacity of the gear set. Typical failure modes affecting the flanks of a bevel gear set are pitting, tooth flank fracture and scuffing. Furthermore, the phenomenon micropitting also affects the load carrying capacity of the tooth flank since it may change the load distribution on the flank. Fig. 2 shows the above-mentioned failure modes on a bevel gear flank without offset and with a shaft angle of 90°.
In order to ensure a reliable operation of a bevel gear set, the load carrying capacity of the tooth flank has to be taken into account within the gear design process. Due to the complex geometry of a bevel gear a virtual cylindrical gear is used as substitute gearing according to the current state of the art [6]. The substitute gearing represents the meshing conditions of a bevel gear with sufficient accuracy. Therefore, it is possible to determine the load carrying capacity of bevel gears regarding pitting, scuffing, micropitting and tooth flank fracture in an early stage of design by use of few, reliable input variables.
The load carrying capacity regarding scuffing, micropitting and tooth flank fracture can be determined by use of the virtual cylindrical gear geometry along the path of contact [5,[7][8][9][10][11]. In contrast to this the determination of the pitting load carrying capacity is only carried out on a single, representative point on the path of contact of the virtual cylindrical gear [12].
Within this paper a new extended calculation method is provided to calculate the pitting load carrying capacity along the entire path of contact of the virtual cylindrical gear, analogous to the current calculation methods for determining scuffing, micropitting and tooth flank fracture load carrying capacities.
In scope of this paper calculation methods comply in assumptions, structure and conditions with the standard calculation method B of ISO 10300:2014 are termed as calculation methods of standard category B. Analogous calculation methods comply with the standard calculation method A of ISO 10300:2014 are termed as calculation methods of standard category A.

Failure mode pitting
The failure mode pitting occurs if the material strength is exceeded by the occurring stress [13,14]. The damage usually starts at or near the gear surface, preferred on areas with negative specific sliding [15][16][17]. The crack is initiated at or near the tooth flank surface and grows into the depth. Consequences are shell-shaped material eruptions [18], which can lead to unfavorable NVH-behavior and additional dy-  [19]. Pitting often shows a progressive damage development and may cause other failure modes such as tooth breakage. Fig. 2 shows a typical pitting on a bevel gear flank.

Calculation methods according to the state of the art
The pitting safety factor SH of gears is defined as the ratio of the permissible contact stress σHP and the contact stresses σH [12,20].
Based on experimental investigations, Wirth [21] developed a local calculation method of standard category A and additionally a calculation approach of standard category B to determine the pitting load carrying capacity of bevel gears (see Fig. 3). By use of the local calculation method, a local safety against pitting on the tooth flank is determined, using the exact description of the macro and micro geometry of a bevel gear set. For calculating the load carrying capacity according the calculation method of standard category B only information about the macro geometry is necessary. The standard-capable method by Wirth [21,22] is to a large extend adopted in ISO 10300-2:2014 method B1 [12].
Klein [8,23] extended the calculation method of standard category B by Wirth [21,22] to the effect of a greater resolution of the stress along the path of contact of the virtual cylindrical gear. Klein [8,23] concentrates his researches on scuffing for bevel and hypoid gears and in 2021 his findings were documented as an international technical specification [9]. Figure 3 illustrates the relationship between the calculation methods described above.
In the following subchapters, the calculation methods that constitute the state of the art and form the basis of the new extended calculation approach are discussed in more detail.

Local calculation method
The local calculation method of standard category A is based on the real bevel gear geometry (macro and micro geometry) and uses local contact stresses to determine the safety factors against pitting. The local contact stresses σH,i can be calculated by the multiplication of the calculated Hertzian stresses σH,lok,i and the square root of the dynamic factor KV according to ISO 10300-1:2014 [6] in order to take load increases due to internal dynamic effects into account. The local contact stresses σH,lok,i are usually determined using a loaded tooth contact analysis (LTCA), BE-CAL [24] was used here for this purpose [25].
The local permissible contact stress σHP,i is calculated on the basis of the strength numbers according to ISO 6336-5:2016 [26] using the corresponding equations according to ISO 10300-2:2014 [12], whereas the bevel specific influence on the permissible contact stress is taken into account by the additional slip factor ZS [25].
The life factor ZNT, the size factor ZX, the lubricant factors ZL ZV ZR and the work hardening factor ZW, determined according to ISO 10300-2:2014 [12], are considered constant over the entire tooth mesh. The hypoid factor ZHyp,i as well as the slip factor ZS,i are local values.
The hypoid factor ZHyp,i takes into account the influence of longitudinal sliding on the pitting load carrying capacity. It is based on the sum of velocities vertical to the contact line v˙; vert;i and the sliding velocity parallel to the contact line v g;par;i . [25].
The slip factor ZS,i considers the influence of the slip on the permissible contact stress. It is determined for pinion and wheel by the slip perpendicular to the contact line ζvert,i [25].
Recalculation of experimental tests showed a very good correlation between the predicted and the real pitting areas on the damaged gear flanks, see Fig. 9. Not only the place of the initial pitting is calculated correctly, also the pitting growth is predicted reasonably [21].
The local calculation method is applicable for case hardened gears which meet the state-of-the-art regarding material and heat treatment and represents the reference for the development of the calculation methods of standard cate-gory B. It can only be applied in combination with a loaded tooth contact analysis based on the real tooth flank geometry (macro and micro-geometry).

Standard calculation method
For practical reasons also a calculation method of standard category B was developed by Wirth [22]. This approach allows an estimation of the pitting load carrying capacity of a bevel or hypoid gear set at a very early design stage of the gear development. Due to the complex geometry of bevel and hypoid gears, the calculation method is based on a virtual gear geometry. The bevel gear geometry is transferred to a virtual tooth system, which represents the tooth contact conditions and thus forms the basis for determining the load carrying capacity regarding different failure modes.
Initially, a virtual cylindrical gear for bevel gears without offset and a virtual crossed helical gear or a virtual bevel gear for hypoid gears were used [27][28][29]. Wirth [21] developed a virtual cylindrical gear, which can be used for bevel and hypoid gears. The virtual cylindrical gear also forms the basis of the international standard ISO 10300 series and represents the state-of-the-art for calculation methods according to the standard category B. The virtual cylindrical gear geometry is structured according to the schematic diagram in Fig. 4. The virtual cylindrical gear is derived from the reference cone of the bevel gear or rather the hypoid gear. Detailed information regarding the representation of a bevel or hypoid gear as a virtual cylindrical gear are available in Annex A of ISO 10300-1:2014 [6].
A comparable plane of action is essential for a reasonable calculation of the load carrying capacity by using the virtual cylindrical gear. Wirth [21] has shown that the zone of action can be described with sufficient accuracy by a parallelogram. Fig. 5 shows the schematic representation of the  [30] zone of action according to Wirth [21]. More detailed information on determining the relevant zone of action for bevel gears can be found in ISO 10300-1:2014 [6]. Based on this, the load distribution between several meshing teeth is determined by the length ratios of the contact lines lb in the zone of action. As shown schematically in Fig. 5, a parabolic load distribution along the line of action and an elliptical load distribution along the contact lines lb are assumed. In total, a maximum of three simultaneously meshing teeth are considered.
Based on the virtual cylindrical gear set and the operating conditions, a global pitting safety factor according to ISO 10300-2:2014 [12] or Wirth [21] can be determined. The nominal contact stress σH0 can be calculated with Eq. 6 according to ISO 10300-2:2014 [12], where Fn is the nominal force of the virtual cylindrical gear at the mean point P. The formula of the nominal contact stress σH0 according to Wirth [21] is given in Eq. 7.
According to ISO 10300-2:2014 [12] and Wirth [21], the pitting load carrying capacity is determined at one decisive point, which lies between the inner point of single contact and the midpoint of the path of contact, depending on the face contact ratio of the virtual cylindrical gear. The mid zone factor ZM-B is used to convert the radii of curvature at the design point P to the radii of curvature at the point which is decisive for the failure mode pitting.
The contact stress can be calculated by use of Eq. 8.
To avoid pitting, the occurring contact stress σH has to be lower than the permissible contact stress σHP, which can be calculated according to ISO 10300-2:2014 [12] as given in Eq. 9 and according to Wirth [21] as specified in Eq. 10.
Analogous to the local calculation method, Wirth [30] uses an additional slip factor ZS1,2, given in Eq. 11, to determine the permissible flank stress σHP, as quoted in Eq. 10.
Finally, the occurring contact stress calculated according to Eq. 8 and the permissible contact stress, calculated according to Eq. 9, can be used for determining a global safety value against pitting according to Eq. 1.
The essential difference between the standard calculation method ISO 10300-2:2014 [12] and the calculation method of standard category B provided by Wirth [21] regards the consideration of the positive and negative sliding zones on the tooth flank. Within ISO 10300-2:2014 [12] the slip factor ZS is not used in the calculation method whereas the bevel factor ZK is used to considerate among other differences the specific sliding conditions for bevel gears in comparison to cylindrical gears. Subsequently for the standard-capable calculation method according to Wirth [21], the permissible contact stress is increased by the slip factor ZS. In the standard calculation method ISO 10300-2:2014 [12] the nominal contact stress is decreased by the bevel gear factor ZK. In sum both considerations of the slip conditions resulting in nearly the same safety factors. On a physical point of view the increasing of the permissible contact stress σHP is more comprehensible than the reduction of the nominal contact stress σH0 due to the slip conditions on the tooth flank. This issue is discussed in more detail in Sect. 6.2.

Extended standard calculation method to determine the Hertzian contact stresses along the path of contact
Based on the virtual cylindrical gear geometry the Hertzian contact stress can be calculated along the path of contact of the virtual cylindrical gear according to ISO/DTS 10300-4:2019 [10], ISO/TS 10300-20:2021 [9] or Klein [8,31]. Therefore, the transverse path of contact between pinion and wheel of the virtual cylindrical gear set is divided in a number of sections (e.g. n = 10) to determine the local safety factors for n + 1 points of contact. To define the contact points the coordinate gY is established, which originates in the pitch point C (= design point P of the bevel gear set). Towards the pinion tip gY is defined as positive, towards the pinion root it is defined as negative. Fig. 6 shows an exemplary path of contact including the coordinate gY. In the Fig. 6 Transverse path of contact [30] boundary points A and E on the transverse path of contact gY is determined as follows [9].
The length of the transverse path of contact can be subdivided in a number of sections i. For a contact point Y on the transverse path of contact of the virtual cylindrical gear the corresponding coordinate gY(Y) can be calculated as follows [9]. Where: g v˛t he length of path of contact of the virtual cylindrical gear.
In the next step the length of each contact line lb,Y can be calculated at each contact point Y of the virtual cylindrical gear by using Eq. 15. Fig. 7 shows the general definition of the length of the contact lines lb,Y according to ISO/TS 10300-20:2021 [9]. The theoretical length of a contact line at a contact point lb0,Y of the virtual cylindrical gear can be calculated by use of the following equation [9].
The contact line lengths are trimmed at the edges by the elliptical correction function Clb [9].
Based on the contact line lengths and the assumption of load distribution shown in Fig. 5, a local load distribution factor ZLS,Y can be determined. For this purpose, the area of the semi-ellipse of the middle contact line is set in relation to the total area consisting of the semi-ellipse areas over all contact lines [9].
Using the curvature factor XY the local equivalent radius of relative curvature ρrel,Y can be calculated for each The conversion of the maximum line load at point M to the contact point Y is calculated by use of the local face load factor KHβ,Y, according Eq. 23-25 [9].
K Fig. 8 Visualization of the connections between the calculation methods regarding the load carrying capacities of bevel gears discussed in this paper The modified Hertzian contact stress at contact point Y H;mod;Y can be calculated based on the Hertzian contact stress σH,Y. The Hertzian contact stress σH,Y is calculated by use of the described load sharing and curvature factors according Eq. 27. Within the determination of the modified Hertzian contact stress σH,mod,Y the profile crowning of the gear set is taken into account by use of the exponent e according to ISO 10300-2:2014 [9, 12] Fig. 9 Comparison of predicted safety factors against pitting on the tooth flank and real pitting damage occurred by experimental investigations [1]. a Pitting on the tooth flank of the pinion after 3 10 6 load cycles at pinion [1]. b Calculation result of the safety factor against pitting according to local calculation method

Description of the methodical approach
In the following the calculation methods used within this paper and their relation between each other are presented to obtain a better overview about the discussed calculation methods.
The suggested extended calculation method of the load capacity against pitting of bevel gears combines the localized calculation method of the contact stresses given in ISO/TS 10300-20:2021 [9] and the calculation method by Wirth [21] enhanced by localized calculation factors. Therefore, within the presented extended calculation method, the load capacity against pitting can be calculated along the path of contact of the virtual cylindrical gear, similar to the given standard of ISO/TS 10300-20:2021 [9]. Fig. 8 shows graphical the connection between the calculation methods.
As visualized in Fig. 8 the local calculation method of standard category A, provided by Wirth [21], is used to validate the new extended calculation method. This approach is reasonable due to the fact, that the local calculation method is based on experimental investigations and is able to predict the areas on the tooth flank with a safety factor against pitting less than 1.0. In Fig. 9 the red dots represent the calculation points with a safety factor against pitting less than 1.0, which match very good with the early pitting damage on the related tooth flank. Additionally, the local calculation method is widely, successfully used with application in practice. Therefore, the validation of the extended calculation method of standard category B based on a proven calculation method of standard category A is considered to be sufficient.

New extended calculation approach for the pitting load carrying capacity of bevel and hypoid gears
Based on extensive experimental studies within the research projects by Wirth [22] and Boiadjiev [3] as well as the findings of Klein [8], a new extended approach of standard category B for calculation of the pitting safety factor for bevel gears was developed by one of the co-authors. This approach is described in the following.
The new extended calculation approach offers the possibility to determine the pitting load carrying capacity of bevel gears over the entire length of the path of contact of the corresponding virtual cylindrical gear. Thus, purely on the basis of the macro geometry, the user has more information about the pitting load carrying capacity of the tooth flank compared to the current standard at an early stage of gear design. For this purpose, the Hertzian contact stresses as well as the permissible contact stresses have to be determined along the path of contact of the virtual cylindrical gear pair.
To determine the modified Hertzian contact stresses σH,Y,mod along the path of contact of the virtual cylindrical gear the calculation method according to Klein [8] and ISO/TS 10300-20:2021 [9] is applied.
Unlike to the known calculation methods of standard category B according to Wirth [21] as well as according to ISO 10300-2:2014 [12], the safety factor against pitting is not determined at a decisive point, but locally for each considered point Y on the path of contact of the virtual cylindrical gear. Based on Wirth [21] the permissible contact stress is determined according to Eq. 28 with the extended, localized hypoid factor ZHyp,Y, slip factor ZS,Y and speed factor ZV,Y for each point of contact Y along the path of contact considered.
The life factor ZNT, the size factor ZX, the work hardening factor ZW, the lubricant factor ZL and the roughness factor ZR represent global values, that are constant along the path of contact. They are calculated according to ISO 10300-2:2014 [12]. The speed factor ZV,Y, hypoid factor ZHyp,Y and slip factor ZS,Y depend on the local velocity conditions at the considered point of contact Y along the path of contact and have to be modified in comparison to the equations given in ISO 10300-2:2014 [12].
The speed factor ZV,Y depends on the tangential speed at mid face width of the reference cone of the wheel vmt2 and can be calculated by using Eq. 29.
The local tangential speed vmt2,Y can be calculated as a function of the coordinate for transverse path of contact gY.  Table 1). The speed factor shows an almost constant behavior along the path of contact of the virtual cylindrical gears. The value of the speed factor depends on gear size and the field of action. Although the speed factor for the test gearings used here show only very small fluctuations in value over the path of contact, the speed factor is used in the new extended calculation method in localized form according to Eq. 29 to establish a conclusive calculation method. For gears of larger sizes, greater fluctuations in the speed factor over the path of contact are expected.
By use of the hypoid factor ZHyp,Y the influence of longitudinal sliding on the pitting load carrying capacity is taken into account. The hypoid factor ZHyp,Y can be calculated for each considered contact point with Eq. 32: With: It is assumed that the sliding velocity parallel to the contact line vg,par shows constant values along the whole path of contact [5,8,9,21]. Therefore, the sliding velocity vg,par is determined at the mean point P by use of the following equations.
The absolute local surface velocities are required for the determination of the local sum of velocities vertical to the contact line v˙; vert;Y . It should be noted that hypoid gears have an additional velocity component in the tooth width direction due to the hypoid offset. However, it can be assumed that this is constant over the line of contact [5,8,9]. The surface velocity in lengthwise direction of the tooth wt1,2s is determined at the mean point P. The surface velocity in profile direction of the tooth is calculated locally over the path of contact of the virtual cylindrical gear.
The total surface velocity can be determined by vectorially addition of the velocity components in lengthwise and profile direction of the tooth. Using the inclination angle of contact line βB, the absolute surface velocity vertical to the contact line can be calculated.
The local sum of velocities vertical to the contact line v˙; vert;Y consists of the local surface velocity of pinion and wheel oriented vertically to the contact line.
v˙; vert;Y = w t 1vert;Y + w t 2vert;Y (43) The localized slip factor ZS,Y considers the influence of the specific slip conditions on the tooth flank on the permissible contact stress. It is defined according to the local calculation method.

Validation of the new extended calculation approach
For validation of the new extended calculation approach of standard category B, a comparison between the local calculation method of standard category A and the new extended calculation approach is presented. For this purpose, the same test gear sets G0, G15, G31.75 and G44 used by Wirth [21,22] are considered. The basic geometry data of these investigated gear sets are given in Table 1. All test gears show a typical ease off design and were ground and blast cleaned after heat treatment. The recalculations were carried out for a torque of 500 Nm and a rotational speed of 4500 min -1 at pinion. As described, the extended calculation approach is based on the virtual cylindrical gear geometry. The virtual cylindrical gear geometry is calculated in the mean section of the bevel gear. In Fig. 11 the calculation points along the Fig. 11 Definition of tooth sections and path of contact gY [32] path of contact are illustrated exemplarily as red dots for the proposed extended calculation approach. Unlike the extended calculation, the local calculation considers the whole flank of the gear by use of tooth sections, which are distributed over the whole flank. The tooth sections are shown as black dashed lines in Fig. 11. The following comparisons between the results of the extended calculation approach, the standard calculation method according to ISO 10300-2:2014 [12] and the local calculation method (method A) are performed along the mean section.

Comparison of radius of relative curvature
Within the newly developed calculation method, the radii of relative curvature are calculated in accordance with ISO/TS 10300-20:2021 [9] at the considered contact points Y of the virtual cylindrical gear. Exemplary values of the local radius of relative curvature along the path of contact are shown in Fig. 12. The dashed lines in Fig. 12 represent the values determined within a tooth contact analysis by use of the software BECAL [24]. The solid lines show the values calculated by the extended approach according to ISO/TS 10300-20:2021 [9], which is based on the approach of Shtipelman [33]. The results of the extended calculation approach based on the virtual cylindrical gear geometry correlate well with the results of the tooth contact analysis based on the real flank geometry. The small deviations can be explained by the fact, that the microgeometry of the bevel gears cannot be represented in detail by the virtual cylindrical gear geometry.

Comparison of Hertzian contact stresses
The modified Hertzian contact stress σH,Y,mod represents a significant factor on the load carrying capacity of gears regarding pitting. Therefore, the correlation between the new extended approach and a LTCA was investigated. In Fig. 13 [12].
All in all, the Hertzian stresses calculated by the extended approach of standard category B and the LTCA of standard category A are in good accordance. However, it can be seen that the geometry of the test gears can only be represented roughly, but sufficiently accurate, by use of a virtual cylindrical gear geometry. The extended calculation approach shows slightly higher stresses than the local calculation method for all considered examples. In summary, it can be stated that the assessment of the risk of pitting failure is conservative if the new extended approach of standard category B is used for the herein investigated test gear sets. For a more accurate consideration of the real flank geometry the described local calculation approach of standard category A should be used. In Fig. 13 it is visible that the Hertzian contact stress according to the standard calculation method ISO 10300-2:2014 [12] is located at a significant lower level in comparison to the extended approach of standard category B and the LTCA. This lower level can be traced back to the consideration of the bevel gear factor ZK within calculation method of ISO 10300-2:2014 [12], which is not considered within other given calculation methods.
With Fig. 13 the given statement in Sect. 3.2.2, that the calculated Hertzian contact stress σH0 by use of the slip factor ZS is a better representation of the reality than by use of the bevel gear factor ZK, can be underlined. This can be seen in the better correlation of the calculated Hertzian contact stresses σH0 by use of the extended calculation method with the results of the LTCA.

Comparison of hypoid factor
In comparison with ISO 10300-2-2014 [12], the newly developed calculation method is based on a hypoid factor determined along the path of contact of the virtual cylindrical gear. In Fig. 14, the hypoid factor determined according to the new extended calculation approach of standard category B is compared with the hypoid factor determined using the local calculation method of standard category A.
The dashed lines represent the values determined within local calculation method. The solid lines show the values calculated by use of the new extended approach. The hypoid factor is ZHyp = 1.0 for bevel gears without hypoid offset. In case of hypoid gear sets, the trend of the increasing hypoid factor along the path of contact, calculated according to the extended calculation method, corresponds to the results of the local calculation method. However, with increasing hypoid offset the gradient of the curve of the hypoid factor differs slightly between the extended and local calculation method. In summary, it can be stated that with increasing hypoid offset, the velocity conditions are more difficult to map using the extended approach that does not consider the microgeometry. Nevertheless, the local hypoid factor can be calculated with sufficient accuracy by use of the extended approach.

Comparison of permissible contact stress
The permissible contact stress σHP,Y is a decisive influence factor within determination of the pitting load carrying capacity. Fig. 15 shows a comparison of the permissible contact stress calculated by use of the new extended calculation method of standard category B and local calculation method of standard category A. The dashed lines represent the results of the local calculation method and the solid lines the results of the extended calculation approach. Since the pinion and the wheel show different permissible contact stresses, the calculation results regarding the pinion are shown in black and regarding the wheel in grey. As a supplement to the given calculation results the permissible contact stress according to the standard calculation method ISO 10300-2:2014 [12] is given in following figures by dot- The application of the new extended calculation approach results in slightly lower permissible contact stresses, compared to the results of the local calculation method. Furthermore, the tendency to increase the permissible contact stress along the path of contact is well reflected by the extended calculation approach, especially for hypoid gears. This behavior can be explained by the local calculation of the hypoid factor, slip factor and the speed factor. However, the presented newly developed extended calculation approach aims at determining the pitting load carrying capacity in an early design stage of the gear set without using a complex tooth contact analysis. From this point of view, the shown accuracy is sufficient for a reliable prediction of pitting of bevel gears in an early design stage and is conservative.
Due to the consideration of the bevel slip factor ZS within calculation methods of the extended approach and the loaded tooth contact analysis (LTCA) the permissible contact stresses of these calculation methods lie mostly above the calculated contact stresses according to the standard calculation method ISO 10300-2:2014 [12]. It can now be seen that the quoted differences between the Hertzian contact stresses and permissible contact stresses, based on the consideration of the bevel gear factor ZK or the bevel slip factor ZS, nearly compensate each other within calculation of safety factor against pitting. The extended calculation approach of standard category B determines slightly lower safety factors against pitting compared to the results of the local calculation method of standard category A for all considered examples and is therefore on the safe side. Thus, the newly developed calculation approach has the character of a standard-capable calculation method of standard category B, which should always guarantee safer evaluation due to lack of input information compared to the local calculation method. All in all, the newly developed calculation approach is able to assess the pitting load carrying capacity of bevel and hypoid gears in an accurate way for the considered gear designs.

Conclusion
In this paper a new extended calculation approach is shown, which offers the possibility to determine the pitting load carrying capacity of bevel and hypoid gears along the path of contact of the corresponding virtual cylindrical gear. Thus, in an early stage of gear design more information about the pitting load carrying capacity along the mid-section of bevel and hypoid gears is generated compared to the current standard, which only enables the calculation of the pitting load carrying capacity at one decisive point on the tooth flank. Additionally, a consistent calculation approach is developed with regard to ISO/TS 10300-20:2021 [9] (scuffing), ISO/DTS 10300-4:2019 [10] (tooth flank fracture) and Hombauer [5] (micropitting). The new extended calculation method was successfully validated on basis of a theoretical study using local calculation results as a reference.
Further researches regarding the phenomenon, that the gear, despite an equal safety factor against pitting as for example given in Fig. 16, does not fail as often as the corresponding pinion were done by Reimann [1] at the Gear Research Center (FZG). Initial experimental investigations see an influence of the driven and drive part of the gear set. The consideration of this investigated influence in calculation methods regarding the load carrying capacity against pitting is content of further releases of the Gear Research Center (FZG).

Parameters
The parameters are shown in Table 2.

Indices
The indices are shown in Table 3. Is the maximum length of path of contact of virtual cylindrical gear in transverse section mm h Is an auxiliary value mm l b Is the length of contact line mm l b0 Is the theoretical length of contact line mm m mn Is the mean normal module mm n Is the rotation speed min -1 p Is the relative peak load for calculating the load sharing factor v g Is the sliding velocity in mean point P m/s v g;par Is the sliding velocity parallel to the contact line m/s v mt Is the tangential speed at mid face width of the reference cone m/s v˙; vert Is the sum of velocities vertical to the contact line m/s w t Is the surface velocity m/s w ts Is the surface velocity in lengthwise direction m/s w th Is the surface velocity in profile direction m/s w t;vert Is the surface velocity vertical to the contact line m/s x Y Is the coordinate of the end of the contact line at contact point Y mm y Y Is the coordinate of the end of the contact line at contact point Y mm z Is the number of teeth z Y Is an auxiliary value at the contact point Y mm A Is the related area for calculation the load sharing factor mm C lb Is the correction factor for the length of contact lines -F n Is the nominal normal force N K A Is the application factor -K H˛I s the transverse load factor for contact stress -K HˇI s the face load factor for contact stress k s Is the auxiliary variable to exclude tip and toot boundary points for bevel gears with mean spiral angle zero .ˇm = 0/ -K V Is the dynamic factor -S H Is the safety factor for contact stress (against pitting) -S Hmin Is the minimum safety factor for contact stress -T D Is the permanently transmittable torque Nm Z LS Is the load distribution factor -Z Y Is the curvature factor -Z E Is the elasticity factor -Z Hyp Is the hypoid factor - Is an auxiliary angle°! wt Is the angle between the sliding velocities in lengthwise and tooth profile direction°T Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4. 0/.