Parameter Sensitivity Analysis of Mounting Pedestals and Multi-Objective Optimization for a Multi-Support Rigid Body System

In this paper, a parameter sensitivity analysis of mounting pedestals and a multi-objective optimization design for vibration reduction in a multi-support rigid body system, taking an aeroengine-lubricating oil tank supported by multiple mounting pedestals as an example, are conducted based on the third version of non-dominated sorting genetic algorithm (NSGA-Ⅲ) combined with Sobol’s sensitivity analysis method (SSAM). An aeroengine-lubricating oil tank with three mounting pedestals is simplified as a three-support dynamic system, and its dynamics model is established. Several structural parameters of mounting pedestals are taken as the design variables, and the system vibration response and the reaction force of the front and rear mounting pedestals are considered as the objective functions. The first-order results and total sensitivity index of different design parameters for each objective function are obtained via SSAM, and the five most sensitive parameters are selected. Based on the above five design parameters, multi-objective optimization designing for vibration reduction in a simplified lubricating oil tank system is conducted based on NSGA-Ⅲ, and the results of the above triple-objective optimization are obtained as a Pareto-front surface with an obvious frontier. It can be observed from the simulation results that the oil tank vibration of the optimized system is effectively suppressed under the unbalanced excitation of two typical engine speeds. The established method and the main results can provide guidance for designers of aeroengine external structural systems, which can help to achieve superior system dynamic performances in engineering applications.


