/KJ*kg-1    

比热比 1.20 薄火药药孔数 1 

挤进压力/MPa 30.0 厚火药质量/kg 1.625 

次要功系数φ 1 1.02 厚火药火药力/KJ*kg-1 980 

火药余容/dm3*kg-1 1.0 厚火药弧厚/mm 1.03 

装药总数 2 厚火药内孔直径/mm 0.50 

薄火药质量/kg 0.34 厚火药长度/mm 12.0 

热损失系数 Kq 0.18 厚火药药孔数 7 

火药力/KJ*kg-1 980.0   

基于 MATLAB 平台，编写程序，完成对内弹道的计算[7-12]。 对于本炮，我们在底凹弹常温全装药、底凹弹低温全装药以及底凹弹高温全装药条件下

各求一次膛压曲线，计算 p-t，结果如图 1.1 所示。p-t 为表格函数,曲线如图 1.1 所示：

由上图可得知，采用常温底凹弹时所受到的平均压力最大。所以在下文的分析计算中， 我们主要针对常温条件下、底凹弹装弹进行后续计算。

本科毕业设计说明书 第 9 页

常温下、底凹弹具体数据详见附录 A。

Pmax=230.7MPa v0=530m/s

2．2 求出炮膛合力曲线

由已有理论得知，炮膛合力分为三个阶段，膛内运动时期、火药气体后效期和惯性后坐 时期[13]。

（1）对于膛内时期有：

Φ 为次要功计算系数，可由下式计算：

k 1 

3 m

K 为经验系数，因为本炮为轻型榴弹炮，所以取 k=1.06

由于 p-t 为表格函数，因此需对 Fpt 进行插值,插值方法选用的是拉格朗日插值法，对函 数进行一元三点不等距插值，计算公式为：

对于拉格朗日插值函数，我们在 MATLAB 中编写如下程序，以便于之后求炮膛合力时方 便调用。详见附录 C。

上式中：b 为后效期时间常数； 为炮口制退器冲量特征量。

几个特征量的计算如下：

当结束后效期时 PK=0.1764MPa。对中小口径榴弹炮，取 B=1300

（3）惯性时期，炮膛合力为 0：来!自~751论-文|网www.751com.cn

带入已知参数： v0=530m/s, A=1.1907e-2m2,

对应炮口制退器效率为 42%, tg=6.3421ms, pg=68.17633MPa,

Fpt=0 （2.43）

可用 matlab 计算出出 Fpt 曲线，编程内容见附录 B，所得结果见图 2.1。