Abstract
As an important parameter to distinguish noncircular gear from cylindrical gear, eccentricity is very important for the meshing characteristics and transmission error of noncircular gear. In order to study the transmission characteristics of the elliptic gear, a pair of elliptic gear in the reversing device of a new type of drum pumping was taken as the research object. Based on the analysis of the transmission pressure angle and instantaneous contact ratio of the elliptic gear, the eccentricity error was introduced into the analysis model of transmission error. The influences of the eccentricity on the transmission pressure angle, instantaneous contact ratio, and transmission error were analyzed, and the analysis accuracy is verified by the finite element method. The results show that the eccentricity has a great influence on the transmission pressure angle, instantaneous contact ratio, and transmission error of the elliptic gear, and the eccentricity error has a significant influence on the transmission error. In order to ensure the normal meshing condition of the elliptic gear, the eccentricity should be less than 0.7071, and the maximum instantaneous contact ratio is 1.809. The research results can provide some guidance for the following noncircular spur gear transmission test and transmission error research.
1. Introduction
As a new type of gear, elliptic gear has the characteristics of strong bearing capacity, compact structure, and variable ratio transmission. It is mainly used in low-speed and large torque occasions, such as agriculture, textile, instruments, navigation ships, hydraulic pump, hydraulic motor, and flow instrument [1]. As is known to all that gear transmission transmits power by the mesh force along the mesh line, the dynamic excitation along the mesh line will then be produced in the transmission process, so that the displacement difference between the driving and driven gears in the direction of the meshing line will be generated. This is what is called transmission error, which can be divided into static and dynamic transmission error [2]. At present, there are mainly two kinds of transmission error research methods: analytical method and numerical method [3], while the numerical calculation methods are mainly divided into conventional finite element method and contact finite element method. The former is mainly represented by gear loaded contact analysis method (LTCA). The contact finite element method is to establish a nonlinear contact element in the meshing area and simulate the dynamic meshing process of the gear and simultaneously obtains the gear deformation, transmission error, and stress [4].
For the gear transmission system, the transmission error is the main factor affecting the accuracy of motion, so it is necessary to study it. For the transmission error of cylindrical gears, a lot of research results have been accumulated. References [5–8] analyzed the transmission error of the helical gear. Hotait and Kahraman [9] studied the relationship between the dynamic transmission error and the dynamic stress parameters of the spur gear pair. Park et al. [10] used the integrated empirical mode decomposition method to analyze the transmission error of the gear experimental device and apply it to the fault diagnosis of the gear. On the basis of studying the modeling of time-varying meshing stiffness and static transmission error of gear, the nonlinear frequency response characteristics of spur gear were analyzed by Yang et al. [11]. The contact model and transmission error of lightweight gear were studied in [12] by means of finite element analysis. The influence of assembly error on transmission error of gear was studied by Chun and Chen [13]. Park [14] analyzed the tooth surface friction and gear transmission error caused by sliding friction of spur gear under quasi-static condition. Based on the consideration of tooth profile modification, manufacturing error, and assembly error, Li et al. [15] proposed a calculation method of no-load transmission error based on measured discrete tooth surfaces. By combining with the principle of triangle intersection judgment, a method of gear tooth contact analysis is proposed to analyze the no-load transmission error, which provides a new method and idea for subsequent transmission error analysis of gear pair. Chin [16] analyzed the influence of average wear depth on gear transmission error. By simulating the meshing process of gear under dry friction condition, the change trend of transmission error was obtained and the efficiency of the proposed method was verified. Based on the establishment of load sharing ratio and quasi-static transmission error, Miguel [17] analyzed the influence on the quasi-static transmission error of symmetric long profile modifications on high contact ratio spur gears and obtained the dynamic matching relationship between the optimal tooth profile repair length and coincidence degree. References [18–22] described the calculation method of different types of gear load transmission error in detail. In [23–27], the analytical calculation model of transmission error was derived under the condition of considering eccentricity error and installation error, and the transmission error of gear was analyzed by finite element method, which provides a new method for the subsequent analysis of gear transmission error.
The above research results have certain guiding significance for the study of transmission errors of elliptic gears. At present, there are few literatures about the transmission characteristics and transmission errors of noncircular gears. Therefore, this paper takes a pair of elliptic gear pairs in the reversing device of a new type of drum pumping as the research object and introduces the eccentric error into the gear transmission error calculation model based on the analysis of the meshing characteristics of the elliptic gear. On this basis, the transmission error of the elliptic gear is studied. Figure 1 shows the elliptic gear reversing device model.

