Next Article in Journal
Human-like Attention-Driven Saliency Object Estimation in Dynamic Driving Scenes
Previous Article in Journal
Acoustic Feature Extraction Method of Rotating Machinery Based on the WPE-LCMV
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization Design of Automotive Body Stiffness Using a Boundary Hybrid Genetic Algorithm

1
Ji Hua Laboratory, Foshan 528200, China
2
State Key Laboratory of Automotive Simulation and Control, Changchun 130025, China
*
Author to whom correspondence should be addressed.
Machines 2022, 10(12), 1171; https://doi.org/10.3390/machines10121171
Submission received: 12 November 2022 / Revised: 3 December 2022 / Accepted: 5 December 2022 / Published: 6 December 2022
(This article belongs to the Section Vehicle Engineering)

Abstract

:
At the conceptual design stage, it is critical to use appropriate structural analysis and optimization methods. The thin-walled beam transfer matrix method (TBTMM) is adopted to establish the mathematical model of the simplified vehicle body-in-white (BIW) structure in this paper and compare it with the results of the finite element method (S-FEM) to verify the approach. In addition, on the basis of the boundary simulation genetic algorithm (BSGA) and local search procedure, a boundary hybrid genetic algorithm (BHGA) is proposed. BHGA is benchmarked on 20 test functions and is compared with current meta-heuristic algorithms to prove its effectiveness and universality. Finally, considering the bending and torsional stiffness constraints, BIW conceptual model is lightweight and designed with an optimizer.

1. Introduction

Compared with other stages, the conceptual design stage has a higher degree of design freedom and various requirements that are easier to meet, which is crucial for vehicle innovation, cost saving, and design cycle shortening [1]. The bending stiffness, torsional stiffness, and NVH performance of the body-in-white (BIW) structure are critical to the safety and comfortability of the vehicle, while the weight performance of the BIW has an effect on cost saving. There is a direct correlation between these two performances; that is, the reduction of vehicle weight requirement will lead to an overall reduction of stiffness and NVH performance. In order to obtain the optimal BIW structure that meets various requirements in the conceptual design, a lightweight design is selected to coordinate vehicle weight (objective function) and its performance (constraints). The performance of bending and torsional stiffness will be considered in this paper.
It is common to establish the mechanical model of the BIW by employing the finite element method (FEM). Bai et al. [2] described a simplified finite element model to provide early-stage predictions of a detailed model. A concept CAE modeling approach based on FE models was presented by Donders et al. [1] to analyze and optimize the structural behaviors of the vehicle BIW. In addition, Mundo et al. [3] proposed a similar approach for replacing beam structures and joints in vehicle BIW. However, few CAD data are available at the conceptual design stage. Meanwhile, most approaches to FEM require a time-consuming FE model and can not obtain clear mathematical relationships between beam section properties and structural performances. Therefore, Qin et al. [4] developed an object-oriented MATLAB toolbox based on the exact transfer stiffness matrix method to calculate the static and dynamic performances. Later, Liu et al. [5] developed a mathematical method based on the reverberation ray matrix method to promote the conceptual design. The vehicle BIW can be regarded as consisting of thin-walled beams, so the warping deformation (longitudinal displacement) in the longitudinal direction of the beam due to torsion should be considered [6]. However, the aforementioned studies are all based on the traditional beam theory, according to which the warping of the beam is not considered. Zhong et al. [7] also pointed out that warping has a great impact on the performances of thin-walled frame structures and developed the thin-walled beam transfer matrix method (TBTMM). Therefore, the TBTMM is adopted in this paper to establish the relationship between the section properties and structural mechanical performances so as to improve the calculation accuracy.
In order to realize the lightweight design of vehicle body-in-white (BIW) at the conceptual design stage, a boundary hybrid genetic algorithm (BHGA) was proposed to solve the engineering-constrained optimization problems. In the next section, BIW mathematical models for bending and torsional stiffness will be established. In Section 3 and Section 4, the development and validation of BHGA will be introduced. Finally, the lightweight design is achieved by using BHGA.

2. Formulation of Vehicle BIW Mathematical Model

2.1. Formulation of Cross-Sectional Properties

In this paper, a simplified BIW model with rectangular thin-walled beams is adopted to analyze the mechanical performances, as shown in Figure 1. Different from the cross-section properties of general beam theory (the inertia moments I y and I z , the cross-sectional area A , and the torsional constants J ), the thin-walled beam theory involves fifteen cross-sectional properties, i.e., the angle of roll ϕ , the cross-sectional area A , the shear areas A s y and A s z in the y- and z-directions, the hybrid shear area of the cross-section A s y z , the sectorial static moments S s y and S s z about the y- and z-axes, the warping torsion moment of inertia I p , the St.Venant torsional constant corresponding to Bredt’s shear stress I B , the St.Venant torsional constant corresponding to St.Venant’s shear stress I s , the inertia moments I y and I z , the sectional moment of inertia I w , and the coordinates of the shear center y s and z s relative to the centroid. Thus, the thin-walled beam cross-section is defined as in Figure 2, in which the right-hand orthogonal coordinate system and the s-coordinate system are adopted, the latter being along the midline of the cross-section. Moreover, the width (b), height (a), and thickness (t) are defined as design variables. The rectangular thin-walled beam is used in this paper; that is, the shear center S coincides with centroid C and ϕ is equal to zero. The calculation formula for section properties has been deduced in detail in the previous paper [7].

2.2. TBTMM Mathematical Model of BIW Structure

