Abstract

A design method for the geometric shape and modification of asymmetric spur gear was proposed, in which the geometric shape and modification of the gear can be obtained directly according to the rack-cutter profile. In the geometric design process of the gear, a rack-cutter with different pressure angles and fillet radius in the driving side and coast side was selected, and the generated asymmetric spur gear profiles also had different pressure angles and fillets accordingly. In the modification design of the gear, the pressure angle modification of rack-cutter was conducted firstly and then the corresponding modified involute gear profile was obtained. The geometric model of spur gears was developed using computer-aided design, and the meshing process was analyzed using finite element simulation method. Furthermore, the transmission error and load sharing ratio of unmodified and modified asymmetric spur gears were investigated. Research results showed that the proposed gear design method was feasible and desired spur gear can be obtained through one time rapid machining by the method. Asymmetric spur gear with better transmission characteristic can be obtained via involute modification.

1. Introduction

Gear drives are one of the most widely used modes in power transmission system with the advantages of large power, high efficiency, and long service life. The involute tooth profile has various advantages, such as good transmission performance, convenient installation, and low cost [1, 2]. Compared with traditional involute spur gears, asymmetric involute gears have larger carrying capacity, smaller volume, lighter weight, and better dynamic characteristics [3]. So the asymmetric spur gear has significant potential applications in the fields of heavy load and high-speed gear transmission system. In the gear transmission process, it is not possible to have the mating gear enter contact in pure involute position because there would be sudden interference corresponding to the elastic deflection and the corner of the tooth tip would gouge into the mating surface [1]. Real meshing tooth profile of asymmetric spur gears will deviate from the ideal location under load and will cause transmission error and the change of load distribution. The transmission error changes over time, which causes displacement excitation in the gear dynamic model and has great influence on the dynamics of meshing gears [46]. Therefore, it is necessary to improve the meshing performance of asymmetric spur gears by the reasonable design and modification.

Until now, many researches have been carried out on the design and modification of symmetric gear. Litvin et al. [7] proposed a modified involute helical gear that matches with a double-crowned pinion and simulated the meshing process between profile-crowned and double-crowned helicoids. Faggioni et al. [8] presented a global optimization method that focused on reducing gear vibration by means of profile modifications. Barbieri et al. [9] used genetic algorithms to research the optimum profile. Zhang et al. [10] performed computer-aided design based on the contact conditions of conjugated gears and the tooth end relief with helix by finite element analysis (FEA). Gonzalez-Perez et al. [11] applied analytical approach in the stress analysis for gear drive with localized bearing contact based on the Hertz theory. Huang and Su [12] presented various approaches about element constructions and dynamic analyses for involute spur/helical gears where meshed elements are constructed directly without CAD geometric model. The gear profile equations were described via C programming language. The effects of modification coefficients and helical angles on the transient meshing performance of the gears were investigated using the method of explicit dynamic FEA from the perspective of energy by Hu et al. [13]. A qualitative analysis about the effects of flank deviation on load distribution for helical gear was performed using FEM by Wei et al. [14]. Wu et al. [15] proposed a tooth profile modification approach for the helical pairs, and the static contact process was analyzed using FEA and experimental tests. Li [16] applied 3D FEM to conduct surface contact stress and root bending stress calculations for a pair of spur gears considering machining errors, assembly errors, and tooth modifications. Simon [17] researched on the influence of tooth modifications induced by machine tool setting and head-cutter profile variations on tooth contact characteristics in face-hobbed spiral bevel gears. Chen and Shao [18] proposed a general analytical meshing stiffness model to establish the relationship between the gear tooth errors and the total mesh stiffness, wherein load sharing among different tooth pairs during meshing and static transmission errors under load were considered. Sankar and Nataraj [19] proposed a modification design approach for increasing the tooth strength in spur gear. Wang et al. [20] proposed a method to simulate residual stress field of planetary gear. Li and Kahraman [21] proposed a tridynamics model for spur gear pairs. The proposed model couples a mixed elastohydrodynamic lubrication model of a spur gear pair with a transverse torsional dynamic model. Kolivand and Kahraman [22] defined new general concept based on the construction of a surface of roll angle to simplify the task of locating instantaneous contact lines of any type of gearing. With the surface of roll angle defined, instantaneous contact lines are determined by a novel approach based on an analogy to parallel-axes gears. Eritenel and Parker [23] provided an analytical solution for the nonlinear vibration of gear pairs that exhibit partial and total contact loss.

