Abstract
To increase quality, reduce cycloidal gear noise, and avoid unnecessary vibration and shock, a compensation of axial geometric errors method is proposed based on the cycloidal gear form grinding. In the process of machining cycloidal gears, the relative position relationship between the grinding wheel and workpiece is affected by geometric errors of the motion axes, which has serious effects on the surface accuracy of the cycloidal gears. Combined with cycloidal gear form grinding kinematic principles, a geometric error model for each axis of a four-axis computer numerical control form grinding machine is established. By changing the compensation value of the geometrical errors on six degrees of freedom, the error of the cycloid gear tooth surface machined is obtained. Based on a sensitivity analysis of geometrical errors of each axis, the corrections are determined through an optimization process that targets the minimization of the tooth flank errors. The geometric errors of each axis of the cycloid gear grinding machine are compensated, and then, the cycloid gears produced by the machine are processed. Through the processing experiment, the error data of the actual processing before and after the compensation are compared, which indicates that the machining accuracy of the cycloid gear grinding machine is obviously improved. It has an important guiding significance in improving the precision and performance of large CNC form gear grinding machines.
1. Introduction
Cycloidal gears are used in rotating vector reducers, instruments and apparatus, planetary gearboxes, and other mechanical transmission devices. A rotary vector reducer is usually referred to as an RV reducer. Compared with mechanical transmissions such as friction, hydraulic pressure, and belts, cycloidal gear transmission efficiency is higher and has a wide range. At the same time, the life of the cycloidal gears is longer, its transmission ratio is stable, and its work reliability is also higher. Therefore, cycloidal gears have necessary and optional applications in various machinery manufacturing industries such as RV reducer, planetary gearbox, and construction machinery. As an important part widely used in mechanical transmission, the quality and accuracy of cycloidal gears during design and processing will directly affect the quality and reliability of mechanical products. To meet the demand for these kinds of gears, the simulation and compensation of axial geometric errors for cycloidal gear grinding are applied to the form grinding machine to increase its load capacity and reduce the transmission error of the RV reducer device. Therefore, how to improve the machining accuracy and optimize the design of cycloidal gears has always been a research hot spot. The problem has been paid close attention to by the business community and academia of various countries [1, 2].
In recent years, with the rising labor costs and the continuous adjustment of the industrial structure, the development trend of the manufacturing industry has become an intelligent manufacturing technology with robots as the core, which has a very important role in promoting the upgrading and transformation of the industry. The RV reducer is a very important mechanical device widely used in the robot industry. At each arm joint of the robot, the main shaft and the arm are connected, which can reduce the rotation speed of the arm and increase the torque. During the process, the precision of the cycloidal gears is closely related to the precision and quality of the RV reducer. One of the most important parts of the RV reducer is the cycloid gear. As an important core component in the RV reducer, the cycloid gears determine the rotation accuracy of the whole reducer, so the machining accuracy of the cycloid gears is also the same as that of the robot important indicators for design and research. In this case, it is very important to manufacture and process the cycloid gear grinding machine based on the forming method. It is of great significance to improve the accuracy of the cycloidal gears [3–5]. Therefore, the research on the precision of the cycloid gear grinding machine is also the focus.
Schultschik [6] puts forward the concept of machine tool space error for the first time, using the closed-loop constraints of the machine tool's mechanical structure. Denavit et al. [7] proposed the ideal dynamic model of multiaxis machine tools based on a standard homogeneous transformation matrix. Paul [8] incorporated robot mathematics, programming, and control into the basic manipulator’s idea of the reference frame and modified the model. Hocken [9] proposed the static positioning error modeling of the measuring machine by the matrix transformation method and measured the three-dimensional measuring machine using the three-dimensional measuring technology. Since the 1980s, scholars have begun to study the machining error caused by the spatial error of the machine tool and established a spatial error model including the manufacturing geometric error [10]. Kim [11] established an error compensation model for a three-axis machine tool, through which the axis motion deviation can be predicted within a certain accuracy range. Soons et al. [12] introduced the model of the rotation axis in the rigid body kinematics model and conducted a comprehensive error model study of multi-axis machine tools including translation and rotation axes. Rahman et al. [13] used the homogeneous coordinate transformation matrix to study the quasi-static error of the machine tool, including the thermal error and geometric error of the rotary axis and the translation axis. Flourssen et al. [14] proposed a new scheme of the three-dimensional hemisphere to identify the error of the rotary axis in a five-axis machine tool by using three-dimensional length measurements. Abbaszadeh et al. [15] classified the geometric errors of multiple motion axes in the machine tool into static errors and motion errors. Tsutsumi et al. [16,17] conducted a mathematical analysis of the influence of installation deviation on the parameter identification of the simplified model of the rotating shaft and simulated the simplified model identification algorithm, which proved the reliability of the simplified method. Zargarbashihe and Mayer [18] analyzed the static and kinematic geometric errors of the first rotary axis of a five-axis machine tool, and mainly studied the five motion errors of the first rotary axis. Lei et al. [19] analyzed the motion error of five-axis machine tools, designed a measurement scheme that is not disturbed by the linear axis error model, and modeled the backlash error and servo mismatch. Sharif and Ibaraki et al. [20] studied the detection scheme of completing all measurements in one installation without considering the geometric error of the translation axis, reducing the uncertain identification results caused by the need for multiple installations and multiple separate modeling. Ibaraki et al. [21] used the probe to measure the deviation between the actual motion position and the theoretical position of the rotary axis of the machine tool during the machining process through the spatial position relationship among the static errors of each axis of the five-axis machine tool. Lee et al. [22] carried out the quintic polynomial modeling of the rotating shaft by limiting some polynomial coefficients, which simplified the identification difficulty and improved the accuracy of error estimation. Tsutsumi et al. [23] proposed two error identification schemes of different coordinate systems for double turntable machine tools, respectively, using a Cartesian coordinate system and cylindrical coordinate system to conduct experiments on the static error of double turntables. Mchichi and Mayer [24] analyzed the inherent errors of the machine tool and identified and analyzed the errors through standard ball arrays and special probes. Fan et al. [25] focused on the unified kinematic modeling method of three-axis machine tools, modularized each kinematic pair model, and studied the commonalities in modeling of three-axis machine tools with four different structural forms. The concept of exponential product and spiral in robotics is used to model the kinematics of the 5-axis machine tool. It is pointed out that the method does not need to establish a fixed coordinate system for each motion pair [26]. An innovative geometric error compensation method for machining with non-rotary cutters is proposed, in which the cutter rotation angle is considered for both modeling and compensation of geometric errors, and the instantaneous ideal contact point is calculated based on the worm grinding of face gears generation process [27–29]. The insufficient machining accuracy of cycloid gear will lead to the increase in transmission error of gear pair and the stress concentration in meshing, which will lead to bending fatigue and contact fatigue, and other problems. These problems will make gear wear and the machine operating results of deviation, unable to achieve the desired effect, and even cause safety accidents in serious cases [30, 31].
This work aims to establish a model of axial geometric errors for cycloidal gear grinding machines based on form grinding. To improve the machining accuracy of machine tools to a large extent, the geometric error compensation technology of machine tools is indispensable, and how to accurately establish the motion error model of machine tools is one of the key issues especially in the cycloidal gear grinding machine. After establishing the numerical control model of the geometric errors of the cycloidal gear grinding machine, the theoretical tooth surface of the cycloidal gears can be compared with the error tooth surface through the actual tooth surface point coordinates. Then, the actual geometric errors of the machine tool are measured by the error measuring tool, and the error is compensated. Finally, experimental results will show that this method has application value in improving the accuracy of machining the cycloidal gears.
The outline of the remainder of this study is as follows. In Section 2, a topological structure analysis of machine tool multibody system based on the gear form grinding is developed. Geometric error modeling of cycloidal gear grinding machine is the main subject of Section 3, Section 4, and Section 5. The cycloidal gears processed and numerical examples are illustrated in Section 6. Finally, some conclusions are drawn in Section 7.
2. Topological Structure Analysis of Machine Tool Multi-Body System
The main research is based on the cycloidal gear form grinding machine. This method of forming a tooth surface can ensure that when grinding two or three cycloidal gears used in the same RV reducer, the parameters such as the tooth shape of the processed cycloidal gears are the same. Compared with the method of grinding gears by the generating grinding method, the form gear grinding method has greater advantages in terms of grinding accuracy, grinding efficiency, cost, and machine tool structure. The correction grinding wheel of the cycloidal gear form grinding machine uses a diamond wheel so that the cycloidal gear within 500 mm can be ground. Gear form grinding is arranged based on the universal CNC gear form grinding machine, which has four digital servo closed-loop controlled axes: three rectilinear motions (axis X, Y, and Z) and one rotational motion (axis C) defined in Figure 1. SP1 and SP2 are spindles of the grinding wheel and diamond wheel, respectively.

