/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。