基于耐撞性的多學(xué)科近似模型預(yù)測精度研究

吳光強(qiáng)　郭建波

摘 要：以整車100%正面碰撞有限元模型為基礎(chǔ)研究了三種近似模型的預(yù)測精度，分析并選取前部結(jié)構(gòu)中對汽車碰撞安全性影響較大的12個部件厚度為變量，利用最優(yōu)拉丁超立方試驗設(shè)計方法生成80個樣本數(shù)據(jù)并進(jìn)行計算，應(yīng)用多學(xué)科優(yōu)化中常用的二次多項式響應(yīng)面（Quadratic Polynomial Response Surface， QPRS）、Kriging以及徑向基函數(shù)（Radial Basis Function，RBF）三種近似方法分別對選取部件的總質(zhì)量、吸收總能量、B柱最大加速度和踏板侵入量建立近似模型。結(jié)果表明：RBF近似方法對部件總質(zhì)量、吸收總能量、B柱最大加速度預(yù)測精度高于其他兩種方法，Kriging近似方法對踏板侵入量預(yù)測模型具有較好的精度，QPRS近似方法適合于部件總質(zhì)量的近似建模。

關(guān)鍵詞：耐撞性；二次多項式響應(yīng)面模型；Kriging模型；RBF模型；預(yù)測精度

中圖分類號：U462.2 文獻(xiàn)標(biāo)志碼：A 文章編號：1005-2550（2014）02-0006-07

隨著社會經(jīng)濟(jì)的發(fā)展，汽車保有量的增加導(dǎo)致交通事故不斷增加，因此汽車安全性受到了前所未有的重視。在整車碰撞安全性設(shè)計過程中，有限元分析方法已經(jīng)得到了廣泛的應(yīng)用。雖然隨著計算機(jī)的硬件與軟件水平的提高，有限元分析方法已經(jīng)成為解決汽車碰撞安全性的重要手段，但是傳統(tǒng)的有限元優(yōu)化方法主要是建立在設(shè)計人員經(jīng)驗之上，難以找到整車碰撞性能的最優(yōu)解[1]。近些年來發(fā)展起來的多學(xué)科優(yōu)化設(shè)計可以通過試驗設(shè)計方法（Design of Experiment， DOE）減少試驗數(shù)量，通過在設(shè)計空間內(nèi)建立響應(yīng)面或者類似的近似模型對整車碰撞性能進(jìn)行優(yōu)化，找到整車碰撞安全性的最優(yōu)解。

在過去十年已有許多研究人員對多學(xué)科近似模型在整車正面碰撞安全性中的應(yīng)用進(jìn)行了研究。其中Shujuan Hou[2]利用二次多項式響應(yīng)面模型對整車正面碰撞進(jìn)行優(yōu)化；王國春[3]利用Kriging模型對汽車前部結(jié)構(gòu)的耐撞性進(jìn)行了優(yōu)化；陳國棟[4]利用徑向基函數(shù)對整車正面碰撞進(jìn)行了多學(xué)科的優(yōu)化；F. FANG[5]對二次多項式響應(yīng)面和RBF模型在整車正面碰撞中的預(yù)測精度進(jìn)行了比較；J. Forsberg[6]利用響應(yīng)面模型和Kriging模型對防撞梁結(jié)構(gòu)進(jìn)行了優(yōu)化。在已有的研究中，大部分研究人員只是將單一或者兩種近似模型應(yīng)用到整車的正面碰撞中進(jìn)行預(yù)測精度的研究，而對于常用的三種近似模型預(yù)測精度以及三種建模方法的適用范圍沒有進(jìn)行詳細(xì)的對比研究。

本文根據(jù)整車正面碰撞過程中零部件的總吸能和吸能密度大小選取了前部結(jié)構(gòu)的12個零件（8個變量），利用最優(yōu)拉丁超立方試驗設(shè)計方法進(jìn)行了80組數(shù)據(jù)采樣，選取了正面碰撞中使用較多的三種多學(xué)科近似模型：二次多項式響應(yīng)面、Kriging以及RBF模型分別對12個部件總質(zhì)量、吸收總能量、B柱最大加速度以及踏板侵入量值構(gòu)造近似模型，比較了三種近似模型的預(yù)測精度，分析得到三種近似模型在整車正面碰撞時的適用范圍。

