功能梯度泡沫薄壁管斜向沖擊耐撞性仿真研究

趙博文，仲衍慧，莊夢(mèng)夢(mèng)

(中國(guó)汽車技術(shù)研究中心汽車試驗(yàn)研究所 天津300300)

功能梯度泡沫薄壁管斜向沖擊耐撞性仿真研究

趙博文，仲衍慧，莊夢(mèng)夢(mèng)

(中國(guó)汽車技術(shù)研究中心汽車試驗(yàn)研究所 天津300300)

采用有限元分析方法，以比吸能和初始碰撞力峰值為評(píng)價(jià)指標(biāo)，通過與等質(zhì)量均值密度泡沫填充管對(duì)比，研究了一種新型功能梯度泡沫填充管在不同沖擊角度下的失效模式和耐撞性能，分析了包括泡沫梯度指數(shù)、壁厚和泡沫密度變化范圍在內(nèi)的結(jié)構(gòu)參數(shù)對(duì)其在斜向沖擊下耐撞特性的影響。結(jié)果表明，功能梯度泡沫填充管在斜向沖擊下的耐撞性能明顯優(yōu)于均值密度泡沫填充管；結(jié)構(gòu)參數(shù)之間相互關(guān)聯(lián)、彼此影響，均對(duì)功能梯度泡沫填充薄壁管的耐撞特性產(chǎn)生顯著影響。

薄壁管 功能梯度泡沫 吸能 耐撞性 斜向沖擊

薄壁管結(jié)構(gòu)在軸向沖擊下容易產(chǎn)生穩(wěn)定的漸進(jìn)式壓潰失效，通過塑性屈曲變形吸能，其作為一種高效吸能結(jié)構(gòu)，被廣泛應(yīng)用于車輛、船舶、航空航天等工業(yè)領(lǐng)域。針對(duì)金屬薄壁結(jié)構(gòu)的耐撞特性，國(guó)內(nèi)外進(jìn)行了大量研究，旨在提高該結(jié)構(gòu)在碰撞過程中的吸能效率，并實(shí)現(xiàn)對(duì)變形模式的有效控制。[1-3]

近年來，泡沫鋁填充薄壁管開始受到關(guān)注。[4-6]相比中空結(jié)構(gòu)，泡沫鋁填充管能更好地保持漸進(jìn)壓潰的變形模式，以此提高吸能穩(wěn)定性。Hassan等進(jìn)行了大量泡沫鋁填充薄壁管在準(zhǔn)靜態(tài)和動(dòng)態(tài)載荷作用下的沖擊實(shí)驗(yàn)，并據(jù)實(shí)驗(yàn)結(jié)果，構(gòu)造出碰撞中平均沖擊力的經(jīng)驗(yàn)公式；[7]Chen 和Wierzbicki構(gòu)造出中空和泡沫填充的多胞管在軸向沖擊下的理論變形值，并發(fā)現(xiàn)泡沫鋁和薄壁管間的相互作用效應(yīng)是提升此類結(jié)構(gòu)吸能特性的主要原因；[8]侯淑娟也以泡沫密度和壁厚等結(jié)構(gòu)參數(shù)為研究變量，對(duì)正方形截面泡沫鋁薄壁管的吸能特性進(jìn)行了優(yōu)化設(shè)計(jì)。[9]

研究表明，泡沫填充管的吸能特性與填充泡沫的相對(duì)密度關(guān)系密切。Reyes等發(fā)現(xiàn)高密度泡沫的使用能有效增加結(jié)構(gòu)總吸能，但會(huì)導(dǎo)致結(jié)構(gòu)比、吸能比相對(duì)中空薄壁管低。為進(jìn)一步提高泡沫填充管的耐撞性能和變形穩(wěn)定性，一種基于填充泡沫密度梯度進(jìn)行優(yōu)化設(shè)計(jì)的新結(jié)構(gòu)——功能梯度泡沫填充薄壁管引起了學(xué)者的研究興趣。