Moreover, more and more attention has been paid to asymmetric gear recently. Litvin et al. [3] proposed the modified geometric model of an asymmetric spur gear drive which integrated an involute gear and a double crowned pinion and simulated the meshing and contact process of the spur gears with proposed geometric parameters using computer-aided design method. Alipiev [24] proposed a new geometric design method for the symmetric and asymmetric involute meshing in which the gear load sharing ratio is equal to its expected value. Pedersen [25] observed the use of two new standard cutting tools for the asymmetric gears since the shape of the gear profile influences its bending strength.

In the gear design process, the traditional method only gives basic parameters of gear structure and it is difficult to obtain precise model for the actual generated gear. In this paper, the equations of gear tooth profile were deduced based on a given rack-cutter, and the accurate geometric model was established in the design process. In the traditional machining process of modified gear, a standard gear without modification is generated at first, and the modification is conducted subsequently. In order to enhance the manufacturing efficiency, the pressure angle modification is carried out on the rack-cutter, so the desired modified spur gear can be obtained through one time rapid machining.

Based on the proposed gear design and modification method, a FEM of asymmetric involute gear was established where computer-aided designs were adopted to establish the geometry model of the spur gears; finally the meshing process was simulated.

The present paper mainly consists of the following three aspects.(1)The equations of a new rack-cutter profile for asymmetric gear were established. The pressure angle modification in rack-cutter design was proposed and the driving side of rack-cutter profile equations was deduced. The given rack-cutter profile has different pressure angles and fillet radius in the driving side and coast side. The generated asymmetric spur gear profiles have different pressure angles and fillets accordingly (Section 2).(2)According to conjugate theory, the equation of asymmetric spur gear was derived by the forming rack-cutter tool. Computer-aided design of spur gears with proposed geometric parameters and the simulation of their meshing were conducted using the proposed methods (Section 3).(3)The transmission error and the load sharing ratio of generated gears were studied. The transmission principles of asymmetric spur gear were discussed (Section 4).

2. Geometric Characteristics of the Rack-Cutter Profile

In order to obtain the precise gear model, the geometric characteristic of the corresponding rack-cutter profile is required. The given rack-cutter has different pressure angles and fillet radius in the driving side and coast side.

2.1. The Rack-Cutter Normal Profile

The typical rack-cutter consists of straight lines and tip circular fillet, which can generate involutes and fillets of typical involute tooth, respectively. As depicted in Figure 1, the tool coordinate system is rigidly fixed to the gear rack-cutter, and the coordinate axis is the reference line of rack-cutter. The basic normal profile of rack-cutter consists of four segments: driving side flank straight line , coast side flank straight line , driving side flank tip fillet curve , and coast side flank tip fillet curve . The gear tool possesses different gradients on either side of the rack-cutter straight tooth profile and has the same addendum coefficient and tip clearance coefficient. The profile of the rack-cutter in can be described with the following equations.(i)Coast side flank straight line is where is an arbitrary point on line , is a variable parameter about the length of straight line , is the pressure angle representing the slope of line , is the radius of fillet curve , is the addendum coefficient factor, and is the normal module.(ii)Driving side flank straight line is where point is an arbitrary point on line , is the variable parameter about the length of straight line , is the pressure angle representing the slope of line , and is the radius of fillet curve .(iii)Coast side flank cutter fillet curve is where is the angle of and is an arbitrary point on the fillet curve .(iv)Driving side flank rack-cutter fillet curve is where is the angle of and is an arbitrary point on the fillet curve .

The radii of fillet curve and meet where is the tip clearance coefficient.

The tangential vector to the generated shape in the is

The unit normal vector in is

The unit normal vectors of each profile in can be derived as where the subscripts , , , and represent the corresponding line and curves.

2.2. The Equations of Pressure Angle Modification

The conventional design of asymmetric spur gear is based on the application of involute profiles [1]. If the lines in the tool tip have different slopes, the generated gear will be modified correspondingly. Thus the driving side of rack-cutter will have two pressure angles, and the generated gear will have two segments of involute profile on the driving side. It is easy and feasible to produce such a rack-cutter.

The pressure angle modification can be represented by two parameters in the tool coordinate: the modified angle and the modified length , as shown in Figure 2. The parameter is called modified coefficient. The modified new gear profiles are shown in Figure 3. is the maximum modified quantity; is the maximum modified length. is the distance between the modified point and the gear tip, which should be within the range of .

The equations of modified tool profiles are listed as follows.(i)Driving side flank straight line is where is the modification coefficient; is the pressure angle after modification.The unit normal vector of line meets (ii)Driving side flank straight line is The unit normal vector of line meets

3. The Gear Profiles Equations

The gear profile equation is the foundation of the parameter modeling. Due to the application of coordinate transformation, we can deduce the generated gear profiles from the equation of rack-cutter profile based on the kinematic relation and conjugate theory. After the pressure angle modification of rack-cutter, the generated gear will be modified with involute curve.