Based on the 3D model of the cycloidal gear grinding machine, the kinematic analysis of the gear grinding system of the research object is carried out according to the theory of multibody system, and the relations between the adjacent low-sequence bodies of the moving parts are determined. The topology diagram of the multibody system and the coordinate system of the machine tool are shown in Figure 2. The machine tool bed is taken as the matrix . The gear grinding system can be divided into two motion chains, for example, the workpiece motion chain extending from the machine tool bed to the worktable and the tool motion chain extending from the fuselage to the spindle, wherein the workpiece motion chain contains the X-axis of the moving body, as C-axis of the moving body, and as the moving parts of the workpiece. On the other hand, the tool motion chain includes as the Z-axis of the moving body and as the moving body of the spindle for the cutting tool.

In the multibody system of gear grinding, it is also necessary to establish a fixed coordinate system for each moving body. The fixed coordinate system is based on the fixed body of the Earth as . The machine tool body coordinate system is taken as the reference coordinate system, and the subcoordinate system of each part is established for other moving bodies according to their own serial numbers in the topological structure. Through the multibody system theory, the kinematic equations of each axis of the machine tool can be established by using the relations among the coordinate systems.
3. Geometric Error Modeling of Cycloidal Gear Grinding Machine
In the process of cycloidal gear grinding, the machining accuracy of cycloidal gears will be directly affected by the motion axis of the machine tool. According to the above analysis of cycloidal wheel gear grinding machine in a total of three axes, respectively, there are two straight axis X and Z and a rotation axis C of the workpiece, so the geometric errors of machine tools include the level of each axis straightness error, vertical straightness error, linear positioning error, deflection angle error, back dip angle error, and the cross error. Then, there is the vertical error between the X-axis and the Z-axis in the Y direction and the vertical error between the X-axis and the Z-axis. Therefore, in the cycloidal gear grinding machine, there are three axes in the six degrees of freedom of 18 geometric errors and three-axis perpendicularity errors, a total of 21 geometric errors.
3.1. Establishing Coordinate System
The space position relation between the coordinate system of the C-axis and the workpiece coordinate system is determined by the transformation matrix . Since the workpiece is fixed on the C-axis turntable and rotates around its axial direction, the coordinate system of the workpiece moving body and the C-axis moving body can be regarded as the whole one body. There is no relative movement between the two bodies, so the position and displacement transformation matrix between the two can be obtained by the unit matrix:
When all the moving bodies are in theoretical motion, there is no motion error between them. Based on the movement of the X-axis, the rotation angle of the C-axis along the Z direction of the gear form grinding machine is expressed by . The theoretical coordinate transformation of rotating the moving body C-axis to the moving body X-axis is determined by the transformation matrix , so it can be expressed as
From above, we found that the second coordinate transformation matrix from the X-axis to the C-axis is obtained and is numerically represented as . The theoretical coordinate transformation from machine tool lathe bed coordinate system to X-axis coordinate system is determined by the transformation matrix . In the ideal state, the X-axis of the moving body moves in the X direction along the machine tool guide rail relative to the fuselage coordinate system. When the moving distance is x, the homogeneous coordinate transformation between the X-axis coordinate system and the fuselage coordinate system can be expressed as
Similarly, the homogeneous coordinate transformation matrix from the fuselage to the X-axis coordinate system is obtained and is numerically represented as . The theoretical coordinate transformation from Z-axis coordinate system to the fuselage coordinate system is determined by the transformation matrix . In the theoretical state, the Z-axis of the moving body moves along the Z direction of the machine tool guide rail relative to the fuselage coordinate system. When the moving distance is z, the homogeneous coordinate transformation between the Z-axis coordinate system and the fuselage coordinate system is expressed as a matrix . That is, it can be expressed as
The theoretical coordinate transformation from the grinding wheel coordinate system to the Z-axis coordinate system is determined by the transformation matrix . Since the grinding wheel is fixedly connected to the skateboard of the Z-axis, it can be regarded as the Z-axis coordinate system of the moving body that coincides with the coordinate system of the grinding wheel of the moving body. There is no relative motion between the two. Therefore, is also the identity matrix:
So far, the theoretical homogeneous coordinate system transformation matrix between each coordinate system of cycloid gear grinding machine gear grinding system has been obtained, and the second coordinate transformation matrix between the coordinate system of moving body grinding wheel and the coordinate system of moving body workpiece in the ideal state can be obtained by multiplying them as :
3.2. Homogeneous Coordinate Transformation with Geometric Errors
The coordinate transformation from X-axis coordinate system to C-axis coordinate system with geometric error is expressed as . In the actual machining process, the actual coordinate transformation from C-axis to X-axis coordinate system is equivalent to the superposition of an error matrix of the C-axis relative to the X-axis coordinate system. On the homogeneous coordinate transformation matrix from the coordinate system to the coordinate system under the state of X-axis movement ,
is the error matrix generated when the moving body moves from the C-axis to the X-axis coordinate system, including the angular displacement error around the Z-axis , straightness error along X-axis, Y-axis, and Z-axis (, , and ), and rotation angle error along Y-axis and X-axis ( and ). So we get the transformation matrix as follows:
The actual coordinate transformation from C-axis to X-axis is obtained by inserting the error matrix into the formula:
Similarly, coordinate transformation from X-axis coordinate system to C-axis coordinate system with geometric error is described as , which is the inverse matrix , so it can be numerically represented as follows:
The actual coordinate system transformation from fuselage coordinate system to X-axis coordinate system with geometric error is expressed as . In the actual processing machine operation, is expressed in terms of the geometric error term order rotation matrix, so
in terms of the geometric error is described by the term order displacement vector, so
Since the X-axis of the moving body moves x distance along the X-axis direction of its coordinate system, then the above vector is expressed as . The motion error of the X-axis of the moving body is described in terms of in the X, Y, and Z directions of the machine tool coordinate system, respectively. The perpendicularity error of the X-axis of a moving body is expressed as in the X, Y, and Z directions.
In the actual state, the coordinate transformation of the geometric error term coordinate system to the fuselage coordinate system can be represented as
The ideal coordinate transformation including the geometric error term X-axis to the fuselage coordinate system can be obtained by substituting the matrix . This includes the linear displacement error along the X-axis, the straightness error and , the yaw angle error and in the X-Y plane and the X-Z plane, and the rotation error on the X-axis, so
Similarly, coordinate transformation from the machine tool coordinate system to the X-axis coordinate system with geometric error is expressed as , that is to say, it is the inverse of . Therefore, it can be represented by an inverse matrix as follows:
The coordinate transformation from the coordinate system to the fuselage coordinate system with geometric error is expressed as . Similarly, when the moving body moves a distance z along the Z direction with the grinding wheel, the homogeneous coordinate transformation matrix from the Z-axis coordinate system to the machine tool coordinate system can be obtained. It includes the linear displacement error along with the z-axis term, straightness error and in the Y-Z plane and X-Z plane, yaw angle error and in the Y-Z plane and X-Z plane, and Z-axis transverse rotation error and can be expressed as
Since there is no relative movement among the workpiece, the turntable, the grinding wheel, and the Z-axis, it can be regarded as the workpiece coordinate system coinciding with the C-axis coordinate system , the grinding wheel coordinate system , and the coordinate system . Therefore, there is no relative position error and rotation angle error between matrix and . Then, it can also be viewed as the identity matrix, which is
So far, the coordinate transformation matrix from the moving body grinding wheel coordinate system containing geometric errors to the moving body workpiece coordinate system can be obtained by multiplying them , so it is represented as
3.3. Geometric Error Model of Gear Grinding Machine
The coordinate transformation matrix from the grinding wheel coordinate system to the workpiece coordinate system , which contains the geometric error, is expressed as . It can also be viewed as the homogeneous coordinate transformation matrix between these two coordinates in the ideal state . An actual matrix containing its rotation angle error and relative position error is superimposed on the grinding wheel coordinate system to the workpiece coordinate system , which is expressed as . So the geometric error can be obtained by the equation as follows:
As mentioned above, when small quantities of second order and above are ignored, the hypothesis principle of minimum error can be used as the basis and is represented aswhere , , and are errors of relative displacement form the grinding wheel coordinate system to the workpiece coordinate system along the X, Y, and Z directions, respectively; , , and are errors of rotation angle between the grinding wheel coordinate system and the workpiece coordinate system along the X direction, Y direction, and Z direction, respectively.
Taking into account the difficulty of calculating errors of relative position and rotation angle, we can use MATLAB software to calculate each matrix between the grinding wheel coordinate system and the workpiece coordinate system in X, Y, and Z directions. The gear form grinding machines involve two processes, namely, gear grinding and wheel dressing, so they must be equipped with four axes to satisfy the demands of these two processes, respectively. So the geometric error model of the cycloid CNC gear grinding machine can be obtained as follows:
4. Simulation of Geometric Error
4.1. Analysis of the Influence of Geometric Error
To reflect the influence of the geometric errors of each axis of the cycloid gear grinding machine, it is also necessary to find the relationship between the geometric errors of the machine tool and the tooth surface errors after the completion of the geometric error model. First, according to the parameters of a cycloidal gear in Table 1, the grid spatial coordinate data of cycloidal gear tooth surface are calculated, and the calculated coordinate data are imported into MATLAB for tooth surface simulation.
Using the control variable method under the condition that the geometric errors of each axis of the machine tool enhance one other and cancel the other ones, it is assumed that only the displacement error of each axis is 0.005 mm or the rotation error 0.01°. Finally, the errors of other directions are zero, and the tooth surface of the error simulation is shown in Figure 3.