功能梯度泡沫填充薄壁管(Functionally Graded Foam-filled Thin-walled Column，F(xiàn)GF管)是指泡沫密度以某種預(yù)定義方式連續(xù)變化的薄壁管結(jié)構(gòu)。[21]與傳統(tǒng)的均值密度泡沫填充薄壁管(Uniform Density Foam-filled Thin-walled Column，UDF管)相比，這種新構(gòu)型耐撞性設(shè)計(jì)空間更加廣闊。Kiernan等人采用實(shí)驗(yàn)和有限元仿真相結(jié)合的方法，比較了梯度泡沫在多種密度梯度下的應(yīng)力波傳遞規(guī)律。研究表明，結(jié)構(gòu)受到的沖擊力峰值和塑性變形吸收的能量是由密度梯度決定的。目前，特定密度梯度的功能梯度泡沫材料已經(jīng)可以在實(shí)驗(yàn)室條件下獲得。

實(shí)車碰撞中，純軸向碰撞發(fā)生幾率很小，大多數(shù)情況下，碰撞方向會(huì)與薄壁管的軸線存在一定夾角。研究表明，發(fā)生斜向沖擊時(shí)，薄壁管的吸能特性隨沖擊角度的改變發(fā)生急劇變化。對(duì)特定結(jié)構(gòu)而言，存在某一臨界角，沖擊角度小于此角度時(shí)，結(jié)構(gòu)產(chǎn)生穩(wěn)定的漸進(jìn)變形；沖擊角度大于此臨界角時(shí)，結(jié)構(gòu)就會(huì)發(fā)生側(cè)向彎曲失效，導(dǎo)致能量吸收迅速下降。[10]因此，針對(duì)薄壁結(jié)構(gòu)在斜向載荷沖擊下的耐撞性研究具有重要的現(xiàn)實(shí)意義，并已得到研究人員的廣泛關(guān)注。

基于以上分析，考慮單位質(zhì)量的能量吸收、初始碰撞力峰值和載荷方向的不確定性，本文采用有限元仿真的方法，分析了功能梯度泡沫填充的方截面薄壁鋁管在斜向載荷沖擊下的變形模式和耐撞特性，并與等質(zhì)量的均值密度泡沫填充管進(jìn)行了全面比較。探討了材料參數(shù)、結(jié)構(gòu)參數(shù)和邊界條件等對(duì)其耐撞性指標(biāo)的影響，旨在為此類結(jié)構(gòu)的優(yōu)化設(shè)計(jì)和工程應(yīng)用提供理論基礎(chǔ)。

1 問題描述

1.1 物理模型

本文研究的結(jié)構(gòu)為軸對(duì)稱的正方形截面直管，內(nèi)部填充功能梯度泡沫鋁材料。如圖1所示，正方形截面邊長(zhǎng)為a，管長(zhǎng)為l，壁厚為t。碰撞過程中，薄壁管沿Z軸正方向撞擊剛性墻，沖擊速度大小為v。剛性墻法向在X-Z平面內(nèi)，與Z軸負(fù)方向夾角(加載角)為α。為模擬實(shí)車碰撞中的汽車的真實(shí)動(dòng)能，薄壁管自由端附加質(zhì)量460,kg。

1.2 斜向沖擊下結(jié)構(gòu)的耐撞性指標(biāo)

薄壁管在碰撞過程中典型的沖擊力-位移曲線如圖2所示。碰撞發(fā)生后，結(jié)構(gòu)受到的沖擊力短時(shí)間內(nèi)急劇增大然后減小，并穩(wěn)定在某一數(shù)值附近振蕩。

圖2 典型薄壁管力-位移曲線Fig.2Impact force versus displacement in a typical thinwalled column

由圖2可知，結(jié)構(gòu)吸收的總能量可表示為：

式中，d為結(jié)構(gòu)的有效變形量。據(jù)文獻(xiàn)[12]，本文中有效變形量取總長(zhǎng)的70%,。

