Multi-Objective Optimization of Square Corrugation Multilayer Nested Structures

: Thin-walled structures, when used for high-speed railways, can effectively mitigate the irreversible destruction when a malfunction occurs. Nested thin-walled tubes, as energy-absorbing structures, possess excellent speciﬁc energy absorption ( SEA ) and crushing force efﬁciency ( CFE ). This paper conducts multi-objective optimization by focusing on a square corrugation nested structure with a double octagon inner wall, namely SCOD, to ameliorate the crashworthiness of the nested structure. The ﬁnite element model of the SCOD is constructed and validated by test data. A set of experimental design points with good spatial distribution are obtained using the optimal Latin hypercube (LHC) method. The polynomial response surface (PRS) method was applied to establish the ﬁtting relationship between design variables and optimization objectives, and validation is accomplished. The DCNSGA-III algorithm is employed for optimization, resulting in a Pareto alternative solution set with good population diversity and convergence. In addition, to observe the optimized performance, a set of optimal solutions considering a single objective value is derived, and a comprehensive optimal solution is obtained by applying the minimum distance selection method (TMDSM). Finally, the proposed optimized system is analyzed and validated. According to the alternative reference solutions, the initial peak force ( IPCF ) reduces by 53.75% and CFE increases by 8.7%. This paper provides some reference for the optimization design in practical engineering.


Introduction
Rail vehicle safety is a continuous and crucial problem that concerns people with the development of railway technology.For the protection of rail vehicle safety, scientists currently focus on two main aspects, namely active safety protection and passive safety protection, to conduct research [1][2][3].
For the aspect of active safety protection, track life prediction and fault diagnosis can effectively and timely maintain railways, thereby preventing railway accidents [4,5].In the study of Ngamkhanong et al., the prediction model was established by Artificial Neural Networks (ANNs) to estimate the buckling temperature of complex tracks [6].A risk analysis method based on fault tree analysis (FTA) was developed by Dindar et al., and this model provides an idea for predicting the risk of railway derailment [7].In addition, measures such as brake control and warning systems can also improve the active safety protection performance of trains [8][9][10].However, it is necessary to investigate passive safety protection, which can mitigate collision force and damage suffered during secondary accidents.
For the aspect of passive safety protection, installing guardrails [11], anti-climbing devices [12], couplers [13], and energy-absorbing devices [14,15] can effectively enhance the crashworthiness of the train.Thin-walled energy absorption tubes as energy-absorbing devices are extensively applied to the transportation industry due to extraordinary crashworthiness characteristics [16][17][18].The thin-walled tubes which are equipped in the front impact zone will first suffer crushing; the reasonable structural design of thin-walled structures can make the crushing process more stable and controllable, thereby protecting the safety of passengers in the survival space.The corrugation configuration was proposed, and researchers confirmed that corrugation tubes have good potential for crashworthiness [19].Experimental studies have shown that the maximum plastic moment at the corrugation groove can first guide the occurrence of collapse to diminish the peak force [20,21], and corrugation tubes with appropriate amplitude possess ordered progressive mode which contributes to more uniform load fluctuations [22].Adding rib plates into thin-walled tubes can generate a stable folding pattern and increase the number of corner elements [23,24].Wang et al. found that the interaction between the foam blocks and the reinforced rib can ameliorate the crash resistance performance [25].In addition, rib plates can affect significantly the energy absorption capacity and deformation mode [23,26].Another method is to set up multi-layer structures or subunits within thin-walled tubes [27,28].The deformation mode usually develops toward a better mode due to the mutual constraints between the inner and outer layers [29,30].Zhu et al. found that the number of nested layers should not be excessive; otherwise, the comprehensive crushing performance will decrease [31].Tran et al. indicated that the deformation mode is more controllable due to the resistance effect of the outer tube [32].The research of Yang et al. showed that the increase in ripple and number of layers will form a more stable and plastic hinge [33].In addition, optimizing structural parameters based on structural design to obtain the most reasonable parameter scheme is crucial for improving vehicle crashworthiness.
For the aspect of structural parameter optimization, non-dominated sorting genetic algorithm II (NSGA-II) became one of the frequently used optimization methods in various fields due to its convenience and adaptability [21,33,34].However, as the number of optimization objectives increases, low-dimensional multi-objective optimization algorithms similar to NSGA-II will lose effectiveness, and the NSGA-III algorithm used to solve highdimensional problems can be applied [3].In addition, other optimization methods, such as multi-objective particle swarm optimization (MOPSO) and Artificial Neural Networks (ANN), are commonly used for crashworthiness structural optimization.Albak adopted MOPSO to conduct multi-objective optimization for optimal structures selected by the decision-making method [35].Qiu et al. utilized the MOPSO algorithm to optimize the front rail and multi-cell tube to enhance the crashworthiness of the vehicle [36].Pirmohammad et al. conducted bi-objective optimization based on the shape of the cross-section by using the ANN algorithm and selected the optimal structure according to the rank results of the technique for ordering preferences by similarity to the ideal solution (TOPSIS) method [37].However, the general evolutionary algorithm ignores the diversity of the population while pursuing the constraint of satisfaction, resulting in the population being too concentrated in a feasible or infeasible region.Jiao et al. developed a DCNSGA-III by combining dynamic-constrained many-objective optimization problems (CMaOPs) and NSGA-III [38].This method can enable the population to span large and infeasible regions, and the proposed environment selection operator can select higher-quality solutions.In addition, few researchers have applied this approach to multi-optimization for energyabsorbing structures.Therefore, this paper intends to combine DCNSGA-III to obtain more convergent and diverse non-dominated solutions.
In summary, a nested structure is proposed inspired by the three elements: corrugation, rib plate, and multi-layer structure.In addition, the octagonal wall can absorb more membrane energy with the multi-corner advantage according to previous researchers.The study of Li et al. indicates that square tubes have better deformation resistance potential compared to circular tubes [39].In this paper, a square corrugation nested tube with a double octagon inner structure (SCOD) is proposed.Distinguishing from previous studies, this research mainly includes three aspects: (1) proposing a novel nested energy-absorbing tube, namely SCOD; (2) developing a multi-objective optimization procedure combining DCNSGA-III algorithm with TMDSM to obtain alternative solutions with better performance; and (3) verifying the effectiveness of the optimization approach proposed and making an evaluation of it.
In this paper, the pattern of organization is described as follows.Section 2 introduces the primordial parameters and finite element modeling of SCOD.Section 3 reveals the methodology and optimization framework.In Section 4, the results of optimization are obtained.In Section 5, discussions and analyses were conducted on the optimization plan.Finally, Section 6 summarizes the entire work.

