肖志 ?k??++劉洪斌??++孔春玉??++莫富灝??++李素雯

摘要：采用有限元仿真技術(shù)研究了泡沫鋁填充薄壁梁不同角度斜向碰撞下的吸能特性。分析了泡沫鋁密度、壁厚、長寬比等設(shè)計(jì)參數(shù)對臨界角和平均碰撞力的影響，并擬合出平均碰撞力和各設(shè)計(jì)參數(shù)及碰撞角度的解析式。結(jié)果表明，泡沫鋁填充薄壁梁的臨界彎折角小于無填充薄壁梁，并且泡沫鋁密度對臨界角的影響最大。隨著橫截面寬度從70 mm增大到90 mm，無填充薄壁梁彎折臨界角從6°增大到9°，而泡沫鋁填充薄壁梁彎折臨界角度則沒有明顯改變。壁厚越厚泡沫鋁填充薄壁梁的彎折臨界角越大，但無填充薄壁梁彎折臨界角不變。平均碰撞力和碰撞角度存在指數(shù)函數(shù)關(guān)系，且反比于長寬比平方、正比于壁厚和泡沫鋁密度的指數(shù)冪。

關(guān)鍵詞：斜向碰撞；薄壁梁；耐撞性；填充泡沫鋁；臨界角；碰撞力

中圖分類號：U467。3文獻(xiàn)標(biāo)識碼：A

金屬薄壁梁作為一種低成本、高效率的緩沖吸能結(jié)構(gòu)，被廣泛應(yīng)用于車輛、船舶、航空航天等幾乎所有交通工具的碰撞吸能裝置中＼[1＼]。現(xiàn)代汽車中大部分框架都設(shè)計(jì)成薄壁梁結(jié)構(gòu)， 通過在碰撞過程中自身的變形來吸收能量。薄壁梁的碰撞主要有軸向潰縮和彎折兩種變形模式。當(dāng)前大量的研究主要從薄壁梁的軸向碰撞進(jìn)行研究，包含材料、厚度和截面形狀等對薄壁梁碰撞特性的影響＼[2-3＼]。

但是，現(xiàn)實(shí)中的薄壁梁吸能結(jié)構(gòu)只承受純軸向碰撞載荷的形態(tài)很少發(fā)生，大多數(shù)情況下，載荷方向會偏離薄壁結(jié)構(gòu)軸線，或者與軸線方向存在一定夾角＼[4＼]。大量的實(shí)驗(yàn)和數(shù)值仿真均表明，薄壁梁在軸向載荷作用下發(fā)生漸進(jìn)壓潰變形能有效地吸收碰撞能量，而斜向碰撞下的能量吸收隨著碰撞角度的增加而急劇減少。研究表明，對特定的薄壁梁存在一個(gè)臨界角度，當(dāng)碰撞角度大于此臨界角度時(shí)，薄壁梁將發(fā)生側(cè)向彎曲，其吸能能力迅速下降＼[5＼]。Nagel等人＼[6＼]研究了變截面管件在斜向碰撞載荷下的力學(xué)響應(yīng)，包括單邊漸弱（single taper）管和雙邊漸弱（double taper）管。研究結(jié)果表明，這種漸弱的設(shè)計(jì)能有效提高薄壁梁件的吸能能力，同時(shí)雙邊漸弱比單邊漸弱更加有效。

湖南大學(xué)學(xué)報(bào)（自然科學(xué)版）2015年

第10期肖志等：泡沫鋁填充薄壁梁斜向碰撞的仿真研究