2. Theoretical Mathematical Models
2.1. Mathematical Model of Elliptic Gear
It can be seen from [28] that the elliptic gear has two focal points. In Figure 2, the equation of elliptic pitch curve with focal point O1 as polar coordinate point is

The tooth profile of elliptic gears can be divided into two parts, the point higher than the pitch curve and the point lower than the pitch curve, and there are different methods for solving the two-part tooth profile equation.
For points above pitch curve profile, the angles between the vector and the polar axis are (right profile angle) and (left profile angle), as shown in Figure 3.

The equation of the right tooth profile is
The equation of the left tooth profile is
For points on the tooth profile below the pitch curve, the angles between the vector and the polar axis are (right profile angle) and (left profile angle), as shown in Figure 4.

The equation of the right tooth profile is
The equation of the left tooth profile is
According to equations (2) to (5), the three-dimensional tooth surface equation of elliptic gear can be obtained, and the right tooth surface equation of elliptic gear is
The left tooth surface equation of elliptic gear iswhere refers to the direction of the tooth line and is equal to the width of the tooth.
2.2. Transmission Pressure Angle of Elliptic Gear
In Figure 2, gear 1 and gear 2 are, respectively, driving gear and driven gear, and the center distance of the elliptic gear is . The pitch curve equation of the driven gear is
The pressure angle of the circular tooth profile at the base circle is equal to the tooth angle of the tool rack profile (), but the transmission pressure angle of the elliptic gear changes due to the time-varying pitch curve radius. Figure 5 shows the transmission pressure angles of the left and right profiles of the elliptic gears.where + and − represent left and right tooth profiles, respectively.

(a)

(b)
According to literature [29], the transmission ratio of elliptic gear is
By simplifying equation (3), we can get
2.3. Instantaneous Contact Ratio of Elliptic Gear
The instantaneous contact ratio of noncircular gear is the same as that of cylindrical gear, which is the ratio of the effective meshing line length to the gear base circular pitch. Figure 6 shows two curves are tangent at point P at a certain moment. The curvature radii of pitch curves of driving and driven gears at point P are, respectively, and . It can be approximated that the cylindrical gear pair with and as the pitch radius is used to calculate the contact ratio instead of the noncircular gear [29]. The curvature radius of each point on the pitch curve of the elliptic gear is different, so the contact ratio in the meshing process is also variable.

The angle between the engagement line and the tangent of the pitch curve at point P is in Figure 6. When letting the intersection of the tooth top and the meshing line of gears 1 and 2 be A and B, respectively, and setting , we can get
Therefore, the contact ratio of elliptic gear is
2.4. Elliptic Gear Transmission Error Mathematical Model
Figure 7 shows the meshing model of elliptic gear pairs with the same parameters. The center of rotation is at the left focal point of the elliptic curve, and gear 1 and gear 2 are, respectively, driving gear and driven gear. In Figure 7, and are the centers of rotation of the gears with eccentric errors, and their trajectories are circles with and as the center of the circle and eccentric errors and as the radius, respectively. If and are rotated to the positions of and , and point is the node at the time of standard meshing. At this time, and are the instantaneous meshing curve radii of the elliptic gear. The connection line of the two gear rotation centers is , and the tangent point of the two instantaneous pitch curves is . The instantaneous meshing line of the two gears intersects the line of the theoretical center of revolution at point . The parameters of the elliptic gear are shown in Table 1.

The coordinate system is established with the rotation center of driving gear 1 as the origin of coordinate axis and as the y-axis. Then the coordinates of points , , and are, respectively,
It can be seen from Figure 8 that the instantaneous transmission ratio of the elliptic gear pair is

The theoretical transmission ratio is
The transmission ratio error is
The slopes of line segments and are, respectively,
Let ; the ordinate of point P can be obtained as follows:
At this time, the transmission ratio error is
For the pitch curve radius of the elliptic gear , equation (20) can be simplified to obtain
The ideal transmission ratio of a pair of elliptic gears is
The actual transmission ratio is
The transmission ratio error and the rotational angle error of driven gear are, respectively,
If the parameters in Table 1 are brought in, the transmission error of the elliptic gear pair can be obtained as follows:
3. Transmission Pressure Angle
According to equation (11), the variation of transmission pressure angle of elliptic gear under different eccentricity is obtained, as shown in Figure 9. When the eccentricity is 0, the transmission pressure angle is which is equivalent to the pressure angle of the spur gear. With the increasing of eccentricity, the transmission pressure angle also increases. The corresponding relationship between different eccentricity and the maximum value of transmission pressure angle is shown in Figure 10.


It is pointed out in [29] that when the pressure angle is large, the force arm between the gear teeth is small, and the force is greater when the same torque is transmitted. When is too large, self-locking may occur, making the gear unable to rotate. In order to avoid this phenomenon, the range of angle is between and ; at this time, the angle is less than . According to the above theory, the derivative of equation (11) is obtained:
Taking equation (26) into equation (11), we get
Therefore, in order to ensure the normal transmission of the elliptic gear pair without self-locking, its eccentricity should be less than 0.7071.
4. Instantaneous Contact Ratios
The variation of the contact ratio of elliptic gears under different eccentricity is shown in Figure 11. If , the contact ratio is always 1.7547. When , the instantaneous contact ratio shows periodic change. There are two peaks in the contact ratio of the elliptic gear when it rotates for one circle. The reason is that the contact ratio of the elliptic gear is related to the curvature radius of the pitch curve. The elliptic pitch curve takes its long axis as the symmetry axis and distributes symmetrically. In this process, the variation of the curvature radius shows two-week distribution. In Figure 11, with the increase of eccentricity, the maximum contact ratio of elliptic gear increases; and when , the contact ratio is 1.8077.

It is the basic condition of gear continuous engagement that the contact ratio is not less than 1. The larger the contact ratio is, the more teeth are engaged in engagement at the same time, and the single tooth load will be reduced, which is conducive to improving the service life of gear. According to the analysis of transmission pressure angle in Chapter 3, the maximum eccentricity is 0.7071, and the maximum contact ratio of elliptic gear is 1.809 at this time.
5. Calculation of Transmission Error
5.1. Transmission Ratio Error
Due to the skew and eccentricity of the gear base circle axis relative to the rotation axis, the error of gear transmission ratio is produced. When the eccentricity is fixed, the change of rotation angle determines the transmission ratio error. When the elliptic gear takes a focal point of its elliptic pitch curve as its rotation center, it is an eccentric elliptic gear. Considering the installation error, the transmission ratio error of the elliptic gear under different eccentricity conditions is shown in Figure 12. In order to reflect the periodicity of transmission ratio error, the horizontal axis is taken as , and the gear rotates for two cycles. With the increase of eccentricity, the fluctuation range of transmission ratio error increases gradually. The mean value and range of transmission ratio errors corresponding to different eccentricities are shown in Table 2. The mean value of transmission ratio error in Table 2 is basically about 0.015, and the range changes significantly with the increase of eccentricity. Transmission ratio error is mainly caused by manufacturing error, gear clearance, and other factors.

5.2. Transmission Error Calculation and Simulation Verification
The variation trend of transmission error of elliptic gear with eccentricity is shown in Figure 13. With the increase of eccentricity, the transmission error also shows an increasing trend. The previous analysis shows that the eccentricity should be less than 0.7071, and the maximum transmission error of the gear is 2.201 μm, which meets the design requirements. In order to verify the rationality of the analysis, the dynamic transmission error of the elliptic gear is obtained by using the finite element analysis, as shown in Figure 14. The black curve shows the no-load transmission error of the gear, and the red curve represents the dynamic transmission error of gear tooth obtained by LS-DYNA finite element software, which considered the load on the teeth.


