許 平， 邵 恒， 嚴佳麗

(中南大學(xué) 交通運輸工程學(xué)院軌道交通安全實驗室， 長沙 410075)

基于自適應(yīng)響應(yīng)面法的膨脹式吸能結(jié)構(gòu)耐撞性優(yōu)化設(shè)計

許 平， 邵 恒， 嚴佳麗

(中南大學(xué) 交通運輸工程學(xué)院軌道交通安全實驗室， 長沙 410075)

借鑒多道次漸進成形方法，提出一種膨脹式吸能新型頂桿。保持錐形頂桿最大直徑不變，在錐形頂桿區(qū)域添加定徑段，將一次高膨脹率的膨脹過程變成多次低膨脹率膨脹過程。以不同膨脹角下的錐形變形區(qū)數(shù)量和定徑段長度為設(shè)計變量，有限行程下的能量吸收為目標，采用全因子試驗設(shè)計方法，基于移動最小二乘法建立的近似模型，進行基于響應(yīng)面法優(yōu)化設(shè)計，得到不同角度下膨脹式吸能頂桿的最佳方案。結(jié)果表明，近似模型的擬合精度較高，采用該模型的優(yōu)化方案可提高膨脹式吸能管的耐撞性能。優(yōu)化的臨界角度為20.782°，誤差小于0.313°。

膨脹吸能； 膨脹角度； 響應(yīng)面法； 優(yōu)化設(shè)計

隨著軌道交通的迅速發(fā)展，列車的運行速度逐步提高，碰撞事故造成的后果趨于嚴重。因而，車輛的被動保護技術(shù)顯得尤為重要，研制耐沖擊吸能列車是眾多學(xué)者關(guān)注的重點。膨脹式吸能裝置由于其良好的耐撞性能越來越受到人們的青睞，目前已運用到車鉤緩沖裝置中[1-2]。膨脹式吸能裝置具有變形穩(wěn)定及易控制的優(yōu)點，一般在速度為15 km/h以上的碰撞事故中發(fā)生不可逆變形，從而吸收沖擊動能[3]。此外，車鉤后置式的膨脹壓潰吸能管也得到了運用，其穩(wěn)定變形力可達2 000 kN[4]。

隨著膨脹式吸能裝置的廣泛應(yīng)用，許多學(xué)者對影響膨脹管耐撞性能的關(guān)鍵參數(shù)進行了研究。羅昌杰等[5]研究了2A12O、2A12T4等鋁管在不同膨脹率及膨脹角度下的吸能特性，并提出對緩沖器不定徑部分進行預(yù)處理可提高理想吸能效率。夏茜等[6]使用有限元仿真方法研究了膨脹角度及管壁厚度對緩沖器容量的影響，提出了從能量要求、撞擊力大小及安裝空間條件綜合考慮的結(jié)構(gòu)設(shè)計方法。羅玗琪[7]建立了B型地鐵車鉤裝置的膨脹管有限元模型，進行了整車碰撞性能研究。馬彥婷[8]利用實驗、理論分析和有限元仿真相結(jié)合的方法對擴徑量、壁厚及錐角等影響因素展開了研究。Jakirahemed等[9]利用有限元仿真方法研究了不同溫度下多種膨脹角度及膨脹率的膨脹管的吸能特性和峰值力，提出膨脹率是影響峰值力及吸能特性的最重要因素。Almeida等[10]用實驗證明了大的膨脹角度和膨脹率容易引起膨脹管的破裂、起皺和局部屈曲。Shakeri等[11]用實驗和有限元相結(jié)合的方法研究了不同膨脹角度及摩擦條件下膨脹管的耐撞性能，實驗結(jié)果還證明了膨脹管膨脹后與膨脹頂桿之間存在一定的間隙。Yang等[12]基于準靜態(tài)實驗研究了不同管壁厚度及膨脹角度對穩(wěn)定變形力和吸能效率的影響，其實驗結(jié)果表明了不同階段膨脹管的變形與變形力的關(guān)系。Choi等[13-14]對膨脹角度對變形能和摩擦耗能及真實膨脹率的影響做了研究，并用沖擊實驗驗證了有限元模型的準確性。隨后，其又對不同沖擊速度和潤滑條件對膨脹吸能特性的影響展開了深入研究，得到了不同沖擊條件下膨脹管形變與沖擊力的關(guān)系。YAN等[15]將膨脹式吸能結(jié)構(gòu)的變形分為以下幾個階段：膨脹管與錐形頂桿開始接觸，變形端發(fā)生彎曲，產(chǎn)生變形抗力；膨脹管變形端在膨脹區(qū)域進行徑向的擴張，變形抗力快速上升；膨脹管變形端完成徑向擴張，進入定徑段，管材發(fā)生反向彎曲，變形抗力再次上升。