Thornton＼[7＼]研究發(fā)現(xiàn)薄壁梁通過填充其他材料的方式可以增強(qiáng)耐撞性。泡沫鋁是其中一種填充功能材料，兼有金屬和泡沫的特性，其壓縮應(yīng)力應(yīng)變曲線有高而寬的應(yīng)力平臺，大量能量在近似恒定應(yīng)力下被吸收。這種特性使泡沫金屬具有很高的吸能效率，可以廣泛應(yīng)用于汽車、軌道列車和航空航天等行業(yè)，因此，泡沫鋁填充薄壁梁引起了廣泛的關(guān)注＼[8＼]。徐雅晨等＼[9＼]研究了泡沫鋁填充鋁合金圓管的軸向壓縮性能；Hanssen等＼[10-11＼]通過大量試驗(yàn)研究了動、靜載荷作用下泡沫鋁填充薄壁梁的軸向壓潰行為；謝中友等＼[12＼]進(jìn)行了泡沫鋁合金填充圓管的三點(diǎn)彎曲實(shí)驗(yàn)研究；Yang等＼[13＼]對無填充薄壁梁和泡沫鋁填充薄壁梁在斜向加載下的吸能進(jìn)行了多目標(biāo)優(yōu)化；Santosa等＼[14＼]進(jìn)行了泡沫鋁填充薄壁梁的軸向靜態(tài)壓縮試驗(yàn)和數(shù)值研究，發(fā)現(xiàn)平均壓潰力增量和泡沫鋁抗壓強(qiáng)度、薄壁梁橫截面積成線性關(guān)系；Ahmad等＼[15＼]研究了泡沫鋁填充圓錐管在斜向動態(tài)載荷下的吸能特性。目前，研究人員從斜向準(zhǔn)靜態(tài)和動態(tài)加載兩方面研究薄壁梁，提出了評價(jià)斜向碰撞的關(guān)鍵指標(biāo)——平均碰撞力＼[2＼]。

本文通過有限元數(shù)值仿真，以平均碰撞力為衡量指標(biāo)研究不同斜向碰撞角度下泡沫鋁密度、薄壁梁壁厚、長寬比等參數(shù)對泡沫鋁填充薄壁梁吸能特性的影響。定量分析不同碰撞角度下平均碰撞力與泡沫鋁密度、長寬比、壁厚等參數(shù)的關(guān)系，為泡沫鋁填充結(jié)構(gòu)優(yōu)化設(shè)計(jì)和工程應(yīng)用提供參考。

1有限元模型

1。1有限元模型描述

考慮到接觸問題的高度非線性、大變形和材料的復(fù)雜性，本文利用顯式非線性有限元程序LSDYNA研究泡沫鋁填充薄壁梁斜向碰撞的響應(yīng)。有限元基準(zhǔn)模型如圖1所示，薄壁梁的橫截面寬度b=80 mm，總長度l=480 mm，薄壁梁壁厚t=1。5 mm，泡沫鋁密度ρ=340 kg/m3，碰撞角度為α?？紤]到在實(shí)際應(yīng)用中吸能結(jié)構(gòu)一般是連接在一個(gè)質(zhì)量相對很大的物體上，因此在薄壁梁末端附加400 kg質(zhì)量。薄壁梁以13。41 m/s的初速度斜向撞擊無摩擦剛性墻。

圖1有限元模型示意圖

Fig。1Finite element model

薄壁梁采用BELYTSCHKOLINTSAY類型殼單元，厚度方向5個(gè)積分點(diǎn)，泡沫鋁采用單點(diǎn)積分8節(jié)點(diǎn)體單元以避免體積鎖定。為了防止零能模式的產(chǎn)生，采用基于剛度的沙漏控制。剛性墻采用殼單元并使用剛性材料（MAT_20）模擬。薄壁梁和泡沫鋁的自接觸采用ATUOMATIC_SINGLE_SURFACE，薄壁梁和泡沫鋁的接觸采用AUTOMATIC_SURFACE_TO_SURFACE，靜態(tài)和動態(tài)摩擦因數(shù)分別為0。3和0。2，薄壁梁和泡沫鋁與剛性墻的接觸采用AUTOMATIC_SURFACE_TO_SURFACE。殼單元和體單元大小分別為4 mm和5 mm。

薄壁梁材料為低碳鋼，密度ρ=7 800 kg/m3，彈性模量E=205 GPa，泊松比ν=0。3，屈服應(yīng)力σy=220 MPa。使用汽車碰撞中廣泛使用的分段線性材料模型（MAT_24），其工程應(yīng)力應(yīng)變曲線參考文獻(xiàn)＼[5＼]，見表1，本文未考慮應(yīng)變率效應(yīng)對碰撞結(jié)果的影響。

表1低碳鋼應(yīng)力應(yīng)變曲線數(shù)據(jù)

Tab。1Stressstrain hardening data for mild steel

等效塑性應(yīng)變

工程應(yīng)力σ/MPa

0。0

220。0

0。005 939 3

250。0

0。016 738

300。0

0。025 750