(a)

(b)

(c)

(d)

(e)

(f)
Above all, we can intuitively see the influence of the error components of the geometric errors of the machine tool on the tooth surface in the six degrees of freedom between the grinding wheel and the workpiece coordinate system. Through the simulation of the tooth surface, we can compensate for the geometric errors of each axis of the machine tool by measuring the parameters of the processed cycloidal gear.
4.2. Relationship between Geometric Errors and Tooth Surface Errors
According to the YK7350 B cycloidal gear grinding machine discussed in this study, the coordinate analysis of the spatial position relationship between the grinding wheel and workpiece is shown in Figure 4.

Coordinate systems and , whose relative positions are described by the machine tool coordinate systems, are rigidly connected to the grinding wheel and work gear, respectively. When the coordinate system of the two moving bodies is without error, the actual processed cycloidal gear coincides with the theoretical cycloidal gear at the same time. Therefore, the gear cross-section z0 is coincident with . The actual tooth surface equation is the profile equation of the grinding wheel in the coordinate system of the machine tool, which can be expressed aswhere and are the radius of the needle wheel and the needle gear sleeve, respectively; and are the number of work gear teeth and needle wheel teeth, respectively; is the short amplitude coefficient; and is the rotation angle of the workpiece.
By rotating the coordinate of the cycloid equation of the first tooth surface, the cycloid equation of the number tooth surface can be determined by
In the processing of cycloidal gears form grinding machine, the movement of each axis determines the trajectory of the grinding wheel. If the dressing error of the grinding wheel is also ignored, then the coordinate vector of the dressed axial wheel profile can be expressed aswhere is the axial parameter of the grinding wheel surface.
The rotary surface of the grinding wheel is the track surface swept by the axial profile of the grinding wheel when it rotates around its axis. In the coordinate system of grinding wheel, the rotation parameter of grinding wheel is expressed by , and then, the coordinate vector of grinding wheel can be expressed as
Among them,
Therefore, the unit normal vector of the wheel surface can be expressed as
In the ideal case, the equation of grinding wheel surface is jointly expressed by the unit normal vector and coordinate vector . According to the homogeneous coordinate transformation matrix between the cycloidal gear coordinate system and the grinding wheel coordinate system in the ideal state, the grinding wheel surface in the workpiece coordinate system can be expressed as
In the actual machining process, the grinding wheel determines the position relationship between its coordinate system and the workpiece coordinate system . The geometric error of each axis of the machine tool also determines the position of the grinding wheel, and all the influence of the geometric error of each axis of the machine tool is used ; then, the actual grinding wheel surface equation can be expressed as
When the parameters of the grinding wheel are known, the grinding contact point is defined as a radial vector from the origin of the workpiece coordinate system to the surface of the grinding wheel. Therefore, the contact conditions in the ideal state and the actual condition are, respectively, as follows:
Here, the machining parameter is represented by constant , and contact conditions between the grinding wheel and workpiece will lead to . While these conditions are true, there remains the only parameter concerning the grinding wheel rotation profile and grinding wheel axial profile parameters in the equation. Here, the range of is known, that is, let . According to the contact conditions between the grinding wheel and the workpiece, we can figure out the corresponding angle . So the coordinate vector of the contact point is obtained, and the unit normal vector is calculated out, too. Because it is through the profile parameters and contact conditions to find out , can be represented as . Therefore, the contact line between the grinding wheel and the workpiece at the end of the k can be expressed as:
In the process of gear processing, the contact line between the grinding wheel and the workpiece together is the tooth surface processed by the machine tool, and the theoretical tooth surface and actual tooth surface composed of discrete contact lines can be, respectively, expressed as
In summary, the relationship between geometric errors of each axis and tooth surface errors of the CNC gear grinder can be established, and its model can be expressed as follows:where and .
4.3. Simulation Analysis of Error Sensitivity
Morris method is used to calculate the mean of the basic effects for each input parameter and the standard deviation as a sensitivity index.
The mean value and standard deviation of the basic effect of each error term of geometric errors can be obtained, and the obtained mean value and standard deviation are normalized as
Finding the normalized mean and the standard deviation can be used as a sensitivity index to identify the sensitivity of geometric errors. Here, the normalized mean value the larger the error, the greater the influence of the error on the tooth surface error. And the standard deviation the larger the error, the greater the coupling effect between the error and other error terms. The sensitivity index of machine tool geometric error on tooth surface error is shown in Figure 5.

