Abstract
Beveloid gears are usually employed to transfer motion of small shaft angle, but the assembly errors or misalignment can affect the gear tooth contact and dynamic performance. A dynamic model and tooth contact analysis model of the intersecting beveloid gear pair are established, the defects of the theoretical tooth surface in the meshing are analyzed, and a microgeometry modification method to improve the meshing performance is proposed. At the same time, the effects of system misalignment and manufacturing errors on its meshing characteristics are considered. By studying the sensitivity of the static and dynamic characteristics of beveloid gear on different misalignment components, the influence laws of misalignment on the contact pattern, transmission error, mesh stiffness, and dynamic excitation force are obtained, which give suggestions for tooth surface modification and support structural stiffness design of beveloid gear with a good NVH performance.
1. Introduction
Beveloid gears are generalized involute gears in which tooth tip and root have a cone angle, tooth thickness is taper, and the tooth flank is involute helicoid. Each transverse section profile has the same base circle, but their profile shift coefficient along the axis has changed linearly. Many AWD vehicles’ transfer cases usually have a beveloid gear pair with a small intersecting shaft angle, the driving pinion is a helical gear, the driven gear is a beveloid gear, and the tooth numbers are equal. The transmission schematic diagrams are shown in Figure 1. Transfer case with beveloid gear can avoid a structural obstacle and split the engine power into front axle power and final drive power. In addition, it has a smaller volume and light weight relative to a single or maybe two-stage cylinder gear transfer case.

(a)

(b)
Merritt and Beam proposed the theory of beveloid gearing, respectively. In 1954, beveloid gear was applied widely until the general theory of the design, processing, and measurement of beveloid gear was elaborated by Mitome [1, 2] in 1983. Beveloid gear pair, which has good backlash compensation characteristics, can be used to replace the cylindrical gear in precision parallel-axis transmission and can also be used in the intersecting or crossed shaft transmission with a small angle. Now beveloid gear is widely used in the transfer case of all-wheel-drive vehicles [3] and in the propulsion system of speedboats, as well as in marine transmission [4, 5].
Many studies on geometric design and machining methods for beveloid gears have been proceeding during the last decades. Liu and Tsay [6] investigated the condition of tooth toe undercutting and put forward two methods to avoid undercutting. Innocenti studied the hobbing kinematics parameters [7] and set up the theory of measurement over spheres of beveloid gears by simplifying its algebra [8]. Brecher [9] introduced a simulation-based method for function-oriented tolerancing of tooth flank modifications of beveloid gears.
In the aspect of tooth contact analysis (TCA) for beveloid gear pair, Mitome [10] analyzed the meshing characteristics of the theoretical tooth flank with an intersecting shaft, and the results show that it does not affect the accuracy of the transfer motion by point contact of the theoretical tooth surface with assembly errors. Meanwhile, Mitome studied the load capacity of beveloid gears based on Hertz theory. Innocenti [11] investigated the dependence of the backlash on the gear geometry and the shaft relative position. Liu and Tsay [12] analyzed the engaging characteristics of the theoretical tooth surface of nonparallel shaft beveloid gears by the TCA method and acquired the same conclusion as reference [10].
The theoretical tooth flank of beveloid gear is an involute helicoid, which is meshing as a point contact; with the increase of loading, the contact point will broaden into an ellipse, but the contact area is small, contact stress is large, and it is easily generated edge contact, which is very undesirable to the gear contact fatigue life. In order to improve these disadvantages, the modification tooth flanks usually are employed instead of the theoretical tooth flanks. Liu and Tsay [13] established a mathematical model of a concave tooth surface and conducted the tooth contact simulation; the tooth surface was a modification to a concave shape lengthwise relative to the theoretical tooth surface so that the instantaneous contact ellipse will be elongated and the contact area will be broader.
Brecher et al. took a simulation for grinding processing and tooth contact analysis of beveloid gears [14]; in addition, they investigated the influence relationships of machining method and motion performance [15]. Zhu and Song have a lot of valuable research on geometry design and tooth contact properties of marine crossed-axis beveloid gears [16–18]. Meanwhile, Song et al. analyzed the meshing characteristics and dynamic response of loaded beveloid gears with small shaft angle by parametric method [19, 20], and further on, they investigated the sliding friction effect on the gear dynamics and transmission error (TE) [21]. Trong Phu Do [22] established the elastic models for different types of beveloid gear and simulated contact force and TE in the process of tooth meshing. Until now, there are few works [23–25] on the static and dynamic mesh characteristics of beveloid gearing, which include structure deformation by load, assembly error, and tooth modifications; especially, there is a scarce study on the mesh properties of modification concave tooth surface under misalignment.
Although there has been a lot of research on the meshing characteristics of marine beveloid gears, the load condition of marine gears in water is not as complex as that of AWD off-road vehicles; because AWD off-road vehicles will bear particularly complex road conditions, especially under harsh conditions such as rapid acceleration, deceleration, and high-speed turning, the gear transmission torque will traverse the whole engine output torque range. In this case, in order to ensure the dynamic stability and reliability of the front and rear wheels of the vehicle, it is necessary to study the static and dynamic characteristics of the transfer case beveloid gear.
Consequently, this paper aims to study the contact characteristics of the modification concave beveloid gear used in the AWD vehicle transfer case with misalignment. The defects of its theoretical tooth surface in practical application will be discussed, a modification tooth surface that can improve the gear meshing performance will be presented, and then the error sensitivity on static meshing characteristics is going to be studied by loaded tooth contact analysis (LTCA, an FEM method) method. Meanwhile, the dynamic response characteristics of the beveloid gear pair under various misalignment errors will be analyzed. Finally, the influence laws of each component of misalignment errors on dynamic meshing performance of the gear system would be obtained.
2. Beveloid Gear Pair Mesh Model and Mesh Characteristics
2.1. Work Pitch Cone and Meshing Equation
Intersecting beveloid gear pair is usually employed to transfer power with a small shaft angle which is less than 20°; especially, if the shaft angle is less than 10°, one member in the gear pair will be designed as a cylinder gear. When the beveloid gear pair is meshing, the pitch cylinder of the helical gear and the pitch cone of the beveloid gear is tangent at the point P, as shown in Figure 2. The origin of the global coordinate system locates in the intersection point of gear shafts, and the origin of the local coordinate system S1(x1, y1, z1) locates in the same position with , but of beveloid gear locates in the pitch cone vertex, and the positive direction of axis z is specified by the toe to the heel of gears. Here, the coordinate system is fixed to the global frame, and the coordinate systems and are fixed to the pinion and gear, respectively.