膨脹管變形吸能屬于金屬塑性大變形，容易產(chǎn)生破裂、屈曲等現(xiàn)象。金屬管材在塑性成形過程中，為提高成形質(zhì)量，采用多道次漸進成形方法。多道次成形的質(zhì)量比單道次更好，且金屬厚度分布均勻；另一方面，多道次成形可以增大金屬材料的塑性變形極限，延遲裂紋的產(chǎn)生[16-17]。

現(xiàn)有針對膨脹式吸能結(jié)構(gòu)的研究，基本是針對膨脹角度、膨脹管壁厚度、膨脹率等影響參數(shù)的單變量研究，未考慮變量之間的相互影響。變量的樣本數(shù)也比較單一，一般為個位數(shù)，沒有呈現(xiàn)變化范圍內(nèi)不同取值水平下變量組合的影響規(guī)律。另外，對提高吸能效率的研究也僅限于對不定徑部分進行預(yù)處理，沒有針對錐形頂桿的形狀提出較大改進的研究?；诖?，本文依據(jù)膨脹管在錐形頂桿作用下的變形力變化規(guī)律，結(jié)合塑性成形中多道次成形理論，提出一種錐形變形區(qū)分段的新型頂桿。在不改變錐形頂桿最大半徑的前提下，在頂桿錐形區(qū)域中增加定徑段，將整個變形過程分多次完成，提高吸能結(jié)構(gòu)變形穩(wěn)定性。本文以膨脹角度、定徑段長度和錐形變形區(qū)個數(shù)為變量，運用全因子法進行試驗設(shè)計，建立400個不同的新型頂桿樣本，基于移動最小二乘法創(chuàng)建各響應(yīng)量的近似模型。考慮綜合幾個因素的影響規(guī)律，基于此近似模型，以最大吸能量為目標，在變量范圍內(nèi)，運用自適應(yīng)響應(yīng)面法求出不同角度下最佳新型頂桿方案。

1 有限元模型

1.1 模型及參數(shù)

車鉤壓潰裝置由膨脹管和膨脹頂桿組成，膨脹管在膨脹頂桿的作用下發(fā)生塑性變形而吸收能量，膨脹頂桿一端為圓臺形狀?；诂F(xiàn)有車鉤壓潰管的形狀與尺寸，結(jié)合膨脹管膨脹行程與變形抗力關(guān)系，提出了一種新型的膨脹頂桿，頂桿結(jié)構(gòu)參數(shù)如圖1所示(本文中不分段頂桿為初始頂桿，分段頂桿為新型頂桿)。在不改變錐形頂桿最大半徑(D=160 mm)和最小半徑(d=140 mm)的前提下，在頂桿錐形區(qū)域中增加定徑段，將錐形變形區(qū)分為均等的n(n=3)段，每段定徑段的長度為l(l=20 mm)，分段后的新型頂桿有3段擴徑段，中間2段定徑段，3段擴徑段總長度與初始頂桿擴徑段長度一致。以此實現(xiàn)“擴徑-定徑”的多次循環(huán)的變形過程，將一次高膨脹率的膨脹過程變成多次低膨脹率的膨脹過程。

圖1 膨脹頂桿幾何結(jié)構(gòu)示意圖

