Abstract
Genetic algorithm (GA) is a common optimization technique that has two fatal limitations: low convergence speed and premature convergence to the local optimum. As an effective method to solve these drawbacks, an adaptive genetic algorithm (AGA) considering adaptive crossover and mutation operators is proposed in this paper. Verified by two test functions, AGA shows higher convergence speed and stronger ability to search the global optimal solutions than GA. To meet the crashworthiness and lightweight demands of automotive bumper design, CFRP material is employed in the bumper beam instead of traditional aluminum. Then, a multiobjective optimization procedure incorporating AGA and the Kriging surrogate model is developed to find the optimal stacking angle sequence of CFRP. Compared with the conventional aluminum bumper, the optimized CFRP bumper exhibits better crashworthiness and achieves 43.19% weight reduction.
1. Introduction
People nowadays have the increasing awareness of energy conservation and emission reduction and focus more on passive safety performance of cars in collisions. Therefore, the crashworthiness and lightweight design of automobile body structure are regarded as one of the most important parts of vehicle design. Substantial researchers devoted themselves to conduct research studies on improvement of crashworthiness and reduction weight of vehicles. To sum up, there are mainly three ways to achieve these goals: adjustment of material usage, improvement of manufacturing processes, and optimization design of vehicle structure. It is proved that the combination of material adjustment and optimized design of body structure is the most effective way to improve the crashworthiness and realize lightweight [1, 2].
In the latest research studies, carbon fiber-reinforced plastic (CFRP) composites have been increasingly utilized instead of conventional metallic materials in the automotive industry due to their advantageous properties like high strength, high stiffness, lightweight, and excellent energy absorption [3–8]. Kiran et al. proposed the method of using composite instead of metallic materials in automobile chassis and analyzed the static structure of equivalent stress and displacement of composite structural parts [9]. Kim et al. designed carbon fiber-reinforced composites and hybrid fiberglass-reinforced composite hood by finite element analysis and proved that composite hood had better performance in crashworthiness and weight reduction than traditional steel and aluminum hood [4].
As the main load-bearing and energy-absorbing part in the low-speed collision of the car, the bumper system plays a vital role in protecting other parts of the car and the safety of the occupants [4, 10–12]. Therefore, it is of great significance to make the optimization design of the car bumper. Combining genetic algorithm, linear perturbation eigenvalue analysis method, and static RIKS analysis technique, Kim et al. designed a carbon fiber-reinforced composite automotive lower arm, which is characterized by twice the stiffness and buckling strength of the steel lower arm and 50% decline in weight [13]. To find out the suitable material of bumper beam, Kim et al. used hybrid glass/carbon fiber-reinforced composite that had 33% lighter weight and better crashworthiness to replace conventional glass mat-reinforced thermoplastic (GMT) [13].
Recently, the integration of finite element analysis and optimization algorithms has been an advanced way to improve the crashworthiness and lightweight performance of the automotive body structure [14–17]. Fang et al. formed a complete optimized design flow that can be utilized as references by introducing DOE (design of experiment), surrogate model, and multiobjective optimization analysis methods into the crashworthiness and lightweight design of body parts [18, 19]. The selection of optimization algorithm will dominate the quality and accuracy of the result, so a variety of intelligent algorithms are proposed by many researchers [11, 20]; i.e., Sun et al. presented a novel optimization algorithm named ICGA (integer-coded genetic algorithm) and applied it on the configurational optimization of multicell topologies [21]. It can be concluded that increasing intelligent algorithms are applied to multiobjective optimization design of the vehicle, such as particle swarm optimization (PSO), genetic algorithms, and simulated annealing algorithms. The genetic algorithm is one of the most common algorithms, but its limitations, i.e., low convergence speed and premature convergence to the local optimum, are often reported in research studies [22, 23]. Therefore, the improvement of GA is necessary and significant.
In this study, design optimization of the carbon fiber-reinforced composite bumper is carried out with the combination of finite element analysis, the Kriging surrogate model, and multiobjective optimization. The carbon fiber stacking angle sequence is determined to be the optimal design variables, with the purpose of improving the crashworthiness and lightweight performance of CFRP composite bumper. For the multiobjective optimization problem, the selection of the algorithm is the most critical, so an adaptive genetic algorithm (AGA) is presented to generate the best possible design of the carbon fiber stacking angle sequence.
2. Adaptive Genetic Algorithm
2.1. Limitations and Improvement of GA
Genetic algorithm (GA) is widely employed to search the global optimal solution in optimization problems; however, it has some defects in solving complicated optimization problem [24], such as the premature convergence to the local optimum and much time consuming. The crossover probability and mutation probability greatly determine the quality of solutions and the convergence speed of GA [25]. The moderately large values of promote the extensive recombination of schemata, while small values of are necessary to prevent the disruption of the solutions. Instead of fixed values of and , AGA utilizes the population information in each iteration and adaptively adjusts the and to maintain the population diversity as well as to sustain the convergence capacity. And the roulette wheel selection method is employed as the selection operator.
It is essential to identify whether the GA is converging to an optimum. The difference in the average and maximum fitness values of the population [26], , is used to be a yardstick for detecting the convergence of GA. And when GA converges to a local optimum, the value of decreases. Therefore, the values of and are supposed to be varied adaptively depending on the value of to prevent premature convergence of the GA to local optimum. In addition, to preserve the “good” solutions of every populations, lower values of and should be set for high fitness solutions and higher values of and should be set for low fitness solutions [27, 28]. While the high fitness solutions aid in the convergence of the GA, the low fitness solutions prevent the GA from getting stuck at a local optimum. So the values of should depend on the fitness values of the solution too. Similarly, the values of should also depend on the fitness values of both the parent solutions. The closer the is to , the smaller the should be. From the above, the adaptive value of and can be expressed as follows:where and are less than 1.0 to constrain and to the range , is the individual fitness, is larger of the fitness values of the solutions to be crossed, and and is the average and maximum fitness values of the population, respectively.
However, the adapting crossover probability and mutation probability defined by equations (1) and (2) have a drawback: if , the values of and are zero that will lead to premature convergence to the local optimum. As an effective method to overcome this drawback, improved equations to calculate and are proposed in this paper; they can be corrected aswhere and are the added parameters that can avoid the condition of , which will avoid premature convergence to the local optimum. The process of adaptive algorithm is depicted in Figure 1.