1 三種近似模型的基本理論

1.1 多項式響應(yīng)面模型

多項式響應(yīng)面（Polynomial Response Surface，PRS）

最初被應(yīng)用于物理試驗的擬合，后來在結(jié)構(gòu)優(yōu)化領(lǐng)域得到應(yīng)用。利用回歸分析和方差分析決定設(shè)計變量和響應(yīng)間的關(guān)系。在多項式響應(yīng)面模型中常用的主要有可以提供設(shè)計變量和響應(yīng)之間的顯式關(guān)系的一次和二次多項式響應(yīng)面模型。本文中將采用二次多項式響應(yīng)面模型，其矩陣表達(dá)式為：

（1）

式中：f=[f1 f2…fn]T為響應(yīng)量組成的向量；

為所有基函數(shù)組成的矩陣；

β=[β1 β2…βm]T為所有基函數(shù)系數(shù)組成的向量，其中m為待定系數(shù)的個數(shù)，n為樣本點的數(shù)量。

利用最小二乘法可求出向量β ：

（2）

1.2 Kriging模型

Kriging是一種距離加權(quán)的插值近似方法，它以已知樣本信息的動態(tài)構(gòu)造為基礎(chǔ)并充分考慮到變量在空間上的相關(guān)特征，通過建立對象問題的近似函數(shù)關(guān)系來預(yù)測未知點。Kriging模型有兩部分組成：其矩陣表達(dá)式為：

（3）

式中：f=[f1 f2…fm]，β=[β1 β2…βm]T，m為基函數(shù)（如常數(shù)項、一階或二階多項式）的數(shù)目；Y是矩陣x的已知函數(shù)；β為基函數(shù)系數(shù)的矢量；Z為隨機(jī)分布的誤差，其均值為0，方差為σ2，協(xié)方差為：

（4）

式（3）中fβ表示設(shè)計空間中的全局模型，類似于響應(yīng)面中的多項式模型。式（3）中的Z用于建立偏離fβ的模型，因此整個模型對試驗設(shè)計產(chǎn)生的樣本點進(jìn)行插值。

1.3 徑向基函數(shù)模型

徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)，簡稱徑向基函數(shù)（Radial Basis Function，RBF），是以徑向函數(shù)通過基函數(shù)線性疊加構(gòu)造的模型。徑向基的矩陣形式為：

（5）

式中：f=[f1 f2…fn]T；Ai，j=φ（‖xi-xj‖）（i=1，2…n， j=1，2…n），‖x-xi‖為歐幾里德（Euclidean）距離，φ為基函數(shù)；λ=[λ1λ2…λn]T 為基函數(shù)系數(shù)，n為樣本點數(shù)量。

徑向基函數(shù)近似模型的特性很大程度上依賴于基函數(shù)，常用的基函數(shù)有薄板樣條函數(shù)（Thin-plate spline）、高斯函數(shù)（Gaussian）、多二次函數(shù)（Multiquadric）、逆多二次函數(shù)（Inverse multiquadric）。當(dāng)采用高斯或逆多二次函數(shù)時，預(yù)測模型具有局部估計的特點，而采用薄板樣條或者多二次函數(shù)時，模型具有全局估計的特點。

2 試驗設(shè)計方法

試驗設(shè)計中樣本點的選取與近似模型的精度具有相關(guān)性，合理的安排試驗才能以最少的試驗次數(shù)最真實的反應(yīng)設(shè)計對象隨設(shè)計變量的變化趨勢。endprint

本文所采用的試驗設(shè)計方法為最優(yōu)拉丁超立方設(shè)計（Optimal Latin hypercube design，OLHD），OLHD是在采樣空間中將每一個設(shè)計變量分成n個間隔（水平），首先在每個變量的n個等間隔中隨機(jī)選取一個不重復(fù)的間隔，按均勻分布隨機(jī)產(chǎn)生一個樣本點，循環(huán)n次，即可得到n個樣本點。與其他傳統(tǒng)的試驗設(shè)計方法相比，最優(yōu)拉丁超立方設(shè)計使所有的試驗點盡量均勻地分布在設(shè)計空間，尤其是對于小樣本設(shè)計具有非常好的的空間填充性和均勻性[7]。圖1為利用拉丁超立方設(shè)計（LHD）與OLHD對兩因素9水平問題的研究，觀察9個采樣點的分布，可以發(fā)現(xiàn)OLHD在設(shè)計空間內(nèi)分布更均勻，在同樣條件下減少試驗設(shè)計的樣本數(shù)量。

