# Analysis of Contact Part of Error Tooth Surface and Dynamic Performance Prediction for Involute Gear

Authors: Xinrong Liu; Zhonghou Wang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Aiming at the problem of constructing digital model of involute gear with error, the method of linear interpolation combined with area weight interpolation is proposed. Based on the non-feature discrete data block technique, the true tooth surface discrete data obtained by the coordinate measuring instrument is divided into blocks, and then the interpolation method is used to interpolate the nonmeasurement area to construct the real tooth surface with errors. The contact part and dynamic performance of the gear are predicted by using the constructed error tooth surface. The contact error of the tooth surface and the transmission error of the gear are verified by the test, and the reliability of the judgment result is judged by measuring the vibration in the direction of the gear meshing line. Compared with the example, this method not only reduces the computational complexity of the interpolation algorithm, but also improves the accuracy of the tooth interpolation data points and the smoothness of the error tooth surface.


---

## Body

## 1. Introduction

It is well known that the contact situation of the gears in the gearbox directly affects the gearbox’s carrying capacity, vibration, and noise. The gear transmission error, the distribution area, and shape of tooth contact part have become an important index to measure the meshing performance of gear [1]. Due to gear manufacturing and installation errors, bearing clearance, load bearing parts’ deformation, and so on, these factors will cause the gears to deviate from the ideal contact state. In order to accurately predict the dynamic performance of involute gears, the construction of error tooth surface is very important.The reconstruction of complex surfaces is one of the most important research techniques in digital construction of tooth surface, and scholars at home and abroad have also done a lot of research about it. Lin et al. [2] assessed the large gear tooth profile deviation by using the NURBS (Nonuniform Rational B-Splines) surface. Zhang et al. [3] obtain digital tooth surface by using double three NURBS surfaces. And the tooth contact analysis of the digital tooth surface was carried according to the spatial meshing theory. Li et al. [4] construct the digital real model of hypoid gear aiming at the wear of tooth surfaces by using non-feature block interpolation technique. Sun et al. [5] propose an algorithm for constructing G1 continuous cubic triangular Bezier Surfaces from triangular mesh surfaces. In these studies, the detection of contact parts is mainly based on the traditional roll test and manual identification method. However, there are few studies about the contact surface of the tooth surface by the gear bearing contact analysis, and also there are studies about the prediction of meshing performance of the error tooth surface by contact part.When the coordinate measurement is performed on the tooth surface, these methods are mainly applied to the smooth tooth surface because of the limited measurement data points of tooth surface. As for irregular tooth surfaces such as error tooth surfaces, the reason why the structural tooth surface can not reflect the true tooth surface is that most of the measurement data is separated from the structural tooth surface. In the above discrete data preprocessing technology, data segmentation is the bottleneck problem of constructing digital real tooth surface.Based on the original research, an algorithm combining linear interpolation method with area weight interpolation method is proposed for constructing error tooth surface [6–10]. Besides, the digital construction of error tooth surface of involute gears is realized on the basis of NURBS surface. In this paper, the traditional smoothing algorithm is improved to ensure accuracy and increase the smoothness of digital tooth surface. The feasibility of the method proposed in this paper is verified by comparing the examples.In contrast to constructing the digital model of single tooth surface, NURBS surface is used to construct the error tooth surface model of modified involute gears in this paper. Combined with the principle of involute gear meshing, the solution of the contact part of involute gear is obtained.


## 2. Synthesis and Reconstruction of Error Tooth Surface

The sum of the deviation values of the normal direction at the corresponding contact point on the tooth surface is called the synthesis error. In this paper, a grinding method with high grinding precision is used to grind the gear. This method is through the grinding wheel to grind the gear, and then get the relevant tooth shape. When the standard gear is ideally engaged, the contact line will cover the entire tooth surface, which is called the geometric contact line. This line also serves as the reference for the contact position. The tooth surface error is defined as the size of the angle in the normal direction when the actual tooth surface corresponds to the reference plane.


### 2.1. Synthesis of Error Tooth Surface

