聚氨酯泡沫填充的碳纖維增強(qiáng)復(fù)合材料錐管吸能性能數(shù)值模擬及試驗(yàn)驗(yàn)證

干年妃 王多華 馮亞楠 張 謙

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

0 引言

硬質(zhì)聚氨酯泡沫(PU泡沫)是密度很小、價格相對低廉、容易合成的多孔介質(zhì)，在絕熱、防振、包裝等場合應(yīng)用比較廣泛。由于PU泡沫具有良好的吸收動能的特性，能夠緩和沖擊、減弱振蕩、減小應(yīng)力幅值，因此目前不少學(xué)者都在研究它。林玉亮等[1]、曹靜等[2]研究了硬質(zhì)聚氨酯泡沫的壓縮性能；PAULINO等[3]將PU泡沫填充在車門中，提高了汽車側(cè)碰安全性；劉文濟(jì)等[4]綜述了車用聚氨酯材料的應(yīng)用進(jìn)展；AVALLE等[5]研究了聚合物結(jié)構(gòu)泡沫在沖擊載荷下的吸能能力。碳纖維增強(qiáng)復(fù)合材料(carbon fiber reinforced polymer，CFRP)作為一種低密度、高比強(qiáng)度、高韌性、高比模量、耐腐蝕、高抗沖擊性能的功能結(jié)構(gòu)材料，已在航空、航天、武器、船舶、汽車[6-7]等領(lǐng)域得到廣泛的研究和應(yīng)用。王榮惠等[8]研究了碳纖維增強(qiáng)復(fù)合材料壓縮力學(xué)性能；蘇尚彬等[9]研究了碳纖維材料保險杠耐撞性及輕量化；JACOB等[10]研究了聚合物復(fù)合材料在汽車碰撞吸能中的應(yīng)用；MAMALIS等[11]研究得出玻璃纖維復(fù)合材料錐形殼結(jié)構(gòu)準(zhǔn)靜態(tài)壓縮過程中的破壞模式主要為漸進(jìn)壓潰和中部失效；BORIA等[12]將錐形碳纖維吸能器用于電動車上，驗(yàn)證了碳纖維吸能器相比金屬吸能器有明顯的優(yōu)勢；李志斌等[13]與張勇等[14]研究了薄壁管及其泡沫金屬填充結(jié)構(gòu)的耐撞性；宋宏偉等[15]研究了多孔材料填充薄壁結(jié)構(gòu)吸能的相互作用效應(yīng)；高偉釗等[16]與PAZ等[17]為得到最優(yōu)的吸能效果，對泡沫填充吸能結(jié)構(gòu)進(jìn)行了優(yōu)化。本文研究了PU泡沫和碳纖維增強(qiáng)復(fù)合材料在準(zhǔn)靜態(tài)壓縮作用下的吸能能力，并通過試驗(yàn)仿真驗(yàn)證了PU泡沫填充的碳纖維錐管的吸能性能。

1 評價指標(biāo)

評價某一結(jié)構(gòu)吸能能力的指標(biāo)是多樣的，主要有峰值載荷、平均載荷、比吸能等[18]。在評價動態(tài)性的碰撞問題中，峰值力、峰值加速度、平均載荷等指標(biāo)都是不可缺少的。由于本次的試驗(yàn)條件是準(zhǔn)靜態(tài)壓縮，所以峰值力只做參考，不做重點(diǎn)對比，主要研究結(jié)構(gòu)的總吸能(energy absorbed，EA)和比吸能(specific absorbed energy，SAE)。

總吸能δEA由下式計算:

(1)

式中,x為壓縮行程;f(x)為吸能結(jié)構(gòu)在壓縮過程中所受到的力;d為總行程。

比吸能δSEA即單位質(zhì)量所吸收的能量，是評價吸能性好壞的重要指標(biāo)，其計算公式為

(2)

式中，m為結(jié)構(gòu)的總質(zhì)量。

δSEA的值越大，結(jié)構(gòu)的吸能性能越好。

2 試驗(yàn)儀器介紹

