基于解耦型GMDH的車身材料參數(shù)反求

殷為洋　王　琥　馮　慧　湯　龍

湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙,410082

基于解耦型GMDH的車身材料參數(shù)反求

殷為洋王琥馮慧湯龍

湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙,410082

動(dòng)態(tài)載荷作用時(shí),具有應(yīng)變率效應(yīng)的材料在碰撞仿真中會(huì)展現(xiàn)出不同于靜態(tài)載荷時(shí)的性能,準(zhǔn)確的材料參數(shù)是獲得可靠仿真結(jié)果的前提。主流的近似模型優(yōu)化方法忽略了對(duì)變量間耦合關(guān)系的判定,造成近似模型中存在不必要的耦合項(xiàng),增大誤差項(xiàng)所占比重,降低模型的效率和泛化能力。為此,提出了解耦型數(shù)據(jù)分組處理方法(GMDH),在建模初期判斷變量之間的耦合關(guān)系,進(jìn)而確定模型的耦合項(xiàng)。在高維非線性函數(shù)測(cè)試中,該方法表現(xiàn)出優(yōu)良的建模性能;將該方法與臺(tái)車試驗(yàn)結(jié)合,反求出兩種材料構(gòu)成的拼焊板參數(shù),與試驗(yàn)結(jié)果的對(duì)比表明,該反求方法具有較高的精度。

應(yīng)變率效應(yīng);耦合關(guān)系;解耦型數(shù)據(jù)分組處理方法(GMDH);近似模型

0　引言

汽車輕量化是涉及到安全、節(jié)能和環(huán)保的綜合性課題。車體減重10%可節(jié)省燃料8%～10%[1]。隨著科學(xué)技術(shù)的進(jìn)步,新型材料不斷涌現(xiàn),這為汽車輕量化開(kāi)辟了廣闊的發(fā)展空間。大部分車身材料具有應(yīng)變率效應(yīng),即在動(dòng)載荷作用時(shí)會(huì)呈現(xiàn)出不同于靜載荷作用時(shí)的力學(xué)特性。汽車碰撞是以大位移、大轉(zhuǎn)動(dòng)及大應(yīng)變?yōu)樘卣鞯膹?fù)雜非線性動(dòng)力學(xué)問(wèn)題。在傳統(tǒng)的材料參數(shù)確定中,材料參數(shù)是在標(biāo)準(zhǔn)試樣件拉伸試驗(yàn)的基礎(chǔ)上通過(guò)測(cè)試解析關(guān)系得到的,不能反映材料在汽車碰撞仿真中的動(dòng)態(tài)性能，在仿真時(shí)會(huì)造成相當(dāng)大的誤差。因此,選擇合適的應(yīng)變率模型并獲得正確的材料參數(shù)是汽車碰撞仿真的前提[2]。

目前,國(guó)際上主流的參數(shù)反求方法是試驗(yàn)與仿真結(jié)合的混合數(shù)值方法，即將材料參數(shù)反求轉(zhuǎn)化為工程優(yōu)化問(wèn)題：以試驗(yàn)測(cè)試結(jié)果為目標(biāo)依據(jù)，利用計(jì)算機(jī)對(duì)真實(shí)試驗(yàn)進(jìn)行仿真，將仿真值與試驗(yàn)值之間的誤差函數(shù)最小化，從而得到材料參數(shù)。曹銀峰等[3]介紹了參數(shù)反求的基本公式和關(guān)鍵算法,并將試驗(yàn)點(diǎn)的活度規(guī)則與試驗(yàn)結(jié)合反求出材料參數(shù);Eggertsen等[4]將周期屈服試驗(yàn)與響應(yīng)面結(jié)合，反求出模擬回彈的材料參數(shù);Sedighi等[5]將Levenberg-Marquardt算法與高速拉伸試驗(yàn)結(jié)合，反求了兩種鋼的材料參數(shù);李恩穎等[6]將最小二乘支持向量機(jī)回歸(least squares support vector regression,LS-SVR)與碰撞試驗(yàn)結(jié)合，反求出了DP600高強(qiáng)鋼的材料參數(shù)。