In this paper, a base section of the gear is divided into 20 equal parts for gear contact analysis. The ideal contact part for the different size of gears is shown in Figure1.Figure 1
The geometrical contact line of the tooth surface during the meshing process.i: the corresponding contact line number during the engagement process; j: split points on the contact line.
(a)
The ideal contact part for spur gears
(b)
The ideal contact part for helical gearsThe interval between the contact lines is one twentieth of the rotation angle of a base section. The contact lines on the tooth surface of the large gear are numbered withi=1,2,…,N (N is the total number of contact lines) in the order of start contact to the end of the contact. The length of a contact line is divided into 19 equal parts, and the corresponding division points are numbered with j=1,2,…,20 in the order of the tooth root to the tooth top. The tooth surfaces meshing with the large gears are divided by the same method. Therefore, the reference plane position on the different size of gears can be expressed by the combination of the serial number of the contact line i and the division number j of the contact line. The shape deviation of the gear can be expressed as data 20×n.The sum of the deviation values of the normal direction at the corresponding contact point on the tooth surface is called the synthesis error. The curve of the 20 synthetic error values on a geometric contact line is called the synthetic error curve. In this paper, the synthetic error surface is expressed in a coordinate system with the split numberj of a geometric contact line being the X-axis and the contact line number i is the plane of the Y axis, which is shown in Figure 2. The direction of the synthetic error coordinate system e is positive and the highest point is zero. The minimum of the synthetic error on the geometric contact line is defined as the ridge. The projection of the ridge to the synthetic error coordinate system Y-e surface is called the ridge shape.Figure 2
The creation of synthetic error surfaces.


### 2.2. Tooth Surface Model Based on NURBS

NURBS surface is the only representation standard of product data model defined in STEP standard (Product data exchange standard), and it has been applied in real tooth surface simulation. However, in the process of real tooth surface reconstruction by using NURBS surface, the precision of the reconstructed tooth surface in the nongrid node is difficult to control due to the limited number of measurement data points. In order to ensure the smoothness of reconstructed tooth surface, this paper proposes a linear interpolation method combined with area weight interpolation method. This method is used to estimate the tooth surface data points which are outside the measurement area to realize the construction of the error tooth surface and to ensure the accuracy of the reconstructed tooth surface to the utmost extent.NURBS surfaces are tensor product forms based on NURBS curves, and the double three NURBS surface can be expressed as(1)Pu,v=∑i=1m∑j=1nBi,kuBj,kvωijVij∑i=1m∑j=1nBi,kuBj,kvωiju,v∈0,1,where m is the number of control vertices for u, n is the number of control vertices for v, and Vij is the control vertex of the surface. ωij is the weight factor of Vij. Bi,k is the kth spline basis function along u. Bj,k is the kth spline basis function along v. Among them, the basis function recursive formula is(2)Bi,0u=1ui≤u≤ui+10otherBi,ku=u-uiui+k-uiBi,k-1u+ui+k+1-uui+k+1-ui+1Bi+1,k-1uk≥1.


### 2.3. Delaunay Triangulation of Discrete Data

In this paper, the scanning method is used to obtain the tooth surface discrete point data. In order to describe the real tooth surface situation more intuitively, this paper deals with the discrete data points by means of Delaunay triangulation principle because the topological relations between the discrete data points are complex. The traditional Delaunay triangulation method often projects the three-dimensional discrete data points onto the two-dimensional plane and then triangles the projection points. However, the triangulation of this method cannot truly reflect its spatial angle, which will affect the quality of the section. On the basis of these methods, this paper presents a method of direct triangulation of three-dimensional discrete data points.Assuming that the scan line isli, k=1,2,…,n, the number of measurement points on the scan line is t, where two adjacent scanning lines are lk,lk+1,1≤k<k+1≤n. The measurement points on each of the two scanning lines are lk(j), j=1,2,…,t;lk+1l,l=1,2,…,t.As shown in Figure3, the starting points for connecting two scan lines are lk(1) and lk+1(1); the end point are lk(t) and lk+1(t). And find another scan line that is closest to the measurement point in the measurement points on lk and lk+1 and connect, where the area between the scanning lines lk and lk+1 is divided into r(r≥t-1) spatial regions. Finally, the triangular and spatial quadrilateral regions are retrieved.Figure 3
Preliminary Delaunay triangulation of adjacent tooth surface scan lines.For the spatial quadrilateral region, the minimum internal angle maximum criterion is used to complete the Delaunay triangulation of the discrete data points in this paper, respectively, as shown in Figure4.Figure 4
Polygon area triangulation.


