亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        正切非線性包裝系統(tǒng)跌落沖擊的MSLP解

        2022-06-17 03:15:30霍銀磊劉彥亨陳無(wú)忌
        振動(dòng)與沖擊 2022年11期
        關(guān)鍵詞:方法系統(tǒng)

        霍銀磊, 劉彥亨, 陳無(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)模型及運(yùn)動(dòng)方程

        圖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)

        2 基于MSLP方法的近似解

        對(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)

        3 算例分析與討論

        為了檢驗(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)

        4 結(jié) 論

        本文基于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)的求解。

        猜你喜歡
        方法系統(tǒng)
        Smartflower POP 一體式光伏系統(tǒng)
        WJ-700無(wú)人機(jī)系統(tǒng)
        ZC系列無(wú)人機(jī)遙感系統(tǒng)
        基于PowerPC+FPGA顯示系統(tǒng)
        學(xué)習(xí)方法
        半沸制皂系統(tǒng)(下)
        連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
        可能是方法不對(duì)
        用對(duì)方法才能瘦
        Coco薇(2016年2期)2016-03-22 02:42:52
        四大方法 教你不再“坐以待病”!
        Coco薇(2015年1期)2015-08-13 02:47:34
        在线视频播放观看免费| 国产精自产拍久久久久久蜜| 精品国产福利在线观看网址2022 | 人人玩人人添人人澡| 一本一本久久久久a久久综合激情| 亚洲av一区二区网址| 亚洲色图视频在线免费看| 日韩中文字幕免费视频| 99热精品成人免费观看| 国内精品熟女一区二区| 久久久国产精品123| 亚洲av麻豆aⅴ无码电影| 久久久久亚洲AV无码专| 日韩男女av中文字幕| 欧美巨鞭大战丰满少妇| 亚洲国产韩国欧美在线| 亚洲av人妖一区二区三区| 亚洲无人区乱码中文字幕动画| 国产极品女主播国产区| 无码国产精品一区二区vr老人| 日韩美女高潮流白浆视频在线观看| 黄片免费观看视频播放| 久久久久国产精品| 亚洲男同志gay 片可播放| 久草国产手机视频在线观看| 午夜福利影院成人影院| 一本一道av中文字幕无码| 亚洲AV无码精品色午夜超碰| 国产一区二区三区小向美奈子| 亚洲av乱码一区二区三区按摩| 无码丰满少妇2在线观看| 亚洲国产精品国自产拍av在线| 亚洲精品午夜久久久九九| 男女车车的车车网站w98免费| 国产偷2018在线观看午夜| 丝袜美腿亚洲综合在线播放| 熟妇熟女乱妇乱女网站| 91视频88av| 性色av一区二区三区四区久久| 欧美成人猛片aaaaaaa| 久久亚洲av成人无码国产|