隨復(fù)雜度的增大,花費(fèi)在仿真計(jì)算的時(shí)間成本大幅度增加，為了提高效率和節(jié)省成本，有效方法是建立近似模型。目前,國(guó)內(nèi)外構(gòu)造近似模型的方法有神經(jīng)網(wǎng)絡(luò)、多項(xiàng)式響應(yīng)面及移動(dòng)最小二乘等[7]。通過(guò)驗(yàn)證,發(fā)現(xiàn)這些近似模型方法應(yīng)用于低維問(wèn)題時(shí)有較好的效果,然而,隨著問(wèn)題維數(shù)的增加,用于構(gòu)造這些近似模型的樣本點(diǎn)數(shù)目及計(jì)算成本會(huì)急劇增加,同時(shí)變量之間不必要的耦合項(xiàng)會(huì)變成冗雜項(xiàng),增加網(wǎng)絡(luò)的冗雜度,大大降低模型的效率和泛化能力。

結(jié)合樣本點(diǎn)成本問(wèn)題與高維變量間耦合性問(wèn)題,本文提出解耦型數(shù)據(jù)分組處理方法(group method of data handling, GMDH)近似模型：在傳統(tǒng)GMDH的基礎(chǔ)上，引入耦合性判斷機(jī)制，根據(jù)判斷結(jié)果決定是否保留耦合項(xiàng)。在高維測(cè)試函數(shù)算例中，算法表現(xiàn)出較高的精度。最后將該算法成功運(yùn)用到車身拼焊板材料參數(shù)的反求中。

1　GMDH近似模型

GMDH的基本思想是用模擬生物有機(jī)體演化的方式來(lái)構(gòu)造數(shù)學(xué)模型,它不要求模型開(kāi)發(fā)者完全掌握系統(tǒng)的內(nèi)部變化機(jī)理,不需要事先確定模型及其參數(shù),而是利用輸入、輸出的實(shí)測(cè)數(shù)據(jù)確定模型輸入變量及其參數(shù)。在建模過(guò)程中并非利用人為事先確定的函數(shù)關(guān)系式，而是利用所提供的建模數(shù)據(jù)通過(guò)計(jì)算機(jī)自動(dòng)尋找數(shù)據(jù)間的函數(shù)關(guān)系。

GMDH近似模型的構(gòu)造過(guò)程[8]如下:

(1)確定輸入變量xi(i=1,2,…,N)和目標(biāo)變量y,輸入變量即構(gòu)成第一層的各個(gè)活動(dòng)節(jié)點(diǎn)。

(2)將輸入變量?jī)蓛膳鋵?duì),組成N(N-1)/2個(gè)活動(dòng)節(jié)點(diǎn),以完全二次多項(xiàng)式作為活動(dòng)節(jié)點(diǎn)的部分描述,并對(duì)該層所有的節(jié)點(diǎn)進(jìn)行精度評(píng)估,保留一定數(shù)目的節(jié)點(diǎn)作為母本進(jìn)入下一層節(jié)點(diǎn)的構(gòu)造中,直至滿足終止準(zhǔn)則。

(3)終止準(zhǔn)則。滿足下列兩個(gè)條件中的任何一個(gè)即終止GMDH網(wǎng)絡(luò)的構(gòu)造:①當(dāng)只有一個(gè)節(jié)點(diǎn)時(shí)；②當(dāng)本層中的最小準(zhǔn)則值大于上一層的最小準(zhǔn)則值時(shí)(此時(shí)選擇上層中最小準(zhǔn)則值對(duì)應(yīng)的節(jié)點(diǎn))。