Introduction
Some devices in industry engineering applications such as aerospace, energy and power, vehicle engineering, and the robotics field, are fixed by multi-point supports, and the composed systems have the typical characteristics of a rigid body system with spatial multi-point supports. Their dynamic characteristics are complex, and they are particularly prone to vibration, or even resonance, under external excitation [1][2][3]. For example, the lubricating oil tank of an aeroengine's external device is fixed on the casing through spatial multi-point anchored supporting structures (usually three mounting pedestals or three supports with a tie rod and two mounting pedestals, and some are with four mounting pedestals as well) at different locations on a typical large-scale rigid body system with spatial multi-point supports. The lubricating oil tank system with multi-point supports has multiple excitations, a wide frequency domain, and a complex vibration response, and it is prone to strong vibration or even damage under aeroengine complex vibration excitations. Owing to the increasing demand for raising the thrust-to-weight ratio of modern aeroengines, higher requirements are placed on the strength, vibration and stability of aeroengine-lubricating oil tanks' mounting pedestals. Therefore, in the design process of aeroengine-lubricating oil tanks and their mounting pedestals, it is of great need to carry out multi-objective optimization design for vibration reduction for a typical multi-support rigid body system.
Structural optimization is currently a frontier subject in engineering structural design, and it has a wide range of applications in engineering practice. It can be divided into size optimization, shape optimization, and topology optimization according to the type of design variables [4]. Due to the importance of structural optimization, the depth and breadth of research in this field have continued to expand, and many achievements have been obtained. Yang et al. [5] proposed a variety of stretcher structure optimization schemes based on the principle of dynamic vibration absorbers (DVAs). By optimizing the stretcher support stiffness and damping coefficient, the performance of the nonlinear vibration damping system of the tracked ambulance was improved, and vibration of the supine human body on the stretcher was reduced. The global sensitivity analysis method was used to calculate the sensitivity of the design parameters, and the non-dominated sorting genetic algorithm (NSGA-II) was used to optimize the nonlinear damping performance of the system. Caresta and Kessissoglou [6] used a resonant converter in the propulsion system to reduce the far-field radiated sound pressure from the submarine, and a combined genetic and pattern search algorithm was used to find the optimum virtual mass, stiffness, and damping parameters of the resonance converter. At present, most topology design studies are focused on single-objective topology optimization (SOTO), and some are also focused on multi-objective topology optimization (MOTO), but only focus on the effect of multi-stiffness on dynamic characteristics or vibration response [7,8]. Tan et al. [9] took stiffness and random vibration response as the optimization objectives, and applied multiobjective topology optimization to reduce the vibration of the primary supporting structure of the video satellite in the frequency domain. The dynamic analysis results showed that the root mean square values in all three spatial directions of the optimal structure by the multi-objective optimization were smaller than that of the single-objective optimization. The multi-objective optimization method significantly improved the vibration response of the base plate, which effectively suppressed the vibration of the satellite. Kapania et al. [10] summarized research on the structural optimization of stiffened plates, and pointed out that the optimization of the layout, thickness and height of the stiffeners was indispensable in practical applications. He et al. [11] adopt a data-driven method to study the multi-objective structural optimization problem, taking the layout, thickness, and height of the stiffeners as design variables to optimize the weight and mean square speed of stiffened plates. The numerical results showed that the proposed method was effective for the optimal design of stiffened plates and satisfied different band gap requirements in engineering. It can be observed that previous research on structural dynamic characteristics optimization was mostly focused on structure stiffness and damping parameters, ignoring structure geometry and mechanical parameters.
For sensitivity analysis (SA), a sensitivity analysis method that is very popular in many fields is the variance-based Sobol method. In general, variance-based sensitivity analysis methods aim to quantify the amount of variance that each parameter contributes to the unconditional variance of the model output [12,13]. Li et al. [14] proposed a novel sparse grid numerical integration method to estimate Sobol's main effect indices, as well as the two decomposed sensitivity measures, and applied the proposed method to a roof truss and a composite beam model to understand the impact of correlation between variables on the output variance explicitly, greatly alleviating the burden of the curse of dimensions on high-dimensional problems. Wu et al. [15] took the uncertain variables of link lengths and joint rotation angles into account, and used Sobol's method to analyze the sensitivity of uncertain variables for the positioning accuracy of industrial robots, by which the sensitive variables were determined to perform the positioning accuracy reliability analysis for industrial robots. Taking the bulb hydro generating unit (BHGU) as the research object, a nonlinear dynamic model of the BHGU with a shaft crack fault was established by Yan et al. [16]. The parameter sensitivity of the BHGU with the shaft crack fault was analyzed via Sobol's method, and ranks of parameter sensitivity for the vibration displacement and unbalanced forces were obtained, respectively.
For optimization algorithm, the non-dominated sorting genetic algorithm (NSGA-II) was proposed by Deb [17] as an effective multi-objective evolutionary algorithm, and has been widely adopted by many researchers. Liu and Jiao [18] presented a multi-objective pipe-routing algorithm considering the vibration performance of an aeroengine based on the Kriging model and NSGA-II, where pipe length, smoothness and natural frequency were taken into account. The proposed algorithm can obtain a set of non-dominant pipeline layout solutions; especially, the natural frequency of the pipeline was far away from the engine excitation frequency, avoiding possible resonance. Chen et al. [19] defined the RMS of weighted acceleration, wheel dynamic load, and suspension dynamic deflection as evaluation indices of vehicle ride comfort, and established a multi-body dynamic rigidflexible coupling model of an in-wheel motor vehicle. The relationship between vehicle ride comfort evaluation indices and suspension parameters was described by the Kriging model, and a multi-objective optimization using NSGA-II was processed using the above model, which improved the ride comfort performance. Zamanifar et al. [20] introduced an optimization technique based on NSGA-II, incorporating DFT analysis to achieve better control the performance of DFIG system stability. The NSGA-II was introduced and applied to solve the optimization problem of the propeller for high-altitude airships by Jiao et al. [21]. Aiming to design an efficient and lightweight propeller for high-altitude airships, the optimization model was established based on the conditions of high-altitude airship propulsion systems, and the effects of various design variables including pitch angle, chord length, diameter, rotational speed, and the number of blades on high-altitude propeller performance were presented and discussed. The third generation of the NSGA-III algorithm has been used for several years, and its performance is better than that of the second-generation algorithm. Based on human electroencephalogram (EEG) signals, Yi et al. [22] tested the performance of three crossover operators of the NSGA-III algorithm to solve a large-scale optimization problem The performance of the proposed NSGA-III variants was verified on six large-scale optimization problems. Aiming for a high-speed automotive turbocharger, Tarlani et al. [23] optimized the dynamic characteristics of the bearings, including the stiffness and damping coefficients, using the NSGA-III multiobjective optimization method to minimize vibrations caused by rotating unbalance. The results demonstrate that using optimized bearings can significantly improve the dynamic behavior of the rotor-bearing system, both at steady and transient operating conditions.
To summarize, the existing research on parameter sensitivity analysis and structure multi-objective topology optimization are separated. Sensitivity analysis can obtain ranks of parameter sensitivity, while subsequent optimization design is a direct step to solve engineering problems. Multi-objective topology optimization can obtain several design variables, while sensitivity analysis is the key step to select the most sensitive design variables. To the best of the author's knowledge, Sobol's sensitivity analysis method combined with the NSGA-III multi-objective optimization design for vibration reduction in an aeroengine-lubricating oil tank has never been reported. In this paper, an aeroenginelubricating oil tank supported by multiple mounting pedestals is firstly simplified as a typical three-support system, and its dynamics model is established. Then, several structural parameters of mounting pedestals are taken as design variables, and the lubricating oil tank vibration response as well as the reaction force of the front and rear mounting pedestals are considered as the objective functions. Furthermore, the first-order results and total sensitivity index of different design parameters for each objective function are obtained, and the five most sensitive parameters are selected. Based on above five design parameters, multi-objective optimization design for vibration reduction in a simplified lubricating oil tank system is conducted based on NSGA-III, and the results of the tripleobjective optimization are obtained with the Pareto-front surface with an obvious frontier. The simulation results show that the oil tank vibration of the optimized system is sup- pressed effectively under the unbalanced excitation of two typical engine speeds. The established method and main results can provide guidance for designers of aeroengine external structural systems, which can help to achieve superior system dynamic performances in engineering applications.