2.2. Performance Test of Adaptive Genetic Algorithm
In this study, two typical test functions are employed to test and evaluate the adaptive genetic algorithm and genetic algorithm and thus verify their optimization performance.
2.2.1. Schaffer Function
The function is a multipeak function with multiple local maxima, which is mainly utilized to test the optimal accuracy of the algorithm [29]. The function can be expressed as equation (4), and then, the optimal problem for the test is constructed as follows:
The global maximum of Schaffer function is , and there are countless local maxima at the value of 0.99028 around the global maximum. The graphical representation of function and contour line is shown in Figure 2.

(a)

(b)
2.2.2. Rosenbrock Function
This function is a nonconvex function, which is widely used as a performance test function for optimization algorithms [30]. The function can be expressed as equation (4), and then, the optimal problem for the test is constructed as follows:
The global minimum is inside a long, narrow, parabolic-shaped flat valley. To find the valley while trivial and to converge to the global minimum are also difficult. The graphical representation of function is shown in Figure 3.

The genetic algorithm and adaptive genetic algorithm are tested by the Schaffer function, respectively, for comparison. To provide a fair comparison and test the robustness of the AGA and GA, the population size is set at 10, 15, and 20, respectively, and the maximum number of iterations for two algorithms is fixed at 100. Referring to the previous research [26], the values of and can be set as 1 and 0.5 for GA, while the values of and are also set as 1 and 0.5 for AGA. And the values of and are 0.5 and 0.01 to constrain the and to the best range [31]. Each test function is tested 30 times in MATLAB, and then, the optimal, worst, and average values of fitness are obtained in AGA and GA. The performance comparison between AGA and GA under different test functions is shown in Table 1, and the fitness curves are plotted in Figures 4 and 5.

(a)

(b)

(c)

(a)

(b)

(c)
From Table 1 and Figures 4 and 5, the performance of AGA and GA can be concluded as follows:(1)In terms of Schaffer function, both AGA and GA can jump to the local optimal solution and reach the global optimal solution at value of 1, but the convergence speed of AGA is faster than GA(2)For Rosenbrock function, the ability of searching global optimal solutions of AGA is better than GA; however, the convergence speed of AGA and GA is close
All in all, compared with the GA, the novel AGA proposed in this paper improves in both convergence speed and the ability to search global optimal solution.
3. Finite Element Model
In the whole process of a car frontal collision, the deformation of the car mainly concentrated in the front cabin while the deformation of the components behind a-pillar is small or almost no deformation occurs. Therefore, the finite element model of the body-in-white car is built as a reduced model, which uses a particle to replace the inessential components behind a-pillar. The size of the mesh is vital to the accuracy of CAE [32]; therefore, the element size of 5 × 5 mm is determined to be sufficient for the front cabin, because this part has a great influence on crashworthiness, and 10 × 10 mm is determined to the other parts. After preprocessing in Hypermesh software, the FE model is modeled by quadrilateral shell elements and the total number of elements is 226503, in which the triangular elements only account for 3.8%, and it is shown in Figure 6.