傳統(tǒng)的GMDH模型沒(méi)有判斷變量間的耦合關(guān)系而直接加入耦合項(xiàng),對(duì)于沒(méi)有耦合關(guān)系的變量會(huì)增加不必要的耦合項(xiàng)即誤差項(xiàng)。隨著構(gòu)造進(jìn)程的推進(jìn)，不必要的耦合項(xiàng)經(jīng)過(guò)多次再耦合，模型中的不必要項(xiàng)所占比重也逐漸增大,導(dǎo)致誤差項(xiàng)增大，大大降低近似模型的效率和精度。

2　解耦型GMDH近似模型

傳統(tǒng)GMDH模型在構(gòu)建初期不判斷耦合關(guān)系就引入耦合項(xiàng),最終影響近似模型精度。解耦型GMDH模型則在構(gòu)造第二層各個(gè)活動(dòng)節(jié)點(diǎn)之前，首先利用高維模型表示(high dimensional model representation,HDMR)方法[9]判斷變量間的耦合關(guān)系，根據(jù)判斷結(jié)果決定是否加入耦合項(xiàng)，其構(gòu)造過(guò)程如下：

(1)所有變量xi構(gòu)成網(wǎng)絡(luò)的第一層節(jié)點(diǎn)。

(2)取問(wèn)題域中心點(diǎn)并計(jì)算樣本值f0,該樣本點(diǎn)作為構(gòu)造單變量非耦合項(xiàng)的評(píng)價(jià)準(zhǔn)則。

(3)在單變量xi∈X(X=(y1,y2,…,xi,…,yN))的定義區(qū)間為構(gòu)造其非耦合項(xiàng)fi(xi)=f(XT)-f0布點(diǎn),然后利用徑向基函數(shù)(radialbasisfunction,RBF)方法[10]構(gòu)造fi(xi)。與其他插值方法相比,RBF插值方法特別適合多維離散點(diǎn)的插值，不會(huì)因?yàn)椴逯迭c(diǎn)的位置關(guān)系而產(chǎn)生病態(tài)的系數(shù)矩陣，而且使用比較少的網(wǎng)格點(diǎn)即可得到同最小二乘方法相當(dāng)?shù)木取?/p>

(4)判斷fi(xi)的線性與非線性。如果滿足收斂條件:

(1)

則單變量的非耦合項(xiàng)構(gòu)造完成;否則,繼續(xù)采樣構(gòu)造單變量的非耦合項(xiàng)直至滿足收斂條件。由于實(shí)際工程問(wèn)題中不確定因素的增多以及系統(tǒng)誤差的存在,在運(yùn)用解耦型GMDH算法解決工程問(wèn)題時(shí),可將收斂條件適當(dāng)放寬。

(5)循環(huán)執(zhí)行步驟(3)和(4)直至所有變量的非耦合項(xiàng)構(gòu)造完畢。

(6)判斷兩兩變量間是否存在耦合關(guān)系。識(shí)別變量之間的耦合性,構(gòu)造新樣本點(diǎn)WT=(y1,y2,…xi,…,xj,…yN)T(i≠j)并計(jì)算其真實(shí)樣本值f(WT)。在精度允許的誤差范圍內(nèi)(即設(shè)定一較小的閾值，如0.0001),若滿足關(guān)系f(WT)=f0+fi(xi)+fj(xj),即常數(shù)項(xiàng)和非耦合項(xiàng)構(gòu)成的近似模型已能夠模擬真實(shí)值,則認(rèn)為輸入變量xi與xj之間不存在耦合關(guān)系,否則兩變量之間存在耦合關(guān)系。

(7)構(gòu)建GMDH模型第二層的活動(dòng)節(jié)點(diǎn)。以兩變量的完全二次多項(xiàng)式為基本表達(dá)式,根據(jù)步驟(6)中的判斷結(jié)果,靈活決定活動(dòng)節(jié)點(diǎn)的部分表達(dá)式是否剔除耦合項(xiàng)a5x1x2。

