非線性系統(tǒng)帶集中質(zhì)量懸臂梁易損件跌落沖擊特性

非線性系統(tǒng)帶集中質(zhì)量懸臂梁易損件跌落沖擊特性

郝蒙1，陳安軍1,2

(1.江南大學(xué)包裝工程系，江蘇無(wú)錫2141222； 2.國(guó)家輕工業(yè)包裝制品質(zhì)量監(jiān)督檢測(cè)中心，江蘇無(wú)錫214122)

摘要：以三次非線性緩沖包裝系統(tǒng)為研究對(duì)象，建立系統(tǒng)跌落沖擊動(dòng)力學(xué)方程；綜合龍格-庫(kù)塔和有限元法設(shè)計(jì)動(dòng)力學(xué)方程的求解算法；探討了集中質(zhì)量和主體振動(dòng)頻率對(duì)易損件響應(yīng)的影響。數(shù)值分析表明：系統(tǒng)跌落沖擊過(guò)程，易損件最大位移和加速度響應(yīng)均位于懸臂梁自由端，易損件與主體連接部?jī)?nèi)應(yīng)力最大；易損件相對(duì)主體質(zhì)量較小時(shí)，系統(tǒng)耦合作用對(duì)易損件響應(yīng)影響不明顯；隨集中質(zhì)量增加或主體振動(dòng)頻率向易損件第一階固有頻率接近時(shí)，易損件內(nèi)應(yīng)力響應(yīng)、自由端相對(duì)位移響應(yīng)幅值顯著增大。在含彈性易損件產(chǎn)品系統(tǒng)包裝設(shè)計(jì)中，連接部的內(nèi)應(yīng)力和易損件相對(duì)位移是關(guān)注的重要參數(shù)。研究結(jié)論可為帶集中質(zhì)量懸臂梁產(chǎn)品緩沖包裝設(shè)計(jì)提供理論依據(jù)。

關(guān)鍵詞：集中質(zhì)量；懸臂梁；易損件；有限元；三次非線性；跌落沖擊

中圖分類號(hào):TB485.3;O322文獻(xiàn)標(biāo)志碼：A

基金項(xiàng)目：國(guó)家自然科學(xué)基金青年基金資助項(xiàng)目(11002084)；上海市教委創(chuàng)新基金資助項(xiàng)目(12YZ092, 12YZ074) ；滬江基金(D14005)；上海理工大學(xué)研究生創(chuàng)新基金項(xiàng)目

收稿日期：2014-03-26修改稿收到日期：2014-07-23

Dropping shock characteristics of a cubic nonlinear system with a cantilever beam type elastic critical component attaching a concentrated tip mass

HAOMeng1,CHENAn-jun1,2(1. Department of Packaging Engineering, Jiangnan University, Wuxi 214122, China; 2. China National Light Industry Packaging Quality Supervision & Inspection Center, Wuxi 214122, China)

Abstract:In order to investigate the dropping shock characteristics of a cubic nonlinear system with a cantilever beam critical component attaching a concentrated tip mass, the dynamic model of the whole system was built. Runge-Kutta method and the finite element method were used to solve the system’s dynamic equations, and the effects of the system parameters, such as, concentrated tip mass and natural frequencies of the main part on the responses of the critical component were studied. The case study showed that the maximum dropping displacement and acceleration responses of the critical component occur at its free end, while the maximum internal stress appears at the joint point between the critical component and the main part; the effects of the coupling action of the system on the responses of the critical component can be ignored when its mass is much less than that of the main component; with increase in concentrated tip mass and/or the natural frequency of the main part close to the 1st natural frequencu of the critical component, the free end of the critical component’s response amplitudes increase obviously; in the packaging design of systems with elastic critical components, the internal stress at the joint point between the critical component and the main part and the relative displacement of the critical component are be parameters to evaluate if products are damaged. The results provided a theoretical foundation for the package cushioning design of this type systems.

Key words:concentrated tip mass; cantilever beam; critical component; finite element; cubic nonlinear; dropping impact

