Multi-Objective Optimization of a Two-Stage Helical Gearbox Using Taguchi Method and Grey Relational Analysis

: This paper presents a novel approach to solve the multi-objective optimization problem designing a two-stage helical gearbox by applying the Taguchi method and the grey relation analysis (GRA). The objective of the study is to identify the optimal main design factors that maximize the gearbox efﬁciency and minimize the gearbox mass. To achieve that, ﬁve main design factors. including the coefﬁcients of wheel face width (CWFW) of the ﬁrst and the second stages, the allowable contact stresses (ACS) of the ﬁrst and the second stages, and the gear ratio of the ﬁrst stage were chosen. Additionally, two single objectives, including the maximum gearbox efﬁciency and minimum gearbox mass, were analyzed. In addition, the multi-objective optimization problem is solved through two phases: Phase 1 solves the single-objective optimization problem in order to close the gap between variable levels, and phase 2 solves the multi-objective optimization problem to determine the optimal main design factors. From the results of the study, optimum values of ﬁve main design parameters for designing a two-stage helical gearbox were ﬁrst introduced.


Introduction
Optimizing gearboxes is a critical aspect of mechanical engineering as it has a direct impact on their efficiency, durability, and reliability, as well as the performance of other machinery and equipment.Multi-objective optimization of gearboxes involves simultaneously optimizing various performance parameters, such as load-carrying capacity, noise, mass, size, and efficiency, making it a complex and challenging task.To address these challenges, various optimization techniques have been developed in recent years.Helical gearboxes, one of many gearbox types, are widely used in industrial applications due to their superior load-carrying capacity, smooth operation [1], simple structure, and low cost.However, designing a helical gearbox involves numerous design parameters, making it difficult to optimize for multiple objectives.
Numerous studies on the optimal design of helical gearboxes have been conducted thus far.This research looked into both single-objective and multi-objective optimization problems.Many authors are also interested in the single-objective optimization problem.I. Römhild and H. Linke [2] presented formulas for calculating gear ratios for two-, three-, and four-stage helical gearboxes in order to obtain the smallest gear mass.Milou et al. [3] presented a practical approach for reducing the mass of a two-stage helical gearbox in their paper.The method entailed analyzing data from gearbox manufacturers.Their findings suggested a center distance ratio (aw 2 /aw 1 ) between the second and first stages in the range of 1.4 to 1.6 to achieve the minimum mass.Once the optimal center distance ratio has been determined, the corresponding partial gear ratios are obtained using a lookup table.Various objective functions are utilized to solve the single-objective optimization problem of determining the optimal gear ratio of helical gearboxes.These objective functions include minimizing gearbox length [4][5][6][7], minimizing gearbox cross section area [6,[8][9][10], minimizing gearbox mass [6,11], minimizing gearbox volume [12], and minimizing gearbox cost [13][14][15][16].All of these studies share a common feature in that only one design parameter, the partial gear ratio, is defined in the form of an explicit model.
The multi-objective optimization problem for designing helical gearboxes has also been of interest to researchers.M. Patil et al. [1] conducted a study in which a twostage helical gearbox was subjected to multi-objective optimization with a broad range of constraints using a specially formulated discrete version of non-dominated sorting genetical algorithm II (NSGA-II).The study formulated two objective functions, namely the minimum gearbox volume and minimum total gearbox power loss.Moreover, the study considered constraints, such as bending stress, pitting stress, and tribological factors.Wu, Y.-R.and V.-T.Tran [17] introduced a new microgeometry modification for helical gear pairs, leading to substantial enhancements in performance with regards to noise and vibration.In their research, C. Gologlu and M. Zeyveli [18] utilized a genetic algorithm (GA) to optimize the volume of a two-stage helical gearbox.To handle the design constraints, such as bending stress, contact stress, number of teeth on pinion and gear, module, and face width of gear, the objective function was subject to static and dynamic penalty functions.The results from the GA were compared to those of a deterministic design procedure, with the GA being found to be the superior method.D. F. Thompson and colleagues [19] presented a generalized optimization technique for reducing the volume of two-stage and three-stage helical gearboxes, while considering tradeoffs with surface fatigue life.In their study, Edmund S. Maputi and Rajesh Arora [20] explored multi-objective optimization by simultaneously considering three objectives: volume, power output, and center-distance.They employed the NSGA-II evolutionary algorithm to generate Pareto frontiers in their research.From the results of the study, insights for the design of compact gearboxes can be gained.The NSGA-II method was also applied by C. Sanghvi et al. [21] to solve the multi-objective optimization problem of a two-stage helical gearbox for minimum volume and maximum load.The results of the multi-objective optimization of the tooth surface in helical gears using the response surface method were presented by Park C.I. [22].
The best or optimal level of the selected criterion is determined by single-objective optimization, which is an absolute optimization.A multi-objective optimization problem is one with two or more simple goals (or criteria).As a result, the solution to the multiobjective optimization problem cannot best satisfy all criteria simultaneously.For instance, it is not possible to meet both the efficiency and cost requirements of the gearbox.In simple terms, determining a solution to a problem that is both "white" and "black" is impossible; only a "gray" solution can be determined.The gray solution is the one that falls between the best and worst solutions, or between "white" and "black" in the multiobjective optimization problem.As a result, it is known as optimization based on gray relation analysis.The original Taguchi method is used to solve the single-objective problem, and the Taguchi method and gray relation analysis are required to solve the multi-objective optimization problem.
While numerous studies have focused on multi-objective optimization for helical gearboxes, the identification of optimal main design factors for such gearboxes has not received adequate attention.Furthermore, previous research on multi-objective optimization for helical gearboxes has not demonstrated the relationship between optimal input factors and total gearbox ratio.This is a critical issue to consider when designing a new gearbox.In this paper, we present a multi-objective optimization study for a two-stage helical gearbox, considering two single objectives: minimizing gearbox mass and maximizing gearbox efficiency.The proposal of five optimal main design factors for the two-stage helical gearbox is the most significant result of this research.These variables include the CWFW for both stages, the ACS for both stages, and the first stage's gear ratio.Furthermore, by combining the Taguchi method and the GRA in a two-stage process not previously described, we present a novel approach to addressing the multi-objective optimization problem in gearbox design.Additionally, a link between optimal input factors and the total gearbox ratio was proposed.

Gearbox Mass Calculation
The mass of the gearbox m gb can be found by the following equation: where m g , m gh , m b , and m s designate the mass of the gears, the shafts, the bearings, and the gearbox housing, respectively.The component mass will be specifically calculated below:

Gearbox Housing Mass Calculation
The gearbox housing mass (m gh ) can be found by: in which ρ gh is the weight density of gearbox housing materials referred; as the material of the gearbox housing is gray cast iron (the most common material for the gearbox housing), ρ gh = 7300 kg/m 3 ; V gh is the gearbox housing volume (m 3 ), which is determined by the following equation (see Figure 1): Appl.Sci.2023, 13, x FOR PEER REVIEW 3 of 18

Gearbox Mass Calculation
The mass of the gearbox   can be found by the following equation: where   ,  ℎ ,   , and   designate the mass of the gears, the shafts, the bearings, and the gearbox housing, respectively.The component mass will be specifically calculated below:

Gearbox Housing Mass Calculation
The gearbox housing mass ( ℎ ) can be found by: in which  ℎ is the weight density of gearbox housing materials referred; as the material of the gearbox housing is gray cast iron (the most common material for the gearbox housing),  ℎ = 7300 (/ 3 );  ℎ is the gearbox housing volume (m 3 ), which is determined by the following equation (see Figure 1): where   ,   , and   are the volumes of side A, B, and C (m 3 ), respectively.
Substituting Equations ( 4) to ( 6) into (3) obtains: In the above equations, L, H, B1, and SG can be calculated by [2]: where V A , V B , and V C are the volumes of side A, B, and C (m 3 ), respectively.

Gear Mass Calculation
The gear mass can be determined by: in which m g1 and m g2 are the gear mass of the first and the second stages (kg), which can be calculated by: where ρ g is the weight density of gear material (kg/m 3 ); as the gear material is steel, ρ g = 7800 kg/m 3 ; e 1 and e 2 are the volume coefficients.Because the pinion has a small diameter, its structural form can be plain whereas the gear has a large diameter and thus requires a hub.As a result, the volume coefficient of the pinion e 1 = 1, and the volume coefficient of the gear is e 2 = 0.6 [13]; b w1 and b w2 are the gear width of the first and the second stages (mm); d w1i and d w2i are the pinion and the gear pitch diameters of the i stage (i = 1 and 2).These parameters are determined by: In the above equations, u 1 is the gear ratio of the first stage; the center distance of i stage a wi is found by the surface fatigue strength [23]: where k a is the material coefficient; k Hβ is the contacting load ratio for pitting resistance; AS i is the allowable contact stress of the i th stage (MPa); X bai is the wheel face width coefficient of ith stage; and T 1i is the drive shaft torque of the ith stage (N.mm), which can be calculated by: After calculating the gear parameters, the bending strength of the i th gear stage must be checked using the following formulas [23]: in which m i is the module of the i th gear stage (mm); K Fi is load factor; Y εi = 1/ε is the contact factor; ε is the contact ratio; Y βi = 1 − β/140 is factor taking into account the helix angle; and Y F1i and Y F2i are the geometry factor of the pinion and gear of i th gear stage.

Calculation of Shaft Mass
The shaft mass can be found by: where m si is the mass of i th shaft of the gearbox (kg), which is determined by: in which d j and l j are the diameter and the length of j th shaft part (mm); d bk and B k are the diameter and the width of k th shaft part on which the bearing is istalled.The values of d j and d bk also are determined by [23]: wherein [σ s ] is the allowable shaft stress (MPa), which can be determined by the material and size of the shaft [23].M e is the equivalent moment (Nmm), which is found by [23]: where M x and M y are the bending moment in x and y directions (Nmm); T is the torque (Nmm).These parameters can be defined based on the diagram for finding shaft dimensions.Figure 2 describes this diagram for calculation of the first shaft of the gearbox.
In Equation ( 22), B k is the bearing width (mm).In this work, radical ball bearings with angular contact were used.From the data in [24], the following regression was proposed to calculate the width of the bearings (with R 2 = 0.9951): and size of the shaft [23].  is the equivalent moment (Nmm), which is found by [23]: where   and   are the bending moment in x and y directions (Nmm); T is the torque (Nmm).These parameters can be defined based on the diagram for finding shaft dimensions.Figure 2 describes this diagram for calculation of the first shaft of the gearbox.

Calculation of Bearing Mass
The bearing mass of the gearbox is calculated by: As mentioned above, radical ball bearings with angular contact were used in this work.From the data about these type of bearings in [24], the mass of the i th bearing can be found by the following regression equation (with R 2 = 0.9833): in which i is the number of bearings (i = 6); d bi is the diameter of the shaft part on which the i th bearing is mounted.

Determination of Gearbox Efficiency
The gearbox efficiency is determined by: where P l is the total power loss in the gear box [25]: where P lg is the power losses in the gears; P lb is the power loss in the bearings; and P ls is the power loss in seals.These factors can be determined as follows: (+) The power losses in the gears: Appl.Sci.2023, 13, 7601 7 of 19 in which P lgi is the the power losses in the gears of the i stage, which is found by: where η gi is the efficieny of the i stage of the gearbox, which can be determined by [26]: where u i is the gear ratio of i stage; f is the friction coefficient; β ai and β ri are the arc of the approach and recess on the i stage, which is calculated by [26]: in which R e1i and R e2i are the outside radius of the pinion and gear, respectively; R 1i and R 2i are the pitch radius of the pinion and gear, respectively; R 01i and R 01i are the base-circle radius of the pinion and gear, respectively; α is the pressure angle.
From the data in [26], the friction coefficient can be determined by the folowing regression equations: -When the sliding velocity is v ≤ 0.424 (m/s), the friction coefficient is calculated by (with R 2 = 0.9958): f = −0.0877•v+ 0.0525 (37) -When the sliding velocity is v > 0.424 (m/s), the friction coefficient is calculated by (with R 2 = 0.9796): (+) The power losses in the bearings [25]: The power losses in rolling bearings can be found by: where f b is the coefficient of the friction of the bearing; as the radial ball bearings with angular contact were used, f b = 0.0011 [25]; F denotes the bearing load in Newtons (N) while v represents the peripheral speed.Additionally, i represents the ordinal number of the bearing, ranging from 1 to 6.
(+) The total power losses in the seals are determined by [25]: in which i is the ordinal number of the seal (i = 1 ÷ 2); P si represents the power loss caused by the sealing for a single seal (w), which can be calculated by: where VG 40 is the ISO viscosity grades number.

Objectives Functions
The multi-objective optimization problem in this study comprises two single objectives: -Minimizing the gearbox mass: -Maximizing the gearbox efficiency: in which X is the design variable vector reflecting variables.In this work, five main design factors, including u 1 , Xba 1 , Xba 2 , AS 1 , and AS 2 were selected as variables, and we have: For a helical gear set, the maximal gear ratio is 9 [23].Additionally, the coefficient of the wheel face width of both gear stages of a two-stage helical gearbox ranges from 0.25 to 0.4 [23].In addition, the gear materials used in this work are steel 40, 45, 40X, and 35XM refining, with teeth hardness on the surface of 350 HB (These are the most commonly used gear materials in gearboxes).As a result of the calculated results, the allowable contact stresses of the first and second stages range from 350 to 420 (MPa).Therefore, the following constraints were derived from these comments:

Methodology
As stated in Section 2.3.1 five main design factors were selected as variables for the multi-objective optimization problem.Table 1 describes these factors and the min and max values of them.In this work, the Taguchi method and grey relation analysis were employed to address the multi-objective optimization problem with five variables.In order to easily determine the solutions of the optimization problem, the larger the number of levels of the variables, the better.To maximize the number of levels for each variable, the L25 (5 5 ) design was selected.However, among the variables mentioned, u 1 has a very wide distribution (u1 ranges from 1 to 9 as stated in Section 2.3.2).As a result, the gap between levels of this variable remained significant even with five levels (in this case, the gap is ((9 − 1)/4 = 2)).To reduce this gap, save time, and improve the accuracy of the solutions, a procedure for solving a multi-objective problem was proposed (Figure 3).This procedure consists of two phases: Phase 1 solves the single-objective optimization problem to close the gap between levels, and phase 2 solves the multi-objective optimization problem to determine the optimal main design factors.consists of two phases: Phase 1 solves the single-objective optimization pr the gap between levels, and phase 2 solves the multi-objective optimizatio determine the optimal main design factors.

Single-Objective Optimization
In this work, the direct search method is used to solve the single-objective optimization problem.Additionally, a computer program has been built using the Matlab language to solve two single-objective problems, including minimizing the gearbox mass and maximizing the gearbox efficiency.From the results of this program, the relation between the optimal value of the gear ratio of the first stage u 1 and the total gearbox ratio u t is shown in Figure 4. Additionally, new constraints for the variable u 1 were found, as shown in Table 2.
In this work, the direct search method is used to solve the single-objective optimization problem.Additionally, a computer program has been built using the Matlab language to solve two single-objective problems, including minimizing the gearbox mass and maximizing the gearbox efficiency.From the results of this program, the relation between the optimal value of the gear ratio of the first stage u1 and the total gearbox ratio ut is shown in Figure 4. Additionally, new constraints for the variable u1 were found, as shown in Table 2.

Multi-Objective Optimization
The multi-objective optimization problem in this research aims to identify the optimal main design factors that satisfy two single-objective functions: minimizing the maximum optimization and maximizing gearbox efficiency in the design of a two-stage helical gearbox with a specific total gearbox ratio.To address this problem, a simulation experiment was conducted.The experiment was designed using the Taguchi method, and the analysis of the results was performed using Minitab R18 software.In addition, as noted above, the design L25 (5 5 ) was chosen for obtaining maximal levels of the variable.A computer program has been developed to perform these experiments.An investigation was conducted to minimize programming intricacy by examining the influence of five key design parameters on gearbox mass.The input first pinion speed of 1480 (rpm) was selected as it is the most common.The steel 45 was selected as the shaft material as it is a very common shaft material.The total gearbox ratios considered for analysis were 10, 15, 20, 25, 30, and 35.Employing a five-level Taguchi design (L25), a total of 25 simulation experiments were carried out for each total gearbox ratio mentioned above.Table 3 describes the main design factors and their levels, and Table 4 presents the experimental plan and the corresponding output results, encompassing the gearbox mass and efficiency, specifically for the total gearbox ratio of 15.

Multi-Objective Optimization
The multi-objective optimization problem in this research aims to identify the optimal main design factors that satisfy two single-objective functions: minimizing the maximum optimization and maximizing gearbox efficiency in the design of a two-stage helical gearbox with a specific total gearbox ratio.To address this problem, a simulation experiment was conducted.The experiment was designed using the Taguchi method, and the analysis of the results was performed using Minitab R18 software.In addition, as noted above, the design L25 (5 5 ) was chosen for obtaining maximal levels of the variable.A computer program has been developed to perform these experiments.An investigation was conducted to minimize programming intricacy by examining the influence of five key design parameters on gearbox mass.The input first pinion speed of 1480 (rpm) was selected as it is the most common.The steel 45 was selected as the shaft material as it is a very common shaft material.The total gearbox ratios considered for analysis were 10, 15, 20, 25, 30, and 35.Employing a five-level Taguchi design (L25), a total of 25 simulation experiments were carried out for each total gearbox ratio mentioned above.Table 3 describes the main design factors and their levels, and Table 4 presents the experimental plan and the corresponding output results, encompassing the gearbox mass and efficiency, specifically for the total gearbox ratio of 15.The multi-optimization optimization problem is solved by applying the Taguchi and GRA methods.The main steps for this process are as follows: (+) Determining the signal to noise ratio (S/N) by the following equations as the object of this work is to reduce the gearbox mass and to increase the gearbox efficiency: -For the gearbox mass objective, the-smaller-is-the-better S/N: -For the gearbox efficiency objective, the-larger-is-the-better S/N: where y i is the ouput response value; m is number of experimental repetitions.In this case, m = 1 because the experiment is a simulation; no repetition is required.The calculated S/N indexes for the two mentioned output targets are presented in Table 5.In fact, the data of the two considered single-objective functions have different dimensions.To ensure comparability, it is essential to normalize the data, bringing them to a standardized scale.The data normalization is performed using the normalization value Zij, which ranges from 0 to 1.This value is determined using the following formula: In the formula, "n" represents the experimental number, which in this case is 25.
(+) Calculating the grey relational coefficient: The grey relational coefficient is calculated by: with i = 1, 2, . . ., n.In the formula, "k" represents the number of objective targets, which is 2 in this case; ∆ j (k) is the absolute value, ∆ j (k) = Z 0 (k) − Z j (k) width; Z 0 (k) and Z j (k)  are the reference and the specific comparison sequences, respectively; ∆ min and ∆ max are the min and max values of ∆i(k); ζ is the characteristic coefficient, 0 ≤ ζ ≤ 1.In this work ζ = 0.5.
(+) Calculating the mean of the grey relational coefficient: The degree of grey relation is determined by calculating the mean of the grey relational coefficients associated with the output objectives.
where y ij is the grey relation value of the j th output targets in the i th experiment.To ensure harmony among the output parameters, a higher average grey relation value is desirable.As a result, the objective function of the multi-objective problem can be transformed into a single-objective optimization problem, with the mean grey relation value serving as the output.
The impact of the main design factors on the average grey relation value (y) was analyzed using the ANOVA method, and the corresponding results are presented in Table 7. From the results in Table 7, AS 2 has the most influence on y (57.35%), followed by the influence of u 1 (25.16%),X ba1 (5.38%), X ba2 2.50%), and AS 1 (0.84%).The order of influence of the main design factors on y through ANOVA analysis is described in Table 8.Theoretically, the set of main design parameters with the levels that have the highest S/N values would be the rational (or optimal) parameter set.Therefore, the impact of the main design factors on the S/N ratio was determined (Figure 5).From Figure 3, the optimal levels and values of main design factors for multi-objective function were found (Table 9).The adequacy of the proposed model is assessed using the Anderson-Darling method, and the results are presented in Figure 6.From the graph, it is evident that the data points corresponding to the experimental observations (represented by blue dots) fall within the region bounded by upper and lower limits with a 95% standard deviation.Furthermore, the p-value of 0.226 significantly exceeds the significance level α = 0.05.These findings indicate that the empirical model employed in this study is appropriate and suitable for the analysis.Continuing from the previous discussion, the optimal values for the main design parameters corresponding to the remaining ut values of 10, 20, 25, 30, and 35 are presented in Table 10.Similarly, the optimal values for AS 1 and AS 2 are also their maximum values.This is because minimizing the gearbox mass requires maximizing the values of AS and AS 2 .By increasing these values, the center distance of gear stage i (as represented by Equation ( 19)) can be minimized.It will lead to the gear widths (as determined by Equations ( 15) and ( 16)), the pinion, and the gear pitch diameters of the ith stage (i = 1 and 2) (as calculated by Equations ( 17) and ( 18)), and therefore, the gear mass (as represented by Equations ( 13) and ( 14)) can be minimized.
Figure 7 depicts an obvious first-order relationship among the optimal values of u 1 and u t .Additionally, the following regression equation (with R 2 = 0.9992) to find the optimal values of u 1 was found: After finding u 1 , the optimum value of u 2 can be determined by u 2 = u t /u 1 .
To assess the effectiveness of the proposed method, the multi-objective optimization problem was solved using the constraints of u 1 , as shown in Table 1 (referred to as the solution by the traditional method).The optimal values for the main design parameters discovered by this method are shown in Table 11. Figure 8 depicts the optimal values of u 1 as determined by the traditional method (data from Table 10) and the new method (data from Table 9).This figure clearly shows that the optimal values of u 1 for the new method are easily determined and obey a very simple first-order function (Equation ( 52)).Furthermore, when determined by traditional methods, these values are distributed randomly rather than according to common rules (Figure 8), and they will almost certainly be less accurate than when determined by the new method.(19)) can be minimized.It will lead to the gear widths (as determined by Equations ( 15) and ( 16)), the pinion, and the gear pitch diameters of the ith stage (i = 1 and 2) (as calculated by Equations ( 17) and ( 18)), and therefore, the gear mass (as represented by Equations ( 13) and ( 14)) can be minimized.Figure 7 depicts an obvious first-order relationship among the optimal values of u1 and ut.Additionally, the following regression equation (with R 2 = 0.9992) to find the optimal values of u1 was found: To assess the effectiveness of the proposed method, the multi-objective optimization problem was solved using the constraints of u1, as shown in Table 1 (referred to as the solution by the traditional method).The optimal values for the main design parameters discovered by this method are shown in Table 11. Figure 8 depicts the optimal values of u1 as determined by the traditional method (data from Table 10) and the new method (data from Table 9).This figure clearly shows that the optimal values of u1 for the new method are easily determined and obey a very simple first-order function (Equation ( 52)).Furthermore, when determined by traditional methods, these values are distributed randomly rather than according to common rules (Figure 8), and they will almost certainly be less accurate than when determined by the new method.