As shown in Figure 3, the mathematical model of the BIW conceptual structure consists of 53 beam members and 36 joints. The ith beam is identified by ‘ ’, and the jth joint is identified by ‘j’. According to thin-walled beam theory, the governing equations of beam members can be expressed as follows
d S ( x ) k d x = H k S ( x ) k
where the S ( x ) k is the state vector of kth member in the place of x containing seven displacement fields and seven force fields. H k is expressed as the state function of kth member.
The Laplace transform is used to solve Equation (1).
S ( x ) k = L 1 ( ( X E H k ) 1 ) S ( 0 ) k
where L 1 ( ) denotes the inverse Laplace transform and X is the Laplacian operator.
Then, the transfer matrix T k can be obtained by substituting the length of beam l k into Equation (2), i.e., T k = L 1 ( ( X E H k ) 1 ) x = l k , and the equation can be rewritten as Equation (3). The T k is a constant matrix containing all fifteen cross-sectional properties.
S ( l k ) k = T k S ( 0 ) k
When the transfer matrices of beams have been determined, the mathematical model for each joint can be assembled as follows
C ¯ i P ¯ i = f ¯ i
where C ¯ i denotes the ith joint coupling matrix; P ¯ i contains all the state vectors defined at ith joint; f ¯ i is the matrix about external forces and moments at the ith joint. In general, the boundary conditions of the structures are homogenous; that is, each joint coupling matrix has half as many rows as columns. Moreover, each state vector contained in P ¯ i is one of the previously defined state vectors S ( l k ) k or S ( 0 ) k .
The global mathematical model of the BIW conceptual structure can be assembled as conventional FEM by combining joint coupling matrices.
C P = f
where C is a 742 × 1484 matrix, and C = d i a g [ C ¯ 1 , C ¯ 2 , C ¯ 36 ] denotes the global joint coupling matrix; P = [ P ¯ 1 , P ¯ 2 , P ¯ 36 ] T contains all the state vectors, and is a 1484 × 1 matrix; f = [ f ¯ 1 , f ¯ 2 , f ¯ 36 ] T is the matrix about external forces and moments, and is a 742 × 1 matrix.
From Equations (3) and (5), P can be written as
P = R T g l o b a l Y
where Y is a 742 × 1 matrix, and Y = [ S ( 0 ) 1 , S ( 0 ) 2 , S ( 0 ) 53 ] T contains all the input state vectors defined at the left ends of the elements; R is a 1484 × 1484 matrix used to rearrange the order of the state vectors; T g l o b a l is the total transfer matrix and can be expressed as follows
T g l o b a l = T 1 T 2 T 53 I 1 I 2 I 53
where the transfer matrix T i   ( i = 1 , 2 , …53 ) is a 14 × 14 matrix and I i   ( i = 1 , 2 , …53 ) is an identity matrix, so the T g l o b a l is a 1484 × 742 matrix.
Then, the complete description of the BIW structure with respect to the state vectors Y and the external forces f can be rewritten as follows
C R T global Y = f
The C R T g l o b a l in Equation (8) is a 742 × 742 square matrix with respect to cross-sectional properties. Moreover, the external forces are set as the known values with respect to the boundary and load conditions. Then, Equation (8) may be formally solved by
Y = [ C R T g l o b a l ] 1 f

2.3. Static Load Cases and Boundary Conditions

In this paper, the loads and boundary conditions of the BIW static analysis are shown in Figure 4. For the bending condition, the external forces F b are set as 1668 N along the Z- direction act on joints 10 and 28, respectively, and the constraint points are set on joints 6, 13, 24, and 31, respectively. The digital 1, 2, and 3 in the triangle region represent that the displacements in the global X-, Y-, and Z-directions are restrained to zero. For the torsion condition, the external forces F t are set as 1668 N act on joints 13 and 31, and the constraint points are set on joints 6, 24, and 36, respectively.
Substituting the load cases and the boundary conditions into Equation (9) would yield the coupling equations of the whole structure. The maximum displacements of load directions Δ l b ( Δ r b ) and Δ l t ( Δ r t ) can be obtained by solving the coupling equations, respectively. And the bending stiffness K b and the torsion stiffness K t can be defined as follows:
K b = 4 F b Δ l b + Δ r b K t = F b L c α
where L c is the distance between the two load points; α is the twist angle which can be calculated as follows:
α = 180 ( Δ l t Δ r t ) π L c

2.4. The Accuracy Verification

In order to verify the accuracy of the developed mathematical model of the BIW structure, the results of bending and torsion stiffness obtained by TBTMM are compared with the results obtained by FEM. The FEM analysis results are obtained by using the program Hypermesh, while the TBTMM analysis results are obtained through the MATLAB code. As listed in Table 1, it is reasonable to establish the vehicle BIW mathematical model by using TBTMM, as the relative error values are much lower than 20% [8].

3. The Development of BHGA

GA is a traditional evolutionary algorithm, which is a stochastic search technique based on a series of possible solutions. Proposed by Holland in 1975 [9], GA has been extensively used in engineering and industry problems related to linear inequality constraints, nonlinear inequality constraints, equality constraints or unconstraints. Thus, it is reasonable to solve the constrained nonlinear optimization problem in equation (9) by using GA. As an unconstrained search technique, constrained problems have traditionally been challenging problems for GA. Several techniques of constraint handling have been developed: special representation and operator methods, penalty methods, separation of objective and constraint methods, repair methods, and hybrid methods. The most common way to introduce constraints in genetic algorithms is the penalty method, which punishes infeasible solutions by reducing their fitness values. Ersavas et al. [10] and Paszkowicz [11] use the static penalty method and dynamic penalty method for constraint optimization, respectively. However, as it is difficult to set the penalty coefficient properly, its performance is not always satisfactory. Considering this obstacle, Lin [12] came up with a rough penalty GA for constraint optimization; nevertheless, it requires another set of parameters to tune the penalty coefficients automatically. The main idea of the special representation and operator methods is to develop special representation schemes to tackle a certain problem for which generic representation schemes might not be appropriate. Koziel and Michalewicz [13] proposed a ‘homomorphous map’, in which they transformed the whole feasible region into a different shape that was easy to optimize. However, the implementation of the algorithm is more complex, and the experiments reported require a large number of fitness function evaluations. Although it is efficient for some intended applications, it can sometimes be difficult or even impossible to develop a special representation. The separation of objective and constraints methods treats the constraints as an objective function so that the original single-objective constraint optimization problem becomes a multi-objective unconstrained optimization problem, to which we can apply any multi-objective optimization techniques. Zhou et al. [14] developed a ranking procedure in accordance with the Pareto strength concept for the bi-objective problem, but with the constraints increasing, the objective function becomes complicated. As pointed out by Runarsson and Yao [15], the multi-objective techniques are difficult to find feasible solutions since most of the time is spent searching infeasible regions. While GA is coupled with another technique (e.g., another heuristic or a mathematical programming approach) to form hybrid methods, the new methods generally require several parameters to work properly, just like penalty methods. The main idea of the repair methods is to transform an infeasible solution into a feasible one, which can reduce the search space by using a repair technique. Furthermore, no special operators or modifications of the fitness function need to be considered in this case. Salcedo et al. [16] proposed a concept of a hybrid genetic algorithm in which the local search (LS) procedure is used as a constraint-handling technique. Later, Li et al. [17] proposed a boundary simulation genetic algorithm (BSGA) to address inequality constraints for GAs and developed a series of genetic operators that would abandon or repair infeasible individuals produced during the search process. However, it was not specified whether the infeasible solutions were abandoned or repaired, and it may not work properly for problems with disconnected feasible regions. Coello [18] has emphasized that a desired constraint-handling technique should be general and incorporate knowledge about the domain, efficiency, etc. Based on the BSGA and hybrid genetic algorithm, a boundary hybrid genetic algorithm (BHGA) that could be applied effectively to engineering is proposed in this paper. In general, the method proposed in this paper has several features, which will be discussed in detail in the following chapters.
  • The BHGA randomly selects individuals from the boundary point set as the feasible initial population and performs a global search (GS) using the GA;
  • Perform the elitist strategy and adaptively tune crossover and mutation operators;
  • The LS procedure is used to handle constraints.