a）隨機(jī)拉丁超立方 b）最優(yōu)拉丁超立方

圖1 隨機(jī)拉丁超立方設(shè)計和最優(yōu)拉丁超立方設(shè)計

3 常用誤差分析指標(biāo)

近似模型的預(yù)測精度直接影響多學(xué)科最優(yōu)解的有效性。模型精度越高，基于近似模型的多學(xué)科優(yōu)化最優(yōu)解的可信度也越高。所以，建立近似模型后需要對近似模型的整體精度進(jìn)行評價，常用評價指標(biāo)有確定性系數(shù)（Coefficient of Determination，R2）、均方根誤差（Root Mean Squared Error，RMSE），其數(shù)學(xué)表達(dá)式如下：

（6）

或：

（7）

（8）

式中：yi為測試樣本點的實際響應(yīng)值；yi為近似模型的預(yù)測值；yi為實際響應(yīng)值的均值；ntest為測試樣本的數(shù)量。

4 實例應(yīng)用及結(jié)果分析

4.1 整車正面碰撞有限元模型

圖2 整車碰撞有限元模型

本文使用整車有限元模型進(jìn)行正面100%重疊剛性壁障碰撞仿真，整車質(zhì)量1.665 t，包含105萬個單元和93萬個節(jié)點。整車以56.5 km/h的速度撞擊剛性墻。在不影響分析結(jié)果的條件下為了減少整車碰撞有限元模型的計算時間將原來的計算時間由150 ms縮短為100 ms。整車有限元模型如圖2所示。

4.2 設(shè)計變量與設(shè)計響應(yīng)量選擇

汽車的正面碰撞必須滿足相關(guān)的國家法規(guī)，而且隨著消費者對汽車安全性的重視，更加嚴(yán)格的新車評價標(biāo)準(zhǔn)（New Car Assessment Program，NCAP）已經(jīng)成為消費者選擇汽車的一個重要參考指標(biāo)。根據(jù)NCAP的評價指標(biāo)，整車在正面碰撞過程中，應(yīng)該使前部結(jié)構(gòu)能充分變形，最大可能的吸收能量，降低B柱的加速度最大值。同時為了滿足汽車輕量化的要求，在設(shè)計過程中不能過多的增加吸能部件總質(zhì)量。

表1 變量的初始值以及上下限

圖3 設(shè)計變量

依據(jù)碰撞中變形與吸能最大的原則，在原車碰撞有限元模型中根據(jù)零部件吸能總能量和吸能密度大小兩種原則選取了車身前部的12個部件的8種板料厚度（對稱部件用一個厚度變量）作為正面碰撞的設(shè)計變量，所選擇的8個變量初始值及上下限如表1所示，變量對應(yīng)的零部件位置如圖3所示。根據(jù)乘員安全性和車身輕量化原則，文中以12個部件總質(zhì)量m、吸收總能量E、B柱最大加速度值aB以及踏板侵入量In為碰撞后的響應(yīng)值。

4.3 模型預(yù)測精度分析

由于整車碰撞有限元模型的計算成本較高，采用最優(yōu)拉丁超立方試驗設(shè)計方法對設(shè)計空間進(jìn)行了10n（小樣本）次采樣，其中n為設(shè)計變量個數(shù)，最后共計選取了80組訓(xùn)練樣本點。

每個訓(xùn)練樣本點應(yīng)用LS-DYNA軟件進(jìn)行計算，每次正面碰撞有限元模型計算仿真耗時約10小時左右，經(jīng)過800小時的計算得到全部訓(xùn)練樣本點的仿真結(jié)果。對80組訓(xùn)練樣本點的仿真結(jié)果提取整車正面碰撞部件總質(zhì)量m，吸收總能量E、B柱最大加速度aB以及踏板侵入量In，利用二次多項式響應(yīng)面、Kriging以及RBF模型對三種響應(yīng)量分別建立近似模型，得到12組近似模型。

