杜群貴　孫智權(quán)　黃崇溪　翟曉晨

華南理工大學(xué),廣州,510640

?

氣缸沖擊動(dòng)力學(xué)建模與仿真分析

杜群貴孫智權(quán)黃崇溪翟曉晨

華南理工大學(xué),廣州,510640

將氣缸在行程末端的沖擊問(wèn)題看成一個(gè)由細(xì)長(zhǎng)活塞桿、質(zhì)量負(fù)載、活塞以及端蓋組成的振動(dòng)系統(tǒng)。首先應(yīng)用彈性體的一維波動(dòng)方程建立細(xì)長(zhǎng)活塞桿縱向振動(dòng)變形的數(shù)學(xué)模型，用模態(tài)疊加方法解析了沖擊力、接觸位移以及接觸時(shí)間。然后利用非線性有限元軟件LS-DYNA對(duì)氣缸沖擊過(guò)程進(jìn)行數(shù)值仿真計(jì)算分析，詳細(xì)研究了當(dāng)驅(qū)動(dòng)負(fù)載質(zhì)量以及活塞初始沖擊速度不同時(shí)，活塞沖擊接觸面接觸合力和接觸作用時(shí)間的變化情況。對(duì)比分析了理論解析解與LS-DYNA數(shù)值解，驗(yàn)證了解析解的準(zhǔn)確性，得出了沖擊力分別隨速度和負(fù)載的變化規(guī)律；進(jìn)一步對(duì)沖擊振動(dòng)響應(yīng)進(jìn)行了頻率分析，并研究了氣缸沖擊力的變化原因。研究結(jié)果為氣缸的設(shè)計(jì)和使用提供了理論和實(shí)際參考。

氣缸沖擊力；解析解；LS-DYNA仿真；頻率分析

0　引言

隨著氣動(dòng)技術(shù)的飛速發(fā)展，氣缸廣泛地應(yīng)用于各種工業(yè)生產(chǎn)線上。但是，由于氣動(dòng)系統(tǒng)中介質(zhì)空氣的高度可壓縮性和低黏度等內(nèi)在固有特性，當(dāng)活塞運(yùn)行到行程末端時(shí)就會(huì)產(chǎn)生沖擊端蓋現(xiàn)象。同時(shí)，氣缸也不可避免地會(huì)由于安裝制造誤差、供氣壓力波動(dòng)、緩沖不良、負(fù)載慣性以及運(yùn)行一段時(shí)間之后緊固螺栓松動(dòng)等外在因素，產(chǎn)生活塞沖擊端蓋現(xiàn)象[1]，并伴隨著不良振動(dòng)與噪聲，這也是導(dǎo)致氣缸失效的最主要原因。

氣缸的沖擊問(wèn)題是一個(gè)經(jīng)典的桿結(jié)構(gòu)彈性碰撞力學(xué)問(wèn)題。近年來(lái)，國(guó)內(nèi)外學(xué)者針對(duì)桿的彈性碰撞問(wèn)題作了很多研究。文獻(xiàn)[2]采用波傳播法給出了質(zhì)點(diǎn)和均勻桿縱向碰撞響應(yīng)解析解；文獻(xiàn)[3-4]也采用波傳播法解析了一般邊界條件下桿的縱向沖擊振動(dòng)問(wèn)題 ，但解析式隨著波傳播次數(shù)的增加而變得過(guò)于復(fù)雜，其工程應(yīng)用難以推廣。文獻(xiàn)[5-6]采用模態(tài)分析方法解析了桿的點(diǎn)彈性碰撞問(wèn)題，但局限于特殊的邊界條件；文獻(xiàn)[7-8]采用直接模態(tài)疊加法將碰撞問(wèn)題轉(zhuǎn)化為具有初速度的振動(dòng)問(wèn)題，解析了不同邊界條件、多種截面形狀的桿類(lèi)沖擊問(wèn)題的頻率方程及動(dòng)力響應(yīng)。但是，對(duì)于類(lèi)似氣缸存在復(fù)雜邊界和初始條件下的沖擊破壞現(xiàn)象的建模研究甚少，對(duì)其沖擊力大小的預(yù)估大都基于經(jīng)驗(yàn)。

