Optimization Design for the Planetary Gear Train of an Electric Vehicle under Uncertainties

: The planetary gear train is often used as the main device for decelerating and increasing the torque of the drive motor of electric vehicles. Considering the lightweight requirement and existing uncertainty in structural design, a multi-objective uncertainty optimization design (MUOD) framework is developed for the planetary gear train of the electric vehicle in this study. The volume and transmission efﬁciency of the planetary gear train are taken into consideration as optimization objectives. The manufacturing size, material, and load input of the planetary gear train are considered as uncertainties. An approximate direct decoupling model, based on subinterval Taylor expansion, is applied to evaluate the propagation of uncertainties. To improve the convergence ability of the multi-objective evolutionary algorithm, the improved non-dominated sorting genetic algorithm II (NSGA-II) is designed by using chaotic and adaptive strategies. The improved NSGA-II has better convergence efﬁciency than classical NSGA-II and multi-objective particle swarm optimization (MOPSO). In addition, the multi-criteria decision making (MCDM) method is applied to choose the most satisfactory solution in Pareto sets from the multi-objective evolutionary algorithm. Compared with the multi-objective deterministic optimization design (MDOD), the proposed MUOD framework has better reliability than MDOD under different uncertainty cases. This MUOD method enables further guidance pertaining to the uncertainty optimization design of transportation equipment, containing gear reduction mechanisms, in order to reduce the failure risk. coefﬁcient is less than 0.8, and the bending stress is less than 0.8 when the uncertainty range is the largest (Case 3). From another perspective, the bending stress constraint and the tooth width coefﬁcient constraint are the most prone to failure types. In the three uncertainty cases, MUOD meets the reliability requirements for all uncertainty constraints. A higher λ indicates that the farther the optimization result is from the boundary of the inequality constraint, the higher its reliability. Compared with the conventional MDOD, the MUOD proposed in this study can design a more reliable planetary gear train and reduce the risk of constraint failure. The results show that MUOD sacriﬁces certain performance, but it is more reliable than MDOD.