準(zhǔn)靜態(tài)試驗(yàn)是在英斯特朗R8505和MTS-647萬能材料試驗(yàn)機(jī)(圖1)上完成的，試驗(yàn)過程中可以直接記錄對試樣施加的力和壓縮位移，MTS647萬能材料試驗(yàn)機(jī)的載荷范圍為0～500 kN，英斯特朗R8505的載荷范圍為0～150 kN。所有試樣放置在試驗(yàn)平臺正中心(圖1a中1位置)，試樣上端與加載裝置(圖1a中2位置)之間放置一厚圓盤，以保證試樣上端受力均勻，試樣兩端不固定。由于試驗(yàn)臺上有同心圓凹環(huán)，且試樣表面粗糙，所以試驗(yàn)過程中沒有發(fā)生試樣相對于試驗(yàn)平臺的相對位移，加載速度為2 mm/min。

(a)MTS 647 (b)英斯特朗R8505圖1 萬能材料試驗(yàn)機(jī)Fig.1 Universal material testing machine

3 材料模型的確定

LS-DYNA有一個非常大的材料庫供用戶選擇，同種類型的材料可能有多種材料模型與其對應(yīng)，仿真前需確定最合適本文材料的材料模型。

3.1 PU泡沫數(shù)值模擬模型的確定

本次試驗(yàn)所用硬質(zhì)PU泡沫的密度為0.2 g/m3，根據(jù)GB/T 8813-2008進(jìn)行硬質(zhì)PU泡沫標(biāo)定試驗(yàn)獲取參數(shù)。

泡沫材料仿真建模一直是研究的難點(diǎn)，目前在LS-DYNA中泡沫材料模型有數(shù)十種，根據(jù)相關(guān)文獻(xiàn)中各種泡沫模型的力學(xué)特性敘述，以及對PU泡沫試驗(yàn)結(jié)果曲線的觀察，擬采用MAT163號材料模型作為本文中PU泡沫的仿真數(shù)學(xué)模型[19]。MAT163號模型在MAT63號模型的基礎(chǔ)上，增加了可選的黏性系數(shù)、拉伸截止應(yīng)力、應(yīng)變速率的影響。模型視為完全彈性卸載，達(dá)到拉伸截止應(yīng)力時，拉伸視為完全彈塑性。該模型需要輸入的參數(shù)比MAT63號模型多了黏性系數(shù)、平均體積應(yīng)變的循環(huán)次數(shù)、拉伸截止應(yīng)力及不同應(yīng)變速率下應(yīng)力和體積應(yīng)變關(guān)系曲線。MAT163材料模型的輸入?yún)?shù)如表1所示。

表1 MAT163材料模型的參數(shù)Tab.1 Performance parameters of MAT163material model

3.2 CFRP材料數(shù)值模擬模型的確定

為了獲得碳纖維管的材料參數(shù)，團(tuán)隊成員對CFRP材料進(jìn)行了拉伸試驗(yàn)。拉伸試樣是由與制造薄壁外皮結(jié)構(gòu)材料一致的碳纖維預(yù)侵布制造而成的，所用纖維布為T300,正交平紋布。試驗(yàn)所使用的設(shè)備為英斯特朗R8505萬能材料試驗(yàn)機(jī)，根據(jù)GB 1447-2005開展試驗(yàn)，進(jìn)行了0°拉伸試驗(yàn)。根據(jù)GB/T 3355-2005，進(jìn)行了45°拉伸試驗(yàn)。通過0°拉伸試驗(yàn)可以獲取材料的泊松比、抗拉強(qiáng)度以及彈性模量。通過45°拉伸試驗(yàn)可以得到材料的剪切模量與剪切強(qiáng)度。根據(jù)GB/T 5258-1995，進(jìn)行了壓縮性能試驗(yàn)，獲得了材料的抗壓強(qiáng)度。

LS-DYNA中復(fù)合材料模型有很多，常用的有MAT54、MAT58、MAT59等，經(jīng)過多次驗(yàn)證，選擇MAT59號材料模型作為碳纖維錐管仿真材料模型最為準(zhǔn)確。MAT59號模型是帶失效的復(fù)合材料模型[19]，適用于體單元?dú)卧疚闹惺褂玫腃FRP材料厚度在2 mm內(nèi)，因此使用殼單元建模，在Hypermesh中鋪層，接觸方式為單面接觸，因此只用輸入兩個主方向的彈性模量、泊松比、剪切模量及密度等。MAT59號材料模型輸入?yún)?shù)如表2所示。