332。0

0。035 750

352。5

0。051 500

370。0

0。105 250

374。0

由于泡沫鋁材料的復(fù)雜性，泡沫鋁模型的建立在有限元仿真中仍然是一個(gè)很大的挑戰(zhàn)。LSDYNA提供了幾種泡沫鋁模型：MAT_26，MAT_63，MAT_75，MAT_83，MAT_126和MAT_154等。文獻(xiàn)＼[16＼]表明MAT_63和MAT_75具有較高的準(zhǔn)確性和計(jì)算效率。本文選擇比較常用的MAT_63號泡沫鋁材料模型。

Hou等＼[17＼]總結(jié)了適用于任意密度泡沫鋁（MAT_63）的應(yīng)力應(yīng)變關(guān)系，但是沒有考慮致密應(yīng)變的影響，本文在此基礎(chǔ)上考慮致密應(yīng)變后泡沫鋁的應(yīng)力應(yīng)變關(guān)系，如表2所示。

表2任意密度泡沫鋁應(yīng)力應(yīng)變關(guān)系

Tab。2Stressstrain corresponding relationship

for arbitrary foam density

應(yīng)變

0

σ*/E*

0。6εD

0。7εD

0。8εD

εD

應(yīng)力 /GPa

0

σ*

σ*

1。35σ*

0。05 E*

1

表中σ*，E*，εD分別為泡沫鋁平臺應(yīng)力、彈性模量、致密應(yīng)變，主要由基體材料的性質(zhì)、泡沫材料的相對密度和胞體結(jié)構(gòu)的幾何性質(zhì)決定。泡沫材料的力學(xué)性能遵守冪法則＼[18＼]：

A（ρfoam）=A0（ρfoamρa(bǔ)l）n。 （1）

式中：A為泡沫鋁材料的性能；ρfoam和ρa(bǔ)l分別為泡沫鋁密度和鋁密度；A0為反映孔壁基體材料性能的因子；n為指數(shù)。A0和n通常使用最小二乘法擬合實(shí)驗(yàn)數(shù)據(jù)得到。本文取值參考文獻(xiàn)＼[19＼]，見表3。

表3 泡沫鋁材料常數(shù)

Tab。3Aluminium foam material parameters

σ*

E*

A0 /GPa

0。59

330

n

2。21

2。45

結(jié)構(gòu)的耐撞性主要考慮其吸能能力，一般表征指標(biāo)是平均碰撞力Fmean。能量的計(jì)算有兩種方式：一種是通過對力位移曲線積分得到，一種是直接從內(nèi)能曲線得到，兩種方式的差別在2%以內(nèi)＼[20＼]。本文直接輸出內(nèi)能數(shù)據(jù)，再計(jì)算得到平均碰撞力。

1。2 模型驗(yàn)證

為了保證有限元模型的有效性，需要對比無填充薄壁梁和泡沫鋁填充薄壁梁在斜向碰撞下數(shù)值仿真和試驗(yàn)的平均碰撞力、初始峰值和壓潰模式。無填充薄壁梁有限元模型通過文獻(xiàn)[5]中的力位移曲線和平均碰撞力結(jié)果來驗(yàn)證。但是關(guān)于泡沫鋁填充薄壁梁斜向動態(tài)加載的實(shí)驗(yàn)數(shù)據(jù)、理論分析或有限元分析結(jié)果未見報(bào)道，泡沫鋁填充薄壁梁有限元模型通過Reyes等＼[19＼]對泡沫鋁填充薄壁梁靜態(tài)斜向加載的結(jié)果來進(jìn)行驗(yàn)證。

圖2和圖3分別為無填充薄壁梁在幾個(gè)典型角度斜向碰撞下的力位移響應(yīng)曲線和在不同角度下斜向碰撞的平均碰撞力。從圖2和圖3可以看出，波形趨勢與文獻(xiàn)[19]中的結(jié)果具有一致性，同時(shí)注意到初始峰值比文獻(xiàn)結(jié)果要高一些，然而對表征本文主要研究的斜向載荷下的平均碰撞力和文獻(xiàn)吻合得很好，因而驗(yàn)證本模型達(dá)到了應(yīng)有的精度，可以用于后續(xù)的分析中。圖4對比了碰撞角度為5°時(shí)，泡沫鋁填充薄壁梁準(zhǔn)靜態(tài)壓縮的有限元仿真和試驗(yàn)結(jié)果。從變形模式來看，有限元仿真結(jié)果與文獻(xiàn)試驗(yàn)結(jié)果符合很好。同時(shí)，表4對比了碰撞角度分別為5°，15°和30°時(shí)的初始峰值力和平均碰撞力，也從另一方面表明了模型的有效性。