Introduction
In recent years, electric vehicle technology has developed rapidly [1,2]. The planetary gear reducer is used in electric vehicles due to its high transmission efficiency and compact structure. Due to the space limitation of electric vehicles, the design of compact planetary gear trains has become a key issue. Numerous optimization methods are involved in the gear train design. For example, Parmar et al. [3] proposed a novel multi-objective optimization method, for planetary gear trains, using NSGA-II. Miler et al. [4] chose transmission volume and power loss as design objectives, and they optimized the parameters of the planetary gear train with multi-objective optimization. Sedak et al. [5] proposed a constrained multi-objective nonlinear optimization problem for planetary gearboxes, based on a hybrid element heuristic algorithm, considering gear volume, center distance, contact ratio, and power loss as optimization objectives. Patil et al. [6] proposed a multi-objective optimization strategy to minimize the total volume and power loss of the two-stage helical gearbox and spur gearbox. Compared to the single-objective optimization method with tribological constraints, the multi-objective optimization results in less power loss. Savsani et al. [7] used the particle swarm optimization algorithm, and the simulated annealing algorithm, to carry out the optimization design of the lightweight spur gear transmission system, and, resultingly, this method is deemed to be suitable for the single-objective or multi-objective optimization design of the multi-stage spur gear transmission. Considering the above research, the main challenge of gear transmission design is in reducing weight and power loss. At present, the optimization design method of planetary gear trains mainly considers the determining system parameters and implements the conventional deterministic optimization method. However, for practical engineering structures, many uncertainties are observed in the material properties, manufacturing, and measurement [8][9][10][11][12]. To obtain a reliable structural design, the uncertainties of the planetary gear train of electric vehicles need to be considered.
Uncertainty optimization in engineering design has gradually attracted attention [13][14][15]. For example, Xian et al. [16] proposed an effective analysis framework for stochastic optimization pertaining to non-linear viscous dampers of energy dissipation structures, which was applied to the uncertainty optimization of non-linear viscous dampers of suspension bridges. Lü et al. [17] proposed an efficient approach for the optimization design of dual uncertain structures, taking into account the dual robust design and the possibility of failure, quickly estimating the dual uncertain target of fuzzy random variables, and equivalently solving the possibility constraints involving fuzzy randomness. Baek et al. [18] developed a design method of a composite microwave absorbing structure using reliability-based optimization (RBO), which considers the failure probability. Compared with the results of deterministic optimization (DO), it was found that the total thickness of the reliability design method increased slightly, but RBO significantly reduced the failure probability. Fang et al. [19] developed an effective multi-objective uncertainty optimization program in order to design car doors. The program analyzed the impact of changing the uncertainty conditions and improving the reliability level, and it provided clear design information for decision-makers. Zhang et al. [20] proposed a reliable uncertainty optimization design route for obtaining optimal energy-absorbing structures. The study found that the solution obtained, by uncertainty optimization, sacrificed certain demand performance, but it was more reliable than deterministic design. The above studies have carried out the uncertainty optimization based on the probability model, which is highly dependent on statistical data. Considering that the distribution of uncertainty requires a lot of data, it is of a high cost to obtain effective probability data from a practical engineering perspective.
To overcome the limitation, of uncertainty optimization, due to the lack of data, some interval uncertainty modes have been gradually developed and applied to engineering optimization [21,22]. The interval uncertainty model mainly focuses on the upper and lower boundaries of uncertainty values, which is easier to implement than the probabilistic uncertainty model. Inuiguchi et al. [23] proposed a linear multi-objective strategy based on maximum and minimum regret criteria to solve the problem of interval uncertainty in the objective function. Fu et al. [24] developed a multi-objective direct structural optimization method for solving interval uncertainty. This method uses the satisfaction value of the interval possibility model to deal with non-linear uncertain constraints, and it judges the feasibility and infeasibility of individual design vectors. Wu et al. [25] proposed a nonprobabilistic robust topology optimization method for interval uncertain structures. The method uses the Chebyshev interval inclusion function to realize the non-invasiveness of the interval algorithm. Wang et al. [26] developed an effective interval uncertain optimization design strategy using Legendre polynomial chaotic expansion, which is more efficient than the conventional method. Hou et al. [27] carried out the uncertainty optimization, pertaining to the energy efficiency of ships in icy areas, considering the interval parameters; the optimization results provided practical guidance for the energy-saving design of ships in the case of uncertainty in the actual environment. Yu et al. [28] regarded friction coefficient, material properties, and wear element thickness as interval uncertainty factors, and proposed an uncertainty optimization method for the noise suppression of the brake system. The above studies have developed a highly effective uncertainty optimization method based on the interval model, and they have applied it to solve practical engineering problems. The interval model has been validated as a highly applicable uncertainty optimization method. Uncertainties in the manufacturing and operation of the planetary gear train of electric vehicles are unavoidable. The process of efficiently solving multi-objective uncertainty problems for the planetary gear train of electric vehicles is still a key issue. Therefore, a multi-objective uncertainty optimization design (MUOD) framework is developed for the planetary gear train of an electric vehicle in this study. Section 2 describes the detailed methodology of MUOD. Section 3 describes the design requirement of the planetary gear train of an electric vehicle. Section 4 shows the optimization results. The main conclusions are drawn in Section 5.