表2 CFRP材料參數(shù)Tab.2 Performance parameters of CFRP material model

4 圓形截面錐形結(jié)構(gòu)準(zhǔn)靜態(tài)壓縮仿真結(jié)果與試驗(yàn)對比分析

4.1 試樣尺寸

為更好地符合汽車吸能的真實(shí)情況，減小壓縮時初始峰值力，有研究表明圓形截面結(jié)構(gòu)吸能性較好[20]，因此本文研究結(jié)構(gòu)為圓形截面錐形結(jié)構(gòu)，樣件如圖2所示，試驗(yàn)所用樣件均用車床加工制作。試驗(yàn)、仿真模型尺寸完全一樣，仿真建模在CATIA中進(jìn)行，使用Hypermesh前處理，在LS-DYNA中計算，為表達(dá)方便，后文中C(circular)表示圓形截面，PU-C表示圓形截面PU錐臺，CFRP-C表示CFRP材料錐形殼； T(test)表示試驗(yàn)結(jié)果，S(simulation)表示仿真結(jié)果；CFRP-C-S表示圓形截面CFRP錐臺仿真所得數(shù)據(jù);CFF-PU-S表示在仿真試驗(yàn)時，填充結(jié)構(gòu)中PU泡沫吸收的能量；CFF-CF-S表示在仿真試驗(yàn)時，填充結(jié)構(gòu)中CFRP外殼吸收的能量；CFF-C-S表示在仿真試驗(yàn)時，填充結(jié)構(gòu)吸收的能量。

(a)截面圖 (b)立體圖圖2 樣件尺寸Fig.2 Size of sample

4.2 圓形截面PU泡沫塑料試驗(yàn)仿真對比分析

(a)仿真

(b)試驗(yàn)圖3 PU泡沫塑料仿真與試驗(yàn)變形模式比較Fig.3 Comparison of deformation model of PU foam plastics between simulation and experimental

PU泡沫塑料仿真與試驗(yàn)變形模式比較如圖3所示。由于泡沫材料在壓縮時會出現(xiàn)大變形而引起過大的沙漏能，且易出現(xiàn)負(fù)體積導(dǎo)致計算錯誤，所以本仿真模型定義了沙漏模式為沙漏控制類型4，沙漏系數(shù)為0.1，單元內(nèi)部接觸類型為CONTTACT_INTERIER。

PU泡沫仿真與試驗(yàn)力-位移曲線及吸能曲線對比如圖4所示。

(a)力-位移曲線

(b)吸能曲線圖4 PU泡沫仿真與試驗(yàn)力-位移曲線及吸能曲線對比Fig.4 Comparison of foam cone platform energy absorption curve between experimental and simulated

對比分析圖3和圖4可知：變形模式幾乎完全一樣，試驗(yàn)和仿真中模型都是從頂部慢慢向下潰縮，力-位移曲線基本吻合，由表3可知，壓縮70 mm總吸能誤差為4.128%。

4.3 圓形截面CFRP材料錐形殼仿真與試驗(yàn)對比分析

圓形截面CFRP材料錐形殼尺寸與截面和上述PU泡沫錐臺對應(yīng)，用熱固性環(huán)氧樹脂和國產(chǎn)3K平紋碳纖維布，先將碳纖維布剪成所需要的大小和形狀，然后用樹脂預(yù)侵，最后將碳纖維布纏繞到之前加工好的模具上，每個模具上纏8層，等待樹脂干后，通過脫模處理就獲取了碳纖維錐形管。為使接縫均勻分布，纏繞一層，把模型旋轉(zhuǎn)90°再纏繞下一層。試件在準(zhǔn)靜態(tài)壓縮過程中受力如圖5所示。由圖5可知，錐形殼上端有向內(nèi)折疊的趨勢，下端有向外折疊的趨勢。仿真與試驗(yàn)對比結(jié)果如圖6和表4所示。