(8)篩選并確定母本。利用AIC準(zhǔn)則(Akaike’s information criterion)[11]對(duì)步驟(7)中構(gòu)造出的所有節(jié)點(diǎn)進(jìn)行評(píng)價(jià),在本層所生成的所有活動(dòng)節(jié)點(diǎn)中,選擇AIC值最小的5個(gè)節(jié)點(diǎn)作為母本繼續(xù)構(gòu)造下一層，直至達(dá)到收斂條件。

AIC準(zhǔn)則是由Akaike在研究信息論特別是在解決時(shí)間序列定階問(wèn)題時(shí)提出來(lái)的[12],其顯著特點(diǎn)之一是“吝嗇原理(principle of parsimony)”的具體化,公式為

PAIC=-2lnP+2Q

(2)

式中,PAIC為AIC值;P為模型的極大似然函數(shù);Q為模型的獨(dú)立參數(shù)個(gè)數(shù)。

當(dāng)從待選擇模型中選擇最佳模型時(shí),AIC值最小的模型是最佳的。當(dāng)兩個(gè)模型之間存在相當(dāng)大的差異時(shí)，這個(gè)差異在式(2)右邊第一項(xiàng)得到表現(xiàn);而當(dāng)兩個(gè)模型之間的差異幾乎沒(méi)有時(shí),則式(2)第二項(xiàng)起作用，從而參數(shù)個(gè)數(shù)小的模型是好的模型。

解耦型GMDH近似模型構(gòu)建流程見(jiàn)圖1。

圖1　解耦型GMDH近似模型構(gòu)建流程圖

3　高維測(cè)試函數(shù)測(cè)試

為驗(yàn)證解耦型GMDH近似模型的精度,采用高維函數(shù)進(jìn)行測(cè)試。

3.1近似模型評(píng)價(jià)指標(biāo)[13]

(1)判定系數(shù)R2為

(3)

R2是一個(gè)整體指標(biāo),其值越接近1,模型越精確。

(2)平均絕對(duì)誤差ERAAE(relative average absolute error)為

(4)

式中,STD為標(biāo)準(zhǔn)差。

ERAAE是一個(gè)整體指標(biāo),其值越小,模型越精確。

(3)最大絕對(duì)誤差ERMAE(relative maximum absolute error)為

(5)

ERMAE是一個(gè)局部指標(biāo),它描述了設(shè)計(jì)空間中的局部誤差,其值越小模型精度越高。

3.2所用測(cè)試函數(shù)及結(jié)果

對(duì)每個(gè)測(cè)試函數(shù)建立其10維問(wèn)題的近似模型,均用隨機(jī)生成的樣本連續(xù)測(cè)試30次,其中De-GMDH代表解耦型GMDH方法,GMDH代表傳統(tǒng)數(shù)據(jù)組處理方法,PR代表多項(xiàng)式方法,RBF代表徑向基函數(shù)方法。

(1)不含耦合項(xiàng)的測(cè)試函數(shù)如下：

(6)

(7)

(8)

上述3個(gè)函數(shù)的各指標(biāo)結(jié)果對(duì)比分別見(jiàn)表1～表3。

表1　式(6)函數(shù)各指標(biāo)結(jié)果對(duì)比

表2　式(7)函數(shù)各指標(biāo)結(jié)果對(duì)比

表3　式(8)函數(shù)各指標(biāo)結(jié)果對(duì)比

以上各個(gè)函數(shù)中變量之間不存在耦合關(guān)系,采用解耦型GMDH構(gòu)建近似模型過(guò)程中均能準(zhǔn)確判斷變量的耦合性，在初始建模階段就剔除了耦合項(xiàng)，不僅簡(jiǎn)化了模型結(jié)構(gòu),還減小了誤差項(xiàng)的比重，提高了精度,即使樣本點(diǎn)數(shù)目減少一半，各個(gè)指標(biāo)值也均優(yōu)于傳統(tǒng)的GMDH。

(2)含有耦合項(xiàng)的測(cè)試函數(shù)如下：

(9)

(10)

(11)