基于以上原因，本文在不考慮緩沖墊和氣緩沖的情況下對(duì)氣缸沖擊現(xiàn)象進(jìn)行簡(jiǎn)化建模研究，求解出沖擊力的解析式，并與有限元仿真數(shù)值解進(jìn)行了對(duì)比，從而驗(yàn)證解析式的可靠性，此外，還分析了沖擊力波動(dòng)變化的規(guī)律及原因。

1　氣缸碰撞系統(tǒng)建模與解析

本文所研究的標(biāo)準(zhǔn)型單桿雙作用氣缸結(jié)構(gòu)簡(jiǎn)圖如圖1所示。

1.質(zhì)量負(fù)載　2.活塞桿　3.端蓋　4.活塞　5.缸筒圖 1　氣缸結(jié)構(gòu)簡(jiǎn)圖

該氣缸缸筒與無(wú)桿側(cè)端蓋壓成一體，缸筒與有桿端蓋3處的連接為螺紋連接，軸向尺寸短，質(zhì)量比其他系列氣缸小10%～15%。氣缸最大工作壓力為1.0 MPa,活塞工作速度為50～1000 mm/s，主要運(yùn)用于裝配生產(chǎn)線上輕型物料輸送。

筆者在該氣缸可靠性試驗(yàn)研究的過(guò)程中發(fā)現(xiàn)氣缸工作一段時(shí)間之后，會(huì)出現(xiàn)端蓋與缸筒螺紋連接松脫失效以及氣缸安裝緊固螺栓松動(dòng)的現(xiàn)象。這是由于行程末端活塞對(duì)端蓋的長(zhǎng)期沖擊造成的，故對(duì)此沖擊現(xiàn)象進(jìn)行建模研究。

對(duì)氣缸碰撞系統(tǒng)數(shù)學(xué)建模時(shí)，為了便于簡(jiǎn)化分析，將質(zhì)量負(fù)載看成細(xì)長(zhǎng)活塞桿左端部的一個(gè)集中質(zhì)量，將活塞看成活塞桿右端部的一個(gè)集中質(zhì)量，集中質(zhì)量隨桿端橫截面一起振動(dòng)。將端蓋與活塞的碰撞接觸部分看成一個(gè)等效接觸彈簧，活塞、活塞桿和負(fù)載以一定的初速度沖擊接觸彈簧，彈簧壓力等效為氣缸沖擊力[9]。并認(rèn)為沖擊過(guò)程中桿的橫截面始終保持為平面，且界面上的應(yīng)力分布均勻，忽略各部分重力的影響。其碰撞系統(tǒng)簡(jiǎn)化模型如圖2所示。

1.質(zhì)量負(fù)載　2.活塞桿　3.等效接觸彈簧　4.活塞圖 2　碰撞系統(tǒng)簡(jiǎn)化模型

1.1碰撞系統(tǒng)數(shù)學(xué)模型的建立

圖2中活塞桿長(zhǎng)為l，橫截面積為A，質(zhì)量密度為ρ，彈性模量為E，等效接觸彈簧剛度系數(shù)為k，負(fù)載質(zhì)量為M，活塞質(zhì)量為m，將質(zhì)量負(fù)載位置處桿端作為坐標(biāo)原點(diǎn)，分別建立桿的縱向碰撞問(wèn)題的控制方程、邊界條件和初始條件[10]。

桿的自由振動(dòng)方程:

(1)

式中，u為桿的縱向振動(dòng)位移；x為桿的某一橫截面距離坐標(biāo)原點(diǎn)的位移；t為振動(dòng)時(shí)間。

碰撞過(guò)程中邊界條件:

(2)

(3)

初始條件:

u(x,0)=0

(4)

(5)

(6)

邊界條件轉(zhuǎn)化為

(7)

(8)

初始條件轉(zhuǎn)化為

