李　鋒,拜云山,馮曉偉

(中國工程物理研究院 總體工程研究所,四川 綿陽 621999)

?

一級氣體炮內(nèi)彈道方程修正及驗證

李鋒,拜云山,馮曉偉

(中國工程物理研究院 總體工程研究所,四川 綿陽 621999)

為了使理論計算結(jié)果更為準(zhǔn)確可靠,對一級氣體炮內(nèi)彈道經(jīng)典方程進(jìn)行修正。分析了一級氣體炮內(nèi)彈道經(jīng)典方程存在的問題及原因,指出氣室容積不能忽略。提出了氣體密度新的表達(dá)方式,并根據(jù)質(zhì)量守恒方程、氣流動量方程、彈丸運動方程、氣體絕熱膨脹方程和理想氣體狀態(tài)方程等,聯(lián)立推導(dǎo)出一級氣體炮內(nèi)彈道修正方程。將修正方程、經(jīng)典方程計算結(jié)果與內(nèi)彈道流體數(shù)值仿真結(jié)果進(jìn)行對比，結(jié)果表明:在容積比較小的情況下三者的偏差很小;但隨著容積比的增大,經(jīng)典方程計算結(jié)果偏離修正方程和數(shù)值仿真計算結(jié)果;經(jīng)典方程僅適用于容積比較小的情況;而修正方程的適用范圍不受容積比大小限制,驗證了修正方程的正確性和適用性。

氣體炮;內(nèi)彈道;炮口速度;容積比

氣體炮作為一種加載試驗裝置,采用壓縮氣體作為動力源驅(qū)動彈丸運動以獲得期望的速度。1946年,美國研制成功第一門利用輕質(zhì)氣體作為工作介質(zhì)的一級氣體炮[1],此后,美國、加拿大、英國和法國以及我國一些單位陸續(xù)建立了氣體炮裝置,使這種類型的彈丸發(fā)射裝置成為一種使用極其廣泛的試驗設(shè)備。

Pidduck和Kent最早進(jìn)行了氣體炮內(nèi)彈道方程的研究,提出采用振蕩函數(shù)近似描述波在彈后的傳播規(guī)律[2],給出了解析模型。Seigel A E采用數(shù)值仿真方法對一級氣體炮的氣體動力學(xué)問題進(jìn)行了研究,并與Pidduck和Kent的解析模型計算結(jié)果進(jìn)行對比,指出當(dāng)壓縮氣體與彈丸質(zhì)量之比小于1/4時,數(shù)值計算結(jié)果與Pidduck的解析解非常接近[3]。王金貴借用火炮內(nèi)彈道理論中虛擬質(zhì)量系數(shù)表達(dá)形式,結(jié)合氣體動力學(xué)方程推導(dǎo)了氣體炮內(nèi)彈道方程[4],獲得了炮口初速的解析表達(dá)式,并在氣體炮設(shè)計上得到了應(yīng)用,實際上已成為國內(nèi)氣體炮設(shè)計的理論依據(jù),本文稱之為經(jīng)典方程。趙俊利等[5]建立了氣體炮內(nèi)彈道數(shù)學(xué)模型,給出了與文獻(xiàn)[4]相同的內(nèi)彈道理論公式,在具體參數(shù)計算方法上進(jìn)行了簡化。Mark Denny[6]將彈丸附加質(zhì)量作為常量處理,建立了考慮彈前氣壓和破膜能量損失影響的內(nèi)彈道解析模型,該模型在去掉彈前氣壓和破膜能量損失項后在形式上與經(jīng)典方程相同。

目前,研究人員在氣體炮內(nèi)彈道方面多針對具體結(jié)構(gòu)型式進(jìn)行數(shù)值仿真。楊均勻等[7]采用Von Neumann Richtmyer差分格式對帶有坡膛、有限藥室長的一級輕氣炮進(jìn)行了內(nèi)彈道一維數(shù)值模擬。仲偉君等[8]針對四腔結(jié)構(gòu)的氣體炮,運用氣體動力學(xué)的相似理論建立了氣體炮的內(nèi)彈道模型,并在計算機(jī)上進(jìn)行了數(shù)值模擬。趙希欣等[9]應(yīng)用氣體動力學(xué)理論及配合間隙的影響建立了考慮氣體泄露的內(nèi)彈道模型,并進(jìn)行了數(shù)值仿真分析。胡靜等[10]建立了三維非定常條件下35 mm等截面輕氣炮數(shù)學(xué)模型,采用ANSYS-CFX軟件對不同氣體工質(zhì)所獲彈丸初速進(jìn)行模擬計算。夏正友[11]采用一種基于推理的神經(jīng)網(wǎng)絡(luò)模型,對氣體炮內(nèi)彈道最大膛壓進(jìn)行了仿真和預(yù)測。

