某火箭炮小行程液壓緩沖裝置研究

李壯壯,王 璐,林家輝,莫宗來(lái),李 軍

(1.南京理工大學(xué) 機(jī)械工程學(xué)院,江蘇 南京 210094;2.中國(guó)電子科技集團(tuán)公司光電研究院,天津 300308)

將電磁彈射技術(shù)應(yīng)用于野戰(zhàn)火箭炮,可以有效解決發(fā)射初期的做功問(wèn)題,節(jié)省推進(jìn)劑質(zhì)量,提高射程和增大有效載荷質(zhì)量[1]。某野戰(zhàn)火箭炮采用新型電磁彈射發(fā)射方式[2],其發(fā)射裝置由定子、動(dòng)子、火箭彈以及導(dǎo)軌等4部分組成。發(fā)射時(shí),定子固定在發(fā)射裝置上,動(dòng)子帶動(dòng)火箭彈沿導(dǎo)軌運(yùn)動(dòng)并賦予其初速。為克服火箭彈與動(dòng)子自身較大的慣性力,往往需要較大的推動(dòng)力。根據(jù)牛頓第三定律,必定存在一個(gè)與推動(dòng)力大小相等、方向相反的反作用力(本文稱為后坐力)作用在導(dǎo)軌上。該力將對(duì)火箭炮的剛性連接部分造成沖擊,使結(jié)構(gòu)發(fā)生較大的變形甚至發(fā)生強(qiáng)度破壞。常規(guī)野戰(zhàn)火箭炮在發(fā)射過(guò)程中受到的最大作用力是燃?xì)馍淞鳑_擊力,且該力僅作用于火箭彈離開(kāi)定向器后的階段;對(duì)于中、大口徑火箭彈,該沖擊力可達(dá)十幾噸[3]。而采用電磁彈射發(fā)射方式時(shí),其推動(dòng)力帶給火箭炮的等值反向作用力可達(dá)50～60噸,且作用于整個(gè)推進(jìn)過(guò)程。因此,與常規(guī)自力發(fā)射相比,電磁彈射發(fā)射方式的作用力大,且作用時(shí)間長(zhǎng)。為防止該作用力對(duì)火箭炮可能造成的發(fā)射系統(tǒng)結(jié)構(gòu)變形和強(qiáng)度破壞,可以采取兩種解決方法:一是增大火箭炮的結(jié)構(gòu)尺寸來(lái)提高剛強(qiáng)度[4];二是通過(guò)延長(zhǎng)作用時(shí)間來(lái)減緩沖擊力的作用強(qiáng)度。對(duì)于野戰(zhàn)火箭武器而言,增大結(jié)構(gòu)尺寸將會(huì)大大增加系統(tǒng)的無(wú)效質(zhì)量,既影響作戰(zhàn)效率,又增加運(yùn)載難度[5]。因此,研究如何通過(guò)延長(zhǎng)作用時(shí)間實(shí)現(xiàn)電磁彈射式火箭炮的后坐力緩沖有重要意義。

在火炮領(lǐng)域已有較成熟的減小巨大后坐力沖擊的裝置[6](稱為反后坐裝置),常見(jiàn)結(jié)構(gòu)形式有彈簧式、氣體式、液壓式以及混合式等,可針對(duì)不同的口徑、射速和應(yīng)用環(huán)境采用相應(yīng)的反后坐裝置[7]。由于火炮的緩沖行程較大,而火箭炮各部分之間是固連或桁架結(jié)構(gòu),在發(fā)射過(guò)程不允許出現(xiàn)發(fā)射系統(tǒng)大的整體后坐位移,因此,火炮的反后坐裝置并不能直接應(yīng)用于火箭炮。文獻(xiàn)[8]建立了新型電磁軌道炮的后坐仿真模型,利用ANSYS軟件計(jì)算了各部位的受力情況,得出了后坐力構(gòu)成及狀態(tài)。依據(jù)粘滯流體從缸桶中流經(jīng)阻尼孔或間隙時(shí),隨著流體截面的變化產(chǎn)生局部阻尼、引起能量損失的原理設(shè)計(jì)了反后坐裝置。該研究利用液體阻力進(jìn)行緩沖,效果較好,但設(shè)計(jì)的反后坐裝置同樣具有較大的緩沖行程,不適于火箭武器。