(a)

(b)
and n denote the unit normal vector of point P on the generatrix of the pitch cone and cylinder, respectively. The position vector of point P can be expressed as
The unit normal vector of point P can be calculated from
Putting formula (1) to (2), then
To obtain the tangential equation of mating pitch cone surfaces, transform the position vector and unit normal vector of point P into the global coordinate system Sg; the condition of tangency for mating pitch surfaces can be represented by
The formulas of coordinate system transformation are given bywhere is the matrix of coordinate transformation. According to the theory of engagement, the equation of meshing tooth surfaces is represented bywhere is the relative velocity vector of the pinion and gear contact point.
2.2. Transmission Error and Mesh Stiffness
The static transmission error (TE) is defined as the deviation of the ideal rotation angle and the actual rotation angle of the driven gear, but for the cylindrical gear pair, particularly, TE is always defined as a linear displacement in the line-of-action direction which is transformed from an angular displacementwhere rb2 is the base radius of the driven gear, Φ1 and Φ2 denote the ideal rotation angle of driving and driven gear, respectively, means the actual rotation angle of the driven gear, i is the transmission ratio, and δl and δ0 denote the loaded and unloaded transmission error in the line-of-action direction, respectively. When the gear tooth is meshing under load the instantaneous contact line can be divided into N microsegments; it is assumed that each microsegment normal vector is ni(nix, niy, niz), and the normal force is fi; as a consequence, the total contact force Fn and its three components of instantaneous contact line can be calculated byand then the actual mesh stiffness is expressed aswhere Ftotal is the total contact force of multiteeth meshing transiently, b is the face width, and M is the number of transient meshing tooth pairs.
By means of the loaded tooth contact analysis (LTCA) method, the distributed contact force and TE can be obtained; furthermore, the total contact force and mesh stiffness also can be calculated by equations (9)–(12).
2.3. Internal Dynamic Excitation of Meshing Gear Pair
In order to study the influence of time-varying mesh stiffness and transmission error on the dynamic mesh performance of the beveloid gear pair, a single-degree-of-freedom discrete rotor system of the beveloid gear is established as shown in Figure 3. The model consists of a pair of intersecting shaft gear pair, power input and output unit, and bearing support unit, drive and driven gear are modeled rigid bodies, and the meshing pair is modeled as a set of damping, stiffness, and error elements.

The equations of gear motion can be represented bywhere ki, ci, and ei, i = 1, 2, denote the mesh stiffness, mesh damping, and error of the i-th tooth pair, respectively, and cm and km indicate the combined mesh damping and mesh stiffness, respectively. It is assumed that the relative displacement of gear along line-of-action, namely, the dynamic transmission error (DTE) x, is expressed asand the excitation force F is expressed aswhere me denotes the equivalent mass and F0 indicates the static circumferential force:and the dynamic contact force Fd can be represented bytherefore, the DTE x can also be expressed as
The discrete steady-state response of the linear time-invariant model can be solved in the frequency domain, and the dynamic transmission error and dynamic mesh force can be obtained from equations (16) and (17), respectively.
2.4. Misalignment Conventions for Beveloid Gear Set
The misalignment of nonparallel shaft gear pair can reference hypoid gears’ definition. As shown in Figure 4, ΔΣ represents the shaft angle Σ deviation from the theoretical position; if ΔΣ is positive, this means the shaft angle Σ increases; ΔE indicates the offset distance E deviation from the theoretical position; if ΔE is positive, this means the offset distance increases; when the design offset is zero, the positive ΔE indicates the pinion deviation along the y direction. Gear mounting distance deviations (ΔXp and ) represent the mounting position error in the axial direction and specify the positive deviation toward the heel end (+z direction).

3. Loaded Tooth Contact Analysis and Tooth Surface Modification under Service Condition
3.1. Beveloid Gear Pair Geometry Parameters in a Transfer Case
A pair of beveloid gears are designed for an AWD vehicle transfer case, and the basic geometry parameters are shown in Table 1.
Take these gear pair parameters in Table 1 into equations (1)–(7); the point set of tooth surfaces which is shown in Figure 5(a) can be solved; furthermore, the 3D geometry model of the beveloid gear pair can be built as shown in Figure 5(b).

(a)

(b)
3.2. The Necessity of Tooth Surface Modification
Beveloid and helical gears’ theoretical tooth surfaces are both spiral involute surfaces, and their transient tooth contact pattern is a point without load, although the contact point will expand into a line when the teeth are loaded, but the contact area maintains narrow still, and the contact stress is high and prone to generate edge contact, which is very unfavorable to fatigue life.
Gear meshing performance has a great relationship with misalignment, which contains system deformation and assembly errors in working. The misalignment usually affects static and dynamic meshing performance, such as premature fatigue failure, scuffing, wear, vibration, and noise (Figure 6).

The beveloid gear pair in this example has a misalignment of △Σ = 0.313 mrad, △E = −0.06 mm, △Xp = 0.162 mm, and △Xw = −0.19 mm under 200 Nm input torque; the results of loaded tooth contact analysis (LTCA) of the theoretical tooth surface are shown in Figure 7(a); there is a serious edge contact on the tooth tip and root; edge contact can induce vibration and impact when the tooth enter or exit meshing and could lead gears to early failure such as pitting or scuffing if edge contact stress value is greater than the allowance value.

(a)

(b)
To improve the undesirable contact pattern, the edge contact, and vibration problem under misalignment, the modification tooth surface must be used to replace the theoretical one tooth surface. A concave parabolic end relief for beveloid gear was designed in the longitudinal direction relative to the theoretical tooth surface; that is, we took the theoretical tooth surface as the baseline; a parabolic end antirelief was used in the tooth toe and heel end; consequently, the instantaneous contact line will be lengthened, and the contact area will be extended. Tip and root relief, maybe plus involute barrelling, can resolve the edge contact problem, but the tip and root relief are difficult for the beveloid gear; therefore, the involute modification is implemented on the mating helical gear.
In the following, the work to optimize the contact pattern and reduce the contact stress and peak-to-peak transmission error (PPTE) by tooth surface microgeometry modification has been carried out. The optimized helical gear microgeometry modifications are shown in Table 2, and the microgeometry modifications of the beveloid gear are shown in Table 3. The modification topographies of the tooth surface are shown in Figure 6.
3.3. Results of Loaded Tooth Contact Analysis
3.3.1. The Improvement of Contact Pattern
After tooth surface micromodification, according to the data of Table 2 and Table 3, the contact stress of the tooth surface has been calculated by LTCA. Figure 7 shows the comparison of tooth surface contact patterns before and after optimization; there are indeed improvements after modification, the maximum contact stress is located in the center of the tooth surface, edge contact does not occur at the tooth tip and root anymore, the contact area is wider, and the load distribution is more reasonable.
3.3.2. The Improvement of Contact Stress
The contact stress curves of the tooth flank by quasistatic calculation are shown in Figure 8; we can find that the contact stress of the optimized tooth flank is generally smaller than the original flank; meanwhile, the maximum stress of the original flank is 1116.4 MPa, but after the optimized maximum stress declined to 736.2 MPa, the maximum value significantly reduced 34%, and the fluctuation degree of stress changed from drastic to stable.

(a)

(b)
3.3.3. The Improvement of Transmission Error
In the design objective of the automobile transfer case, the peak-to-peak transmission error (PPTE) should be less than 0.5 μm at least in the NVH testing conditions; in order to satisfy a low noise and vibration for the gear transmission system, it is well known that PPTE is a significant internal excitation. The TE curves calculated by quasistatic LTCA are shown in Figure 9; after tooth modification, the PPTE reduced from 0.51 μm to 0.39 μm, and the harmonics of TE after fast Fourier transform had also reduced.

(a)

(b)
4. Error Sensitivity Analysis of Beveloid Gear Meshing Characteristics under Misalignment
As we all know, the torque of a vehicle is constantly changing during driving, and the gear pair misalignment is changed with torque; therefore, whether a matched tooth surface optimization parameters can adapt to a variety of load conditions and whether it has a satisfactory meshing performance, practices have proved that this is a great challenge. Otherwise, if we can let the gear pair works in a well-supported environment so that the misalignment is in a reasonable range under the working load spectrum, then a set of optimized tooth surface parameters can be well adapted to different working conditions. Consequently, we should first study the sensitivity of gear meshing performance under different misalignment and then reduce the misalignment in the highly sensitive direction during structural design and assembly so that the gear system can acquire a good performance. Therefore, three types of misalignment of the modification beveloid gear pairs were studied separately, and their effects on the static and dynamic transmission performance were analyzed.
4.1. Error Sensitivity Analysis of Static Meshing Characteristics
4.1.1. Contact Patterns
The contact patterns under diverse shaft angle deviation ΔΣ are shown in Figure 10 (top). When ΔΣ is positive, the contact pattern tends to move to the toe end; otherwise, when ΔΣ is negative, the contact pattern will move to the heel end; if the absolute value of ΔΣ is large enough, the stress concentration will appear on the heel or toe edge of the tooth, and the contact area will be smaller as long as there exists an error of shaft angle. Consequently, the contact area is influenced significantly by the shaft angle deviation.

The contact patterns under various offset deviations ΔE are shown in Figure 10 (middle). The effect of ΔE on tooth contact pattern is similar to the ΔΣ; likewise, if the absolute value of ΔE is large enough, the stress concentration will appear on the heel or toe end.
The contact patterns under different mounting distance deviation ΔX are shown in Figure 10 (bottom). As a rule, cylindrical gears are insensitive to axial direction displacement; the cylindrical gear has neither a positive ΔXp nor a negative ΔXp; the contact pattern will hardly change, but the contact pattern of the beveloid gear has an obvious change if is not zero; when is positive, the contact pattern moves to the toe end slightly; in contrast, the contact pattern moves to the heel end slightly when is negative. Furthermore, we can use of the beveloid gear to adjust the backlash of the gear pair and to modify the transmission accuracy.
4.1.2. Transmission Error
The TE curves under various misalignments are shown in Figure 11. It can be seen that the ΔΣ has a few influences on the TE; the negative ΔΣ makes the average TE and PPTE decline, but the positive ΔΣ changes the average TE and PPTE indistinctively. ΔE can reduce the average TE, but it will increase the PPTE value. ΔXp has a hard influence on the TE but is similar to ΔE, which has a signal influence on TE and PPTE.

(a)

(b)

(c)

(d)
4.1.3. Time-Varying Mesh Stiffness
Figure 12 shows the time-varying mesh stiffness curve of the gear set under diverse misalignment; it can be seen that the shaft angle deviation ΔΣ has slight effects on the mesh stiffness; specifically, the negative ΔΣ strengthens the average mesh stiffness, and the positive ΔΣ changes the average mesh stiffness indistinctively; with the increase of ΔΣ, the mesh stiffness peak-to-peak value first increases and then decreases. The offset distance deviation ΔE will increase the mesh stiffness clearly; both positive and negative ΔE will strengthen mesh stiffness and increase peak-to-peak value. The wheel mounting distance deviation has a regular influence on mesh stiffness as same as ΔE, and pinion mounting distance deviation ΔXp has a great effect on the mesh stiffness.

(a)

(b)

(c)

(d)
4.2. Error Sensitivity Analysis of Gear Dynamic Characteristics
The misalignment also has a certain influence on the dynamic performance of the gear system, such as dynamic transmission error (DTE) and dynamic excitation force, which are two excitation sources of vibration and noise in the gear system.
4.2.1. Dynamic Transmission Error
The DTE in the frequency domain curves under various misalignments is shown in Figure 13. It can be seen that the negative ΔΣ makes the DTE decline, but the positive ΔΣ changes the DTE indistinctively. The plus or minus ΔE will increase the DTE value; ΔXp and have little influence on the DTE. As a result, the negative shaft angle deviation and absolute offset distance deviation have a great effect on the DTE.

4.2.2. Dynamic Excitation Force
The dynamic excitation force (DEF) in the frequency domain curves under diverse misalignment is shown in Figure 14. It can be seen that whether ΔΣ, ΔE, ΔXp, or has a regular influence on the DEF, the same as DTE above. As a consequence, the negative shaft angle deviation and absolute offset distance deviation have a great effect on the DEF.

4.2.3. Vibration Response Acceleration
The bearing vibration response acceleration (VRA) curves of the support structure response in the frequency domain under diverse misalignment are shown in Figure 15. It can be seen that if ΔΣ is unequal to zero, the VRA of support will decrease; when ΔΣ = 0.3 mrad, VRA has reduced about 3.9%, and the frequencies of peak points have not changed, and when ΔΣ = −0.3 mrad, VRA has reduced about 15.5%, but their frequencies of peak points have all increased. Whether ΔE is positive or negative, VRA will increase about 8%. Otherwise, ΔXp has a little influence on VRA, and with a range of −0.2 mm to +0.2 mm can make no more than 12% increment in VRA.

As a result, the significant influence of misalignment components on the vibration of the transmission system is negative ΔΣ, absolute ΔE, and ; the former tends to improve the dynamic performance, but the latter two are inclined to worsen the dynamic performance.
5. Conclusions and Discussion
In this paper, an intersecting shaft beveloid gear pair has been taken as the study object. The defects of its theoretical tooth surface in practical application were analyzed. On this base, a modification tooth surface that can improve the gear meshing performance was put forward, and then the sensitivity of the beveloid gear pair under various misalignment errors was analyzed. Finally, the influence law of each component of misalignment errors on the dynamic meshing performance of the gear system has been obtained. Through the above studies, some main conclusions are listed as follows:(1)The contact pattern of the theoretical tooth surface of the beveloid gear pair is a point, which will expand into line contact after loaded, but the length of the contact line is short, which often leads to large contact stresses and serious edge contact. Therefore, gear modification is recommended, widens the instantaneous contact line length, and eliminates the edge contact to improve the contact stress distribution and fatigue life.(2)For the intersecting shaft beveloid gear pair, the contact pattern is more sensitive to the shaft angle and offset error but less sensitive to the beveloid gear mounting distance error and almost not sensitive to the helical gear mounting distance error. The contact pattern will move to the toe end if the shaft angle and offset error are positive values and will move to the heel end if they are negative values. TE and PPTE are more sensitive to the offset error and beveloid gear distance error, more sensitive to the negative shaft angle error than the positive value, and almost not sensitive to the helical gear mounting distance error. Except the shaft angle error, the others error will increase the PPTE value. Mesh stiffness is more sensitive to the offset error and beveloid gear mounting distance error; the peak-peak amplitude of mesh stiffness will increase if there exist errors, whether the error is positive or negative. Besides, the shaft angle is not sensitive to meshing stiffness.(3)The influence of misalignment on the dynamic characteristics of the beveloid gear system is different from the static characteristics. For the dynamic exciting force, the offset error will increase the DEF, but the negative intersection angle error will decrease the DEF. For the vibration response, the offset error will increase VRA, but the shaft angle error will decrease VRA; this is an interesting phenomenon for beveloid gear transmission. In addition, the influence of beveloid gear mounting distance error on the dynamic performance of the gear system is significant. Therefore, the offset and beveloid gear mounting distance deviation have the greatest impact on the system dynamic performance.
Consequently, in the design process of the gear support structure of beveloid gear transmission, enough stiffness should be designed in the offset direction, and also the assembly error in the offset direction and beveloid gear mounting direction should be as small as possible; these are good for system NVH performance.
Data Availability
The data cannot be shared due to nondisclosure agreement with the industrial partner.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the National Key Research and Development Project of China (no. TC210H02Q, “Development and Industrialization of High Power and High Efficiency Electric Drivetrain, 2021”).