一級氣體炮經(jīng)典方程針對較為理想的結(jié)構(gòu)型式,借鑒火炮內(nèi)彈道[12]的相關(guān)理論,通過理論分析和推導(dǎo),獲得了炮口初速與內(nèi)彈道參數(shù)的解析表達(dá)式,與數(shù)值模擬相比在工程應(yīng)用上更為方便。但由于簡化較多,實際應(yīng)用中發(fā)現(xiàn)經(jīng)典方程在大氣室容積的情況下存在一定的局限性。本文分析了經(jīng)典內(nèi)彈道方程存在的問題及原因,對炮管內(nèi)氣體密度表達(dá)式進(jìn)行了修正,并推導(dǎo)了新的氣體炮內(nèi)彈道方程,稱之為修正方程,進(jìn)一步將修正方程和經(jīng)典方程與數(shù)值仿真結(jié)果進(jìn)行對比,給出了各自的適用范圍,完成了對修正方程的驗證。

1　氣體炮內(nèi)彈道經(jīng)典方程分析

根據(jù)一般的一級氣體炮結(jié)構(gòu)形式,建立氣體炮內(nèi)彈道物理模型,如圖1所示,圖中,S為炮管截面積,L為彈丸在膛內(nèi)的行程,lp為t時刻彈丸在膛內(nèi)的運動位移。采用如下基本假設(shè):

①氣室氣體為理想氣體,整個發(fā)射過程視為理想氣體絕熱膨脹過程;

②不考慮彈前氣體壓力的影響;

③不考慮各種損耗的具體形式,在計算過程中引入虛擬質(zhì)量系數(shù)φ來計入能量損失因素,φ的具體形式與火炮內(nèi)彈道方程相同[4]。

圖1　一級氣體炮內(nèi)彈道物理模型

王金貴[4]基于上述假設(shè),通過聯(lián)立求解理想氣體狀態(tài)方程、氣體絕熱過程方程、彈丸運動方程等,得到彈丸炮口初速與氣室容積、氣室壓力、炮彈質(zhì)量、炮管尺寸等的關(guān)系:

(1)

式中:vm為彈丸炮口初速,p0為氣室初始壓力,V0為氣室容積,m為彈丸質(zhì)量,φ為虛擬質(zhì)量系數(shù),c為絕熱指數(shù)。

本文稱式(1)為一級氣體炮內(nèi)彈道經(jīng)典方程。根據(jù)式(1),炮口初速隨氣室容積的變化規(guī)律如圖2所示(炮管直徑為80 mm,彈丸行程為1 m,彈丸質(zhì)量為2 kg,氣室壓力為5 MPa)。從圖2中可以看出,炮口初速隨著氣室容積的增大先增大后減小,而內(nèi)彈道數(shù)值仿真結(jié)果表明炮口初速隨著氣室容積增大而增大。

圖2　炮口初速隨氣室容積的變化規(guī)律對比

進(jìn)一步分析理論計算與數(shù)值仿真結(jié)果不一致的原因,對一級氣體炮內(nèi)彈道經(jīng)典方程的推導(dǎo)過程進(jìn)行梳理,發(fā)現(xiàn)經(jīng)典方程中氣體密度ρ的表達(dá)式為[12]

(2)

式中:m0為氣體質(zhì)量。

式(2)中將氣體密度ρ定義為氣體質(zhì)量與炮管容積之比,沒有考慮氣室容積的大小。當(dāng)氣室容積相對于炮管容積較大時,若忽略氣室容積就會導(dǎo)致氣體密度顯著變大,從而使得氣體的附加質(zhì)量效應(yīng)人為變大,導(dǎo)致炮口初速變小,甚至出現(xiàn)氣室容積越大炮口初速反而越小的情況。故分析認(rèn)為,當(dāng)氣室容積相對于炮管容積較大時,一級氣體炮內(nèi)彈道經(jīng)典方程的氣體密度表達(dá)式中忽略氣體容積是不合理的。