基于HyperMesh建立了膨脹式吸能管的有限元模型，采用顯式有限元LS-DYNA進行碰撞仿真，計算選用八節(jié)點實體單元。綜合考慮計算的效率和精度的要求，采用2 mm×2 mm×2 mm的網(wǎng)格，因此，管壁厚度方向劃分為4個單元，便于分析管材厚度方向的金屬流動趨勢。膨脹管的擴徑變形屬于大變形，材料為具有良好強度和塑性的材料，而頂桿在整個變形過程中幾乎無變化，因此膨脹管選用彈塑性材料模型，膨脹頂桿選用剛性材料模型。其中，膨脹管材料的真實應(yīng)力應(yīng)變參數(shù)如表1所示。膨脹管與頂桿之間添加固體潤滑劑，動摩擦因數(shù)為0.14。采用“AUTOMATIC_SURFACE_TO_SURFACE”接觸算法，使用移動剛性墻模擬沖擊小車，剛性墻質(zhì)量為30 t，初始速度為15 km/h。采用初始頂桿和新型錐形頂桿的膨脹式吸能結(jié)構(gòu)有限元模型如圖2所示。其中，在初始頂桿的改進過程中，膨脹管的結(jié)構(gòu)尺寸保持不變：初始總長度L為350 mm，初始內(nèi)徑d1為142 mm，管壁厚度t為7 mm。

本文所采用的計算有限元模型方法與文獻[15]中一致，此種有限元計算建模方法模擬膨脹式吸能結(jié)構(gòu)的沖擊性能已在特定的沖擊實驗下證明具有較高的可信度，實驗與有限元仿真變形比較結(jié)果如圖3所示。

表1 膨脹管材料的真實應(yīng)力應(yīng)變關(guān)系

圖2 有限元模型

1.2 耐撞性指標

吸能結(jié)構(gòu)的耐撞性指標包括：① 碰撞過程中所能吸收的能量；② 碰撞過程中的初始峰值力；③ 材料的利用效率；④ 結(jié)構(gòu)的有效壓縮距離。因本文中膨脹管不作改變，其質(zhì)量保持不變。因此，本文以吸能量為目標，平均變形力為約束，對該膨脹式吸能管進行耐撞性優(yōu)化。

吸能結(jié)構(gòu)在碰撞過程中所吸收的總能量EA可表示為

(1)

平均變形力是吸能結(jié)構(gòu)在整個碰撞過程中載荷力的平均值：

(2)

2 吸能結(jié)構(gòu)耐撞性分析

對頂桿優(yōu)化前后的仿真計算結(jié)果進行分析，圖4是膨脹管在膨脹過程中的力-位移曲線。從圖中可以看出，在初始頂桿的作用下，膨脹管變形大致可分為三個階段：初始時刻的彎曲，變形區(qū)內(nèi)的膨脹和定徑區(qū)內(nèi)的反向彎曲。在新型頂桿的作用下，膨脹管變形可分為五個階段。與初始頂桿相比，由于變形區(qū)之間的定徑段長度略短，其間的彎曲變形和反向彎曲變形均不明顯。從圖4中d到e的變化可知，反向彎曲行程至少需要40 mm的行程。因此其五個階段為：初始時刻彎曲，第一段變形區(qū)內(nèi)膨脹，第二段變形區(qū)內(nèi)膨脹，第三段變形區(qū)內(nèi)膨脹，第三段定徑區(qū)內(nèi)反向彎曲。在每段彎曲及膨脹后，載荷都稍有回落，是因為膨脹管初始變形端處于自由狀態(tài)。初始頂桿作用下，膨脹管膨脹的最終穩(wěn)定膨脹變形力為343.48 kN。新型頂桿作用下，膨脹管膨脹的最終穩(wěn)定膨脹變形力為419.70 kN。后者與前者相比，增加了22.19%。

圖4 力-位移曲線

圖5是膨脹管在初始頂桿及新型頂桿作用下的變形圖，變形序列分別對應(yīng)圖3中字母標注的對應(yīng)狀態(tài)，并且與此前所描述的變形過程相對應(yīng)。從圖中可以看出，在兩種不同頂桿的作用下，膨脹管都有序地進行徑向擴張。圖左顯示的是膨脹管的最大等效應(yīng)變，比較兩者可以發(fā)現(xiàn)，最大等效應(yīng)變變化不大，后者比前者略有增加。

表2顯示了兩種不同頂桿作用下膨脹管的能量吸收。從表中可以看出，初始頂桿作用下膨脹管的吸能量為104.236 kJ，新型頂桿作用下膨脹管的吸能量為117.684 kJ，比前者增加了12.90%。膨脹式的吸能結(jié)構(gòu)耗散能量主有兩個途徑：塑性變形和摩擦耗能。與初始頂桿相比，新型頂桿的變形區(qū)長度較長，因此，當(dāng)350 mm行程結(jié)束時，新型頂桿下膨脹管的摩擦耗能量比前者多，而塑性變形能較少。