Multi-Objective Uncertainty Optimization Problem
In general, the multi-objective deterministic optimization design (MDOD) model can be expressed as follows [29,30]: In the formula, f 1 , . . . , f q are the objective functions and q is the number of objectives. G i (X) is the inequality constraint and l is the number of its constraints; h j (X) is the equality constraint, and g is the number of its constraints; and {S} is the design space. Different from the conventional deterministic optimization, the uncertainties of optimization variables and other relevant design parameters need to be considered during actual processing. Stochastic probability models are often used to construct uncertainty models, but the distribution information of uncertainties is unknown due to the lack of test samples. Therefore, the interval uncertainty model is employed in this study [31]. The multi-objective deterministic optimization can be transformed into the interval uncertainty problem, as follows: In the formula, x I and d I are interval design variables and other relevant design parameters, respectively. The superscripts IC and IR represent the nominal value and interval radius, respectively. The interval radius of an interval value reflects its fluctuation range and can be expressed as uncertainty deviation. When the design variables and other relevant parameters are interval values, the relationship of reliability-based possibility degree P d can be used to transform the interval uncertainty models into general noninterval models [31]. For the interval values A 1 , A 2 and A 1 ≤ A 2 , The superscripts IL and IU represent the lower and the upper values, respectively. The reliability-based possibility degree of the interval level should be given beforehand based on the actual reliable problem. Therefore, the multi-objective uncertainty optimization model can be expressed as: In the formula, λ i is the requirement of reliability-based possibility degree, and it also represents the equivalent reliability with different constraints. The main optimization goals and constraints have been described in Section 3.
Nested optimization design is often used in interval uncertainty optimization, which treats the uncertainty analysis problem as an internal optimization problem. The purpose of inner optimization is to evaluate the propagation of uncertainty and feed it back to the outer optimization route. It is worth considering that adding a new optimization solver will cause low computational efficiency. Therefore, Taylor expansion, as an effective decoupling method, is applied to analyze the propagation of uncertainty in this study [32].
The constraint function G i x I , d I can be approximately constructed by first-order Taylor expansion, that is: Therefore, the lower and upper bounds of the constraint function can be expressed as follows: Generally, the Taylor formula can achieve the best approximation in the case of a small interval uncertainty. Further, the calculation accuracy can be improved by establishing a subinterval to compensate for the nonlinear approximation error. For the uncertainty values U, In the formula, U s and S_n are the s th subinterval and the subinterval number, respectively. The subinterval number can be determined by referring to the number of uncertain parameters. The interval range of constraint function G i (U) is expressed as follows: Through the above interval uncertainty analysis method, the uncertain information of constraint function G i (U) can be solved by using the approximate direct decoupling method.

Improved Evolutionary Algorithm
The classical non-dominated sorting genetic algorithm (NSGA-II) generally uses the random function to generate the initial population [33], and its population uniformity is poor. The crossover probability and mutation probability of classical NSGA-II are set to a fixed value, respectively, and the optimization algorithm falls into the premature problem. Therefore, this paper adopts the improved NSGA-II designed by using chaotic and adaptive evolutionary strategies in order to obtain the multi-objective solution set.
Here, a chaotic strategy is used to generate the initial population of a multi-objective evolutionary algorithm, which can improve the diversity of the population. Tent map is one of the most commonly used mapping functions for generating chaotic sequences [34]. Here, the main steps of population chaos initialization and assignment, using the Tent mapping method, are as follows: Step 1: Randomly generate an N-dimensional random number vector, 1], where N v is the number of optimization variables.
Step 2: The improved Tent mapping method is used to calculate the chaotic component of each optimized variable, as follows: In the formula, i = 1, 2, . . . , P s , and P s is the population size; j = 1, 2, . . . , N v .
Step 3: Substitute each chaotic component obtained in Step 2 into the real range of each optimization variable, as follows: In the formula, x j lower and x j upper are the lower and upper bounds of the j th optimized variable respectively.
Here, the adaptive evolutionary strategy mainly improves the crossover and mutation operators. The adaptive crossover probability and mutation probability are generated according to the number of iterations, which is helpful to accelerate the convergence of optimization. In this study, the exponential function is applied to the adaptive adjustment mode of crossover probability and mutation probability. The calculation formula is described as follows: In the formula, pc(n i ) and pm(n i ) are the crossover probability and mutation probability at the n i th iteration; pc(0) and pm(0) are the initial crossover probability and mutation probability respectively; n t is the total evolutionary generation.