(9)

(10)

1.2碰撞模型的求解

(11)

(12)

(13)

(14)

其中，λ為無(wú)因次頻率，由式(11)得到主振型通解為

U(ξ)=acosλξ+bsinλξ

(15)

式中，a、b分別為主振型中余弦、正弦信號(hào)的幅值。

將式(15)代入式(13)、式(14)得

(16)

由式(16)得到頻率方程：

(17)

由式(15)、式(17)得到各階主振型函數(shù)：

(18)

由式(12)解得各階時(shí)間函數(shù)：

Ti(τ)=Cicosλiτ+Disinλiτ

(19)

式中，Ci、Di分別為第i階時(shí)間函數(shù)中余弦信號(hào)和正弦信號(hào)的幅值。

根據(jù)模態(tài)疊加法，桿的振動(dòng)位移為無(wú)窮個(gè)主振動(dòng)的疊加：

(20)

1.3主振型的正交性求解

為了求得各階時(shí)間函數(shù)中系數(shù)Ci、Di，首先應(yīng)求得主振型之間的正交性。設(shè)Ui(ξ)、Uj(ξ)分別是對(duì)應(yīng)于無(wú)因次頻率λi及λj的主振型，由式(11)有

(21)

(22)

對(duì)式(21)等號(hào)兩邊同時(shí)乘以Uj(ξ)并對(duì)其從0到1求積分，得

(23)

利用分部積分法可將式(23)等號(hào)左邊寫(xiě)為

(24)

聯(lián)立式(23)、式(24)并化簡(jiǎn)，得

(25)

同理，對(duì)式(23)化簡(jiǎn)得

(26)

(27)

為了便于連續(xù)性求解，引入δ函數(shù)，得到主振型的正交性關(guān)系式：

(28)

(29)

1.4沖擊接觸求解

由初始條件可得

(30)

(31)

利用主振型的正交性，將式(30)、式(31)等號(hào)兩邊同時(shí)乘以[α+δ(ξ)+αηδ(ξ-1)]Ui(ξ)，然后從0到1區(qū)間求定積分，可以求解得到:

Ci=0

(32)

(33)

則桿的振動(dòng)位移解析式為

(34)

在活塞沖擊端蓋處，有x=l，將其代入式(34)，可得到?jīng)_擊處桿端橫截面振動(dòng)位移：

(35)

用等效接觸彈簧所受的壓力大小來(lái)表征碰撞過(guò)程中沖擊力大小，得到?jīng)_擊力表達(dá)式：

(36)

2　基于LS-DYNA的有限元仿真分析

以圖1所示的典型氣缸結(jié)構(gòu)為基礎(chǔ)，忽略結(jié)構(gòu)中倒角、圓角、氣孔等因素，根據(jù)該氣缸實(shí)際結(jié)構(gòu)尺寸，建立簡(jiǎn)化的有限元仿真模型。由于氣缸為軸對(duì)稱(chēng)性結(jié)構(gòu)，利用其對(duì)稱(chēng)性可以簡(jiǎn)化模擬過(guò)程和減少自由度，并縮短分析所需計(jì)算時(shí)間。本文所采用的1/2模型如圖3所示。

圖 3　氣缸活塞沖擊端蓋有限元模型

活塞和活塞桿的材料為結(jié)構(gòu)鋼，端蓋和缸筒的材料為鋁合金，有限元仿真時(shí)均采用線彈性材料模型。其材料性能參數(shù)見(jiàn)表1。

表 1　材料性能參數(shù)