The Geometric Parameter of the SCOD Tube
The square corrugation nested tube with a double octagonal inner wall (SCOD) is shown in Figure 1.The SCOD can be divided into four sections, including the outer corrugation wall, the mid-octagon wall, the inner octagon wall, and rib plates.The four different colored lines correspond to the four thickness parameters of the sections.
membrane energy with the multi-corner advantage according to previous researchers.The study of Li et al. indicates that square tubes have better deformation resistance potential compared to circular tubes [39].In this paper, a square corrugation nested tube with a double octagon inner structure (SCOD) is proposed.Distinguishing from previous studies, this research mainly includes three aspects: (1) proposing a novel nested energy-absorbing tube, namely SCOD; (2) developing a multi-objective optimization procedure combining DCNSGA-III algorithm with TMDSM to obtain alternative solutions with better performance; and (3) verifying the effectiveness of the optimization approach proposed and making an evaluation of it.
In this paper, the pattern of organization is described as follows.Section 2 introduces the primordial parameters and finite element modeling of SCOD.Section 3 reveals the methodology and optimization framework.In Section 4, the results of optimization are obtained.In Section 5, discussions and analyses were conducted on the optimization plan.Finally, Section 6 summarizes the entire work.

The Geometric Parameter of the SCOD Tube
The square corrugation nested tube with a double octagonal inner wall (SCOD) is shown in Figure 1.The SCOD can be divided into four sections, including the outer corrugation wall, the mid-octagon wall, the inner octagon wall, and rib plates.The four different colored lines correspond to the four thickness parameters of the sections.The design of the outer corrugation wall is based on a square shape with a side length of 40 mm.The corrugation on the four side lines is a cosine curve which can be represented in the Cartesian coordinate system as follows: where L0 denotes the length of the side; the parameter A denotes the amplitude of the cosine curve.The initial thickness t of all sections is designed to be 1 mm.To ensure adequate compression distance, the height of the nested tubes is set to 100 mm [40].Two octagonal sub-thin-walled structures, namely the mid-wall and inner wall, are embedded within the corrugated outer wall.The mid-wall is circumscribed from a circular shape with a diameter of 22√2 mm; The inner wall is circumscribed from a circular shape The design of the outer corrugation wall is based on a square shape with a side length of 40 mm.The corrugation on the four side lines is a cosine curve which can be represented in the Cartesian coordinate system as follows: where L 0 denotes the length of the side; the parameter A denotes the amplitude of the cosine curve.The initial thickness t of all sections is designed to be 1 mm.To ensure adequate compression distance, the height of the nested tubes is set to 100 mm [40].Two octagonal sub-thin-walled structures, namely the mid-wall and inner wall, are embedded within the corrugated outer wall.The mid-wall is circumscribed from a circular shape with a diameter of 22 √ 2 mm; The inner wall is circumscribed from a circular shape with a diameter of 10 √ 2 mm.For the design of the size of octagonal sub-thin-walled structures, we refer readers to the study of Albak [35].

The Establish of the FE Model
The finite element (FE) method has been proven to have sufficient precision and is often utilized for the crashworthiness analysis of tubes [41].In this paper, the characteristic is analyzed by using the nonlinear explicit FE analysis program LS_DYNA.The simulation model of SCOD is shown in Figure 2.
with a diameter of 10√2 mm.For the design of the size of octagonal sub-thin-walled struc tures, we refer readers to the study of Albak [35].