表3 泡沫錐臺試驗(yàn)與仿真吸能對比Tab.3 Comparison of foam cone platform energyabsorption between experimental and simulated

圖5 圓形截面CFRP材料錐形殼受力圖Fig.5 Force analysis of circular cross section CFRP composite conical shell

(a)仿真

(b)試驗(yàn)圖6 圓形截面CFRP仿真與試驗(yàn)變形模式對比Fig.6 Comparison of deformation model of CFRPbetween simulation and experimental

總吸能(kJ)CFRP?C?T3．323CFRP?C?S4．299相對誤差(%)29．4

對比分析圓形截面CFRP材料準(zhǔn)靜態(tài)壓縮變形模式(圖6)、力-位移曲線(圖7a)、吸能曲線(圖7b)可知：在壓縮前段時間，試驗(yàn)、仿真變形模式基本一樣，都是從頂部往下部潰縮，試驗(yàn)中頂部碳纖維復(fù)合材料都是往內(nèi)折疊、基本無往外翻現(xiàn)象，纖維能否往內(nèi)外折疊主要由圓錐殼的半錐頂角決定[21-22]，如圖8所示。由圖6變形模式可知，本文試驗(yàn)是圖8中形式三的變形模式，隨著壓縮的進(jìn)行，上端復(fù)合材料全部向內(nèi)側(cè)彎曲碎裂失效。

(a)力-位移曲線

(b)吸能曲線圖7 圓形截面CFRP錐管準(zhǔn)靜態(tài)壓縮仿真與試驗(yàn)力-位移曲線及吸能曲線對比Fig.7 Comparison of CFRP energy absorption curve between experimental and simulated

(a)形式一 (b)形式二 (c)形式三圖8 直管及圓錐管復(fù)合材料結(jié)構(gòu)準(zhǔn)靜態(tài)壓縮不同破壞形式[21]Fig.8 Different failure forms of quasi static compressionof the composite structure of straight pipe and conical pipe[21]

對于復(fù)合材料管在壓縮過程中的吸能方式，GUPTA等[23]和BORIA等[24]都進(jìn)行了分析。能量耗散途徑包括：基體剪切變形能、層束彎曲能、纖維拉伸變形能、碎片與壓縮載荷接觸面的摩擦能。

本文試驗(yàn)過程中，在前段穩(wěn)定壓潰階段，出現(xiàn)的是圖8c的情況，符合BORIA等[24]所描述的過程，當(dāng)壓縮位移達(dá)到40 mm左右時，CFRP錐管失穩(wěn)，靠近底部20 mm的地方出現(xiàn)橫向裂紋，同時載荷在此刻開始減小，吸能能力大幅降低，由于部分纖維之間的相互作用，壓縮力-位移曲線緩慢下降。造成前半段穩(wěn)定，后半段失穩(wěn)的原因，可能是試件內(nèi)部有缺陷，承受一段時間載荷后，試件出現(xiàn)斷裂的現(xiàn)象，缺陷是碳纖維的致命傷[25]。仿真過程中，在壓縮位移前70 mm，結(jié)構(gòu)都比較穩(wěn)定，由上至下層層壓潰，到壓縮位移接近80 mm的時候，下部出現(xiàn)翹曲，由于向內(nèi)彎曲的材料堆積，錐管被壓實(shí)，此時的力-位移曲線迅速上升。試驗(yàn)過程中，當(dāng)壓縮位移到達(dá)80 mm附近時，錐管被壓實(shí)，力-位移曲線也是迅速上升。從力-位移曲線和吸能曲線看，在穩(wěn)定變形階段，試驗(yàn)結(jié)果與仿真結(jié)果吻合比較好，說明仿真可以用來預(yù)測穩(wěn)定狀態(tài)下復(fù)合材料錐管的吸能情況。

4.4 圓形截面PU泡沫填充的CFRP結(jié)構(gòu)仿真與試驗(yàn)對比分析