研究包裝系統(tǒng)的動(dòng)力學(xué)模型及跌落條件下系統(tǒng)的動(dòng)力學(xué)特性是產(chǎn)品緩沖包裝設(shè)計(jì)的理論基礎(chǔ)。Mindlin[1]首次提出易損件的概念，建立兩種包裝系統(tǒng)的動(dòng)力學(xué)模型(單自由度和二自由度集中質(zhì)量-彈簧系統(tǒng)模型)，研究線性和非線性系統(tǒng)跌落沖擊動(dòng)力學(xué)特性。此后，考慮易損件的二自由度集中質(zhì)量-彈簧模型的動(dòng)力學(xué)特性研究得到不斷豐富和發(fā)展[2-5]。孫勇[6]、王軍[7]和高德[8]分別對(duì)多自由度包裝系統(tǒng)振動(dòng)和沖擊動(dòng)力學(xué)特性進(jìn)行了分析，Suhir[9]和Wong[10-12]針對(duì)帶彈性板易損件的產(chǎn)品包裝系統(tǒng)，建立動(dòng)力學(xué)模型，分析彈性板的跌落沖擊特性，結(jié)果表明對(duì)于帶彈性易損件的產(chǎn)品，由于易損件為無(wú)限自由度，傳統(tǒng)集中質(zhì)量-彈簧模型應(yīng)用有失準(zhǔn)確性；奚德昌等[13-14]針對(duì)具有簡(jiǎn)支梁彈性部件工藝品包裝系統(tǒng)，建立剛體-彈性體模型，探究了緩沖材料參數(shù)對(duì)易損件和系統(tǒng)響應(yīng)的影響；高德等[15-17]分別以桿式、簡(jiǎn)支梁和懸臂梁3種不同類型易損件產(chǎn)品與發(fā)泡聚乙烯緩沖材料組成的耦合系統(tǒng)為研究對(duì)象，建立了系統(tǒng)的動(dòng)力學(xué)方程，利用顯式有限差分法研究了跌落沖擊過(guò)程中易損件位移響應(yīng)、加速度響應(yīng)和應(yīng)力分布的特點(diǎn)。對(duì)包含集中質(zhì)量懸臂梁彈性易損件產(chǎn)品包裝系統(tǒng)的研究現(xiàn)未見(jiàn)報(bào)道，且此類易損件在機(jī)電產(chǎn)品中較為常見(jiàn)，其包裝防護(hù)理論分析有待進(jìn)一步探討。

本文以三次非線性包裝系統(tǒng)為研究對(duì)象，產(chǎn)品易損件為自由端有集中質(zhì)量的懸臂梁，建立系統(tǒng)跌落沖擊下的動(dòng)力學(xué)模型；龍格-庫(kù)塔法與有限元法結(jié)合，設(shè)計(jì)系統(tǒng)動(dòng)力學(xué)方程的數(shù)值分析算法。探討易損件位移、加速度響應(yīng)以及內(nèi)應(yīng)力變化規(guī)律，并分析其影響因素。研究方法和結(jié)論可為此類產(chǎn)品緩沖包裝設(shè)計(jì)提供理論依據(jù)。

1系統(tǒng)動(dòng)力學(xué)模型

圖1　跌落條件下考慮易 損件包裝系統(tǒng)力學(xué)模型 Fig.1 The dropping model of packaging system with critical component

圖2　帶集中質(zhì)量懸臂梁易損件有限單元?jiǎng)澐?Fig.2 The finite element division of cantilever beam component with concentrated tip mass

(1a)

(1b)

式中：E、I、ρ和A分別為易損件彈性模量、橫截面慣性矩、材料密度和截面面積。

易損件動(dòng)能T和應(yīng)變能U為

(2)

以u(píng)、w分別表示易損件節(jié)點(diǎn)絕對(duì)位移和相對(duì)位移向量，即

單元節(jié)點(diǎn)位移向量與總節(jié)點(diǎn)位移向量關(guān)系為

ui=Ciu，wi=Ciw

(3)

式中:Ci為4×(2n+2)的矩陣，記為第i個(gè)單元位移提取矩陣，具體為

