霍銀磊, 劉彥亨, 陳無(wú)忌
(河南科技大學(xué) 包裝工程系, 河南 洛陽(yáng) 471000)
包裝中的很多緩沖材料壓縮時(shí)體現(xiàn)為明顯的正切型彈性特性,例如緩沖氣柱、泡沫橡膠、預(yù)壓后的泡沫塑料等。此類(lèi)材料組成的包裝系統(tǒng)具有較強(qiáng)的非線性,基于小參數(shù)的傳統(tǒng)Lindstedt-Poincare攝動(dòng)法(L-P)對(duì)此類(lèi)問(wèn)題的求解一般不再有效。
由于初始條件的差別,上述針對(duì)非線性系統(tǒng)振動(dòng)響應(yīng)的求解結(jié)果并不能直接應(yīng)用于跌落沖擊問(wèn)題。目前,針對(duì)正切型非線性系統(tǒng)的求解的研究還比較少,而且主要針對(duì)于無(wú)阻尼系統(tǒng):宋浩等[12]基于何氏PEM法分析了無(wú)阻尼正切型非線性包裝系統(tǒng)跌落沖擊響應(yīng),利用沖擊過(guò)程中的能量關(guān)系修正了何氏PEM法的解析表達(dá),得到較高精度的近似解析解;郭蓓蓓等[13-14]分別基于同倫攝動(dòng)法和Li-He氏修正同倫攝動(dòng)法分析了正切型非線性包裝系統(tǒng)跌落沖擊響應(yīng)并對(duì)所得解析解基于能量關(guān)系進(jìn)行了修正;Song[15]還利用修正的同倫攝動(dòng)法分析了雙曲正切型無(wú)阻尼非線性包裝系統(tǒng)的跌落沖擊響應(yīng),得到了系統(tǒng)響應(yīng)的級(jí)數(shù)解;趙曉兵等[16]基于牛頓諧波平衡法(NHB)分析了正切型緩沖系統(tǒng)跌落沖擊響應(yīng)。最近,仲晨等[17]也基于NHB法分析了包含線性阻尼的正切型包裝系統(tǒng)跌落沖擊響應(yīng)。
上述沖擊響應(yīng)求解方法大多需要基于沖擊過(guò)程中的能量變化對(duì)所得解析解的幅值和頻率進(jìn)行修正以提高沖擊響應(yīng)的計(jì)算精度。對(duì)于有阻尼的包裝系統(tǒng),能量修正會(huì)變得困難。因此,考慮到跌落沖擊的瞬時(shí)特性,本文嘗試將結(jié)合多尺度法和改進(jìn)的L-P攝動(dòng)法來(lái)討論正切型含有線型阻尼的非線性包裝系統(tǒng)在發(fā)生跌落沖擊時(shí)的響應(yīng)問(wèn)題。
圖1 包裝系統(tǒng)跌落動(dòng)力學(xué)模型Fig.1 The dropping dynamics model of packaging system
由牛頓第二定律可知系統(tǒng)的運(yùn)動(dòng)方程
(1)
由于變形過(guò)程中總有x (2) 利用式(2),式(1)寫(xiě)為 (3) 其中: 對(duì)于跌落沖擊,系統(tǒng)運(yùn)動(dòng)的初始條件為 (4) 對(duì)于α,β不為小量的非線性系統(tǒng),借鑒改進(jìn)的L-P法,令τ=ωt,對(duì)式(3)進(jìn)行時(shí)間變換得 (5) 式中,“′” 表示為對(duì)新時(shí)間變量τ的微分。 利用多尺度法,分別定義快、慢時(shí)間尺度 T0=τ,T1=ετ,T2=ε2τ (6) 變換后的時(shí)間導(dǎo)數(shù)和因變量的多尺度展開(kāi)分別為 x=x0(T0,T1,T2)+εx1(T0,T1,T2)+ε2x2(T0,T1, T2)+… (7) 定義系統(tǒng)基頻ω,有: (8) ε2x2)+ε2μω[D0+εD1+ε2D2](x0+εx1+ε2x2)+(ω2-εω1-ε2ω2)(x0+εx1+ε2x2)+εα(x0+εx1+ε2x2)3+ε2β(x0+εx1+ε2x2)5=0 (9) 分離參數(shù)ε各階,得到: (10) (11) (12) 系統(tǒng)運(yùn)動(dòng)的初始條件 (13) 由式(10)可方便得到系統(tǒng)的0次近似解 x0=AeiT0+cc (14) 式中:A為復(fù)振幅,包含振幅與相位信息;cc表示前面各項(xiàng)的復(fù)共軛。 將式(14)代入式(11)得 (15) 令D1A=0,根據(jù)消除久期項(xiàng)條件可得 (16) 考慮初始條件,式(11)有如下形式解 (17) 式中:B為復(fù)振幅;cc表示前面各項(xiàng)的復(fù)共軛。 進(jìn)一步地,將式(14)和(17)代入式(12)得 (18) 令D2A=0可知ω2為復(fù)數(shù),方程的解不可得[18]。因此,根據(jù)消除久期項(xiàng)條件,令ω2=0,得 (19) 考慮初始條件,式(12)有如下形式的解 (20) (21) 整理可得式(5)的一次漸進(jìn)解 3 sin(ωt+θ)] (22) 及二次漸進(jìn)解 (23) 其中: 由式(8)及(16)可得 (24) (25) 為了檢驗(yàn)?zāi)P偷臏?zhǔn)確性,不考慮阻尼的影響,利用文獻(xiàn)[13-14]中的系統(tǒng)參數(shù):m=10 kg,H=0.6 m,k0=600 N/cm,μ=0,εα=72 N/cm3,ε2β=0.000 1N /cm5,對(duì)比MSLP法(c=1×10-10N·s/m)與同倫攝動(dòng)解 (HPM)、修正的同倫攝動(dòng)解(CHPM)、Li-He修正同倫解(LHHPM)及修正的Li-He同倫解(MLHHPM)的計(jì)算結(jié)果見(jiàn)表1??梢钥闯觯簩?duì)于無(wú)阻尼系統(tǒng),考慮未修正的沖擊響應(yīng),同倫攝動(dòng)法(HPM)的對(duì)應(yīng)誤差分別為6.14%和11.8%;同樣,未修正的Li-He同倫解(LHHPM)誤差也較大,分別為5.1%和8.5%。相比之下,修正的同倫攝動(dòng)法(CHPM)和修正Li-He同倫解(MLHHPM)的最大位移及加速度的相對(duì)誤差都降低了很多,精度有了很大的提高,顯示了能量修正的優(yōu)越性。本文MSLP計(jì)算的最大位移和加速度響應(yīng)的一階近似解與R-K數(shù)值結(jié)果對(duì)比的相對(duì)誤差分別為1.35%和3.28%,二次近似階的相對(duì)誤差分別為0.62%和1.84%。雖然相對(duì)于修正的同倫攝動(dòng)解精度稍低,考慮到阻尼系統(tǒng)能量修正的復(fù)雜性,本文方法具有更大的實(shí)用價(jià)值。與期望的一樣,相比于一次近似,取MSLP的二次近似解有助于提高系統(tǒng)響應(yīng)的計(jì)算精度。 表1 不同方法計(jì)算的近似解與數(shù)值解的對(duì)比Tab.1 Comparison of solutions obtained by different methods for the cubic-nonlinear system 考慮到文獻(xiàn)[13-14]中針對(duì)式(3)進(jìn)行的參數(shù)選擇,其中參數(shù)的取值并不嚴(yán)格滿足正切型彈性關(guān)系。為更嚴(yán)格的考察正切非線性MSLP解的準(zhǔn)確性, 若不做特殊說(shuō)明,以下分析均基于文獻(xiàn)[17]的系統(tǒng)參數(shù):m=20 kg,g=9.8 m/s2,H=0.6 m,c=1 200 N·s/m,k0=2.5×105N/m,db=0.05 m。系統(tǒng)沖擊響應(yīng)與文獻(xiàn)[17]結(jié)果對(duì)比見(jiàn)圖2??梢钥闯觯罕疚腗SLP法計(jì)算的系統(tǒng)位移響應(yīng)與數(shù)值法得到的結(jié)果吻合很好,在其峰值處誤差最大,其一、二次位移響應(yīng)最大值分別為0.020 55 m和0.020 62 m,較數(shù)值解0.020 76 m相對(duì)誤差僅分別為1.0%和0.67%,較NHB法所得結(jié)果的5.58%的平均誤差精確度更高;同樣,本文方法計(jì)算的系統(tǒng)加速度響應(yīng)與數(shù)值法得到的結(jié)果吻合也較好,在其峰值處誤差最大,其一、二次加速度響應(yīng)最大值分別為-346 m/s2、-344 m/s2,較數(shù)值解-334.2 m/s2的相對(duì)誤差僅有3.53%和2.96%,較NHB法所得結(jié)果的5.61%的平均誤差精確度更高,證明了本文方法對(duì)阻尼正切型非線性系統(tǒng)的有效性。二次MSLP解較一次解精度有所提高,但總體差別不大,考慮計(jì)算的簡(jiǎn)便性和工程實(shí)際,認(rèn)為一次MSLP近似解即可滿足實(shí)際計(jì)算要求。后續(xù)的討論將基于一次MSLP近似解進(jìn)行。 圖2 不同方法得到的系統(tǒng)跌落沖擊響應(yīng)結(jié)果對(duì)比Fig.2 Comparison of dropping shock responses of the packaging system based on different methods 圖3 不同阻尼系統(tǒng)跌落沖擊響應(yīng)(k0=2.5×105 N/m)Fig.3 Dropping shock responses of the damped packaging system for different damping c. (k0=2.5×105 N/m) 對(duì)于小阻尼(c=500 N·s/m)正切型非線性系統(tǒng),不同初始彈性系數(shù)k0對(duì)系統(tǒng)響應(yīng)的影響見(jiàn)圖4。結(jié)果表明:隨著初始彈性系數(shù)逐漸增大,響應(yīng)幅值逐漸減小,響應(yīng)周期逐漸減小。這也與線性系統(tǒng)的結(jié)論一致,較大的初始彈性系數(shù)意味著較大的系統(tǒng)剛度和較大的響應(yīng)頻率;另外,隨著k0的逐漸增大,系統(tǒng)的最大位移響應(yīng)與最大加速度響應(yīng)的相對(duì)誤差也逐漸減小。 圖4 不同初始彈性系數(shù)的系統(tǒng)跌落沖擊響應(yīng)(c=500 N·s/m)Fig.4 Dropping shock responses of the damped packaging system for different initial elastic constant k0.(c=500 N·s/m) 本文基于MSLP法分析了含有阻尼的正切型非線性包裝系統(tǒng)的跌落沖擊響應(yīng)的一、二次近似表達(dá),并與數(shù)值解進(jìn)行對(duì)比,主要結(jié)論如下: (1) MSLP方法可有效用于含有阻尼的正切型非線性包裝系統(tǒng)的跌落沖擊響應(yīng)的近似解析求解,其一、二次近似解都具有較高的精度并且無(wú)需額外的基于能量關(guān)系的幅值與頻率修正,二次近似解的精度相對(duì)于一次近似解有所提高,但其形式也更復(fù)雜。 (2) 隨著系統(tǒng)阻尼的增大,一次MSLP近似解的誤差逐漸增大。對(duì)于實(shí)際的小阻尼緩沖材料系統(tǒng),一次MSLP近似解具有較高的精度。 (3) 隨著系統(tǒng)初始彈性系數(shù)的增大,基于MSLP方法的一次近似解的誤差逐漸減小,解析解的精度提高。 (4) 考慮到雙曲正切彈性材料力-變形關(guān)系有類(lèi)似式(2)形式的泰勒展開(kāi)式,本文方法可方便的應(yīng)用于雙曲正切系統(tǒng)的求解。2 基于MSLP方法的近似解
3 算例分析與討論
4 結(jié) 論