Outsize straight bevel gears are indispensable components in heavy-duty transmission systems, particularly in industries such as power generation, construction materials, marine engineering, and mining machinery. In my research, I focus on a new envelope shaping method for outsize straight bevel gear tooth surfaces, aiming to overcome the limitations of existing machining technologies. The proposed method is based on the geometric properties of the theoretical tooth surface of a straight bevel gear. I derive the tooth surface model, analyze its convexity, establish the machining principle, develop the machine motion model, control the envelope accuracy, and implement a dedicated calculation software system. Through numerical simulation and calculation examples, I verify the correctness and feasibility of the proposed method. This research is intended to provide a theoretical foundation for the design of CNC planning machines for outsize straight bevel gears and to promote the progress of large gear manufacturing technology.
1. Introduction
A straight bevel gear is used to transmit motion and power between intersecting axes. Compared with spiral bevel gears and hypoid gears, the straight bevel gear is easier to design, manufacture, and adjust. Therefore, it remains one of the most widely used gear forms in modern mechanical transmission systems. The outsize straight bevel gear, usually defined as a gear with a large-end diameter larger than 3000 mm, is critical in large-scale industrial equipment. However, the machining of outsize straight bevel gears is much more difficult than that of small and medium-sized gears because of the huge workpiece size, high inertia, and strict requirements for tooth surface accuracy.
Traditional machining methods for straight bevel gears are mainly based on the generating principle with a virtual crown gear. In such methods, the workpiece and the virtual generating gear rotate in a constrained relationship, and the cutting edges of reciprocating planing tools generate the tooth flank through enveloping motions. The generating method works well for small and medium straight bevel gears, but it is not economical for outsize straight bevel gears. Since the generating gear must be much larger than the workpiece, the machine tool becomes extremely large, expensive, and inefficient in terms of floor space. Some manufacturers try to machine large straight bevel gears on modified shapers or milling machines using form cutters. These approaches often suffer from low accuracy, poor efficiency, poor cutter universality, and high cost. In addition, traditional calculation methods based on the equivalent spur gear at the large end inevitably introduce theoretical errors when computing tooth thickness and cutter settings.
In my work, I address these problems by proposing a new envelope shaping method for outsize straight bevel gear tooth surfaces. Instead of using a virtual generating gear, the machining principle is derived directly from the geometric property that the theoretical tooth surface of a straight bevel gear is a convex ruled surface. Multiple planes are generated by the straight cutting edge of a planing tool, and these planes envelope the theoretical tooth surface to form a polyhedral surface. By controlling the number of enveloping planes, the machined surface can approximate the theoretical tooth surface within a specified accuracy. This method avoids the traditional equivalent gear calculation, reduces theoretical errors, and allows the workpiece to rotate continuously during machining. Consequently, the processing efficiency is significantly improved and the requirement for the workpiece driving motor is reduced.
2. Theoretical Tooth Surface of Straight Bevel Gear
To develop a reliable machining method, I first establish the mathematical model of the theoretical tooth surface of a straight bevel gear. The generation principle is described as follows. A circular plane, called the generating plane, rolls without sliding on a base cone. The tangent line between the generating plane and the base cone sweeps a spatial surface, namely the involute cone surface. The intersection of this surface with a sphere centered at the cone apex is a spherical involute. The whole tooth surface of a straight bevel gear can be considered as a family of spherical involutes corresponding to different spherical radii.
From the viewpoint of differential geometry, the theoretical tooth surface of a straight bevel gear is a ruled surface. I therefore express the tooth surface in the following form:
$$ \mathbf{S}(r,\phi) = (1-r)\mathbf{W}(\phi) + r\mathbf{Q}(\phi), \quad 0\le r \le 1, $$
where \(\mathbf{Q}(\phi)\) represents the spherical involute at the large end, \(\mathbf{W}(\phi)\) represents the spherical involute at the small end, \(r\) is the parameter along the straight generatrix, and \(\phi\) is the parameter characterizing the rolling angle of the generating plane.
3. Mathematical Model of the Large-End and Small-End Tooth Profiles
In order to obtain the analytical expressions of \(\mathbf{Q}(\phi)\) and \(\mathbf{W}(\phi)\), I introduce two coordinate systems. The fixed coordinate system \(S(x,y,z)\) is attached to the workpiece, with the origin at the base cone apex, the \(z\)-axis coinciding with the base cone axis, and the \(x\)-axis parallel to the radius corresponding to the starting point of the spherical involute. An auxiliary moving coordinate system \(S_1(x_1,y_1,z_1)\) is attached to the generating plane. At any instant, the \(z_1\)-axis coincides with the instantaneous contact line between the generating plane and the base cone, the \(x_1\)-axis lies in the generating plane, and the \(y_1\)-axis is determined by the right-hand rule.
For a point \(K\) on the spherical involute, its coordinates in the moving coordinate system are given by:
$$ \begin{bmatrix} x_1 \\ y_1 \\ z_1 \end{bmatrix} = \begin{bmatrix} R\sin\psi \\ 0 \\ R\cos\psi \end{bmatrix}. $$
Because the generating plane rolls without sliding, the pure rolling condition yields:
$$ \psi = \phi \sin\delta_b, $$
where \(R\) is the cone distance, \(\delta_b\) is the base cone angle, and \(\phi\) is the angle between the initial radius and the instantaneous contact radius on the base cone bottom plane. Using the coordinate transformation from \(S_1\) to \(S\), I obtain the large-end tooth profile expression:
$$
\mathbf{Q}(\phi) =
\begin{bmatrix}
Q_x(\phi) \\
Q_y(\phi) \\
Q_z(\phi)
\end{bmatrix}
=
\begin{bmatrix}
R\left[\cos(\phi\sin\delta_b)\sin\delta_b\cos\phi+\sin(\phi\sin\delta_b)\sin\phi\right] \\
R\left[\cos(\phi\sin\delta_b)\sin\delta_b\sin\phi-\sin(\phi\sin\delta_b)\cos\phi\right] \\
R\cos(\phi\sin\delta_b)\cos\delta_b
\end{bmatrix}.
$$
Similarly, the small-end tooth profile is obtained by replacing \(R\) with \(R-B\), where \(B\) is the face width:
$$
\mathbf{W}(\phi) =
\begin{bmatrix}
(R-B)\left[\cos(\phi\sin\delta_b)\sin\delta_b\cos\phi+\sin(\phi\sin\delta_b)\sin\phi\right] \\
(R-B)\left[\cos(\phi\sin\delta_b)\sin\delta_b\sin\phi-\sin(\phi\sin\delta_b)\cos\phi\right] \\
(R-B)\cos(\phi\sin\delta_b)\cos\delta_b
\end{bmatrix}.
$$
Therefore, the theoretical tooth surface of the straight bevel gear is completely determined by the two base curves \(\mathbf{Q}(\phi)\), \(\mathbf{W}(\phi)\), and the straight lines connecting the corresponding points with the same parameter \(\phi\). This ruled-surface model is the foundation of the subsequent machining principle and motion analysis.
4. Geometric Characteristics of the Tooth Surface
In my research, the geometric convexity of the theoretical tooth surface of a straight bevel gear is essential for the envelope machining method. By differentiating the tooth surface model with respect to the parameters \(r\) and \(\phi\), I obtain:
$$ \mathbf{S}_r = \frac{\partial \mathbf{S}(r,\phi)}{\partial r} = \mathbf{Q}(\phi)-\mathbf{W}(\phi), $$
$$ \mathbf{S}_\phi = \frac{\partial \mathbf{S}(r,\phi)}{\partial \phi} = (1-r)\mathbf{W}'(\phi)+r\mathbf{Q}'(\phi). $$
Using the relationship between the large-end and small-end tooth profiles,
$$ \mathbf{W}(\phi)=\frac{R-B}{R}\mathbf{Q}(\phi), $$
the \(\phi\)-direction tangent vector can be simplified as:
$$ \mathbf{S}_\phi = \left[(1-r)\frac{R-B}{R}+r\right]\mathbf{Q}'(\phi). $$
This expression shows that, for a fixed \(\phi\), the tangent vector \(\mathbf{S}_\phi\) has the same direction at all points along the straight generatrix. In other words, the tangent plane of the tooth surface at every point on the same cone generatrix is identical. Consequently, the convexity of the tooth surface is determined by the change of \(\mathbf{Q}'(\phi)\).
I project \(\mathbf{S}_\phi\) onto the \(xOy\) plane of the fixed coordinate system and define the angle \(\eta\) between the projection and the \(x\)-axis. The cosine of this angle is:
$$ \cos\eta(\phi)=\frac{Q_x'(\phi)}{\sqrt{Q_x’^2(\phi)+Q_y’^2(\phi)}}. $$
Differentiating with respect to \(\phi\), I obtain:
$$ \frac{d\cos\eta}{d\phi}=
\frac{Q_y'(\phi)\left[Q_y'(\phi)Q_x”(\phi)-Q_x'(\phi)Q_y”(\phi)\right]}
{\left(Q_x’^2(\phi)+Q_y’^2(\phi)\right)^{3/2}}<0.
$$
Because the derivative is negative, the angle \(\eta\) monotonically increases with \(\phi\). This fact indirectly indicates that the spatial tangent vector \(\mathbf{S}_\phi\) rotates continuously in one direction. Therefore, the theoretical tooth surface of a straight bevel gear is a convex ruled surface. This convex property guarantees that all tangent planes along the cone generatrices lie on the same side of the tooth surface and thus a polyhedral surface formed by a finite number of tangent planes can approximate the theoretical tooth surface from the outside.
5. Envelope Shaping Principle
Based on the convex property of the tooth surface, I propose an envelope shaping principle for the outsize straight bevel gear. The basic idea is to use a sequence of planes tangent to the theoretical tooth surface along different cone generatrices. Adjacent planes intersect in straight lines, forming a polyhedral surface. When the number of planes is sufficiently large, the polyhedral surface approximates the theoretical tooth surface with acceptable accuracy.
In the machining process, the straight cutting edge of a planing tool sweeps a plane while moving along a cone generatrix of the work gear. During each cutting stroke, the cutting edge is kept tangent to the theoretical tooth surface. Multiple cutting strokes generate multiple planes, and these planes together form the final machined tooth flank. In my method, I choose a rhombic indexable insert with a straight cutting edge. The cutter is mounted on a CNC planing machine, and the workpiece is rotated continuously to realize efficient machining.
This envelope shaping principle is different from traditional generating methods because it does not rely on an imaginary crown gear. Instead, the tooth surface is generated directly by enveloping tangent planes. This makes the machining process more intuitive and reduces the size of the required machine tool. In addition, the calculation of machining parameters is based on the theoretical tooth surface of the straight bevel gear itself, rather than on the approximated equivalent spur gear at the large end. Therefore, theoretical calculation errors are avoided.
6. Machining Method and Tool-Workpiece Relative Position
In the proposed machining method, I choose to cut from the large end toward the small end of the straight bevel gear. This choice avoids unnecessary conversion of gear parameters because the modulus and tooth profile parameters are defined at the large end. The pitch cone generatrix is placed horizontally. The workpiece axis and the pitch cone generatrix lie in a vertical plane, and the machining coordinate system is defined accordingly.
The relative position of the tool and the workpiece is shown in the conceptual arrangement of the machine. The projection plane is the plane that contains the cone apex and is parallel to the plane of the cutting edge. In this projection plane, the workpiece axis projection forms an angle \(\theta\) with the \(x\)-axis of the workpiece coordinate system. The angle \(\theta\) is determined from two parts. The first part is the angle between the radius vector corresponding to the starting point of the spherical involute and the radius vector corresponding to the pitch cone point. The second part is the half-angle corresponding to the pitch-cone tooth thickness. The total projection angle is:
$$
\theta = \theta_1+\theta_2
=
\arccos\left[\frac{Q_x(\phi_\delta)}{R\sin\delta}\right]
+
\frac{\pi m}{4R\sin\delta},
$$
where \(\phi_\delta\) is the value of \(\phi\) corresponding to the pitch cone generatrix, \(m\) is the large-end modulus, \(\delta\) is the pitch cone angle, and \(\theta_1\), \(\theta_2\) have the following meanings:
$$ \theta_1=\arccos\left[\frac{Q_x(\phi_\delta)}{R\sin\delta}\right], \quad
\theta_2=\frac{\pi m}{4R\sin\delta}.
$$
The projection plane is therefore defined by the normal vector \(\mathbf{n}\) and the cone apex:
$$ x\sin\delta\cos\theta + y\sin\delta\sin\theta + z\cos\delta = 0. $$
7. Cutting-In Position Adjustment Motions
Before each cutting stroke, the tool cutting edge must be positioned so that it is tangent to the large-end tooth profile at the current cut-in point. I divide the angular range between the face cone and the root cone into \(k-1\) equal intervals, where \(k\) is the number of enveloping planes. For the \(n\)-th cut-in point, the corresponding cone apex angle is:
$$ \delta_n = \delta_a – \frac{n(\delta_a-\delta_f)}{k-1}, \quad n=0,1,\ldots,k-1, $$
where \(\delta_a\) is the face cone angle and \(\delta_f\) is the root cone angle. For each cone apex angle \(\delta_n\), the parameter \(\phi_n\) of the large-end spherical involute is obtained from the geometry:
$$ \phi_n = \frac{\arccos\left(\dfrac{\cos\delta_n}{\cos\delta_b}\right)}{\sin\delta_b}. $$
Substituting \(\phi_n\) into \(\mathbf{Q}(\phi)\) gives the coordinates of the cut-in point on the large-end tooth profile:
$$ \mathbf{P}_n = \mathbf{Q}(\phi_n). $$
When the cone generatrix at the cut-in point is brought into coincidence with the workpiece axis projection in the projection plane, the workpiece must rotate by an additional angle \(\Delta\theta_n\). In traditional calculations, the chordal tooth thickness at the large end is approximated by the equivalent spur gear parameters. In my research, I derive the chordal tooth thickness directly from the spherical involute model. I introduce an auxiliary coordinate system \(S_2(x_2,y_2,z_2)\) whose \(z_2\)-axis is the same as the workpiece \(z\)-axis and whose \(x_2\)-axis is rotated by \(\theta\) from the workpiece \(x\)-axis. The large-end spherical involute is transformed into this auxiliary system, and the \(y_2\)-coordinate gives the lateral distance from the tooth-profile center. The chordal tooth thickness at the cut-in point is:
$$ S_n = 2\left[Q_y(\phi_n)\cos\theta – Q_x(\phi_n)\sin\theta\right]. $$
Then the workpiece rotation angle is:
$$ \Delta\theta_n = \arcsin\left(\frac{S_n}{2R\sin\delta_n}\right). $$
After the workpiece rotation, the large-end tooth profile in the workpiece coordinate system is expressed as \(\mathbf{P}(\phi)\):
$$
\mathbf{P}(\phi) = \begin{bmatrix}
\cos\Delta\theta_n & \sin\Delta\theta_n & 0 \\
-\sin\Delta\theta_n & \cos\Delta\theta_n & 0 \\
0 & 0 & 1
\end{bmatrix}
\mathbf{Q}(\phi).
$$
By projecting this spatial curve onto the projection plane, I obtain the projected tooth profile. For an arbitrary point on the tooth profile, the projection is calculated by drawing a perpendicular line to the projection plane:
$$
\begin{bmatrix} x_t \\ y_t \\ z_t \end{bmatrix}
=
\begin{bmatrix}
P_x(\phi) – s_0 R\sin\delta\cos\theta \\
P_y(\phi) – s_0 R\sin\delta\sin\theta \\
P_z(\phi) – s_0 R\cos\delta
\end{bmatrix},
$$
where
$$ s_0 = \frac{P_x(\phi)\sin\delta\cos\theta + P_y(\phi)\sin\delta\sin\theta + P_z(\phi)\cos\delta}{R}. $$
The tangent direction of the projected tooth profile at the cut-in point is obtained by differentiating the above expression with respect to \(\phi\). The angle \(\lambda_n\) between this tangent and the workpiece axis projection is:
$$
\lambda_n = \arccos\left[
\frac{(x_t’\cos\theta+y_t’\sin\theta)\cos\delta+z_t’\sin\delta}
{\sqrt{x_t’^2+y_t’^2+z_t’^2}}
\right].
$$
The initial swing angle of the tool is related to \(\lambda_n\) and to the blade included angle \(\gamma\). For the first cut-in point, I set:
$$ \beta_0 = \frac{\gamma}{2}-\lambda_0, $$
and for the subsequent cut-in points, the swing angle is determined by the change of the tangent direction:
$$ \beta_n = \left|\lambda_n-\lambda_{n-1}\right|. $$
The tool swing motion is realized about the tool shank center. The horizontal position of the tool in the machine \(X\)-direction is adjusted so that the blade center lies on the workpiece axis projection. The adjustment value is:
$$ H_{X,n} = l’\sin\beta_n, $$
where \(l’\) is the swing radius of the blade center. The vertical position in the machine \(Y\)-direction is also adjusted according to the distance between the blade center and the cone apex projection in the projection plane. These four motions, namely workpiece rotation, tool swing, \(X\)-axis adjustment, and \(Y\)-axis adjustment, constitute the cutting-in position adjustment motions.
8. Tooth Surface Forming Motions
During the forming process, the workpiece is continuously rotating, and the tool cutting edge must follow the cone generatrix while remaining tangent to the theoretical tooth surface. Therefore, the forming process requires five coordinated motions: the ram primary motion in the \(Z\)-direction, the continuous workpiece rotation, the tool swing following the ram motion, the \(X\)-axis following motion, and the \(Y\)-axis following motion.
I divide the cone generatrix at the cut-in point into \(N\) discrete points. When the workpiece rotates by an angle \(\omega_1 t\), a point on the cone generatrix has moved in the \(Z\)-direction. The total ram displacement at the \(N\)-th point is:
$$ L_Z(N) = L_Z'(N) + r_{n(N)}\left(1-\cos\omega_1 t\right)\sin\delta, $$
where \(L_Z'(N)\) is the ram displacement when the workpiece is fixed, and \(r_{n(N)}\) is the radius of the cone cross-section at that point. In practice, the workpiece rotational speed \(\omega_1\) must be selected based on the maximum ram stroke and the allowable cutting speed of the machine tool.
The tool must swing continuously as it moves along the tooth surface. The instantaneous swing angle \(\beta_{n(N)}\) is computed from the projected tooth profile at the corresponding point. The projected coordinates of a point on the theoretical tooth surface are obtained by transforming the spherical involute at an arbitrary cone distance into the workpiece coordinate system and then projecting it onto the projection plane. The instantaneous tangent direction is used to calculate \(\lambda_{n(N)}\), and the swing angle becomes:
$$ \beta_{n(N)} = \lambda_{n(N)}-\lambda_n. $$
To keep the contact point of the cutting edge on the same cone generatrix during workpiece rotation, the tool must also move in the machine \(Y\)-direction. The relationship between the \(Y\)-axis following motion and the ram motion depends on whether the tool swing angle is larger than the blade contact angle \(\eta_n\). In the case \(\beta_{n(N)}\ge\eta_n\), the \(Y\)-direction displacement is:
$$
\begin{aligned}
L_Y(N) = \pm\Big[&L_Z(N)\tan(\delta-\delta_n) – r_{n(N)}(1-\cos\omega_1 t)\cos\delta \\
&- l’\left(\cos(\beta_{n(N)}-\eta_n)-\cos(\beta_{n(N)}+\eta_n)\right)\Big],
\end{aligned}
$$
where the positive sign is used when cutting the tooth flank from the addendum to the pitch cone, and the negative sign is used when cutting from the pitch cone to the dedendum. For the case \(\beta_{n(N)}<\eta_n\), a similar expression is applied with the appropriate sign changes. The \(X\)-direction following motion is:
$$
L_X(N) = r_{n(N)}\sin\omega_1 t
+ l’_{n(N)}\left[\sin(\beta_{n(N)}+\eta_n)-\sin(\beta_{n(N)}-\eta_n)\right].
$$
These motion equations fully describe the tooth surface forming process. The cutting edge of the tool sweeps a plane while the workpiece rotates continuously. By synthesizing the five motions, I can generate a series of planes tangent to the theoretical tooth surface along different cone generatrices, and the resulting polyhedral surface is the machined tooth flank of the straight bevel gear.
9. Mathematical Model of the Enveloping Plane
In order to control the machining accuracy, I need to obtain the mathematical equation of each enveloping plane. A plane is uniquely determined by the cutting edge line and the cone generatrix along which the tool moves. In the workpiece coordinate system, the direction vector of the cutting edge is \(\mathbf{T}_n\), and the direction vector of the cone generatrix is \(\mathbf{F}_n = \mathbf{P}_n\). The normal vector of the generated plane is:
$$ \mathbf{n}_n = \mathbf{T}_n \times \mathbf{F}_n. $$
Then the equation of the \(n\)-th plane is:
$$ n_{nx}\left(x-P_{nx}(\phi_n)\right) + n_{ny}\left(y-P_{ny}(\phi_n)\right) + n_{nz}\left(z-P_{nz}(\phi_n)\right)=0. $$
After all \(k\) planes are generated, the adjacent planes intersect in straight lines. The intersection line between the \(i\)-th and \((i+1)\)-th planes is:
$$
L_i:\quad \frac{x-x_0}{A_i}=\frac{y-y_0}{B_i}=\frac{z-z_0}{C_i}, \quad i=0,1,\ldots,k-2,
$$
where the direction vector \((A_i,B_i,C_i)\) of the intersection line is the cross product of the two plane normal vectors:
$$
A_i = n_{iy}n_{(i+1)z}-n_{iz}n_{(i+1)y}, \quad
B_i = n_{iz}n_{(i+1)x}-n_{ix}n_{(i+1)z}, \quad
C_i = n_{ix}n_{(i+1)y}-n_{iy}n_{(i+1)x}.
$$
10. Envelope Accuracy Control
The polyhedral surface formed by the enveloping planes is an approximation of the theoretical tooth surface of the straight bevel gear. To evaluate the approximation error, I define the envelope accuracy as the maximum distance between the intersection lines of adjacent planes and the theoretical tooth surface. This distance is a reasonable measure because the largest deviation usually occurs along the intersection line.
In my computation, I use a grid method combined with iterative refinement. The main steps are as follows:
- Discretize each intersection line \(L_i\) into \(M\) points.
- Discretize the theoretical tooth surface \(\mathbf{S}(r,\phi)\) into grid points with a given step size.
- For each discrete point on the intersection line, find the nearest grid point on the theoretical tooth surface, and use it as the initial iteration point.
- Reduce the grid step by half and search for a new nearest grid point in the neighborhood of the previous point.
- Repeat the refinement until the grid step is smaller than the prescribed tolerance \(\varepsilon\).
- The obtained distance is the minimum distance from the discrete point to the tooth surface.
After processing all discrete points on all intersection lines, I obtain a distance set:
$$ \{d_1,d_2,\ldots,d_j\}, \quad j=1,2,\ldots,M(k-1). $$
The envelope accuracy \(\mu\) is defined as:
$$ \mu = \max\{d_1,d_2,\ldots,d_j\}. $$
If \(\mu\) is larger than the specified tolerance, I increase the enveloping number \(k\) and repeat the whole calculation. The smallest \(k\) that satisfies the accuracy requirement is selected as the optimal enveloping number. This approach gives a direct relationship between the machining efficiency and the required tooth surface accuracy, which is useful for practical production planning.
11. Machine Motion Layout
Based on the required motions, I propose a rational machine motion layout for the outsize straight bevel gear planing machine. The machine coordinate system is defined as follows: the \(Z\)-axis is the ram primary motion direction, the \(X\)-axis is horizontal and perpendicular to the ram direction, and the \(Y\)-axis is vertical. The workpiece is mounted on a rotary table, and its axis is adjusted according to the pitch cone angle of the straight bevel gear.
Because the workpiece is very large and heavy, it is not desirable to impose fast following motions on the workpiece table. Therefore, I assign the \(X\)-axis and \(Y\)-axis following motions to the tool holder. The workpiece only rotates continuously at a low speed. The tool holder carries the ram motion, the tool swing motion, the \(X\)-axis following motion, and the \(Y\)-axis following motion. This arrangement reduces the inertia of the moving parts and improves the dynamic response of the machine.
The motion allocation is summarized in the following table:
| Axis / Element | Motion | Function |
|---|---|---|
| Workpiece table | Continuous rotation | Indexing and tooth surface generation |
| Tool holder | Ram reciprocating motion in \(Z\)-direction | Primary cutting motion |
| Tool holder | Tool swing about a vertical axis | Keeping the cutting edge tangent to the tooth surface |
| Tool holder | \(X\)-axis following motion | Compensating the lateral movement of the cone generatrix |
| Tool holder | \(Y\)-axis following motion | Compensating the vertical movement of the cone generatrix |
This machine layout can be adapted to different sizes of straight bevel gears by adjusting the angle of the worktable. The proposed layout provides a practical basis for the design of a special CNC planing machine for outsize straight bevel gears.
12. Software Development for Calculation
The calculation of the machining parameters and envelope accuracy is complicated and time-consuming. In my research, I develop a special calculation software system for the envelope shaping method of outsize straight bevel gears. The software combines the numerical computation capability of MATLAB with the convenient graphical user interface design of Visual Basic. The main program is written in Visual Basic, while the core numerical algorithms are written as MATLAB M-files. These M-files are compiled into COM components using MATLAB Builder for COM, and then called from the Visual Basic environment.
The software system includes the following functional modules:
| Module | Description |
|---|---|
| Parameter input module | Input of workpiece geometric parameters, machine parameters, and cutting tool parameters |
| Tooth surface calculation module | Calculation of the theoretical tooth surface and large-end tooth profile |
| Cut-in point module | Calculation of the cone apex angles and corresponding cut-in point positions |
| Projection plane module | Calculation of the projection plane normal vector and plane equation |
| Motion calculation module | Calculation of cutting-in adjustment motions and forming motions |
| Envelope accuracy module | Calculation of plane equations, intersection lines, and maximum tooth surface error |
| Database management module | Storage, loading, and deletion of part parameters and machining results |
| Result display module | Display of the optimal enveloping number and all associated machining parameters |
The main calculation procedure is shown in the following flow. First, the user inputs the gear parameters, tool parameters, and machine parameters. Then the software calculates the theoretical tooth surface and the large-end tooth profile. An initial enveloping number \(k\) is selected, and the cut-in points are determined. The projection plane and the projected tooth profile are calculated. Afterward, the cutting-in position adjustment motions and the tooth surface forming motions are derived. The software then computes the enveloping planes, their intersection lines, and the envelope accuracy \(\mu\). If \(\mu\) is larger than the required accuracy, the enveloping number is increased and the calculation is repeated. Finally, the software outputs the optimal enveloping number and the complete set of machining parameters.
The use of MATLAB and Visual Basic hybrid programming greatly reduces the software development cycle. It also provides a user-friendly interface for shop-floor engineers who need to use the proposed machining method in practice.
13. Calculation Example
To verify the correctness and feasibility of the proposed method, I perform a numerical calculation example for an outsize straight bevel gear. The main geometric parameters of the workpiece are listed in the following table.
| Symbol | Parameter | Value |
|---|---|---|
| \(z\) | Number of teeth | 200 |
| \(m\) | Large-end modulus | 40 mm |
| \(\alpha\) | Pressure angle | 20° |
| \(\delta\) | Pitch cone angle | 84°43′ |
| \(B\) | Face width | 400 mm |
| \(d\) | Large-end pitch diameter | 8000 mm |
| \(R\) | Cone distance | 4017.07 mm |
| \(\delta_a\) | Face cone angle | 85°17’26” |
| \(\delta_f\) | Root cone angle | 84°02’07” |
I select a rhombic indexable insert with included angle \(\gamma=35^\circ\). The required envelope accuracy is 0.03 mm. For the initial calculation, I choose the enveloping number \(k=3\). The corresponding cut-in point parameters are calculated and listed in the following table.
| Cut-in point \(n\) | Cone apex angle \(\delta_n\) | Parameter \(\phi_n\) | Initial tool swing angle \(\beta_n\) |
|---|---|---|---|
| 0 | 85°17’26” | 81°48’05” | 2°38’22” |
| 1 | 84°39’47” | 79°50’47” | 2°28’44” |
| 2 | 84°02’07” | 77°52’34” | 2°15’47” |
Using the equations derived in the previous sections, I calculate the projection plane equation and the projected large-end tooth profile. The cutting-in adjustment motions and the forming motions are then obtained. A portion of the calculated following-motion displacements is shown in the following table.
| Discrete point number | \(L_X\) for first plane (mm) | \(L_X\) for second plane (mm) | \(L_Z\) for first plane (mm) | \(L_Z\) for second plane (mm) |
|---|---|---|---|---|
| 1 | 0 | 0 | 0 | 0 |
| 5 | 1.644 | 1.643 | 84.207 | 84.211 |
| 10 | 3.600 | 3.597 | 189.466 | 189.475 |
| 15 | 5.446 | 5.441 | 294.726 | 294.741 |
| 20 | 7.182 | 7.175 | 399.987 | 400.007 |
The tool swing angle \(\beta_{n(N)}\) at different discrete points is also calculated. The values increase monotonically as the tool moves from the large end to the small end, which is necessary to compensate for the continuous rotation of the workpiece.
After calculating the three enveloping planes for \(k=3\), I obtain their mathematical equations. The intersection lines between adjacent planes are then derived. The envelope accuracy \(\mu\) is calculated using the grid refinement method. For \(k=3\), the maximum deviation between the polyhedral surface and the theoretical tooth surface is:
$$ \mu = 0.033\ \text{mm}. $$
This value is larger than the required accuracy of 0.03 mm. Therefore, I increase the enveloping number to \(k=4\) and repeat the entire calculation. The new envelope accuracy is:
$$ \mu = 0.009\ \text{mm}. $$
This result satisfies the accuracy requirement and leaving enough margin for practical cutting process variations. Therefore, the optimal enveloping number is \(k=4\). The corresponding machining parameters are used as the basis for subsequent CNC programming.
The calculation example demonstrates that the proposed envelope shaping method can accurately generate the tooth surface of an outsize straight bevel gear. The workpiece is continuously rotated, the tool path is smooth, and the required accuracy can be achieved by simply adjusting the number of enveloping planes.
14. Discussion of the Proposed Method
The proposed envelope shaping method has several advantages over traditional machining methods for outsize straight bevel gears. First, the method avoids the use of a large imaginary generating gear. This significantly reduces the required machine size and cost. Second, the calculation of the machining parameters is based directly on the spherical involute tooth surface model, so the theoretical errors introduced by equivalent spur gear approximation are eliminated. Third, the workpiece rotates continuously during the entire cutting process. This continuous motion not only improves machining efficiency but also reduces the acceleration and deceleration requirements of the workpiece driving motor. Fourth, the use of indexable rhombic inserts with straight cutting edges is economical and convenient. The inserts can be easily replaced when worn, and the same tool holder can be used for different gear sizes.
There are also some limitations that should be addressed in future work. Because the machine tool requires five coordinated motions, the control system must be carefully designed to achieve the desired synchronization. The envelope accuracy depends on the number of planes, but increasing the number of planes also increases the machining time. Therefore, an optimal choice must be made according to the required tooth surface accuracy and production efficiency. In addition, the current research is mainly based on numerical simulation and theoretical calculation. Experimental validation on a real outsize straight bevel gear machine should be carried out in the future.
15. Conclusion
In my research, I propose a complete envelope shaping method for the tooth surface of an outsize straight bevel gear. The main conclusions are as follows:
First, the theoretical tooth surface of a straight bevel gear is a convex ruled surface. I establish the mathematical model using the spherical involute curves at the large end and small end as the two base curves and the cone generatrices as the straight lines. The convexity analysis proves that the tangential direction of the tooth surface rotates monotonically, which is the theoretical basis for the plane-envelope machining principle.
Second, I propose an envelope planing principle in which a straight cutting edge sweeps a plane along a cone generatrix. Multiple planes are generated sequentially, and adjacent planes intersect to form a polyhedral surface. This polyhedral surface approximates the theoretical tooth surface of the straight bevel gear. The machining process starts from the large end and proceeds toward the small end, and the pitch cone generatrix is oriented horizontally to simplify the calculation and the machine setup.
Third, the tooth surface forming process requires five coordinated motions: workpiece rotation, ram primary motion, tool swing, \(X\)-axis following motion, and \(Y\)-axis following motion. I derive the mathematical model of each motion and give the corresponding formulas. Based on these motion models, I propose a machine motion layout in which the workpiece only rotates continuously while all following motions are allocated to the tool holder. This layout reduces the effect of workpiece inertia and improves the dynamic performance of the machine.
Fourth, the envelope accuracy is controlled by the number of enveloping planes. I calculate the intersection lines between adjacent planes and compute the maximum distance between the intersection lines and the theoretical tooth surface using a grid refinement method. The optimal enveloping number is determined by comparing the calculated accuracy with the specified tolerance.
Fifth, I develop a calculation software system for the envelope shaping method using MATLAB and Visual Basic hybrid programming. The software greatly simplifies the calculation process and provides a user-friendly interface for practical application.
Finally, the calculation example shows that for an outsize straight bevel gear with a pitch diameter of 8000 mm, the initial enveloping number \(k=3\) gives an envelope accuracy of 0.033 mm, while increasing the enveloping number to \(k=4\) gives an accuracy of 0.009 mm, which satisfies the requirement of 0.03 mm. The example confirms the correctness and feasibility of the proposed envelope shaping method for outsize straight bevel gears.
Future work will include the construction of a prototype planing machine, actual cutting experiments, vibration analysis of the machine tool, optimization of the cutting parameters, and further improvement of the calculation software. I believe that the proposed method will provide a new and effective approach for machining outsize straight bevel gears and will contribute to the advancement of large gear manufacturing technology.