The Establish of the FE Model
The finite element (FE) method has been proven to have sufficient precision and i often utilized for the crashworthiness analysis of tubes [41].In this paper, the characteris tic is analyzed by using the nonlinear explicit FE analysis program LS_DYNA.The simu lation model of SCOD is shown in Figure 2. From Figure 2, the simulation construction of SCOD can be divided into three com ponents: motional plate, SCOD tube, and stationary plate.To save computational costs the moving plate imposes a load at a constant speed of 3 m/s [39].The fixed plate is re stricted to six degrees of freedom, and the SCOD tube is placed in the fixed plate.When the moving plate moves at a constant speed, the fixed plate limits the longitudinal dis placement of the SCOD in the longitudinal direction to ensure successful crushing.To guarantee that the tube is not crushed and that the endergonic particularity of the mode can be observed, the tube is compressed by 60% of its original length, and loading i stopped.
The Belytschko-Tsay shell element with one center integration point and fiv through-thickness integration points is employed to establish the FE model of SCOD.Th convergence analysis of mesh is obtained in Figure 3.As shown in the analysis results, th variation of MCF has a slight change with the increase in mesh size, but the solving tim is greatly reduced.To ensure sufficient accuracy and computational efficiency, we decided to adopt the mesh with the size of 1 mm × 1 mm to construct numerical models.Becaus the plate hardly deforms during the collision which can be regarded as a rigid body, both of the plates are simulated with MAT_RIGID (Mat_20), and all nodes will maintain a con stant relative displacement in this condition.The SCOD tubes are simulated with MAT_PIECEWISE_LINEAR_PLASTICITY (Mat_24) which is a model that reflects th elastic-plastic mechanical properties of materials.This type of material has an arbitrary stress-strain curve, and any strain rate dependence can be defined.In addition, the strain rate can be explained by the Cowper-Symonds model.The "AUTO MATIC_NODE_TO_SURFACE" is applied to define the contact between SCOD and th plate, and the dynamic and static friction coefficients are both 0.3.And the "AUTO MATIC_SINGLE_SURFACE" contact algorithm is employed to define the self-contact o the SCOD to avoid the occurrence of a pressure breakdown phenomenon.The dynami and static friction coefficients of self-contact are both 0.3 as well.The material of the SCOD tube is aluminum alloy AA6061O, and the constitutive model and material characteriza tion are shown in Table 1.From Figure 2, the simulation construction of SCOD can be divided into three components: motional plate, SCOD tube, and stationary plate.To save computational costs, the moving plate imposes a load at a constant speed of 3 m/s [39].The fixed plate is restricted to six degrees of freedom, and the SCOD tube is placed in the fixed plate.When the moving plate moves at a constant speed, the fixed plate limits the longitudinal displacement of the SCOD in the longitudinal direction to ensure successful crushing.To guarantee that the tube is not crushed and that the endergonic particularity of the model can be observed, the tube is compressed by 60% of its original length, and loading is stopped.
The Belytschko-Tsay shell element with one center integration point and five throughthickness integration points is employed to establish the FE model of SCOD.The convergence analysis of mesh is obtained in Figure 3.As shown in the analysis results, the variation of MCF has a slight change with the increase in mesh size, but the solving time is greatly reduced.To ensure sufficient accuracy and computational efficiency, we decided to adopt the mesh with the size of 1 mm × 1 mm to construct numerical models.Because the plate hardly deforms during the collision which can be regarded as a rigid body, both of the plates are simulated with MAT_RIGID (Mat_20), and all nodes will maintain a constant relative displacement in this condition.The SCOD tubes are simulated with MAT_PIECEWISE_LINEAR_PLASTICITY (Mat_24) which is a model that reflects the elastic-plastic mechanical properties of materials.This type of material has an arbitrary stress-strain curve, and any strain rate dependence can be defined.In addition, the strain rate can be explained by the Cowper-Symonds model.The "AU-TOMATIC_NODE_TO_SURFACE" is applied to define the contact between SCOD and the plate, and the dynamic and static friction coefficients are both 0.3.And the "AUTO-MATIC_SINGLE_SURFACE" contact algorithm is employed to define the self-contact of the SCOD to avoid the occurrence of a pressure breakdown phenomenon.The dynamic and static friction coefficients of self-contact are both 0.3 as well.The material of the SCOD tube is aluminum alloy AA6061O, and the constitutive model and material characterization are shown in Table 1.From Figure 1, it can be easily seen that the SCOD is formed based on the primary cosine outline; therefore, in this study, the square corrugation square tube, namely SC will be utilized to check the exactitude of the FE simulation.The SC used for verification has a normal side length with 40 mm and an amplitude with 3 mm, which is the same as the geometric parameters of SCOD.The SC was compressed under quasi-static compres sion by using the test machine MTS Landmark, and the loading speed was set as 2 mm/min.The quasi-static compression test is shown in Figure 4.The tubes used to test and FE simulation have the same height and thickness, which are 40 mm and 100 mm, respectively.Furthermore, the test and numerical tube have iden tical material properties.To ensure the consistency of experimental conditions, SC was compressed to 60% of the initial distance in both the testing machine and the finite elemen numerical simulation.In addition, the collapse modes and force-displacement of SC un der an experiment test and FE simulation are compared in Figure 5.  From Figure 1, it can be easily seen that the SCOD is formed based on the primary cosine outline; therefore, in this study, the square corrugation square tube, namely SC, will be utilized to check the exactitude of the FE simulation.The SC used for verification has a normal side length with 40 mm and an amplitude with 3 mm, which is the same as the geometric parameters of SCOD.The SC was compressed under quasi-static compression by using the test machine MTS Landmark, and the loading speed was set as 2 mm/min.The quasi-static compression test is shown in Figure 4.  From Figure 1, it can be easily seen that the SCOD is formed based on the primary cosine outline; therefore, in this study, the square corrugation square tube, namely SC, will be utilized to check the exactitude of the FE simulation.The SC used for verification has a normal side length with 40 mm and an amplitude with 3 mm, which is the same as the geometric parameters of SCOD.The SC was compressed under quasi-static compression by using the test machine MTS Landmark, and the loading speed was set as 2 mm/min.The quasi-static compression test is shown in Figure 4.The tubes used to test and FE simulation have the same height and thickness, which are 40 mm and 100 mm, respectively.Furthermore, the test and numerical tube have identical material properties.To ensure the consistency of experimental conditions, SC was compressed to 60% of the initial distance in both the testing machine and the finite element numerical simulation.In addition, the collapse modes and force-displacement of SC under an experiment test and FE simulation are compared in Figure 5.The tubes used to test and FE simulation have the same height and thickness, which are 40 mm and 100 mm, respectively.Furthermore, the test and numerical tube have identical material properties.To ensure the consistency of experimental conditions, SC was compressed to 60% of the initial distance in both the testing machine and the finite element numerical simulation.In addition, the collapse modes and force-displacement of SC under an experiment test and FE simulation are compared in Figure 5. From Figure 5, the tendency of force change over time of SC under test and numerical simulation coincide in principle, and the deformation modes are accordant.In addition, the crashworthiness indexes obtained by experiment and FE modeling are displayed in Table 2.

The Mathematical Model Based on SCOD
Based on the Simplified Super Folding Element (SSFE) theory, a mathematical model is proposed to predict the mean crushing force (MCF) of SCOD [42][43][44].During the collision, the energy balance equation of SCOD can be obtained as shown in Equation ( 2): where Fm is MCF; Eb and Em are the bending energy and membrane energy, which are primary forms of energy dissipation; 2H represents the wavelengths of the fold; and k is the effective compression factor, which is 0.7 in this paper [45].
A. The prediction of Eb.
The bending energy can be calculated as follows [46]: where Lc is the perimeter of the cross-section, and M0 is the full plastic bending moment per unit length, which can be obtained as Equation ( 4): where t is the thickness of the wall and σ0 is the flow stress of the structural material: B. The prediction of Em.From Figure 5, the tendency of force change over time of SC under test and numerical simulation coincide in principle, and the deformation modes are accordant.In addition, the crashworthiness indexes obtained by experiment and FE modeling are displayed in Table 2. From Table 2, the IPCF of the test and FE model are 18.34 kN and 17.95 kN, respectively, which are in reasonable agreement.It is obvious that the error between the test and FE modeling is less than 2.5%, which also indicates that the FE numerical simulation has adequate accuracy.