氣缸裝配體模型均采用sweep掃掠方式劃分網(wǎng)格，網(wǎng)格均為八節(jié)點(diǎn)六面體單元。將缸筒與端蓋、活塞與活塞桿以及質(zhì)量負(fù)載與活塞桿的連接設(shè)置成固連綁定接觸，活塞桿與端蓋導(dǎo)向環(huán)的相對(duì)滑動(dòng)以及活塞與端蓋的碰撞接觸采用自動(dòng)面面雙向接觸(ASTS)。將缸筒右端面和端蓋左端面設(shè)置成固定約束，將對(duì)稱(chēng)面設(shè)置成對(duì)稱(chēng)約束，保證活動(dòng)部件可在對(duì)稱(chēng)平面內(nèi)平動(dòng)?；钊麠U與導(dǎo)向環(huán)、活塞與端蓋的摩擦因數(shù)均為0.1。在不考慮重力作用下，給活塞、活塞桿以及質(zhì)量負(fù)載一定的初速度去沖擊端蓋，仿真得到碰撞沖擊響應(yīng)。

為了有效地處理氣缸沖擊過(guò)程中的接觸問(wèn)題、非線性問(wèn)題和隨時(shí)間變化很快的沖擊力學(xué)過(guò)程，LS-DYNA求解器采用顯式直接積分方法，采用罰函數(shù)法處理沖擊過(guò)程中的動(dòng)態(tài)接觸，并指定輸出RCFORC界面合力ASCII文件，獲得沖擊力的時(shí)間歷程曲線[11]。

氣缸在實(shí)際工作過(guò)程中，因應(yīng)用場(chǎng)合不同，氣缸搭載的負(fù)載和運(yùn)行速度就不一樣，運(yùn)行到行程末端時(shí)，必然產(chǎn)生不同的沖擊情況。本文分析了氣缸在不同速度和不同質(zhì)量下的沖擊狀況，得到了沖擊力響應(yīng)曲線。

3　結(jié)果分析

本文研究的典型氣缸活塞桿直徑d=12 mm，活塞桿長(zhǎng)l=193 mm，氣缸內(nèi)直徑D=32 mm，端蓋長(zhǎng)度l1=29 mm，活塞寬度l2=17 mm，活塞質(zhì)量m=0.08 kg，仿真質(zhì)量負(fù)載長(zhǎng)度l3=50 mm，活塞與活塞桿質(zhì)量比η=0.5，工作行程s=100 mm。

氣缸的碰撞為活塞與端蓋的短暫接觸，由于活塞長(zhǎng)度較端蓋短，而其彈性模量比端蓋大得多，使得活塞的抗壓剛度比端蓋的抗壓剛度大得多，物理建模時(shí)將活塞看成一個(gè)集中質(zhì)量，而將端蓋看成接觸彈簧；此外，由于活塞運(yùn)行到行程末端時(shí)，有桿腔內(nèi)氣體迅速排空，故忽略氣體對(duì)接觸剛度的耦合影響，其等效接觸剛度為

1.12×109(N/m)

其中，E1為端蓋材料彈性模量，Acon為活塞與端蓋碰撞有效接觸面積，碰撞處端蓋內(nèi)半徑r=10.3 mm。

將上述參數(shù)代入解析解中，利用MATLAB編程計(jì)算，取前100階(已趨于穩(wěn)定)模態(tài)疊加，得到理論解析曲線。

3.1沖擊力與速度的關(guān)系

為了驗(yàn)證解析解表達(dá)式中沖擊力與沖擊速度的線性關(guān)系，在負(fù)載為20 N的情況下，分別取活塞沖擊速度為0.1 m/s、0.2 m/s、0.3 m/s、0.4 m/s、0.5 m/s進(jìn)行有限元仿真計(jì)算，得到在不同速度下的沖擊力-時(shí)間曲線，如圖4所示。并將0.5 m/s沖擊速度下的仿真曲線與理論曲線進(jìn)行對(duì)比，如圖5所示。

圖 4　不同速度下的沖擊力-時(shí)間曲線

圖 5　仿真曲線與理論曲線對(duì)比(負(fù)載為20 N，沖擊速度為0.5 m/s)

從圖4可看出，在不同沖擊速度下，有限元仿真得到的沖擊力波形隨時(shí)間同步變化，且沖擊力隨速度線性增大；同時(shí)，還可以看出，在同一負(fù)載下，沖擊作用時(shí)間并不隨著沖擊速度的變化而變化，而是保持不變。從圖5可看出，同一沖擊條件下，LS-DYNA仿真曲線與理論解析曲線十分吻合，表明仿真結(jié)果與解析表達(dá)式一致，從而驗(yàn)證了解析表達(dá)式中沖擊力與速度的關(guān)系。