### 2.4. Linear Interpolation Method and Area Weight Interpolation Method

In order to describe the error of the tooth surface quickly and accurately, it is necessary to block the discrete data points of the tooth surface in the construction of the digital real tooth surface. And use the interpolation method to accurately describe the error surface to reduce the amount of tooth interpolation calculation. The traditional data block method is often based on the geometric characteristics of the part, which does not apply to free surfaces that describe the error tooth surface. In this paper, an interpolation method based on linear interpolation combined with area weight interpolation is proposed to accurately identify tooth surface errors.The error tooth surface shape is calculated by using the area weight method according to the following interpolation method. In Figure5, assume that A,B,C,D are the measurement points and the tooth surfaces of these points have been determined, and the data of point O can be obtained. Firstly, calculate the 4 points around point for linear interpolation. The value of the point O surrounded by point A,B,C,D is calculated by using the area weights of the interpolation values of the four triangles ABC,BCD,ABD,ACD. If point O is close to a region, the area weighting factor is close to one. If point O is far from the area, the area weighting factor is close to zero. When the interpolation point coincides with the measurement point, the interpolation value is the measured value itself.Figure 5
Four-point interpolation algorithm.As shown in Figure5, assumptions are as follows: the linear interpolation value where the triangles ABC pass through the interpolation point O is a1; the area of triangle ABC is SACOB, The interpolation value of triangle BCD is a2; the area of triangle BCD is SOCDB. Thus the interpolation value is b1.(3)b1=SACOBSACDBa2+SOCDBSACDBa1,where SACDB is the area of quadrilateral ABCD.On the other hand, the triangularABD and ACD’s interpolation values for point O are a3,a4 and the triangular ABD and ACD’s areas are SOABD,SACDO, respectively. Thus the interpolation value is b2.(4)b2=SOABDSACDBa4+SACDOSACDBa3.The final interpolation value b is(5)b=b1+b22.


## 2.1. Synthesis of Error Tooth Surface

In this paper, a base section of the gear is divided into 20 equal parts for gear contact analysis. The ideal contact part for the different size of gears is shown in Figure1.Figure 1
The geometrical contact line of the tooth surface during the meshing process.i: the corresponding contact line number during the engagement process; j: split points on the contact line.
(a)
The ideal contact part for spur gears
(b)
The ideal contact part for helical gearsThe interval between the contact lines is one twentieth of the rotation angle of a base section. The contact lines on the tooth surface of the large gear are numbered withi=1,2,…,N (N is the total number of contact lines) in the order of start contact to the end of the contact. The length of a contact line is divided into 19 equal parts, and the corresponding division points are numbered with j=1,2,…,20 in the order of the tooth root to the tooth top. The tooth surfaces meshing with the large gears are divided by the same method. Therefore, the reference plane position on the different size of gears can be expressed by the combination of the serial number of the contact line i and the division number j of the contact line. The shape deviation of the gear can be expressed as data 20×n.The sum of the deviation values of the normal direction at the corresponding contact point on the tooth surface is called the synthesis error. The curve of the 20 synthetic error values on a geometric contact line is called the synthetic error curve. In this paper, the synthetic error surface is expressed in a coordinate system with the split numberj of a geometric contact line being the X-axis and the contact line number i is the plane of the Y axis, which is shown in Figure 2. The direction of the synthetic error coordinate system e is positive and the highest point is zero. The minimum of the synthetic error on the geometric contact line is defined as the ridge. The projection of the ridge to the synthetic error coordinate system Y-e surface is called the ridge shape.Figure 2
The creation of synthetic error surfaces.


## 2.2. Tooth Surface Model Based on NURBS