位移d/mm（a） 碰撞角度α=2°

位移d/mm（b） 碰撞角度α=6°

位移d/mm（c） 碰撞角度α=7°

位移d/mm（d） 碰撞角度α=11°

圖2不同碰撞角度時(shí)薄壁梁碰撞力位移曲線

Fig。2Impact forcedisplacement

curve in different angles

碰撞角度α/（°）

圖3薄壁梁平均碰撞力碰撞角度關(guān)系

Fig。3 Relationship between average

impact force and angle

圖4泡沫鋁填充薄壁梁斜向加載變形模式對比

Fig。4Comparison of oblique load deformation modes

表4不同碰撞角度時(shí)的峰值力和碰撞力

Tab。4Peak force and impact force of different

impact angleskN

5°

15°

30°

初始峰

值力

平均碰

撞力

初始峰

值力

平均碰

撞力

初始峰

值力

平均碰

撞力

本文結(jié)果

66。39

25。96

35。82

21。07

22。28

15。64

試驗(yàn)結(jié)果

66。52

27。58

36。49

20。05

22。32

15。57

相對誤差/%

1。9

5。8

1。8

4。1

1。8

4。5

2參數(shù)研究

在基準(zhǔn)模型基礎(chǔ)上，以平均碰撞力為評價(jià)指標(biāo)，研究薄壁梁壁厚、薄壁梁長寬比和泡沫鋁密度對泡沫鋁填充薄壁梁在斜向載荷下碰撞吸能特性的影響。平均碰撞力取薄壁梁末端節(jié)點(diǎn)位移為240 mm。為方便起見，薄壁梁長度不變，僅改變橫截面寬度以表示不同的長寬比。各參數(shù)取值如表5所示。

表5泡沫鋁填充薄壁梁研究參數(shù)及其取值

Tab。5Study parameters and values

泡沫鋁密度

/（kg·m-3）

橫截面寬度

/mm

長寬比

壁厚

/mm

170

70

6。86

1。3

250

75

6。40

1。5

340

80

6。00

1。7

420

85

5。65

1。9

510

90

5。33

-

2。1泡沫鋁密度的影響

基于基準(zhǔn)模型，對不同密度泡沫鋁填充薄壁梁結(jié)構(gòu)的斜向碰撞進(jìn)行研究。圖5比較了不同密度泡

沫鋁填充薄壁梁和無填充薄壁梁斜向碰撞下平均碰撞力碰撞角度的關(guān)系。隨著角度的增加，薄壁管會從疊縮變形過渡到彎曲變形，首次發(fā)生彎曲變形的角度稱為臨界角度。由圖5可知，泡沫鋁填充薄壁梁在斜向碰撞載荷下，其發(fā)生彎折的臨界角度αc小于無填充薄壁梁臨界角度（αc=7°），同時(shí)泡沫鋁密度對臨界角度的影響也很明顯。另一方面，彎折情況下泡沫鋁填充薄壁梁的平均碰撞力明顯高于無填充薄壁梁，甚至當(dāng)泡沫鋁密度較大時(shí)和無填充薄壁梁軸向壓潰時(shí)平均碰撞力相當(dāng)，因此，泡沫鋁填充薄壁梁在大斜向碰撞角度下具有良好的吸能特性。圖6為不同密度泡沫鋁填充薄壁梁發(fā)生彎折時(shí)，平均碰撞

碰撞角度α/（°）

圖5不同密度薄壁梁平均碰撞力碰撞角度關(guān)系

Fig。5Comparison of average impact

force and angle of different foam density

力和斜向碰撞角度的關(guān)系?？梢钥闯觯煌菽X密度下，平均碰撞力均隨著碰撞角度的增加而減小，其關(guān)系可由下式擬合：

Fmean=a+b×e-cα。（2）