薄壁結(jié)構(gòu)耐撞性設(shè)計(jì)通常關(guān)注兩個(gè)方面：①結(jié)構(gòu)是否具有良好的吸能效率；②沖擊載荷極值是否低于容許值。因此，本文采用比吸能(Specific Energy Absorption，SEA)和初始碰撞力峰值(Peak Crushing Force，PCF)作為評(píng)判結(jié)構(gòu)耐撞性的指標(biāo)。

比吸能即結(jié)構(gòu)單位質(zhì)量吸收的沖擊能量，據(jù)公式(1)，可表示為：

式中，M為結(jié)構(gòu)總質(zhì)量，E為結(jié)構(gòu)吸收的總能量。

考慮碰撞角度的不確定性，為全面準(zhǔn)確評(píng)估薄壁結(jié)構(gòu)在軸向和斜向碰撞下的耐撞性能，本文采用亓昌等提出的耐撞性能綜合評(píng)價(jià)準(zhǔn)則。[32]其中，斜向沖擊下結(jié)構(gòu)的比吸能定義為：

式中，SEAαi為沖擊角度為時(shí)的結(jié)構(gòu)比吸能；為此角度下的比吸能權(quán)重，與該角度下碰撞概率有關(guān)。

本文中，將初始碰撞力峰值(PCF)定義為碰撞全過程中前20%,階段內(nèi)沖擊力的最大值。考慮載荷方向不確定性的初始碰撞力峰值，定義為：

2 仿真模型

2.1 功能梯度泡沫材料

本文中，薄壁管中填有密度沿軸向變化的功能梯度泡沫材料，其密度梯度可用如下指數(shù)函數(shù)描述：

式中，ρfa和ρfb分別為薄壁管碰撞端和自由端泡沫的相對(duì)密度(相對(duì)于金屬鋁)。l代表管長(zhǎng)，y代表距碰撞端的距離，m是決定密度梯度變化的指數(shù)參數(shù)。

在參數(shù)m作用下，泡沫密度自碰撞端沿軸向逐漸變化。m＜1時(shí)，函數(shù)為凸函數(shù)；m≥1時(shí)，函數(shù)為凹函數(shù)，如圖3所示。碰撞中，薄壁管理想的變形模式是由碰撞端向自由端的漸進(jìn)式壓潰，因此設(shè)定即泡沫密度從碰撞端開始沿軸向遞增。

圖3 泡沫密度梯度隨m值變化示意圖(0表示碰撞端，1表示自由端)Fig.3 Variation in density versus normalized distance (0 representing impact end，1 representing free end)

現(xiàn)存有限元軟件中沒有能夠有效模擬泡沫密度連續(xù)變化的材料模型。據(jù)此，本文采用分層建模方法，以求實(shí)現(xiàn)準(zhǔn)確仿真。將薄壁管模型沿軸向均勻劃分為10層，每層填充各向同性的單一密度泡沫。其中，第1層(碰撞端)泡沫的相對(duì)密度設(shè)為ρf1，第10層(自由端)的相對(duì)密度為ρf10。各層厚度相同，均為泡沫芯總長(zhǎng)度的1/10。結(jié)合公式(5)得到離散化的各層泡沫材料的相對(duì)密度公式：

式中，ρfn為從碰撞端開始第n層泡沫的相對(duì)密度，m為控制密度梯度變化的指數(shù)參數(shù)。

2.2 有限元建模

FGF管有限元模型如圖4所示。結(jié)構(gòu)的建模在有限元前處理軟件HYPERMESH中完成。薄壁管外壁采用Belytschko-Tsay四節(jié)點(diǎn)薄殼單元建模，沿厚度方向取4個(gè)積分點(diǎn)，面內(nèi)采用縮減積分。泡沫材料使用8節(jié)點(diǎn)實(shí)體單元建模，以保證有限元仿真的可靠性。泡沫和管壁的單元特征長(zhǎng)度均取4,mm。采用非線性顯式有限元求解器LS-DYNA對(duì)薄壁結(jié)構(gòu)的力學(xué)響應(yīng)進(jìn)行數(shù)值仿真。采用自動(dòng)單面接觸算法(CONTACT_AUTOMATIC_SINGLE_SURFACE)模擬薄壁管與剛性墻之間，及管壁自身變形過程中可能產(chǎn)生的自接觸。采用自動(dòng)面面接觸算法(CONTACT_AUTOMATIC_SURFACE_TO_SURFACE)模擬管壁與泡沫間的接觸。采用內(nèi)部接觸算法(CONTACT_INTERIOR)模擬碰撞中泡沫內(nèi)部產(chǎn)生的自接觸。計(jì)算中考慮各部件間摩擦，摩擦系數(shù)為0.3。剛性墻用關(guān)鍵字RIGID WALL_RWPlanar定義。

