Abstract

The manufacturing error of spiral bevel gear tooth surface has a great influence on transmission efficiency and gear life. The error of the gear tooth surface needs to be measured accurately and fed back to the machine tool to adjust the parameters. When measuring the spiral bevel gear using a gear measuring machine, combined with the measurement theory of the tooth flank of spiral bevel gear, this paper proposed a method to compensate the inclination error in the measurement result precisely. Based on the iterative search method, a precision matching method for the theoretical and the measured tooth surface of the spiral bevel gear was designed to calculate the compensation results. The experimental results show that the inclination errors included in tooth surface measurement results reduced from more than 3 μm to less than 0.5 μm, and more than 70% of the errors are compensated by the proposed method. The accuracy of the measurement results improved significantly after compensation, and furthermore, it can provide a more accurate basis for the adjustment of machine tool parameters in the manufacturing process.

1. Introduction

Spiral bevel gear is widely used in the field of mechanical transmissions, such as automobile, aviation, wind turbines, and so on. When manufacturing the gear by machine tools, the manufacturing error affects the accuracy of the gear tooth surface and then directly affects the performance of gear transmission efficiency, noise, lifetime, and reliability [14]. In the closed-loop manufacturing system, in order to feedback and adjust the machine tool parameters, the gear tooth surface requires to improve its measuring accuracy. At present, the coordinate measurement method is one of the most widely used methods to measure the tooth surface error of the spiral bevel gear [57].

When measuring the tooth surface of spiral bevel gear by gear measuring machine (GMM), the manufacturing and assembly error of the gear is inevitable. The inclination error caused by the nonparallel between the gear axis and the vertical rotation axis of GMM will affect the measurement results of spiral bevel gear [8, 9]. Therefore, to improve the measurement accuracy of the spiral bevel gear tooth surface, first of all, the measured tooth surface of the gear and its theoretical tooth surface must be accurately matched to eliminate the inclination error in the measurement result of the gear tooth surface [10].

Many works have been done on the measurement error and the compensation method of spiral bevel gear tooth surface. In order to predict the relative errors of alignment between spiral bevel gears, Fuentes et al. proposed a procedure of determination of the relative spatial position of spiral bevel gear supporting shafts during torque transmission [1113]. The misalignment problem for the spiral bevel and hypoid gears usually exists in the actual manufacturing and transmission, and Ding et al. presented an automatic data-driven operation and optimization to determine the uncertain misalignment [14]. Considering residual tooth flank form error, Shao et al. presented an accurate systematic CMM measurement method to prescribe and data-driven control the tooth flank form error and get a flexible compensation of the error [15]. Fu et al. proposed a new online detection method of flatness of spiral bevel gears, and the data of flatness is evaluated by the least square method after error separation [16]. When manufacturing the spiral bevel gears, Shih and You developed an on-machine measurement system on five-axis machines using a quasi-3D probe and verified by comparing the evaluation results with the same gear measured on a dedicated GMM [17]. By using the virtual conjugate reference tooth surface, Li et al. presented a novel scanning measurement method to measure the tooth flank form of hypoid gears [18, 19]. Wang et al. established the on-machine measuring and the data processing method for spiral bevel gears, and the method was validated by comparing simulation results with the actual measuring result [20, 21]. Simon researched the influence of misalignments of the mating members and tooth errors on mesh performances of the spiral bevel gear and found that the misalignments and the tooth spacing error could increase the angular position error of the driven gear [22]. By considering spiral bevel gear geometry, Shunmugam reported a method of determining the normal deviation from the theoretical surface of the bevel gear and then gave the validation of the method [23]. Peng et al. researched the effects of eccentricity error on the gear dynamic responses and proposed a simple method using translational kinematic transmission error modification to reduce the computational time [24]. Using the truncated singular value decomposition and the L curve method, Chen and Yan solved the identification equation of the tooth surface deviation and compared it with the least square method to verify the validation of this method [25]. Xie et al. presented a coordinate measuring method for the variable ratio noncircular bevel gear by coordinate measuring machine [26]. Combined with the theoretical gear and coordinate measuring machine, Cao and Deng provided a method for measuring the errors between the actual surface and the theoretical surface of the bevel gear and then developed the corrective machine setting [27]. Zhou et al. investigated a higher-accuracy fitting method for the form error tooth flank, and the universal machine tool settings are exploited for the identification of the real tooth flank form error [28]. Mo and Zhang gained the corresponding digitized true tooth surface by changing the machining adjustment parameters, which can lay a solid basis for the subsequent finite element analysis of gear contact and transmission error analysis [29]. However, few works have been focused on the measurement error of the spiral bevel gear caused by the inclination error, so it is necessary to research the compensation method for the inclination errors.

Combined with the measurement theory of the tooth flank of spiral bevel gear, this paper proposed a precision matching method for the theoretical tooth surface and the measured tooth surface of the spiral bevel gear, and then compensated the inclination error in the measurement result precisely. The experimental results show that the accuracy of the measurement results improved significantly after compensation.

2. Establishment of the Measuring Coordinate System