3.1. Generate the Initial Population

In this paper, we mainly concentrate on the problems of inequality constraints. Thus, the mathematical form of the optimization problem could be formulated as the following minimum optimization problem:
Minimize: f ( x )
Subject to:
g i ( x ) 0 ,                     i = 1 , 2 , , m x j l     x j x j u ,     j = 1 , 2 , , n                      
where x = ( x 1 , x 2 , x n ) is the vector of design variables; m and n are the number of constraints and design variables, respectively; x j l and x j u denote the lower and upper bounds of x j , respectively.

3.1.1. The Calculation of Boundary Points

Isaacs et al. [19] indicated that the optimal solutions to the constrained optimization problems are usually spread along the constraint boundary. The BSGA [18] is proposed to solve the constrained optimization problem based on the binary search method, but the binary search method would not search the maximum constraint boundary and may fall into the local optimum in the constraint boundary. Therefore, the reverse binary search method and LS strategy are developed in the BHGA, and it will be discussed in detail later. The flowchart for the main process of the generation of feasible regions is illustrated in Figure 5.
The first step is to generate FP feasible points. Compared with unconstrained and simple constrained problems, it becomes relatively complicated to generate feasible points for complex constraints. In order to improve the search efficiency, the GA method is still used to generate the feasible points, and the main process is illustrated in Figure 6. Moreover, a real coding representation scheme is adopted in this paper. Initially, an empty set is developed for storing the feasible points and then generates the initial population, including PS individuals, randomly. The initial individuals satisfying the constraints are put into the feasible set, and the number of feasible points NP is recorded. The new populations are generated by selection, crossover, and mutation operators of GA [20], and the feasible points set is updated based on the feasible points in the new populations until the number of feasible points reaches FP.
Since the GA is used to obtain a set of feasible points instead of an optimal solution, the fitness function needs to be modified. The aim is to reduce the differences between the individuals but not convergence, and the constraint violation status for each individual is expressed as Equation (13). The fitness value of every individual can be calculated by summing the constraint violation degree for each of them.
q i ( x ) = g i ( x )     i f   g i ( x ) > 0 0                       i f   g i ( x ) 0      
where q i ( x ) represents the ith constraint violation degree, while g i ( x ) represents the ith constraint.
The following step is to generate the infeasible points. Define the expansion region and transform the feasible region ( S f e a s i b l e = { x R n | | x j l x j x j u   and   g i ( x ) 0 ,     f o r   j = 1 , 2 , , n } ) into infeasible region S :
S = { x R n | x ˜ j l x j x ˜ j u ,   f o r   j = 1 , 2 , , n }
The x ˜ j l and x ˜ j u are the lower and upper extension region bounds of the jth design variable and are defined as Equation (15).
x ˜ j l = x j min ( x j u x j l ) x ˜ j u = x j max + ( x j u x j l )
where x j min and x j max are the minimum and maximum values of jth variable about the feasible point, respectively.
As shown in Figure 7, generate IP individuals in the extension region and set one of the variables in each individual to its lower or upper bound.
Figure 8 illustrates the feasible point a moves to the constraint boundary by a large step size in the binary search method of BSGA, and the process ignores the possibility that there is a constraint boundary between the infeasible points c and d. Thus, in order to increase the search range of the feasible points (i.e., to find the maximum boundary of the feasible points) and reduce the step size, the reverse binary search method is proposed to find the boundary of the feasible region shown in Figure 8. Compared with the binary search method, the latter method requires ε and Δ parameters to control the search accuracy. Moreover, the specific steps of the reverse binary search method are as follows:
  • Step 1: Choose a feasible point a and an infeasible point b, then go to step 2.
  • Step 2: Calculate the middle point c between a and b, then go to step 3.
  • Step 3: If c is feasible, a = c, then go to step 6; otherwise, calculate the middle point d between c and b, then go to step 4.
  • Step 4: If d is feasible, a = d, then go to step 6; otherwise, go to the next step.
  • Step 5: Calculate the distance between b and d; if it Δ , b = d, then go to step 2; otherwise, c = d, then go to step 3.
  • Step 6: Calculate the distance between a and b; if it ε , terminate; otherwise, go to step 2.

3.1.2. The Initial Population

A set of boundary points have been generated in the previous section, and then the initial population is randomly selected from the boundary points in the following steps. Initially, popsize individuals are randomly generated in the extension region, and the distances from each individual to all boundary points are calculated. Then the nearest boundary points are selected as the initial population individuals.

3.2. The Definition of GA Operators

The GA operators are mainly selection, crossover, and mutation operators. The selection operator is utilized to select some good individuals from the current population as parents to generate offspring. In this paper, the well-known roulette wheel selection operator is adopted [20]. Due to the randomness of the selection strategy, the elitist strategy is implemented by reinsertion to prevent the good individuals from being abandoned.
In the crossover and mutation operations, the setting of the crossover probability (Pc) and the mutation probability (Pm) is crucial to the generation of new individuals. In the BSGA method, the use of fixed crossover and mutation probabilities may lead to local optimum. The higher the value of Pc, the quicker the new individuals will be introduced into the population. However, as Pc increases, individuals can be disrupted faster than selection can exploit them. Similarly, if the values of Pm are too small, it is not easy to generate new individuals. However, the larger values of Pm transform the GA into a pure random search algorithm. In order to improve the performance of the BHGA, the fixed values of Pc and Pm can no longer meet the dynamic performance of the algorithm, so the adaptive crossover and mutation probabilities are proposed in this paper.
P c = P c 1 ( P c 1 P c 2 ) ( f f a v g ) f max f a v g ,     f f a v g P c 1                                                                 ,       f < f a v g
P m = P m 1 ( P m 1 P m 2 ) ( f f a v g ) f max f a v g ,     f f a v g P m 1                                                                   ,       f < f a v g
where f max is the maximum fitness value of the current population, while f avg is the average fitness value of the population and f is the individual fitness value; P c 1 and P c 2 represent the upper and lower limits of the crossover probability, while P m 1 and P m 2 represent the upper and lower limits of the mutation probability.