NURBS surface is the only representation standard of product data model defined in STEP standard (Product data exchange standard), and it has been applied in real tooth surface simulation. However, in the process of real tooth surface reconstruction by using NURBS surface, the precision of the reconstructed tooth surface in the nongrid node is difficult to control due to the limited number of measurement data points. In order to ensure the smoothness of reconstructed tooth surface, this paper proposes a linear interpolation method combined with area weight interpolation method. This method is used to estimate the tooth surface data points which are outside the measurement area to realize the construction of the error tooth surface and to ensure the accuracy of the reconstructed tooth surface to the utmost extent.NURBS surfaces are tensor product forms based on NURBS curves, and the double three NURBS surface can be expressed as(1)Pu,v=∑i=1m∑j=1nBi,kuBj,kvωijVij∑i=1m∑j=1nBi,kuBj,kvωiju,v∈0,1,where m is the number of control vertices for u, n is the number of control vertices for v, and Vij is the control vertex of the surface. ωij is the weight factor of Vij. Bi,k is the kth spline basis function along u. Bj,k is the kth spline basis function along v. Among them, the basis function recursive formula is(2)Bi,0u=1ui≤u≤ui+10otherBi,ku=u-uiui+k-uiBi,k-1u+ui+k+1-uui+k+1-ui+1Bi+1,k-1uk≥1.


## 2.3. Delaunay Triangulation of Discrete Data

In this paper, the scanning method is used to obtain the tooth surface discrete point data. In order to describe the real tooth surface situation more intuitively, this paper deals with the discrete data points by means of Delaunay triangulation principle because the topological relations between the discrete data points are complex. The traditional Delaunay triangulation method often projects the three-dimensional discrete data points onto the two-dimensional plane and then triangles the projection points. However, the triangulation of this method cannot truly reflect its spatial angle, which will affect the quality of the section. On the basis of these methods, this paper presents a method of direct triangulation of three-dimensional discrete data points.Assuming that the scan line isli, k=1,2,…,n, the number of measurement points on the scan line is t, where two adjacent scanning lines are lk,lk+1,1≤k<k+1≤n. The measurement points on each of the two scanning lines are lk(j), j=1,2,…,t;lk+1l,l=1,2,…,t.As shown in Figure3, the starting points for connecting two scan lines are lk(1) and lk+1(1); the end point are lk(t) and lk+1(t). And find another scan line that is closest to the measurement point in the measurement points on lk and lk+1 and connect, where the area between the scanning lines lk and lk+1 is divided into r(r≥t-1) spatial regions. Finally, the triangular and spatial quadrilateral regions are retrieved.Figure 3
Preliminary Delaunay triangulation of adjacent tooth surface scan lines.For the spatial quadrilateral region, the minimum internal angle maximum criterion is used to complete the Delaunay triangulation of the discrete data points in this paper, respectively, as shown in Figure4.Figure 4
Polygon area triangulation.


## 2.4. Linear Interpolation Method and Area Weight Interpolation Method

In order to describe the error of the tooth surface quickly and accurately, it is necessary to block the discrete data points of the tooth surface in the construction of the digital real tooth surface. And use the interpolation method to accurately describe the error surface to reduce the amount of tooth interpolation calculation. The traditional data block method is often based on the geometric characteristics of the part, which does not apply to free surfaces that describe the error tooth surface. In this paper, an interpolation method based on linear interpolation combined with area weight interpolation is proposed to accurately identify tooth surface errors.The error tooth surface shape is calculated by using the area weight method according to the following interpolation method. In Figure5, assume that A,B,C,D are the measurement points and the tooth surfaces of these points have been determined, and the data of point O can be obtained. Firstly, calculate the 4 points around point for linear interpolation. The value of the point O surrounded by point A,B,C,D is calculated by using the area weights of the interpolation values of the four triangles ABC,BCD,ABD,ACD. If point O is close to a region, the area weighting factor is close to one. If point O is far from the area, the area weighting factor is close to zero. When the interpolation point coincides with the measurement point, the interpolation value is the measured value itself.Figure 5
Four-point interpolation algorithm.As shown in Figure5, assumptions are as follows: the linear interpolation value where the triangles ABC pass through the interpolation point O is a1; the area of triangle ABC is SACOB, The interpolation value of triangle BCD is a2; the area of triangle BCD is SOCDB. Thus the interpolation value is b1.(3)b1=SACOBSACDBa2+SOCDBSACDBa1,where SACDB is the area of quadrilateral ABCD.On the other hand, the triangularABD and ACD’s interpolation values for point O are a3,a4 and the triangular ABD and ACD’s areas are SOABD,SACDO, respectively. Thus the interpolation value is b2.(4)b2=SOABDSACDBa4+SACDOSACDBa3.The final interpolation value b is(5)b=b1+b22.


