In modern mechanical transmission systems, gears play a pivotal role, and with the advancement of technology, high-end equipment demands increasingly higher precision in gear manufacturing. Tooth profile modification, particularly in gear grinding processes, is essential for improving load distribution and enhancing the load-bearing capacity of gears. However, during the CNC form grinding of internal helical gears, multi-axis additional motions introduce machining errors, leading to phenomena such as tooth surface distortion and even grinding cracks. This study focuses on analyzing the mechanisms behind tooth profile errors in internal helical gears during form grinding and proposes effective compensation methods to enhance machining accuracy. We establish a comprehensive error model based on the kinematic relationships in CNC gear profile grinding machines and develop a software-based compensation approach to mitigate these errors. Through theoretical analysis and experimental validation, we demonstrate the effectiveness of our method in achieving higher precision in tooth profile modification.
The foundation of this research lies in understanding the spatial engagement between the grinding wheel and the internal helical gear during the gear grinding process. We begin by setting up a coordinate system for the CNC form grinding machine, which includes three linear motion axes (x, y, z) and three rotational motion axes (A, B, C). The grinding wheel is mounted on a wheel arm that rotates about the A-axis, while the workpiece is installed on the B-axis, allowing for helical motion through the interpolation of B and z axes. The spatial engagement relationship is derived using homogeneous transformation matrices. Let $s_1$ and $s_t$ represent the coordinate systems fixed to the gear and grinding wheel, respectively. The position vector of the tooth surface in $s_1$ is given by $\mathbf{r}_1(u, \theta) = [x_1, y_1, z_1, 1]^T$, where $u$ and $\theta$ are surface parameters. The transformation to the grinding wheel coordinate system $s_t$ involves rotation and translation operations:
$$\mathbf{r}_t(u, \theta, \phi_1) = \mathbf{M}_{tc} \mathbf{M}_{ca} \mathbf{M}_{a1} \mathbf{r}_1(u, \theta)$$
Here, $\phi_1$ is the rotation angle of the workpiece, $\mathbf{M}_{a1}$ is the rotation matrix about the $z_1$-axis, $\mathbf{M}_{ca}$ accounts for the center distance $E_{tp}^{(0)}$ and axial movement $L_t$, and $\mathbf{M}_{tc}$ represents the grinding wheel installation angle $\gamma_m$. The contact condition between the grinding wheel and the tooth surface is determined by the conjugate surface theory, where the vector $\mathbf{r}_t$, unit normal vector $\mathbf{n}_t$, and grinding wheel axis $\mathbf{y}_t$ are coplanar. The contact equation is expressed as:
$$f_t = \mathbf{n}_t(u_1, \theta_1, \phi_1) \cdot [\mathbf{y}_t \times \mathbf{r}_t(u_1, \theta_1, \phi_1)] = 0$$
The unit normal vector $\mathbf{n}_t$ is derived from the partial derivatives of $\mathbf{r}_t$ with respect to $u_1$ and $\theta_1$. The contact line, which is a spatial curve, is projected onto the grinding wheel’s axial section to obtain the wheel profile. The coordinates in the axial plane are given by:
$$x_w(u_1, \theta_1, \phi_1) = \sqrt{x_t^2(u_1, \theta_1, \phi_1) + z_t^2(u_1, \theta_1, \phi_1)}$$
$$y_w(u_1, \theta_1, \phi_1) = y_t(u_1, \theta_1, \phi_1)$$
To achieve tooth profile modification, additional radial motion (x-axis) and rotational motion (B-axis) are superimposed on the grinding wheel’s path. The modification amount $\Delta E$ follows a parabolic form along the tooth width $l$, with the gear rotation angle $\theta$ varying based on workpiece parameters. The radial feed amount is defined as:
$$\Delta E = \begin{cases}
\frac{a_{ml} l_p}{2} (\theta – \theta_b)^2 & \theta_a \leq \theta \leq \theta_b \\
0 & \theta_b \leq \theta \leq \theta_c \\
\frac{a_{ml} l_p}{2} (\theta – \theta_c)^2 & \theta_c \leq \theta \leq \theta_d
\end{cases}$$
where $a_{ml}$ is the tooth profile modification coefficient, and $\theta_i$ (i = a, b, c, d) represents the gear rotation angles at specific points along the tooth profile.
The multi-axis additional motions during gear grinding introduce errors that affect tooth surface accuracy. Specifically, the x-axis additional motion causes a radial displacement $\Delta x$, leading to varying normal offsets at different points on the involute profile due to changing pressure angles. This results in tooth profile slope deviations. The deviations at the tooth tip and root are calculated as:
$$\Delta x_a = \Delta x \sin \lambda_a$$
$$\Delta x_f = \Delta x \sin \lambda_f$$
$$f_H = \Delta x_f – \Delta x_a$$
where $\lambda_a$ and $\lambda_f$ are the pressure angles at the tip and root, respectively. The B-axis additional rotation $\Delta b$ causes a phase shift in the tooth profile, resulting in deviations proportional to the radii:
$$\Delta b_a = r_a \Delta b$$
$$\Delta b_f = r_f \Delta b$$
$$f_H = \Delta b_f – \Delta b_a$$
The helical motion in gear grinding is achieved through the interpolation of z-axis movement and workpiece rotation, governed by the relationship $C_B = Z_B \times \tan \beta / r$, where $Z_B$ is the axial movement distance, $\beta$ is the helix angle, and $r$ is the pitch radius. The grinding wheel installation angle $\gamma_m$ also influences the contact line morphology. By adjusting $\gamma_m$ within a small range (e.g., ±1° to 2°), the contact line length and curvature can be optimized, thereby reducing errors in tooth profile modification. For instance, decreasing the installation angle shortens the contact line and increases its deviation on the left and right tooth surfaces, while increasing it lengthens and curves the contact line more prominently.
To quantify and compensate for these errors, we develop a tooth surface error model based on the actual grinding process. The model compares the theoretical tooth surface with the actual ground surface, identifying distortions such as over-modification or under-modification on one side of the tooth. For internal helical gears, the left and right tooth surfaces are ground simultaneously, so compensation must balance errors on both sides. We propose a method that neutralizes the errors by adjusting the multi-axis motions. The adjustment amounts for the helix angle $\beta_g$ and the swing angle $C_{Bg}$ are derived from the tooth profile slope deviations $f_{H\beta L}$ and $f_{H\beta R}$ measured over the effective tooth length $D_1$:
$$\beta_g = \arctan \left( \frac{(f_{H\beta L} – f_{H\beta R}) / 2}{D_1} \right)$$
$$C_{Bg} = \frac{D_1 \times \tan \beta_g}{r}$$
Based on this, we create a software tool for tooth profile accuracy compensation, which calculates the necessary adjustments to the machine axes. The software inputs include measurement positions along the tooth profile and the slope deviations for left and right surfaces, outputting optimized grinding parameters.
In our experimental validation, we use a Siemens-controlled YK7350 CNC grinding machine to process an internal helical gear with parameters as listed in Table 1. The gear undergoes initial grinding, followed by precision measurement using a Gleason 650GMS inspection center. The results, summarized in Table 2, show that before compensation, the left tooth surface has an average total profile deviation $F_\beta$ of 8.6 μm and slope deviation $f_{H\beta}$ of -2.3 μm, corresponding to grade 7 accuracy, while the right surface has $F_\beta$ of 5.4 μm and $f_{H\beta}$ of -4.7 μm, at grade 6. After applying compensation through our software, the left surface improves to $F_\beta$ of 4.3 μm and $f_{H\beta}$ of -2.6 μm, and the right surface to $F_\beta$ of 10.0 μm and $f_{H\beta}$ of -6.0 μm, both achieving grade 6 accuracy. This demonstrates the effectiveness of our error compensation method in enhancing tooth profile precision during gear profile grinding.
Throughout the gear grinding process, issues such as grinding cracks can arise due to excessive thermal loads or improper wheel-workpiece engagement. These cracks are often initiated by high grinding forces and inadequate cooling, leading to micro-fractures on the tooth surface. In severe cases, grinding cracks can compromise the gear’s fatigue life and operational reliability. Therefore, controlling the grinding parameters and implementing error compensation are crucial not only for accuracy but also for preventing defects like grinding cracks. The following image illustrates typical grinding cracks observed in gear manufacturing, highlighting the importance of optimized grinding strategies:

To further analyze the impact of multi-axis motions, we conduct numerical simulations using MATLAB, modeling the theoretical tooth surface and grinding contact lines. The simulation results, presented in Table 3, show how variations in x-axis and B-axis motions affect the tooth profile deviations. For example, a 10 μm additional x-axis motion can lead to a slope deviation of up to 3 μm, while a 0.01° B-axis rotation causes deviations proportional to the tooth radius. These findings underscore the need for precise control of machine axes during gear profile grinding.
In conclusion, our research establishes a comprehensive framework for error compensation in tooth profile modification of internal helical gears using CNC form grinding. By modeling the spatial engagement, analyzing error sources, and developing a practical compensation method, we achieve significant improvements in gear accuracy. The experimental results validate our approach, showing that balanced compensation of left and right tooth surfaces can elevate the gear grinding quality to higher precision grades. Future work will focus on integrating real-time monitoring and adaptive control to further reduce errors and prevent grinding cracks, advancing the capabilities of gear manufacturing in high-performance applications.
| Parameter | Value |
|---|---|
| Normal module $m_n$ (mm) | 2 |
| Number of teeth $z$ | 79 |
| Normal pressure angle $\alpha_n$ (°) | 20 |
| Helix angle $\beta$ (°) | 15 |
| Hand of spiral | Right-hand |
| Face width $d$ (mm) | 65 |
| Normal shift coefficient $x_n$ | 0.4987 |
| Tooth profile modification amount (μm) | 5 ± 4 |
| Evaluation start position (mm) | 4.5 |
| Evaluation end position (mm) | 40.5 |
| Tooth Surface | Tooth Number | $F_\beta$ (μm) Before | $f_{H\beta}$ (μm) Before | $F_\beta$ (μm) After | $f_{H\beta}$ (μm) After |
|---|---|---|---|---|---|
| Left | 1 | 6.8 | 6.9 | 5.1 | -4.6 |
| 20 | 10.0 | -10.2 | 6.2 | -5.2 | |
| 40 | 6.2 | 5.4 | 2.9 | 0.5 | |
| 60 | 11.5 | -11.4 | 3.1 | -1.0 | |
| Right | 1 | 6.6 | -6.2 | 5.4 | -2.4 |
| 20 | 4.0 | -3.7 | 9.0 | -4.4 | |
| 40 | 8.0 | -8.3 | 14.0 | -9.6 | |
| 60 | 2.9 | -0.4 | 11.5 | -7.4 |
| Axis Motion | Motion Amount | Slope Deviation $f_H$ (μm) | Total Deviation $F_\beta$ (μm) |
|---|---|---|---|
| X-axis | 5 μm | 1.5 | 3.2 |
| 10 μm | 3.0 | 6.5 | |
| 15 μm | 4.5 | 9.8 | |
| B-axis | 0.005° | 2.1 | 4.3 |
| 0.01° | 4.2 | 8.6 | |
| 0.015° | 6.3 | 12.9 |
The gear grinding process involves complex interactions between the grinding wheel and workpiece, where parameters such as wheel speed, feed rate, and cooling conditions must be optimized to minimize errors and defects. In gear profile grinding, the form of the grinding wheel is critical, as it directly imprints the desired profile onto the tooth surface. Any deviations in the wheel profile or machine kinematics can lead to inaccuracies, emphasizing the need for precise modeling and compensation. Our approach not only addresses geometric errors but also contributes to reducing the risk of grinding cracks by ensuring smoother engagement and controlled material removal. As the demand for high-precision gears grows in industries like aerospace and automotive, advancements in gear grinding techniques will continue to play a vital role in meeting these challenges.