3.3. The LS Strategy

There are two main operators proposed to handle the infeasible individuals during the search process of BSGA: one is to regenerate it again until a feasible individual is obtained, but this will change the characteristics of the original individual; the other is to repair the infeasible individuals with the binary search method, but this will also cause the algorithm to converge to the boundary local optimum.
Thus, the LS procedure occurs as a method of constraint handling in the process of infeasible offspring approaching the feasible boundary in the BHGA. The current parent Pr and the infeasible individuals Or generated by the crossover or mutation operation are selected to perform a nonlinear search along the P r O r direction and are defined as Equation (18). As shown in Figure 9, the dashed blue box represents the process of feasible individual Pr searching for the boundary, while the red dashed box is the process of searching for the discrete feasible domain by the reverse binary search method. The main LS procedure is as follows:
Step 1: Select the current parent Pr and the infeasible individual Or, then go to step 2.
Step 2: Generate the repaired individual P o by Equation (18), and operate as follows
  • If the P o is feasible, calculate the objective function value of P o and the distance between P o and O r . If the current objective function value of P o is optimal compared to other individuals or the distance the set value Δ , then P o is output as the repaired individual; otherwise, P r = P o , and then go to step 2.
  • If the P o is infeasible, then repair the infeasible individual with the reverse binary search method and go to step 2.
    P O = P r + c r ( O r P r )
    where c is used to tune the LS step and r is the random number which is uniformly distributed in the interval [ 0 , 1 ] .

4. Experimental Study and Discussion

In order to verify the numerical efficiency of the BHGA, 20 benchmark functions (6 unconstrained, 11 constrained, and 3 engineering-constrained problems, all of which are minimization problems) are adopted. The algorithm is running under Windows 7 Ultimate, and the code is programmed and compiled in MATLAB R2014a. All the details of the benchmark functions are listed in Table 2, where N is the number of design variables; LI and NI are the number of linear inequality and nonlinear inequality constraints, respectively; ρ is the estimated ratio between the feasible region and the search space; BO represents the ratio between the number of constraints at the boundary with the total number of constraints. It can be seen from Table 2 that the optimal solution is easily located at the constraint boundary, so it is meaningful to propose the BHGA and start searching from the boundary point. All the used corresponding parameters of the BHGA are listed in Table 3.

4.1. Benchmark of Unconstrained Functions

For the unconstrained problem, six composite benchmark functions with several randomly located global and deep local optima of CEC’2005 [21] are considered in this paper. Considering these unconstrained optimization problems is mainly to verify whether the BHGA is more advantageous than the BSGA for optimal local problems and whether it can be used to deal with unconstrained optimization problems.
The statistic of optimization results of the BHGA and the comparison algorithms are plotted in Table 4. 10 comparison algorithms are adopted, which are BSGA [17], IGA [22], SSA [23], GOA [24], WOA [25], GWO [26], PSO [27], GSA [28], MVO [29], and HS [30]. For all the algorithms, the same population size and iteration number equal to 30 and 1000 have been utilized, and all algorithms run 30 times. The optimization results of the comparison algorithms use the results of the literature [22].
From the ranking of the average results in Table 4, it can be seen that although the BHGA has not obtained an optimal solution, the stability of the solution for different unconstrained problems is generally acceptable. Moreover, the BSGA is more likely to fall into the optimal local solution for composite functions, and the proposed BHGA is more effective.

4.2. Benchmark of Constrained Functions

The 11 constrained problems (G01, G02, G04, G06-G10, G12, G19, G24) selected from the CEC’2006 [31] are adopted to verify the constrained optimization ability of the BHGA. The optimal results (best, worst, average, and standard deviation) obtained by the BHGA are listed in Table 5, and the algorithm runs 30 times with 240,000 function evaluations. The exact/near-optimal results calculated by the BHGA are highlighted in boldface, and 8 of the 11 benchmark functions (G01, G02, G04, G06, G08, G09, G12, G24) obtained the known optimal results. Only G10 did not obtain the exact/near-optimal result, but it is also hard for G10 to get the optimal result through other algorithms. Table 6 shows the comparative results of benchmark functions obtained by the BHGA and other comparison algorithms (d-DS [32], BSGA [17], HTS [33], BBO [33], TLBO [33], GA [34], PSO [35], DE [36], and ABC [37]). As an algorithm that mainly solves constraint problems, the BHGA has a great improvement in computing ability compared with GA, and it has better computing stability than the BSGA.

4.3. Benchmark of Constrained Engineering Functions

In this section, the BHGA was also tested with three constrained engineering problems [26]: a tension/compression spring (Eg01), a welded beam (Eg02), and a pressure vessel (Eg03).
The comparison of the BHGA optimization results with literature for the engineering problem is listed in Table 7, Table 8 and Table 9, where the bold number indicates the best results. Inspecting the results of the algorithms on those problems makes it evident that the BHGA managed to show very competitive results compared to IGA [22], BSGA [17], GA [38], GA [18], and TLBO [39], and obtains a better result on the pressure vessel problem. As an improved GA, the BHGA obtains better results than GA with fewer evaluations and also has a great improvement in computing ability compared with the BSGA. Taken together, the BHGA is efficient as a constrained handling method, especially for engineering constraint problems and constrained problems.

5. Lightweight Design Based on BHGA

In this chapter, the BHGA method is adopted to lighten the BIW mass. As the symmetry of the BIW structure, the design variables of beam members (1, 2, 3, 4, 5, 8, 9, 10, 48) are defined to have the same properties as beam members (17, 18, 19, 20, 21, 24, 25, 26, 49). Take the bending stiffness and the torsion stiffness as the constraints, and set the allowable limit values of constraint condition with respect to bending stiffness and torsion stiffness are 6000 N/mm and 2500 Nm/deg, respectively. The configurations of the BHGA are listed in Table 10.
The mass convergence curve for the optimization process is obtained, as depicted in Figure 10. Moreover, the optimized values of the cross-sections are listed in Table 11, where the initial and bounds values of design variables are also listed. According to the optimized and initial values of the cross-sections, it can be seen that the mass of the auto-body decreases by about 14.8 kg (from 111.6 kg to 96.8 kg). Consequently, the results indicate that it is effective to use the BHGA for the lightweight design of the BIW structure.

6. Conclusions