將式(3)和Ci代入式(2)，易損件動(dòng)能和應(yīng)變能分別為

(4)

式中:M和K均為2(n+1)×2(n+1)的矩陣，分別表示易損件質(zhì)量矩陣和剛度矩陣，即

(5)

考慮自由端的集中質(zhì)量影響，將m2添加到總質(zhì)量矩陣的第(2n+1)×(2n+1)元素上。

由拉格朗日方程建立易損件的動(dòng)力學(xué)方程

(6)

式中:w=u-y1(t)v，v=[1010…1010]T。

系統(tǒng)跌落過(guò)程中，可將易損件在主體位移激勵(lì)下的響應(yīng)分析轉(zhuǎn)化為求解易損件在慣性力作用下的相對(duì)位移響應(yīng)問(wèn)題。將u=w+y1(t)v代入方程(6)得

(7)

考慮易損件對(duì)主體的影響，跌落沖擊條件下主體的動(dòng)力學(xué)方程為

(8)

若不考慮易損件對(duì)主體的影響，初始條件不變，主體動(dòng)力學(xué)方程為

(9)

同時(shí)，由變分迭代法[19]可得的主體頻率近似表達(dá)式

(10)

2動(dòng)力學(xué)方程求解

對(duì)二階常微分方程，龍格-庫(kù)塔數(shù)值分析方法具有很好的精度[20]。方程(9)可直接采用龍格-庫(kù)塔法求解，由于方程(8)包含易損件對(duì)主體的影響，需有限差分離散后再使用龍格-庫(kù)塔法求解。有限差分離散方法如下。

(11)

(12)

式中：B=EIτ2/(ρAl4)。進(jìn)一步

(13)

(14)

應(yīng)用Newmark法[23]求解；對(duì)應(yīng)廣義特征值方程可表示為

K*φ=ω2M*φ

(15)

采用子空間迭代法[23]求解。

(16)

式中：J(n×2n)為橫向相對(duì)位移提取矩陣，即

易損件橫向絕對(duì)位移向量可表示為

u2=w2+y1(:,t)

(17)

易損件內(nèi)應(yīng)力為

(18)

式中:r為易損件橫截面半徑。

3數(shù)值計(jì)算結(jié)果及其分析

某小型機(jī)電產(chǎn)品，主體質(zhì)量m1=10 kg；懸臂梁(45號(hào)鋼)橫截面半徑r=0.003 m，長(zhǎng)L=0.1 m，材料密度ρ=7 850 kg/m3，彈性模量E=200 GPa，彈性極限σe=180 MPa，材料安全系數(shù)取1.2，許用內(nèi)應(yīng)力為150 MPa；集中質(zhì)量m2=0.02 kg。緩沖材料初始彈性常數(shù)k0=100 N/cm，非線性常數(shù)e=72 N/cm3，系統(tǒng)跌落高度H=0.6 m。

圖3　易損件跌落沖擊響應(yīng) Fig.3 The dropping shock response of critical component

表1　易損件跌落沖擊響應(yīng)最值

考慮易損件與主體耦合作用，易損件位移、加速度響應(yīng)和懸臂梁內(nèi)應(yīng)力分布如圖3所示，相對(duì)位移、絕對(duì)加速度及內(nèi)應(yīng)力最值如表1所示。易損件絕對(duì)位移響應(yīng)和絕對(duì)加速度響應(yīng)沿懸臂梁長(zhǎng)度方向變化不明顯；相對(duì)位移響應(yīng)和絕對(duì)加速度響應(yīng)最值出現(xiàn)在自由端，分別為0.41 mm和-560.02 m/s2；內(nèi)應(yīng)力最大值81.60 MPa出現(xiàn)在與主體連接部。傳統(tǒng)脆值理論以加速度是否超過(guò)許用值作為判斷產(chǎn)品失效標(biāo)準(zhǔn)，而對(duì)于含彈性易損件產(chǎn)品，根據(jù)彈性材料失效準(zhǔn)則，當(dāng)最大應(yīng)力超過(guò)其彈性極限，材料將發(fā)生塑性變形，產(chǎn)品失效，同時(shí)考慮到產(chǎn)品內(nèi)部構(gòu)件間隙，避免碰撞，應(yīng)綜合考慮最大應(yīng)力是否超過(guò)許用值和易損件相對(duì)位移是否超過(guò)許用值作為判定該產(chǎn)品失效的準(zhǔn)則。