試驗(yàn)與仿真模型尺寸完全一樣，在制作填充管時，PU泡沫既可以做填料也可以充當(dāng)模具材料，可以直接把碳纖維布纏繞到PU泡沫塑料上，樹脂干了以后就得到了碳纖維填充管樣件，建模與仿真方式同上文一樣，試驗(yàn)與仿真變形模式及結(jié)果如圖9、圖10、表5所示。

(a)仿真

(b)試驗(yàn)圖9 圓形截面填充結(jié)構(gòu)仿真與試驗(yàn)變形模式對比Fig.9 Comparison of deformation model of thecomposite filled structure between simulation and experimental

比較圖9所示的圓形截面填充結(jié)構(gòu)變形模式可以看到，試驗(yàn)中填充結(jié)構(gòu)從頂部開始往下潰縮，在向下潰縮過程中由于填充了PU泡沫，外層復(fù)合材料殼不能向內(nèi)彎折，所以碎裂后的復(fù)合材料結(jié)構(gòu)外溢。隨著壓縮的進(jìn)行，上部結(jié)構(gòu)慢慢被壓實(shí)，壓縮力增大，結(jié)構(gòu)從下部開始潰縮，下部壓實(shí)后，又開始從上部潰縮，交替進(jìn)行，所以得到的力-位移曲線是上下波動的折線。仿真模擬中同樣出現(xiàn)了這種現(xiàn)象，壓縮填充結(jié)構(gòu)未出現(xiàn)像壓縮空殼結(jié)構(gòu)斷裂導(dǎo)致力-位移曲線突然下降的現(xiàn)象，說明填充結(jié)構(gòu)吸能穩(wěn)定性比空殼結(jié)構(gòu)吸能穩(wěn)定性要好。由表6與圖11可以看出，對于CFRP外殼,填充后有效吸能行程由80 mm降為70 mm，但是填充后CFRP外殼與PU泡沫的總吸能都提高了，這說明填充結(jié)構(gòu)具有很好的吸能性能，這是因?yàn)樘畛浣Y(jié)構(gòu)中，外壁與內(nèi)部PU泡沫之間的相互作用也消耗了一部分能量。由變形模式(圖9)與力-位移曲線和吸能曲線(圖10)可知，試驗(yàn)與仿真結(jié)果吻合得很好。

(a)力-位移曲線

(b)吸能曲線圖10 圓形截面填充結(jié)構(gòu)準(zhǔn)靜態(tài)壓縮仿真與試驗(yàn)力-位移曲線及吸能曲線對比Fig.10 Comparison of the composite filled structure energy absorption curve between experimental and simulated

總吸能(kJ)CFF?C?T6．897CFF?C?S6．829相對誤差(%)0．986

5 結(jié)論

(1)通過材料性能試驗(yàn)和復(fù)合材料準(zhǔn)靜態(tài)壓潰試驗(yàn)，獲得了復(fù)合材料的力學(xué)性能參數(shù)及其壓潰吸能數(shù)據(jù)。仿真與試驗(yàn)在能量吸收方面誤差比較小，表明仿真模型可以很好地模擬準(zhǔn)靜態(tài)壓潰過程，試驗(yàn)與仿真相互驗(yàn)證，有助于縮短結(jié)構(gòu)試驗(yàn)準(zhǔn)備周期、減少試驗(yàn)次數(shù)、降低試驗(yàn)費(fèi)用，而且還可以幫助試驗(yàn)人員確定優(yōu)化試驗(yàn)方案和試驗(yàn)流程。

表6 填充前與填充后各種結(jié)構(gòu)吸能變化

(a)填充PU泡沫與未填充PU泡沫吸能對比

(b)填充CFRP與未填充CFRP吸能對比圖11 填充前與填充后同一結(jié)構(gòu)仿真吸能對比Tab.11 Simulation of energy absorption curve of the same structure before filling and after filling

(2)試驗(yàn)和仿真結(jié)果均表明PU泡沫填充的碳纖維錐形結(jié)構(gòu)有較好的吸能性能，雖然填充結(jié)構(gòu)的有效吸能行程縮短了，但是吸能能力提升了，填充結(jié)構(gòu)比單獨(dú)使用一種結(jié)構(gòu)吸能性好，填充結(jié)構(gòu)吸收的能量大于兩種結(jié)構(gòu)單獨(dú)使用時吸收能量之和，同時從壓縮力-位移曲線可知，填充結(jié)構(gòu)具有很好的穩(wěn)定性，將填充結(jié)構(gòu)用于汽車前部吸能結(jié)構(gòu)是一個可選方案。