3.2沖擊力與負(fù)載的關(guān)系

通過(guò)改變程序中質(zhì)量比α的大小，得到一系列不同負(fù)載下的沖擊力響應(yīng)曲線。同時(shí)，為了進(jìn)一步研究活塞沖擊端蓋的沖擊情況，利用LS-DYNA仿真分析了負(fù)載分別為10 N、20 N(圖5)、30 N、40 N、50 N在沖擊速度為0.5 m/s時(shí)的沖擊過(guò)程，經(jīng)過(guò)后處理得到?jīng)_擊力-時(shí)間歷程曲線。將相應(yīng)的理論解析曲線與LS-DYNA仿真曲線進(jìn)行對(duì)比分析，結(jié)果對(duì)比如圖6～圖9所示。

圖 6　仿真曲線與理論曲線對(duì)比(負(fù)載為10 N，沖擊速度為0.5 m/s)

圖 7　仿真曲線與理論曲線對(duì)比(負(fù)載為30 N，沖擊速度為0.5 m/s)

圖 8　仿真曲線與理論曲線對(duì)比(負(fù)載為40 N，沖擊速度為0.5 m/s)

圖 9　仿真曲線與理論曲線對(duì)比(負(fù)載為50 N，沖擊速度為0.5 m/s)

結(jié)合圖5～圖9，對(duì)比LS-DYNA仿真數(shù)值解與理論解析解可知，兩者得出的活塞沖擊面沖擊力時(shí)間歷程曲線在變化趨勢(shì)、峰值大小、峰值數(shù)量以及沖擊作用時(shí)間等方面，理論解析與有限元仿真十分吻合，進(jìn)一步說(shuō)明了上述物理模型建立以及模型簡(jiǎn)化處理的正確性和可靠性。

同時(shí)，從圖6～圖9可以看出，沖擊力隨時(shí)間的變化曲線并不具有單峰值和完全對(duì)稱(chēng)性，也并不是先增大后減小，而是表現(xiàn)為多峰值且與時(shí)間不對(duì)稱(chēng)的波動(dòng)形式。隨著負(fù)載質(zhì)量的不斷增加，峰值和峰值個(gè)數(shù)也不斷增加，沖擊力曲線波動(dòng)越來(lái)越“劇烈”，沖擊作用時(shí)間也不斷變長(zhǎng)，這是由于振動(dòng)頻率的變化所致。

3.3振動(dòng)頻率分析

沖擊接觸面的振動(dòng)為細(xì)長(zhǎng)活塞桿各階主振動(dòng)在桿碰撞端x=l處的疊加，其振動(dòng)位移響應(yīng)是由具有各階固有頻率的諧波信號(hào)組成的。不同負(fù)載質(zhì)量下的前6階無(wú)因次固有頻率及相應(yīng)的振幅見(jiàn)表2。

表2不同負(fù)載下沖擊處振動(dòng)位移中前六階振動(dòng)頻率與幅值

負(fù)載(N)頻率/幅值一階二階三階四階五階六階0λi1.407883.741245.493988.1835111.2031314.29036Ai0.158060.137720.040710.003420.000561.50581×10-410λi0.379982.756514.563756.776329.6956312.7555Ai0.267550.126060.100410.010890.001272.76765×10-420λi0.272852.733524.554136.765439.6872112.74897Ai0.360110.126880.101330.011020.001272.77676×10-4

從表2可以看出，隨著頻率階次的增加，其對(duì)應(yīng)的幅值不斷衰減，且從第四階頻率開(kāi)始，衰減得很快。同時(shí)，隨著負(fù)載質(zhì)量的不斷增加，同一階次的頻率值有所減小，但同一階次的幅值有所增大，且第一階頻率和幅值對(duì)負(fù)載質(zhì)量的變化最敏感。但是，隨著階數(shù)的增加，各階頻率值和幅值對(duì)負(fù)載質(zhì)量的變化越來(lái)越不敏感，使得不同負(fù)載質(zhì)量下的各階頻率值和幅值都十分接近。