Conclusions
The Taguchi method and the GRA are used in this paper to solve the multi-objective optimization problem in designing a two-stage helical gearbox.The goal of the study is to discover the optimal main design parameters that maximize gearbox efficiency while minimizing gearbox mass.To accomplish this, five major design factors were chosen: the CWFW for the first and second stages, the ACS for the first and second stages, and the first-stage gear ratio.In addition, the multi-objective optimization problem is solved in two stages.Phase 1 is concerned with solving the single-objective optimization problem of closing the gap between variable levels, while Phase 2 is concerned with determining the optimal main design factors.The following conclusions were proposed as a result of this work: -A novel approach to handling the multi-objective optimization problem in the gearbox design was presented by combining the Taguchi method and the GRA in a twostage process.The distance between the values of the lower and upper bounds of the constants of u1 is shortened as a result of this approach, which leads to an easier and a more accurate determination of optimal values.- The solution of the single-objective optimization problem bridges the gap between variable levels, making the solution of the multi-objective optimization problem easier and more accurate.-From the results of the study, optimal values for the five main design factors in the design of a two-stage helical gear gearbox were proposed (Equation (52) and Table 10).- The effect of the main design parameters on y ̅ was analyzed using the ANOVA method.The results revealed that AS2 had the highest influence on y ̅ (57.35%), followed by u1 (25.16%),Xba1 (5.38%), Xba2 (2.50%), and AS1 (0.84%).- The proposed model of u1 demonstrates a high level of consistency with the experimental data, validating their reliability.This model can be effectively utilized for multi-objective optimization of a two-stage helical gearbox, providing a valuable ap-

