The pursuit of manufacturing excellence in gear production necessitates continuous advancements in process efficiency, product quality, and environmental sustainability. Among various gear manufacturing techniques, gear hobbing stands out due to its versatility and productivity. The advent of high-speed dry gear hobbing, which eliminates the need for cutting fluids, presents a significant green manufacturing alternative. However, the absence of coolant exacerbates thermal loads and tool wear, making the optimal selection of process parameters more critical than ever. This work focuses on developing and validating a robust multi-objective optimization framework for high-speed gear hobbing parameters, aiming to simultaneously minimize machining energy consumption and maximize gear quality.
The core challenge lies in the inherent conflict between these objectives. Parameters that enhance surface finish or dimensional accuracy (e.g., reduced feed, specific tool geometry) may increase machining time and energy, while parameters favoring high productivity might compromise quality. Traditional single-objective optimization falls short. Therefore, this article proposes a hybrid methodology integrating an improved Multi-Objective Grey Wolf Optimizer (MOGWO) for generating a Pareto-optimal solution set, followed by the Technique for Order Preference by Similarity to Ideal Solution (TOPSIS) for multi-attribute decision-making to select the final optimal parameters.
Mathematical Modeling of the Gear Hobbing Optimization Problem
The optimization model is constructed with defined decision variables, objective functions, and constraints specific to the high-speed dry gear hobbing process.
Decision Variables
For small-module gear batch processing, the cutting depth has a negligible impact on machine tool deformation and spindle dynamics. Therefore, the primary decision variables chosen for optimization are:
- Axial feed rate, $v_a$ (mm/min) – A continuous variable.
- Spindle speed (hob rotational speed), $n_0$ (r/min) – A continuous variable.
- Hob diameter, $d_0$ (mm) – Treated as a continuous variable within a discrete practical range.
- Number of hob threads (starts), $z_0$ – A discrete variable (e.g., 2 or 3).
The hob gullet number $k$ is considered fixed. Thus, the solution vector is represented as $\mathbf{X} = [v_a, n_0, d_0, z_0]^T$.
Objective Functions
Two conflicting objectives are formulated:
- Minimize Machining Energy Consumption ($E$): The total energy consumed during the active gear hobbing process. A simplified model can relate energy to machining power and time. The specific power $P_c$ can be estimated using empirical models related to cutting forces in gear hobbing.
$$
P_c(v_a, n_0, d_0, z_0) = F_c \cdot v_c
$$
where $F_c$ is the tangential cutting force and $v_c$ is the cutting speed ($v_c = \pi d_0 n_0 / 1000$). The total energy $E$ is:
$$
E = \int P_c(t) \, dt \approx P_c \cdot t_m
$$
where $t_m$ is the machining time for one gear, inversely related to $v_a$ and $z_0$. - Minimize Gear Quality Error ($Q$): Represents a composite error metric, potentially combining deviations in tooth profile, lead, and pitch. This is modeled as a function influenced by process parameters, typically increasing with higher thermal loads (high $n_0$) and improper chip load. A surrogate model or empirical relation is used:
$$
Q = f_Q(v_a, n_0, d_0, z_0)
$$
The function aims to quantify the degradation in gear precision, where a lower $Q$ value indicates higher quality.
The multi-objective optimization problem is formally stated as:
$$
\begin{aligned}
& \text{Minimize} \quad \mathbf{F}(\mathbf{X}) = [E(\mathbf{X}), \, Q(\mathbf{X})]^T \\
& \text{subject to} \\
& v_{a,\text{min}} \leq v_a \leq v_{a,\text{max}} \\
& n_{0,\text{min}} \leq n_0 \leq n_{0,\text{max}} \\
& d_0 \in \{d_{0,1}, d_{0,2}, …, d_{0,m}\} \\
& z_0 \in \{z_{0,1}, z_{0,2}\} \\
& \text{(Other process constraints)}
\end{aligned}
$$
Improved Multi-Objective Grey Wolf Optimizer (MOGWO) for Solution Generation
The standard Grey Wolf Optimizer (GWO) mimics the social hierarchy and hunting behavior of grey wolves. In MOGWO, this metaphor is extended to handle multiple objectives by incorporating an external archive to store non-dominated Pareto solutions and a grid mechanism for archive maintenance. The key improvements in this implementation focus on enhancing exploration and handling mixed variable types.
Algorithm Fundamentals and Enhancement
The wolf pack is categorized into four levels: alpha ($\alpha$), beta ($\beta$), delta ($\delta$), and omega ($\omega$). The $\alpha$, $\beta$, and $\delta$ wolves represent the best, second-best, and third-best solutions, respectively, and guide the hunting (optimization). The position update for a wolf (candidate solution) is guided by these leaders:
$$
\begin{aligned}
\vec{D}_{\alpha} &= |\vec{C}_1 \cdot \vec{X}_{\alpha} – \vec{X}(t)| \\
\vec{D}_{\beta} &= |\vec{C}_2 \cdot \vec{X}_{\beta} – \vec{X}(t)| \\
\vec{D}_{\delta} &= |\vec{C}_3 \cdot \vec{X}_{\delta} – \vec{X}(t)| \\
\vec{X}_1 &= \vec{X}_{\alpha} – \vec{A}_1 \cdot \vec{D}_{\alpha} \\
\vec{X}_2 &= \vec{X}_{\beta} – \vec{A}_2 \cdot \vec{D}_{\beta} \\
\vec{X}_3 &= \vec{X}_{\delta} – \vec{A}_3 \cdot \vec{D}_{\delta} \\
\vec{X}(t+1) &= \frac{\vec{X}_1 + \vec{X}_2 + \vec{X}_3}{3}
\end{aligned}
$$
where $\vec{A}$ and $\vec{C}$ are coefficient vectors calculated as:
$$
\vec{A} = 2\vec{a} \cdot \vec{r}_1 – \vec{a}, \quad \vec{C} = 2 \cdot \vec{r}_2
$$
Here, $\vec{r}_1, \vec{r}_2$ are random vectors in [0,1], and $\vec{a}$ is a convergence factor that decreases linearly from 2 to 0 over iterations. To improve global search capability, the decay of $\vec{a}$ is modified using a non-linear strategy, allowing for more extensive exploration in the early stages and faster convergence later. Additionally, a special mutation operator is applied to discrete variables ($z_0$) to ensure they explore their defined set effectively.
Algorithm Workflow for Gear Hobbing Optimization
The step-by-step workflow of the improved MOGWO for gear hobbing parameter optimization is as follows:
- Initialization: Define population size $N$, maximum iterations $T$, and parameter bounds. Initialize the wolf population $\mathbf{X}_i$ randomly within feasible ranges.
- Evaluation: For each wolf (parameter set), calculate the objective functions: machining energy $E_i$ and gear quality error $Q_i$.
- Archive Management: Initialize an external archive to store non-dominated solutions. Update the archive after each iteration using Pareto dominance principles and a grid-based mechanism to ensure diversity.
- Hierarchy Assignment: Select the alpha, beta, and delta wolves from the current archive using a roulette wheel selection based on grid crowding distance to promote diversity.
- Position Update: Calculate $\vec{A}$ and $\vec{C}$ vectors. Update the position of each omega wolf using the equations guided by $\alpha$, $\beta$, and $\delta$. Apply boundary control and handle discrete variables.
- Iteration: Repeat steps 2-5 until the maximum iteration $T$ is reached.
- Output: The final external archive contains the Pareto-optimal set of gear hobbing parameter solutions.
| Parameter | Symbol | Value/Range |
|---|---|---|
| Population Size | $N$ | 20 |
| Maximum Iterations | $T$ | 200 |
| Number of Grids | – | 10 |
| Inflation Rate | – | 0.1 |
| Variable Bounds: $v_a$ | – | 68 – 92 mm/min |
| Variable Bounds: $n_0$ | – | 650 – 780 r/min |
| Variable Bounds: $d_0$ | – | 78 – 89 mm (discrete steps) |
| Variable Bounds: $z_0$ | – | {2, 3} |
TOPSIS for Multi-Attribute Decision Making on Pareto Solutions
The MOGWO algorithm yields a set of non-dominated solutions (Pareto front), each representing a different trade-off between energy $E$ and quality $Q$. To recommend a single best parameter set for implementation, the TOPSIS method is employed. TOPSIS selects the alternative that has the shortest geometric distance from the positive ideal solution (PIS) and the longest distance from the negative ideal solution (NIS).
Let the Pareto set have $m$ solutions and $n=2$ criteria ($E$, $Q$). The decision matrix $\mathbf{D}$ is constructed:
$$
\mathbf{D} =
\begin{bmatrix}
E_1 & Q_1 \\
E_2 & Q_2 \\
\vdots & \vdots \\
E_m & Q_m
\end{bmatrix}
$$
Step 1: Normalize the decision matrix.
$$
r_{ij} = \frac{x_{ij}}{\sqrt{\sum_{i=1}^{m} x_{ij}^2}}, \quad i=1,\ldots,m; \, j=1,2
$$
Step 2: Calculate the weighted normalized matrix. Assuming equal importance, weights are set as $w_E = w_Q = 0.5$.
$$
v_{ij} = w_j \cdot r_{ij}
$$
Step 3: Determine the ideal and negative-ideal solutions.
$$
\begin{aligned}
A^+ &= \{ v_1^+, v_2^+ \} = \{ (\min_i v_{i1}), (\min_i v_{i2}) \} = \{ v_E^+, v_Q^+ \} \\
A^- &= \{ v_1^-, v_2^- \} = \{ (\max_i v_{i1}), (\max_i v_{i2}) \} = \{ v_E^-, v_Q^- \}
\end{aligned}
$$
Step 4: Calculate the separation measures.
$$
\begin{aligned}
S_i^+ &= \sqrt{ \sum_{j=1}^{2} (v_{ij} – v_j^+)^2 }, \quad i=1,\ldots,m \\
S_i^- &= \sqrt{ \sum_{j=1}^{2} (v_{ij} – v_j^-)^2 }, \quad i=1,\ldots,m
\end{aligned}
$$
Step 5: Calculate the relative closeness to the ideal solution.
$$
C_i = \frac{S_i^-}{S_i^+ + S_i^-}, \quad 0 \leq C_i \leq 1
$$
The solution with the highest $C_i$ value is the best compromise solution according to TOPSIS.
Experimental Validation and Results Analysis
The proposed methodology was validated using a YKS3112CNC7 high-speed gear hobbing machine. The workpiece was a small-module gear with the parameters shown below. Dry hobbing conditions were maintained throughout the experiments.
| Parameter | Value |
|---|---|
| Module | 6.5 mm |
| Number of Teeth | 48 |
| Pressure Angle | 22° |
| Face Width | 32 mm |
The improved MOGWO algorithm was implemented in MATLAB with the parameters from Table 1. After 200 iterations, the algorithm converged to a well-distributed Pareto front.
| Solution ID | $v_a$ (mm/min) | $n_0$ (r/min) | $d_0$ (mm) | $z_0$ | $Q$ (mm) | $E$ (kJ) |
|---|---|---|---|---|---|---|
| P1 | 91.00 | 782.01 | 89 | 2 | 0.035 | 0.16 |
| P2 | 91.00 | 772.39 | 88 | 3 | 0.085 | 0.81 |
| P3 | 91.00 | 781.98 | 89 | 3 | 0.432 | 0.63 |
| P4 | 69.00 | 649.92 | 78 | 2 | 0.143 | 0.72 |
| P8 | 69.00 | 780.00 | 89 | 3 | 0.015 | 10.96 |
| P9 | 69.00 | 780.00 | 89 | 3 | 0.014 | 7.68 |
Analysis of Pareto Solutions:
- Trade-off Trend: A clear trade-off is observed. Solutions with very low gear error (e.g., P8, P9 with $Q$ ~0.015 mm) require significantly higher energy ($E$ > 7 kJ). Conversely, solutions with very low energy (e.g., P1, P2, P3 with $E$ < 1 kJ) result in substantially higher quality error.
- Impact of Hob Diameter ($d_0$): A significant number of Pareto solutions, especially those favoring quality, feature the maximum hob diameter (89 mm). A larger $d_0$ increases tool rigidity, reduces deflection, improves heat dissipation, and can allow for higher metal removal rates, contributing to both quality and productivity in the gear hobbing process.
- Impact of Number of Threads ($z_0$): Solutions with $z_0=3$ are prevalent across the front. Increasing the hob threads reduces machining time per gear ($t_m \propto 1/z_0$), which directly lowers energy consumption for a given power, explaining its dominance in low-energy solutions.
- Impact of Spindle Speed ($n_0$): Higher $n_0$ generally increases cutting temperature. The results show that the best quality solutions (low $Q$) often use high $n_0$ but paired with a low $v_a$, suggesting a strategy to control chip load and thermal effects while maintaining a high cutting speed for better surface finish.
- Impact of Axial Feed ($v_a$): No simple monotonic relationship exists. Optimal $v_a$ is highly dependent on the combination with other parameters. Very high $v_a$ can lead to high energy but poor quality due to increased load, while very low $v_a$ with high $n_0$ can achieve excellent quality at the cost of higher energy due to longer contact time.
TOPSIS Decision and Final Parameter Selection
Applying TOPSIS with equal weights ($w_E=w_Q=0.5$) to the full Pareto set (16 solutions) provides a ranking. The relative closeness $C_i$ is calculated for each.
| Rank | TOPSIS Score ($C_i$) | Corresponding Solution | Primary Characteristic |
|---|---|---|---|
| 1 | 0.1345 | P8 | Best Quality Focus |
| 2 | 0.1503 | P3 | Best Energy Focus |
| 3 | 0.1436 | P6 | Good Compromise |
| 4 | 0.1194 | P9 | High Quality |
| … | … | … | … |
According to TOPSIS, solution P8 ($v_a$=69 mm/min, $n_0$=780 r/min, $d_0$=89 mm, $z_0$=3) is the best compromise, being closest to the ideal of low energy and low error. If the primary goal is supreme gear quality, P8 or P9 should be chosen. If minimizing energy is paramount, P3 is the best choice.
Performance Comparison
To validate the superiority of the proposed improved MOGWO-TOPSIS framework, its recommended solution (P8) was compared against the optimal solution derived from a standard MOGWO algorithm coupled with TOPSIS, under identical constraints and objectives.
| Method | Optimal Parameters ($v_a$, $n_0$, $d_0$, $z_0$) | Gear Error $Q$ (mm) | Energy $E$ (kJ) |
|---|---|---|---|
| Proposed Improved MOGWO-TOPSIS | (69.00, 780.00, 89, 3) | 0.015 | 10.96 |
| Standard MOGWO-TOPSIS | (73.56, 735.29, 89, 3) | 0.042 | 152.62 |
| Improvement | – | -24.6% | -6.3% |
The proposed algorithm successfully identified a parameter set that simultaneously reduced gear error by 24.6% and energy consumption by 6.3% compared to the standard algorithm’s recommendation. This demonstrates the effectiveness of the enhancements made to the MOGWO in navigating the complex, constrained search space of the gear hobbing optimization problem.
Conclusion
This article presented a comprehensive methodology for the multi-objective optimization of high-speed dry gear hobbing process parameters. The integrated approach of an improved MOGWO and TOPSIS provides a powerful tool for manufacturing engineers. The key findings are:
- The inherent conflict between machining energy consumption and gear quality necessitates a multi-objective framework. The improved MOGWO algorithm effectively generates a diverse Pareto-optimal set, revealing the complex trade-offs in gear hobbing.
- Process parameters have distinct and interconnected effects. Larger hob diameters and multi-thread hobs generally promote productivity and can benefit quality or energy depending on other parameter combinations. The optimal axial feed and spindle speed are highly interdependent.
- The TOPSIS method offers a rational and flexible way to select a final implementation point from the Pareto set based on customizable weightings of the objectives, aligning the choice with specific production priorities.
- Experimental validation confirms the practical efficacy of the proposed algorithm. The recommended parameters achieved a significant simultaneous improvement in both gear quality (24.6% error reduction) and energy efficiency (6.3% consumption reduction) compared to a standard optimization approach.
Future work will focus on incorporating more precise physical models for tool wear and dynamic machine tool behavior into the objective functions. Furthermore, exploring real-time adaptation of these parameters based on in-process monitoring signals could pave the way for adaptive, intelligent gear hobbing systems, pushing the boundaries of sustainable and high-precision gear manufacturing.