負(fù)載質(zhì)量的增加，使得振動(dòng)位移響應(yīng)中低頻(尤其是第一階頻率)成分所占比重不斷增加，從而導(dǎo)致沖擊力幅值的不斷增大；各階頻率的減小，使得沖擊作用時(shí)間不斷增大，沖擊力波動(dòng)變化越來(lái)越“劇烈”。

在氣缸短暫劇烈的沖擊過(guò)程中，高頻響應(yīng)衰減得十分迅速。為了使解析式進(jìn)一步推廣到實(shí)際運(yùn)用，在不增大負(fù)載質(zhì)量的情況下，通過(guò)取前5階模態(tài)和前100階模態(tài)所得到的沖擊力解析曲線與LS-DYNA仿真曲線對(duì)比分析，如圖10所示。

圖10　仿真曲線與理論曲線對(duì)比(負(fù)載為0,沖擊速度為0.5 m/s)

圖10結(jié)果顯示,三條曲線十分吻合；同時(shí),模態(tài)階數(shù)的減少,使得沖擊力幅值和沖擊作用時(shí)間稍有延長(zhǎng)。但隨著負(fù)載質(zhì)量的增加，兩種階次取值所得到的沖擊力響應(yīng)曲線趨于重合，差距完全可以忽略不計(jì)。這是由于隨著負(fù)載質(zhì)量的增加，高頻成分衰減得越來(lái)越迅速，高階振型中基本不引起振動(dòng)。因此，實(shí)際工程運(yùn)用(如螺紋連接設(shè)計(jì)所考慮的最大沖擊力、端蓋螺紋孔的大小和深度以及安裝方式)中，只需要考慮前面幾階振型分析就足夠準(zhǔn)確了。

4　結(jié)論

(1)本文對(duì)氣缸行程末端活塞沖擊端蓋的現(xiàn)象進(jìn)行動(dòng)力學(xué)建模，將其看成一個(gè)振動(dòng)系統(tǒng)，采用模態(tài)分析方法推導(dǎo)出沖擊力響應(yīng)解析解，計(jì)算結(jié)果與LS-DYNA仿真分析結(jié)果吻合，證明了本文簡(jiǎn)化建模及其解析解的正確性。

(2)沖擊力與活塞沖擊速度成線性關(guān)系；隨著負(fù)載質(zhì)量的增加，沖擊力幅值不斷增大，且波動(dòng)得更“劇烈”。沖擊作用時(shí)間隨負(fù)載質(zhì)量的增加而變長(zhǎng)，但與與沖擊速度無(wú)關(guān)。

(4)沖擊響應(yīng)中第一階模態(tài)對(duì)負(fù)載變化最敏感，但隨著階次的增加，高階模態(tài)對(duì)負(fù)載變化越來(lái)越不敏感；同時(shí)，高階頻率成分衰減得很快，實(shí)際運(yùn)用中只需要取幾階低頻分析即可。

(5)研究結(jié)果為氣缸的可靠性研究提供了理論依據(jù)，為氣缸的結(jié)構(gòu)(尤其是緩沖結(jié)構(gòu))設(shè)計(jì)、活塞螺紋連接、使用安裝、維護(hù)保養(yǎng)提供了理論依據(jù)和數(shù)據(jù)支持。

[1]廖理.沖擊負(fù)載下氣缸性能測(cè)試平臺(tái)的研究與實(shí)現(xiàn)[D].成都:電子科技大學(xué),2013.

[2]Love A E H. A Treatise on the Mathematical Theory of Elasticity [M].New York: Dover, 1994.

[3]張繼業(yè),曾京,舒仲周.桿的縱向沖擊振動(dòng)[J].振動(dòng)與沖擊, 1999,18(3):58-63.

