劉 勝,李 鵬,楊翊仁
(西南交通大學(xué)力學(xué)與工程學(xué)院,成都 610031)
圖1為一種新型水力緩沖器結(jié)構(gòu)模型,主要由套筒、活塞、碰撞部件、緩沖頭和緩沖彈簧組成。套筒與活塞桿之間安裝有預(yù)緊緩沖彈簧,使得相對(duì)靜止?fàn)顟B(tài)下套筒底部與活塞頭部緊密貼合。運(yùn)動(dòng)部件安裝在緩沖器上部,并置于靜水中。工作時(shí),運(yùn)動(dòng)部件隨緩沖器一起從一定高度落下。當(dāng)進(jìn)入下落行程末段時(shí),緩沖頭與碰撞部件接觸,并以一定速度和加速度沖擊碰撞部件。沖擊過程中,緩沖頭壓縮碰撞部件彈性件,同時(shí)自身速度不斷減小,直至彈性件被完全壓縮,而后緩沖頭由剛性件提供支撐;另一方面,沖擊時(shí)緩沖頭壓縮活塞桿,從而擠壓套筒內(nèi)的冷卻劑和緩沖彈簧,起到對(duì)運(yùn)動(dòng)部件緩沖的效果。運(yùn)動(dòng)部件完全靜止后,緩沖彈簧推動(dòng)活塞桿和緩沖頭復(fù)位至安裝初始狀態(tài)。
本文依據(jù)水力緩沖器的結(jié)構(gòu)和工作原理,建立下落過程和緩沖過程的結(jié)構(gòu)控制方程和流體控制方程,計(jì)算緩沖過程中緩沖頭對(duì)碰撞部件的沖擊力、套筒內(nèi)部壓強(qiáng)、套筒和活塞位移等參數(shù)的動(dòng)態(tài)變化,分析各節(jié)流孔對(duì)沖擊力和壓強(qiáng)影響,對(duì)緩沖器結(jié)構(gòu)進(jìn)行優(yōu)化。
圖1 水力緩沖器結(jié)構(gòu)
為建立緩沖過程的控制方程,首先給出以下基本假設(shè):①在下落過程中,套筒與活塞為無相對(duì)運(yùn)動(dòng),可視為一個(gè)整體;②各部件之間的碰撞為完全彈性碰撞,忽略碰撞阻尼,碰撞等效為彈簧作用,組件之間的碰撞剛度為組件拉壓剛度的串、并聯(lián)值[1-2];③ 流體不可壓縮,其能量損失僅為局部能量損失;④僅考慮外部流場對(duì)緩沖器的流體形狀阻力。
緩沖器在下落和緩沖過程中的結(jié)構(gòu)示意圖如圖2所示。M1,M2分別為套筒和活塞桿的質(zhì)量;Mf1,Mf2分別為二者排開水的質(zhì)量;Cf,Cs分別為黏性阻力系數(shù)和形狀阻力系數(shù);D1,D2,H1,H2分別為套筒與活塞桿的外直徑和高度;Sht為緩沖頭截面積;K,F(xiàn)0分別為緩沖彈簧剛度和預(yù)緊力;K1,K2,K3分別為組件各部分的碰撞剛度;H為活塞有效工作高度;He為碰撞部件彈性件最大壓縮行程;P1,P'1,P2,P'2分別為作用與套筒和活塞上的流體壓力。
圖2 緩沖器結(jié)構(gòu)簡化示意圖
下落過程結(jié)構(gòu)的運(yùn)動(dòng)方程可寫為[3]:
緩沖過程結(jié)構(gòu)控制方程為
其中:
緩沖頭與碰撞部件之間的沖擊力為
其中δ為Dirichlet函數(shù)。
圖3給出了緩沖過程中流量的計(jì)算模型。套筒與活塞之間的相對(duì)運(yùn)動(dòng)會(huì)擠出套筒主腔內(nèi)的流體或吸入外部流體,因此各節(jié)流孔的流體流向及流量與此相對(duì)運(yùn)動(dòng)相關(guān)。
圖3 流量計(jì)算示意圖
首先假定當(dāng)主腔內(nèi)流體被擠出,且擠出流體均是從活塞頭部流出,而后經(jīng)過各節(jié)流孔流出套筒;當(dāng)流體被吸入主腔時(shí),認(rèn)為所有流體均流向活塞頭部。各節(jié)流孔流量分別記為αQ,βQ,γQ,ηQ,ξQ。忽略非定常項(xiàng)和位置水頭對(duì)總水頭的影響,可由Bernoulli方程[4]沿流線建立能量方程:
流體從節(jié)流孔流出主腔時(shí)
流體從節(jié)流孔流進(jìn)主腔時(shí)
其中:i=1,2,…,5,分別對(duì)應(yīng)5個(gè)不同節(jié)流孔;p,U分別為活塞頭部壓強(qiáng)與流速;pi,Ui分別為流出i號(hào)節(jié)流孔后的流體壓強(qiáng)與流速;∑hm,∑h'm分別為沿流線的所有局部能量損失,流出及流進(jìn)二者不同。由于套筒側(cè)壁開有節(jié)流孔,各條流線上的能量方程與下落過程中活塞頭與各節(jié)流孔的相對(duì)位置相關(guān)。
忽略外部流場的流動(dòng),外部流場滿足
流動(dòng)中的流量方程為
其中:λi= ±1,i=1,2,3,4,5為流向迭代指標(biāo);ΔV為副腔體積變化;λV=±1為副腔體積增大和縮小指標(biāo)。
聯(lián)立式(5)~(8),可得流量控制方程:
流量控制方程中,活塞頭部流體流速應(yīng)等于套筒與活塞的相對(duì)運(yùn)動(dòng)流速,即
采用 Newton迭代方法[5]求解式(9)。由式(9)獲得節(jié)流孔流量后,代入式(5)或(6)即可得到緩沖器內(nèi)部壓強(qiáng)分布。
采用Newmark積分格式[6]對(duì)結(jié)構(gòu)控制方程求解。基于Fortran程序[7]編制整個(gè)求解過程數(shù)值計(jì)算程序。圖4給出了程序計(jì)算流程。圖5給出了當(dāng)下落高度為1.75 m時(shí)沖擊力和活塞頭部壓力的時(shí)程曲線。顯然,在第1次沖擊之后,沖擊力和壓強(qiáng)都迅速降低。此外,當(dāng)活塞經(jīng)過套筒側(cè)壁節(jié)流孔時(shí),壓強(qiáng)會(huì)出現(xiàn)顯著突變。
圖4 緩沖分析的數(shù)值模擬流程
圖6給出了不同下落高度時(shí)最大沖擊力和最大壓強(qiáng)的變化曲線,可見隨著下落高度的增加,最大沖擊力和最大壓強(qiáng)隨之增大。
圖7給出了節(jié)流孔尺寸對(duì)最大沖擊力的影響曲線,節(jié)流孔編號(hào)與圖3中對(duì)應(yīng)。由圖7可以看出:隨著#4、5節(jié)流孔直徑的增大,最大沖擊力減小,#5節(jié)流孔直徑的影響更加明顯。#1節(jié)流孔直徑較小時(shí)對(duì)最大沖擊力影響也很小,當(dāng)直徑大于3 mm后,沖擊力隨之先增大后減小。
圖5 沖擊力、壓強(qiáng)時(shí)程曲線
圖6 最大沖擊力和最大壓強(qiáng)隨下落高度變化
圖7 最大沖擊力隨節(jié)流孔直徑變化
針對(duì)一新型水力緩沖器建立了理論模型。該模型包括結(jié)構(gòu)控制方程、流量控制方程和結(jié)構(gòu)-流量耦合控制方程。并根據(jù)該模型,進(jìn)行了緩沖器的程序?qū)崿F(xiàn)與數(shù)值模擬,得到以下結(jié)論:
1)最大沖擊力與最大壓強(qiáng)出現(xiàn)在沖擊開始的短暫時(shí)間內(nèi),之后,二者均顯著下降。
2)隨著初始下落高度的增加,最大沖擊力和最大壓強(qiáng)隨著增大。
3)套筒頂部節(jié)流孔在尺寸達(dá)到一定程度后對(duì)緩沖影響較大,在側(cè)壁節(jié)流孔中,上方節(jié)流孔對(duì)緩沖影響更大。
[1]孫訓(xùn)方,方孝淑,關(guān)來泰.材料力學(xué)(I)[M].4版.北京:高等教育出版社,2002.
[2]倪振華.振動(dòng)力學(xué)[M].西安:西安交通大學(xué),1989.
[3]金棟平,胡海巖.碰撞振動(dòng)與控制[M].北京:科技出版社,2005.
[4]張鳴遠(yuǎn).流體力學(xué)[M].北京:高等教育出版社,2010.
[5]John H M,Kurtis D F.數(shù)值方法(MATLAB 版)[M].4版.陳渝,錢方,譯.北京:電子工業(yè)出版社,2005.
[6]王勖成.有限單元法[M].北京:清華大學(xué)出版社,2003.
[7]彭國倫.Fortran95程序設(shè)計(jì)[M].北京:中國電力出版社,2006.