Multi-Criteria Decision Making (MCDM) Method
Usually, the Pareto solution set in multi-objective optimization can provide decisionmakers with numerous feasible design schemes at the early stage of design, but it cannot directly obtain the most satisfactory solution. In addition, the weight method aggregates multi-objective optimization into a single comprehensive objective to obtain the ideal optimal solution. However, although some decision-makers are full of engineering experience, it is nonetheless difficult to assign the optimal weight to each optimization objective. Therefore, as a multi-criteria decision making (MCDM) model, grey relational analysis (GRA) will be applied to select the most satisfactory scheme in Pareto sets [35][36][37]. Here, the GRA with entropy weight method is proposed to identify the most satisfactory solution. The normalization method can be adopted in the grey relation analysis, depending on the characteristics of the original sequence. When the target value of the original sequence is "the larger the better", the original sequence can be normalized as: In the formula, y * i (k) is a new sequence after normalization; max[y i (k)] is the maximum value of the original sequence; and min[y i (k)] is the minimum value of the original sequence. When the target value of the original sequence is "the smaller the better", the original sequence can be normalized as: After normalization, the grey relational coefficient γ i (k), which is used to quantify the relationship between the target and actual normalized results, can be formulated as [26]: In the formula, ∇ oi (k) is the deviation between reference sequence ∇ oi (k) and the compared sequence x * i (k), as follows: ρ is the distinguishing coefficient, ρ ∈ [0, 1], and ρ = 0.5 in this study. After obtaining the grey relational coefficient, the grey relational grade c i is presented in a weighted sum of the grey relational coefficients, as follows: For the actual engineering requirements, the effect of each criterion on the design objectives is not exactly the same; resultingly, Equation (25) can be modified to In the formula, w k is a weight of k th criterion. In this study, w j is determined by the entropy weight method. The weight is calculated by using the entropy weight method according to the variation degree of each criterion.
Different from the analytic hierarchy process (AHP) [38], the entropy weight method can objectively obtain the weight of each criterion according to the amount of information provided by each criterion and the correlation between the criteria, which overcomes the subjectivity in determining the weight of the criterion. Assuming O ik is the i th alternative value of the k th evaluation criterion, and the initial evaluation matrix is The proportion of the i th alternative value of the k th evaluation criterion is [39]: The entropy e k of the k th criterion is: When p ik is equal to 0, to ensure that ln(p ik ) is meaningful, Equation (23) can be modified to: Therefore, the entropy weight of the k th criterion can be expressed as follows:

Main Processes of MUOD
The main steps of the multi-objective uncertainty optimization design (MUOD) framework are as shown in Figure 1: Step 1: Define multi-objective optimization problems, optimization variables, objectives, and constraint functions. This step is similar to conventional deterministic multiobjective optimization.
Step 2: The interval optimization problem is transformed into a deterministic optimization problem using the relationship of reliability-based possibility degree. It should be noted that the interval uncertainty transformation mainly aims at the inequality constraints in the multi-objective optimization model. The lower and upper values of constraint functions can be solved directly by using the approximate direct decoupling method.
Step 3: The execution process of the multi-objective evolutionary algorithm. An improved multi-objective evolutionary algorithm is applied to calculate the transformed mathematical model in Step 2. The improved NSGA-II is designed in Section 2.2. The initial crossover probability and mutation probability are 0.8 and 0.1.
Step 4: The execution process of the MCDM method. The MCDM method is described in Section 2.3. The GRA with entropy weight method is applied to choose the most satisfactory solution in Pareto sets. The objective weight is calculated by using the entropy weight method, and the grey relational grade is calculated by Equation (26).

Main Design Variables and Optimization Objectives
The electric drive system and its planetary gear train are shown in Figure 2. The planetary gear train is used to reduce the speed and increase the output torque of the motor.