Dynamics Model of a Rigid Body System with Multi-Point Supports
The specific structural object of an aeroengine-lubricating oil tank supported by multiple mounting pedestals can be simplified as a rigid body system with multi-point supports. The lubricating oil-tank is a typical three-support rigid body system that is mainly composed of an oil tank body, a casing, two fixed mounting pedestals, and an articulated mounting pedestal (tie rod), as shown in Figure 1. estals are considered as the objective functions. Furthermore, the first-order results and total sensitivity index of different design parameters for each objective function are obtained, and the five most sensitive parameters are selected. Based on above five design parameters, multi-objective optimization design for vibration reduction in a simplified lubricating oil tank system is conducted based on NSGA-III, and the results of the tripleobjective optimization are obtained with the Pareto-front surface with an obvious frontier. The simulation results show that the oil tank vibration of the optimized system is suppressed effectively under the unbalanced excitation of two typical engine speeds. The established method and main results can provide guidance for designers of aeroengine external structural systems, which can help to achieve superior system dynamic performances in engineering applications.

Dynamics Model of a Rigid Body System with Multi-Point Supports
The specific structural object of an aeroengine-lubricating oil tank supported by multiple mounting pedestals can be simplified as a rigid body system with multi-point supports. The lubricating oil-tank is a typical three-support rigid body system that is mainly composed of an oil tank body, a casing, two fixed mounting pedestals, and an articulated mounting pedestal (tie rod), as shown in Figure 1. The mounting location, mounting angle, elastic modulus of the mounting pedestals, cross-sectional area, initial length of the tie rod, and the damper damping coefficient of the mounting pedestals have a certain influence on the system's vibration response and the reaction force of the mounting pedestals. A system dynamics model is established considering the above-mentioned mounting pedestal design parameters.
In order to the reduce vibration of the oil tank with a thin-walled structure, the internal supporting structure of the oil tank is mainly composed of a truss welded by a variety of square steel tubes. In the present study, the influence of the structural parameters of The mounting location, mounting angle, elastic modulus of the mounting pedestals, cross-sectional area, initial length of the tie rod, and the damper damping coefficient of the mounting pedestals have a certain influence on the system's vibration response and the reaction force of the mounting pedestals. A system dynamics model is established considering the above-mentioned mounting pedestal design parameters.
In order to the reduce vibration of the oil tank with a thin-walled structure, the internal supporting structure of the oil tank is mainly composed of a truss welded by a variety of square steel tubes. In the present study, the influence of the structural parameters of the mounting pedestals on the oil tank body vibration response is identified without considering the complex vibration of oil tank's thin-walled structure, so the concentrated mass simplification of the oil tank is processed. Considering that the oil tank adopts a three-support type of tie rod and two mounting pedestals, the main vibration behavior of the oil tank is the rotation along the two fixed mounting pedestals, and herein the oil tank is simplified as a rigid rod and the tie rod is simplified as a flexible rod. The simplified mechanical model of the oil tank supported by multiple mounting pedestals is illustrated in Figure 2. In the above model, the initial length of the tie rod, the location of the rear mounting pedestals, and the distance from the center of mass to the rear support are denoted by l 0 , d, and e, respectively. θ 0 and β 0 denote the mounting angles of the mounting pedestals at two ends. The oil tank rotation angle is θ. The oil tank is subjected to the torque of the tie rod M C , the external excitation of the casing Hsin(ωt), and its own gravity mg. The pulling force of the front mount by the tie rod is F 1 , and the reaction force of the rear mounting pedestals is F 2 ; these are decomposed into F n and F τ along the normal and tangential directions of the simplified oil tank, respectively. mounting pedestals, and the distance from the center of mass to the rear support are denoted by l0, d, and e, respectively. θ0 and β0 denote the mounting angles of the mounting pedestals at two ends. The oil tank rotation angle is θ. The oil tank is subjected to the torque of the tie rod MC, the external excitation of the casing Hsin(ωt), and its own gravity mg. The pulling force of the front mount by the tie rod is F1, and the reaction force of the rear mounting pedestals is F2; these are decomposed into Fn and Fτ along the normal and tangential directions of the simplified oil tank, respectively. A force analysis of the simplified oil tank system is carried out and the forces are, respectively, decomposed along the normal and tangential directions of the oil tank's rotational motion, obtaining Normal direction: Tangential direction: The reaction force of the rear mounting pedestals can be calculated as Considering that the vibration response of the oil tank θ is small, the deformation of the tie rod can be expressed as where E represents the elastic modulus of tie rod and A represents the cross-sectional area of the tie rod. The moment of the oil tank applied by the tie rod can be calculated as A force analysis of the simplified oil tank system is carried out and the forces are, respectively, decomposed along the normal and tangential directions of the oil tank's rotational motion, obtaining The reaction force of the rear mounting pedestals can be calculated as Considering that the vibration response of the oil tank θ is small, the deformation of the tie rod can be expressed as Introducing the mounting angles of the oil tank α = θ 0 + β 0 , the tie rod tension can be obtained as where E represents the elastic modulus of tie rod and A represents the cross-sectional area of the tie rod. The moment of the oil tank applied by the tie rod can be calculated as The oil tank torque balance equation can be described as where I represents the moment of inertia of the oil tank rotating along the mounting pedestals and γ represents the damping coefficient of the mounting pedestals. Substituting Formulas (3) and (5) into Formula (6) yields The oil tank vibration differential equation can be written as The natural frequency of oil tank rotation is obtained as Assuming that the oil tank is subjected to external excitation H sin(ωt) from the casing, etc., and the torque amplitude can be calculated as where m represents mass of the empty oil tank, a represents external acceleration excitation, and e is the distance from the center of mass of the oil tank to the rear mounting pedestals. The forced vibration differential equation of the simplified oil tank system can be expressed as Supposing that the equation is solved as Substituting Equation (12) into Equation (11) yields The vibration response of the oil tank can be obtained as