為解決某電磁彈射火箭炮后坐力巨大的問(wèn)題,本文應(yīng)用粘滯流體流經(jīng)小孔產(chǎn)生阻尼力原理,提出了一種小行程的液壓緩沖裝置,通過(guò)延長(zhǎng)作用時(shí)間,使火箭炮在一定后坐行程下的受力得到緩沖,從而可有效克服后坐力巨大對(duì)發(fā)射裝置結(jié)構(gòu)和強(qiáng)度的破壞,保證電磁彈射技術(shù)在野戰(zhàn)火箭發(fā)射中的應(yīng)用。

1 緩沖裝置工作原理及緩沖行程確定

1.1 緩沖裝置工作原理

本文設(shè)計(jì)的緩沖裝置的工作原理如圖1所示。由圖1可見(jiàn),緩沖裝置的一端與導(dǎo)軌定子由鋼架連接,另一端連接地面。依據(jù)動(dòng)量定理,發(fā)射時(shí)后坐力傳遞至緩沖裝置,液體粘滯阻力延長(zhǎng)了作用時(shí)間,減小了沖擊力。

緩沖裝置必須在一定的緩沖行程下才能發(fā)揮作用,因此緩沖行程是本文的關(guān)鍵設(shè)計(jì)指標(biāo)。野戰(zhàn)火箭武器通常采用由45鋼、Q345等普通鋼材焊接而成的桁架結(jié)構(gòu),其自身結(jié)構(gòu)強(qiáng)度有限,在發(fā)射過(guò)程不允許出現(xiàn)大的整體結(jié)構(gòu)后坐位移,只允許在一定的強(qiáng)度條件下發(fā)生小的彈性變形。因此,如何將工作行程限定在結(jié)構(gòu)的彈性變形范圍內(nèi),通過(guò)校核負(fù)載作用下結(jié)構(gòu)允許的彈性變形,來(lái)確定緩沖裝置的工作行程,是本文的首要設(shè)計(jì)任務(wù)。

1.2 緩沖行程確定

本文采用仿真分析法,利用有限元軟件Abaqus對(duì)結(jié)構(gòu)在負(fù)載作用下的受力情況進(jìn)行分析,從而在保證強(qiáng)度的情況下確定緩沖行程。針對(duì)本文研究問(wèn)題,對(duì)發(fā)射裝置模型進(jìn)行合理簡(jiǎn)化,如圖2所示。其中,導(dǎo)軌簡(jiǎn)化為長(zhǎng)方形箱體,并對(duì)其材料屬性進(jìn)行自定義,使其質(zhì)量、質(zhì)心位置以及轉(zhuǎn)動(dòng)慣量均與導(dǎo)軌真實(shí)情況相符。除導(dǎo)軌外,發(fā)射裝置各部分材料均設(shè)置為45鋼。高低機(jī)與起落架之間的關(guān)節(jié)軸承設(shè)置為接觸關(guān)系;耳軸與起落架、回轉(zhuǎn)體之間設(shè)為轉(zhuǎn)動(dòng)副;鋼架結(jié)構(gòu)采用綁定設(shè)置。邊界條件的設(shè)定如圖3所示,在不考慮車體彈性變形及自身重力的情況下,將座圈與大地進(jìn)行固連。

圖4(a)、(b)分別給出了后坐位移為10 mm、15 mm時(shí)回轉(zhuǎn)座圈在彈性變形下的有限元計(jì)算應(yīng)力云圖。通過(guò)圖4結(jié)果對(duì)比發(fā)現(xiàn),位移為10 mm時(shí),座圈應(yīng)力較大位置滿足材料的強(qiáng)度要求;位移為15 mm時(shí),模型中很多部位的應(yīng)力已大于400 MPa,局部焊接點(diǎn)應(yīng)力近700 MPa,易發(fā)生危險(xiǎn)。因此,本文將緩沖行程設(shè)計(jì)為10 mm。