(a)

(b)

(c)

(d)

(e)

(f)
5. Geometric Error Measurement and Compensation
A laser interferometer is a length and angle measurement instrument. Using the laser wavelength as the standard, two beams of coherent light in the process of transmission will form interference and get the light and dark stripes on the laser. The wavelength of the laser is relatively stable, and the distance between each fringe is one wavelength of the laser. Through these stripes, the laser head can obtain the moving distance of the measured object.
According to the formula,
Due to the and the , , so
The measured distance l is obtained and numerically represented aswhere N is the number of the cumulative pulse; is the laser wavelength; c is the speed of light; and v is the moving speed of the moving mirror.
In actual machine tool error measurement work, a Renishaw dual-frequency laser interferometer is used. In addition to the principle of light interference, it also uses the frequency change generated by the Doppler effect to detect the position of each axis of the machine tool. After the laser head is electrified, a beam of light is divided into two beams of coherent light through the spectroscope. The polarization plane is perpendicular to the incident plane. A beam of light is reflected through the spectroscope and radiated to the fixed mirror to form F1, which is the reference light. Another beam of light in the incident plane is transmitted through the spectroscope to the moving mirror to form F2, which is the measurement light. The actual measurement is shown in Figure 6.

5.1. Measuring Principle and Measurement of Rotating Axis
The principle of measuring the error of the rotation axis of the CNC machine tool is the same as that of measuring the error of the linear axis. It only needs to change that the linear length change of the moving mirror is changed to the angle change of the laser angle mirror by using the trigonometric function. The Renishaw XR20-W nonlinear rotary shaft calibration device is used in the actual measurement. The high-precision angle mirror on the device also has a high-precision small motor. The XR20-W is mounted at the center of the machine tool rotary shaft turntable, and the angle mirror small turntable can be precisely rotated through the internal high-precision grating system and drive.
Similarly, using the frequency difference generated by light interference and the Doppler effect, the displacement difference of the mirror can be obtained. According to the optical path principle of measurement, the relative optical path change of measuring arm 2 and measuring arm 1 is . If the nominal spacing s between two mirrors in an angular mirror is yes, then the measuring system software can calculate the angle error by arcsine as follows:where s is the parameter of the nominal spacing between two mirrors.
5.2. Measurement Geometric Error of Each Axis
To investigate the response of each axis to the geometric error that is independently changed and the corresponding geometric error that is obtained, various error data of the machine tool are obtained through the actual measurement, which provide a data basis for the subsequent geometric error compensation of the CNC machine tool. Figure 7 shows the geometric error of actual measurement by a laser interferometer.

