沈曉軍，王雨時，聞泉，徐鷺林，劉延中

(1. 南京理工大學(xué) 機械工程學(xué)院，江蘇 南京 210094； 2. 吉林市江機民科實業(yè)有限公司，吉林 吉林 132021)

用Vallier公式近似估算槍炮引信膛內(nèi)后坐過載規(guī)律

沈曉軍1，王雨時1，聞泉1，徐鷺林2，劉延中2

(1. 南京理工大學(xué) 機械工程學(xué)院，江蘇 南京 210094； 2. 吉林市江機民科實業(yè)有限公司，吉林 吉林 132021)

為了得出火藥氣體壓力與時間關(guān)系的解析式，在前人得出的內(nèi)彈道半實驗解法的基礎(chǔ)上，對火藥氣體壓力與時間的關(guān)系近似地計算，以前人計算得出的火藥氣體壓力與時間關(guān)系為依據(jù)，論證用Vallier公式簡化計算火藥氣體壓力的可行性，并用中間值迭代方法對給定一些條件后的火藥氣體壓力曲線進行逼近，得出逼近后的Vallier公式曲線，將各平臺火炮內(nèi)彈道修正系數(shù)的平均值作為經(jīng)驗值來計算火藥氣體壓力。結(jié)果表明，Vallier公式能夠較好地計算出一般火炮的火藥氣體壓力和無后坐火炮的火藥氣體壓力上升段。

火炮；內(nèi)彈道；后坐過載；估算；火藥氣體壓力；Vallier公式

0 引言

在求解引信后坐機構(gòu)運動微分方程時，需要知道膛內(nèi)火藥氣體與時間的關(guān)系式，但對精度的要求不是很高?；鹚帤怏w壓力曲線主要通過實驗和傳統(tǒng)內(nèi)彈道學(xué)計算兩種方法得到，如今已有許多精確解法，并且編制了各種內(nèi)彈道表。然而隨著裝藥結(jié)構(gòu)的差異、火藥燃燒規(guī)律的多樣性，這些解法還需要根據(jù)實驗修正[1]，工作量很大。法國瓦里耶(E.Vallier)提出以時間為自變量表示的火藥氣體壓力的實驗式[2]。文獻[3]介紹的海登賴希解法，是瓦里耶解法的發(fā)展，但計算還是比較復(fù)雜，考慮進了火藥參數(shù)，達不到簡化的目的。文獻[4]在前人工作的基礎(chǔ)上，給出了一種用函數(shù)式表示內(nèi)彈道時期膛壓變化規(guī)律的方法，但此函數(shù)中待定系數(shù)的求解有一定困難，不夠簡化。

文中應(yīng)用Vallier公式對不同火炮發(fā)射時的火藥氣體壓力進行計算，并與實際火藥氣體壓力進行對比，論證Vallier公式的可行性和適用范圍。提出用迭代的方式計算實際情況下給出某些已知條件時的火藥氣體壓力，對比火藥參數(shù)計算法和迭代法計算得出的火藥平均壓力代入Vallier公式得出的曲線，證明迭代法的準(zhǔn)確性。計算出不同平臺不同口徑的內(nèi)彈道修正系數(shù)最優(yōu)值的平均值作為經(jīng)驗值。

1 Vallier公式

法國瓦里耶(E.Vallier)提出的Vallier公式[2]：

式中:Pm——膛內(nèi)最大壓力；

tm——最大膛壓時刻的時間值；

β1——內(nèi)彈道修正系數(shù)。

按照Vallier公式，t=0時，Pb也為0。而實際情況下，膛內(nèi)存在擠進壓力P0，所以需要對瓦里耶公式進行修正。修正后的公式為[2]：

修正前和修正后的Vallier公式計算的XX年式76.2 mm加農(nóng)炮火藥氣體壓力曲線如圖1所示。

圖1 修正前后的Vallier公式

2 Vallier公式近似估算火藥氣體壓力的可行性分析

文獻[4]在MZ-80B型微機上使用BASIC語言在經(jīng)典彈道學(xué)基礎(chǔ)上編制了數(shù)值模型，彈道方程的求解采用四階龍格-庫塔法，特殊點數(shù)值求解時，最大壓力點采用二次插值法或黃金分割法；燃燒結(jié)束點、分裂點、炮口點采用變量更換法。其中的某些火炮內(nèi)彈道參數(shù)如表 1所示。 各種火炮的實際火藥氣體壓力與Vallier公式計算值對比情況如圖2所示。

表1 不同火炮的火藥氣體壓力與時間關(guān)系[4]

圖2 各種火炮的火藥氣體壓力與Vallier計算值對比圖

由圖2可知，對于一般火炮修正后的Vallier公式計算值與實際火藥氣體壓力較接近，而對于迫擊炮彈，兩者相差較大，對于無后坐炮，在壓力上升段較接近，在壓力下降段相差較大。總之，當(dāng)對火藥氣體壓力要求不是很精確時，一般火炮的火藥氣體壓力曲線和無后坐炮火藥氣體壓力的上升段都可以用Vallier公式計算。

3 Vallier公式的逼近

圖3 計算流程圖

4 算例

圖4 Vallier公式第一次逼近

圖5 Vallier公式第二次逼近

圖6 Vallier公式第二次逼近

圖7 最后曲線與實際火藥氣體壓力的對比

由圖7可見，根據(jù)Vallier公式計算出的火藥氣體的壓力與實際情況較為接近，說明Vallier公式計算火藥氣體壓力的準(zhǔn)確性。