Conclusions
The Taguchi method and the GRA are used in this paper to solve the multi-objective optimization problem in designing a two-stage helical gearbox.The goal of the study is to discover the optimal main design parameters that maximize gearbox efficiency while minimizing gearbox mass.To accomplish this, five major design factors were chosen: the CWFW for the first and second stages, the ACS for the first and second stages, and the first-stage gear ratio.In addition, the multi-objective optimization problem is solved in two stages.Phase 1 is concerned with solving the single-objective optimization problem of closing the gap between variable levels, while Phase 2 is concerned with determining the optimal main design factors.The following conclusions were proposed as a result of this work: -A novel approach to handling the multi-objective optimization problem in the gearbox design was presented by combining the Taguchi method and the GRA in a two-stage process.The distance between the values of the lower and upper bounds of the constants of u 1 is shortened as a result of this approach, which leads to an easier and a more accurate determination of optimal values.- The solution of the single-objective optimization problem bridges the gap between variable levels, making the solution of the multi-objective optimization problem easier and more accurate.-From the results of the study, optimal values for the five main design factors in the design of a two-stage helical gear gearbox were proposed (Equation (52) and Table 10).- The effect of the main design parameters on y was analyzed using the ANOVA method.The results revealed that AS 2 had the highest influence on y (57.35%), followed by u 1 (25.16%),X ba1 (5.38%), X ba2 (2.50%), and AS 1 (0.84%).