The TBTMM is used to establish the relationship between the section properties and structural mechanical performances to improve the calculation accuracy. The bending and torsional stiffness errors of the mathematical simulation model and the finite element model are 6.9% and 2.1%, respectively, which are within the reasonable error range. Moreover, as the optimal solutions of the constrained optimization problems are usually distributed along the constraint boundary, a more general and simpler constrained optimization algorithm BHGA is proposed based on a hybrid genetic algorithm and LS. Twenty problems (six unconstrained, eleven constrained, and three engineering-constrained problems) are benchmark tested and compared with well-regarded algorithms, which proves the effectiveness of this algorithm. Finally, BHGA is used in the lightweight design of vehicle BIW, providing initial design parameters for the vehicle conceptual design stage.
This paper only considers the static performances of BIW with single-cell thin-walled beams and describes its lightweight design as a single objective optimization problem. In further studies, we will focus on two aspects. On the one hand, we will improve the mechanical model of the BIW, such as establishing the static transfer matrix and dynamic transfer matrix of thin-walled beams with arbitrary cross-sections, improving the joint mechanical transfer model of the thin-walled beam, and establishing the transfer matrix of the curved thin-walled beam. On the other hand, we will combine other genetic operators and diversity-maintaining strategies to improve the BHGA and develop its parallel processing capability.

Author Contributions

Conceptualization, T.X. and M.S.; methodology, H.Z. and M.S.; software, H.Z.; validation, T.X.; M.S., H.Z. and J.Y.; formal analysis, H.Z. and J.Y.; investigation, H.Z.; resources, T.X.; data curation, H.Z.; writing—original draft preparation, H.Z.; writing—review and editing, T.X. and F.G.; visualization, M.S.; supervision, T.X.; project administration, T.X.; funding acquisition, T.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “Research and development of energy-saving and environment-friendly high-performance non-pneumatic tire” and “Research and development of a new type of non pneumatic tire applied to micro vehicles”. The grant numbers are X220091TL220 and X201011XQ200 respectively.

Data Availability Statement

Data are available upon request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Donders, S.; Takahashi, Y.; Hadjit, R.; Van Langenhove, T.; Brughmans, M.; Van Genechten, B.; Desmet, W. A reduced beam and joint concept modeling approach to optimize global vehicle body dynamics. Finite Elem. Anal. Des. 2009, 45, 439–455. [Google Scholar] [CrossRef]
  2. Bai, J.; Zhao, Y.; Meng, G.; Zuo, W. Bridging Topological Results and Thin-Walled Frame Structures Considering Manufacturability. J. Mech. Des. 2021, 143, 091706. [Google Scholar] [CrossRef]
  3. Mundo, D.; Hadjit, R.; Donders, S.; Brughmans, M.; Mas, P.; Desmet, W. Simplified modelling of joints and beam-like structures for BIW optimization in a concept phase of the vehicle design process. Finite Elem. Anal. Des. 2009, 45, 456–462. [Google Scholar] [CrossRef]
  4. Qin, H.; Liu, Z.; Liu, Y.; Zhong, H. An object-oriented MATLAB toolbox for automotive body conceptual design using distributed parallel optimization. Adv. Eng. Softw. 2017, 106, 19–32. [Google Scholar] [CrossRef]
  5. Liu, Y.; Liu, Z.; Qin, H.; Zhong, H.; Lv, C. An efficient structural optimization approach for the modular automotive body conceptual design. Struct. Multidiscip. Optim. 2018, 58, 1275–1289. [Google Scholar] [CrossRef]
  6. Nguyen, N.-L.; Jang, G.-W. Joint modeling using nonrigid cross-sections for beam-based analysis of a car body. Comput. Struct. 2021, 257, 106648. [Google Scholar] [CrossRef]
  7. Zhong, H.; Liu, Z.; Qin, H.; Liu, Y. Static analysis of thin-walled space frame structures with arbitrary closed cross-sections using transfer matrix method. Thin-Walled Struct. 2018, 123, 255–269. [Google Scholar] [CrossRef]
  8. Hou, W.B.; Zhang, H.Z.; Chi, R.F.; Hu, P. Development of an intelligent CAE system for auto-body concept design. Int. J. Automot. Technol. 2009, 10, 175–180. [Google Scholar] [CrossRef]
  9. Maghawry, A.; Hodhod, R.; Omar, Y.; Kholief, M. An approach for optimizing multi-objective problems using hybrid genetic algorithms. Soft Comput. 2020, 25, 389–405. [Google Scholar] [CrossRef]
  10. Ersavas, C.; Karatepe, E. Optimum allocation of FACTS devices under load uncertainty based on penalty functions with genetic algorithm. Electr. Eng. 2016, 99, 73–84. [Google Scholar] [CrossRef]
  11. Paszkowicz, W. Properties of a genetic algorithm equipped with a dynamic penalty function. Comput. Mater. Sci. 2009, 45, 77–83. [Google Scholar] [CrossRef]
  12. Lin, C.-H. A rough penalty genetic algorithm for constrained optimization. Inf. Sci. 2013, 241, 119–137. [Google Scholar] [CrossRef]
  13. Koziel, S.; Michalewicz, Z. Evolutionary Algorithms, Homomorphous Mappings, and Constrained Parameter Optimization. Evol. Comput. 1999, 7, 19–44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zhou, Y.; Li, Y.; He, J.; Kang, L. Multi-objective and MGG evolutionary algorithm for constrained optimization. In Proceedings of the Congress on Evolutionary Computing 2003 (CEC’2003), Canberra, Australia, 8–12 December 2003; pp. 1–5. [Google Scholar]
  15. Runarsson, T.P.; Yao, X. Search Biases in Constrained Evolutionary Optimization. IEEE Trans. Syst. Man Cybern. 2005, 35, 233–243. [Google Scholar] [CrossRef] [Green Version]
  16. Salcedo-Sanz, S. A survey of repair methods used as constraint handling techniques in evolutionary algorithms. Comput. Sci. Rev. 2009, 3, 175–192. [Google Scholar] [CrossRef]
  17. Li, X.; Du, G. Inequality constraint handling in genetic algorithms using a boundary simulation method. Comput. Oper. Res. 2012, 39, 521–540. [Google Scholar] [CrossRef]
  18. Coello, C.A.C.; Montes, E.M. Constraint-handling in genetic algorithms through the use of dominance-based tournament selection. Adv. Eng. Inform. 2002, 16, 193–203. [Google Scholar] [CrossRef]
  19. Isaacs, A.; Ray, T.; Smith, W. Blessings of maintaining infeasible solutions for constrained multi-objective optimization problems. In Proceedings of the 2008 IEEE Congress on Evolutionary Computation (IEEE World Congress on Computational Intelligence), Hong Kong, China, 1–6 June 2008; pp. 2780–2787. [Google Scholar]
  20. Genetic Algorithms Toolbox 2001. Available online: http://codem.group.shef.ac.uk/index.php/ga-toolbox (accessed on 2 March 2021).
  21. Liang, J.J.; Suganthan, P.N.; Deb, K. Novel composition test functions for numerical global optimization. In Proceedings of the IEEE Swarm Intelligence Symposium, Pasadena, CA, USA, 8–10 June 2005; pp. 68–75. [Google Scholar]
  22. Qin, H.; Guo, Y.; Liu, Z.; Liu, Y.; Zhong, H. Shape optimization of automotive body frame using an improved genetic algorithm optimizer. Adv. Eng. Softw. 2018, 121, 235–249. [Google Scholar] [CrossRef]
  23. Mirjalili, S.; Gandomi, A.H.; Mirjalili, S.Z.; Saremi, S.; Faris, H.; Mirjalili, S.M. Salp swarm algorithm: A bio-inspired optimizer for engineering design problems. Adv. Eng. Softw. 2017, 114, 163–191. [Google Scholar] [CrossRef]
  24. Saremi, S.; Mirjalili, S.; Lewis, A. Grasshopper Optimisation Algorithm: Theory and application. Adv. Eng. Softw. 2017, 105, 30–47. [Google Scholar] [CrossRef]
  25. Mirjalili, S.; Lewis, A. The whale optimization algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  26. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  27. Kennedy, J. Particle swarm optimization. In Encyclopedia of Machine Learning; Springer: New York, NY, USA, 2011; pp. 760–766. [Google Scholar]
  28. Rashedi, E.; Nezamabadi-Pour, H.; Saryazdi, S. GSA: A Gravitational Search Algorithm. Inf. Sci. 2009, 179, 2232–2248. [Google Scholar] [CrossRef]
  29. Mirjalili, S.; Mirjalili, S.M.; Hatamlou, A. Multi-verse optimizer: A nature-inspired algorithm for global optimization. Neural Comput. Appl. 2016, 27, 495–513. [Google Scholar] [CrossRef]
  30. Geem, Z.W.; Kim, J.H.; Loganathan, G.V. A new heuristic optimization algorithm: Harmony search. Simulation 2001, 76, 60–68. [Google Scholar] [CrossRef]
  31. Liang, J.; Runarsson, T.P.; Mezura-Montes, E.; Clerc, M.; Suganthan, P.N.; Coello, C.A.C.; Deb, K. Problem Definitions and Evaluation Criteria for the CEC 2006 Special Session on Constrained Real-Parameter Optimization; Technical report; Nanyang Technological University: Singapore, 2006. [Google Scholar]
  32. Liu, J.; Teo, K.L.; Wang, X.; Wu, C. An exact penalty function-based differential search algorithm for constrained global optimization. Soft Comput. 2015, 20, 1305–1313. [Google Scholar] [CrossRef]
  33. Patel, V.K.; Savsani, V.J. Heat transfer search (HTS): A novel optimization algorithm. Inf. Sci. 2015, 324, 217–246. [Google Scholar] [CrossRef]
  34. Mezura-Montes, E.; Coello, C.A.C. A Simple Multimembered Evolution Strategy to Solve Constrained Optimization Problems. IEEE Trans. Evol. Comput. 2005, 9, 1–17. [Google Scholar] [CrossRef]
  35. Zavala, A.E.M.; Aguirre, A.H.; Diharce, E.R.V. Constrained optimization via particle evolutionary swarm optimization algorithm (PESO). In Proceedings of the GECCO’05, Washington, DC, USA, 25–29 June 2005; pp. 209–216. [Google Scholar] [CrossRef] [Green Version]
  36. Mezura-Montes, E.; Miranda-Varela, M.E.; Gómez-Ramón, R.D.C. Differential evolution in constrained numerical optimization: An empirical study. Inf. Sci. 2010, 180, 4223–4262. [Google Scholar] [CrossRef]
  37. Karaboga, D.; Akay, B. A modified Artificial Bee Colony (ABC) algorithm for constrained optimization problems. Appl. Soft Comput. 2011, 11, 3021–3031. [Google Scholar] [CrossRef]
  38. Coello, C.A.C. Use of a self-adaptive penalty approach for engineering optimization problems. Comput. Ind. 2000, 41, 113–127. [Google Scholar] [CrossRef]
  39. Rao, R.V.; Savsani, V.J.; Vakharia, D.P. Teaching–learning-based optimization: A novel method for constrained mechanical design optimization problems. Comput. Aided Des. 2011, 43, 303–315. [Google Scholar] [CrossRef]