上述3個(gè)函數(shù)的各指標(biāo)結(jié)果對(duì)比分別見(jiàn)表4～表6。

表4　式(9)函數(shù)各指標(biāo)結(jié)果對(duì)比

表5　式(10)函數(shù)各指標(biāo)結(jié)果對(duì)比

表6　式(11)函數(shù)各指標(biāo)結(jié)果對(duì)比

以上各個(gè)函數(shù)中變量之間存在耦合關(guān)系,解耦型GMDH保留了部分表達(dá)式中的耦合項(xiàng)，整個(gè)建模過(guò)程與傳統(tǒng)GMDH幾乎一樣,即樣本數(shù)目與精度基本保持一致。

根據(jù)上述兩種類型的高維測(cè)試函數(shù)測(cè)試結(jié)果可以得到：在構(gòu)建上述近似模型過(guò)程中，解耦型GMDH能夠準(zhǔn)確判定出各個(gè)變量之間的耦合關(guān)系，進(jìn)而在構(gòu)建第二層活動(dòng)節(jié)點(diǎn)時(shí)只需決定是否添加耦合項(xiàng)，從而達(dá)到了減小誤差項(xiàng)比重、提高擬合精度的效果。

4　建立薄壁吸能筒近似模型

4.1吸能筒的幾何參數(shù)

圖2　四部分構(gòu)成的薄壁吸能梁撞擊剛性墻

圖3　橫截面尺寸

根據(jù)實(shí)際的吸能筒臺(tái)車碰撞試驗(yàn)中所用吸能筒的幾何尺寸及材料匹配情況,建立雙吸能筒有限元模型，如圖2所示。其中，前后兩個(gè)吸能筒的材料分別為DP590和DP780,厚度分別為1.0 mm和1.5 mm,每個(gè)吸能筒由上下兩部分焊接組合而成,具體的橫截面幾何參數(shù)如圖3所示。圖3中各部分長(zhǎng)度如表7所示。此外每個(gè)吸能筒的長(zhǎng)度為200 mm。

表7　吸能筒橫截面尺寸　mm

4.2吸能筒臺(tái)車碰撞有限元模型

參照試驗(yàn)條件下的臺(tái)車結(jié)構(gòu)，在Hypermesh中建立八邊形主機(jī)吸能前縱梁-臺(tái)車結(jié)構(gòu)的有限元模型(圖4)。吸能筒的單元尺寸設(shè)置為5 mm, 臺(tái)車的速度設(shè)置為30 km/h。臺(tái)車部分的屬性設(shè)置為不變形的剛體,采用第20號(hào)剛性殼體,調(diào)整模型質(zhì)心位置,使之與實(shí)際臺(tái)車質(zhì)心位置一致。

圖4　吸能筒臺(tái)車碰撞的有限元仿真模型

在臺(tái)車碰撞仿真中,輪胎是臺(tái)車模型中比較重要的部件。為了提高模型的精度,在臺(tái)車模型中考慮車輪模型高速旋轉(zhuǎn)運(yùn)動(dòng)和摩擦特性，使車輪由高速旋轉(zhuǎn)逐漸停止。為實(shí)現(xiàn)仿真與試驗(yàn)中輪胎剛度的擬合,綜合考慮了輪胎單元屬性、材料參數(shù)、輪胎氣壓等，本文中的輪胎厚度與材料參數(shù)均為實(shí)際輪胎數(shù)據(jù)，輪胎剛度擬合參照文獻(xiàn)[14]。

4.3吸能筒的本構(gòu)模型

常用的率相關(guān)本構(gòu)模型主要有Johnson-Cook(JC)模型、Zerilli-Armstrong(ZA)模型、Cowper-Symonds(CS)模型及Piecewise-Linear(PL)模型等。其中JC模型引入了材料的應(yīng)變率強(qiáng)化及熱軟化參數(shù)，適用于大多數(shù)金屬材料的準(zhǔn)靜態(tài)變形到高速變形的仿真，如彈道穿透、碰撞等。本文選擇JC模型：