當(dāng)考慮車體與起落架變形時(shí),允許有更大的位移,可以適當(dāng)增大緩沖行程。由于具體的變形量值難以獲得,本文將緩沖行程的取值拓展到15 mm、20 mm和30 mm。

2 緩沖裝置設(shè)計(jì)

緩沖裝置結(jié)構(gòu)如圖5所示。圖5中,S和P分別表示兩腔的活塞面積與壓強(qiáng),D為活塞直徑。工作腔通過(guò)管道連接節(jié)流閥。在后坐力F的作用下,推動(dòng)活塞向后運(yùn)動(dòng),液體在壓力作用下流經(jīng)節(jié)流孔產(chǎn)生阻力。通過(guò)對(duì)節(jié)流孔開(kāi)口面積進(jìn)行控制,可實(shí)現(xiàn)對(duì)系統(tǒng)流量的控制,從而使后坐運(yùn)動(dòng)按設(shè)計(jì)規(guī)律進(jìn)行。為使后坐運(yùn)動(dòng)行程達(dá)到設(shè)計(jì)要求,對(duì)后坐運(yùn)動(dòng)軌跡進(jìn)行如下設(shè)計(jì)。

2.1 后坐運(yùn)動(dòng)軌跡設(shè)計(jì)

將緩沖裝置的后坐運(yùn)動(dòng)方程假設(shè)為

(1)

式中:a為加速度,v為速度,x為行程,t為緩沖時(shí)間,A、B、C為未知常數(shù)。

緩沖裝置所受負(fù)載作用力與推進(jìn)力大小相等,方向相反,其加載曲線如圖6所示。根據(jù)圖6中的負(fù)載作用時(shí)間,將緩沖裝置作用時(shí)間設(shè)計(jì)為200 ms。認(rèn)為作用結(jié)束時(shí)無(wú)沖擊,即后坐部分速度為0。在緩沖行程確定后,得到兩組后坐運(yùn)動(dòng)方程邊界條件,如表1所示。

表1 后坐運(yùn)動(dòng)方程邊界條件

將表1中的邊界條件代入式(1),得到不同緩沖行程下的后坐運(yùn)動(dòng)方程。以緩沖行程為10 mm為例,求解得到后坐運(yùn)動(dòng)方程為

(2)

根據(jù)后坐運(yùn)動(dòng)方程,得到緩沖裝置運(yùn)動(dòng)的位移、速度曲線,如圖7所示。

2.2 節(jié)流閥選取與節(jié)流孔開(kāi)口面積計(jì)算

為使緩沖裝置的后坐運(yùn)動(dòng)按照設(shè)計(jì)的軌跡進(jìn)行,需要通過(guò)節(jié)流孔開(kāi)口面積的變化對(duì)系統(tǒng)流量進(jìn)行調(diào)節(jié),從而控制缸內(nèi)壓力變化。通過(guò)液壓系統(tǒng)力的作用關(guān)系,應(yīng)用牛頓第二定律可求解出節(jié)流孔面積變化曲線。

(1)節(jié)流閥的選取。

不同結(jié)構(gòu)的節(jié)流孔,其節(jié)流效果差異明顯,如圖8所示的細(xì)長(zhǎng)孔和薄壁孔(孔長(zhǎng)lk通常為孔徑dk的4倍以上,即lk≥4dk)是兩種常見(jiàn)的典型節(jié)流孔結(jié)構(gòu)。根據(jù)不同的節(jié)流孔結(jié)構(gòu)調(diào)節(jié)流閥,通過(guò)改變節(jié)流面積或節(jié)流長(zhǎng)度來(lái)控制流體流量[9]。

對(duì)于細(xì)長(zhǎng)孔式節(jié)流閥,其節(jié)流孔流量計(jì)算公式為

(3)

式中:qv為節(jié)流孔流量,dk為細(xì)長(zhǎng)孔直徑,Δpk為孔進(jìn)出口壓差,μ為動(dòng)力粘度。

由于式(3)中含動(dòng)力粘度μ,在其他條件不變的情況下,節(jié)流孔流量受溫度影響較大,且隨細(xì)長(zhǎng)孔進(jìn)出口壓差呈線性變化,說(shuō)明節(jié)流孔的流量隨負(fù)載變化較大,在液壓調(diào)速系統(tǒng)中,體現(xiàn)為執(zhí)行端速度變化靈敏,穩(wěn)定性相對(duì)較差[10]。

對(duì)于薄壁孔式節(jié)流閥,其節(jié)流孔流量計(jì)算公式為

(4)

式中:Cd為流量系數(shù),S為小孔等效面積,ρ為液壓油密度。

與細(xì)長(zhǎng)孔情況不同的是,薄壁孔的流量與進(jìn)出口壓差的平方根成正比,而與動(dòng)力粘度無(wú)直接關(guān)系,說(shuō)明薄壁孔穩(wěn)定性較好,在液壓調(diào)速系統(tǒng)中,體現(xiàn)為執(zhí)行端速度穩(wěn)定。因此,本文選用薄壁孔電磁節(jié)流閥。

(2)節(jié)流孔開(kāi)口面積計(jì)算。

流量系數(shù)Cd的近似計(jì)算公式[11]為

(5)

式中:Cdmax根據(jù)經(jīng)驗(yàn)取0.7;λ為流數(shù);λcrit為臨界流數(shù),即流量系數(shù)達(dá)到常數(shù)時(shí)的λ值,一般取為1 000;μ為流體的動(dòng)力粘度。

節(jié)流孔流量與后坐速度的關(guān)系為

qv=vS1

(6)

液壓系統(tǒng)力的作用關(guān)系(不考慮沿程損失與局部損失)為

(7)

式中:F為負(fù)載;Ft為液壓推力;Ff為密封圈產(chǎn)生的密封阻力[12];M為后坐部分總質(zhì)量;a為后坐運(yùn)動(dòng)的加速度。

利用MATLAB編寫(xiě)計(jì)算程序,對(duì)由式(4)-(7)組成的方程組進(jìn)行求解,得到節(jié)流孔開(kāi)口面積曲線。以10 mm緩沖行程為例,求解得到節(jié)流孔開(kāi)口面積隨時(shí)間的變化曲線,如圖9所示。

根據(jù)后坐過(guò)程中液壓系統(tǒng)各力的作用關(guān)系,應(yīng)用牛頓第二定律求解節(jié)流孔開(kāi)口面積變化曲線。通過(guò)控制開(kāi)口面積,可實(shí)現(xiàn)對(duì)緩沖裝置系統(tǒng)流量的控制,從而使其后坐運(yùn)動(dòng)按照設(shè)計(jì)軌跡運(yùn)行,在確保緩沖行程的情況下實(shí)現(xiàn)緩沖。

對(duì)緩沖裝置各部分構(gòu)件進(jìn)行設(shè)計(jì)及強(qiáng)度校核后,設(shè)計(jì)如圖10所示的緩沖裝置。

當(dāng)緩沖行程為30 mm時(shí),末端反向加速度最大,所需提供的液壓推力最大。按照負(fù)載均勻分配到兩個(gè)裝置估算,單個(gè)液壓缸所提供的最大推力約為350 kN,經(jīng)計(jì)算設(shè)定最高工作壓力為20 MPa。根據(jù)流速與流量的關(guān)系,計(jì)算出管道內(nèi)最大流量為321.9 L/min,管道內(nèi)徑33.7 mm,圓整為40 mm,重新選取推薦流量為400 L/min。假設(shè)液壓油不可壓縮,整個(gè)液壓系統(tǒng)的任何部位流量相同,由設(shè)計(jì)手冊(cè)選取最大流速為6 m/s。

3 緩沖裝置液固聯(lián)合仿真

為驗(yàn)證所設(shè)計(jì)的緩沖裝置是否達(dá)到設(shè)計(jì)要求,利用剛體動(dòng)力學(xué)軟件Adams和復(fù)雜系統(tǒng)仿真平臺(tái)AMEsim進(jìn)行液固聯(lián)合仿真。Adams的坐標(biāo)系定義和AMEsim仿真模型分別如圖11(a)、(b)所示。

圖11(b)Adams模塊13中,f1和f2為通過(guò)AMEsim計(jì)算出的液壓推力的輸出信號(hào)接口。在Adams中計(jì)算出的位移和速度信號(hào)通過(guò)V1、V2、X1和X2輸出到AMEsim中進(jìn)行實(shí)時(shí)交互數(shù)據(jù)傳遞。ct1與ct2分別為Adams計(jì)算出的兩個(gè)基座的接觸力,a為活塞加速度,jiao1x、jiao1y、jiao1z、jiao2x、jiao2y和jiao2z為對(duì)應(yīng)基座的3個(gè)方向位移。

在固液聯(lián)合仿真過(guò)程中,Adams中計(jì)算出的活塞運(yùn)動(dòng)的速度、位移作為輸出傳遞到AMEsim中,AMEsim通過(guò)速度位移計(jì)算出液壓系統(tǒng)的壓強(qiáng)、流量等參數(shù),并將壓強(qiáng)產(chǎn)生的力作為輸出傳遞到Adams中,從而實(shí)現(xiàn)了數(shù)據(jù)實(shí)時(shí)傳遞[13]。

4 仿真結(jié)果分析

4.1 10 mm緩沖行程仿真結(jié)果

圖12(a)、(b)分別給出了行程為10 mm時(shí),導(dǎo)軌的后坐位移和速度的設(shè)計(jì)曲線與軟件仿真曲線。由圖12(a)可知,在導(dǎo)軌后坐過(guò)程的中間階段,仿真得到的位移比設(shè)計(jì)位移大0.15 mm左右。在仿真結(jié)果中,總的后坐位移為10.01 mm,誤差為0.1%。說(shuō)明緩沖裝置能夠較好滿足緩沖行程設(shè)計(jì)要求。由圖12(b)后坐速度的結(jié)果對(duì)比發(fā)現(xiàn),仿真結(jié)果在設(shè)計(jì)曲線上下波動(dòng)。這是由于在仿真中,基座與地面的連接關(guān)系設(shè)置為接觸,以監(jiān)測(cè)基座與地面之間的受力。而在本文研究中,導(dǎo)軌的受力狀態(tài)是時(shí)刻改變的,因接觸力而產(chǎn)生的速度波動(dòng)很難消除。但仍可看出,后坐速度的變化趨勢(shì)和幅值與設(shè)計(jì)值較吻合。整個(gè)后坐過(guò)程中,后坐速度的幅值為90 mm/s,負(fù)載作用結(jié)束時(shí)后坐速度為0,說(shuō)明本文設(shè)計(jì)能夠達(dá)到較好的緩沖效果。

液壓缸工作腔壓強(qiáng)隨時(shí)間變化曲線的仿真結(jié)果與設(shè)計(jì)計(jì)算結(jié)果如圖13所示。由圖13可以看出,仿真結(jié)果圍繞計(jì)算結(jié)果上下振蕩,這是因?yàn)橛?jì)算曲線為理想狀態(tài),而在仿真計(jì)算中,前100 ms負(fù)載線性增加,由節(jié)流孔開(kāi)口面積變化引起的缸內(nèi)壓力變化來(lái)不及與負(fù)載時(shí)時(shí)匹配,因而工作腔壓強(qiáng)出現(xiàn)較大起伏。而在后100 ms負(fù)載維持恒定,此后開(kāi)口面積變化相對(duì)平穩(wěn),缸內(nèi)壓強(qiáng)變化較小,總體壓強(qiáng)變化趨勢(shì)與設(shè)計(jì)曲線相符合。設(shè)計(jì)曲線的最大壓強(qiáng)值為15.90 MPa,仿真曲線得到的最大壓強(qiáng)值為14.97 MPa,均小于缸體所能承受的最大壓強(qiáng)20 MPa,滿足缸內(nèi)壓力設(shè)計(jì)要求。

緩沖裝置與地面的接觸力、液壓推力和負(fù)載隨時(shí)間的變化曲線如圖14所示。

由圖14可以看出,在緩沖行程結(jié)束時(shí),緩沖裝置的液壓推力大于負(fù)載作用力,說(shuō)明液壓缸能夠承受后坐力且具有反向運(yùn)動(dòng)的趨勢(shì)。在炮身的彈性變形恢復(fù)過(guò)程中,活塞將反向運(yùn)動(dòng)。