3.1. The Coordinate Transformation

Coordinate systems , , and were rigidly fixed to the rack-cutter, the spur gear to be generated, and the frame (see Figure 4) [7]. A point in coordinate system was represented by position vector . A vector in coordinate system was signed as . is the rotate angle of gear blank. The coordinate transformation from to was based on the matrix equation The transformation matrix can be determined by In this way, we can obtain the equation of gear profile from that of rack-cutter profile as follows:

3.2. Meshing Equations

Suppose that the rack moves at a velocity of and the gear rotates at an angular velocity of (Figure 4). The instantaneous center of rotation is located on line that is perpendicular to . The location of satisfies the vector equation The relative motion of centrode is pure rolling around and the displacements of rack and the rotate angle of gear blank can be correlated by where is the pitch circle of gear.

The vector radius of a moving point in coordinates is The sliding velocity at point is where is the velocity of point in coordinates given by where From (8) through (20), we can derive

Any meshing point on the gear tooth profile needs to meet the following conditions: the relative motion velocity vector of the point should be located in a common tangent of the two tooth profiles along a direction; namely, where is the unit normal vector of meshing point (equation (8)) and is the relative velocity vector.

According to the normal vector of different line segments and the relative velocity in (23), meshing equation can be obtained by

3.3. Conjugate Gear Profile Equations

The tooth shape of the generated gear without modification is represented by (15) and (24). Each segment of gear profile curves is listed as follows.(i)Coast side flank involute gear profile is (ii)Driving side flank involute gear profile is (iii)Coast side flank gear fillet curve is (iv)Driving side flank gear fillet curve is

We can deduce the general equation of asymmetric involute gears profile with double pressure angles in from (1)–(4) and (25)–(28) as follows.(i)Coast side flank involute gear profile is (ii)Driving side flank involute gear profile is (iii)Coast side flank gear fillet curve is (iv)Driving side flank gear fillet curve is

When the modification is taken into consideration, the modified profile of driving side can be derived as follows.

In Figure 2, the meshing process can be described by (33), according to (10) and (22)-(23): From (10) and (33), we can deduce the gear profile equation of profile : where .

Similarly, the deduced meshing equation of profile is The corresponding gear profiles satisfy the following equations: where .

3.4. Simulation of Meshing and Contact

The gear profile with and without modification was generated according to the derived equations in Section 3.3, and the corresponding finite element models were also established for subsequent numerical simulation.

The solid models are constructed with the following design parameters. The number of the teeth , normal module , face width , normal pressure angle of coast side flank , normal pressure angle of driving side flank , addendum coefficient , and tip clearance coefficient . The solid model was discretized with 8-node hexahedron isoparametric elements. The material was steel with the properties of Young’s modulus and Poisson’s ratio . The motor power , the rated speed of driving gear was , and the applied torque Nm according to the equation .

The considerations for the finite element models are as follows.(1)The 3D models are adopted due to the differences of tooth surface contact stress along the tooth width direction. Since the overlap ratio of gear pair is 1.463, it is suitable to have the solid model with five gear teeth and inner ring. The simplification is sufficient and highly efficient for the simulation. The two-tooth geometric model is shown in Figure 5(a).(2)Gear meshing process is a dynamic problem that the meshing position changes with the change of time. So dynamic meshing simulation can be discrete into a series of quasistatic processes. After establishing the corresponding quasistatic mathematical model, using the existing computer technology fulfills the iterative calculation repeatedly to get suitable results. Finite element models of the gear drive can be automatically obtained for any meshing position. Through rotating the gear drive, we can achieve the simulation in the whole meshing cycle.(3)The gear solid model was discretized with 8-node hexahedron isoparametric elements, as shown in Figure 5(b).(4)The tooth contact is performed by contact elements. A contact pair consists of a target surface and a contact surface. One contact pair is shown in Figure 5(c).(5)Boundary conditions were imposed automatically in each cylindrical coordinate system. The whole five teeth mating gear finite element models are shown in Figure 5(d).(6)Nodes of the inner rim of driven gear were fixed by restricting all degrees of freedom; nodes of the inner rim of driving gear are fixed by restricting the radial and axial degree of freedom and releasing the rotational degree of freedom; the transmission torque was exerted to each node in the form of circumferential force. The relationship between the number of nodes in the driving gear rim and the radius of rim meets the equation

Assumption of load distribution in the contact area is not required since the contact algorithm of the general computer program is used to get the contact area and stresses by application of the torque to the driving gear while the driven gear is considered at rest. Both finite element models of driving and driven gears are rotated to simulate whole meshing cycle, and nonlinear contact analysis was accomplished by the contact algorithm in the finite element program.

Through the numerical simulation, we can predict the characteristic of the gear drive comprehensively and can reveal the principle of gear mating in the design process.

4. Calculation Results and Discussions

The purposes of modification are to improve meshing performance and to realize load distribution balance. The two aspects can be evaluated by transmission errors and load sharing ratio.

The parameter (Figure 3) was defined as the modified ratio,

The parameter was defined as the modified angle meeting The pitch point where the gears began to mate with each other was defined as the valid starting point. The driven gear rotated counter-clockwise was defined as the positive direction, and the clockwise direction was chosen as the negative direction.

The theoretic angle at the start of active profile (SAP) is , and theoretic angle at the end of active profile (EAP) is according to the formula [26]. Therefore the actual mating angle of one gear tooth should be within the range of to . The theoretic angles of the lowest point of single tooth contact (LPSTC) and the highest point for single tooth contact (HPSTC) are and calculated by formula, respectively [2, 26]. The investigators chose the rotation angles of symmetric spur gear pairs ranging from to with an interval of .

4.1. Influence of the Transmission Errors

Transmission error (TE) is one of the most important factors causing gear noise and vibration. TE refers to the difference between the displacements of gears in reality and those in theory. TE is calculated as the relative displacement at the pitch point in the line of action between the two engaged gears, using small angle approximations for any relative rotations as shown in where and are the displacements of engaged gears at the pitch point in the line of action. The resulting transmission error is specified as a distance. The displacement of each gear can be obtained by the simulation results.

4.1.1. The Transmission Errors under Different Relief Lengths

When the torque () was applied, the modified angle , and the modified ratio values were 1/3, 1/2, 2/3, and 1. The TE curves of driven gear in a meshing cycle under different relief lengths are shown in Figure 6. TE peak amplitudes and percentage increments of double tooth meshing area compared with single tooth meshing area for the unmodified and modified gears under different relief lengths are given in Table 1. Through analysis based on Figure 6, the following results can be drawn.(1)Mating zones and are the single tooth meshing areas, and mating zone is double tooth meshing area. In the unmodified model, there exist convex line in zone and concave lines in zones and . The gear has lower mating stiffness in zone ; therefore the amplitude of TE in zone is greater than that in zones and . The lowest point in zones and is called the negative TE peak point. The TE amplitude in zones and is 39.782%, lower than in zone (Table 1).(2)Compared with the unmodified model, the relief length is a critical factor for determining TE curve shape. When , zone has wider space. Zones and have same concave shape as the unmodified model but have a sharp point at the lowest position. When , zone also has wider range than the unmodified model. The curves in zones and are rough and uneven. When , there is obvious difference between single tooth meshing and double tooth meshing areas. The curves in the zones and also change from concave to convex shape, and its range decreases with increasing zone . When , the relief starts from the pitch point, so zone changes from a line to a point. The curve in zones and changes from concave to convex shape dramatically.(3)The relief length has great effect on the TE amplitude, especially the double tooth meshing areas. When , the curve in zones and has sharper negative peak but even amplitude compared with the unmodified model. The negative peak amplitude in zones and is 6.57885, 40.0133% lower than the peak amplitude in zone . When , the negative peak amplitude in zones and is 11.64246 , 16.07% higher than the peak amplitude in zone . When , the peak amplitude in zones and is 15.71421, which is 55.98% higher than the peak amplitude in zone . When , the negative peak amplitude in zones and is 24.89898, 142.6177% higher than the peak amplitude in zone (Table 1).(4)The TE peak amplitude of modified model in zone changes little when compared with the unmodified one with the increasing modified factor . The TE peak amplitude of the modified model in zones and increases as the relief factor increases. It was found that the relief factor is relatively optimal to make the curve of TE level and smooth. When or , the relief length is sufficient, and when , the relief length is excessive.

4.1.2. The Transmission Errors under Different Relief Pressure Angles

When the applied torque was  N·m, the modified angle values were °, °, or °, and the modified ratio values were 1/3, 1/2, 2/3, or 1. The TE curves of driven gear in a meshing cycle under different relief pressure angles are shown in Figure 7. TE peak amplitudes of modified gears under different relief lengths are given in Table 2. Through the analysis based on Figure 7, the following results can be obtained.(1)The relief pressure angle has great effect on the variation range of the TE when the modified ratio is a constant. As the increases, the value of TE increases.(2)When , the maximum value of TE in single tooth meshing area increases with increasing . The negative peak amplitudes in double tooth meshing area change very little (Table 2). The width of double tooth meshing areas becomes narrower with increasing . When , the shape of curve in double tooth meshing areas changes from concave with two peaks into convex with one peak. The peak amplitudes increase with increasing . When or , the TE values increase with increasing .(3)The modification length is the main factor that determines the shape of the transmission error curve. Appropriate modification length can make the gear transmission error curve in the process of alternating single and double teeth smooth, thereby reducing the meshing dynamics. Modification angle is the main factor influencing the peak transmission error. Suitable pressure angle can make the transmission error fluctuation decrease. In present study, it was found that the TE curve would be level and smooth when for the assumed torque.