Main Design Variables and Optimization Objectives
The electric drive system and its planetary gear train are shown in Figure 2. The planetary gear train is used to reduce the speed and increase the output torque of the motor. Since the helical gear has the advantages of good meshing, stable transmission, and low noise, the helical planetary gear train is designed in this study. The main parameters of an electric commercial vehicle are shown in Table 1, which are provided by a vehicle company. The transmission ratio of the gear and the final output torque can be calculated according to the vehicle parameters. The planetary gear train should be compact; that is, the overall volume of the planetary gear train should be small enough to facilitate the arrangement of the electric drive system in the chassis of the electric vehicle. Smaller gear volume corresponds to lighter weight, which is more conducive to the improvement of energy efficiency. Therefore, the volume of the helical planetary gear train is used as the optimization objective function. To simplify the calculation, the volume of the ring gear is chosen as the design objective, as follows: In the formula, is the pitch circle diameter of the ring gear; is the tooth width. In general, the power loss in gear transmission mainly includes the friction loss caused by gear tooth surface meshing, the bearing loss, and the stirring loss of lubricating oil. The meshing friction power loss is the main reason for the gear transmission power loss. Therefore, the transmission efficiency, considering the meshing friction power loss, is regarded as the second design objective. The planetary gear train mainly includes external and internal meshing of the gear, as shown in Figure 3a   The planetary gear train should be compact; that is, the overall volume of the planetary gear train should be small enough to facilitate the arrangement of the electric drive system in the chassis of the electric vehicle. Smaller gear volume corresponds to lighter weight, which is more conducive to the improvement of energy efficiency. Therefore, the volume of the helical planetary gear train is used as the optimization objective function. To simplify the calculation, the volume v of the ring gear is chosen as the design objective, as follows: In the formula, d r is the pitch circle diameter of the ring gear; b is the tooth width. In general, the power loss in gear transmission mainly includes the friction loss caused by gear tooth surface meshing, the bearing loss, and the stirring loss of lubricating oil. The meshing friction power loss is the main reason for the gear transmission power loss. Therefore, the transmission efficiency, considering the meshing friction power loss, is regarded as the In the formula, is the pitch circle diameter of the ring gear; is the tooth width. In general, the power loss in gear transmission mainly includes the friction loss caused by gear tooth surface meshing, the bearing loss, and the stirring loss of lubricating oil. The meshing friction power loss is the main reason for the gear transmission power loss. Therefore, the transmission efficiency, considering the meshing friction power loss, is regarded as the second design objective. The planetary gear train mainly includes external and internal meshing of the gear, as shown in Figure 3a, and 3b. Here, the transmission efficiency of external meshing can be expressed as [40]: Here, the transmission efficiency η ex of external meshing can be expressed as [40]: In the formula, µ is the friction coefficient; Z s is the number of teeth of the sun gear; Z p is the number of teeth of the planet gear; α t is the transverse pressure angle; p t is the transverse circular pitch; R bp and R p are the base circle radius and pitch circle radius of planet gear respectively; and ε ei and ε eo are the meshing in and meshing out contact ratio of the gear external meshing, respectively, as follows: In the formula, α is the working pressure angle; α as and α ap are the tooth top pressure angle of sun gear and planet gear respectively. As shown in Figure 3c, assuming that the load distribution coefficient in the regions "a-b" and "c-d" is 0.5, f ex (ε) can be expressed as follows: Similarity, the transmission efficiency η in of internal meshing can be expressed as follows [40]: In the formula, Z r is the is the number of teeth of the ring gear; ε ii and ε io are the meshing in and meshing out contact ratio of the gear external meshing, respectively, as follows: Same as Equation (31), assuming that the load distribution coefficient is 0.5, f in (ε) can be expressed as follows: Therefore, the transmission efficiency η p of the planetary reduction gear train can be expressed as follows [41]: (36) In this study, six main parameters are considered as design variables to find the minimum volume and maximum transmission efficiency. These design variables X v include the teeth number of the sun gear Z s , the teeth number of the planet gear Z p , the teeth number of the ring gear Z r , helix angle β, face width b, and normal module m n . Z s , Z p , and Z r are integers, β and b are continuous, and m n is discrete. The alternative modulus m n is shown in Equation (38). Table 2 shows the detailed information of all design variables.

Main Design Constraints
The gear design should meet the specified constraints to meet the actual geometric, load, and material requirements. The main constraints of the planetary gear train in this study are as follows.

Equally Spaced Planets
To prevent the gear teeth from interfering with the mating gear, the gear teeth of all gears must mesh with the center gear teeth at the same time. The installation requirement needs to meet the following conditions: In the formula, n p is the number of planet gears.

Equally Spaced Planets
According to the actual power requirement of an electric vehicle, the transmission ratio r i of the planetary gear train needs to meet the following conditions:

Tooth Width Coefficient
The size of the tooth width is related to the strength of the gear; the larger the tooth width, the higher the strength. However, it should be noted that, if the tooth width is too large, there will be a larger number of tooth contact errors as well as a more uneven load distribution in the tooth direction. Therefore, it is critical to select an applicable tooth width. Here, the tooth width coefficient Φ d is the primary indicator of the tooth width design, which needs to meet the following constraints [42]:

Minimum Teeth of No-Undercut
Gear undercutting not only weakens the root of gear teeth while reducing the bending strength, but it also reduces the coincidence degree. Therefore, undercutting should be avoided in the gear design stage. The minimum number of teeth without undercutting of the helical cylindrical gear needs to meet the following constraints:

Concentric Constraint
The center distance between sun gear, ring gear, and planetary gear should be equal. The concentric constraint is:

Adjacency Constraint
To prevent the planet gears from colliding with each other, it is necessary to ensure that the planet gears have a certain clearance on their connecting lines; that is, the sum of the tooth top circle radius of two adjacent planetary gears shall be less than the center distance of two adjacent planetary gears. The adjacency constraint is: In the formula, d ap is the addendum circle diameter of planet gear; l sp is the center distance between the sun gear and the planet gear.