(a)

(b)

(c)

(d)
5.3. Compensation Result of Geometric Error
From the measurement analysis curve of the C-axis after compensation, we can see that the angle error and repeated angle error of the machine tool are significantly reduced (as shown in Figure 8), and the performance of the machine tool is significantly improved. By measuring the geometric error data of each axis of the gear grinding machine and analyzing the compensation values of the error compensation points of each motion axis, a machine tool error compensation table is established. Using software compensation, the error compensation value of each compensation point is written into the program through Siemens 840D system to compensate for the error of each axis of the machine tool. After the error compensation, the machine tool is run again and the error value of each motion axis of the machine tool is measured again. Geometric errors are significantly reduced.

(a)

(b)

(c)

(d)
6. Experimental
6.1. Cycloidal Gear Grinding
The cycloidal gears processed by the gear grinder are experimentally ground, and the cycloidal gear processing by the CNC gear grinder before and after the geometric error compensation is carried out. Grinding is performed on the gear form grinding according to the parameters of the cycloidal gear (in Table 2). The actual grinding process is shown in Figure 9.

6.2. Detection and Measurement
The installation and measurement processes of each part in the actual measurement of cycloidal gear are shown in Figure 10.

The measuring equipment uses the Gleason 650GMS gear detection analyzer, which can measure gears up to 650 mm in diameter, and the measuring speed is increased by 25%. The gear detection analyzer is also equipped with a Renishaw SP80H 3D scanning probe, which can carry out combined measurement of measuring needles of various sizes and extended lengths, with high speed and high-precision measurement capability. The measurement results are shown in Figures 11 and 12.


The test results are shown in Table 2, indicating that the machining accuracy of the cycloid gear grinding machine is significantly improved.
7. Conclusions
In this study, the geometric error model of the CNC gear grinding machine is established by a multibody motion system. The theoretical tooth surface and error tooth surface are simulated, and the compensation of the axial geometric error method is used to research the sensitivity of the influence of geometric error of each axis on the tooth surface. Then, the geometric error of each axis of the machine tool is measured and compensated. The accuracy of the machine tool is verified through the gear processing experiment, which provided a reliable guarantee for the manufacture of high-precision cycloid gear. As the numerical examples show, compensation of the axial geometric error method can improve tooth quality. An additional advantage of the approach is to provide the user with valuable information about the grinding machine and machining capacity.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors are grateful to the Natural Science Foundation of China Youth Fund and the Postdoctoral Science Foundation of China for their financial support. Part of this work was performed under contract no. 52005157 and no. 2021M690051.