2　氣體炮內(nèi)彈道修正方程推導(dǎo)

根據(jù)前文的分析,氣體密度表達(dá)式應(yīng)考慮氣室容積,將氣體密度ρ修正為

(3)

基于修正的氣體密度表達(dá)式,對內(nèi)彈道方程進(jìn)行推導(dǎo),具體過程如下。

根據(jù)彈后氣體質(zhì)量守恒方程,并假設(shè)任一時刻彈底和膛底之間的氣體密度ρ為常量[3],則彈后氣體速度vx為

(4)

式中:v為彈丸速度,x為氣體微團(tuán)距炮管左端面的距離。

根據(jù)氣流動量方程,彈后氣體的壓力分布為

(5)

式中：pd為彈底壓力。則膛內(nèi)平均壓力p為

(6)

根據(jù)彈丸運動方程,并引入阻力系數(shù)φ1和動態(tài)虛擬質(zhì)量系數(shù)φ′,則彈丸的加速度為

(7)

由式(7)有:

(8)

將式(3)和式(5)代入式(6),并與式(8)聯(lián)立求解,可得虛擬質(zhì)量系數(shù)φ′的表達(dá)式:

(9)

氣體炮內(nèi)彈道經(jīng)典方程[4]中虛擬質(zhì)量系數(shù)φ表達(dá)式為

(10)

對比式(9)和式(10)可以看出,修正方程的虛擬質(zhì)量系數(shù)表達(dá)式與經(jīng)典方程明顯不同。

根據(jù)熱力學(xué)基本方程,平均膛壓為

(11)

將式(9)和式(11)代入式(7)并積分,可以得到:

(12)

式(12)即為一級氣體炮內(nèi)彈道修正方程。

3　氣體炮內(nèi)彈道修正方程驗證

為了驗證修正方程的正確性和適用性,將修正方程的計算結(jié)果與數(shù)值仿真結(jié)果進(jìn)行對比。數(shù)值仿真在流體仿真軟件Fluent平臺下進(jìn)行,采用了與修正方程相同的假設(shè),進(jìn)行了不同工作氣體(空氣、氦氣、氫氣)、不同氣室容積下的彈丸膛內(nèi)運動過程仿真。

具體仿真方法:采用二維軸對稱模型,使用結(jié)構(gòu)化網(wǎng)格對計算域進(jìn)行離散。由于彈丸發(fā)射過程中流場計算區(qū)域網(wǎng)格隨著彈丸的運動會不斷變化,數(shù)值仿真應(yīng)用彈簧近似法和動態(tài)分層法相結(jié)合的域動分層網(wǎng)格更新技術(shù)對計算模型進(jìn)行處理。彈丸受力通過Fluent求解器對彈表面壓力積分獲得(不考慮彈前壓力),加速度根據(jù)彈丸受到的氣體壓力由牛頓第二定律計算獲得,彈的速度和位置通過對加速度進(jìn)行積分求解獲得。

編寫自定義彈丸運動UDF程序嵌入Fluent求解器控制彈丸運動,采用密度基Couple算法對彈丸發(fā)射過程中的流場進(jìn)行計算,并得到彈丸在膛內(nèi)任意位置的速度及炮口初速。

主要仿真參數(shù):氣室容積分別為4.71 L,18.84 L,75.36 L;氣室壓力為5 MPa;彈丸質(zhì)量為2 kg;炮管直徑為80 mm;彈丸膛內(nèi)運動行程為1 m。

當(dāng)氣室容積為18.84 L,工作氣體為空氣時,彈丸出炮口時的流場壓力p分布如圖3所示。

圖3　彈丸出炮口時的流場壓力分布

不同工作氣體、不同容積比r情況下的數(shù)值仿真結(jié)果與修正方程計算結(jié)果、經(jīng)典方程計算結(jié)果的對比如表1和圖4所示,其中,e為修正方程與數(shù)值仿真結(jié)果炮口速度的偏差,容積比r為

(13)

圖4　炮口初速隨容積比的變化規(guī)律對比

工況2下彈丸速度與其在炮膛內(nèi)的運動行程關(guān)系如圖5所示。

圖5　彈丸速度隨其行程的變化規(guī)律(空氣)

從圖4、圖5可以看出:

①當(dāng)容積比較小時,經(jīng)典方程和修正方程的炮口初速曲線基本上是重合的,且與數(shù)值仿真結(jié)果吻合;但隨著容積比的增大,修正方程的炮口初速逐漸增大,與數(shù)值仿真結(jié)果吻合,而經(jīng)典方程的炮口初速是先增大后減小。

②氣體的摩爾質(zhì)量越小,經(jīng)典方程和修正方程的差別越小,這是因為氣體摩爾質(zhì)量小時式(9)和式(10)的差別小,引起的炮口初速偏差也小。

③經(jīng)典方程、修正方程和數(shù)值仿真得到的彈丸在膛內(nèi)運動速度的變化規(guī)律相同,但修正方程與數(shù)值仿真結(jié)果在數(shù)值上更為接近,說明修正方程的計算精度更高。

4　討論

一級氣體炮內(nèi)彈道經(jīng)典方程提出較早,且也在實踐中得到了應(yīng)用,從一些公開文獻(xiàn)[2]中看,在容積比較小的情況下,經(jīng)典方程的計算結(jié)果與實測結(jié)果的偏差并不大,說明經(jīng)典方程在容積比較小時是適用的,可以滿足工程設(shè)計精度的要求。但一級氣體炮內(nèi)彈道經(jīng)典方程的密度表達(dá)式忽略了氣室容積大小,導(dǎo)致在容積比較大時炮口初速的變化規(guī)律不正確,說明經(jīng)典方程不適用于容積比較大的情況。

由于氣體炮的氣室容積往往比較大,容積比也較大,若忽略氣室容積大小的影響,則可能引入較大的誤差甚至錯誤。因此修正方程基于氣室容積不能忽略的考慮,通過修正氣體密度表達(dá)式,重新推導(dǎo)出炮口初速的表達(dá)式。修正方程與經(jīng)典方程和數(shù)值仿真結(jié)果的對比表明:在容積比較小時,三者的計算結(jié)果基本一致;但在容積比較大時,修正方程的炮口初速變化規(guī)律與數(shù)值仿真結(jié)果一致,而與經(jīng)典方程明顯不同,進(jìn)一步驗證了修正方程的正確性和適用性。

5　結(jié)論

本文基于氣室容積不能忽略的考慮,對炮管內(nèi)氣體密度表達(dá)式進(jìn)行了修正,聯(lián)立推導(dǎo)出一級氣體炮內(nèi)彈道修正方程,并將修正方程與經(jīng)典方程和數(shù)值仿真結(jié)果進(jìn)行對比,分析結(jié)果表明,經(jīng)典方程僅適用于容積比較小的情況,而修正方程的適用范圍不受容積比大小的限制。修正方程在一定程度上更準(zhǔn)確地揭示了炮口初速隨氣室容積的變化規(guī)律,可以用于指導(dǎo)設(shè)計,并提高炮口初速在設(shè)計階段的估算精度。根據(jù)修正方程進(jìn)行理論設(shè)計,設(shè)計可行域較經(jīng)典方程大,可增加設(shè)計的靈活性。

[1]CROZIER W D,HUME W.High-velocity light-gas gun[J].Journal of Applied Physics,1957,28(8):892-894.

[2]BOURNE N K.A 50 mm bore gas gun for dynamic loading of materials and structures[J].Measurement Science and Technology,2003,14:273-278.

[3]SEIGEL A E.The theory of high speed guns，AD475660XAB[R].1965.

[4]王金貴.氣體炮原理及技術(shù)[M].北京:國防工業(yè)出版社,2001.

WANG Jin-gui.Principle and technology of gas gun[M].Beijing:National Defense Industry Press,2001.(in Chinese)

[5]趙俊利,曹鋒.氣體炮實用內(nèi)彈道方程及應(yīng)用[J].火炮發(fā)射與控制學(xué)報,2003(3):48-51.