圖4 功能梯度泡沫填充薄壁管有限元模型Fig.4 Finite element models

薄壁管外壁采用AA6061-T4鋁合金材料，密度ρ=2.7×10?6kg/mm3，楊氏模量E=64.8GPa ，泊松比為ν=0.28，材料屈服應(yīng)力σs=110.3MPa，強(qiáng)度極限為σu=213.0MPa，延伸率ε=19%?？紤]載荷條件，采用顯式有限元軟件LS-DYNA中第24號(hào)材料模型(MAT_PIECEWISE_LINEAR_PLASTICITY)建模。鋁合金對(duì)應(yīng)變率不敏感，故不考慮材料應(yīng)變率效應(yīng)。

考慮求解精度和計(jì)算效率之間的平衡，本文選用LS-DYNA中63號(hào)(MAT_CRUSHABLE_FOAM)材料進(jìn)行泡沫鋁的建模。其中E=66GPa，泊松比ν=0.01，拉伸應(yīng)力臨界值為1.11，阻尼系數(shù)為0.05。

2.3 有限元模型驗(yàn)證

現(xiàn)有文獻(xiàn)中沒有FGF管在斜向沖擊下力學(xué)響應(yīng)的準(zhǔn)確實(shí)驗(yàn)數(shù)據(jù)，難以進(jìn)行直接驗(yàn)證。m＝0時(shí)，可將UDF管看作FGF管的特殊形式，而針對(duì)UDF管在不同方向沖擊下動(dòng)態(tài)響應(yīng)的驗(yàn)證是切實(shí)可行的。

圖5 泡沫填充薄壁管軸向碰撞仿真值與理論值對(duì)比Fig.5 Simulation and theory results comparison regarding foam-filled columns under axial loadings

圖5為泡沫填充薄壁管受軸向沖擊時(shí)平均碰撞力的理論值[11]和仿真值的對(duì)比。由圖可知，本文有限元模型準(zhǔn)確反映了結(jié)構(gòu)在軸向沖擊下的力學(xué)響應(yīng)。

據(jù)文獻(xiàn)[26]的實(shí)驗(yàn)參數(shù)，對(duì)泡沫填充管進(jìn)行準(zhǔn)靜態(tài)斜向壓潰仿真。圖6對(duì)比了仿真值和實(shí)驗(yàn)值?？紤]碰撞力峰值和變形模式，仿真值均與試驗(yàn)值吻合良好，驗(yàn)證了模型在斜向沖擊下的正確性和有效性。

圖6 泡沫填充薄壁管準(zhǔn)靜態(tài)斜向壓潰仿真結(jié)果與實(shí)驗(yàn)對(duì)比Fig.6 Simulation and experimental results comparison regarding foam-filled columns under quasi-static oblique loadings

3 結(jié)果與分析

3.1 FGF管與UDF管斜向沖擊下耐撞性比較

圖7為FGF管和UDF管碰撞過程中的典型力學(xué)響應(yīng)比較。UDF管受到的沖擊力在碰撞初期迅速增大，后穩(wěn)定在某數(shù)值附近震蕩；作為一種新型吸能結(jié)構(gòu)，F(xiàn)GF管能準(zhǔn)確控制泡沫填充物的密度分布，其受到的沖擊力隨碰撞過程逐漸增大?？梢钥闯?，泡沫鋁材料的密度梯度對(duì)薄壁管吸能特性和變形模式影響很大?？紤]載荷方向的不確定性，比較不同碰撞角度下兩結(jié)構(gòu)的耐撞性能。