The design basis of spiral bevel gear is the apex of the pitch cone of the bevel gear, and its design coordinate system is set as {O: X, Y, Z}. When measuring spiral bevel gear with GMM, the inner end of the gear was fixed upward, and the outer end was fixed downward on the platform, and the measurement basis is the rotation axis of the gear and the outer end surface of the gear.

As shown in Figure 1, the measuring coordinate system of spiral bevel gear {Om: Xm, Ym, Zm} is set to coincide with the coordinate system {O: X, Y, Z}. According to the above coordinate systems, the position of the theoretical gear pitch cone apex in the Z-axis direction can be determined by the Z-axis coordinate value of the mounting platform and the mounting distance of the gear, and the coordinate relationship between the measured pitch cone apex and the theoretical pitch cone apex can be solved.

3. Precision Compensation Method for the Measurement Errors of Gear Tooth Surface

When measuring the tooth surface of spiral bevel gear by GMM, the theoretical data point cloud and the corresponding unit normal vector of the gear tooth surface measuring points should be input GMM first. According to the coordinate value of the theoretical data point cloud, the GMM measures the actual coordinates of the points on the tooth surface and finally obtains the actual measurement error. In order to compensate the inclination error of the rotating axis of the gear, matching the measured tooth surface with the theoretical tooth surface of the gear accurately is the first step.

The spatial relative positions of two curved surfaces can be represented by six relative position parameters of the corresponding points P and P1: {Δx, Δy, Δz, Δθx, Δθy, Δθz}. When matching the two curved surfaces, the above six parameters should be calculated first, and then one of the two surfaces translated and rotated according to the value of the six parameters to match the other surface in spatial. This paper mainly focuses on the translation deviations Δx and Δy along the X and Y axis, and the angle deviations Δθx and Δθy around the X and Y axis, and finally, find a compensation method for the above four parameters.

Assuming that the point on the measured tooth surface is recorded as point Pm (xm, ym, zm), and the corresponding point on the theoretical tooth surface is named as point Pt (xt, yt, zt). After the surface matching process, the measured tooth surface moves to a new position and the corresponding point of Pm is recorded as point Pc (xc, yc, zc). The distance from the point Pm to the theoretical tooth surface is Δd, and the intersection between the normal vector of the theoretical tooth surface that passing through the point Pc and the theoretical tooth surface is point Ps (xs, ys, zs). The objective function of the matching method is shown in the following:

Due to the inclination error, the position relationship between the measured tooth surface and the theoretical tooth surface is uncertain. Therefore, the position of the point Ps corresponding to the point Pm needs to be obtained by the search calculation method. As the error of the actual tooth surface is generally small, the position of the point Ps is considered nearby the point Pt, as shown in Figure 2.

If the corresponding normal vector of point Ps passes through the point Pm, the point Ps is the closest point from the point Pm to the theoretical tooth surface, and the cosine value of the angle between the vectors and is 1. In order to find the position of the nearest point Ps by searching, eight measured points Pi (xi, yi, zi) (i = 1∼8) around the point Pt and their corresponding unit normal vectors are selected to form the search region, and the objective function is established as in equation (2). The position of the point Ps could be searched by iteratively calculating and comparing the value of the function f of the point Pi.

This paper researches the single direction matching method with the example of tooth surface matched by moving along the X-axis and rotating around the X-axis.

3.1. Matching Method for the Measured and Its Theoretical Tooth Surface in Single Direction

The iterative search calculation diagram is shown in Figure 3. The point Pm is taken as the origin of the coordinate system, the moving distance and rotating angle are taken as the X-axis and Y-axis, and the optimal matching point Pc of the two surfaces is searched by the iteration method in this figure.

First, set the initial search step Δx0 and Δθ0 to ensure that the initial search region could enclose the point Pc. Set the matching parameters Δx = 0, and Δθx = 0. The specific process is as follows:(1)Move the point Pm with a distance Δx0 along with the positive and negative directions of the X-axis, a distance Δθ0 along with the positive and negative directions of the θx axis, respectively. Together with the point Pm, nine new matched points are obtained and named P0P8. According to equation (1), the objective function value Fi (i = 0 − 8) could be calculated. Compare the value Fi (i = 0 − 8), and the new search region is judged according to the minimum value. For example, if F1 is the minimum value, the new search region is the first quadrant shown in Figure 3.(2)The new search direction is determined according to the new search region, and the step size is reduced to half of the original step, that is and . If the point Pc is in the first quadrant, the coefficients corresponding to the search direction are recorded as α1 = 1 and β1 = 1. If the point Pc is in the second quadrant, the coefficients is recorded as α1 = -1 and β1 = 1, and so on. In order to ensure that the optimization point is always in the search region, when a new search process is carried out, the search region is expanded with a distance of and , where r is the ratio. Generally, the value of r is 0.5. The new step size can be calculated by the following:(3)Along the search direction, move the point Pm with a distance and to get a new matched point Pc1. The distance could be calculated by the following:Take the point Pc1 as the center, and Δx1 and Δθx1 as the step size, repeat steps (1) and (2) to get the new search region, the search direction coefficients α2 and β2, the new step size Δx2 and Δθx2, and the moving distance and .(4)After the k-th iteration, the moving distance of the matched point is calculated. Move the point according to the moving distance and to get the new matched point Pck, and then repeat step (1).(5)After the k-th iteration, if the steps Δxk along the X-axis and Δθxk around the X-axis are less than ε, the iteration search results meet the accuracy requirements, and the iteration ended.