Nonlinear finite element Ls-Dyna is utilized to simulate the crash behavior of automobile in this study. In the initial model, aluminum is used as the material of the bumper beam. Based on the latest C-NCAP (China-New Car Assessment Program) [33], the collision velocity of the car is determined to be 50 km/h, the acceleration is 9.8 m/s2, and the total time is 0.12 s.
Then, the material of bumper beam is replaced by the GFRP composite. The bumper beam is the main frontal part of automotive structure during a frontal crash and is generally subjected to bending load. Therefore, stiffness of the initial bumper beam is supposed to be the same as the GFRP bumper beam [34], and the equal stiffness equation can be expressed as equation (9):where and are the inertia moment and elastic modulus of CFRP composite bumper, respectively, while the and are the inertia moment and elastic modulus of the aluminum bumper. And I can be calculated by
Then, the thickness of CFRP composite bumper can be calculated bywhere is the thickness of the CFRP composite bumper and is the thickness of the aluminum bumper.
The thickness of the aluminum bumper beam is 1.80 mm; therefore, the thickness of the CFRP bumper beam is 1.92 mm calculated by using equation (10). In Hypermesh software, the material model 54 (CFRP) is selected instead of aluminum in the bumper beam. The failure criteria of CFRP were introduced in the research [35]. After replacing the aluminum bumper beam by the CFRP bumper beam, the weight of the bumper beam is reduced from to that achieves 43.19% weight reduction. And the material properties of CFRP are listed in Table 2 [36].
4. Multiobjective Optimization of CFRP Composite Bumper
4.1. Definition of the Optimization Problem
The structures of composite have great flexibility in design, because they can be customized in different ways to meet specific design requirements. By selecting the best parameters such as fiber and matrix material, fiber orientation, or layer thickness, the fatigue strength of the laminated composite can be significantly improved, thereby improving the structural properties [37, 38]. The CFRP composite studied in this research can be regarded as a multidirectional laminate plate composed by a single-layer plate that is laminated in different directions. The microstructure of CFRP can be shown in Figure 7.

The stiffness of CFRP is greatly affected by the different stacking angle sequences [39, 40]; therefore, further optimization is carried out to find the optimal stacking angle sequence, with the purpose of improving the crashworthiness and reducing the quality of bumper.
4.1.1. Design Variables
The thickness of the bumper beam of the above model is 1.92 mm, and CFRP laminates are designed as symmetric based on the center of the laminate and composed of ten layers in total, so the thickness of each ply is 0.192 mm. Then, the stacking angle of each layer is determined to be the design variables; therefore, five independent design variables can be considered as [x1/x2/x3/x4/x5]s, which represent the centrosymmetric stacking angle sequence of CFRP laminates. And the value of design variable x can be taken as [0/±45/90] based on the global coordinates. The schematic diagram of 10-plies CFRP based on the global coordinates is vividly depicted in Figure 8.

4.1.2. Constraints
The peak impact force F, the intrusion L, and the energy absorption Q of bumper beam are considered to be the constraints of optimization design, and the values of them are supposed to be lower than those of initial aluminum model.
4.1.3. Objectives
According to C-NCAP, upward displacement automotive steering column S1 and rearward displacement automotive steering column S2 can be used to evaluate the crashworthiness and they are shown in Figure 9. And the purpose of optimization design is to minimize them.