(12)

(13)

在不考慮溫度效應(yīng)的情況下忽略式(13)，JC模型簡(jiǎn)化為

(14)

4.4吸能筒的材料參數(shù)反求

在該焊接雙吸能筒的材料反求中,以兩種材料的JC模型參數(shù)為自變量(DP590的A、B、n、C和DP780的A、B、n、C),均勻選取21個(gè)加速度仿真值與對(duì)應(yīng)試驗(yàn)值的平均相對(duì)誤差作為目標(biāo)變量：

(15)

式中,ati為第i個(gè)試驗(yàn)數(shù)據(jù)加速度值;asi為第i個(gè)仿真數(shù)據(jù)加速度值；M為加速度值取樣點(diǎn)數(shù)。

各個(gè)加速度值是加速度曲線經(jīng)濾波處理后均勻截取的。目標(biāo)函數(shù)與各個(gè)待匹配參數(shù)的關(guān)系為

F=min(S(A1i,B1i,n1i,C1i,A2i,B2i,

n2i,C2i)-atest)=S(Xbest)-atest

(16)

(A1i,B1i,n1i,C1i,A2i,B2i,n2i,C2i)∈Di=1,2,…,M

式中,A1、B1、n1、C1、A2、B2、n2、C2為待求參數(shù);S(A1,B1,n1,C1,A2,B2,n2,C2)為該組參數(shù)對(duì)應(yīng)的仿真計(jì)算值;atest為試驗(yàn)值;Xbest為參數(shù)的最佳組合值；D為參數(shù)區(qū)間。

以R2為收斂條件,逐漸增加樣本點(diǎn),當(dāng)滿足R2>0.9時(shí)即停止構(gòu)造，并結(jié)合遺傳算法完成材料參數(shù)的反求。

DP590和DP780的JC模型參數(shù)取值區(qū)間設(shè)置分別如表8、表9所示。

表8　DP590的JC模型參數(shù)取值區(qū)間

表9　DP780的JC模型參數(shù)取值區(qū)間

采用解耦型GMDH和傳統(tǒng)GMDH兩種方法分別構(gòu)造近似模型,并與遺傳算法結(jié)合來(lái)反求材料參數(shù)。設(shè)置種群大小為100,迭代次數(shù)為50,樣本點(diǎn)成本及收斂結(jié)果對(duì)比如表10所示。

表10　兩種方法收斂結(jié)果對(duì)比

反求得到的8個(gè)參數(shù)對(duì)比結(jié)果如表11所示。

表11　參數(shù)反求對(duì)比結(jié)果

5　試驗(yàn)驗(yàn)證

為了檢驗(yàn)參數(shù)的可靠度,將反求的8個(gè)參數(shù)輸入有限元模型中，并將仿真結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比，圖5為碰撞試驗(yàn)的變形結(jié)果圖。

圖5　雙焊接吸能筒的臺(tái)車碰撞試驗(yàn)

分別將有限元模型的加速度與試驗(yàn)過(guò)程中的加速度在Original中繪制成圖形,對(duì)比結(jié)果如圖6和表12所示。

圖6　吸能筒碰撞仿真與試驗(yàn)對(duì)比加速度曲線

試驗(yàn)值(g)DeGMDHGMDH仿真值(g)相對(duì)誤差(%)仿真值(g)相對(duì)誤差(%)最大值28.243228.24030.010330.21526.9822平均值13.964814.79485.943516.137715.5598

進(jìn)一步探究仿真加速度結(jié)果與試驗(yàn)加速度結(jié)果之間的誤差,由近似模型技術(shù)反求得到的參數(shù)計(jì)算得到的仿真加速度與試驗(yàn)加速度之間的誤差見(jiàn)表12。

