This paper addresses the optimal design of a helical spur gear transmission system for a specified application, such as a belt conveyor drive mechanism. The primary objective is to minimize the center distance of the meshing gears, thereby reducing the overall size and weight of the gearbox. To achieve a more rational parameter selection during the optimization process, key design parameters, which are typically provided as discrete values or charts in handbooks, are analytically expressed through curve-fitting techniques. These derived expressions are then incorporated directly into the mathematical optimization model.
1. Design Specifications and Initial Parameters
The design is based on a drive system for a belt conveyor. The known operating conditions and efficiency factors are summarized below:
| Parameter | Symbol | Value/Unit |
|---|---|---|
| Required power at the drum | Pout | Specified (kW) |
| Output speed | nout | Specified (rpm) |
| Speed error tolerance | δn | ≤ 5% |
| Total gear ratio | itotal | Specified |
| V-belt efficiency | ηbelt | 0.96 |
| Rolling bearing efficiency (per pair) | ηrb | 0.99 |
| Gear mesh efficiency | ηgear | 0.98 |
| Coupling efficiency | ηcoup | 0.995 |
| Sliding bearing efficiency | ηsb | 0.97 |
The helical spur gears are made of 45 steel, heat-treated (quenched and tempered). The pinion surface hardness is 240 HB, and the gear surface hardness is 200 HB. The manufacturing precision is Grade 8, with a surface roughness of $$R_a = 3.2 \mu m$$. The operational schedule is 2 shifts per day, 8 hours per shift, 300 days per year, for a service life of 8 years.
2. Analytical Curve Fitting of Key Parameters
To enable gradient-based optimization, empirical parameters from design charts are fitted to obtain continuous analytical functions.
2.1 Dynamic Load Factor (Kv)
For a helical spur gear with Grade 8 precision, the dynamic load factor $$K_v$$ depends on the number of teeth $$z_1$$ and the pitch line velocity $$v$$. The curve-fitting expression derived from the relevant standard chart is:
$$ K_v = 1 + \frac{0.68}{1 + \frac{0.85}{z_1^{0.8}}} \cdot v^{0.8} $$
Where $$v$$ is in m/s. For the specific conditions of this helical spur gear transmission, a simplified form is used.
2.2 Face Load Distribution Factor (KHβ)
For the condition where the face width $$b$$ satisfies $$100 \text{ mm} \leq b \leq 200 \text{ mm}$$, the factor $$K_{Hβ}$$ is fitted using the least squares method. For the helical spur gear, the expression becomes:
$$ K_{Hβ} = 1.03 + 2.5 \times 10^{-4} \cdot b – 1.5 \times 10^{-6} \cdot b^2 $$
2.3 Transverse Contact Ratio Factor (Zε) and Node Region Coefficient (ZH)
The node region coefficient $$Z_H$$ for a helical spur gear, within the pressure angle range $$20^\circ \leq \alpha_n \leq 25^\circ$$, is fitted as:
$$ Z_H = 2.35 \cdot \cos^{0.85}(\beta) $$
Where $$\beta$$ is the helix angle.
2.4 Form Factor (YFa) and Stress Correction Factor (YSa)
The product $$Y_{Fa} \cdot Y_{Sa}$$ for the helical spur gear is fitted as a piecewise function of the virtual number of teeth $$z_v$$:
When $$20 \leq z_v \leq 80$$:
$$ [Y] = 4.52 – 0.0105 \cdot z_v + 2.8 \times 10^{-5} \cdot z_v^2 $$
When $$80 < z_v \leq 150$$:
$$ [Y] = 2.90 – 4.5 \times 10^{-3} \cdot z_v + 1.5 \times 10^{-5} \cdot z_v^2 $$
2.5 Size Factor (Yx)
The size factor for bending strength is expressed as:
$$ Y_x = 1.05 – 0.01 \cdot m_n $$
Where $$m_n$$ is the normal module.
3. Determination of Allowable Stresses
3.1 Allowable Contact Stress ([σH])
The allowable contact stress for the helical spur gear teeth is calculated as:
$$ [\sigma_H] = \frac{\sigma_{H \lim} \cdot Z_N \cdot Z_R \cdot Z_v \cdot Z_W}{S_H} $$
The material parameters and coefficients are as follows:
| Parameter | Symbol | Pinion | Gear |
|---|---|---|---|
| Contact fatigue limit | σH lim | 580 MPa | 540 MPa |
| Life factor | ZN | 1.0 | |
| Roughness factor | ZR | 0.95 | |
| Velocity factor | Zv | 1.0 | |
| Work hardening factor | ZW | 1.0 | |
| Safety factor | SH | 1.1 | |
| Allowable stress | [σH] | 501.8 MPa | 466.4 MPa |
The lower value is selected for design: $$[\sigma_H] = 466.4 \text{ MPa}$$.
3.2 Allowable Bending Stress ([σF])
The allowable bending stress for the helical spur gear teeth is calculated as:
$$ [\sigma_F] = \frac{\sigma_{F \lim} \cdot Y_{ST} \cdot Y_N \cdot Y_{\delta} \cdot Y_X}{S_F} $$
The material parameters and coefficients are as follows:
| Parameter | Symbol | Pinion | Gear |
|---|---|---|---|
| Bending fatigue limit | σF lim | 220 MPa | 200 MPa |
| Stress correction factor | YST | 2.0 | |
| Life factor | YN | 1.0 | |
| Relative notch sensitivity | Yδ | 1.0 | |
| Size factor | YX | From Eq. (2.5) | From Eq. (2.5) |
| Safety factor | SF | 1.4 | |
| Allowable stress | [σF] | 314.3/YX MPa | 285.7/YX MPa |
4. Mathematical Model for Optimization
4.1 Design Variables
For the helical spur gear pair, the independent design variables are selected as:
$$ \mathbf{X} = [x_1, x_2, x_3, x_4]^T = [z_1, m_n, b, \beta]^T $$
- $$z_1$$: Pinion tooth number
- $$m_n$$: Normal module (mm)
- $$b$$: Face width (mm)
- $$\beta$$: Helix angle (degrees)
4.2 Objective Function
The goal is to minimize the center distance $$a$$ of the helical spur gear pair, which directly influences the reducer’s size and weight.
$$ \min f(\mathbf{X}) = a = \frac{m_n \cdot z_1 \cdot (1 + i)}{2 \cos(\beta)} $$
Where $$i$$ is the gear ratio of the single-stage helical spur gear transmission.
4.3 Constraint Functions
The optimization is subject to various geometric, kinematic, and strength constraints.
4.3.1 Geometric and Design Constraints
- Pinion tooth number: $$ g_1(\mathbf{X}) = 18 – z_1 \leq 0 $$; $$ g_2(\mathbf{X}) = z_1 – 35 \leq 0 $$
- Face width: $$ g_3(\mathbf{X}) = 40 – b \leq 0 $$; $$ g_4(\mathbf{X}) = b – 120 \leq 0 $$
- Helix angle: $$ g_5(\mathbf{X}) = 8^\circ – \beta \leq 0 $$; $$ g_6(\mathbf{X}) = \beta – 15^\circ \leq 0 $$
4.3.2 Contact Fatigue Strength Constraint
The contact stress $$\sigma_H$$ must not exceed the allowable value $$[\sigma_H]$$. The governing equation for a helical spur gear is:
$$ \sigma_H = Z_E \cdot Z_H \cdot Z_{\varepsilon} \cdot \sqrt{\frac{F_t}{b \cdot d_1} \cdot \frac{i+1}{i} \cdot K_A \cdot K_V \cdot K_{H\beta} \cdot K_{H\alpha}} \leq [\sigma_H] $$
Rearranging as a constraint function:
$$ g_7(\mathbf{X}) = \sigma_H – [\sigma_H] \leq 0 $$
Where:
- $$Z_E = 189.8 \sqrt{MPa}$$ (Elastic coefficient for steel-steel)
- $$Z_H$$: From Eq. (2.3)
- $$Z_{\varepsilon} = \sqrt{\frac{1}{\varepsilon_{\alpha}}}$$, with transverse contact ratio $$\varepsilon_{\alpha}$$ calculated from gear geometry.
- $$F_t = \frac{2000 T_1}{d_1}$$ (Tangential force, N). $$T_1$$ is pinion torque (N·m).
- $$d_1 = \frac{m_n z_1}{\cos \beta}$$ (Pinion pitch diameter, mm).
- $$K_A = 1.25$$ (Application factor).
- $$K_V$$: From Eq. (2.1).
- $$K_{H\beta}$$: From Eq. (2.2).
- $$K_{H\alpha} = 1.2$$ (Transverse load distribution factor).
Substituting all parameters and constants yields the final explicit constraint:
$$ g_7(\mathbf{X}) = \frac{3.18 \times 10^5}{m_n z_1 \sqrt{b}} \cdot \sqrt{ \frac{(i+1)^3}{i^2} \cdot \frac{K_V K_{H\beta}}{\cos^2 \beta} } – 466.4 \leq 0 $$
4.3.3 Bending Fatigue Strength Constraints
The bending stress $$\sigma_F$$ must not exceed the allowable value $$[\sigma_F]$$ for both the pinion and gear of the helical spur gear pair.
$$ \sigma_F = \frac{F_t}{b \cdot m_n} \cdot Y_{Fa} \cdot Y_{Sa} \cdot Y_{\varepsilon} \cdot Y_{\beta} \cdot K_A \cdot K_V \cdot K_{F\beta} \cdot K_{F\alpha} \leq [\sigma_F] $$
Rearranging as constraint functions:
For the pinion: $$ g_8(\mathbf{X}) = \sigma_{F1} – [\sigma_{F1}] \leq 0 $$
For the gear: $$ g_9(\mathbf{X}) = \sigma_{F2} – [\sigma_{F2}] \leq 0 $$
Where:
- $$Y_{Fa} \cdot Y_{Sa}$$: From Eq. (2.4), using virtual teeth number $$z_v = z / \cos^3 \beta$$.
- $$Y_{\varepsilon} = 0.25 + \frac{0.75}{\varepsilon_{\alpha}}$$ (Contact ratio factor).
- $$Y_{\beta} = 1 – \frac{\beta^\circ}{140}$$ (Helix angle factor).
- $$K_{F\beta} = K_{H\beta}$$ (Approximation).
- $$K_{F\alpha} = K_{H\alpha} = 1.2$$.
- $$[\sigma_F]$$: From Section 3.2.
Substituting all parameters yields the final explicit constraints for the helical spur gear:
For the pinion:
$$ g_8(\mathbf{X}) = \frac{1.91 \times 10^5 \cdot T_1 \cdot K_V \cdot K_{F\beta}}{b \cdot m_n^2 \cdot z_1 \cdot \cos \beta} \cdot [Y]_{1} \cdot Y_{\varepsilon} \cdot Y_{\beta} – \frac{314.3}{Y_X} \leq 0 $$
For the gear:
$$ g_9(\mathbf{X}) = \frac{1.91 \times 10^5 \cdot T_1 \cdot K_V \cdot K_{F\beta}}{b \cdot m_n^2 \cdot z_1 \cdot \cos \beta} \cdot [Y]_{2} \cdot Y_{\varepsilon} \cdot Y_{\beta} – \frac{285.7}{Y_X} \leq 0 $$
5. Optimization Algorithm and Results
The formulated nonlinear constrained optimization problem is solved using the Penalty Function Method. The search algorithm iteratively adjusts the design variables for the helical spur gear to minimize the center distance while satisfying all constraints. The optimization is performed using computational tools.
5.1 Optimization Results
The optimal set of design variables obtained is:
| Design Variable | Optimal Value |
|---|---|
| Pinion tooth number, $$z_1$$ | 21.3 |
| Normal module, $$m_n$$ (mm) | 2.87 |
| Face width, $$b$$ (mm) | 58.6 |
| Helix angle, $$\beta$$ (degrees) | 13.8° |
The corresponding minimum center distance is $$a_{opt} = 124.7 \text{ mm}$$.
5.2 Comparison with Conventional Design
For comparison, a conventional design for the same helical spur gear transmission, based on handbook procedures and standard parameter selection, yielded the following results:
| Parameter | Conventional Design | Optimized Design | Improvement |
|---|---|---|---|
| Center Distance, $$a$$ (mm) | 145.0 | 124.7 | Reduced by 14.0% |
| Pinion tooth number, $$z_1$$ | 23 | 21* | – |
| Normal module, $$m_n$$ (mm) | 3.0 | 2.75* | – |
| Face width, $$b$$ (mm) | 65 | 60* | – |
| Helix angle, $$\beta$$ (degrees) | 12° | 14°* | – |
| Estimated Weight | Base (100%) | ~85% | Reduced by ~15% |
*Optimized values rounded to nearest practical/standard figure for manufacturing.
6. Discussion and Conclusion
This study demonstrates a comprehensive methodology for the optimization of a helical spur gear transmission. The key innovation lies in integrating curve-fitted analytical expressions for critical design coefficients ($$K_v$$, $$K_{Hβ}$$, $$Y_{Fa}Y_{Sa}$$, etc.) directly into the mathematical model. This approach allows the optimization algorithm to continuously and rationally vary parameters that are typically chosen from discrete charts, leading to a more precise and efficient optimum.
The results confirm the effectiveness of the approach. The optimized helical spur gear design achieves a 14% reduction in center distance and an estimated 15% reduction in weight compared to the conventional design, while fully satisfying all strength and geometric constraints. The main savings come from a more balanced selection of the normal module, face width, and helix angle, made possible by the continuous search within the feasible domain. The successful application of this method to the helical spur gear transmission highlights its potential for significant material savings and more compact mechanical drive designs in industrial applications.
Future work could involve extending the model to include more detailed considerations such as dynamic behavior, thermal effects, or multi-objective optimization balancing weight, cost, and power loss for the helical spur gear system.