圖7 FGF管和UDF管在碰撞過程中的典型力學(xué)響應(yīng)比較Fig.7Typical mechanical response of FGF and UDF columns

參照實(shí)車前縱梁尺寸，本文基準(zhǔn)模型的幾何參數(shù)如下：a＝80,mm；l＝240,mm；t＝1.5,mm。作為常量，F(xiàn)GF管兩端的相對(duì)密度(相對(duì)于金屬鋁)分別為：ρf1=0.08(碰撞端)和ρf10=0.20(自由端)。薄壁管沖擊剛性墻的初速度為v=14m/s(50,km/h)。

圖8 斜向沖擊下FGF管和UDF管比吸能對(duì)比Fig.8 SEA of FGF and UDF columns under oblique impact

圖8比較了斜向沖擊下兩結(jié)構(gòu)的比吸能大小?；鶞?zhǔn)模型基礎(chǔ)上，取m＝1(關(guān)于壁厚和泡沫密度范圍等參數(shù)的討論中，均取m＝1)?？汕宄乜闯?，兩種結(jié)構(gòu)的SEA均隨沖擊角度的增大單調(diào)減小，小角度時(shí)，兩結(jié)構(gòu)SEA相差不大；大角度時(shí)，F(xiàn)GF管優(yōu)勢(shì)明顯。總體來看，F(xiàn)GF管的SEA始終高于UDF管，并在α=25°時(shí)優(yōu)勢(shì)最大。與泡沫密度均勻分布的UDF管相比，F(xiàn)GF管填充泡沫的密度(剛度)自碰撞端沿軸向遞增的特點(diǎn)，更容易促使結(jié)構(gòu)產(chǎn)生對(duì)稱塑性鉸，發(fā)生穩(wěn)定高效的漸進(jìn)壓潰。這種泡沫密度分布方式在大角度沖擊下的優(yōu)勢(shì)尤其明顯。

另外，兩種結(jié)構(gòu)吸能特性的差異與薄壁管的失效模式有關(guān)。臨界加載角是指結(jié)構(gòu)的失效模式隨沖擊角度的增大發(fā)生改變時(shí)的角度，是評(píng)判薄壁管在斜向沖擊下吸能性能優(yōu)劣的重要依據(jù)。沖擊角度小于臨界值時(shí)，結(jié)構(gòu)發(fā)生高效漸進(jìn)壓潰；沖擊角度大于臨界值時(shí)，結(jié)構(gòu)發(fā)生吸能很少的側(cè)向彎曲失效。因此，臨界加載角越大，結(jié)構(gòu)吸能穩(wěn)定性越好，斜向沖擊下耐撞性越好。由圖8可知，UDF管的臨界加載角介于20～25,°之間，F(xiàn)GF管的臨界加載角介于25～30,°之間。

圖9 α＝25 °時(shí)FGF管和UDF管變形模式對(duì)比Fig.9 Failure modes of FGF and UDF columns atα＝25 °

圖9為α=25,°時(shí)，兩種構(gòu)型的變形模式對(duì)比。可以看出，UDF管已明顯側(cè)向彎曲失效，而FGF管仍能較好保持漸進(jìn)壓潰變形，故此角度下吸能優(yōu)勢(shì)最明顯。

考慮載荷方向的不確定性，斜向沖擊下兩種結(jié)構(gòu)的比吸能SEAα(設(shè)各角度下比吸能所占比重均勻分布，即權(quán)重系數(shù)均相同，45,°)分別為：SEAα(FGF)=9.9 kJ/kg ，SEAα(UDF)= 7.9 kJ/kg。與等質(zhì)量的UDF管相比，F(xiàn)GF管在斜向沖擊下比吸能增大25.3%,，吸能效率提升顯著。

圖10 斜向沖擊下FGF管和UDF管初始碰撞力峰值對(duì)比Fig.10 PCF of FGF and UDF columns under oblique impact