Contact Stress Requirement
The planetary gear train should be able to resist material failure (deformation and fracture) during contact behavior. The real contact stress requirement σ c needs to meet the following constraint [43]: In the formula, F n is the normal load; E is the elastic modulus of the material; [σ c ] is the allowable contact stress, and the gear material is 40Cr in this study; and ρ a and ρ p denote the radius of curvature at the nodes of the driving and driven gears, respectively. In the planetary gear system, the contact stress mainly occurs between the sun gear and the planetary gear, and between the planetary gear and the ring gear. However, considering that the dangerous position is usually present between the sun gear and the planetary gear, this study will focus on the contact stress between the sun gear and the planetary gear.

Bending Stress Requirement
The bending stress requirement σ w_s needs to meet the following constraint [42]: In the formula, T t is the transmitted torque; d is the pitch circle diameter; K σ is the stress concentration factor of contact ratio, K σ = 1.5; K ε is the influence factor of contact ratio, K ε = 2.0; and Y is the tooth profile coefficient. According to the empirical formula, the tooth profile coefficient of the sun gear can be expressed as follows: In the formula, Z v is a virtual number of teeth pertaining to the helical gear, Z v = Z/ cos 3 β.

Optimization Results and Discussions
In this study, the helix angle β, face width b, elastic modulus E, and input torque T are considered uncertain. The uncertain helix angle β and face width b are regarded as the uncertainty of manufacturing size. The uncertain elastic modulus E is regarded as the uncertainty of material. The uncertain input torque T is regarded as the uncertainty of load input. This study defines three uncertainty cases with different uncertainty deviations, which correspond to different degrees of uncertainty deviations, as shown in Table 3. Therefore, the constraint functions related to the above uncertain values can be regarded as uncertainty constraints. Uncertainty constraints mainly include tooth width coefficient, minimum teeth of no-undercut, adjacency constraint, contact stress constraint, and bending stress constraint. The nominal values of elastic modulus E IC and input load T IC t are 210 GPa and 3936 N·m, respectively. All requirements of reliability-based possibility degree λ are defined as 0.8. Here, the classical NSGA-II and multi-objective particle swarm optimization (MOPSO) are implemented in order to explore the feasibility of the improved NSGA-II. The initial population size is 400, the maximum number of iterations is 200, and the objective number of non-dominated solutions is 200. Figure 4 shows the iterative history of MDOD by using MOPSO, NSGA-II, and improved NSGA-II. The number of non-dominated solutions obtained by improved NSGA-II increases steadily, and improved NSGA-II can obtain nondominated solutions more efficiently than MOPSO and NSGA-II. Therefore, the improved NSGA-II designed in this paper is effective, and it contains better optimization potential than the classical NSGA-II and MOPSO. The improved NSGA-II will be implemented for MUOD. Figure 5 shows the optimal Pareto solution sets of MDOD and MUOD. There is an intense conflict between volume and transmission efficiency, which cannot achieve the common optimization; that is, the further improvement of one objective will inevitably worsen the other objective. The Pareto solution set of deterministic optimization design is lower than that of uncertainty optimization design, and the optimization objective of deterministic optimization design is better than uncertainty optimization design. In general, the inequality constraint of deterministic optimization is mainly concentrated near the constraint boundary, so its Pareto solution set has more loose space, and it is easier to obtain the better solution. It should be noted that, with the increase of uncertainty, the optimization results of MUOD tend to be conservative. To obtain the optimal solutions in different multi-objective optimization models, this study makes a trade-off analysis on the Pareto solution set by using the MCDM method. The optimal results of MDOD and MUOD are shown in Appendix A. Table A1 presents the optimization results of MDOD, and Tables A2-A4 present the optimization results of MUOD. The alternatives of different optimization methods are sorted according to the grey correlation degree, and the optimal solutions of all optimization methods are shown in bold. It is found that all moduli are 2, which means that the optimization space for modulus is small. Here, three uncertainty cases are substituted into all optimization results, and the obtained constraints are shown in Table 4. The optimal result of MDOD shows that the λ of the upper bound of tooth width coefficient is less than 0.8, and the bending stress is less than 0.8 when the uncertainty range is the largest (Case 3). From another perspective, the bending stress constraint and the tooth width coefficient constraint are the most prone to failure types. In the three uncertainty cases, MUOD meets the reliability requirements for all uncertainty constraints. A higher λ indicates that the farther the optimization result is from the boundary of the inequality constraint, the higher its reliability. Compared with the conventional MDOD, the MUOD proposed in this study can design a more reliable planetary gear train and reduce the risk of constraint failure. The results show that MUOD sacrifices certain performance, but it is more reliable than MDOD.
(N m) 393.6 590.4 Here, the classical NSGA-II and multi-objective particle swar (MOPSO) are implemented in order to explore the feasibility of the imp The initial population size is 400, the maximum number of iterations is 20 tive number of non-dominated solutions is 200. Figure 4 shows the ite MDOD by using MOPSO, NSGA-II, and improved NSGA-II. The numb nated solutions obtained by improved NSGA-II increases steadily, and im II can obtain non-dominated solutions more efficiently than MOPSO and fore, the improved NSGA-II designed in this paper is effective, and it co timization potential than the classical NSGA-II and MOPSO. The improv be implemented for MUOD.  Figure 5 shows the optimal Pareto solution sets of MDOD and MU intense conflict between volume and transmission efficiency, which can common optimization; that is, the further improvement of one objectiv worsen the other objective. The Pareto solution set of deterministic optim lower than that of uncertainty optimization design, and the optimization terministic optimization design is better than uncertainty optimization de the inequality constraint of deterministic optimization is mainly conce constraint boundary, so its Pareto solution set has more loose space, a obtain the better solution. It should be noted that, with the increase of ctuators 2022, 11,49 optimization results of MUOD tend to be conservative. To obtain the optimal so different multi-objective optimization models, this study makes a trade-off analy Pareto solution set by using the MCDM method. The optimal results of M MUOD are shown in Appendix A. Table A1 presents the optimization results o  and Table A2 ~ A4 present the optimization results of MUOD. The alternatives o optimization methods are sorted according to the grey correlation degree, and th solutions of all optimization methods are shown in bold. It is found that all mo which means that the optimization space for modulus is small. Here, three u cases are substituted into all optimization results, and the obtained constraints a in Table 4. The optimal result of MDOD shows that the of the upper boun width coefficient is less than 0.8, and the bending stress is less than 0.8 when tainty range is the largest (Case 3). From another perspective, the bending stress and the tooth width coefficient constraint are the most prone to failure types. In uncertainty cases, MUOD meets the reliability requirements for all uncerta straints. A higher indicates that the farther the optimization result is from the of the inequality constraint, the higher its reliability. Compared with the con MDOD, the MUOD proposed in this study can design a more reliable planetary and reduce the risk of constraint failure. The results show that MUOD sacrific performance, but it is more reliable than MDOD.