式中：a，b，c都是與泡沫鋁密度成正相關(guān)的常數(shù)，泡沫鋁密度越大，相應(yīng)的平均碰撞力也越大，表明其吸能能力越大。擬合得到不同泡沫鋁密度的參數(shù)值見表6。

碰撞角度α/（°）

圖6不同密度薄壁梁發(fā)生彎折時(shí)

平均碰撞力碰撞角度擬合曲線

Fig。6Fitting curve of the average impact force

and the impact angle under different density

表6曲線擬合參數(shù)值

Tab。6Fitting Curve parameters

ρ/（kg·m-3）

a/kN

b/kN

c/（°）-1

170

5。792 3

15。940 2

0。076 7

250

8。079 2

18。191 8

0。111 0

340

10。545 2

20。909 4

0。147 5

420

12。120 2

25。883 4

0。174 8

510

14。722 0

26。161 6

0。180 8

2。2 寬度（長寬比）的影響

在基準(zhǔn)模型基礎(chǔ)上，研究泡沫鋁填充薄壁梁截面寬度（長寬比）對斜向碰撞吸能的影響。橫截面寬度w分別為70，75，80，85，90 mm，代表長寬比分別為6。86，6。40，6。00，5。65，5。33。圖7為不同橫截面寬度（長寬比）時(shí)，泡沫鋁填充薄壁梁和無填充薄壁梁斜向碰撞下平均碰撞力碰撞角度的關(guān)系。結(jié)果表明，泡沫鋁填充薄壁梁在斜向碰撞下吸能整體明顯高于無填充薄壁梁。由于較早地發(fā)生彎折，在小碰撞角度情況下泡沫鋁填充薄壁梁優(yōu)勢不明顯，碰撞角度越大其吸能優(yōu)勢越明顯。從圖中可以看出，隨著橫截面寬度的增加，無填充薄壁梁彎折臨界角從6°增大到9°，而泡沫鋁填充薄壁梁彎折臨界角度則沒有明顯改變。

如圖8所示，在基準(zhǔn)模型上，如果僅改變截面寬度（長寬比），平均碰撞力可表示為長寬比和碰撞角度的函數(shù)：

Fmean=（468。3+1 433。1e-0。3α）（wL）2。（3）

碰撞角度α/（°）（a） 泡沫鋁填充薄壁梁

碰撞角度α/（°） （b） 無填充薄壁梁

圖7不同寬度時(shí)薄壁梁平均碰撞力對比

Fig。7Comparison of average impact

force of different breadth

碰撞角度α/（°）

圖8不同寬度時(shí)薄壁梁彎折時(shí)平均

碰撞力碰撞角度擬合曲線

Fig。8Fitting curve of the average impact force

and the impact angle of different breadth

2。3壁厚的影響

在基準(zhǔn)模型基礎(chǔ)上，研究泡沫鋁填充薄壁梁壁厚對斜向碰撞的影響。圖9比較了不同壁厚泡沫鋁填充薄壁梁和無填充薄壁梁斜向碰撞下平均碰撞力碰撞角度的關(guān)系。無填充薄壁梁彎折臨界角保持不變，與壁厚無關(guān)；壁厚對泡沫鋁填充薄壁梁彎折臨界角有明顯的影響，且壁厚越厚彎折臨界角越大。

碰撞角度/（°） （a） 泡沫鋁填充薄壁梁

碰撞角度/（°） （b） 無填充薄壁梁

圖9不同壁厚時(shí)薄壁梁平均碰撞力對比

Fig。9Comparison of average impact

force under different wall thickness

圖10為不同碰撞角度下，泡沫鋁填充薄壁梁和無填充薄壁梁平均碰撞力的比值。在大碰撞角度時(shí)的彎折吸能和小碰撞角度時(shí)的軸向吸能，泡沫鋁填充薄壁梁都高于無填充薄壁梁，彎折時(shí)泡沫鋁填充薄壁梁平均碰撞力約是無填充薄壁梁的1。5倍。但是填充泡沫鋁薄壁梁臨界角小于無填充薄壁梁，當(dāng)碰撞角度介于泡沫鋁填充薄壁梁臨界角和無填充薄壁梁臨界角之間時(shí)，泡沫鋁填充薄壁梁吸能小于無填充薄壁梁。

碰撞角度/（°）

圖10不同壁厚時(shí)薄壁梁平均碰撞力比值