圖4　不同集中質(zhì)量易損件跌落沖擊響應(yīng) Fig.4 The dropping shock response of critical component with different concentrated tip mass

考慮易損件和主體耦合作用，集中質(zhì)量取不同值(0.01、0.02、0.03和0.04 kg)時(shí)易損件自由端相對(duì)位移、相對(duì)加速度和根部?jī)?nèi)應(yīng)力變化如圖4。隨著集中質(zhì)量增加，易損件響應(yīng)顯著增加；集中質(zhì)量為0.04 kg時(shí)，內(nèi)應(yīng)力及相對(duì)位移最大值分別為138.44 MPa和0.73 mm，分別是集中質(zhì)量為0.01 kg時(shí)的2.53倍和2.7倍。最大值出現(xiàn)時(shí)間和響應(yīng)周期基本不變。

考慮耦合作用，易損件響應(yīng)最值及出現(xiàn)時(shí)間如表2所示。集中質(zhì)量為0.02 kg時(shí)，易損件總質(zhì)量(0.042 3 kg)遠(yuǎn)小于主體質(zhì)量(10 kg)，耦合作用對(duì)易損件的響應(yīng)影響不明顯；當(dāng)集中質(zhì)量為0.08 kg時(shí)，相對(duì)位移及絕對(duì)加速度最大值誤差達(dá)到5%以上，耦合作用影響不能忽略。

表2　考慮和不考慮耦合作用易損件跌落沖擊響應(yīng)最值

表3　主體振動(dòng)頻率

圖5　不同主體頻率易損件跌落沖擊響應(yīng) Fig.5 The dropping shock response of critical component with different main component frequency

初始條件和產(chǎn)品主體質(zhì)量不變，緩沖材料參數(shù)e取不同的值，相應(yīng)主體的響應(yīng)頻率通過(guò)式(10)求得如表3所示；集中質(zhì)量m2=0.02 kg時(shí)，由式(14)模態(tài)分析可求得懸臂梁易損件第一階頻率為195.63 Hz。圖5為主體運(yùn)動(dòng)頻率取不同值時(shí)，易損件自由端相對(duì)位移、相對(duì)加速度和根部?jī)?nèi)應(yīng)力變化。當(dāng)主體振動(dòng)頻率為26.91時(shí)，內(nèi)應(yīng)力最大值為185.46 MPa已超過(guò)許用應(yīng)力150 MPa，緩沖材料需重新設(shè)計(jì)；且自由端相對(duì)位移最大值為0.94 mm，為主體振動(dòng)頻率14.32 Hz時(shí)的2.29倍。主體振動(dòng)頻率向易損件一階頻率靠近時(shí)，易損件自由端相對(duì)位移、相對(duì)加速度響應(yīng)以及根部?jī)?nèi)應(yīng)力幅值均顯著增大，易損件響應(yīng)對(duì)主體振動(dòng)頻率變化敏感。

4結(jié)論

以易損件為帶集中質(zhì)量的懸臂梁產(chǎn)品包裝系統(tǒng)為研究對(duì)象，建立了多自由度剛體-彈性體動(dòng)力學(xué)模型和運(yùn)動(dòng)方程；結(jié)合龍格-庫(kù)塔法和有限元法提出解決此類問(wèn)題的有效數(shù)值算法，主要結(jié)論如下：

(1)帶集中質(zhì)量懸臂梁易損件的包裝系統(tǒng)，易損件自由端位移和加速度響應(yīng)最大，內(nèi)應(yīng)力最??；而與主體連接部?jī)?nèi)應(yīng)力最大。