The Mathematical Model Based on SCOD
Based on the Simplified Super Folding Element (SSFE) theory, a mathematical model is proposed to predict the mean crushing force (MCF) of SCOD [42][43][44].During the collision, the energy balance equation of SCOD can be obtained as shown in Equation ( 2): where F m is MCF; E b and E m are the bending energy and membrane energy, which are primary forms of energy dissipation; 2H represents the wavelengths of the fold; and k is the effective compression factor, which is 0.7 in this paper [45].
A. The prediction of E b .
The bending energy can be calculated as follows [46]: where L c is the perimeter of the cross-section, and M 0 is the full plastic bending moment per unit length, which can be obtained as Equation ( 4): where t is the thickness of the wall and σ 0 is the flow stress of the structural material: B. The prediction of E m .
According to SSFE, the SCOD can be divided into several constituent elements in order to analyze theoretically the E m, , and then the constituent elements will be simplified [40,46].The constituent elements of SCOD are defined as five types, including the shell corner element, T-shape element, Y-shape element, 3-plane corner element, and criss-cross element.The specific situation is shown in Figure 6.
Appl.Sci.2023, 13, x FOR PEER REVIEW 7 of 23 According to SSFE, the SCOD can be divided into several constituent elements in order to analyze theoretically the Em,, and then the constituent elements will be simplified [40,46].The constituent elements of SCOD are defined as five types, including the shell corner element, T-shape element, Y-shape element, 3-plane corner element, and criss-cross element.The specific situation is shown in Figure 6.As shown as Figure 6, the shell corner element can be simplified as a plane corner element (180°) [47]: where θ = 180°.The T-shape element is equivalent to the three-plane corner element (30°) [35], which can be obtained by Equation (7): The Y-shape element can be regarded as the three-plane corner element with the α = 90° [47,48], which can be obtained by Equation ( 8): The Em of the criss-cross element is [42]: As mentioned above, the total membrane energy can be summarized as: C. The prediction of MCF.By uniting Equations ( 2), ( 3) and (10), the predicted MCF of the SCOD can be expressed as: In addition, parameter H can be calculated by: As shown as Figure 6, the shell corner element can be simplified as a plane corner element (180 • ) [47]: where θ = 180 • .The T-shape element is equivalent to the three-plane corner element (30 • ) [35], which can be obtained by Equation ( 7): The Y-shape element can be regarded as the three-plane corner element with the α = 90 • [47,48], which can be obtained by Equation ( 8): The E m of the criss-cross element is [42]: As mentioned above, the total membrane energy can be summarized as: C. The prediction of MCF.
By uniting Equations ( 2), ( 3) and ( 10), the predicted MCF of the SCOD can be expressed as: In addition, parameter H can be calculated by: By calculating Formulas ( 11) and ( 12), F m SCOD can be solved, and the results predicted by the mathematical model are compared with the FE model, as shown in Table 3.In terms of Table 3, the relative error between the FE model and the mathematical model is 2.24%, which indicates that FE modeling satisfies the accuracy requirements and can be employed in further study.

The Sensitivity Analysis of SCOD
To investigate the effect of structural parameters of SCOD on energy-absorbing performance, the sensitivity analyses are conducted in this section, which is shown in Figures 7 and 8.

R PEER REVIEW 8 of 23
By calculating Formulas ( 11) and ( 12), Fm SCOD can be solved, and the results predicted by the mathematical model are compared with the FE model, as shown in Table 3.In terms of Table 3, the relative error between the FE model and the mathematical model is 2.24%, which indicates that FE modeling satisfies the accuracy requirements and can be employed in further study.

The Sensitivity Analysis of SCOD
To investigate the effect of structural parameters of SCOD on energy-absorbing performance, the sensitivity analyses are conducted in this section, which is shown in Figures 7 and 8.According to Figure 7, with the increase in thickness tint, tmid, tout, and trib, the three crashworthiness indicators, namely IPCF, SEA, and CFE, all show an upward trend.It can be concluded that the wall thickness of the SCOD and IPCF basically exhibit a linear func- By calculating Formulas ( 11) and ( 12), Fm SCOD can be solved, and the results predicted by the mathematical model are compared with the FE model, as shown in Table 3.In terms of Table 3, the relative error between the FE model and the mathematica model is 2.24%, which indicates that FE modeling satisfies the accuracy requirements and can be employed in further study.

The Sensitivity Analysis of SCOD
To investigate the effect of structural parameters of SCOD on energy-absorbing per formance, the sensitivity analyses are conducted in this section, which is shown in Figures 7 and 8.According to Figure 7, with the increase in thickness tint, tmid, tout, and trib, the three crashworthiness indicators, namely IPCF, SEA, and CFE, all show an upward trend.It can be concluded that the wall thickness of the SCOD and IPCF basically exhibit a linear func tion relationship.Among the four parameters related to thickness, tout has the greates According to Figure 7, with the increase in thickness t int , t mid , t out , and t rib , the three crashworthiness indicators, namely IPCF, SEA, and CFE, all show an upward trend.It can be concluded that the wall thickness of the SCOD and IPCF basically exhibit a linear function relationship.Among the four parameters related to thickness, tout has the greatest impact on peak force, while t int has the smallest impact.For the sensitivity on SEA, as the thickness of the outer wall increases, the values will sacrifice slowly.In addition, when the thickness varies within the range of 1.5~2.0mm, CFE exhibits the most drastic change.
As shown in Figure 8, IPCF always rises monotonically.When D mid changes within the range of 25~35 mm, both the SEA and CFE values first increase and then decrease, and the peak appears around D mid = 30 mm.Establishing a surrogate model is regarded as an efficient method to handle complex multi-objective optimization problems (MOOPs), which can save cost and calculation time, and it has been applied extensively in the multi-objective optimization of endergonic tubes [49,50].The DoE method is a primary and indispensable step to establish an accurate surrogate model, which is also an approach to reflect integral design space by training minimum points [51].According to existing studies, the DoE generally consists of a Latin hypercube design (LHD) [48], optimal LHD [52], full factorial design, orthogonal arrays, central composite design, and Modified Extensible Lattice Sequence (MELS) [28].During optimization, the optimal LHD approach is adopted to conduct DoE.This method not only reserves excellent projective properties but also ameliorates the space-filling properties by combining novel indexes combining the Kullback-Leibler information [53,54].It has been concluded that the optimal LHC method maintains robustness as the criteria vary and is adequate to the requirement of changes in design parameters [55].

The Establish of Surrogate Model
The MOOP of crashworthiness requires extensive FE simulation or experiments.To conserve costs, proxy models are developed and used to express the response correlation between design variables and optimization goals regarded [56,57].The commonly used approaches to establish a surrogate model include the polynomial response surface (PRS) [58], kriging [51] and radial basis function (RBF) [21].As a useful method, the PRS has been proven to be credible and efficient when dealing with small group issues with noise data and a non-massive DoE sample scale [55].Therefore, the PRS establishes a surrogate model in this paper.The core idea of the PRS is fabricating nonlinear polynomials by fitting simple points to express approximately objective responses in the entire design space.The polynomial usually is expressed as follows: where α j (j = 1, 2, . .., M) in vector α represents the j-th unknown coefficient of the j-th basis function f j (x).M represents the quantity of basis functions; c represents the coefficients of each basis function by using the multiple regression method of least squares.
In addition to satisfying the requirements of the polynomial degree, four criteria for evaluating the precision of the fitting results, maximum value (MAX), mean relative error (MRE), root mean square error (RMSE) and the square value (R 2 ) are proposed [52].The formula for the four indicators is depicted as follows: 16)