3 膨脹頂桿參數(shù)研究

(a) 初始頂桿

(b) 新型頂桿

表2 吸能量比較

從上一節(jié)的分析中可知，此種新型頂桿有利于提高膨脹管的耐撞性能。因此，在保持頂桿最大和最小直徑(D和d)不變的前提下，在新型頂桿現(xiàn)有形狀的基礎(chǔ)上做以下設(shè)計：① 改變膨脹頂桿錐角ɑ；② 改變定徑段的長度l：5 mm～50 mm；③ 改變錐形變形區(qū)數(shù)量n。分別研究這三個變量對變形力及能量吸收的影響?？紤]到要使定徑段的段數(shù)為整數(shù)，膨脹管的總長度為350 mm，膨脹頂桿不宜太長?？傻玫饺绫?所示的膨脹角度、定徑段長度及錐形變形區(qū)數(shù)量。

表3 新型膨脹頂桿各變量的取值范圍

從表3中看出，一共有三個設(shè)計變量，其中膨脹頂桿的錐角ɑ有8個取值水平，定徑段長度有10個取值水平，錐形變形區(qū)數(shù)量有5個取值水平。對三個變量進行不同水平的全因子正交計算，得到400種不同結(jié)構(gòu)形態(tài)的新型頂桿。此外，還有8種不同角度的初始頂桿。

3.1 不同角度初始頂桿作用下膨脹管的耐撞性能

從圖4膨脹管的變形序列圖中可以看出，初始變形端為自由端。在自由端越過頂桿錐形變形區(qū)時，由于具有徑向速度，還保持向外擴張的趨勢，當(dāng)徑向速度減為零時，其與頂桿之間的間隙最大，變形力實現(xiàn)單個較完整的上升過程。

表4是在當(dāng)前膨脹管參數(shù)及潤滑條件下，在不同角度的初始頂桿作用下的能量吸收及最大間隙，便于與不同的新型頂桿方案結(jié)果進行比較。錐形變形區(qū)的長度隨著角度的增大而減小，因此膨脹管與頂桿的接觸面積減小。因此，隨著角度的增大，膨脹管的塑性變形能增大，摩擦耗能減小，而總吸能量先減小后增大。

表4 初始頂桿作用下的膨脹管耐撞性能

3.2 不同角度新型頂桿作用下的耐撞性能

3.2.1 基于移動最小二乘法的總吸收能量近似模型

移動最小二乘法是基于樣本點構(gòu)造的方法，是對傳統(tǒng)加權(quán)最小二乘法技術(shù)的推廣，其與傳統(tǒng)法本質(zhì)的區(qū)別是在試驗設(shè)計中各個采樣點處的加權(quán)系數(shù)是隨著樣本點到取值點的距離改變而變化的。在最小二乘法中，求解最小二乘法擬合的問題等價于求解正則方程：

式中:A是由試驗設(shè)計(DOE)輸入值構(gòu)建的(n×p)矩陣；β為回歸系數(shù)向量(p×1)；W是一個n階對角陣，主元素為加權(quán)系數(shù)，各個主元素值的大小代表了各個采樣點的重要程度。

最小二乘法的權(quán)函數(shù)ω(x)不是定值，其會隨著取值點的遠離而衰減。權(quán)函數(shù)的幅值隨著x移動而改變，權(quán)值ωi是x的函數(shù)，因此最終擬合得到的多項式系數(shù)β也與x相關(guān)，因此采用最小移動二乘法擬合時，無法獲得具體近似函數(shù)的解析形式。

因此在本節(jié)中，基于全因子試驗設(shè)計所得的400個樣本，采用三維響應(yīng)面對近似模型的變化規(guī)律進行呈現(xiàn)，如圖6所示。從圖中可以發(fā)現(xiàn)，小角度及大角度下的響應(yīng)面形態(tài)有較大的差別。當(dāng)膨脹角度低于20°時，隨著角度的增大，總吸收能量最大值向較小n值和較大l值的位置移動。當(dāng)膨脹角度達到25°以上時，總吸能量最大值在較小n值和較小l值的位置。在膨脹管為350 mm時，膨脹管的總吸能量與分段段數(shù)(n)及定徑段長度(l)不成正比。