Design Parameters and Objective Functions of a Rigid Body System with Multi-Point Supports
In the design process of the oil tank and its mounting pedestals, the system vibration response amplitude |θ|, the front mount reaction force amplitude |F 1 | and the rear mount reaction force amplitude |F 2 | are defined as three objective functions. The mounting location of rear mounting pedestals d, mounting angle α, the elastic modulus E, crosssectional area A, initial length l 0 of the tie rod, the damping coefficient of the mounting pedestals γ, and the external acceleration excitation α are defined as decision variables. Sobol's sensitivity analysis method (SSAM) is used to investigate the design parameters' sensitivity to each objective function. Based on the above results, the NSGA-III is employed to optimize system objective functions.

Objective Functions
The three objective functions of the simplified oil tank system are defined as

Domain Intervals of Design Parameters
The design parameters of mounting pedestals are given in Table 1. The structural parameters, material parameters, and load parameters of mounting pedestals have a sufficiently wide range of values, which helps to obtain optimal results. Table 1. Domain intervals of design parameters.

Sensitivity Analysis of Structural Parameters of Mounting Pedestals
In this section, the mounting location of rear mounting pedestals d, the mounting angle α, the elastic modulus E, cross-sectional area A, initial length l 0 of the tie rod, the damping coefficient of the mounting pedestals γ, and the external acceleration excitation a are defined as decision variables, and SSAM is employed to investigate the design parameters' sensitivity to each objective function.