The DCNSGA-III Algorithm
According to Jiao et al. [38], the DCNSGA-III method is proposed, which is created by employing reference-point-based NSGA-III to resolve constrained multi-objective optimization problems (CMaOPs).The NSGA-III method is advantageous to handle constraints and optimize objectives concurrently, which can make a CMaOP dynamic [59].The DCNSGA-III method performs great competitiveness compared with other sophisticated CMaOEAs and is competent for solving CMaOPs.The general process of DCNSGA-III is described as follows.
The DCNSGA-III, which integrates widely used NSGA-III and a problem conversion means, is applied to resolve CMaOP.Algorithm 1 displays the trunk of how the DCNSGA-III solves the CMaOP.Search incipient population P 0 ; 2 Form the incipient relaxed constrained margin δ = δ (0) , and make the iteration state t = 0; 3 while the termination requirement is not satisfied do 4 Shrink the boundary: δ = δ (t + 1) ; 5 Renew population δ-feasibility, t = t + 1; 6 Obtain the offspring population O t in terms of Algorithm 2; 7 Renew the parent population P t+1 in terms of Algorithm 3; 8 end 9 return non-dominated solutions.
At the beginning of the framework, the population P 0 is generated randomly by DCNSGA-III from unitary search space X.The non-constraint degree of every solution is computed subsequently, and the maximal non-constraint degree is regarded as the incipient constraint bounds δ.Thus, any solution in population P 0 is δ-feasible.In subsequent iterations, the constraint bounds will decrease over the evolution of the population.In addition, a modified mating and environment selection mechanism is utilized to obtain highly qualified δ-feasible solutions.

B. Select mating parents.
It is inevitable that δ-infeasible solutions are generated.In this connection, a binary tournament selection is employed to generate the offspring operator, which picks two individuals from P t .To be specific, the tournament selection intends to select an δ-feasible solution instead of an δ-infeasible solution or an δ-infeasible with a minor non-constraint level instead of an δ-infeasible with a major non-constraint level so that more δ-feasible offspring solutions are obtained.Algorithm 2 describes the main procedure of tournament selection.To screen for more feasible solutions, the following two scenarios were considered:

C. Update the parent population
For t generation, the coalition of parent and offspring population is classified into two sets, the δ-feasible set S 1 = {x ∈ U t | g(x) ≤ δ(t)} and δ-infeasible set S 2 = {x ∈ U t | g j (x) > δ j (t), j ∈ {1, . .., q}}.This process is discussed in Algorithm 3.

Algorithm 3: Renew the parent population
Input: parent population Pt, offspring population Ot To screen for more feasible solutions, the following two scenarios were considered: To screen for more feasible solutions, the following two scenarios were considered: Case 1.When |S 1 |≥ N, the non-dominated ranking and prime selection based on the reference point will proceed to choose N solutions from S 1 , and then, the N solutions will be appended to P t + 1 .
Case 2. When |S 1 |< N, there are still blank elements in P t + 1 , after the solutions S 1 are introduced.The N-S 1 rest of the solutions are selected from δ-infeasible set S 2 .Therefore, the solutions in S 2 are sorted in terms of the non-constraint level, and we append the top N − |S 1 |solutions inside the ranked S 2 into P t + 1 .D. Non-dominated ranking based on reference point method.
When the case |S 1 |≥ N occurs, non-dominated ranking based on a reference point will proceed as the environmental selection.Figure 9 displays the environmental selection of DCNSGA-III, which contains two parts: (1) the non-dominated ranking aims to obtain nondominated solutions inside the δ-feasible set S 1 to approach its PF; (2) the prime selection based on a reference point aims to retain the population multiplicity through providing a group of reference points with good spatial distribution characteristics.As shown in Figure 9, the first section is that the non-dominated ranking given m + 1 objectives: (f1(x), …, fm(x), cv(x)) is employed to classify into different non-dominated levels (F1, F2, …) and select prime solutions which have higher objective values and lower nonconstraint levels.This approach helps to optimize goals and handle constraints simultaneously.
In the second section, an interception approach is employed to intercept the optimal N solutions.The optimal k levels are determined, and k denotes the maximum which satisfies │F1 ∪ F2 ∪…∪ Fk│ < N. The k + 1 solutions │F1 ∪ F2 ∪…∪ Fk ∪Fk + 1│ − N are determined in terms of the admixture of solutions and the pre-defined reference points, which is depicted in Algorithm 4. In addition, the non-constraint objective will not be referred for the normalization of targets and the associated employment.The phenomenon is because infeasible solutions which have lower non-constraint levels are selected during the nondominated ranking.The purpose of the selection based on a reference point is to pick solutions that are close to the basis edge so that all reference points maintain the great distributed feature of solutions.Otherwise, the solution will erratically distribute within the target space.F. Evaluate computational complicacy It can be observed that the DCNSGA-III method does not requires a very large amount of calculation.The additional complicacy of computational time is reflected in the non-dominated ranking in O((m + 1) N 2 ), and m denotes the targets quantity.Since the complicacy of the original NSGA-III is O(mN 2 ), the extra complicacy of the computational time this approach brings is tolerable.As shown in Figure 9, the first section is that the non-dominated ranking given m + 1 objectives: (f 1 (x), . .., f m (x), cv(x)) is employed to classify into different non-dominated levels (F 1 , F 2 , . ..) and select prime solutions which have higher objective values and lower nonconstraint levels.This approach helps to optimize goals and handle constraints simultaneously.
In the second section, an interception approach is employed to intercept the optimal N solutions.The optimal k levels are determined, and k denotes the maximum which satisfies are determined in terms of the admixture of solutions and the pre-defined reference points, which is depicted in Algorithm 4. In addition, the non-constraint objective will not be referred for the normalization of targets and the associated employment.The phenomenon is because infeasible solutions which have lower non-constraint levels are selected during the non-dominated ranking.The purpose of the selection based on a reference point is to pick solutions that are close to the basis edge so that all reference points maintain the great distributed feature of solutions.Otherwise, the solution will erratically distribute within the target space.

E. Evaluate computational complicacy.
It can be observed that the DCNSGA-III method does not requires a very large amount of calculation.The additional complicacy of computational time is reflected in the nondominated ranking in O((m + 1) N 2 ), and m denotes the targets quantity.Since the complicacy of the original NSGA-III is O(mN 2 ), the extra complicacy of the computational time this approach brings is tolerable.
complicacy of the original NSGA-III is O(mN 2 ), the extra complicacy of the computational time this approach brings is tolerable.

The Minimum Distance Selection Method (TMDSM)
Through the above optimization method, the Pareto optimal front is obtained so that the designer can consult and make a determination according to the amount of candidate solutions.For all that, sifting the most satisfactory individual solution is essential to evaluate the final optimization result [60].The most commonly used methods for determining the optimal solution include the weight decision-making method and the minimum distance selection method (TMDSM).The former assigns weights to the optimization objectives based on their relative importance and determines the optimal solution.However, considering different actual operating conditions and preferences, it is difficult to determine appropriate weights [61].In this paper, TMDSM, which searches the "keen point" by calculating the distance between the ideal solution and each candidate solution, is applied to conduct the selection of the optimal solution.The formulation below exhibits the principle of TMDSM:

The Minimum Distance Selection Method (TMDSM)
Through the above optimization method, the Pareto optimal front is obtained so that the designer can consult and make a determination according to the amount of candidate solutions.For all that, sifting the most satisfactory individual solution is essential to evaluate the final optimization result [60].The most commonly used methods for determining the optimal solution include the weight decision-making method and the minimum distance selection method (TMDSM).The former assigns weights to the optimization objectives based on their relative importance and determines the optimal solution.However, considering different actual operating conditions and preferences, it is difficult to determine appropriate weights [61].In this paper, TMDSM, which searches the "keen point" by calculating the distance between the ideal solution and each candidate solution, is applied to conduct the selection of the optimal solution.The formulation below exhibits the principle of TMDSM: where D is the distance between the ideal solution and the "keen point".And K denotes the quantity of the objective elements; f cτ represents τ-th objective value corresponding to the c-th Pareto solution.

The Framework of the Optimization Process
In the previous study, the TOPSIS method combining multiple operating conditions is employed to determine the ideal structures from 20 nested structures proposed in order to consider comprehensively the preference of experts for various performances.To obtain the most satisfactory alternative, the DCNSGA-III method is referred to handle the multi-objective optimization issue.Therefore, the entire process combines the theories of multi-attribute decision-making (MADM) and multi-objective optimization, which can be accomplished by the following frame shown in Figure 10.

The Framework of the Optimization Process
In the previous study, the TOPSIS method combining multiple operating conditions is employed to determine the ideal structures from 20 nested structures proposed in order to consider comprehensively the preference of experts for various performances.To obtain the most satisfactory alternative, the DCNSGA-III method is referred to handle the multiobjective optimization issue.Therefore, the entire process combines the theories of multiattribute decision-making (MADM) and multi-objective optimization, which can be accomplished by the following frame shown in Figure 10.As shown in Figure 10, the optimization is roughly divided into three phases: Phase 1: The final scores of the crashworthiness criteria of each structure proposed are recorded through numerical simulation.And then, the TOPSIS based on 10 groups of weight is applied to determine the ideal structure.
Phase 2: The optimization problem is defined, and the optimal LHC is selected to conduct DoE.The surrogate model is constructed by using PRS and the DCNSGA-III method is used to acquire the candidate Pareto solution.
Phase 3: The TMDSM is utilized to solve the ideal solution in the set of alternative solutions.Finally, the results of optimization are verified to assess the exactitude of the surrogate model and the effect of the optimization model.

The Determination of the Optimization Problem
To successfully establish an optimization model, it is the priority to determine the optimization problem [37,62].In the study of Albak, the thickness of the inner wall, middle wall, outer wall and ribs of nested tubes are selected as the design variables [35].However, to explore the enhancement of the size of the inside structure, the circumscribed diameter of the middle wall of SCOD, namely Dmid, is quantified, as shown in Figure 1b.Therefore, As shown in Figure 10, the optimization is roughly divided into three phases: Phase 1: The final scores of the crashworthiness criteria of each structure proposed are recorded through numerical simulation.And then, the TOPSIS based on 10 groups of weight is applied to determine the ideal structure.
Phase 2: The optimization problem is defined, and the optimal LHC is selected to conduct DoE.The surrogate model is constructed by using PRS and the DCNSGA-III method is used to acquire the candidate Pareto solution.
Phase 3: The TMDSM is utilized to solve the ideal solution in the set of alternative solutions.Finally, the results of optimization are verified to assess the exactitude of the surrogate model and the effect of the optimization model.

The Determination of the Optimization Problem
To successfully establish an optimization model, it is the priority to determine the optimization problem [37,62].In the study of Albak, the thickness of the inner wall, middle wall, outer wall and ribs of nested tubes are selected as the design variables [35].However, to explore the enhancement of the size of the inside structure, the circumscribed diameter of the middle wall of SCOD, namely D mid , is quantified, as shown in Figure 1b.Therefore, the geometric parameters which include the thickness of the inner wall tint, the thickness of the middle wall t mid , the thickness of the outer wall tout and the cross-section size of the middle wall D mid are determined as design variables.According to previous studies, the t int , t mid , t out , t rib and D mid generated a high sensitivity impact on the crashworthiness with the range of 0.5 to 2.0 mm and 25 to 35 mm, respectively.Hence, the same interval is selected as the range of variation in this study.
In the field of practical engineering, higher SEA and CFE values and lower IPCF values are the orientation commonly pursued [29,63].However, the increase in SEA values typically accompanies an increase in IPCF, and a higher IPCF is not acceptable.In addition, the CFE represents the performance of energy absorbing, which is also a significant criterion.Hence, the SEA, CFE and IPCF are taken as optimization objectives.The mathematical model of the optimization problem can be represented as the following equation:

The Design of the Experiment (DoE)
In this section, an optimal LHC is used to contact for DoE, which has been proven to be an effective method for estimating the overall average [64,65].A total of 60 simple points are obtained by using the optimal LHC method, and all of the points are established by FE modeling.The 60 groups of SEA, CFE and IPCF values of DoE are displayed in Table A1.

The Construction of the Surrogate Model
It has been proven that the polynomial response surface (PRS) method has such bodacious effectiveness and rationality that PRS has been widely applied in multi-objective optimization design studies of thin-walled structures [36].In this study, the nonlinear quintic polynomial is used to perform the fitting relationship between the obtained parameter variables of the sample points, i.e., t int , t mid , t out , t rib and D mid , and the optimization objectives, i.e., IPCF, SEA and CFE.The established RS values, which reflect the corresponding relationship between each two parameter variables and a target value, are shown in Figures A1-A3.In addition, error analysis is applied to verify the assessment of the surrogate model in this study.To assure reasonable accuracy, four commonly used key indicators, including MAX, MRE, RMSE and R 2 are computed, which are shown in Table 4.As shown in Table 4, all the final scores of the four criteria satisfy the tolerance value of error.It should be noted that the lower MAX, MRE and RMSE are and the higher R 2 is, the better the fitting effect and the higher the accuracy.In general, the value of R 2 should not be more than 1 and should not be less than 0.88.It can be concluded that the proxy model has satisfactory exactitude.In Figure 11, the Pareto solutions obtained by the DCNSGA-III algorithm can be observed in the 3D space and the projections are also marked in planes that reflect the relationships between every two target values.As the value of IPCF increases, the SEA also increases.The changes between IPCF and CFE reflect the same trend.It is worth noting that as the SEA value increases, the CFE also increases, indicating that while reducing the cost indicator, the benefit indicator will also be sacrificed.In addition to considering the optimal solution obtained through TMDSM, the other three optimal single performance solutions generated by referencing a single target value are also marked in the figure.As shown in Figure 11, the optimal solution is consistently located in the lower left quarter in each plane.It can be expressed by the fact that the optimal solution searched by TMDSM prefers to possess smaller IPCF values and the highest possible SEA and CFE values.In addition, the other three single best performance values are all located at the extreme positions of the Pareto alternative solution set.It can be seen that the alternative solutions implement the brightest endergonic performance and highest IPCF, which confirms that this optimization algorithm achieves satisfactory diversity and convergence of the population.Due to the conflicting objectives and indicators, it is necessary to balance In Figure 11, the Pareto solutions obtained by the DCNSGA-III algorithm can be observed in the 3D space and the projections are also marked in planes that reflect the relationships between every two target values.As the value of IPCF increases, the SEA also increases.The changes between IPCF and CFE reflect the same trend.It is worth noting that as the SEA value increases, the CFE also increases, indicating that while reducing the cost indicator, the benefit indicator will also be sacrificed.In addition to considering the optimal solution obtained through TMDSM, the other three optimal single performance solutions generated by referencing a single target value are also marked in the figure.As shown in Figure 11, the optimal solution is consistently located in the lower left quarter in each plane.It can be expressed by the fact that the optimal solution searched by TMDSM prefers to possess smaller IPCF values and the highest possible SEA and CFE values.In addition, the other three single best performance values are all located at the extreme positions of the Pareto alternative solution set.It can be seen that the alternative solutions implement the brightest endergonic performance and highest IPCF, which confirms that this optimization algorithm achieves satisfactory diversity and convergence of the population.Due to the conflicting objectives and indicators, it is necessary to balance the relationship between the three in practical engineering.The set of optimal solutions obtained by the methods used in this paper only provides some reference significance for engineers.

The Validation of Optimization Results
To verify the solutions selected are reasonable and reliable, it is effective to demonstrate that the errors between optimization results and the FE model are acceptable.In this section, the numerical model of the obtained optimal solution is reconstructed according to the parameter design variables of Pareto solutions.Table 5 exhibits the alternative solution objective values of Pareto solutions and numerical simulation in order to estimate the proposed optimization system.As shown in Table 5, the relative error between the optimization results and the numerical simulation is within 12%, and the max error occurs in CFE.The comparison dates show that the developed optimization model is reasonable and meets the accuracy requirements.

The Discussion of Optimization Results
In addition, the initial pattern is introduced to compare with the optimal solutions to evaluate the effectiveness of the optimization system.After optimization, the IPCF value decreases by 53.36% and the CFE value improves by 8.68% while the SEA value decreases by 17.88%, which is shown in Figure 12.Although the SEA value has been greatly weakened, the increase in IPCF value and CFE value has compensated for this.It is obvious that tremendous room for amelioration exists in IPCF for SCOD.To sum up, the final result evidences that the optimization method combining the DCNSGA-III algorithm and TMDSM can be designed for specific projects, making it convenient for engineers to refer to.
the relationship between the three in practical engineering.The set of optimal solutions obtained by the methods used in this paper only provides some reference significance for engineers.

The Validation of Optimization Results
To verify the solutions selected are reasonable and reliable, it is effective to demonstrate that the errors between optimization results and the FE model are acceptable.In this section, the numerical model of the obtained optimal solution is reconstructed according to the parameter design variables of Pareto solutions.Table 5 exhibits the alternative solution objective values of Pareto solutions and numerical simulation in order to estimate the proposed optimization system.As shown in Table 5, the relative error between the optimization results and the numerical simulation is within 12%, and the max error occurs in CFE.The comparison dates show that the developed optimization model is reasonable and meets the accuracy requirements.

The Discussion of Optimization Results
In addition, the initial pattern is introduced to compare with the optimal solutions to evaluate the effectiveness of the optimization system.After optimization, the IPCF value decreases by 53.36% and the CFE value improves by 8.68% while the SEA value decreases by 17.88%, which is shown in Figure 12.Although the SEA value has been greatly weakened, the increase in IPCF value and CFE value has compensated for this.It is obvious that tremendous room for amelioration exists in IPCF for SCOD.To sum up, the final result evidences that the optimization method combining the DCNSGA-III algorithm and TMDSM can be designed for specific projects, making it convenient for engineers to refer to.

Conclusions
This study conducts multi-objective optimization of multi-layer nested thin-walled corrugated pipes, namely SCOD, combing the DCNSGA-III algorithm and TMDSM.The DCNSGA-III algorithm is applied to solve the non-dominated Pareto alternative to the design of optimization.The optimal design point is searched by utilizing TMDSM from the alternatives to determine the definitive combination of structural parameters.
The FE model of the SCOD is established based on the basic parameters, and validation is completed.Referring to the error of experiment and simulation results, the accuracy of the FE model is proved as reasonable and effective.In the early stage of optimization work, the design of experimental points is completed using the optimal LHD method, and then variables and target values are fitted according to the obtained experimental point data by utilizing PRS.By computing the four measurements, the surrogate model is validated.An optimal with comprehensive consideration of multiple objectives is obtained by retrieving from the Pareto front.The ultimate optimization result indicates that the IPCF reduces by 53.4% and the CFE increases by 8.7% after optimization.It can be concluded that the optimization method combining the DCNSGA-III algorithm and TMDSM used is effective.
In the future investigation, we will focus on three aspects: (1) developing further experiments and providing more adequate test data to support the validation; (2) applying fuzzy methods and other theories to solve the problem of weight uncertainty to find a more suitable Pareto solution set; (3) studying the impact of ripple amplitude and number on structural crashworthiness and optimize for them.

Figure 1 .
Figure 1.The parameter and structural configuration of SCOD: (a) the sketch for thickness parameters; (b) the sketch for dimension of middle wall.

Figure 1 .
Figure 1.The parameter and structural configuration of SCOD: (a) the sketch for thickness parameters; (b) the sketch for dimension of middle wall.

Figure 3 .
Figure 3.The result of mesh convergence analysis.

Figure 4 .
Figure 4.The schematic diagram of experiment and finite element simulation: (a) the test tube i placed in the test machine; (b) the test and FE model tubes; (c) FE Simulation of SC.

Figure 3 .
Figure 3.The result of mesh convergence analysis.

23 Figure 3 .
Figure 3.The result of mesh convergence analysis.

Figure 4 .
Figure 4.The schematic diagram of experiment and finite element simulation: (a) the test tube is placed in the test machine; (b) the test and FE model tubes; (c) FE Simulation of SC.

Figure 4 .
Figure 4.The schematic diagram of experiment and finite element simulation: (a) the test tube is placed in the test machine; (b) the test and FE model tubes; (c) FE Simulation of SC.

Figure 5 .
Figure 5.The force-displacement curves and collapse modes of SC under experiment test and FE simulation.

Figure 5 .
Figure 5.The force-displacement curves and collapse modes of SC under experiment test and FE simulation.

Figure 6 .
Figure 6.The constituent elements of SCOD and simplified elements.

Figure 6 .
Figure 6.The constituent elements of SCOD and simplified elements.

Figure 7 .
Figure 7.The sensitivity analyses of thickness t on crashworthiness: (a) effect of t on IPCF; (b) effect of t on SEA; (c) effect of t on CFE.

Figure 8 .
Figure 8.The sensitivity analysis of Dmid on crashworthiness.

Figure 7 .
Figure 7.The sensitivity analyses of thickness t on crashworthiness: (a) effect of t on IPCF; (b) effect of t on SEA; (c) effect of t on CFE.

Figure 7 .
Figure 7.The sensitivity analyses of thickness t on crashworthiness: (a) effect of t on IPCF; (b) effec of t on SEA; (c) effect of t on CFE.

Figure 8 .
Figure 8.The sensitivity analysis of Dmid on crashworthiness.

Figure 8 .
Figure 8.The sensitivity analysis of D mid on crashworthiness.

Algorithm 1 :
The DCNSGA-III to solve the framework of CMaOP 1 Appl.Sci.2023, 13, x FOR PEER REVIEW 12 of 23 non-dominated solutions inside the δ-feasible set S1 to approach its PF; (2) the prime selection based on a reference point aims to retain the population multiplicity through providing a group of reference points with good spatial distribution characteristics.

Algorithm 4 :
Select N solutions by non-dominated ranking based on the reference point Input: δ-feasible set S1, population size N, reference points Z Output: Pt + 1

Figure 10 .
Figure 10.Schematic diagram of the entire optimization process.

Figure 10 .
Figure 10.Schematic diagram of the entire optimization process.

4. 4 . 23 4. 4 .
The Pareto Solution of Multi-Objective Optimization A group of Pareto solutions of multi-objective optimization with better convergence, diversity and non-constraint degree is obtained by using the DCNSGA-III algorithm.It should be noticed that SEA and CFE are benefit indicators and IPCF is a cost indicator.In engineering practice, the structure with a larger SEA and CFE and smaller IPCF are preferentially considered.The 3D graphic of Pareto solutions corresponding to objective values is shown in Figure 11, which delivers the relationships of objective values.Appl.Sci.2023, 13, x FOR PEER REVIEW 16 of The Pareto Solution of Multi-Objective Optimization A group of Pareto solutions of multi-objective optimization with better convergence, diversity and non-constraint degree is obtained by using the DCNSGA-III algorithm.It should be noticed that SEA and CFE are benefit indicators and IPCF is a cost indicator.In engineering practice, the structure with a larger SEA and CFE and smaller IPCF are preferentially considered.The 3D graphic of Pareto solutions corresponding to objective values is shown in Figure 11, which delivers the relationships of objective values.

Figure 11 .
Figure 11.The Pareto front in 3D space and its projection in 2D planes: (a) the Pareto solutions are reflected in 3D space; (b) the Pareto solutions are reflected in the plane of IPCF-SEA; (c) the Pareto solutions are reflected in the plane of IPCF-CFE; (d) the Pareto solutions are reflected in the plane of SEA-CFE.

Figure 11 .
Figure 11.The Pareto front in 3D space and its projection in 2D planes: (a) the Pareto solutions are reflected in 3D space; (b) the Pareto solutions are reflected in the plane of IPCF-SEA; (c) the Pareto solutions are reflected in the plane of IPCF-CFE; (d) the Pareto solutions are reflected in the plane of SEA-CFE.

Figure 12 .
Figure 12.The parametric comparison between initial and optimized structures.

Figure A1 .
Figure A1.Second-order interaction response surface of CFE.
Figure A2.Second-order interaction response surface of IPCF.
Figure A3.Second-order interaction response surface of SEA.

Figure A1 .
Figure A1.Second-order interaction response surface of CFE.
Figure A3.Second-order interaction response surface of SEA.

Figure A3 .
Figure A3.Second-order interaction response surface of SEA.

Table 2 .
The comparative result between the experiment test and FE modeling.

Table 2 ,
the IPCF of the test and FE model are 18.34 kN and 17.95 kN, respectively, which are in reasonable agreement.It is obvious that the error between the test and FE modeling is less than 2.5%, which also indicates that the FE numerical simulation has adequate accuracy.

Table 2 .
The comparative result between the experiment test and FE modeling.

Table 3 .
The comparison of SCOD between theoretical prediction and FE model.

Table 3 .
The comparison of SCOD between theoretical prediction and FE model.

Table 3 .
The comparison of SCOD between theoretical prediction and FE model.

Algorithm 2 :
The procedure of tournament selection This process is discussed in Algorithm 3.

Table 4 .
Tolerance and scores of error criteria.

Table 5 .
Comparison of optimization results and FE model.

Table 5 .
Comparison of optimization results and FE model.

Table A1 .
The 60 groups of SEA, CFE and IPCF values of DoE.

Table A1 .
The 60 groups of SEA, CFE and IPCF values of DoE.

Table A1 .
The 60 groups of SEA, CFE and IPCF values of DoE.