Figure 1. The conceptual BIW model.
Figure 1. The conceptual BIW model.
Machines 10 01171 g001
Figure 2. The details of the cross-section.
Figure 2. The details of the cross-section.
Machines 10 01171 g002
Figure 3. The mathematic simulation model of BIW conceptual structure.
Figure 3. The mathematic simulation model of BIW conceptual structure.
Machines 10 01171 g003
Figure 4. The static load cases and boundary conditions of BIW structure: (a) bending condition, (b) torsion condition.
Figure 4. The static load cases and boundary conditions of BIW structure: (a) bending condition, (b) torsion condition.
Machines 10 01171 g004
Figure 5. The flowchart for the calculation of boundary region.
Figure 5. The flowchart for the calculation of boundary region.
Machines 10 01171 g005
Figure 6. The flowchart for the generation of feasible points.
Figure 6. The flowchart for the generation of feasible points.
Machines 10 01171 g006
Figure 7. The distribution of point sets in the design domain.
Figure 7. The distribution of point sets in the design domain.
Machines 10 01171 g007
Figure 8. The difference between the binary search method and the reverse binary search method.
Figure 8. The difference between the binary search method and the reverse binary search method.
Machines 10 01171 g008
Figure 9. The process of the Local search.
Figure 9. The process of the Local search.
Machines 10 01171 g009
Figure 10. The convergence of object function.
Figure 10. The convergence of object function.
Machines 10 01171 g010
Table 1. Material property of tread rubber.
Table 1. Material property of tread rubber.
Stiffness TypeFEM AnalysisTBTMM AnalysisError
Bending stiffness (N/mm)6515.38465963.5323−6.9%
Torsion stiffness (N.m/°)2642.91292697.23382.1%
Table 2. Details of benchmark functions.
Table 2. Details of benchmark functions.
Prob.NType of FunctionLINI ρ (%) BO
CF110nonlinear00100-
CF210nonlinear00100-
CF310nonlinear00100-
CF410nonlinear00100-
CF510nonlinear00100-
CF610nonlinear00100-
G0113quadratic900.01116/9
G0220nonlinear0299.99711/2
G045quadratic0652.12302/6
G062cubic020.00662/2
G0710quadratic350.00036/8
G082nonlinear020.85600/2
G097polynomial040.51212/4
G108linear330.00106/6
G123quadratic014.77130/1
G1915nonlinear0533.47615/5
G242linear0279.65562/2
Eg013nonlinear130.75142/4
Eg024cubic252.66274/7
Eg034cubic3175.91502/4
Table 3. The configurations of optimization.
Table 3. The configurations of optimization.
Prob. F P P c 1 P c 2 P m 1 P m 2 c ε Δ P o p s i z e MaxGen
CF1400.950.850.20.0110.11301000
CF2400.950.850.20.0110.11301000
CF3400.950.850.20.0110.11301000
CF4400.950.850.20.0110.11301000
CF5400.950.850.20.0110.11301000
CF6400.950.850.20.0110.11301000
G01400.950.850.20.0110.00010.12001200
G02400.950.850.20.0110.00010.12001200
G04400.950.850.20.0110.00010.12001200
G06400.950.850.20.0110.00010.12001200
G07400.950.850.20.0110.00010.12001200
G08400.950.850.20.0110.00010.12001200
G09400.950.850.20.0110.00010.12001200
G10400.950.850.20.0110.00010.12001200
G12400.950.850.20.0110.00010.12001200
G19400.950.850.20.0110.00010.12001200
G24400.950.850.20.0110.00010.12001200
Eg01400.950.850.20.0110.00010.1100200
Eg02400.950.850.20.0110.00010.1100200
Eg03400.950.850.20.0110.00010.1100200
Table 4. Statistical results of composite benchmark functions.
Table 4. Statistical results of composite benchmark functions.
AlgorithmStatisticF1F2F3F4F5F6
BHGAAve.43.6796.27198.42417.1364.99535.21
Std.68.2084.2058.1378.6575.1891.25
Rank334621
BSGAAve.66.67123.75205.02572.80410.51819.30
Std.118.42150.4675.75140.7483.57162.49
Rank54610118
IGAAve.46.6789.10191.15344.5590.21544.26
Std.77.614.1061.3454.4757.4168.40
Rank422242
SSAAve.36.67303.35240.79334.9928.18608.64
Std.55.56373.2586.5130.4533.2818.42
Rank2117113
GOAAve.120.00261.49341.83516.50195.43846.43
Std.121.48121.63132.22166.93191.4113.79
Rank101010899
WOAAve.120.92173.93417.50610.51141.62673.17
Std.129.6491.30157.50138.71122.15196.52
Rank117111174
GWOAve.82.19146.49200.72430.5693.48860.55
Std.115.1793.0672.71125.21102.17123.24
Rank6657511
PSOAve.113.33124.15191.89346.67131.98751.69
Std.107.4296.5679.37102.83106.71194.32
Rank953365
GSAAve.3.33186.67157.14410.00195.43814.97
Std.18.2650.7455.12156.06191.41113.47
Rank1915106
MVOAve.86.68177.41299.11392.9386.15815.65
Std.81.72118.18154.94126.35112.3816.17
Rank789437
HSAve.93.2575.29273.79524.60193.47846.86
Std.44.60255.75101.48138.01128.0612.77
Rank8189810
Table 5. Results obtained by the BHGA algorithm for 11 benchmark functions over 30 independent runs with 240,000 function evaluations.
Table 5. Results obtained by the BHGA algorithm for 11 benchmark functions over 30 independent runs with 240,000 function evaluations.
Prob.Opti.BestWorstMeanSD
G01−1515.0000−15.0000−15.00001.9 × 10−6
G02−0.8036−0.8036−0.7926−0.80084.5 × 10−3
G04−30,665.5387−30,665.5387−30,665.5387−30,665.53874.5 × 10−6
G06−6961.8139−6961.8139−6961.8139−6961.81394.5 × 10−6
G0724.306224.307824.941924.48641.9 × 10−1
G08−0.095825−0.095825−0.095825−0.0958253.1 × 10−11
G09680.6301680.6301680.6626680.64088.1 × 10−3
G107049.24807114.83057795.38017354.55701.2 × 10−2
G12−1−1.0000−1.0000−1.00006.1 × 10−13
G1932.655632.891238.855035.30451.5
G24−5.5080−5.5080−5.5080−5.50803.7 × 10−8
Table 6. Comparative results of benchmark functions obtained by different algorithms.
Table 6. Comparative results of benchmark functions obtained by different algorithms.
Prob. BHGAd-DSBSGAHTSGAPSODEABCBBOTLBO
g01Best−15.0000−15−14.9999−1514.44−15−15−15−14.977−15
Worst−15.0000−13−14.9996−15 −13−11.828−15−14.5882−6
Mean−15.0000−12.3−14.9997−15−14.236−14.71−14.555−15−14.7698−10.782
SD1.9 × 10−61.8 × 10−26.7 × 10−5-------
g02Best−0.8036−0.803518−0.8036−0.7515−0.796321−0.669158−0.472−0.803598−0.7821−0.7835
Worst−0.7926−0.7743−0.7215−0.5482−0.299426-−0.749797−0.7389−0.5518
Mean−0.8008−0.7880−0.7669−0.6437−0.788588−0.41996−0.655−0.792412−0.7642−0.6705
SD4.5 × 10−37.0 × 10−42.3 × 10−2-------
g04Best−30,665.5387−30,665.539−30,665.5385−30,665.5387−30626.053−30,665.539−30,665.539−30,665.539−30,665.539−30,665.539
Worst−30,665.5387−30,665.6475−30,665.5380−30,665.5387−30,665.539−30,665.539−30,665.539−29942.3−30,665.539
Mean−30,665.5387−30,665.8862−30,665.5383−30,665.5387−30590.455−30,665.539−30,665.539−30,665.539−30411.865−30,665.539
SD4.5 × 10−61.2 × 10−11.3 × 10−4-------
G06Best−6961.8139−6961.8139−6961.6025−6961.814−6952.472−6961.814−6954.434−6961.814−6961.814−6961.814
Worst−6961.81393.6128E+07−6959.5077−6961.814−6961.814−6954.434−6961.805−5404.4941−6961.814
Mean−6961.81391.8436E+06−6961.1706−6961.814−6872.204−6961.814−6954.434−6961.813−6181.7461−6961.814
SD4.5 × 10−68.1 × 1064.5 × 10−1-------
g07Best24.307824.31524.325024.310431.09724.3724.30624.3325.664524.3103
Worst24.941925.533636.381025.0083-56.05524.3325.1937.691227.6106
Mean24.486424.715325.312624.494534.9832.40724.3124.47329.82924.837
SD1.9 × 10−13.1 × 10−22.2-------
g08Best−0.095825−0.095825−0.095825−0.095825−0.095825−0.095825−0.095825−0.095825−0.095825−0.095825
Worst−0.095825−0.09582−0.095825−0.095825-−0.095825−0.095825−0.095825−0.095817−0.095825
Mean−0.095825−0.0958−0.095825−0.095825−0.095799−0.095825−0.095825−0.095825−0.095824−0.095825
SD3.1 × 10−1107.2 × 10−11-------
g09Best680.6301680.630680.6321680.6301685.994680.63680.63680.634680.6301680.6301
Worst680.6626681.1324680.7393680.644 680.631680.631680.653721.0795680.6456
Mean680.6408680.7132680.6587680.6329692.064680.63680.63680.64692.7162680.6336
SD8.1 × 10−31.1 × 10−32.6 × 10−2-------
g10Best7114.83057056.767479.55477049.48369079.777049.4817049.5487053.9047679.06817250.9704
Worst7795.38017846.789810074.69067252.0546-7894.8129264.8867604.1329570.57147291.3779
Mean7354.55707350.34498945.58457119.701510003.2257205.57147.3347224.4078764.98647257.0927
SD1.2 × 10−22.0 × 1018.0 × 102------
g12Best−1.0000−1−1−1−1−1−1−1−1−1
Worst−1.0000−1−1−1 −0.994−1−1−1−1
Mean−1.0000−1−1−1−1−0.998875−1−1−1−1
SD6.1 × 10−1306.0 × 10−13-------
g19Best32.891232.655633.536432.7132-33.535832.685133.332539.147132.7916
Worst38.855046.165847.206233.2140-39.844332.907838.561471.310636.1935
Mean35.304532.804737.458532.7903-36.617232.768036.007851.876934.0792
SD1.52.83.3--2.046.28 × 10−21.831.12 × 1019.33 × 10−1
g24Best−5.5080−5.5080−5.5080−5.5080-−5.5080−5.5080−5.5080−5.5080−5.5080
Worst−5.5080−5.4661−5.5080−5.5080-−5.5080−5.5080−5.5080−5.4857−5.5080
Mean−5.5080−5.5080−5.5080−5.5080-−5.5080−5.5080−5.5080−5.4982−5.5080
SD3.7 × 10−83.4 × 10−63.5 × 10−6--9.36 × 10−169.36 × 10−169.36 × 10−166.75 × 10−39.36 × 10−16
Table 7. Comparison of the BHGA optimization results with literature for the tension/compression spring problem.
Table 7. Comparison of the BHGA optimization results with literature for the tension/compression spring problem.
AlgorithmDesign VariablesOptimum
Result
Max.
Eval.
x1x2x3
BHGA0.0517020.35703411.2704550.01266520,000
IGA0.0517600.35842111.1910340.01266750,000
BSGA0.0524990.37650510.2166920.01267720,000
GA (2000)---0.012822900,000
GA (2002)0.0519890.36396510.8905220.01297380,000
TLBO---0.01266510,000
Table 8. Comparison of the BHGA optimization results with literature for the welded beam problem.
Table 8. Comparison of the BHGA optimization results with literature for the welded beam problem.
AlgorithmDesign VariablesOptimum
Result
Max.
Eval.
x1x2x3x4
BHGA0.2057113.4708419.0367810.2057291.72489320,000
IGA0.2052183.4815379.0368230.2057311.72559750,000
BSGA0.1918423.8023799.0234410.2063321.74919320,000
GA (2000)----1.748309900,000
GA (2002)0.2059863.4713289.0202240.2064801.72822680,000
TLBO----1.72485210,000
Table 9. Comparison of the BHGA optimization results with literature for the pressure vessel problem.
Table 9. Comparison of the BHGA optimization results with literature for the pressure vessel problem.
AlgorithmDesign VariablesOptimum
Result
Max.
Eval.
x1x2x3x4
BHGA0.7899380.39053040.9293191.6831315905.963320,000
IGA0.8157520.40393242.248583174.8147125957.989850,000
BSGA0.80740.399041.8153180.17745939.185720,000
GA (2000)0.8125000.43450040.323900200.0000006288.7445900,000
GA (2002)0.8125000.43750042.097398176.6540506059.946380,000
TLBO----6059.71433510,000
Table 10. The configurations of BIW optimization.
Table 10. The configurations of BIW optimization.
Prob. F P P c 1 P c 2 P m 1 P m 2 c ε Δ P o p s i z e MaxGen
BIW200.950.850.20.0110.00010.540200
Table 11. The initial bounds and optimized values of the design variables.
Table 11. The initial bounds and optimized values of the design variables.
No.Design Variable
a (mm)b (mm)t (mm)
InitialLBUBOptimumInitialLBUBOptimumInitialLBUBOptimum
180709070.009150406040.32962131.2210
250406040.000260507050.02982131.1240
350406058.928360507050.30712131.5336
450406043.04011009011090.09022131.4500
550406040.0000140130150130.00872131.1200
650406059.054280709070.38452131.6700
750406050.497980709077.09802131.7085
850406040.001380709070.77702132.2068
950406050.835960507050.00002131.1000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhong, H.; Xu, T.; Yang, J.; Sun, M.; Gao, F. Optimization Design of Automotive Body Stiffness Using a Boundary Hybrid Genetic Algorithm. Machines 2022, 10, 1171. https://doi.org/10.3390/machines10121171

AMA Style

Zhong H, Xu T, Yang J, Sun M, Gao F. Optimization Design of Automotive Body Stiffness Using a Boundary Hybrid Genetic Algorithm. Machines. 2022; 10(12):1171. https://doi.org/10.3390/machines10121171

Chicago/Turabian Style

Zhong, Haolong, Ting Xu, Jianglin Yang, Meng Sun, and Fei Gao. 2022. "Optimization Design of Automotive Body Stiffness Using a Boundary Hybrid Genetic Algorithm" Machines 10, no. 12: 1171. https://doi.org/10.3390/machines10121171

APA Style

Zhong, H., Xu, T., Yang, J., Sun, M., & Gao, F. (2022). Optimization Design of Automotive Body Stiffness Using a Boundary Hybrid Genetic Algorithm. Machines, 10(12), 1171. https://doi.org/10.3390/machines10121171

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop