制導(dǎo)炮彈膛內(nèi)轉(zhuǎn)速特性建模及仿真研究

齊竹昌，張鵬飛，柴 勁，張 意

(西安現(xiàn)代控制技術(shù)研究所，西安 710065)

0 引言

火炮以射速快、射程遠、威力大的優(yōu)勢在以往的戰(zhàn)爭中發(fā)揮著巨大的作用，而隨著制導(dǎo)技術(shù)的發(fā)展，高精度制導(dǎo)彈藥在現(xiàn)代化戰(zhàn)爭中已經(jīng)占據(jù)了主導(dǎo)作用，常規(guī)炮彈已經(jīng)無法滿足當(dāng)前戰(zhàn)爭的需要，制導(dǎo)炮彈應(yīng)運而生。目前世界上主流的制導(dǎo)炮彈為激光末制導(dǎo)炮彈、衛(wèi)星制導(dǎo)炮彈及二維修正制導(dǎo)炮彈[1]。與采用高速旋轉(zhuǎn)飛行穩(wěn)定的榴彈不同，制導(dǎo)炮彈通常采用折疊式尾翼和鴨式氣動布局[2]，為保證制導(dǎo)炮彈有足夠的過載能力，翼展通常設(shè)計較大，尾翼在火炮膛內(nèi)處于折疊狀態(tài)，出炮口后尾翼張開，彈體成為靜穩(wěn)定狀態(tài)。

制導(dǎo)炮彈出炮口時轉(zhuǎn)速過高會導(dǎo)致尾翼張開過程中離心力過大，翼片張開角速度增加，使翼片無法正常鎖定；同時翼片與尾翼基座發(fā)生高速碰撞，會使尾翼結(jié)構(gòu)受損，影響飛行彈道。為了避免尾翼在炮口高速旋轉(zhuǎn)時張開結(jié)構(gòu)受損，制導(dǎo)炮彈多采用金屬滑動彈帶設(shè)計，以減小制導(dǎo)炮彈膛內(nèi)及出炮口的轉(zhuǎn)速，制導(dǎo)炮彈炮口轉(zhuǎn)速通常設(shè)計到6～10 r/s。

慈明森[3]針對滑動彈帶彈丸炮口轉(zhuǎn)速的計算方法進行了研究，推導(dǎo)了炮口轉(zhuǎn)速計算模型，該方法未引入火炮內(nèi)彈道模型，僅能在已知膛壓條件下求出炮口轉(zhuǎn)速，無法得到制導(dǎo)炮彈在膛內(nèi)運動過程中的轉(zhuǎn)速變化情況，同時其采用膛底壓力的處理方法與真實結(jié)果存在一定偏差。文中針對制導(dǎo)炮彈膛內(nèi)運動過程，通過建立火炮內(nèi)彈道數(shù)學(xué)仿真模型和彈丸在膛內(nèi)的受力分析模型，對制導(dǎo)炮彈在火炮膛內(nèi)運動過程中的轉(zhuǎn)速變化進行仿真計算，分析了不同因素對制導(dǎo)炮彈膛內(nèi)轉(zhuǎn)速的影響，同時通過與某型制導(dǎo)炮彈膛內(nèi)轉(zhuǎn)動過程實測結(jié)果對比，對計算模型的正確性進行驗證，為制導(dǎo)炮彈膛內(nèi)轉(zhuǎn)速設(shè)計提供參考。

1 彈丸火炮膛內(nèi)建模

1.1 火炮零維內(nèi)彈道模型

制導(dǎo)炮彈在火炮發(fā)射過程中存在前期、熱力學(xué)第一時期、熱力學(xué)第二時期和后效期4個階段[4]。在前期點火發(fā)生后，火炮發(fā)射裝藥開始燃燒，膛內(nèi)壓力增加，但此時壓力小于彈丸擠進壓力P0，彈丸在膛內(nèi)不運動。從膛內(nèi)壓力上升至P0，彈丸開始運動，至裝藥燃燒結(jié)束時刻點tk為熱力學(xué)第一時期。從火藥燃燒結(jié)束點tk開始，一直持續(xù)到彈丸運動到彈底與炮口重合時刻te為熱力學(xué)第二時期。從te開始一直到平均彈道壓力等于臨界壓力Pcr時結(jié)束為后效期，對于火藥氣體流出到空氣中的情況，臨界壓力Pcr約為0.18 MPa。

根據(jù)發(fā)射裝藥在膛內(nèi)燃燒過程和彈丸運動過程，建立火炮零維內(nèi)彈道方程[4-5]為：

(1)

式中：zi為第i種火藥燃燒過程中的相對厚度，zi≤1表示火藥分裂前，zi>1表示火藥分裂后;Iki表示第i種火藥壓力全沖量;P為火炮膛內(nèi)平均壓力;ν為火藥壓力燃速指數(shù);ψi為第i種火藥燃燒過程中生成相對氣體質(zhì)量;χ1i、χ2i、λ1i、λ2i、μ1i為火藥藥粒分裂前后的與相對弧厚、相對周長及相對燃燒面積形狀函數(shù)系數(shù);v為彈丸在膛內(nèi)的運動速度;l為彈丸在膛內(nèi)運動距離;m為彈丸重量;S為炮管截面積;Pd為彈底壓力;Pt為膛底壓力;ωm為裝藥總質(zhì)量;ωm,i為第i種火藥質(zhì)量;θ為氣體比熱比;Kq為平均熱損失系數(shù);φ1為次要功系數(shù);W0為藥室初始容積;αi第i種火藥余容;δi為第i種火藥密度。

1.2 彈丸膛內(nèi)受力分析建模

制導(dǎo)炮彈在發(fā)射過程中，火藥燃燒產(chǎn)生的高溫高壓氣體推動彈丸沿著身管運動。由于制導(dǎo)炮彈控制系統(tǒng)響應(yīng)頻率的限制，制導(dǎo)炮彈基本采用滑動彈帶使彈丸在膛內(nèi)實現(xiàn)閉氣和減旋。當(dāng)膛內(nèi)壓力大于彈帶擠進壓力P0時，彈丸開始運動，滑動彈帶嵌入膛線槽內(nèi)，沿膛線高速旋轉(zhuǎn)，彈丸本體在滑動彈帶摩擦接觸面的帶動下低速旋轉(zhuǎn)。

為研究膛線對滑動彈帶的作用力，將炮膛表面展開成Oxy平面，圖1(a)為真實火炮膛線，圖1(b)為彈丸在膛內(nèi)的受力關(guān)系圖，其中Ox軸與炮膛軸線平行，Oy軸為展開后的徑向方向，OO曲線為漸速膛線，膛線纏繞角為α，A點為導(dǎo)轉(zhuǎn)側(cè)中心，在該點膛線對滑動彈帶的作用力為正壓力N和摩擦力μN，μ為膛動彈帶與膛線的摩擦力系數(shù)，將該力沿Ox，Oy軸方向進行分解，則沿Ox軸方向的分力R和沿Oy軸方向的分力T可表示為：

圖1 膛線展開平面彈丸受力分析圖

(2)

彈丸在膛內(nèi)運動過程中，滑動彈帶與彈體相互作用力存在兩種狀態(tài)，第一種狀態(tài)為滑動彈帶與壓螺接觸，如圖2(a)所示，此時彈體推動滑動彈帶沿炮膛軸線向前運動，第二種狀態(tài)為滑動彈帶與彈體接觸，如圖2(b)所示，此時滑運彈帶推動彈體沿炮膛軸線向前運動。因為滑動彈帶與彈體間隙約1～2 mm，滑動彈帶運動速度可視為與彈體相同。

圖2 滑動彈帶與彈體相互作用力關(guān)系圖

若滑動彈帶與彈體相互作用力為第一種狀態(tài)，滑動彈帶沿膛線的直線運動、旋轉(zhuǎn)運動和彈體的旋動運動可表示為方程(3)所示，滑動彈帶與壓螺作用力為Q1。

(3)

若滑運彈帶與彈體相互作用力為第二種狀態(tài)，滑動彈帶沿膛線的直線運動、旋轉(zhuǎn)運動和彈體的旋動運動可表示為方程(4)所示，滑動彈帶與彈體作用力為Q2。

(4)