Fig。10Ratio of the average impact

force under different wall thickness

圖11為不同壁厚泡沫鋁填充薄壁梁在彎折時(shí)不同碰撞角度下的平均碰撞力，平均碰撞力與碰撞角度及壁厚之間存在如下關(guān)系：

Fmean=6。49t1。5+12。85t2e-0。146tα。（4）

碰撞角度/（°）

圖11不同壁厚薄壁梁在彎折時(shí)

平均碰撞力碰撞角度擬合曲線

Fig。11Fitting curve of the average impact force

and the impact angle under different wall thickness

3結(jié)論

本文研究了泡沫鋁填充薄壁梁在不同角度斜向碰撞時(shí)的吸能特性。采用有限元仿真技術(shù)研究了不同泡沫鋁密度、長寬比和壁厚的泡沫鋁填充薄壁梁在不同角度斜向碰撞時(shí)的平均碰撞力，通過和無填充薄壁梁對比發(fā)現(xiàn)泡沫鋁填充薄壁梁能明顯提高吸能量。本文定量地給出了泡沫鋁填充薄壁梁在不同斜向碰撞角度下平均碰撞力和泡沫鋁密度、長寬比和壁厚的公式。研究結(jié)果表明：

1）隨著碰撞角度的增加，泡沫鋁填充薄壁梁吸能量減小。在臨界角附近泡沫鋁填充薄壁梁吸能與無填充薄壁梁相當(dāng)或更小，在大碰撞角度下泡沫鋁填充薄壁梁吸能遠(yuǎn)高于無填充薄壁梁。

2）平均碰撞力隨著泡沫鋁密度的增加而增加，而且是碰撞角度的指數(shù)函數(shù)，各系數(shù)正相關(guān)于泡沫鋁密度。不同長寬比泡沫鋁填充薄壁梁斜向碰撞平均碰撞力和碰撞角度滿足同一個(gè)指數(shù)關(guān)系，各系數(shù)反比于長寬比平方。不同壁厚泡沫鋁填充薄壁梁斜向碰撞平均碰撞力和碰撞角度是指數(shù)關(guān)系，各系數(shù)正比于壁厚的指數(shù)冪。

3）泡沫鋁填充薄壁梁彎折臨界角小于無填充薄壁梁，且隨著泡沫鋁密度增加薄壁梁彎折臨界角顯著減小；長寬比對泡沫鋁填充薄壁梁彎折臨界角沒有明顯影響，但對無填充薄壁梁彎折臨界角影響明顯，壁厚則與之相反。

同時(shí)本文僅以平均碰撞力為耐撞性評價(jià)指標(biāo)，未全面考慮初始峰值力和結(jié)構(gòu)質(zhì)量的影響，具有一定的局限性。

參考文獻(xiàn)

[1]BANHART J。 Manufacture， characterization and application of cellular metals and metal foams＼[J＼]。Progress in Materials Science， 2001， 46（6）： 559-632。

[2]張維剛，黃棟。汽車薄壁梁斜向碰撞性能仿真研究＼[J＼]。 汽車工程，2010，32（6）：515-518。

ZHANG Weigang， HUANG Dong。 Simulation study of automatic thinwalled beam crashworthiness under oblique impact loading ＼[J＼]。Automotive Engineering， 2010， 32（6）： 515-518。（In Chinese）

[3]孟志強(qiáng)，何濤，尹望吾，等。六角彈簧管汽車碰撞吸能裝置的效能研究＼[J＼]。湖南大學(xué)學(xué)報(bào)： 自然科學(xué)版，2011， 38（11）：43-46。

MENG Zhiqiang， HE Tao， YIN Wangwu， et al。 Research on the performances of hexagonal spring tube energyabsorbing device for car crash＼[J＼]。Journal of Hunan University： Natural Sciences， 2011， 38（11）： 43-46。 （In Chinese）

[4]亓昌，董方亮，楊姝，等。錐形多胞薄壁梁斜向碰撞吸能特性仿真研究＼[J＼]。振動與沖擊，2012，31（24）：102-107。

QI Chang， DONG Fangliang， YANG Shu， et al。 Energyabsorbing characteristics of a tapered multicell thinwalled tube under oblique impact＼[J＼]。 Journal of Vibration and Shock， 2012， 31（24）：102-107。（In Chinese）