參考文獻(xiàn)：

[1] 林玉亮,盧芳云,王曉燕,等. 低密度聚氨酯泡沫壓縮行為實(shí)驗(yàn)研究[J]. 高壓物理學(xué)報,2006,20(1):88-92.

SUN Yuliang, LU Fangyun, WANG Xiaoyan, et al. Experimental Study on Compressive Behavior of Low Density Polyurethane Foam [J]. Journal of High Pressure Physics,2006,20(1):88-92.

[2] 曹靜,閆澍旺,馮現(xiàn)洪,等. 聚氨酯泡沫材料的壓縮試驗(yàn)研究[J]. 中國造船2005,46(b11):232-238.

CAO Jing, YAN Shuwang, FENG Xianhong, et al. Compression Test of Polyurethane Foam[J]. Shipbuilding of China,2005,46(b11):232-238.

[3] PAULINO M, TEIXEIRA-DIAS F. On the Use of Polyurethane Foam Paddings to Improve Passive Safety in Crashworthiness Applications[M]// Polyurethane. London:INTECH,2012:337-354.

[4] 劉文濟(jì),涂偉萍,王鋒,等. 汽車用聚氨酯材料的應(yīng)用與研究進(jìn)展[J]. 聚氨酯工業(yè),2015(3):26-30.

LIU Wenji, TU Weiping, WANG Feng, et al. Application and Research Progress of Polyurethane Material for Automobile[J]. Polyurethane Industry,2015(3):26-30.

[5] AVALLE M, BELINGARDI G, MONTANINI R. Characterization of Polymeric Structural Foams under Compressive Impact Loading by Means of Energy-absorption Diagram[J]. International Journal of Impact Engineering,2001,25(5):455-472.

[6] 張定金,陳虹,于今. 汽車輕量化對碳纖維的需求[J]. 中國石油和化工經(jīng)濟(jì)分析,2014(10):30-33.

CHEN Dingjin, CHEN Hong, YU Jin. Automobile Lightweight Needs Carbon Fiber[J]. Economic Analysis of Petroleum and Chemical Industry of China,2014(10):30-33.

[7] 馮瑞華. 碳纖維在汽車領(lǐng)域的應(yīng)用[J]. 新材料產(chǎn)業(yè),2014(8):15-18.

FENG Ruihua. Application of Carbon Fiber in Automobile Industry[J]. New Material Industry,2014(8):15-18.

[8] 王榮惠, 李樹虎, 王丹勇,等. 碳纖維增強(qiáng)復(fù)合材料壓縮力學(xué)性能實(shí)驗(yàn)研究[C]// 全國沖擊動力學(xué)學(xué)術(shù)會議.咸陽, 2013:1-6.

WANG Ronghui, LI Shuhu, WANG Danyong, et al. Experimental Study on Compressive Mechanical Properties of Carbon Fiber Reinforced Composites[C]// National Conference on Impact Dynamics. Xianyang,2013:1-6.

[9] 蘇尚彬, 焦學(xué)健, 薛遠(yuǎn)水,等. 碳纖維材料保險杠耐撞性分析及輕量化[J]. 山東理工大學(xué)學(xué)報(自然科學(xué)版),2015,29(3):49-52.

SU Shangbin, JIAO Xuejian, XUE Yuanshui, et al. Crashworthiness Analysis and Lightweight of Carbon Fiber Bumper[J]. Journal of Shandong University of Technology (Natural Science Edition),2015,29(3):49-52.

[10] JACOB G C, FELLERS J F, SIMUNOVIC S, et al. Energy Absorption in Polymer Composites for Automotive Crashworthiness[J]. Journal of Composite Materials,2002,36(7):813-850.