表5列出了近似模型與有限元模型計算之間的誤差。其中，只列出了400個樣本中誤差絕對值高于1%的模型。從表中可以得到，近似模型的最大誤差為2.316%，近似模型擬合效果良好。

表5 近似模型的誤差

3.2.2 基于自適應(yīng)響應(yīng)面法的最佳方案優(yōu)化

在自適應(yīng)響應(yīng)面法中，目標函數(shù)和約束函數(shù)按照以下二階多項式進行擬合：

式中:m為約束的個數(shù)；n為設(shè)計變量的個數(shù)；aj0、an、am分別為二次項的系數(shù)。

二次項系數(shù)是由之前試驗設(shè)計所得樣本點通過最小二乘法擬合得到，擬合過程原理如下：

(1) 分析初始的設(shè)計變量生成的n個新的設(shè)計點。

(2) 使用最小二乘法確定目標函數(shù)和各約束函數(shù)的多項式系數(shù)。

(3) 用數(shù)學(xué)規(guī)劃法對基于響應(yīng)面的問題進行優(yōu)化求解。

(4) 對模型近似最優(yōu)解進行分析。

(5) 如已經(jīng)收斂，則停止求解。

(6) 如果不收斂，則返回步驟(2)。

基于上一節(jié)樣本所建立的近似模型，以總吸收能量最大化為目標，以變形力及變形中的最大等效應(yīng)變?yōu)榧s束，得到不同角度下的最佳設(shè)計方案，如表6所示，并將其與初始頂桿方案進行優(yōu)劣對比。其中優(yōu)化比例是指新型頂桿的總吸收能量相較于初始頂桿的增加比例。

從表6中可以看出，當(dāng)膨脹角度小于20°時，膨脹管在新型頂桿的作用下耐撞性能有明顯改進。其中當(dāng)膨脹角度為10°時，總吸能量增加24.59%。在總能量增加的最佳方案中，膨脹管膨脹時的最大變形力均大于初始頂桿作用下的變形力。但是，隨著角度的改變，最大變形力變化較小。當(dāng)膨脹角度大于20°時，最佳新型頂桿方案作用下的膨脹管的耐撞性能與初始頂桿相比無明顯改進，且所吸收的能量比初始頂桿作用下的小。這是由于在膨脹管長度有限的前提下，膨脹管后部大部分都處于較低膨脹率下的膨脹狀態(tài)，只有前端變形端達到預(yù)定膨脹率。

表6 不同角度下的最佳方案

在較小的膨脹角度下，當(dāng)錐形變形區(qū)分為4段時，方案最佳；當(dāng)膨脹角度增大時，錐形變形區(qū)分為3段時，方案最佳；當(dāng)在較大的膨脹角度下，錐形變形區(qū)分段不能改進膨脹管的耐撞性能。錐形變形區(qū)的分散，使得大角度膨脹下的總彎曲變形減小。從表4中可以看出，當(dāng)角度越大時，膨脹管自由端與頂桿之間的最大徑向距離也越大。然而此優(yōu)化模型中，頂桿的最小與最大直徑是定值，他們之間的差值是10 mm。當(dāng)膨脹管經(jīng)過一次膨脹后，與頂桿之間的徑向間距過大時，會使得第二錐形變形段不起作用，膨脹管徑向擴大效應(yīng)明顯減弱，從而大大降低膨脹管的實際膨脹率。

3.2.3 臨界角度

臨界角度是對于已知尺寸的膨脹管，隨著膨脹角度的增大，最佳新型頂桿的總吸能量開始小于初始頂桿總吸收能量的角度。圖7是基于本文的初始頂桿，不同膨脹角度下的的新型頂桿吸能量優(yōu)化比例。從圖中可以看出，曲線的零點位于20°～25°之間，因此，本文采用二分法求臨界角度。當(dāng)膨脹角度固定時，以錐形變形區(qū)(n)和定徑段長度(l)為變量，以最大吸收能量為目標，基于3.1節(jié)中所得近似模型，使用自適應(yīng)響應(yīng)面法，可求出不同角度下的最佳方案。建立最佳方案的有限元模型，進一步求出優(yōu)化比例。此外，為了使曲線的駐點附近更加平滑，使用相同的方法計算駐點附近的優(yōu)化比例。當(dāng)膨脹角度為11.25°時，新型頂桿膨脹管吸能量比初始頂桿提高了24.92%。