(2)集中質(zhì)量對(duì)易損件位移、加速度響應(yīng)及內(nèi)應(yīng)力峰值影響顯著。當(dāng)集中質(zhì)量為0.04 kg時(shí)，內(nèi)應(yīng)力及相對(duì)位移最大值分別為138.44 MPa和0.73 mm，是集中質(zhì)量為0.01 kg時(shí)的2.53倍和2.7倍。

(3)考慮動(dòng)力學(xué)方程耦合作用具有實(shí)際意義。集中質(zhì)量為0.02 kg時(shí)，耦合作用對(duì)易損件的響應(yīng)影響不明顯；當(dāng)集中質(zhì)量為0.08 kg時(shí)，相對(duì)位移及絕對(duì)加速度最大值誤差達(dá)到5%以上，耦合作用影響不能忽略。

(4)易損件響應(yīng)峰值與周期對(duì)主體振動(dòng)頻率敏感。當(dāng)主體頻率向易損件的第一階固有頻率靠近時(shí)，易損件相對(duì)位移和內(nèi)應(yīng)力峰值增加尤為顯著，峰值出現(xiàn)時(shí)間及響應(yīng)周期同樣變化明顯。

對(duì)含彈性易損件產(chǎn)品包裝設(shè)計(jì)中，需根據(jù)彈性材料失效準(zhǔn)則，將彈性易損件的最大內(nèi)應(yīng)力是否超過(guò)許用值，作為判斷該類產(chǎn)品是否失效的標(biāo)準(zhǔn)之一，同時(shí)考慮相對(duì)位移是否超過(guò)允許值，作為產(chǎn)品是否失效的判定標(biāo)準(zhǔn)，使得產(chǎn)品包裝系統(tǒng)設(shè)計(jì)及評(píng)價(jià)更具有針對(duì)性。緩沖包裝設(shè)計(jì)中，緩沖材料的選擇應(yīng)使主體的振動(dòng)頻率遠(yuǎn)離易損件的固有頻率。

參考文獻(xiàn)

[1]Mindlin R D. Dynamics of packaging cushioning[J]. Bell System Technical Journal, 1945, 25:353-445．

[2]王志偉, 胡長(zhǎng)鷹. 洗衣機(jī)包裝系統(tǒng)在矩形波沖擊下的特性研究[J]. 包裝工程, 1997, 18(2,3):14-16.

WANG Zhi-wei, HU Chang-ying. The Shock characteristics of washing machine under the action of rectangular pulse[J]. Packaging Engineering, 1997, 18(2,3):14-16.

[3]Wang Zhi-wei, Hu Chang-ying. Shock spectra and damage boundary curves for nonlinear package cushioning system[J]. Packaging Technology and Science, 1999, 12(5): 207-217.

[4]王軍,王志偉.半正弦脈沖激勵(lì)下考慮易損件的正切型包裝系統(tǒng)沖擊特性研究[J].振動(dòng)與沖擊,2008,27(1):167-168.

WANG Jun, WANG Zhi-wei.3-Dimensional shock response spectra characterizing shock response of a tangent packaging system with critical component[J]. Journal of Vibration and Shock, 2008, 27(1):167-168.

[5]段寧寧,陳安軍.矩形脈沖激勵(lì)下斜支承系統(tǒng)易損件的沖擊特性研究[J].包裝工程,2012,34(7):21-24.

DUAN Ning-ning, CHEN An-jun. Shock characteristics of vulnerable components of tilted support system under rectangular pulse excitation[J]. Packaging Engineering, 2012, 34(7):21-24.

[6]孫勇, 熊彪, 王石. 多自由度緩沖包裝系統(tǒng)的平均能量設(shè)計(jì)法[J]. 包裝工程, 1997, 18(1):17-20.

SUN Yong, XIONG Biao, WANG Shi. The averaging energy design method used for a multi-freedom-degree cushioning packaging system[J]. Packaging Engineering, 1997, 18(1):17-20.

[7]王軍, 王志偉. 多層堆碼包裝系統(tǒng)沖擊動(dòng)力學(xué)特性研究(I): 沖擊譜[J]. 振動(dòng)與沖擊, 2008, 27(8):106-107.