## 3. The Simulation of Contact Pattern

### 3.1. The Basic Theory of Contact Parts

When the tooth surface is transferring load, the contact tooth surface will be elastically deformed so that the area around the contact part will make contact. As the teeth continue meshing, the area left of the tooth surface is the actual tooth contact area, as shown in Figure6.Figure 6
Schematic diagram of contact parts.


### 3.2. Analytical Method of Contact Part of Involute Gear

As shown in Figure7, σ0i(i=1,2) represents two coordinate systems whose origin is O1,O2. koi coincides with i axis of the gear coordinate system, O1O2¯=E is the gear offset distance, and joi coincides with the vertical axis of the two axes. The two gears are rotated φ1,φ2 angles, respectively. The coordinate systems are σ1,σ2, respectively. A coordinate system σM(M;ας,αη,n) is established at the contact point M of the tooth surface Σ2, where n is the unit normal vector at point M on the tooth surface Σ2. ας and αη are two vertical unit vectors at the contact point M of the tooth face Σ2. Take any node Q2 on the tooth face Σ2. And draw a straight line from this point Q2 along a normal direction n to the tooth surface Σ. Then get the intersection Q1 of the straight line and the tooth surface Σ1.Figure 7
Conversion diagram of tooth surface meshing coordinate system.Coordinate transformation matrixσ2→σM is(6)T2M=a11a12a13a21a22a23a31a32a33,where a11=ας·i2, a12=ας·j2, a13=ας·k2, a21=αη·i2, a22=αη·j2, a23=αη·k2, a31=n·i2, a32=n·j2, a33=n·k2.According to Figure7(7)T2M=-cos⁡Σ0-sin⁡Σ010sin⁡Σ0-cos⁡Σ.Thus Coordinate transformation matrix σ2→σM is(8)T1M=T2M·Rk2,-φ2·RjO2,180∘+Σ·RkO1,φ1=T2Mcos⁡φ2sin⁡φ20-sin⁡φ2cos⁡φ20001-cos⁡Σ0-sin⁡Σ010sin⁡Σ0-cos⁡Σcos⁡φ1-sin⁡φ10sin⁡φ1cos⁡φ10001,where Rk2,-φ2 is the coordinate transformation matrix of σ2→σM; RjO2,180∘+Σ is the coordinate transformation matrix of σO1→σO2.R
k
O
1
,
φ
1 is the coordinate transformation matrix of σ1→σO1.The vector diameter of pointQ2 in the coordinate system σM is(9)rO2M=r2M+MO2M=T2Mr22+T2MMO22.The vector diameter of point Q1 in the coordinate system σM is(10)rO1M=r1M+MO1M=r1M+MO2M+O1O2M,where (r2)2 is the r2 vector in coordinate system σ2 and (MO2)2 is a vector of point O2 to point M in coordinate system σ2. (r1)M is the r1 vector in coordinate system σM. (MO2)M is a vector of point O2 to point M in coordinate system σM. (O1O2)M is a vector of point O1 to point O2 in coordinate system σM. The expression of vector (O1O2)M is(11)O1O2M=T2MRk2,-φ20-O1O20,where Rk2,-φ2 is the coordinate transformation matrix of σ2→σM.In coordinate systemσM, the projection for difference between MQ2 and MQ1 in the normal vector direction n1 at the point M on the tooth surface Σ2 is the clearance in the direction of the engagement of the two tooth surfaces at the point Q1 and point Q2. When the pinion turns through angle φ1, the gear turns through angle φ2:(12)rQ1M=xQ1,yQ1,zQ1rQ2M=xQ2,yQ2,zQ2.In order to ensure that the pointQ2 on the tooth surface Σ2 of the large gear is still the point Q1 on the tooth surface Σ1 of the pinion in the direction of n, the vectors rQ1 and rQ2 must satisfy the following relationships:(13)xQ2-xQ1≤10-4yQ2-yQ1≤10-4.The clearance between the two engaging tooth surfacesΣ1 and Σ2 in the direction of n at points Q1 and point Q2 is(14)Δz=Q1Q2¯=zQ2-zQ1.According to GB/Z18620.4-2008, the thickness of the detection coating for contact part is set to 0.01 mm. The light load contact part exists in the area of contact point separation amountδ1=0.01 mm. When Δz≤0.01 mm, it is considered that the two tooth surfaces Σ1 and Σ2 are in contact at point Q1 and point Q2. By dividing the tooth surface into the grid, the clearance in the direction of n the corresponding grid node is calculated, so is the tooth contact part.