Zhang Jiye,Zeng Jing,Shu Zhongzhou.Longitudinal Vibration of Prismatic Bar during Impact[J].Journal of Vibration and Shock,1999,18(3):58-63.

[4]Gan Chunbiao, Wei Yimin, Yang Shixi. Longitudinal Wave Propagation in a Rod with Variable Cross-section[J]. Journal of Sound and Vibration, 2014,333(2) :434-445.

[5]諸德超,邢譽(yù)峰. 點(diǎn)彈性碰撞問(wèn)題之解析解[J].力學(xué)學(xué)報(bào), 1996,28(1):99-103.

Zhu Dechao, Xing Yufeng. Analytical Solution of Point Elastic Impact between Structures[J].Chinese Journal of Theoretical and Applied Mechanics,1996, 28(1):99-103.

[6]吳家強(qiáng),王宏志. 桿的縱向沖擊問(wèn)題全過(guò)程分析[J]. 振動(dòng)與沖擊,2004,23(1):103-106.

Wu Jiaqiang, Wang Hongzhi. General Procedure for Analysis of Axial Response of Prismatic Bar during Impact[J].Journal of Vibration and Shock, 2004, 23(1):103-106.

[7]鮑四元. 彈性碰撞問(wèn)題的模態(tài)疊加法及哈密頓體系下的計(jì)算研究[D].西安:西北工業(yè)大學(xué),2005.

[8]Li Q S,Wu J R, Xu J. Longitudinal Vibration of Multi-step Non-uniform Structures with Lumped Masses and Spring Supports[J]. Applied Acoustics,2002,63(3):333-350.

[9]孫龍.潛孔錘及聲波鉆柔體振動(dòng)沖擊的理論與建模[D].北京:中國(guó)地質(zhì)大學(xué),2014.

[10]倪振華.振動(dòng)力學(xué)[M].西安:西安交通大學(xué)出版社,1990.

[11]王澤鵬,胡仁喜,康士廷,等.ANSYS14.5/LS-DYNA非線性有限元分析實(shí)例指導(dǎo)教程[M].北京:機(jī)械工業(yè)出版社,2013.

(編輯陳勇)

Impact Dynamics Modeling and Simulation Analyses of Pneumatic Cylinders

Du QunguiSun ZhiquanHuang ChongxiZhai Xiaochen

South China University of Technology,Guangzhou,510640

Impact of pneumatic cylinder at the stroke terminal was considered to be an integrated system of vibrations, including slender piston rod, mass load, piston and end cover.Based on one-dimensional wave equation, a mathematical model of longitudinal vibration of slender piston rod was established firstly. Impact force,contact displacement and contact time were analyzed by means of mode superposition method.Then impact process of pneumatic cylinder was numerically simulated by nonlinear FEM software LS-DYNA.The variation of impact interface force and contact time were studied in detail under different masses of driven load and impact velocities. Compared with numerical solution of LS-DYNA,the accuracy of analytical solution was verified, and how impact force varies as velocity and load was concluded. Vibration frequencies of impact response were furtherly analyzed, and reasons for the changes of impact forces were studied. Research conclusions provide reference for design and applications of pneumatic cylinders.

impact force of pneumatic cylinder; analytical solution; LS-DYNA simulation; frequency analysis

杜群貴,男,1965年生。華南理工大學(xué)機(jī)械與汽車(chē)工程學(xué)院教授、博士研究生導(dǎo)師。主要研究方向?yàn)楝F(xiàn)代機(jī)電工程設(shè)計(jì)與控制。孫智權(quán)，男，1991年生。華南理工大學(xué)機(jī)械與汽車(chē)工程學(xué)院碩士研究生。黃崇溪，男，1990年生。華南理工大學(xué)機(jī)械與汽車(chē)工程學(xué)院碩士研究生。翟曉晨，男，1988年生。華南理工大學(xué)機(jī)械與汽車(chē)工程學(xué)院博士研究生。

2015-04-10

TH138.51;O326

10.3969/j.issn.1004-132X.2016.08.011