圖10比較了兩種結(jié)構(gòu)斜向沖擊下的初始碰撞力峰值。可以看出，隨碰撞角度的增大，兩種結(jié)構(gòu)的初始碰撞力峰值均單調(diào)減小，F(xiàn)GF管的PCF始終小于等質(zhì)量UDF管，且在碰撞力較大的小角度載荷條件下優(yōu)勢(shì)最為明顯。

綜上所述，斜向沖擊下，考慮小角度沖擊下的吸能效率和大角度沖擊下的穩(wěn)定性，F(xiàn)GF管的綜合耐撞性能明顯優(yōu)于等質(zhì)量的UDF管。

3.2 結(jié)構(gòu)參數(shù)對(duì)耐撞性的影響分析

泡沫填充薄壁管的吸能主要由3方面組成，即管壁吸能、泡沫吸能和管壁與泡沫相互作用產(chǎn)生的交互能。因此，本文將密度梯度指數(shù)m、結(jié)構(gòu)自由端泡沫相對(duì)密度ρf10和壁厚t作為3個(gè)結(jié)構(gòu)參數(shù)。各參數(shù)取值如表1所示。

表1 功能梯度泡沫填充薄壁管結(jié)構(gòu)參數(shù)取值Tab.1 Structural parameters of FGF columns

3.2.1 泡沫密度梯度的影響

圖11 斜向沖擊下梯度參數(shù)對(duì)比吸能的影響Fig.11 Effect of parameter m on SEA under oblique impact

圖11為斜向沖擊下梯度參數(shù)m對(duì)結(jié)構(gòu)SEA的影響?？梢钥闯觯S著沖擊角度的增大，不同密度梯度薄壁管SEA均單調(diào)減小，但結(jié)構(gòu)的比吸能并不與梯度參數(shù)m之間單調(diào)相關(guān)。α＜15°時(shí)，各管SEA值隨m值的增大先增后減，在m＝0.2時(shí)取到最大值；15°≤α≤30°時(shí)，泡沫平均密度(剛度)較大的管子(m＝0.1、m＝0.2)的SEA迅速減小，其余各管變化較平緩；α＞30°時(shí)，各結(jié)構(gòu)均已發(fā)生側(cè)向彎曲失效，比吸能變化趨勢(shì)類似，差距逐漸縮小。

圖12 斜向沖擊下薄壁管比吸能對(duì)比Fig.12 SEAαof different columns under oblique impact

小角度時(shí)，各管均能保持高效的漸進(jìn)式變形模式，產(chǎn)生塑性鉸的個(gè)數(shù)相近，平均密度較大的管子(如m＝0.1、m＝0.2)平均剛度也較大，因此吸能更多；大角度時(shí)，剛度越大的結(jié)構(gòu)越容易發(fā)生側(cè)向彎曲失效，無法產(chǎn)生足夠多的塑性鉸，因此吸能表現(xiàn)反而不如剛度小的結(jié)構(gòu)。

圖12為不同m值的薄壁管在斜向沖擊下的比吸能(SEAα)比較?？梢悦黠@看出，m≤1時(shí)，薄壁管在斜向沖擊下的吸能效率較高，并在m＝0.5時(shí)達(dá)到最大。

綜上所述，m值的大小對(duì)薄壁管耐撞性能影響很大。綜合考慮小角度沖擊下的吸能效率和大角度沖擊下的耐撞穩(wěn)定性，m＝0.5的薄壁管斜向沖擊下耐撞性能最佳。

圖13 斜向碰撞下梯度參數(shù)對(duì)初始碰撞力峰值的影響Fig.13 Effect of parameter m on PCF under oblique impact

圖13為斜向沖擊下參數(shù)m對(duì)PCF的影響。可以看出，隨著沖擊角度的增大，各種結(jié)構(gòu)的PCF均單調(diào)減小。本組對(duì)比中，各管碰撞端泡沫密度均相同，沖擊角度一定時(shí)，PCF差距不明顯。

3.2.2 泡沫密度變化范圍的影響

圖14 斜向碰撞下泡沫密度范圍對(duì)比吸能的影響Fig.14 Effect of the range of foam density on SEA under oblique impact