WANG Jun, WANG Zhi-wei. Combined shock spectrum of linear stacking packaging system[J]. Journal of Vibration and Shock, 2008, 27(8):106-107.

[8]高德, 盧富德. 聚乙烯緩沖材料多自由度跌落包裝系統(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.

[9]Suhir E, Burke R. Dynamic response of a rectangular plate to a shock load, with application to portable electronic products[J]. Components, Packaging and Manufacturing Technology, Part B: Advanced Packaging, IEEE Transactions on, 1994,17(3):449-460.

[10]Wong E H. Dynamics of board-level drop impact[J]. Journal of Electronic Packaging, 2005, 127(3): 200-207.

[11]Wong E H, Mai Y W, Seah S K. Board level drop impact-fundamental and parametric analysis[J]. Journal of Electronic Packaging, 2005, 127(4):496-502.

[12]Wong E H, Mai Y W. New insights into board level drop impact[J]. Microelectronics Reliability, 2006, 46(5):930-938.

[13]奚德昌, 彭南陵. 具有簡(jiǎn)支梁彈性部件包裝物品的防震[J]. 包裝工程, 1989, 10(1):42-48.

XI De-chang, PENG Nan-ling. Shock proof of the packaging for elastic parts with simple beam[J]. Packaging Engineering, 1989, 10(1):42-48.

[14]彭南陵, 奚德昌. 具有簡(jiǎn)支梁彈性部件包裝物品的緩沖[J]. 南昌大學(xué)學(xué)報(bào)(工程技術(shù)版), 1994, 16(1):43-48.

PENG Na-ling, XI De-chang. Cushioning of packaged articles with an elastic element of simply supported beam[J]. Journal of Nanchang University, 1994, 16(1):43-48.

[15]盧富德,陶偉明,高德.具有簡(jiǎn)支梁式易損部件的產(chǎn)品包裝系統(tǒng)跌落沖擊研究[J].振動(dòng)與沖擊,2012,31(5):79-81.

LU Fu-de, TAO Wei-ming, GAO De.Drop impact analysis on item packaging system with beam type elastic critical component[J]. Journal of Vibration and Shock, 2012, 31(15):79-81.

[16]高德,盧富德.基于桿式彈性易損部件的非線性系統(tǒng)跌落沖擊研究[J].振動(dòng)與沖擊,2012,31(15):7-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.

[17]Gao D, Lu F D, Chen S J. Drop impact analysis of cushioning system with an elastic critical component of cantilever beam type[J]. Mathematical Problems in Engineering,2013, Article(379068):1-5.

[18]Rao S S. The finite element method in engineering[M]. Butterworth-heinemann, 2005.

[19]陳安軍. 非線性包裝系統(tǒng)跌落沖擊問(wèn)題變分迭代法[J]. 振動(dòng)與沖擊, 2013, 32(18):105-107.

CHEN An-jun. Variational iteration method for droppingshock problem of a cubic non-linear packaging system[J]. Journal of Vibration and Shock, 2013, 32(18):105-107.

[20]陳安軍. 矩形脈沖激勵(lì)下斜支承彈簧系統(tǒng)沖擊特性的研究[J]. 振動(dòng)與沖擊, 2010, 29(10):225-227.

CHEN An-jun. Shock characteristics of a tilted support spring system under action of a rectangular pulse[J]. Journal of Vibration and Shock, 2010, 29(10):225-227.

[21]杜紹洪. 梁自由橫振動(dòng)方程的有限差分方法[J]. 重慶交通大學(xué)學(xué)報(bào):自然科學(xué)版, 2012, 31(1):6-10.

DU Shao-hong. Finite difference method of free transverse vibration equations of beams[J]. Journal of Chongqing Jiaotong University:Natural Science, 2012, 31(1):6-10.

[22]丁科, 殷水平. 有限單元法[M]. 2版.北京:北京大學(xué)出版社, 2012.

[23]王勖成, 邵敏. 有限單元法基本原理和數(shù)值方法[M]. 2版.北京:清華大學(xué)出版社, 1991.

第一作者胡育佳男，副教授，1979年生