The optimization problem can be defined mathematically as follows:
4.2. Kriging Model
Design optimization of crashworthiness and lightweight is a nonlinear problem that usually costs much time. In the optimization problem above, the objective function and constraint function are complicated nonlinear function. As an effective approach, the Kriging surrogate model is employed to represent multimodal and nonlinear functions in this paper [41–43]. It can be expressed aswhere is a constant approximate model, and it is similar to the polynomial response surface that provides design space for a global model and is a stochastic process with mean zero and variance and it creates localized deviations to assist the Kriging model generating interpolation among sample points. The covariance matrix of is expressed aswhere is a correlation matrix of Gauss–Markov theorem.
In this study, the Kriging surrogate model is used to approximate F, L, Q, S1, and S2 functions for each case. Design of experiment (DOE), as a widely used technique to address how to select the minimum number of training points for mapping the entire design space properly, is the first step of surrogate modeling technique. The optimal Latin hypercube sampling (OLHS) [44] can efficiently generate uniformly distributed sampling points, so it is applied to generate 120 baseline designs. Then, the F, L, Q, S1, and S2 at these training points are calculated by Ls-Dyna software. All the calculation results are based on the global coordinate. Based on the simulation results, the Kriging models of the above optimal problem are constructed. To calculate easily, 1, 2, 3, and 4 are substituted for 0°, −45°, 45°, and 90°. The 120 baseline designs are shown in Table 3.
The fitting accuracy of the Kriging surrogate model is crucial to the result of an optimal problem, and the R2, RMAE (relative maximum absolute error), and RMSE (Root-mean-square error) are employed as the criterion to evaluate the accuracy of it [45]. They can be calculated bywhere is the mean response value of test points, is the predicted response value and is the real test response value and N is the number of test samples. The ideal surrogate model is supposed to meet the requirements that , , and . After evaluating the Kriging surrogate model in I-sight software, it was obtained that all the values of R2, RMAE, and RMSE meet the requirements, which indicates the Kriging model in this paper has high accuracy and reliability.
4.3. Implement of Adaptive Genetic Algorithm
The adaptive genetic algorithm is employed in the optimization problem mentioned in Section 4.2. The stacking angles 0°, −45°, 45°, and 90° of each layer in CFRP are the genes, and they are encoded as 1, 2, 3, and 4. The stacking angle sequence [x1/x2/x3/x4/x5] is the chromosomes. The Kriging surrogate models of S1 and S2 are the fitness functions while the Kriging surrogate models of F, L, and Q are the constraint functions. The values of and are 1 and 0.5 for GA, and the values of , , , and are 1, 0.5, 0.5, and 0.01, respectively, for AGA. To compare the optimal result between AGA and GA, the population size is set as 20 and maximum iteration is 2000 in both GA and AGA.
In terms of AGA, 24 Pareto solutions are obtained and 14 Pareto solutions are obtained using GA. The Pareto-optimal fronts of AGA and GA are shown as Figure 10. In practical engineering requirement of crashworthiness, the upward displacement automotive steering column S1 is as important as the rearward displacement automotive steering column S2. Therefore, in this paper, the Pareto solution that has the minimum sum of S1 and S2 is selected as the optimal solution. Then, the optimal solutions of AGA and GA and the comparison to the initial aluminum model are shown in Table 4.

It can be concluded from Figure 10 that a more widespread and even distribution nondominated solution set is obtained by AGA. Both the number and the distributivity of Pareto solutions of AGA are better than those of GA. And the results also indicate that the crashworthiness indicators of the AGA-optimized bumper have slightly improvement than those of GA.
4.4. Verification of the Results to FE Model
To verify the optimal result of AGA in practical application, the best stacking angle consequence is applied in the FE model and then calculated in Ls-Dyna. The whole process of collision can be considered as the conversion of kinetic energy to internal energy, and the C curve of total energy basically maintains a horizontal shape, which represents that the total energy satisfies energy conservation. The maximum of the hourglass energy and slip energy should be less than 5% of the maximum kinetic energy, which means that they are within the acceptable range, which is depicted in Figure 11. It is a necessary precondition of collision simulation, and the optimized FE model meets the requirements. And the optimal results of the Kriging surrogate model are compared to the simulation results by FEA. The accuracy of surrogate model is shown in Table 5.

It can be concluded from the above results that the Kriging surrogate model is of high accuracy and the optimization algorithm is effective. And the comparison of crashworthiness performance in Ls-Dyna is shown in Figure 12.

(a)

(b)

(c)
5. Conclusion
In this study, a novel multiobjective optimization method considering the adaptive genetic algorithm and Kriging surrogate model is proposed for improving crashworthiness and reducing weight of the automotive bumper. The CFRP material is employed in the bumper beam instead of aluminum, and the optimal stacking angle sequence of each layer in CFRP is obtained by AGA. The following conclusion can be summarized:(1)Compared with the genetic algorithm, the adaptive genetic algorithm proposed in this paper has higher convergence speed and stronger ability to search global optimal solutions. And applying them to the optimization problem in crashworthiness, the comparison results indicate the crashworthiness of the CFRP bumper optimized by AGA is better than that of GA.(2)After the use of CFRP material in the bumper beam and optimization by AGA, the new bumper achieves 43.19% weight reduction and exhibits better crashworthiness, compared to the initial aluminum bumper.
Furthermore, better optimization algorithm and more design variables of material will be considered into the research of the optimization of automotive bumper in the future.
Data Availability
No data were used to support this study.
Disclosure
The views expressed in this paper are those of the authors and do not necessarily reflect the official views of affiliations.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the research fund of the Shanghai Automotive Industry Technology Development Foundation (1744), National Natural Science Foundation of China (Grant no. 51809168), and Young Teacher Initiation Program of Shanghai Jiao Tong University (17X100040060).