圖14為斜向碰撞下泡沫密度范圍對(duì)比吸能的影響。可以看出，各管SEA均隨沖擊角度的增大單調(diào)減小，變化趨勢(shì)類似。沖擊角度一定時(shí)，密度變化范圍更大的結(jié)構(gòu)，SEA越大。碰撞端泡沫密度固定時(shí)，密度變化范圍大的管子平均剛度更大，能在小角度碰撞時(shí)保持較高的吸能效率；大角度下，密度變化范圍大的管子，吸能泡沫的剛度變化區(qū)間更廣，更有利于誘導(dǎo)漸進(jìn)壓潰變形的發(fā)生。但須注意的是，泡沫鋁相對(duì)密度過大，亦會(huì)引起結(jié)構(gòu)比吸能的下降。綜上所述，一定范圍內(nèi)增大填充泡沫密度的變化范圍，能有效提升結(jié)構(gòu)的耐撞性能。

圖15為斜向碰撞下泡沫密度范圍對(duì)初始碰撞力峰值的影響。可以看出，圖15與圖13變化趨勢(shì)非常類似，說明碰撞端泡沫密度是影響FGF管初始碰撞力峰值的重要因素，而泡沫的梯度和變化范圍對(duì)結(jié)構(gòu)的PCF影響不大。

圖15 斜向碰撞下泡沫密度范圍對(duì)初始碰撞力峰值的影響Fig.15Effect of the range of foam density on PCF under oblique impact

3.2.3 壁厚的影響

圖16為斜向沖擊下壁厚t對(duì)SEA的影響。各管SEA均隨沖擊角度的增大而減小。同一角度下，各管的SEA隨壁厚的增加單調(diào)增大。沖擊角度為25 °(臨界加載角附近)時(shí)，壁厚不同的各結(jié)構(gòu)比吸能差距最小。這是因?yàn)?，管壁厚度與薄壁結(jié)構(gòu)剛度成正比，在臨界加載角附近時(shí)，剛度小的結(jié)構(gòu)更易發(fā)生局部屈服，促使剩余結(jié)構(gòu)產(chǎn)生穩(wěn)定的逐層壓潰變形。因此，一定范圍內(nèi)增大壁厚，可有效提高結(jié)構(gòu)的比吸能，但在大角度沖擊下的耐撞穩(wěn)定性有所下降。

圖16 斜向碰撞壁厚對(duì)比吸能的影響Fig.16Effect of the wall thickness on SEA under oblique impact

圖17給出了壁厚對(duì)PCF的影響。可以看出，同一沖擊角度下，碰撞力峰值隨壁厚單調(diào)增大。壁厚的增加使管子的剛度和碰撞接觸面面積均增大。因此，減小壁厚可有效降低結(jié)構(gòu)的初始碰撞力峰值。

圖17 斜向碰撞下壁厚對(duì)初始碰撞力峰值的影響Fig.17 Effect of the wall thickness on PCF under oblique impact

4 結(jié) 論

本文研究了一種功能梯度泡沫填充管在斜向碰撞下的失效模式和耐撞性能，分析了包括泡沫梯度分布指數(shù)、泡沫密度變化范圍和壁厚在內(nèi)的結(jié)構(gòu)參數(shù)對(duì)其在不同沖擊角度下耐撞特性的影響，得到了一系列指導(dǎo)工程實(shí)踐的有用結(jié)論。研究結(jié)果表明：①與等質(zhì)量的UDF管相比，F(xiàn)GF管在斜向沖擊載荷作用下SEA更高，PCF更小，小角度沖擊時(shí)吸能效率更高，大角度沖擊時(shí)吸能穩(wěn)定性更好，綜合耐撞性能優(yōu)勢(shì)明顯；②一定范圍內(nèi)，調(diào)整指數(shù)參數(shù)m使其接近m＝0.5、增大泡沫密度變化范圍以及減小壁厚，可有效提高FGF管在斜向沖擊下的臨界加載角，提高結(jié)構(gòu)在斜向碰撞下的耐撞穩(wěn)定性；③影響功能梯度泡沫填充管PCF的兩大因素為管壁的厚度和碰撞端泡沫密度。特定沖擊角度下，壁厚越薄，碰撞端泡沫密度越小，結(jié)構(gòu)的初始碰撞力峰值就越小。梯度泡沫的分布方式和變化范圍對(duì)結(jié)構(gòu)的PCF影響不大。需要指出的是，各種結(jié)構(gòu)參數(shù)的變化都會(huì)引起功能梯度泡沫填充薄壁管在斜向沖擊下耐撞性能的改變，且相互影響，彼此相關(guān)，有必要進(jìn)行進(jìn)一步優(yōu)化設(shè)計(jì)，以找到符合要求的最優(yōu)結(jié)構(gòu)?！?/p>