After the iteration, the matched results of the two tooth surfaces could be obtained, and the compensation value for the tooth surface measurement errors of the spiral bevel gear could be calculated by the following:

The specific iteration process is shown in Figure 4.

3.2. Matching Method for the Measured and Its Theoretical Tooth Surface in Multiple Directions

In this iteration process, four parameters Δx, Δy, Δθx and Δθy should be calculated. In order to reduce the calculation work, the initial iteration (coarse matching) is carried out first, and then the accurate iteration (precision matching) in two directions is carried out. The specific process is shown in Figure 5.

The coarse matching process is as follows: According to the matching method for a single direction, match the tooth surface along and around X-axis for n times and then match the tooth surface along and around Y-axis for n times. Judge whether the tooth surface deviation is less than the end criterion of the iteration. If the deviation is less than the criterion, skip the cycle and turn it into the accurate iteration process; otherwise, continue the iteration.

The accurate iteration process is as follows: first, match the tooth surface along and around X-axis for 1 time according to the matching method for a single direction, and then match the surface along and around Y-axis for 1 time. Judge whether the deviation meets the end criterion of the iteration. If it does not meet the criterion, return and continue the iteration until the criterion is met, and then calculate the matching results of the tooth surface according to the matching method.

Finally, the compensation value of the tooth surface measurement errors of the spiral bevel gear in two directions could be calculated by equation (5).

4. Experiment Results and Analysis of the Compensation Method

4.1. Basic Data of the Measured Gear

The parameters of the measured spiral bevel gear are shown in Table 1.

The theoretical points of the tooth surface are selected along three profile lines named PL7, PL15, and PL23, and one tooth trace line named TL5, as shown in Figure 6. Each line has 113 theoretical points. The theoretical data both for the convex and the concave surfaces contain the coordinates and unit normal vectors of these theoretical points. According to these theoretical data, actual measurement experiments are done and the measured tooth surface are the convex and concave surfaces of Tooth 1, Tooth 15, Tooth 29, and Tooth 43 of the gear. In the actual measurement experiment, some assembly errors are added artificially.

According to the compensation method proposed in the paper, the point cloud of the measured tooth surface could be matched with the theoretical tooth surface, and the compensation values could also be calculated by equation (5). Compared with the errors of the measured value, the accuracy of the compensation could be calculated.

4.2. The Compensation Results and Analysis

Figure 7 shows the comparisons of the measurement errors of the measured points along the tooth trace line both on the concave and convex tooth surface of Tooth 1 of the gear before and after compensation. Figures 7(a) and 7(b) are the compensation results of the convex tooth surface in the X direction and Y direction, respectively, and Figures 7(c) and 7(d) are the compensation results of the concave tooth surface in the X direction and Y direction, respectively.

After matching and compensation between the measured tooth surface and the theoretical tooth surface, the average measurement errors of the tooth surface along the X-axis and Y-axis of all measured tooth surfaces are shown in Tables 2 and 3, and the corresponding figures of the average errors of four teeth of the gear along the X-axis and Y-axis before and after compensation are shown in Figures 8 and 9, respectively.

From Tables 2 and 3, Figures 8 and 9, it can be easily seen that, before the compensation process, the average errors of some tooth flank are more than 3 μm. After matching and the compensation process, the average errors of the measurement results of all tooth surfaces are less than 0.5 μm, and more than 70% of the measurement errors are compensated by the proposed method.

5. Conclusions

In order to improve the accuracy of the measurement results of the spiral bevel gear tooth surface, this paper presents a precision compensation method for the inclination error of the gear tooth surface. Based on the iterative search method, this paper established an objective function to match the theoretical tooth surface and actual tooth surface along and around X-axis and Y-axis accurately and then compensated the measurement errors of the tooth surface according to the matching method.

The experimental results show that the errors included in the tooth surface measurement results along and around the X-axis, Y-axis are reduced from more than 3 μm to less than 0.5 μm, and more than 70% of the errors are compensated by the proposed compensation method. The results show that the measurement errors in the above directions can be eliminated obviously, and this method can provide a more accurate basis for the adjustment of machine tool parameters in the manufacturing process.

Data Availability

The processed data cannot be shared at this time as the data also form part of an ongoing study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported in part by the Natural Science Foundation of Shaanxi Province (2019JQ-556), in part by the Open Research Project of Key Laboratory of Advanced Manufacture Technology for Automobile Parts (Chongqing University of Technology), MOE (2017KLMT06), in part by the Major Science and Technology Projects of Shaanxi Province (2018zdzx01-01-01), and in part by the Fundamental Research Funds for the Central Universities, CHD (300102250111).