## 3.1. The Basic Theory of Contact Parts

When the tooth surface is transferring load, the contact tooth surface will be elastically deformed so that the area around the contact part will make contact. As the teeth continue meshing, the area left of the tooth surface is the actual tooth contact area, as shown in Figure6.Figure 6
Schematic diagram of contact parts.


## 3.2. Analytical Method of Contact Part of Involute Gear

As shown in Figure7, σ0i(i=1,2) represents two coordinate systems whose origin is O1,O2. koi coincides with i axis of the gear coordinate system, O1O2¯=E is the gear offset distance, and joi coincides with the vertical axis of the two axes. The two gears are rotated φ1,φ2 angles, respectively. The coordinate systems are σ1,σ2, respectively. A coordinate system σM(M;ας,αη,n) is established at the contact point M of the tooth surface Σ2, where n is the unit normal vector at point M on the tooth surface Σ2. ας and αη are two vertical unit vectors at the contact point M of the tooth face Σ2. Take any node Q2 on the tooth face Σ2. And draw a straight line from this point Q2 along a normal direction n to the tooth surface Σ. Then get the intersection Q1 of the straight line and the tooth surface Σ1.Figure 7
Conversion diagram of tooth surface meshing coordinate system.Coordinate transformation matrixσ2→σM is(6)T2M=a11a12a13a21a22a23a31a32a33,where a11=ας·i2, a12=ας·j2, a13=ας·k2, a21=αη·i2, a22=αη·j2, a23=αη·k2, a31=n·i2, a32=n·j2, a33=n·k2.According to Figure7(7)T2M=-cos⁡Σ0-sin⁡Σ010sin⁡Σ0-cos⁡Σ.Thus Coordinate transformation matrix σ2→σM is(8)T1M=T2M·Rk2,-φ2·RjO2,180∘+Σ·RkO1,φ1=T2Mcos⁡φ2sin⁡φ20-sin⁡φ2cos⁡φ20001-cos⁡Σ0-sin⁡Σ010sin⁡Σ0-cos⁡Σcos⁡φ1-sin⁡φ10sin⁡φ1cos⁡φ10001,where Rk2,-φ2 is the coordinate transformation matrix of σ2→σM; RjO2,180∘+Σ is the coordinate transformation matrix of σO1→σO2.R
k
O
1
,
φ
1 is the coordinate transformation matrix of σ1→σO1.The vector diameter of pointQ2 in the coordinate system σM is(9)rO2M=r2M+MO2M=T2Mr22+T2MMO22.The vector diameter of point Q1 in the coordinate system σM is(10)rO1M=r1M+MO1M=r1M+MO2M+O1O2M,where (r2)2 is the r2 vector in coordinate system σ2 and (MO2)2 is a vector of point O2 to point M in coordinate system σ2. (r1)M is the r1 vector in coordinate system σM. (MO2)M is a vector of point O2 to point M in coordinate system σM. (O1O2)M is a vector of point O1 to point O2 in coordinate system σM. The expression of vector (O1O2)M is(11)O1O2M=T2MRk2,-φ20-O1O20,where Rk2,-φ2 is the coordinate transformation matrix of σ2→σM.In coordinate systemσM, the projection for difference between MQ2 and MQ1 in the normal vector direction n1 at the point M on the tooth surface Σ2 is the clearance in the direction of the engagement of the two tooth surfaces at the point Q1 and point Q2. When the pinion turns through angle φ1, the gear turns through angle φ2:(12)rQ1M=xQ1,yQ1,zQ1rQ2M=xQ2,yQ2,zQ2.In order to ensure that the pointQ2 on the tooth surface Σ2 of the large gear is still the point Q1 on the tooth surface Σ1 of the pinion in the direction of n, the vectors rQ1 and rQ2 must satisfy the following relationships:(13)xQ2-xQ1≤10-4yQ2-yQ1≤10-4.The clearance between the two engaging tooth surfacesΣ1 and Σ2 in the direction of n at points Q1 and point Q2 is(14)Δz=Q1Q2¯=zQ2-zQ1.According to GB/Z18620.4-2008, the thickness of the detection coating for contact part is set to 0.01 mm. The light load contact part exists in the area of contact point separation amountδ1=0.01 mm. When Δz≤0.01 mm, it is considered that the two tooth surfaces Σ1 and Σ2 are in contact at point Q1 and point Q2. By dividing the tooth surface into the grid, the clearance in the direction of n the corresponding grid node is calculated, so is the tooth contact part.