［1］ Mamalis A G，Manolakos D E，Ioannidis M B，et al. Finite element simulation of the axial collapse of metallic thin-walled tubes with octagonal cross-section [J]. Thin-Walled Structures，2003(41)：891-900.

［2］ 馬琳琳，楊娜，趙桂范. 轎車結(jié)構(gòu)耐撞性分析與改進(jìn)[J]. 機(jī)械設(shè)計(jì)與制造，2008(6)：106-108.

［3］ Nagel G M，Thambiratnam D P. A numerical study on the impact response and energy absorption of tapered thin-walled tubes [J]. International Journal of Mechanical Sciences，2004(46)：201-216.

［4］ Nagel G M，Thambiratnam D P. Computer simulation and energy absorption of tapered thin-walled rectangular tubes [J]. Thin-Walled Structures，2005(43)：1225-1242.

［5］ Zhang X，Cheng G，Zhang H. Theoretical prediction and numerical simulation of multi-cell square thin-walled structures [J]. Thin-Walled Structures，2006(44)：1185-1191.

［6］ Reyes A，Hopperstad O S，Berstad T，et al. Constitutive modeling of aluminum foam including fracture and statistical variation of density [J]. European Journal of Mechanics-A/Solids，2003(22)：815-835.

［7］ Hanssen A G，Langseth M，Hopperstad O S. Static and dynamic crushing of square aluminium extrusions with aluminium foam filler [J]. International Journal of Impact Engineering，2000(24)：347-383.

［8］ Chen W G，Wierzbicki T. Relative merits of singlecell，multi-cell and foam-filled thin-walled structures in energy [J]. Thin-Walled Structures，2001，39(4)：287-306.

［9］ Hou S，Li Q，Long S，et al. Crashworthiness design for foam filled thin-wall structures [J]. Materials & Design，2009(30)：2024-2032.

［10］ 亓昌，董方亮，楊姝. 錐形多胞薄壁管斜向沖擊吸能特性仿真研究[J]. 振動(dòng)與沖擊，2012，31(24)：102-106.

Crashworthiness Analysis for Functionally Graded Foam-filled Thin-walled Column Under Oblique Impact

ZHAO Bowen，ZHONG Yanhui，ZHUANG Mengmeng
(China Automotive Technology & Research Center，Automobile Testing Research Institute，Tianjin 300300，China)

The failure modes and crashworthiness of a new functionally graded foam-filled thin-walled columnunder oblique impact were investigated using finite element analysis(FEA)．The specific energy absorption(SEA)and the initial peak force were adopted as performance indices．The effects of structural parameters including index parameter of graded foam，wall thickness and range of foam density on the crashworthiness characteristics of the columns were analyzed by comparing with the uniform foam-filled counterparts．It is found that the functionally graded foam-filled thin-walled columns are superior to their uniform counterparts considering energy-absorbing characteristics under oblique impact．It should be noted that structural parameters influence with each other，and all of them significantly influence the crashworthiness characteristics of the structures.

thin-walled column；functionally graded foam；energy absorption；crashworthiness；oblique impact

O313.4；U463.8

A

1006-8945(2016)03-0068-07

中國(guó)汽車技術(shù)研究中心青年創(chuàng)新基金。

2016-01-26