圖7 不同膨脹角度下的的新型頂桿吸能量優(yōu)化比例

角度區(qū)間/°中間角度/(°)總能量EA/kJ新型頂桿初始頂桿優(yōu)化比例/%精確度/(°)(20．000，25．000)22．500135．994143．823-5．445．000(20．000，22．500)21．250135．732137．588-1．352．500(20．000，21．250)20．625135．659134．7400．681．250(20．625，21．250)20．938135．664136．172-0．370．625(20．625，20．938)20．782135．686135．4710．160．313

從表7中可得，膨脹頂桿的臨界角度約為20.782°，誤差小于0.313°。從而可以得到，當(dāng)角度大于20.782°時，對于當(dāng)前的膨脹管，定徑段分段優(yōu)化無法使膨脹管的耐撞性能提高。

4 結(jié) 論

(1) 借鑒多道次漸進塑性成形理論，膨脹頂桿錐形變形區(qū)的分段可改進膨脹管的耐撞性能，在一定程度使膨脹管的變形質(zhì)量提升，吸能更加穩(wěn)定。

(2) 膨脹變形力隨著分段數(shù)量及定徑段長度的增加而增大。不同角度下膨脹管總吸能量與錐形變形區(qū)數(shù)量和定徑段長度的響應(yīng)面形態(tài)不一致，且不成正比?；谝苿幼钚《朔ǖ慕颇Ｐ妥畲笳`差為2.316%。

(3) 當(dāng)膨脹頂桿錐角為11.25°時，最佳新型頂桿方案(n=3,l=45 mm)作用下的膨脹管，總吸能量的比初始頂桿提高了24.92%。膨脹頂桿的臨界角度約為20.782°，誤差小于0.313°。此種優(yōu)化設(shè)計方法為膨脹式吸能結(jié)構(gòu)設(shè)計提供了新的思路。

[1] 王學(xué)亮，姜曉東. 車鉤緩沖裝置用過載保護吸能裝置綜述[J]. 機車車輛工藝，2011(3)：14-16.

WANG Xueliang, JIANG Xiaodong. The review of overload protection energy absorption of coupling buffer device[J]. Locomotive & Rolling Stock Technology, 2011(3): 14-16.

[2] 鄭偉，魯青君，陳書翔. 柴田式車鉤緩沖裝置過載保護功能研究[J]. 鐵道機車車輛，2014, 34(5)：28-34.

ZHENG Wei, LU Qingjun, CHEN Shuxiang. The research of the function of energy absorption coupling buffer device[J]. Railway Locomotive & CAR, 2014,34(5): 28-34.

[3] 馬青春. 城鐵用吸能裝置結(jié)構(gòu)淺析[J]. 鐵道車輛，2013,51(1)：32-34.

MA Qingchun. Analysis of energy absorption device of commuter trains[J]. Rolling Stock,2013,51(1): 32-34.

[4] 楊杰，劉鳳剛. 后置式壓潰裝置的發(fā)開研究[J]. 鐵道車輛，2014, 52(6)：31-34.

YANG Jie, LIU Fenggang. The development research of rear-mounted energy absorption device[J]. Rolling Stock, 2014,52(6): 31-34.

[5] 羅昌杰，劉榮強，鄧宗全，等. 薄壁金屬管塑性變形緩沖器吸能特性的試驗研究[J]. 振動與沖擊，2010,29(4)：101-106.

LUO Changjie, LIU Rongqiang, DENG Zongquan, et al. The experimental study of energy absorption characteristics of thin walled metal tube buffer[J]. Journal of Vibration and Shock, 2010, 29(4): 101-106.

[6] 夏茜,劉國偉,羅玗琪. 鼓脹式壓潰吸能結(jié)構(gòu)的碰撞性能研究[J]. 城市軌道交通研究，2013,11:20-24.

XIA Xi, LIU Guowei, LUO Yuqi. On crashworthiness of bulgy type energy absorb structure[J]. Urban Mass Transit, 2013, 11:20-24.

[7] 羅玗琪. B型地鐵車輛間鼓脹與誘導(dǎo)組合式吸能結(jié)構(gòu)碰撞性能研究[D]. 長沙：中南大學(xué)，2011.

[8] 馬彥婷. 脹管式緩沖器吸能特性及優(yōu)化設(shè)計研究[D]. 秦皇島：燕山大學(xué)，2014.

[9] VENKATESWARLU G. A study on effect of process parameters on the expansion of thin walled aluminum 7075 tubes[J]. Int. J. Adv. Sci. Technol, 2011,36:83-93.

[10] ALMEIDA B P P, ALVES M L, ROSA P A R, et al. Expansion and reduction of thin-walled tubes using a die: experimental and theoretical investigation[J]. Int. J. Mach Tools Manuf 2006,46:1643-52.

[11] SHAKERI M, SALEHGHAFFARI S, MIRZAEIFAR R. Expansion of circular tubes by rigid tubes as impact energy absorbers: experimental and theoretical investigation[J]. International Journal of Crashworthiness, 2007,12(5):493-501.

[12] YANG J L, LOU M, HUA Y, et al. Energy absorption of expansion tubes using a conical-cylindrical die: experiments and numerical simulation[J]. Int. J. Mech. Sci., 2010,52(5):716-725.

[13] CHOI W M, KIM J S, JUNG H S, et al. Effect of punch angle on energy absorbing characteristics of tube-type crash elements[J]. Int. J. Automot Technol, 2011, 12(3): 383-389.

[14] CHOI W M, KWON T S, JUNG H S,et al. Influence of impact velocity on energy absorption characteristics and friction coefficient of expansion tube[J]. International Journal of Crashworthiness, 2012, 17(6): 1-9.

[15] YAN J L, YAO S G, XU P, et al. Theoretical prediction and numerical studies of expanding circular tubes as energy absorbers[J]. International Journal of Mechanical Sciences, 2015, 105:206-214.

[16] 金玲玲，高霖，史曉凡,等. 多道次單點漸進成形變形規(guī)律實驗研究[J]. 中國機械工程，2014, 25(2)：262-266.

JIN Lingling, GAO Lin, SHI Xiaofan, et al. Experimental study on deformation law of multi-stage single point incremental forming[J]. China Mechanical Engineering, 2014, 25(2): 262-266.

[17] 李軍超，耿佩，潘俊杰. 金屬板材多道次漸進成形的模擬分析[J]. 熱加工工藝，2014,41(13): 79-84.

LI Junchao, GENG Pei, PAN Junjie. Numerical simulation of multistage incremental forming in sheet metal[J]. Hot Working Technology, 2014, 41(13): 79-84.

Crashworthiness optimization design of expanding type energy absorption devices based on adaptive response surface method

XU Ping, SHAO Heng, YAN Jiali

(Key Laboratory of Traffic Safety on Track of Ministry of Education, School of Traffic and Transportation Engineering, Central South University, Changsha 410075, China)

A new expanding type energy absorption device was proposed based on the multistage incremental forming method. The device’s conical mandrel was divided into several sections with certain diameters and the largest diameter of the mandrel kept unchanged. Thus, the expansion process with a high expansion ratio was changed into several expansion processes with lower expansion ratios. Taking the number (n) of conical deformation zones and the length (l) of sections with certain diameters under different expansion angles as design variables, taking the energy absorption under a finite stroke as an objective function, the full factorial test design method and an approximate model built based on the moving least squares method (MLAM) were adopted to conduct an optimization design based on the adaptive response surface method (ARSM). The optimal design schemes of expanding type energy absorption mandrels under different expansion angles were obtained. The results showed that the approximation model has a higher fitting accuracy, the optimization design schemes based on this approximate model can improve the crashworthiness of expanding type energy absorption tubes; the optimal critical expansion angle is 20.782°, its error is less than 0.313°.

expanding type energy absorption; expansion angle; adaptive response surface method; optimization design

教育部重點項目(113051A)；國家支撐項目(2015BAG12B01)

2015-12-22 修改稿收到日期：2016-04-08

許平 男，博士，副教授，1971年

U270.34

A

10.13465/j.cnki.jvs.2017.11.018