盧富德，高 德

(浙江大學(xué) 寧波理工學(xué)院，寧波 315100)

在緩沖包裝或頭盔領(lǐng)域，一些學(xué)者研究了層狀緩沖結(jié)構(gòu)的壓縮響應(yīng)及吸收能量特性。Sek[1]試驗(yàn)方法研究多層瓦楞紙板結(jié)構(gòu)的緩沖性能，拓展了瓦楞紙板作為緩沖包裝材料的可能。盧富德[2-3]提出虛擬質(zhì)量方法，求解多層瓦楞紙板靜態(tài)、動(dòng)態(tài)壓縮響應(yīng),Lu等[4]研究發(fā)泡聚苯乙烯與發(fā)泡聚乙烯串聯(lián)結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)。Coelho等[5]運(yùn)用顯式有限元研究EPS與軟木串聯(lián)應(yīng)用的頭盔結(jié)構(gòu)，文獻(xiàn)[6]考慮外包裝箱作用，研究瓦楞紙板與發(fā)泡聚乙烯的沖擊響應(yīng)及結(jié)構(gòu)的優(yōu)化設(shè)計(jì)。

以上文獻(xiàn)研究對(duì)象都是層狀結(jié)構(gòu)，最近出現(xiàn)功能梯度緩沖材料，此功能梯度材料的密度隨厚度方向按照一定的函數(shù)變化，它克服了層狀結(jié)構(gòu)界面易開裂等缺點(diǎn)[7]。Cui等[8-9]用有限元方法研究層狀泡沫與梯度泡沫在頭盔中的優(yōu)化,把梯度泡沫按照近似的層狀結(jié)構(gòu)進(jìn)行處理。報(bào)道的文獻(xiàn)對(duì)梯度泡沫的壓縮響應(yīng)研究還比較少，僅限于有限元方法，還需要尋求梯度泡沫的響應(yīng)簡(jiǎn)便方法。

鑒于此，本文運(yùn)用Runge-Kutta法研究功能梯度泡沫緩沖結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)規(guī)律，為工程優(yōu)化設(shè)計(jì)提供有效方法。

1 功能梯度泡沫一維本構(gòu)

典型泡沫結(jié)構(gòu)應(yīng)力-應(yīng)變曲線由線彈性區(qū)、屈服平臺(tái)區(qū)及壓實(shí)區(qū)組成。泡沫彈性模量、屈服強(qiáng)度與實(shí)體塑料彈性模量、屈服強(qiáng)度有關(guān)[8]，即

(1a)

(1b)

(1c)

式中：Ef為泡沫彈性模量，Es為塑料實(shí)體的彈性模量，σf為泡沫平臺(tái)應(yīng)力，σs為塑料實(shí)體的屈服應(yīng)力，ρf為泡沫的密度，ρs為塑料實(shí)體的密度，Φ為單位體積上邊上固體部分，1-Φ為單位體積上面上固體部分。

如圖1所示，雙曲正切函數(shù)與正切函數(shù)的特點(diǎn)，典型泡沫的本構(gòu)關(guān)系可由下式表示[10]：

(2)

圖1 雙曲函數(shù)與正切函數(shù)的特征

泡沫結(jié)構(gòu)線彈性、屈服段及壓實(shí)段的本構(gòu)可以由下式表示[10]

(3)

式中：α為控制應(yīng)力-應(yīng)變曲線平臺(tái)的斜率。

由聚苯乙烯的Es=3.3 GPa,σys=15 MPa,ρs=1 050 kg/m3,Φ=0.52，α=4 kPa等參數(shù)[9]，由式(3)可得近似發(fā)泡聚苯乙烯EPS的應(yīng)力-應(yīng)變曲線，如圖2所示。圖2所示，隨泡沫密度的增加，泡沫屈服平臺(tái)應(yīng)力增加，其壓實(shí)應(yīng)變卻減小。

圖2 應(yīng)力與應(yīng)變、相對(duì)密度的關(guān)系

2 功能梯度泡沫系統(tǒng)模型

質(zhì)量為m的物體，以一定的速度沖擊梯度泡沫緩沖結(jié)構(gòu)，為了建模的方便，把泡沫分成n層，近似模擬連續(xù)的功能梯度函數(shù)變化，如圖3所示為功能梯度泡沫系統(tǒng)動(dòng)力學(xué)模型。

梯度函數(shù)為指數(shù)函數(shù)，下式所示[7]

(4)

由式(4)得到梯度泡沫平均密度ρa(bǔ)ve

(5)

圖3 功能梯度泡沫系統(tǒng)沖擊模型

式中：ρfi(i=1,2,…n)為第i層的密度，ρ1為上層的密度，ρ2為最下層的密度，M為指數(shù)。

由牛頓定律得系統(tǒng)動(dòng)力學(xué)方程：

式中：m1為每一層的等效質(zhì)量，A為泡沫的橫截面積，h為每層的厚度，

初始條件為：

y1(0)=y2(0)=…=yn(0)

(7)

3 數(shù)值結(jié)果

m=6 kg,m1=0.01 kg,n=11,h=0.003 64 m,Es=3.3 GPa,σys=15 MPa,ρs=1 050 kg/m3,Φ=0.52，α=4 kPa。

圖曲線(ρa(bǔ)ve=40 kg/m3 ,Δρ=50 kg/m3)實(shí)線：梯度泡沫系統(tǒng)響應(yīng),點(diǎn)劃線：均勻泡沫響應(yīng)響應(yīng)

取V=1.98,2.8,3.42,3.95,4.42 m/s，求質(zhì)量塊的加速度響應(yīng)，加速度是衡量緩沖沖擊響應(yīng)的一個(gè)標(biāo)準(zhǔn)[10-11]，本文研究加速度響應(yīng)與一些參數(shù)的關(guān)系，本文僅考慮壓縮階段的沖擊響應(yīng)，不考慮泡沫回彈階段的結(jié)果。

如圖4所示，具有均勻密度為40 kg/m3的泡沫，當(dāng)m=6 kg的物體沖擊時(shí)，其最大響應(yīng)加速度為49.34 g(g=9.8 m/s2),不隨沖擊速度的變化而變化；當(dāng)沖擊速度1.98,2.8,3.42,3.95與4.42 m/s時(shí)，其最大加速度分別為32.15,40.67,48.02 g,55.84和61.77 g,沖擊能量較小時(shí)，梯度泡沫響應(yīng)加速度小于均勻泡沫，故梯度泡沫對(duì)較小能量有利。

取M=0.2,0.5,1,2,3,4,5和10，分析梯度函數(shù)對(duì)梯度泡沫結(jié)構(gòu)沖擊響應(yīng)的影響,在V=1.4～5.4 m/s范圍,也就是初始輸入能量為5.88～88.2 J時(shí), 取Δρ分別為25,50和75 kg/m3,對(duì)平均密度為70 kg/m3的梯度泡沫進(jìn)行分析,以說明梯度Δρ對(duì)泡沫沖擊響應(yīng)的影響。以梯度泡沫系統(tǒng)響應(yīng)加速度除以均勻泡沫沖擊響應(yīng)為歸一化的值作為應(yīng)變量，以指數(shù)M及輸入能量為自變量，得到結(jié)果如圖5所示。

(a) Δρ=25 kg/m3 (b) Δρ=50 kg/m3 (c) Δρ=75 kg/m3

(a) ρ1=50 kg/m3 (b) ρ1=60 kg/m3 (c) ρ1=70 kg/m3

在圖5中，當(dāng)M=1，即為線性梯度的泡沫系統(tǒng)，在輸入能量為5.88時(shí)，其均一化加速度值最小，隨Δρ=25、50和75 kg/m3變化時(shí)，均一化加速度值分別為0.84、0.75及0.56，可以看出，梯度增大，小沖擊能量沖擊系統(tǒng)時(shí)，梯度泡沫系統(tǒng)響應(yīng)會(huì)變小。

圖6所示，ρa(bǔ)ve=100 kg/m3,ρ1分別取50、60與70 kg /m3時(shí)，均一化加速度與梯度函數(shù)的關(guān)系，在所給的例子中，均一化加速度隨M的增大而減小。

4 結(jié) 論

基于泡沫的特點(diǎn)，建立了有效的本構(gòu)關(guān)系。根據(jù)對(duì)泡沫系統(tǒng)動(dòng)力學(xué)響應(yīng)分析的結(jié)果，得到梯度泡沫系統(tǒng)對(duì)小能量的沖擊響應(yīng)比均勻泡沫密度系統(tǒng)要小，克服了均勻泡沫系統(tǒng)對(duì)一定沖擊能量范圍響應(yīng)幾乎不變化的特點(diǎn)，增加了泡沫的使用范圍，為梯度泡沫的運(yùn)用提供有效方法。

[1]Sek M, Rouillard V,Tarash H, et al. Enhancement of cushioning performance with paperboard crumple inserts[J]. Packaging Technology and Science,2005, 18(5): 273-278.

[2]盧富德，陶偉明，高德.串聯(lián)緩沖結(jié)構(gòu)壓縮響應(yīng)虛擬質(zhì)量分析[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版)，2012,46(8):1431-1436.

LU Fu-de,TAO Wei-ming, GAO De. Compression responses of series cushioning structures by a virtual mass method[J]. Journal of Zhejiang University (Engineering Science), 2012,46(8):1431-1436.

[3]盧富德，陶偉明，高德.瓦楞紙板串聯(lián)緩沖系統(tǒng)動(dòng)力學(xué)響應(yīng)[J]. 振動(dòng)與沖擊，2012 ，31 (21): 30-32.

LU Fu-de,TAO Wei-ming, GAO De. Dynamic response of a series cushioning packaging system made of multi-layer corrugated paperboard[J]. Journal of Vibration and Shock, 2012 ，31 (21): 30-32.

[4]Lu F D,Tao W M,Gao D.Virtual mass method for solution of dynamic response of composite cushion packaging system[J]. Packaging technology and Science,2013,26(S1):32-42.

[5]Coelho R M,Alves D, Sousa R J, et al. New composite liners for energy absorption purposes [J]. Materials & Design,2013, 43 :384-392.

[6]盧富德，陶偉明，高德. 串聯(lián)緩沖系統(tǒng)沖擊響應(yīng)與結(jié)構(gòu)優(yōu)化分析[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版),2012,46(10):1773-1777.

LU Fu-de,TAO Wei-ming, GAO De. Impact response of series cushioning system and structure optimization analysis[J]. Journal of Zhejiang University (Engineering Science),2012,46(10):1773-1777.

[7]Kieback B, Neubrand A, Riedel H. Processing techniques for functionally graded materials [J]. Materials Science and Engineering A,2003,362(1-2): 81-105.

[8]Cui L, Rueda M A F,Gilchrist M D. Optimisation of energy absorbing liner for equestrian helmets. Part II: Functionally graded foam liner[J]. Materials & Design,2009,30(9):3414-3419.

[9]Cui L, Kiernan S, Gilchrist M D. Designing the energy absorption capacity of functionally graded foam materials[J]. Materials Science and Engineering A,2009, 507(1-2): 215-225.

[10]高德，盧富德.考慮轉(zhuǎn)動(dòng)的雙曲正切與正切組合模型緩沖系統(tǒng)沖擊響應(yīng)研究[J]. 振動(dòng)工程學(xué)報(bào)，2012,25(1)：6-11.

GAO De, LU Fu-de. The shock response of hyperbolic tangent and tangent comprehensive model on cushion system considering rotary motion[J]. Journal of Vibration Engineering, 2012,25(1)：6-11.

[11]高德，盧富德.基于桿式彈性易損部件的非線性系統(tǒng)跌落沖擊研究[J]. 振動(dòng)與沖擊，2012 ，31 (15): 47-49 ,59.

GAO De, LU Fu-de. Drop impact analysis of packaging system with bar type elastic critical components[J]. Journal of Vibration and Shock, 2012 ，31 (15): 47-49 ,59.

[12]高德，盧富德.聚乙烯緩沖材料多自由度跌落包裝系統(tǒng)優(yōu)化設(shè)計(jì)[J].振動(dòng)與沖擊，2012,31(3)：69-72.

GAO De, LU Fu-de. Optimization design of MDOF package cushioning system made of polyethylene[J]. Journal of Vibration and Shock, 2012,31(3)：69-72.