式中：mb為滑動彈帶的質(zhì)量;Jxb為滑動彈帶繞彈軸的轉(zhuǎn)動慣量;ωb為滑動彈帶環(huán)繞彈軸的轉(zhuǎn)動角速度;n為膛線數(shù);rcp為考慮膛線時炮膛端面的換算半徑;rh為壓螺的外半徑;ra為滑動彈帶內(nèi)端面半徑;r為彈丸半徑;f1為彈帶環(huán)后端面與壓螺之間的摩擦力系數(shù);f2為彈帶環(huán)前端面與彈體之間的摩擦力系數(shù);Jx為彈丸繞縱軸的轉(zhuǎn)動慣量;ω為彈丸在膛內(nèi)的轉(zhuǎn)動角速度。

當(dāng)Q1>0，Q2<0時，表明滑動彈帶后端面與壓螺接觸，滑動彈帶后端面通過摩擦帶動彈體轉(zhuǎn)動；當(dāng)Q1<0，Q2>0時，表明滑動彈帶前端面與彈體接觸，滑動彈帶前端面通過摩擦帶動彈體轉(zhuǎn)動。

2 彈丸膛內(nèi)轉(zhuǎn)速影響因素分析

2.1 彈丸膛內(nèi)轉(zhuǎn)速求解

通過聯(lián)合求解式(1)、式(3)和式(1)、式(4)可求解出火炮發(fā)射過程中制導(dǎo)炮彈彈丸的運動速度、運動距離及轉(zhuǎn)速。

如今，因為擔(dān)心夫人摔著，潘際鑾院士不再騎車載她了，而是改為手牽手，去菜場、去學(xué)校、去訪學(xué)……時時刻刻不分離。潘際鑾是南昌大學(xué)的老校長，南昌大學(xué)有一條“際鑾路”，老兩口回南昌大學(xué)故地重游，牽手走在這條路上，也被學(xué)生拍到了。

若滑動彈帶與壓螺接觸，則膛線作用在滑動彈帶上的正壓力N、滑動彈帶與壓螺作用力Q1及彈丸轉(zhuǎn)速ω可表示為：

(5)

(6)

(7)

若滑動彈帶與彈體接觸，則膛線作用在滑動彈帶上的正壓力N、滑動彈帶與彈體作用力Q2及彈丸轉(zhuǎn)速ω可表示為：

(8)

(9)

(10)

某122 mm榴彈炮制導(dǎo)炮彈火炮內(nèi)彈道計算參數(shù)見表1，彈體與滑動彈帶幾何物理參數(shù)見表2，對其內(nèi)彈道及膛內(nèi)轉(zhuǎn)速進行仿真計算，膛壓、炮口速度和膛內(nèi)轉(zhuǎn)速仿真結(jié)果與火炮試驗實測結(jié)果對比分別見圖3～圖5，仿真結(jié)果與實測結(jié)果基本一致，表明火炮內(nèi)彈道模型和膛內(nèi)受力分析模型正確。彈丸在膛內(nèi)運動過程中，滑動彈帶前后接觸力如圖6所示，Q1<0，Q2>0，表明彈丸在膛內(nèi)運動過程中，滑動彈帶前端面始終與彈體接觸，滑動彈帶前端面通過摩擦帶動彈體轉(zhuǎn)動。

圖3 仿真膛壓與實測膛底壓力對比圖

圖4 仿真膛內(nèi)速度與實測初速對比圖

圖5 仿真彈丸膛內(nèi)轉(zhuǎn)速與實測轉(zhuǎn)速對比圖

圖6 滑動彈帶前后接觸力

表1 某122 mm火炮內(nèi)彈道計算原始參數(shù)表

表2 某122 mm制導(dǎo)炮彈彈體與滑動彈帶幾何物理參數(shù)

2.2 發(fā)射裝藥號對膛內(nèi)轉(zhuǎn)速的影響

發(fā)射裝藥號不同，對應(yīng)的膛壓和彈丸膛內(nèi)運動速度不同，某122 mm榴彈炮1#、2#、3#裝藥的膛壓、速度和轉(zhuǎn)速曲線分別如圖7～圖9所示，隨著裝藥量減小，膛壓和初速降低，炮口轉(zhuǎn)速下降，下降規(guī)律與速度變化基本一致。

圖7 1#裝藥的膛壓、速度和轉(zhuǎn)速曲線

圖8 2#裝藥的膛壓、速度和轉(zhuǎn)速曲線

圖9 3#裝藥的膛壓、速度和轉(zhuǎn)速曲線

2.3 滑動彈帶摩擦系數(shù)對膛內(nèi)轉(zhuǎn)速的影響

由于彈丸發(fā)射時，滑動彈帶會嵌入膛線槽，滑動彈帶與膛線的摩擦力系數(shù)只與其材料有關(guān)，制導(dǎo)炮彈滑動彈帶大多采用紫銅材料，其與炮管的摩擦力系數(shù)基本固定，無法調(diào)整，而前后表面的摩擦力系數(shù)可通過增加不同的潤滑工藝進行調(diào)整，前后端面取不同摩擦力系數(shù)時，彈丸膛內(nèi)轉(zhuǎn)速變化情況如圖10所示。隨著摩擦力系數(shù)的增大，彈丸在膛內(nèi)的轉(zhuǎn)速急劇增大。

圖10 不同摩擦力系數(shù)對轉(zhuǎn)速影響

2.4 滑動彈帶幾何參數(shù)對膛內(nèi)轉(zhuǎn)速的影響

滑動彈帶外端半徑rcp和彈丸半徑r均與火炮身管尺寸相關(guān)，尺寸參數(shù)無法調(diào)整，彈丸運動過程中，只有滑動彈帶前端面與彈體接觸，壓螺外半徑rh不影響膛內(nèi)轉(zhuǎn)速。滑動彈帶內(nèi)端半徑ra的變化會明顯改變滑動彈帶與彈體的摩擦力矩，影響轉(zhuǎn)速變化。當(dāng)ra取不同值時，彈體轉(zhuǎn)速變化情況如圖11所示，隨著半徑的增大，摩擦力矩增大，轉(zhuǎn)速增大，但轉(zhuǎn)速增加量值較小，當(dāng)滑動彈帶內(nèi)端半徑增大13 mm時，轉(zhuǎn)速僅增加0.9 r/s。

圖11 不同滑動彈帶內(nèi)端半徑對轉(zhuǎn)速影響

3 彈丸膛內(nèi)轉(zhuǎn)速改善方法

通過對彈丸膛內(nèi)轉(zhuǎn)速影響因素分析，影響轉(zhuǎn)速變化的主要因素為滑動彈帶摩擦力系數(shù)，發(fā)射裝藥和滑動彈帶幾何參數(shù)對膛內(nèi)轉(zhuǎn)速的影響均較小。因此，要改善彈丸炮口轉(zhuǎn)速，可通過調(diào)整滑動彈帶摩擦力系數(shù)實現(xiàn)。

目前常用的潤滑方式分為流體潤滑和固體潤滑，流體潤滑以潤滑油、潤滑脂等液體材料作為潤滑劑，但是流體潤滑對溫度非常敏感，如油脂類潤滑材料的最高使用溫度不超過200 ℃，溫度過高將導(dǎo)致潤滑劑因氧化變質(zhì)而失效。固體潤滑劑種類有很多，如鋅、鉛等軟金屬，二硫化鉬、硫化亞鐵等金屬化合物，以及鎳基合金、銅基合金等合金材料[6-8]。其中，MoS2對高溫、低溫、高負載、高速等條件下的設(shè)備有優(yōu)異的潤滑功效。文獻[9]研究發(fā)現(xiàn)，MoS2涂層在全接觸運動條件下，摩擦力系數(shù)基本在0.067～0.081。根據(jù)前述仿真計算結(jié)果，該摩擦力系數(shù)符合炮口轉(zhuǎn)速6～9 r/s的需求。此外，MoS2還具有在火炮膛內(nèi)高溫高壓的氮氧燃氣環(huán)境下不易分解、經(jīng)濟性好等優(yōu)點。因此MoS2十分適合作為炮彈閉氣環(huán)的潤滑材料。

4 結(jié)論

結(jié)合火炮零維內(nèi)彈道模型，完成了彈丸膛內(nèi)運動受力建模分析，對制導(dǎo)炮彈膛內(nèi)轉(zhuǎn)速特性進行了研究，通過與試驗結(jié)果對比，驗證了模型的正確性。分析了發(fā)射裝藥、滑動彈帶摩擦力系數(shù)和滑動彈帶幾何參數(shù)參對彈丸膛內(nèi)轉(zhuǎn)速的影響，提出了改善制導(dǎo)炮彈膛內(nèi)轉(zhuǎn)速的方法，為制導(dǎo)炮彈膛內(nèi)轉(zhuǎn)速設(shè)計提供了參考依據(jù)。