Conclusions
In order to design a reasonable planetary reduction gear system, matching the electric vehicle motor, this study proposes a multi-objective uncertainty optimization design (MUOD) framework for the planetary gear train of an electric vehicle. An approximate direct decoupling model, based on subinterval Taylor expansion, is applied to evaluate the propagation of uncertainties; the improved evolutionary algorithm is designed by using chaotic and adaptive evolutionary strategies. The volume and transmission efficiency of the planetary gear are optimization objectives. The optimization results of MUOD show that the Pareto front gradually moves to the upper right corner with the uncertainty increases. The most satisfactory solutions (improving lightweight and improving transmission efficiency) pertaining to different multi-objective optimization models can be obtained by the MCDM method. Compared with the conventional multi-objective deterministic optimization design (MDOD) method, the uncertainty optimization design of the planetary gear train sacrifices certain performance. When the reliability-based possibility degree λ is defined as 0.8, the optimization results of MUOD always meet this requirement, but at least one constraint violates this requirement in MDOD. As the degree of uncertainty increases, the optimization results of MUOD tend to be conservative, but MUOD is more reliable than MDOD. The uncertainty of planetary gear trains of electric vehicles is very complex in actual working conditions. The MUOD framework proposed in this study is able to continue carrying out optimization design with more complex high dimension uncertainties and objectives in order to ensure that the structure has better potential to resist the risk of failure.
Funding: This research was funded by "Project of Shanghai Science and Technology Committee, grant number 20511104602" and "Prospective Technology Project of Nanchang Intelligent New Energy Vehicle Research Institute, grant number 17092380013".

Conflicts of Interest:
The authors declare no conflict of interest.