Research on CNC System for Large Cylindrical Gear Milling Machining

The development of high-performance CNC systems for large cylindrical gear milling is critical for industries like wind power, mining, and shipbuilding, where traditional hobbing solutions are cost-prohibitive for small-to-medium enterprises. We present an economical PC-based CNC system leveraging motion control cards and digital gear modeling to enable precision milling of large gears (modules >12mm, diameters >250mm). Our approach integrates parametric gear profile generation, real-time interpolation, and modular software architecture to achieve flexible and cost-effective manufacturing.

Gear milling requires precise involute curve generation. The involute profile is mathematically defined by the base circle radius $r_b$ and pressure angle $\alpha_k$ at any point $K$:

$$
r_k = \frac{r_b}{\cos\alpha_k}, \quad \theta_k = \tan\alpha_k – \alpha_k = \text{inv}\alpha_k
$$

In Cartesian coordinates relative to gear center $O$:

$$
\begin{cases}
x_k = r_b(\cos\theta_k + \theta_k \sin\theta_k) \\
y_k = r_b(\sin\theta_k – \theta_k \cos\theta_k)
\end{cases}
$$

Transformation to CNC machining coordinates $X’O’Y’$ involves rotation by $\omega_b$, the base circle half-thickness angle:

$$
\omega_b = \frac{180^\circ}{\pi} \left( \frac{\pi}{2z} + \text{inv}\alpha \right)
$$

Resulting in transformed coordinates:

$$
\begin{cases}
x’_k = \sqrt{x_k^2 + y_k^2} \sin\left( \tan^{-1}\frac{y_k}{x_k} + \omega_b \right) \\
y’_k = \sqrt{x_k^2 + y_k^2} \cos\left( \tan^{-1}\frac{y_k}{x_k} + \omega_b \right)
\end{cases}
$$

Standard Gear Module Series (GB/T1357-1987)
Series Module (mm)
I 0.12, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5, 0.6, 0.8, 1, 1.25
II 2, 2.5, 3, 4, 5, 6, 8, 10, 12, 16, 20, 25, 32
III 0.35, 0.7, 0.9, 1.75, 2.25, 2.75, (3.25), 3.5, (3.75), 4.5, 5.5, (6.5), 7, 9, (11), 14, 18, 22, 28, (30), 36

G60 Gear Milling Machine

Our CNC hardware architecture uses a GT-400-SV motion controller with PCI bus communication to an IPC-810 industrial computer. Servo control employs PID+Kvff+Kaff algorithms:

$$
U_n = K_p E_n + K_i \sum_{i=1}^n E_i + K_d (E_n – E_{n-1}) + K_{vff} (P_{\text{target}} – P_{\text{actual}}) + K_{aff} \frac{d^2P}{dt^2}
$$

where $E_n$ is position error, and $P_{\text{target}}$, $P_{\text{actual}}$ are target/actual positions.

G/M Code Definitions
Group Code Function
I G00/G01 Rapid traverse/Linear interpolation
G02/G03 Clockwise/Counter-clockwise arc
G90/G91 Absolute/Incremental coordinates
II G17/G18/G19 XY/XZ/YZ plane selection
III G41/G42 Left/Right cutter compensation

Software modules were developed in C++ using MFC classes:

  1. CCodeParser: NC program lexical/syntactic analysis
  2. CNCDriver: Motion control (GotoXYZ(), GotoARC())
  3. GEAR: Involute interpolation point calculation

Interpolation points $n$ for involute segment AB satisfy surface tolerance constraints:

$$
| x'(i+1) – x'(i) | + f_{pt} \leq 8 \mu m, \quad | y'(i+1) – y'(i) | + f_{pt} \leq 8 \mu m
$$

where $f_{pt}$ is pitch tolerance (e.g., 0.063mm for Grade 8 gears). For a module 25, 40-tooth gear, key points were calculated:

Key Gear Profile Coordinates (mm)
Point A B C D E
X 9.5078 25.4411 43.5595 29.9881 43.5595
Y 524.9139 469.1570 539.000 467.7898 466.7317

Testing confirmed 0.005mm axis repeatability and successful simulation of full gear milling cycles. The system enables high-accuracy gear milling at 60% lower cost than traditional hobbing for small-batch production.

Scroll to Top