[11] MAMALIS A G, MANOLAKOS D E, VIEGEL-AHN G L, et al. On the Axial Crumpling of Fibre-reinforced Composite Thin-walled Conical Shells[J]. International Journal of Vehicle Design,1991,12(4):061692.

[12] BORIA S, PETTINARI S. Mathematical Design of Electric Vehicle Impact Attenuators: Metallic vs Composite Material[J]. Composite Structures,2014,115(18):51-59.

[13] 李志斌, 虞吉林, 鄭志軍,等. 薄壁管及其泡沫金屬填充結(jié)構(gòu)耐撞性的實(shí)驗(yàn)研究[J]. 實(shí)驗(yàn)力學(xué),2012,27(1):77-86.

LI Zhibin, YU Jilin, ZHENG Zhijun, et al. Experimental Study on the Crashworthiness of Foam Metal Filled Thin Walled Tubes[J]. Experimental Mechanics,2012,27(1):77-86.

[14] 張勇, 林福泳. 鋁泡沫填充薄壁結(jié)構(gòu)耐撞可靠性優(yōu)化設(shè)計[J]. 機(jī)械工程學(xué)報,2011,47(22):93-99.

ZHANG Yong, LIN Fuyong. Crashworthiness Optimization Design of Aluminum Foam Filled Thin-walled Structures[J]. Chinese Journal of Mechanical Engineering,2011,47(22):93-99.

[15] 宋宏偉, 虞鋼, 范子杰,等. 多孔材料填充薄壁結(jié)構(gòu)吸能的相互作用效應(yīng)[J]. 力學(xué)學(xué)報,2005,37(6):697-703.

SONG Hongwei, YU Gang, FAN Zijie, et al. Interaction Effect of Energy Absorption of Porous Material Filled Thin-walled Structures[J]. Journal of Mechanics,2005,37(6):697-703.

[16] 高偉釗,莫旭輝,付銳,等. 基于Kriging的泡沫填充錐形薄壁結(jié)構(gòu)耐撞性6σ穩(wěn)健性優(yōu)化設(shè)計[J]. 固體力學(xué)學(xué)報,2012,33(4):370-378.

GAO Weizhao, MO Xuhui, FU Rui, et al. Foam Filled Conical Thin-walled Structure Crashworthiness 6 Sigma Robust Optimization Based on Kriging[J]. Journal of Solid Mechanics,2012,33(4):370-378.

[18] 柳艷杰. 汽車低速碰撞吸能部件的抗撞性能研究[D]. 哈爾濱:哈爾濱工程大學(xué),2012.

LIU Yanjie. Study on the Crashworthiness of Automobile Low Speed Collision Energy Absorbing Components[D]. Harbin: Harbin Engineering University,2012.

[19] HALLQUIST J. LS-DYNA Keyword User’s Manual[M]. Livermore: Livermore Software Technology Corporation,2006.

[20] WIERZBICKI T. Energy Absorption of Structures and Materials[J]. International Journal of Impact Engineering,2004,30(7):881-882.

[21] MAMALIS A G, MANOLAKOS D E, DEMOSTHENOUS G A, et al. Energy Absorption Capability of Fibreglass Composite Square Frusta Subjected to Static and Dynamic Axial Collapse[J]. Thin-walled Structures,1996,25:269-295.

[22] SHARGHI H, SHAKOURI M, KOUCHAKZADEH M A. An Analytical Approach for Buckling Analysis of Generally Laminated Conical Shells under Axial Compression[J]. Acta Mechanica,2016,227(4):1-18.

[23] GUPTA N K, VELMURUGAN R, GUPTA S K. An Analysis of Axial Crushing of Composite Tubes[J]. Journal of Composite Materials,1997,31:1262-1286.

[24] BORIA S, PETTINARI S, GIANNONI F. Theoretical Analysis on the Collapse Mechanisms of Thin-walled Composite Tubes[J]. Composite Structures,2013,103(9):43-49.

[25] 賀福. 缺陷是碳纖維的致命傷[J]. 高科技纖維與應(yīng)用,2010(4):25-31.

HE Fu. Defect is the Fatal Wound of Carbonfiber[J]. High Tech Fibers and Applications,2010(4):25-31.