Computers and Structures 83 (2005) 2121–2136 www.elsevier.com/locate/compstruc
A comparative study of metamodeling methods for multiobjective crashworthiness optimization
H. Fang a,*, M. Rais-Rohani b, Z. Liu c, M.F. Horstemeyer a,d
a Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS 39762, USA b Department of Aerospace Engineering, Mississippi State University, Mississippi State, MS 39762, USA c Department of Civil Engineering, Mississippi State University, Mississippi State, MS 39762, USA d Department of Mechanical Engineering, Mississippi State University, Mississippi State, MS 39762, USA
Received 28 January 2004; accepted 28 February 2005 Available online 24 June 2005
Abstract
The response surface methodology (RSM), which typically uses quadratic polynomials, is predominantly used for metamodeling in crashworthiness optimization because of the high computational cost of vehicle crash simulations. Research shows, however, that RSM may not be suitable for modeling highly nonlinear responses that can often be found in impact related problems, especially when using limited quantity of response samples. The radial basis func- tions (RBF) have been shown to be promising for highly nonlinear problems, but no application to crashworthiness problems has been found in the literature. In this study, metamodels by RSM and RBF are used for multiobjective opti- mization of a vehicle body in frontal collision, with validations by finite element simulations using the full-scale vehicle model. The results show that RSM is able to produce good approximation models for energy absorption, and the model appropriateness can be well predicted by ANOVA. However, in the case of peak acceleration, RBF is found to generate better models than RSM based on the same number of response samples, with the multiquadric function identified to be the most stable RBF. Although RBF models are computationally more expensive, the optimization results of RBF models are found to be more accurate. 2005 Elsevier Ltd. All rights reserved.
Keywords: Metamodeling; Crashworthiness; Multiobjective optimization; Radial basis function; Response surface methodology
1. Introduction permanent injuries. In recent years, vehicle safety (crash- worthiness) has drawn more public attention and In 2002, 42,815 people were killed in vehicle crashes. research on vehicle crashworthiness has gained momen- That represents an increase of 1.5% over 2001 and the tum in both academe and the automotive industry. highest level since 1990 [1]. Also in 2002, 2,926,000 Vehicle safety can be measured by parameters such as people were injured from vehicle crashes, many with the contact forces exerted on the occupants and/or the resulting accelerations during a vehicle crash [2]. Both safety parameters (i.e., contact force and acceleration) * Corresponding author. Tel.: +1 662 325 6696; fax: +1 662 are closely related to the amount of energy absorbed 325 5433. by the vehicle before the impact wave reaches the occu- E-mail address: [email protected] (H. Fang). pants. With the aid of finite element (FE) analysis
0045-7949/$ - see front matter 2005 Elsevier Ltd. All rights reserved. doi:10.1016/j.compstruc.2005.02.025 2122 H. Fang et al. / Computers and Structures 83 (2005) 2121–2136 programs designed for dynamic contact problems, such low-order RS models. Although these techniques are as LS-DYNA and PAM-CRASH, it is possible to per- very efficient for single objective optimization problems, form a crash simulation and evaluate these parameters. they may not be appropriate for problems involving Furthermore, by coupling such simulation tools with multiple objectives. Yang et al. [15], in their survey of lit- nonlinear mathematical programming procedures, we erature on local RSM approaches, have concluded that can optimize a vehicle body to improve its crashworthi- in multiobjective optimization problems these ap- ness characteristics while considering other important proaches could be ineffective because the response re- factors such as manufacturing cost and vehicle perfor- gion of interest would never be reduced to a small mance. However, there are many challenges with crash- neighborhood that was good for all the objectives that worthiness optimization. Among them, we can point to typically conflict with each other. the implicit relationships between safety parameters and Jin et al. [16], compared RSM, Kriging method structural design attributes as well as the computational (KM), radial basis functions (RBF), and multivariate cost associated with repeated transient FE analyses. adaptive regression splines (MARS) using fourteen dif- Such concerns have led to the widespread use of meta- ferent problems, with one representing a complex engi- modeling approaches in crashworthiness optimization. neering application. They showed that RBF was the Despite their efficiency, metamodeling techniques best for both large-scale and small-scale problems based could still require a significant number of crash simula- on evaluations of the coefficient of multiple determina- tions, especially when the number of design variables is tion (R2), relative average absolute error (RAAE), and large. For this reason, reduced models (i.e., component relative maximum absolute error (RMAE). RBF was FE models or full-vehicle FE models with fewer degrees found to be the best for overall performance on accu- of freedom) have been developed and used in crash sim- racy, robustness, problem types, sample size, efficiency, ulations [3–6]. Although helpful in understanding the and simplicity. By contrast, they showed RSM to be mechanism of crash and improving vehicle designs, the the worst, in fact not suitable, for modeling highly non- reduced models have two major limitations. Firstly, it linear problems. However, they used only a linear RBF is difficult to accurately match the exact structural con- that, in general, is not a very accurate metamodel for straint and loading conditions in the reduced model with modeling highly nonlinear problems. Krishnamurthy those in a full-scale model during impact accompanied [17] compared augmented and compactly supported with large deformations. Secondly, a reduced model RBF with KM, local moving least square (MLS), and can typically be used for only one impact scenario and global least square (GLS) using one mathematical func- as such it would not be appropriate for crashworthiness tion and one FE based problem. The MLS and GLS optimization involving multiple impacts. However, the models in his study were basically quadratic polynomi- recent advances in computer technology have made it als; therefore, they represented local and global RS mod- possible to use full-scale FE models in high fidelity crash els, respectively. He showed that RBF, KM, and MLS simulations. Such models, developed by the National produced comparable and accurate results, and that Crash Analysis Center (NCAC) and other agencies in GLS performed poorly. However, the example problems the US, have been successfully used for crash simulation in that study did not require complex engineering anal- with the results compared with real vehicle tests [7–10]. ysis and thus relatively large sample sizes were generated Among the commonly used metamodeling tech- and used. niques in crashworthiness optimization, the response Despite its simplicity, RSM provides efficient yet surface methodology (RSM), particularly the use of sec- accurate solutions to many engineering problems [18] ond-order polynomials, has been the predominant meth- and analysis of variance (ANOVA) can be used to pre- od mainly because of its simplicity and efficiency. dict model appropriateness or fitness before the model However, the drawback of using second-order response is used in design optimization. RBF, on the other hand, surface (RS) models is that they may not be appropriate is more expensive than RSM, because it uses a series of for creating global models over the entire design space computationally expensive functions for a single model; for highly nonlinear problems. Although it is possible therefore, it is less efficient in performing function eval- to develop higher-order RS models, they may not be uations. This drawback becomes apparent when solving effective or appropriate for crashworthiness optimiza- multiobjective design optimization problems in which tion, partly due to the high computational cost in exten- millions sometimes even billions of solutions need to sive sampling of the design space. be found in order to develop the Pareto Frontier. An- Recent innovations to improve both the accuracy other disadvantage of using RBF is that model fit- and efficiency of RSM include the development and ness cannot be checked using ANOVA, because by application of the sequential local RSM [11,12], adaptive definition an RBF passes exactly through all the design RSM [13], and trust-region-based RSM [14]. All these points. approaches partition the feasible design space into mul- The focus of this study is to compare RSM and RBF tiple small regions that can be accurately represented by using limited samplings from both a nonlinear mathe- H. Fang et al. / Computers and Structures 83 (2005) 2121–2136 2123 matical function and full-scale crash simulations, and to Eq. (2) may be expressed in matrix form as assess the appropriateness of the resulting metamodels f ¼ Xb^; ð3Þ in both problems. The remaining portion of the paper is organized as follows. In Section 2 a brief overview where the vector of unknown coefficients b^ represents of RSM and RBF is provided followed by the metamod- the least-square estimator of the true coefficient vector eling of a nonlinear mathematical function in Section 3. and is solved using the method of least squares as In Section 4, we describe the crash simulation model, b^ ¼ðXTXÞ 1ðX Tf Þ. ð4Þ multiobjective optimization problem, and corresponding results. This is followed by the concluding remarks in Statistical analysis techniques such as ANOVA can Section 5. be used to check the fitness of an RS model and to iden- tify the main effects of design variables. However, main effect analysis is not the focus of this study and will not 2. Metamodeling methodologies be discussed here. The major statistical parameters used for evaluating model fitness are the F statistic, R2, ad- 2 2 The basic idea of metamodeling is to construct an justed R ðRadjÞ, and root mean square error (RMSE). approximate model using function values at some sam- Note that these parameters are not totally independent pling points, which are typically determined using exper- of each other and are calculated as imental design methods such as factorial design, Latin ðSST SSEÞ=p F ¼ ; ð5Þ hypercube, central composite design, or Taguchi ortho- SSE=ðn p 1Þ gonal array. Model fitness is subsequently checked using various statistical methods. In this section, we give a R2 ¼ 1 SSE=SST; ð6Þ brief overview of the two methods of interest in our n 1 research. R2 ¼ 1 ð1 R2Þ ; ð7Þ adj n p 1 2.1. Response surface methodology (RSM) sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi SSE RMSE ¼ ; ð8Þ In RSM, we typically use first- or second-order mod- n p 1 els in the form of linear or quadratic polynomial func- tions to develop an approximate model that provides where p is the number of non-constant terms in the RS an explicit relationship between design variables and model, SSE is the sum of square errors, and SST is the the response of interest. The unknown coefficients in total sum of squares. SSE and SST are calculated as the model are approximated using the method of least Xn 0 2 squares. SSE ¼ ðfi fi Þ ; ð9Þ Let f(x) be the true objective or response function i¼1 and f 0(x) its approximation obtained using the second- Xn 2 order polynomial in the form SST ¼ ðfi f Þ ; ð10Þ i¼1 Xm Xm Xm 1 Xm 0 2 where f is the measured function value at the ith design f ðxÞ¼b0 þ bixi þ biixi þ bijxixj; ð1Þ i 0 i¼1 i¼1 i¼1 j¼iþ1 point, fi is the function value calculated from the poly- nomial at the ith design point, and f is the mean value of where m is the total number of design variables, x is i f . the ith design variable, and the bs are the unknown i Generally speaking, the larger the values of R2 and coefficients. For n sampling of design variables R2 , and the smaller the value of RMSE, the better the x (k =1,2 ...n, i =1,2,...m) and the corresponding adj ki fit. In situations where the number of design variables function values f (k =1,2...n), Eq. (1) leads to n linear k is large, it is more appropriate to look at R2 , because equations expressed as adj R2 always increases as the number of terms in the model 2 Xm Xm Xm 1 Xm is increased while Radj actually decreases if unnecessary ^ ^ ^ 2 ^ terms are added to the model [19]. In addition to these f1 ¼ b0 þ bix1i þ biix1i þ bijx1ix1j; i¼1 i¼1 i¼1 j¼iþ1 statistics, the accuracy of the RS model can also be mea- Xm Xm Xm 1 Xm sured by checking its predictability of response using the ^ ^ ^ 2 ^ 2 f2 ¼ b0 þ bix2i þ biix2i þ bijx2ix2j; prediction error sum of squares (PRESS) and R for pre- i¼1 i¼1 i¼1 j¼iþ1 ð2Þ 2 2 diction (Rprediction). The PRESS statistic and Rprediction are calculated as Xm Xm Xm 1 Xm Xn f b^ b^ x b^ x2 b^ x x . 0 2 n ¼ 0 þ i ni þ ii ni þ ij ni nj PRESS ¼ ½fi fðiÞ ; ð11Þ i¼1 i¼1 i¼1 j¼iþ1 i¼1 2124 H. Fang et al. / Computers and Structures 83 (2005) 2121–2136
2 Rprediction ¼ 1 PRESS=SST; ð12Þ An RBF using the aforementioned highly nonlinear functions does not work well for linear responses [17]. 0 where fðiÞ is the predicted value at the ith design point To solve this problem, we can augment an RBF by using the model created by (n 1) design points that ex- including a polynomial function such that clude the ith point. Xn Xm 0 f ðxÞ¼ ki/ðkx xikÞ þ cjpjðxÞ; ð16Þ 2.2. Radial basis functions (RBF) i¼1 j¼1 where m is the total number of terms in the polynomial, An RBF uses a series of basic functions that are sym- and c (j =1,2,..., m) is the corresponding coefficient. A metric and centered at each sampling point, and it was j detailed discussion on the polynomial functions that originally developed for scattered multivariate data may be used can be found in Ref. [17]. interpolation [20]. Applications of RBF include ocean It can be seen that Eq. (16) is underdetermined be- depth measurement, altitude measurement, rainfall cause there are more parameters to be solved than the interpolation, surveying, mapping, geographics and number of equations created with available sampling geology, image warping, and medical imaging. points. Therefore, the orthogonality condition is further Let f(x) be the true objective or response function imposed on coefficients k as and f0(x) its approximation obtained from a classical Xn RBF with the general form kipjðxiÞ¼0; for j ¼ 1; 2 ...m. ð17Þ Xn i¼1 f 0ðxÞ¼ k /ðkx x kÞ; ð13Þ i i Combining Eqs. (16) and (17) in matrix form gives i¼1 AP k f where n is the number of sampling points, x is the vector ¼ ; ð18Þ PT 0 c 0 of design variables, xi is the vector of design variables at the ith sampling point, kx xik is the Euclidean dis- where Ai,j = /(kxi xjk)(i =1,2...n,j =1,2,...n), T tance, / is a basis function, and ki is the unknown Pi,j = pj(xi)(i =1,2,...n, j =1,2...m), k =[k1k2...kn] , weighting coefficient. Therefore, an RBF is actually a T 0 0 0 T c =[c1c2 cm] and f =[f (x1)f (x2) f (xn)] . Eq. (18) linear combination of n basis functions with weighted consists of (n + m) equations and its solution gives coef- coefficients. Some of the most commonly used basis ficients k and c for the RBF in the form of Eq. (16). functions include: It should be noted that an RBF passes through all the sampling points exactly. This means that function values 2 2 6 • Thin-plate spline: /(r)=r log(cr ), 0 < c 1; from the approximate function are equal to the true cr2 6 • Gaussian: /ðrÞ¼e p,0<ffiffiffiffiffiffiffiffiffiffiffiffiffiffic 1; function values at the sampling points. This can be seen • Multiquadric: /ðrÞ¼ r2 þ c2,0