In order to simulate the actual contact situation during the tooth engagement process, the following boundary conditions should be set: the inner ring of the rigid body shaft hole drives the gear body to rotate, the gear material is Solid-164 flexible body, and the inner hole of the shaft hole is Shell-163 rigid body. The driving and driven gears are limited to X, Y, Z three-direction moving degree of freedom and X, Y rotation degrees of freedom, and the rotation speed of the driving gear is 600 r/min. In the process of solving the tooth meshing model, the time step and the scale factor of the calculation time step are too large to interrupt the simulation, while the generation of negative volume is mostly caused by grid distortion, which is related to mesh quality and material and load conditions. Therefore, the appropriate time step should be taken to avoid the negative volume. The value of the debug time step scale factor TSSFAC is taken to be 0.5, and the time step DT2MS values can complete the analog tooth engagement. After gear meshing, the number of driving gear nodes is 205716, the number of units is 210211, the number of driven gear nodes is 181740, and the number of units is 186180.
The variation of the two curves is the same, but the dynamic transmission error curve has certain fluctuation. The reason is that the dynamic transmission error mainly considers the impact vibration, material damping, and other factors of the gear. In the process of analyzing gear tooth meshing, a load is applied and the penetration amount of gear tooth meshing is set. Also, it is assumed that the gear teeth have a certain rigidity, which makes the curve between the above two have a certain difference. The existence of dynamic backlash makes the curve of dynamic transmission error different from that of no-load transmission. The NLTE curve of theoretical analysis is relatively smooth, and there are violent fluctuations in the dynamic transmission error curve, which further shows that the calculation of dynamic transmission error by the finite element method truly reflects the meshing impact in the process of gear transmission.
5.3. Influence of Eccentricity Error on Gear Transmission Error
When manufacturing and installing the gear, it is inevitable to produce tooth shape error, pitch error, tooth direction error, and so forth, which will affect the transmission accuracy, transmission stability, and load distribution uniformity of the gear. The manufacturing error and installation error are reflected by the eccentric error of the gear. Figure 15 shows the effect of the eccentric error of the elliptic gear on the transmission error. With the increase of eccentricity error, the transmission error of elliptic gear is gradually increasing. Therefore, in order to ensure good transmission performance, the transmission error of gear can be reduced by improving the installation and manufacturing accuracy.

6. Conclusions
The eccentricity has a great influence on the meshing characteristics of the elliptic gear. By introducing the eccentricity error into the transmission error calculation model, the variation trend of transmission pressure angle, instantaneous contact ratio, and transmission error with eccentricity is obtained. The specific conclusions are as follows:(1)In order to ensure the normal meshing between the elliptic gears without self-locking, the eccentricity of the elliptic gears should be less than 0.7071.(2)The maximum instantaneous contact ratio of elliptic gear will increase with the increase of eccentricity. With the increase of contact ratio, the number of teeth participating in meshing at the same time will be increased, and the load of single tooth will be reduced, which is conducive to improving the service life of gears. Under the normal meshing condition, the maximum instantaneous contact ratio of the elliptic gear can reach 1.809.(3)The transmission ratio error and transmission error of the elliptic gear will increase with the increase of the eccentricity. In order to ensure the accuracy of the gear transmission, a small eccentricity can reduce the transmission error of the elliptic gear under the condition of satisfying the given motion law.(4)The eccentricity error has a significant effect on the transmission error of the elliptic gear. The transmission error will increase dramatically with the increase of the eccentricity error. The eccentricity error can be reduced by improving the installation and manufacturing accuracy of the elliptic gear, and then the transmission error of the elliptic gear can be reduced.
Symbols
r1, r2: | Radii of pitch curves of driving and driven gears |
A: | Long axis radius of ellipse |
c: | Half focal length of ellipse |
e: | Eccentricity of ellipse |
t1, t2: | The angle of the ellipse around the rotation center |
a: | Center distance of elliptic gear |
: | Profile angle of tool rack |
: | Transmission pressure angle |
n-n: | Tangent pitch curve at point P |
s: | Normal pitch curve at point P |
: | Angle between tangent and normal |
, : | Curvature radius of pitch curve of driving and driven gears at point P |
, : | Rotation speed of driving and driven gears |
: | Transmission ratio |
: | Contact ratio |
O1X1Y1: | Coordinate system built on the driving gear |
: | The center of rotation of gear teeth under the condition of eccentricity error |
P0: | Node in standard engagement |
R01, R02: | The instantaneous curve radius of elliptic gear |
e1, e2: | Eccentricity error of driving and driven gears |
: | The instantaneous curve radius of elliptic gear under the condition of eccentricity error |
, : | Initial phase angle of driving and driven gears |
, : | Ideal rotation angle of driving and driven gears |
x1, y1: | Coordinate of |
x2, y2: | Coordinate of P′ |
x3, y3: | Coordinate of |
, : | Actual rotation angles of driving and driven gears. |
Data Availability
The data used to support the findings of this study are included within the article, and the results of the article can be obtained by simply entering the various parameters into the formula derived in the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The research was supported by the National Natural Science Foundation of China (no. 51765032).