吸能筒的碰撞涉及眾多復(fù)雜的因素,計(jì)算不可避免地產(chǎn)生誤差。在綜合各個(gè)影響因素的基礎(chǔ)上,從上述結(jié)果比較可以看出,仿真加速度與試驗(yàn)加速度相差不大，且加速度曲線變化趨勢(shì)一致,加速度最大值及加速度平均值的誤差較小:將基于解耦型GMDH方法得到的材料參數(shù)輸入有限元模型并計(jì)算,與試驗(yàn)結(jié)果對(duì)比,最大加速度的相對(duì)誤差為0.0103%;加速度平均值的相對(duì)誤差為5.9435%。這說(shuō)明，材料參數(shù)精度較高，滿足工程應(yīng)用，同時(shí)驗(yàn)證了該近似模型方法在實(shí)際工程問(wèn)題中的可靠度。

6　結(jié)論

(1)在傳統(tǒng)GMDH算法基礎(chǔ)上，提出了基于高維模型表示(HDMR)方法判斷變量間耦合關(guān)系的解耦型GMDH模型：在構(gòu)造第二層活動(dòng)節(jié)點(diǎn)時(shí)利用HDMR方法判斷變量間耦合關(guān)系，根據(jù)結(jié)果智能選取耦合項(xiàng)。該模型不但繼承了GMDH自組織的進(jìn)化建模思想，還能在建模初期即可明確變量間的耦合性，提前剔除不必要的耦合項(xiàng)。

(2)由高維測(cè)試函數(shù)評(píng)價(jià)結(jié)果可以看出：對(duì)于不含耦合項(xiàng)的測(cè)試函數(shù)，解耦型GMDH在準(zhǔn)確判定兩變量之間不存在耦合性的基礎(chǔ)上剔除了耦合項(xiàng)，與傳統(tǒng)GMDH相比，可以在樣本點(diǎn)減少一半的情況下各個(gè)指標(biāo)仍能優(yōu)于傳統(tǒng)GMDH所得結(jié)果；對(duì)于含有耦合項(xiàng)的測(cè)試函數(shù)，解耦型GMDH在準(zhǔn)確判斷耦合項(xiàng)存在的基礎(chǔ)上保留了耦合項(xiàng)，此時(shí)的近似模型構(gòu)造過(guò)程與傳統(tǒng)GMDH基本一致。

(3)試驗(yàn)操作、試驗(yàn)環(huán)境、試驗(yàn)設(shè)備、試件制作及數(shù)據(jù)處理等存在大量的不確定因素,在綜合各個(gè)實(shí)際因素的基礎(chǔ)上，與試驗(yàn)結(jié)果對(duì)比表明，反求得到的參數(shù)滿足實(shí)際工程精度要求，證明了基于解耦型GMDH近似模型的材料反求方法的可行性。

[1]Zhang Y,Lai X,Zhu P,et al.Lightweight Design of Automobile Component Using High Strength Steel Based on Dent Resistance[J].Materials & Design,2006,27(1):64-68.

[2]Smerd R,Winkler S,Salisbury C,et al.High Strain Rate Tensile Testing of Automotive Aluminum Alloy Sheet[J].International Journal of Impact Engineering,2005,32(1):541-560.

[3]曹銀鋒,李光耀,鐘志華.金屬成型材料參數(shù)的反求技術(shù)[J].計(jì)算力學(xué)學(xué)報(bào),2004,6(3):291-297.

Cao Yinfeng,Li Guangyao,Zhong Zhihua.Reverse Technology of Metal Forming Material Parameters[J].Chinese Journal of Computational Mechanics,2004,6(3):291-297.

[4]Eggertsen P A,Mattiasson K.An Efficient Inverse Approach for Material Hardening Parameter Identification from a Three-point Bending Test[J].Engineering with Computers,2010,26(2):159-170.

[5]Sedighi M, Khandaei M, Shokrollahi H.An Approach in Parametric Identification of High Strain Rate Constitutive Model Using Hopkinson Pressure Bar Test Results[J].Materials Science and Engineering:A,2010,527(15):3521-3528.