由于Kriging是一種精確的插值技術(shù)，模型精確通過樣本點[3]，而根據(jù)RBF的定義，RBF也具有相同的特性[5]，所以這兩種模型在訓(xùn)練樣本點處的擬合誤差基本為0。故若對三種模型進(jìn)行預(yù)測精度的評估，需額外選取檢測點對模型的預(yù)測精度進(jìn)行評價。本文在80組訓(xùn)練樣本點之外隨機(jī)選取了15組數(shù)據(jù)進(jìn)行近似模型預(yù)測精度的評價。

通過15組檢測點得到12組近似模型的確定性系數(shù)R2以及均方根誤差RMSE，12組近似模型的預(yù)測值與有限元模型計算出來的實際值的分布如圖4所示。對于確定性系數(shù)R2的值越接近1，近似模型的預(yù)測精度越高，而對于均方根誤差RMSE的值越小，模型的精度越高。在表2中列出了通過15個檢測點得到的近似模型對四種響應(yīng)量的R2和RMSE值。

從圖4中的樣本點分布可以看出，對于三種類型的近似模型，部件總質(zhì)量m的預(yù)測精度最高，吸收總能量E與B柱最大加速度aB次之，踏板侵入量In精度最差，從表2的數(shù)據(jù)上也可以得出相同的結(jié)論。分析原因主要是因為前部結(jié)構(gòu)質(zhì)量m完全取決于12個部件的厚度，所以與變量之間主要成線性關(guān)系；而對于吸收總能量E不僅與零部件的厚度有關(guān)系，且與部件的形狀有一定的關(guān)系，所以其具有弱非線性的特征；而影響車身B柱最大加速度aB和踏板侵入量In的因素最多，呈現(xiàn)出了強(qiáng)非線性的特性，所以可以得到三種近似模型對于響應(yīng)量的預(yù)測精度會隨著非線性程度的增強(qiáng)而出現(xiàn)下降的趨勢，在圖4中的表現(xiàn)為隨著非線性程度的增加，檢測樣本點沿著直線兩側(cè)的分布更加分散。

從三種近似模型對于部件總質(zhì)量m的預(yù)測精度中可以看出，RBF與二次多項式響應(yīng)面的預(yù)測精度高于Kriging模型，從圖4中可以看到RBF和二次多項式響應(yīng)面模型樣本點基本分布在直線之上。分析原因主要是由于RBF模型對于樣本點具有較強(qiáng)的容錯功能，即使樣本中含有“噪聲”輸入，也不會影響網(wǎng)絡(luò)的整體性能，而二次多項式響應(yīng)面模型近似過程中并不是通過每個樣本點，也同樣具有過濾噪點的功能，但是Kriging模型需要通過每個樣本點，所以對于樣本中的噪點不具有過濾功能。endprint

對于零部件吸收總能量E，隨著非線性程度的增加，三種近似模型的預(yù)測精度都有所降低，但是Kriging模型的預(yù)測精度相對于其他兩種近似模型預(yù)測精度相對會增加。

B柱最大加速度aB和踏板侵入量In具有強(qiáng)非線性的特點，從預(yù)測精度的值可以看出，RBF預(yù)測精度最好，Kriging其次，而且Kriging模型的預(yù)測精度已經(jīng)更加接近RBF模型的精度，尤其在踏板侵入量的預(yù)測精度中已超越RBF的預(yù)測精度，但二次多項式響應(yīng)面的預(yù)測精度最差，從圖4的分布上也可以得出這樣的結(jié)論。

通過以上分析，可以得出以下相應(yīng)的結(jié)論，在多維度小樣本的條件下，文中選擇的三種近似模型中，RBF模型在線性、弱非線性和強(qiáng)非線性的情況下都具有較好的預(yù)測精度，這與R. Jin[8]所得到的結(jié)論相同；Kriging模型在弱非線性和強(qiáng)非線性的預(yù)測精度次于RBF模型，二次多項式響應(yīng)面模型在線性情況下具有較好的預(yù)測精度。

表3 三種近似模型適用范圍

考慮到構(gòu)建三種近似模型所需的時間：RBF最長，Kriging模型次之，二次多項式模型時間最短，對于B柱最大加速度aB近似建模二次多項式時間約為0.5 s，Kriging約為1.7 s，RBF約為35 s。所以在構(gòu)建部件總質(zhì)量m近似模型時可以考慮使用二次多項式模型，在建立吸收總能量E和B柱最大加速度aB近似模型時可以考慮使用RBF近似模型來提高預(yù)測的精度。三種模型的適用范圍如表3所示。

5 結(jié)論

本文結(jié)合最優(yōu)拉丁超立方試驗設(shè)計方法、三種近似模型構(gòu)建方法以及整車正面碰撞有限元模型，分析了三種近似模型在整車正面碰撞中的預(yù)測精度使用范圍，得出以下結(jié)論：

（1）隨著響應(yīng)量非線性程度的增加，三種近似模型的預(yù)測精度都會降低，QPRS近似方法的降低程度最大；

（2）在討論的三種近似方法中，RBF方法在三種響應(yīng)量的近似建模中預(yù)測精度都比較好，Kriging方法在吸收總能量E和B柱最大加速度aB近似中精度小于RBF方法但是高于QRPS方法，Kriging模型適用于踏板侵入量In建模，QRPS方法只有在建立部件總質(zhì)量m時具有較高的精度；

（3）綜合考慮三種方法的建模精度和建模時間兩個因素，得到在整車100%正面碰撞工況下RBF方法適合于建立B柱最大加速度aB和吸收總能量E模型，Kriging適用于建立踏板侵入量模型，QPRS適合于建立部件總質(zhì)量m模型。

參考文獻(xiàn)：

[1]張勇，李光耀，鐘志華.基于移動最小二乘響應(yīng)面法的整車輕量化設(shè)計優(yōu)化[J].機(jī)械工程學(xué)報，2008， 44（11）：192-196.

[2]Shujuan Hou， Duo Dong， Lili Ren， et al. Multivariable crashworthiness optimization of vehicle body by unreplicated saturated factorial design[J]. Structural and Multidisciplinary Optimization， 2012， 46（6）：891-905.

[3]王國春，成艾國，胡朝輝，等.基于Kriging模型的汽車前部結(jié)構(gòu)的耐撞性優(yōu)化[J].汽車工程，2011， 33（3）：208-212.

[4]陳國棟，韓 旭，劉桂萍，等.基于自適應(yīng)徑向基函數(shù)的整車耐撞性多目標(biāo)優(yōu)化[J].中國機(jī)械工程，2011，22（40）：488-493.

[5]H Fang， M Rais-Rohani， Z Liu， M F Horstemeyer. A comparative study of metamodeling methods for multiobjective crashworthiness optimization[J]. Computers & Structures， 2005， 83（25-26）：2121-2136.

[6]J.Forsberg， L.Nilsson. On polynomial response surfaces and Kriging for use in structural optimization of crashworthiness[J]. Structural and Multidisciplinary Optimization， 2005， 29（3）：232-243.

[7]Xingtao Liao， Qing Li， Xujing Yang， Weigang Zhang， Wei Li. Multiobjective optimization for crash safety design of vehicles using stepwise regression model [J]. Structural and Multidisciplinary Optimization， 2008， 35（6）：561-569.

[8]R Jin， W Chen， T W Simpson. Comparative studies of metamodelling techniques under multiple modeling criteria[J]. Structural and Multidisciplinary Optimization， 2001， 23（1）：1-13.endprint

對于零部件吸收總能量E，隨著非線性程度的增加，三種近似模型的預(yù)測精度都有所降低，但是Kriging模型的預(yù)測精度相對于其他兩種近似模型預(yù)測精度相對會增加。

B柱最大加速度aB和踏板侵入量In具有強(qiáng)非線性的特點，從預(yù)測精度的值可以看出，RBF預(yù)測精度最好，Kriging其次，而且Kriging模型的預(yù)測精度已經(jīng)更加接近RBF模型的精度，尤其在踏板侵入量的預(yù)測精度中已超越RBF的預(yù)測精度，但二次多項式響應(yīng)面的預(yù)測精度最差，從圖4的分布上也可以得出這樣的結(jié)論。

通過以上分析，可以得出以下相應(yīng)的結(jié)論，在多維度小樣本的條件下，文中選擇的三種近似模型中，RBF模型在線性、弱非線性和強(qiáng)非線性的情況下都具有較好的預(yù)測精度，這與R. Jin[8]所得到的結(jié)論相同；Kriging模型在弱非線性和強(qiáng)非線性的預(yù)測精度次于RBF模型，二次多項式響應(yīng)面模型在線性情況下具有較好的預(yù)測精度。

表3 三種近似模型適用范圍

考慮到構(gòu)建三種近似模型所需的時間：RBF最長，Kriging模型次之，二次多項式模型時間最短，對于B柱最大加速度aB近似建模二次多項式時間約為0.5 s，Kriging約為1.7 s，RBF約為35 s。所以在構(gòu)建部件總質(zhì)量m近似模型時可以考慮使用二次多項式模型，在建立吸收總能量E和B柱最大加速度aB近似模型時可以考慮使用RBF近似模型來提高預(yù)測的精度。三種模型的適用范圍如表3所示。

5 結(jié)論

本文結(jié)合最優(yōu)拉丁超立方試驗設(shè)計方法、三種近似模型構(gòu)建方法以及整車正面碰撞有限元模型，分析了三種近似模型在整車正面碰撞中的預(yù)測精度使用范圍，得出以下結(jié)論：

（1）隨著響應(yīng)量非線性程度的增加，三種近似模型的預(yù)測精度都會降低，QPRS近似方法的降低程度最大；

（2）在討論的三種近似方法中，RBF方法在三種響應(yīng)量的近似建模中預(yù)測精度都比較好，Kriging方法在吸收總能量E和B柱最大加速度aB近似中精度小于RBF方法但是高于QRPS方法，Kriging模型適用于踏板侵入量In建模，QRPS方法只有在建立部件總質(zhì)量m時具有較高的精度；

（3）綜合考慮三種方法的建模精度和建模時間兩個因素，得到在整車100%正面碰撞工況下RBF方法適合于建立B柱最大加速度aB和吸收總能量E模型，Kriging適用于建立踏板侵入量模型，QPRS適合于建立部件總質(zhì)量m模型。

參考文獻(xiàn)：

[1]張勇，李光耀，鐘志華.基于移動最小二乘響應(yīng)面法的整車輕量化設(shè)計優(yōu)化[J].機(jī)械工程學(xué)報，2008， 44（11）：192-196.

[2]Shujuan Hou， Duo Dong， Lili Ren， et al. Multivariable crashworthiness optimization of vehicle body by unreplicated saturated factorial design[J]. Structural and Multidisciplinary Optimization， 2012， 46（6）：891-905.

[3]王國春，成艾國，胡朝輝，等.基于Kriging模型的汽車前部結(jié)構(gòu)的耐撞性優(yōu)化[J].汽車工程，2011， 33（3）：208-212.

[4]陳國棟，韓 旭，劉桂萍，等.基于自適應(yīng)徑向基函數(shù)的整車耐撞性多目標(biāo)優(yōu)化[J].中國機(jī)械工程，2011，22（40）：488-493.

[5]H Fang， M Rais-Rohani， Z Liu， M F Horstemeyer. A comparative study of metamodeling methods for multiobjective crashworthiness optimization[J]. Computers & Structures， 2005， 83（25-26）：2121-2136.

[6]J.Forsberg， L.Nilsson. On polynomial response surfaces and Kriging for use in structural optimization of crashworthiness[J]. Structural and Multidisciplinary Optimization， 2005， 29（3）：232-243.

[7]Xingtao Liao， Qing Li， Xujing Yang， Weigang Zhang， Wei Li. Multiobjective optimization for crash safety design of vehicles using stepwise regression model [J]. Structural and Multidisciplinary Optimization， 2008， 35（6）：561-569.

[8]R Jin， W Chen， T W Simpson. Comparative studies of metamodelling techniques under multiple modeling criteria[J]. Structural and Multidisciplinary Optimization， 2001， 23（1）：1-13.endprint

對于零部件吸收總能量E，隨著非線性程度的增加，三種近似模型的預(yù)測精度都有所降低，但是Kriging模型的預(yù)測精度相對于其他兩種近似模型預(yù)測精度相對會增加。

B柱最大加速度aB和踏板侵入量In具有強(qiáng)非線性的特點，從預(yù)測精度的值可以看出，RBF預(yù)測精度最好，Kriging其次，而且Kriging模型的預(yù)測精度已經(jīng)更加接近RBF模型的精度，尤其在踏板侵入量的預(yù)測精度中已超越RBF的預(yù)測精度，但二次多項式響應(yīng)面的預(yù)測精度最差，從圖4的分布上也可以得出這樣的結(jié)論。

通過以上分析，可以得出以下相應(yīng)的結(jié)論，在多維度小樣本的條件下，文中選擇的三種近似模型中，RBF模型在線性、弱非線性和強(qiáng)非線性的情況下都具有較好的預(yù)測精度，這與R. Jin[8]所得到的結(jié)論相同；Kriging模型在弱非線性和強(qiáng)非線性的預(yù)測精度次于RBF模型，二次多項式響應(yīng)面模型在線性情況下具有較好的預(yù)測精度。

表3 三種近似模型適用范圍

考慮到構(gòu)建三種近似模型所需的時間：RBF最長，Kriging模型次之，二次多項式模型時間最短，對于B柱最大加速度aB近似建模二次多項式時間約為0.5 s，Kriging約為1.7 s，RBF約為35 s。所以在構(gòu)建部件總質(zhì)量m近似模型時可以考慮使用二次多項式模型，在建立吸收總能量E和B柱最大加速度aB近似模型時可以考慮使用RBF近似模型來提高預(yù)測的精度。三種模型的適用范圍如表3所示。

5 結(jié)論

本文結(jié)合最優(yōu)拉丁超立方試驗設(shè)計方法、三種近似模型構(gòu)建方法以及整車正面碰撞有限元模型，分析了三種近似模型在整車正面碰撞中的預(yù)測精度使用范圍，得出以下結(jié)論：

（1）隨著響應(yīng)量非線性程度的增加，三種近似模型的預(yù)測精度都會降低，QPRS近似方法的降低程度最大；

（2）在討論的三種近似方法中，RBF方法在三種響應(yīng)量的近似建模中預(yù)測精度都比較好，Kriging方法在吸收總能量E和B柱最大加速度aB近似中精度小于RBF方法但是高于QRPS方法，Kriging模型適用于踏板侵入量In建模，QRPS方法只有在建立部件總質(zhì)量m時具有較高的精度；

（3）綜合考慮三種方法的建模精度和建模時間兩個因素，得到在整車100%正面碰撞工況下RBF方法適合于建立B柱最大加速度aB和吸收總能量E模型，Kriging適用于建立踏板侵入量模型，QPRS適合于建立部件總質(zhì)量m模型。

參考文獻(xiàn)：

[1]張勇，李光耀，鐘志華.基于移動最小二乘響應(yīng)面法的整車輕量化設(shè)計優(yōu)化[J].機(jī)械工程學(xué)報，2008， 44（11）：192-196.

[2]Shujuan Hou， Duo Dong， Lili Ren， et al. Multivariable crashworthiness optimization of vehicle body by unreplicated saturated factorial design[J]. Structural and Multidisciplinary Optimization， 2012， 46（6）：891-905.

[3]王國春，成艾國，胡朝輝，等.基于Kriging模型的汽車前部結(jié)構(gòu)的耐撞性優(yōu)化[J].汽車工程，2011， 33（3）：208-212.

[4]陳國棟，韓 旭，劉桂萍，等.基于自適應(yīng)徑向基函數(shù)的整車耐撞性多目標(biāo)優(yōu)化[J].中國機(jī)械工程，2011，22（40）：488-493.

[5]H Fang， M Rais-Rohani， Z Liu， M F Horstemeyer. A comparative study of metamodeling methods for multiobjective crashworthiness optimization[J]. Computers & Structures， 2005， 83（25-26）：2121-2136.

[6]J.Forsberg， L.Nilsson. On polynomial response surfaces and Kriging for use in structural optimization of crashworthiness[J]. Structural and Multidisciplinary Optimization， 2005， 29（3）：232-243.

[7]Xingtao Liao， Qing Li， Xujing Yang， Weigang Zhang， Wei Li. Multiobjective optimization for crash safety design of vehicles using stepwise regression model [J]. Structural and Multidisciplinary Optimization， 2008， 35（6）：561-569.

[8]R Jin， W Chen， T W Simpson. Comparative studies of metamodelling techniques under multiple modeling criteria[J]. Structural and Multidisciplinary Optimization， 2001， 23（1）：1-13.endprint