## 4. Prediction of Dynamic Performance of Error Tooth Surface

The main causes of gear vibration are gear manufacturing error, installation error, gear time varying meshing stiffness and load gear, shaft, gearbox, and other deformation caused by the gear position error. The gear contact part and transmission error is the visualization of gear meshing performance. In the example, the reliability of the prediction results is judged by measuring the vibration of the gears in the direction of the meshing line.


### 4.1. Gear Transmission Error

The factors that affect the dynamic excitation of the gear are not only the structural form of the gear itself, the geometric characteristics, and the error, but also the other components in the system. The cylindrical gear system can be simplified as a torsional vibration system for the gear pair when the stiffness of the support member such as the drive shaft, bearing, and housing is large. The dynamics model is shown in Figure8.Figure 8
Torsional vibration model of gear pair.Suppose the gear pair of coincidence between one and two. According to the mechanical dynamics theory and the mechanical vibration theory, the torsional vibration equation of a pair of gear pair is(15)Md2xdt2+cmdxdt+kmx=W0+Wd,where(16)x=r1θ1-r2θ2M=I1I2/I1r22+I2r12W0=T1r1=T2r2Wd=k1e1+k2e2+c1de1dt+c2de2dt.In the above formula,θi(i=1,2) is the torsional vibration displacement of gear and pinion.Ii(i=1,2) is the inertia of the gear and pinion.ri(i=1,2) is the base radius of the gear and pinion.ki(i=1,2) is the integrated stiffness of the pair of ith gears.ci(i=1,2) is the damping coefficient of the ith gear.k(t) is the combined meshing stiffness of gear pair.c is meshing damping of gear pair.ei(i=1,2) is the error of the ith gear.Ti(i=1,2) is the external load moment on the driving and driven gears.When the dynamic term such as inertia force and damping force term is ignored, the static transmission error is(17)xs=W0kt+k1e1+k2e2kt.In the quasi-static transmission error test of the gear pair, it is necessary to know exactly the rotation angle of the driving and driven gears in the transmission process.


## 4.1. Gear Transmission Error

The factors that affect the dynamic excitation of the gear are not only the structural form of the gear itself, the geometric characteristics, and the error, but also the other components in the system. The cylindrical gear system can be simplified as a torsional vibration system for the gear pair when the stiffness of the support member such as the drive shaft, bearing, and housing is large. The dynamics model is shown in Figure8.Figure 8
Torsional vibration model of gear pair.Suppose the gear pair of coincidence between one and two. According to the mechanical dynamics theory and the mechanical vibration theory, the torsional vibration equation of a pair of gear pair is(15)Md2xdt2+cmdxdt+kmx=W0+Wd,where(16)x=r1θ1-r2θ2M=I1I2/I1r22+I2r12W0=T1r1=T2r2Wd=k1e1+k2e2+c1de1dt+c2de2dt.In the above formula,θi(i=1,2) is the torsional vibration displacement of gear and pinion.Ii(i=1,2) is the inertia of the gear and pinion.ri(i=1,2) is the base radius of the gear and pinion.ki(i=1,2) is the integrated stiffness of the pair of ith gears.ci(i=1,2) is the damping coefficient of the ith gear.k(t) is the combined meshing stiffness of gear pair.c is meshing damping of gear pair.ei(i=1,2) is the error of the ith gear.Ti(i=1,2) is the external load moment on the driving and driven gears.When the dynamic term such as inertia force and damping force term is ignored, the static transmission error is(17)xs=W0kt+k1e1+k2e2kt.In the quasi-static transmission error test of the gear pair, it is necessary to know exactly the rotation angle of the driving and driven gears in the transmission process.


## 5. Experiment Study

This paper chooses the gears of a certain type of solar generator gearbox to carry on the research. According to the support stiffness of the gear and other factors, refer to optimizing the gear modification [11–27]. This paper would not study the best trimming tooth surface; that is, select two measurements modification gears; it is short for group A and group B in the following paper, and its dynamic performance in analysis.The transmission mechanism in the gearbox is composed of three spur gears, and the structure of the gearbox is shown in Figure9. The two gears with 35 teeth are the pinion gear, and the gears with the number of teeth are 42. The speed of the drive wheel is 1800 r/min and the rated power is 25 KW. The basic parameters of the gear are shown in Table 1.Table 1
Factors and their levels.
Parameters
Letter symbols
Gear modulus
m
2.5
Pressure angle
α
20
Addendum
h
a
2.5
Dedendum
h
f
3.125
Teeth
Z
35/42Figure 9
Gear distribution diagram.Using the measurement method of ZEISS coordinate measuring instrument, the topological measurement of tooth width and tooth height is taken. According to the actual measurement data, the surface of the tooth surface deviation is directly constructed, and the shape deviation surface in the measurement area of the two sets of gear driving wheels is shown in Figure10. Through the calculation of the points of the measurement area through the interpolation algorithm, tooth surface with the shape deviation on the different mesh nodes is obtained, which are shown in Figure 11.Figure 10
The tooth surface with shape deviation in the measurement area.Figure 11
The tooth surface with shape deviation in the measurement area after interpolation.On the basis of the analysis of the contact part of the error tooth surface after interpolation, the contact surface of the tooth surface is shown in Figure12.Figure 12
Tooth surface contact parts (theoretical calculation).The theoretical calculation of the transmission error is shown in Figure13.Figure 13
Gear drive error (theoretical calculation).A gear loading test is performed for gears whose tooth surfaces are coated with steel red after steel red is dried, and the gear contact parts are shown in Figure14, The transmission error of the gear is shown in Figure 15.Figure 14
Contact part of test gear.Figure 15
Transmission error of test gear.In order to test the accuracy of the structural error surface, the vibration measuring instrument is used to measure the time-domain diagram of the vibration in the direction of the gear meshing line, and its FFT is shown in Figure16.Figure 16
Gearbox vibration spectrum.It can be seen from the figure that the predictive gear contact part and gear transmission error is basically the same as test results. The vibration measurement of the gearbox along the gear meshing line shows that the vibration condition of the gearbox is basically the same as the transmission error of the gear. This shows that the model of error tooth surface is accurate.


## 6. Conclusions

(1)
An algorithm combining linear interpolation with area weight interpolation is proposed to realize the construction of error tooth surface. This method not only reduces the computational complexity of the tooth surface interpolation algorithm, but also improves the efficiency of the digital tooth surface construction.
(2)
The contact part and the transmission error of the gear are analyzed by using the constructed error tooth surface, which provides a new idea for the prediction of the dynamic performance of the gear.
(3)
Through the gear loading test, the contact part and gear transmission error of the gear are obtained. The analysis result and the test result are basically the same, which shows the accuracy of the constructed error tooth surface.
(4)
The construction of the error tooth surface in this paper lays the foundation for further improving the grinding accuracy of the tooth surface and the grinding of the high performance gear.



---

*Source: 6143054-2017-12-12.xml*