[6]李恩穎,王琥,李光耀.基于支持向量機(jī)回歸的材料參數(shù)反求方法[J]. 機(jī)械工程學(xué)報(bào)，2012, 48(6):90-96.

Li Enying,Wang Hu,Li Guangyao.Material Parameter Inverse Technique Based on Support Vector Regression[J].Journal of Mechanical Engineering,2012,48(6):90-96.

[7]Simpson T W,Peplinski J D,Koch P N,et al.Metamodels for Computer-based Engineering Design:Survey and Recommendations[J].Engineering with Computer,2001,17(2):129-150.[8]Ivakhnenko A G. Heuristic Self-organization on Problems of Engineering Cybemetics[J]. Automatic,1970,6(3):207-219.[9]Wang H,Tang L,Li G Y.Adaptive MLS-HDMR Metamodeling Techniques for High Dimensional Problems[J].Expert Systems with Applications,2011,38(10):14117-14126.[10]Kondo T,Abhijit S P,Nagashino H.GMDH-type Neural Network Algorithm with a Feedback Loop for Structural Identification of RBF Neural Network[J].International Journal of Knowledge-based and Intelligent Engineering Systems.2007,11:157-168.[11]Bozdgan H.Model Selection and Akaike’s Information Criterion(AIC):the General Theory and Its Analytical Extensions[J].Psychometrika,1987,52:345-370.

[12]Bozdgan H. Akaike’s Information Criterion and Recent Developments in Information Complexity[J].Journal of Mathematical Psychology,2000,44:62-91.[13]Marcelo J C,George S D,Sahoo Debasis.A Response Surface Method-based Hybrid Optimizer[J].Inverse Problems in Science and Engineering,2008,16(6):717-741.

[14]白中浩,陳天治,曹立波,等.基于正交設(shè)計(jì)的汽車前縱梁吸能結(jié)構(gòu)的優(yōu)化[J].汽車工程,2010,32(10):935-940.

Bai Zhonghao,Chen Tianzhi,Cao Libo,et al.Optimization of Energy-absorbing Structure of Vehicle front Side Rail Based on Orthogonal Design[J].Automotive Engineering,2010,32(10):935-940.

(編輯蘇衛(wèi)國(guó))

Inverse Determination of Material Parameters Based on Decoupled GMDH Algorithm

Yin WeiyangWang HuFeng HuiTang Long

State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body,Hunan University,Changsha,410082

Since materials with strain rate effect showed different performances in crash simulation compared with they were under static conditions,accurate constitutive model parameters were the premise of reliable simulation results.The determination for the coupling relationship among the input variables was ignored during the construction program in traditional metamodel method,there would be coupled terms which were not necessary,on the contrary the proportion of the error term could be increased and efficiency and generalization ability of the metamodel were reduced.To deal with this problem,a decoupled GMDH algorithm was proposed:there would be a judgment process for the coupling relationship of all the input variables before the construction process,and the results are taken to determine whether choose the coupled term.In the test program of high dimension functions,the algorithm could provide accurate metamodel.Based on the trolley test,the method was implemented into match of the material parameters for two kinds advanced high-strength steel,and the results show that this algorithm owns good modeling performance.

strain rate effect;coupling relationship;decoupled GMDH(group method of data handling);metamodel

2013-12-20

國(guó)家自然科學(xué)基金資助項(xiàng)目(11172097,61232014);新世紀(jì)優(yōu)秀人才支持計(jì)劃資助項(xiàng)目(NCET-11-0131);湖南省自然科學(xué)基金資助項(xiàng)目(11JJA001);國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)資助項(xiàng)目(2012AA111802)

U467DOI：10.3969/j.issn.1004-132X.2015.09.016

殷為洋,男,1989年生。湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室碩士研究生。研究方向?yàn)槠嚤粍?dòng)安全與工程優(yōu)化。王琥，男,1975年生。湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室副教授、博士。馮慧,女,1989年生。湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室博士研究生。湯龍,男,1986年生。湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室博士研究生。