ZHAO Jun-li,CAO Feng.Practical interior ballistic equation of the gas gun and its application[J].Journal of Gun Launch and Control,2003(3):48-51.(in Chinese)

[6]DENNY M.Gas gun dynamics[J].European Jouranal of Physics,2013,34:1 327-1 336.

[7]楊均勻,金志明,張小兵.一級輕氣炮內(nèi)彈道數(shù)值模擬[J].彈道學(xué)報,1998,10(4):21-25.

YANG Jun-yun,JIN Zhi-ming,ZHANG Xiao-bing.Numerical simulation of the interior ballistic process in one-stage light-gas gun[J].Journal of Ballistics,1998,10(4):21-25.(in Chinese)

[8]仲偉君,趙曉利,齊杏林.氣體炮內(nèi)彈道建模與發(fā)射環(huán)境模擬研究[J].動力學(xué)與控制學(xué)報,2005,3(1):61-65.

ZHONG Wei-jun,ZHAO Xiao-li,QI Xing-lin.Study on the modeling of the internal ballistic of a gas gun and launching environment simulation[J].Journal of Dynamics and Control,2005,3(1):61-65.(in Chinese)

[9]趙希欣,高元樓.一種氣體炮的建模與分析[J].液壓氣動與密封,2012(10):43-44.

ZHAO Xi-xin,GAO Yuan-lou.Modeling and analysis of a gas gun[J].Hydraulics Pneumatics and Seals,2012(10):43-44.(in Chinese)

[10]胡靜,張明安,岳文龍,等.內(nèi)彈道過程中氣體工質(zhì)做功能力計算分析[J].火炮發(fā)射與控制學(xué)報,2011(2):82-85.HU Jing,ZHANG Ming-an,YUE Wen-long,et al.Calculation and analysis of gas power capability in process of interior ballistics[J].Journal of Gun Launch and Control,2011(2):82-85.(in Chinese)

[11]夏正友,李建華,褚永勝.基于神經(jīng)網(wǎng)絡(luò)的氣體炮內(nèi)彈道仿真與預(yù)測[J].現(xiàn)代引信,1998(4):36-40.XIA Zheng-you,LI Jian-hua,CHU Yong-sheng.Imitation and prediction of the gas gun interior ballistics based on neural network[J].Modern Fuze,1998(4):36-40.(in Chinese)

[12]錢林方,侯保林,徐亞棟.火炮彈道學(xué)[M].北京:北京理工大學(xué)出版社,2009.

QIAN Lin-fang,HOU Bao-lin,XU Ya-dong.Artillery ballistics[M].Beijing:Beijing Institute of Technology Press,2009.(in Chinese)

Modification and Validation of Interior Ballistic Equation for One-stage Gas Gun

LI Feng,BAI Yun-shan,FENG Xiao-wei

(Institute of System Engineering,China Academy of Engineering Physics,Mianyang 621999,China)

In order to improve the accuracy and reliability of the theoretical calculation result,the classical equation for interior ballistics of one-stage gas gun was modified.The problem of classical equation was analyzed.The result shows that the influence of air chamber volume on bullet velocity can not be ignored.A new expression of gas density was presented.According to the equation of mass conservation,the momentum equation of flow,the motion equation of projectile,the adiabatic expansion equation of gas and ideal-gas equation,the modified equation of one-stage gas gun were deduced.The calculation results of muzzle velocity by classical equation and amendatory equation were compared to the result of numerical simulation of interior ballistics.The results show that the muzzle velocity obtained by above mentioned three methods is almost same when the volume ratio of air chamber to gun barrel is very small,but the differences between the results obtained by classical equation and the results obtained by the modified equation and numerical simulation increase markedly with the increase of volume ratio.The classical equation is only suitable for calculating velocity of gas gun with small volume ratio,and the modified equation can be used without volume ratio limit.

gas gun;interior ballistics;muzzle velocity;volume ratio

2014-09-02

中國工程物理研究院高新裝備專項

李鋒(1980- ),男,工程師,碩士研究生,研究方向為氣體炮設(shè)計及氣體驅(qū)動技術(shù)。E-mail:15983633005@163.com。

TJ012.1

A

1004-499X(2016)01-0014-05