4.2. Influence of the Contact Force
4.2.1. The Load Sharing Ratio under Different Relief Lengths

Due to the material having flexibility, the deformed teeth can rotate, bend, and shear, and the tooth flanks have been flattened at the contact points. In this situation, the load varies outside the normal contact path of asymmetric spur gears. The curves of tooth load sharing ratio can clearly show the complex variation from the tooth running in, to running out of contact. The circumference force and radial force of meshing gear satisfy the following equations, respectively:

The theoretic total meshing force is given by where . Then the total meshing force can be determined to be by substituting the involved parameters. A single gear meshing force can be obtained by simulation results. The load sharing ratio is given by Taking one gear tooth for consideration in the whole meshing cycle, distribution of the contact force under different relief lengths is given in Figure 8. Through analysis based on Figure 8, the following results can be obtained.(1)The load sharing ratio is 100% in single tooth meshing area and 35%–46% in double tooth meshing areas of unmodified model. There is a sharp drop at the LPSTC and the HPSTC.(2)The load sharing ratio in double tooth meshing areas is 0–100%. The slopes of the line increase and the space becomes narrower with the increasing relief length in double tooth meshing areas.(3)When the rotation angles are within the range of 15°–9.6° and −9.5°–15°, the change amplitude of load sharing ratio decreases as increases. When the rotation angle is in the range of 9.6°–4° and from −4° to 9.6°, the change amplitude of load sharing ratio increases as increases.

4.2.2. The Contact Force under Different Relief Pressure Angles

When the applied torque was , the modified angle values were , , or , and the modified ratio values were 1/3, 1/2, 2/3, or 1. Through analysis based on Figure 9, the following results can be obtained.(1)When the rotation angles are in the range of 15°–9.6° and from −9.5° to 15°, the amplitude of load sharing ratio declines as increases. When the rotation angles are in the range of 9.6°–4° and −4°–9.6°, the change amplitude of load sharing ratio increases as increases.(2)The load sharing ratio is 100% in single tooth meshing area. The slope of the line increases and the space range becomes narrower with the increasing relief length in double tooth meshing areas.(3)Length of modification and modification angle determines the shape of the curve of load distribution and proportion. In the analysis of load distribution, is reasonable.

5. Conclusions

Based on the research performed, the following conclusions can be drawn.(1)A design method for the geometric profile of asymmetric spur gear was proposed, in which the geometric profile of the gear can be obtained directly according to the rack-cutter profile in the design process. The geometric equations of rack-cutter were conducted, and the corresponding equations of generated gear were deduced based on the conjugated theory.(2)A design method for the modification of asymmetric spur gear was proposed, in which the modification of the gear can be obtained via the pressure angle modification of the rack-cutter. The geometric equations of the modified driving side of rack-cutter were conducted, and the corresponding equations of the generated gear were deduced based on the conjugated theory.(3)The designed and modified gear models can be adopted to carry out the subsequent mating analysis. The research results indicated that the involute modification of the generated gear contributed to better transmission performance.

Nomenclature

:Metric module, mm
:Tooth number
:Pressure angle for asymmetric spur gear drive: for coast side profile, for driving side profile, for modified driving side profile
:Radius of rack-cutter fillet: for coast side profile, for driving side profile
:Addendum coefficient
:Tip clearance coefficient
:End point of coast side straight line
:Start point of coast side straight line
:End point of driving side straight line
:Start point of driving side straight line
:Bottom point of rack-cutter fillet curve
:Center point of rack-cutter fillet: for coast side profile, for driving side profile
:Arbitrary point in straight line , ,
:Arbitrary point in rack-cutter fillet curve ,
:Length of straight line , ,
:Angle of ,
:Modified coefficient
:Modified ratio
:Modified angle
:Coordinate systems: for rack-cutter, for conjugate gear, for fixed frame
:Radius of pitch cylinder of gear.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the Natural Science Foundation of China (no. 51105287) and the Key Research and Development Project of New Product and New Technology of Hubei Province (no. 2012BAA08003) for the support given to this research.