4.2 4種緩沖行程結(jié)果對(duì)比分析

考慮到輪胎、車體、發(fā)射裝置及千斤頂?shù)炔课豢赡馨l(fā)生的變形,設(shè)計(jì)時(shí)可選取更大的緩沖行程來(lái)滿足結(jié)構(gòu)強(qiáng)度要求,以保證設(shè)計(jì)的穩(wěn)定性與安全性。除上文10 mm緩沖行程外,本文還設(shè)計(jì)了15 mm、20 mm與30 mm 3種緩沖行程,將4種緩沖行程下的計(jì)算結(jié)果進(jìn)行對(duì)比分析。

圖15給出了4種緩沖行程下節(jié)流孔開(kāi)口面積隨時(shí)間的變化。由圖15對(duì)比發(fā)現(xiàn),緩沖行程越大,節(jié)流孔開(kāi)口面積也越大,且面積的變化具有一致的規(guī)律。4種緩沖行程的最大開(kāi)口面積分別為14.07 mm2、21.20 mm2、28.193 mm2、42.42 mm2。當(dāng)其他條件不變時(shí),緩沖行程每增加5 mm,最大開(kāi)口面積增大約7 mm2,且最大開(kāi)口面積的比值與緩沖行程的比值相同,所有開(kāi)口面積均小于節(jié)流孔最大開(kāi)口面積。

4種緩沖行程下的后坐位移對(duì)比如圖16所示。圖16顯示,10 mm、15 mm、20 mm和30 mm 4種緩沖行程的最終行程分別為10.01 mm、15.01 mm、20.05 mm和30.20 mm,誤差分別為0.1%、0.067%、0.25%和0.67%。所設(shè)計(jì)行程的結(jié)果誤差均在1%的可接受范圍。

4種緩沖行程下的工作腔壓強(qiáng)對(duì)比如圖17所示。從圖17可以看出,緩沖行程越大時(shí),缸內(nèi)壓強(qiáng)越平穩(wěn),4種緩沖行程下的最終壓強(qiáng)分別為14.97 MPa、15.26 MPa、15.64 MPa和16.44 MPa,均小于缸體所能承受的最大設(shè)計(jì)壓強(qiáng)20 MPa。

4種緩沖行程下的節(jié)流孔流量對(duì)比如圖18所示。由圖18可見(jiàn),節(jié)流孔流量以相同規(guī)律成遞增趨勢(shì),緩沖行程越大,流量越平穩(wěn)。在緩沖行程為30 mm時(shí),節(jié)流孔流量最大(325.1 L/min),小于管道的推薦流量400 L/min,滿足設(shè)計(jì)要求。

5 結(jié)論

針對(duì)新型電磁彈射火箭炮帶來(lái)的巨大后坐力問(wèn)題,本文設(shè)計(jì)了一種小行程液壓緩沖裝置,通過(guò)仿真計(jì)算與對(duì)比分析,得到如下主要結(jié)論:

(1)對(duì)10 mm緩沖行程的設(shè)計(jì)計(jì)算與仿真結(jié)果的對(duì)比分析發(fā)現(xiàn),本文設(shè)計(jì)的緩沖裝置的緩沖行程、液壓缸壓強(qiáng)均滿足設(shè)計(jì)要求;

(2)緩沖行程為10 mm時(shí),整個(gè)后坐過(guò)程最大后坐速度不超過(guò)90 mm/s,緩沖效果較好;

(3)考慮車體變形,還設(shè)計(jì)了15 mm、20 mm和30 mm 3種較大的緩沖行程。對(duì)4種緩沖行程下節(jié)流孔開(kāi)口面積、后坐位移、工作腔壓強(qiáng)和節(jié)流孔流量等進(jìn)行了對(duì)比,表明不同緩沖行程時(shí)各結(jié)構(gòu)均可滿足要求。緩沖行程越大,系統(tǒng)穩(wěn)定性越好,但同時(shí)對(duì)系統(tǒng)結(jié)構(gòu)強(qiáng)度要求更高。