火藥氣體壓力平均值還可以通過火藥參數(shù)直接計算得出，計算公式為：

式中：q——彈丸質(zhì)量；

ω——裝藥質(zhì)量；

v——初速；

g——重力加速度；

S——身管內(nèi)膛橫斷面面積；

l——彈丸全行程長。

用火藥氣體參數(shù)計算出的平均壓力代入Vallier公式進行計算，得出的曲線如圖8所示。

圖8 各火炮火藥氣體壓力

在上升段迭代法和火藥參數(shù)計算法得出的曲線相差很小，迭代法甚至更接近于實際火藥氣體壓力曲線；在下降段火藥參數(shù)計算法的炮口壓力更接近于實際壓力曲線。在引信設(shè)計過程中，主要考慮火藥氣體上升段對彈丸的作用，因此在引信設(shè)計過程中，可以使用迭代法用Vallier公式計算火藥氣體壓力。

5 內(nèi)彈道修正系數(shù)β1的估算

文獻[3]中提出內(nèi)彈道修正系數(shù)與火藥性能有關(guān)，例如黑火藥的β1=1，速燃藥β1=2～3，緩燃藥的β1=4。

運用Vallier公式對不同平臺不同口徑的火炮進行迭代計算，得出最優(yōu)內(nèi)彈道修正系數(shù)，結(jié)果如表2-表7所示。對不同火炮取平均值作為以后估算膛內(nèi)火藥壓力的計算值。

表2 加農(nóng)炮的內(nèi)彈道修正系數(shù)

表3 高射炮的內(nèi)彈道修正系數(shù)

表4 坦克炮的內(nèi)彈道修正系數(shù)

表5 艦炮的內(nèi)彈道修正系數(shù)

表6 航空炮的內(nèi)彈道修正系數(shù)

表7 無后坐炮的內(nèi)彈道修正系數(shù)

用內(nèi)彈道修正系數(shù)的平均值和最優(yōu)值分別代入Vallier公式進行計算，得出的XX年式76.2 mm口徑加農(nóng)炮火藥氣體壓力曲線與時間的關(guān)系如圖8(a)所示，由圖8(a)可知，平均值和最優(yōu)值得出的火藥氣體壓力曲線相差不多，可以用平均值作為估算火藥氣體壓力，如圖9所示。

圖9 選用不同內(nèi)彈道修正系數(shù)的Vallier公式 表示的火藥氣體壓力曲線

6 結(jié)語

1) 對于一般火炮的火藥氣體壓力和無后坐炮火藥氣體壓力的上升段，可以用Vallier公式進行估算，對于迫擊炮，Vallier公式的估算結(jié)果相差較大。

3) 在引信設(shè)計過程中，可以使用迭代法用Vallier公式計算火藥氣體壓力。在上升段迭代法和火藥參數(shù)計算法得出的曲線相差很小，迭代法甚至更接近于實際火藥氣體壓力曲線；在下降段火藥參數(shù)計算法的炮口壓力更接近于實際壓力曲線。

4) 對于不同的平臺火炮可以取不同內(nèi)彈道修正系數(shù)，選用這些內(nèi)彈道修正系數(shù)的Vallier公式能夠較快地表達出火藥氣體壓力。加農(nóng)炮取3.39，高射炮取3.17，坦克炮取3.58，艦炮取3.36，航空炮取2.18，無后坐炮取3.09。

[1] 苪筱亭，楊啟仁. 彈丸發(fā)射過程理論[M]. 南京: 東南大學(xué)出版社，1992.

[2] 苪筱亭. 彈丸在膛內(nèi)運動的研究(一)[J]. 兵工學(xué)報彈箭分冊，1987：4-6.

[3] 蔡偉妹，易連軍，徐萬和. 一種新的內(nèi)彈道膛壓計算方法[J]. 四川兵工學(xué)報，2012, 33(2)：16-17.

[4] 李晝堂. 火藥與內(nèi)彈道[M]. 北京: 兵器工業(yè)出版社，2001.

[5] 王連榮，張佩勤. 火炮內(nèi)彈道計算手冊[M]. 北京: 國防工業(yè)出版社，1987.

Approximate Estimate of Law of Breech Pressure of Gun Fuze with Vallier Formula

SHEN Xiaojun1, WANG Yushi1, WEN Quan1, XU Lulin2, LIU Yanzhong2

(1. School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China;2. Jilin Jiangjiminke Industry Co., Ltd., Jilin 132021, China )

In order to draw the analytic relationship between gunpowder gas pressure and time, based on previous results of the semi-empirical formula of ballistic solution, the relationship between gunpowder gas pressure and time is approximately calculated. The feasibility of simplifying the gunpowder gas pressure calculation by Vallier formula is demonstrated. The powder gas pressure curve is approximated after using the intermediate value iterative method to give out some conditions. The average value of trajectory correction coefficient with each platform artillery is calculated as the experience value of the gunpowder gas pressure. The results show that Vallier formula can be used to calculate the gas pressure of general artillery gunpowder and ascent of recoilless artillery gunpowder gas pressure well.

artillery; ballistic; recoil overload; estimation; gunpowder gas pressure; Vallier formula

沈曉軍(1989-)，男，江蘇如皋人，碩士研究生，主要研究方向為引信系統(tǒng)設(shè)計。

TJ303

B

1671-5276(2015)05-0192-05

2014-02-26