Basic Principle of Sobol's Sensitivity Analysis
The SSAM can not only quantify the influence degree of a single system parameter on the system, but can also analyze the influence degree of the simultaneous change of multiple parameters and the coupling effect between parameters on the system [12,13]. The basic process of SSAM is: (1) construct a mathematical model, which is the basis for the calculation of the sensitivity index; (2) clarify the domain interval of each design parameter; (3) sample the domain space of the design parameters and generate the corresponding sampling matrix; (4) calculate the parameter sensitivity index.
The mathematical model and parameter domain intervals are given in Sections 3.1 and 3.2, respectively. The next step is to sample in the domain space of design parameters and generate the corresponding sampling matrix. Domain intervals of seven structural parameters are given, and we need to perform quasi Monte Carlo sampling on the design parameters via the Sobol method to generate two mutually independent m×n sampling matrices, A and B, where m is the number of samples for each parameter and the number is taken as m = N, n is the number of design parameters (n = 7), and the dimensions of the sampling matrices A and B are N×7. Meanwhile, another two matrices A   m represents the m-th row of the corresponding sampling matrix, and the row vector in the matrix can be expressed as where x 1 represents the mounting location of the rear mounting pedestals d; x 2 represents mounting angle α; x 3 , x 4 , and x 5 represent the elastic modulus E, cross-sectional area A, and the initial length l 0 of the tie rod, respectively; x 6 represents the damping coefficient of the mounting pedestals γ; and x 7 represents the external acceleration excitation a.
The estimation of the first-order influence index S i and total sensitivity index ST i can be expressed as The first-order influence index S i reflects the contribution degree of a single design parameter X i to the output result. The total sensitivity index ST i reflects the influence of a single design parameter X i on the output result and the influence of the interaction effect between X i and the other parameters on the output result.
Finally, the above-mentioned calculation process is simulated via MATLAB software, and the influence weight of each structural parameter of the simplified lubricating oil tank system on each objective function can be obtained. Furthermore, the total sensitivity index of each parameter can be calculated.

Sensitivity Analysis Results
The more samples are taken, the more the calculation results tend to converge. Numbers of samples are taken as m= (7, 507, 1007, 1507, . . . , 9507), and the first-order influence index S i and total sensitivity index ST i of each parameter are calculated.
For the objective function f 1 (X), the total sensitivity index of system parameters is much larger than the first-order sensitivity index of system parameters, indicating that the coupling effect of different parameters is significant. It can be observed from the total sensitivity index that the mounting angle α has the greatest influence on the function, followed by the initial length of the tie rod l 0 and the mounting location of the oil tank d, as shown in Figure 3 and Table 2.

Sensitivity Analysis Results
The more samples are taken, the more the calculation results tend to converge. Numbers of samples are taken as m= (7, 507, 1007, 1507, …, 9507), and the first-order influence index Si and total sensitivity index STi of each parameter are calculated.
For the objective function f1(X), the total sensitivity index of system parameters is much larger than the first-order sensitivity index of system parameters, indicating that the coupling effect of different parameters is significant. It can be observed from the total sensitivity index that the mounting angle α has the greatest influence on the function, followed by the initial length of the tie rod l0 and the mounting location of the oil tank d, as shown in Figure 3 and Table 2. For the objective function f2(X), it can be observed from the total sensitivity index that the mounting angle α has the greatest influence on the function, and the coupling effect of other design parameters is not obvious. It can be observed from the first-order sensitivity index that the mounting location d, the mounting angle α of the oil tank, and the initial length of the tie rod l0 are the three most sensitive parameters, in turn, as shown in Figure  4 and Table 2.  For the objective function f 2 (X), it can be observed from the total sensitivity index that the mounting angle α has the greatest influence on the function, and the coupling effect of other design parameters is not obvious. It can be observed from the first-order sensitivity index that the mounting location d, the mounting angle α of the oil tank, and the initial length of the tie rod l 0 are the three most sensitive parameters, in turn, as shown in Figure 4 and Table 2.
shown in Figure 3 and Table 2. For the objective function f2(X), it can be observed from the total sensitivity index that the mounting angle α has the greatest influence on the function, and the coupling effect of other design parameters is not obvious. It can be observed from the first-order sensitivity index that the mounting location d, the mounting angle α of the oil tank, and the initial length of the tie rod l0 are the three most sensitive parameters, in turn, as shown in Figure  4 and Table 2. For the objective function f3(X), the coupling effect between the various system parameters is significant. It can be observed from the total sensitivity index that the mounting angle α and mounting location d of the oil tank have the greatest influence on the function, and they are almost equal, followed by cross-sectional area A and the elastic modulus E of the tie rod, as shown in Figure 5 and Table 2. For the objective function f 3 (X), the coupling effect between the various system parameters is significant. It can be observed from the total sensitivity index that the mounting angle α and mounting location d of the oil tank have the greatest influence on the function, and they are almost equal, followed by cross-sectional area A and the elastic modulus E of the tie rod, as shown in Figure 5 and Table 2. Furthermore, the specific first-order sensitivity index and total sensitivity index of each parameter are calculated, as listed in Table 2. Table 2. First-order and total sensitive index of each parameter for different objective functions. Furthermore, the specific first-order sensitivity index and total sensitivity index of each parameter are calculated, as listed in Table 2.

Multi-Objective Optimization of a Rigid Body System with Multi-Point Supports
Based on the sensitivity analysis results in Section 4.2, the mounting location of the rear mounting pedestals d, the mounting angle α, the elastic modulus E, the cross-sectional area A, and the initial length l 0 of the tie rod are defined as decision variables, and NSGA-III is employed to optimize the system's objective functions.

Definition of Pareto Optimal Solution
Generally, the mathematical model of a multi-objective optimization problem can be expressed as where x = (x 1 , . . . , x n ) ∈ X ⊂ R n is an n-dimensional decision vector, X is defined as an n-dimensional decision space, y = (y 1 , . . . , y n ) ∈ Y ⊂ R n is an m-dimensional object vector, and Y is an m-dimensional object space. Additionally, g l (x) ≤ 0, i = 1, 2, . . . , q represents q inequality constraints and h l (x) = 0, j = 1, 2, . . . , p represents p equality constraints. For multi-objective optimization problems there is usually a set of solutions, and it is impossible to improve any objective function without weakening at least one other objective function. Such solutions are called nondominated solutions, or Pareto optimal solutions, and are defined as follows:

Definition 2 (feasible solution set)
: the set of all feasible solutions in X is called the feasible solution set, denoted by X f , X f ⊆ X.

Definition 3 (Pareto dominated)
: assuming x A , x B ∈ X f are two feasible solutions, if then it is called x A dominant x B , denoted by x A ≺ x B .

Definition 4 (Pareto optimal solution): a solution, x * ∈ X f , is called a non-dominated solution if and only if
Definition 5 (Pareto optimal solution set): the Pareto optimal solution set is the set of all non-dominated solutions, which is defined as follows Definition 6 (Pareto-front surface): defining the Pareto optimal solution set as P * , the Paretofront surface PF * can be defined as

Steps of NSGA-III Algorithm
The traditional multi-objective optimization of evolutionary algorithms (MOEAs) cannot satisfy the requirements for online data processing because of its high computational costs. This drawback makes it difficult for MOEAs to solve multi-objective, large-scale optimization problems. Among different evolutionary algorithms, the NSGA-III is a fairly new method capable of solving large-scale optimization problems with a good ability to search for Pareto optimal solution sets. The main idea of NSGA-III is to introduce a reference point mechanism on the basis of NSGA-II and retain those individuals that are non-dominated and close to the reference points. NSGA-III and NSGA-II have similar frameworks. The main difference between the two algorithms is the selection mechanism. NSGA-II is mainly sorted by crowding degree, and its effect is not obvious in high-dimensional objective function spaces. NSGA-III makes a drastic adaptation of the crowding degree, maintaining population diversity by widely introducing distributed reference points. The flow chart of the proposed scheme is presented in Figure 6.

Optimization Results and Verification
Based on the proposed mathematic model, tri-objective optimization was carried out using the NSGA-II and NSGA-III algorithms. Then, vibration responses under the unbalanced excitation of two typical engine speeds (n1 = 10,200 r/min and n2 = 13,200 r/min) were calculated to validate the optimization results. Figure 7 is the Pareto front obtained by NSGA-III optimization with tri-objective functions under engine rotation speed n1 = 10,200 r/min (ω1 = 170 Hz). The results of points A1 and points A2 obtained by NSGA-II and NSGA-III optimization are listed in Table A1.

Optimization Results and Verification
Based on the proposed mathematic model, tri-objective optimization was carried out using the NSGA-II and NSGA-III algorithms. Then, vibration responses under the unbalanced excitation of two typical engine speeds (n 1 = 10,200 r/min and n 2 = 13,200 r/min) were calculated to validate the optimization results. Figure 7 is the Pareto front obtained by NSGA-III optimization with tri-objective functions under engine rotation speed n 1 = 10,200 r/min (ω 1 = 170 Hz). The results of points A1 and points A2 obtained by NSGA-II and NSGA-III optimization are listed in Table A1. It can be seen that A1 has the lowest vibration response and reaction force of the rear mount, while it has the largest reaction force of the front mount. A2 has the lowest vibration response and reaction force of the front mount, while it has the largest reaction force of the rear mount. The optimized results are similar under the unbalanced excitation of engine speed n 2 = 13,200 r/min (ω 2 = 220 Hz), as shown in Figure 8.  In order to validate the optimization results and compare the two optimization methods, a reference point Pref.A with mounting pedestal parameters given in Table A1 is selected, and its corresponding objective functions under the unbalanced excitation of engine speed ω1 = 170 Hz were calculated. It can be observed from the optimization results in Table A1 that the vibration of point A1 obtained via NSGA-II dropped 35% and the reaction force of the rear mount dropped 80%, compared with those of point Pref.A. Additionally, the vibration of point A1 obtained via NSGA-III dropped 44% and the reaction force of the rear mount dropped 94%, compared with those of point Pref.A. The vibration of point A2 dropped 41% and the reaction force of the front mount dropped 94%, com-  In order to validate the optimization results and compare the two optimization methods, a reference point Pref.A with mounting pedestal parameters given in Table A1 is selected, and its corresponding objective functions under the unbalanced excitation of engine speed ω1 = 170 Hz were calculated. It can be observed from the optimization results in Table A1 that the vibration of point A1 obtained via NSGA-II dropped 35% and the reaction force of the rear mount dropped 80%, compared with those of point Pref.A. Additionally, the vibration of point A1 obtained via NSGA-III dropped 44% and the reaction force of the rear mount dropped 94%, compared with those of point Pref.A. The vibration of point A2 dropped 41% and the reaction force of the front mount dropped 94%, com- In order to validate the optimization results and compare the two optimization methods, a reference point P ref.A with mounting pedestal parameters given in Table A1 is selected, and its corresponding objective functions under the unbalanced excitation of engine speed ω 1 = 170 Hz were calculated. It can be observed from the optimization results in When the simplified lubricating oil tank was subjected to the unbalanced excitation of engine speed ω 2 = 220 Hz, the optimization results were similar to those of the unbalanced excitation of engine speed ω 1 = 170 Hz.
The vibration response of the two optimized results (points A1 and Pref.A, points B1 and P ref.B ) under different unbalanced excitation engine speeds are shown in Figure 9. It can be seen that as the speed increases, the vibration response of the oil tank increases. Compared with the vibration response of the unoptimized points (P ref.A and P ref.B ), that of the optimized points (points A1 and points B1) drops significantly, and the optimization results obtained via NSGA-III are better than those of NSGA-II.

Conclusions
In the current study, the non-dominated sorting genetic algorithm (NSGA-III) combined with Sobol's sensitivity analysis method has been employed to optimize an aeroengine simplified lubricating oil tank supported by multiple mounting pedestals. In the optimization process, several structural parameters of mounting pedestals are taken as the design variables, and the system vibration response as well as the reaction forces of the front and rear mounting pedestals are considered as the objective functions. The results of the triple-objective optimization are obtained as Pareto-front surfaces. The following conclusions are obtained: • By Sobol's sensitivity analysis method, the importance ranking of mounting pedestals parameters can be obtained. Many other sensitivity analysis problems in engineering applications can adopt the method used in this paper. • By using NSGA-III, multi-objective optimization problems of vibration reduction for aeroengines with simplified lubricating oil tanks supported by multiple mounting pedestals can be solved efficiently. The Pareto optimal solution set can be obtained, by which the designer can consider different objective functions simultaneously and make final decisions according to design demands. The optimization results are verified, and the system vibration is suppressed effectively under the unbalanced excitation of two typical engine speeds (n1 = 10,200 r/min and n2 = 13,200 r/min). Many

Conclusions
In the current study, the non-dominated sorting genetic algorithm (NSGA-III) combined with Sobol's sensitivity analysis method has been employed to optimize an aeroengine simplified lubricating oil tank supported by multiple mounting pedestals. In the optimization process, several structural parameters of mounting pedestals are taken as the design variables, and the system vibration response as well as the reaction forces of the front and rear mounting pedestals are considered as the objective functions. The results of the triple-objective optimization are obtained as Pareto-front surfaces. The following conclusions are obtained: • By Sobol's sensitivity analysis method, the importance ranking of mounting pedestals parameters can be obtained. Many other sensitivity analysis problems in engineering applications can adopt the method used in this paper. • By using NSGA-III, multi-objective optimization problems of vibration reduction for aeroengines with simplified lubricating oil tanks supported by multiple mounting pedestals can be solved efficiently. The Pareto optimal solution set can be obtained, by which the designer can consider different objective functions simultaneously and make final decisions according to design demands. The optimization results are verified, and the system vibration is suppressed effectively under the unbalanced excitation of two typical engine speeds (n 1 = 10,200 r/min and n 2 = 13,200 r/min). Many other multi-objective optimization problems in engineering applications can adopt the method used in this paper.

•
The optimization results of mounting pedestal parameter design for vibration reduction in a multi-support rigid body system can help to achieve better system dynamic performances according to design requirements in engineering applications. They can contribute to the practical application of structural optimization in engineering applications in order to reduce vibration and reaction force level, and improve the service time of key parts, such as mounting pedestals, in aeroengine-lubricating oil tanks.