[5]HAN D C， PARK S H。 Collapse behavior of square thinwalled columns subjected to oblique loads＼[J＼]。 Thinwalled Structures， 1999， 35（3）：167-184。

[6]NAGEL G M， THAMBIRATNAM D P。 Computer simulation and energy absorption of tapered thinwalled rectangular tubes ＼[J＼]。 ThinWalled Structures， 2005，43（8）：1225-1242。

[7]THORNTON P H。 Energy absorption by foam filled structures＼[R＼]。 Detroit：SAE，1980。

[8]趙廣寧，王社鋒，趙洪倫。泡沫鋁填充薄壁圓管吸能特性優(yōu)化＼[J＼]。計(jì)算機(jī)輔助工程，2009，18（4）：62-66。

ZHAO Guangning， WANG Shefeng， ZHAO Honglun。 Optimization on energy absorption performance of thinwalled pipe filled with aluminum foam＼[J＼]。Computer Aided Engineering， 2009，18（4）：62-66。（In Chinese）

[9]徐雅晨，鳳儀，湯靖婧，等。 泡沫鋁填充薄壁鋁合金圓管軸向壓縮性能的數(shù)值模擬＼[J＼]。稀有金屬材料與工程，2011，40（5）：875-879。

XU Yachen， FENG Yi， TANG Jingjing， et al。 Numerical simulation of vertical compressing experiments of thinwalled cylindrical aluminum alloy tubes filled with aluminum foam＼[J＼]。Rare Metal Materials and Engineering， 2011，40（5）：875-879。（In Chinese）

[10]HANSSEN A G， LANGSETH M， HOPPERSTAD O S。 Static and dynamic crushing of circular aluminum extrusions with aluminum foam filled＼[J＼]。International Journal of Impact Engineering， 2000，24（5）：475-507。

[11]HANSSEN A G， HOPPERSTAD O S， LANGSETH M，et al。 Validation of constitutive models applicable to aluminum foam＼[J＼]。 International Journal of Mechanical Sciences， 2002， 44（2）：359-406。

[12]謝中友，虞吉林，李劍榮。泡沫鋁合金填充圓管三點(diǎn)彎曲實(shí)驗(yàn)研究＼[J＼]。實(shí)驗(yàn)力學(xué)，2007，22（2）：104-110。

XIE Zhongyou， YU Jilin， LI Jianrong。 An experimental study on threepoint bending of aluminum alloy foamfilled cylindrical aluminum alloy pipe＼[J＼]。Journal of Experimental Mechanics，2007，22（2）：104-110。（In Chinese）

[13]YANG S， QI C。 Multiobjective optimization for empty and foamfilled square columns under oblique impact loading ＼[J＼]。 International Journal of Impact Engineering， 2013， 54（4）：177-191。

[14]SANTOSA S P， WIERZBICKI T， HANSSEN A G， et al。 Experimental and numerical studies of foamfilled sections ＼[J＼]。International Journal of Impact Engineering， 2000， 24（5）：509-534。

[15]AHMAD Z， THAMBIRATNAM D P， TAN A C。 Dynamic energy absorption characteristics of foamfilled conical tubes under oblique impact loading＼[J＼]。International Journal of Impact Engineering，2010，37（5）：475-488。

[16]HANSSEN A G， LANGSETH M， HOPPERSTAD O S。 Static and dynamic crushing of square aluminum extrusions with aluminum foam filler＼[J＼]。International Journal of Impact Engineering， 2000， 24（4）：347-383。

[17]HOU S， LI Q， LONG S， et al。 Crashworthiness design for foam filled thinwall structures ＼[J＼]。Materials & Design， 2009， 30（6）：2024-2032。

[18]GIBSON L J， ASHBY M F。 Cellular solids： structure and properties ＼[M＼]。Cambridge： Cambridge University Press， 1999：152-205。

[19]REYES A， HOPPERSTAD O S， LANGSETH M。 Aluminum foamfilled extrusions subjected to oblique loading： experimental and numerical study ＼[J＼]。International Journal of Solids and Structures， 2004， 41（5）：1645-1675。

[20]REID J D。 Crashworthiness of automotive steel midrails： Thickness and material sensitivity ＼[J＼]。Thinwalled Structures， 1996， 26（2）：83-103。