Figure 2 .
Figure 2. Diagram for determining shaft dimensions.Figure 2. Diagram for determining shaft dimensions.

Figure 2 .
Figure 2. Diagram for determining shaft dimensions.Figure 2. Diagram for determining shaft dimensions.

Figure 3 .
Figure 3.The procedure for solving a multi-objective problem.

Figure 3 .
Figure 3.The procedure for solving a multi-objective problem.

Figure 4 .
Figure 4. Gear ratio of the first stage versus total gearbox ratio.

Figure 4 .
Figure 4. Gear ratio of the first stage versus total gearbox ratio.

Figure 7 .
Figure 7. Optimal gear ratio of the first stage versus total gearbox ratio.

Figure 7 . 18 Figure 8 .
Figure 7. Optimal gear ratio of the first stage versus total gearbox ratio.Appl.Sci.2023, 13, x FOR PEER REVIEW 16 of 18

Figure 8 .
Figure 8. Optimal gear ratio of first stage of traditional and new method.

Table 2 .
New constraints of u 1 .

Table 3 .
Main design factors and their levels for u t = 15.

Table 4 .
Experimental plans and output responses for u t = 15.

Table 5 .
Values of S/N of each experimental run of u t = 15.
Table 6 displays the calculated results of the grey relation value (y i ) and the average grey relation value of all experiments.

Table 6 .
Values of ∆ i (k) and y i .

Table 7 .
Factor effect on y.

Table 8 .
Order of main design factor effect on y.

Table 9 .
Optimal levels and values of main design factors.

Table 10 .
Optimum main design factors finding by traditional method.

Table 11 .
Optimum main design factors finding by new method.

Table 11 .
Optimum main design factors finding by new method.