短时间工作发动机内弹道计算方法

李春艳; 王立武; 刘沛; 林志远

doi:10.15892/j.cnki.djzdxb.2019.05.028

弹箭与制导学报 >

2019 , Vol. 39 >Issue 5: 119 - 122

DOI: https://doi.org/10.15892/j.cnki.djzdxb.2019.05.028

短时间工作发动机内弹道计算方法

李春艳 ¹^,² ,
王立武 ¹^,² ,
刘沛 ¹^,² ,
林志远 ¹^,²

展开

¹ 中国航天科技集团有限公司第四研究院第41研究所, 西安 710025
² 固体火箭发动机燃烧、热结构与内流场国防科技重点实验室, 西安 710025

李春艳(1990-),女,湖南常宁人,工程师,研究方向:固体火箭发动机。

收稿日期: 2018-10-08

网络出版日期: 2025-05-12

收起

Investigation on the Calculating Method for Short Time Working Solid Rocket Motor

LI Chunyan ¹^,² ,
WANG Liwu ¹^,² ,
LIU Pei ¹^,² ,
LIN Zhiyuan ¹^,²

Expand

¹ The 41st Institute of the Fourth Academy,CASC, Xi’an 710025, China
² National Key Laboratory of Combustion Flow and Thermo-structure, Xi’an 710025, China

Received date: 2018-10-08

Online published: 2025-05-12

Fold

摘要

为了提高短时间工作固体火箭发动机的内弹道计算精度,对现有固体发动机内弹道性能计算方法进行适应性修改,建立了适用于短时间工作固体发动机的高精度内弹道计算方法。通过预示内弹道性能与实测发动机试验数据的对比,修改后计算方法的预示精度在±5%以内,计算精度显著提高。计算结果表明,在发动机平衡段,喷管效率和理论推力系数基本保持不变,可以为短时间工作发动机的设计提供依据。

关键词： 短时间工作; 发动机; 内弹道; 计算方法

本文引用格式

李春艳 , 王立武 , 刘沛 , 林志远 . 短时间工作发动机内弹道计算方法[J]. 弹箭与制导学报, 2019 , 39(5) : 119 -122 . DOI: 10.15892/j.cnki.djzdxb.2019.05.028

Abstract

In order to increase the calculation accuracy of performance for short time working solid rocket motor, a method is established by modifying the conventional one. The analysis of calculation and test results show that the calculation accuracy is not more than ±5 percent, which means that the accuracy of calculation is improved significantly. The calculation results show that nozzle efficiency and theoretical thrust coefficient will remain unchanged when working in high pressure. This calculation method can provide an important reference for the design of short time working solid rocket motor.

Key words： short time working; solid rocket motor; performance; calculation method

0 引言

在导弹的级间分离环节,为提高分离可靠性,广泛采用短时间、大推力固体火箭发动机为分离提供动力^[1]。此类固体发动机工作时间很短,一般小于1 s,在内弹道预示过程中存在几个难点:1)推进剂燃速较高,燃速偏差较大;2)标准发动机与全尺寸发动机燃速相关性不确定;3)在进行内弹道反算过程中,处理发动机的试车数据时,难以准确地确定推进剂的燃烧时间。短时间工作小型固体发动机的内弹道预示精度较低,不得不增加试验数量以摸清发动机的内弹道性能,增加了研制成本。

目前,关于固体发动机内弹道性能预示及分析的书籍^[2-5]并不适用于短时间工作固体发动机的内弹道性能分析。文中建立了一种适用于短时间工作固体发动机的高精度内弹道计算方法,介绍了内弹道计算方法的理论模型,并通过预示内弹道性能与实测发动机试车数据的对比分析,阐明了该方法的适用性。

1 理论模型

对短时间工作发动机影响最大的是燃烧时间,其次为燃速相关性和推进剂燃速偏差。短时间工作固体发动机内弹道计算的流程见图1,利用内弹道反算获得推进剂燃烧规律、喷管效率及推力系数随肉厚的变化关系,以此为基础进行内弹道性能计算,预示发动机的工作压强和推力,并进行计算精度分析。

显示原图|下载原图ZIP|生成PPT

图1 短时间工作小型固体发动机内弹道计算方法流程图

1.1 发动机燃烧时间的确定

短时间工作发动机药柱肉厚较小、燃速较高,能否准确判断发动机燃烧时间直接影响内弹道计算的精度。根据大量实验分析,发动机燃烧时间压强冲量与工作时间压强冲量存在如下关系^[1]:

(1)

K = ∫ 0 t b P c · d t ∫ 0 t a P c · d t

式中:K为与发动机相关的试验修正系数;t_b为发动机燃烧时间;t_a为发动机工作时间;P_c为发动机试验工作压强。

发动机燃烧时间按照式(1)迭代求解,可以较为准确地确定发动机的燃烧时间,从而提高内弹道预示的精度。

1.2 推进剂燃面计算

短时间工作发动机喷管喉部直径一般较小,发动机工作结束后喉部烧蚀极小,一般难以测量。因此,在进行内弹道燃面反算时认为喷管喉部直径不变,将喷管喉部直径的变化因素折合到推进剂的燃烧规律,按式(2)计算发动机工作过程中的推进剂燃烧规律,并采用最小二乘法获得多台发动机平均燃烧规律^[6]。

(2)

A b = P c 1 - n · A t ρ · a · c

式中:A_b为推进剂反算燃烧规律;P_c为燃烧室压强;A_t为喷管喉径;n为压强指数;ρ为推进剂密度;a为燃速系数;c为特征速度。

1.3 喷管效率和实际推力系数计算

根据小型固体发动机试车数据,利用式(3)计算喷管效率,利用式(4)计算实际推力系数。

(3)

η = F C f · P c · A t

(4)

C f r e a l = F P c · A t

式中:η为喷管效率;F为发动机实测推力;C_f为理论推力系数;C_freal为实际推力系数。

1.4 发动机工作压强预示

对于短时间工作发动机,由于不考虑喷管喉部烧蚀,可以根据上述计算以及推进剂性能参数,按照内弹道计算公式(5)计算发动机的工作压强^[2-4]。

(5)

P c 预 示 = ρ · a · c A b A t 1 1 - n

1.5 发动机推力预示

发动机试车过程中,环境压强P_a为当地大气压强,存在一个使燃气完全膨胀的工作压强P_L。短时间工作发动机工作压强较高,一般大于P_L,但是在发动机工作上升段和下降段,存在工作压强小于P_L的时间段,此时理论推力系数计算误差已失真。因此,为了提高内弹道预示精度,根据喷管扩张比计算公式(6),采用二分法迭代求解燃气处于完全膨胀状态下的工作压强P_L,当预示压强P_c预示≥P_L时,采用式(7)计算发动机的预示推力;当预示压强P_c预示<P_L时,采用式(8)计算发动机的预示推力^[2-5]。

(6)

A e A t = 2 k + 1 1 k - 1 k - 1 k + 1 P e P c 2 k - P e P c k + 1 k

(7)F=η·C_f·P_c_预示·A_t

(8)F=C_freal·P_c_预示·A_t

1.6 工作时间预示

发动机工作压强和工作时间是一一对应的,可以通过燃面-肉厚将二者对应起来,从而可以得到压强-燃面-肉厚-时间的对应关系。计算过程中采用基于肉厚步长的算法,每个肉厚对应的发动机工作时间可以按式(9)和式(10)计算。

(9)t_i+₁=t_i+Δt

(10)

Δ t = Δ w a · P c, i + P c, i + 1 2 n

式中:t_i为第i时刻发动机的工作时间;Δt为发动机第i+1时刻与第i时刻的工作时间差;Δw为发动机第i+1时刻与第i时刻之间的肉厚差;P_c,_i为第i时刻发动机预示工作压强。

2 计算结果及分析

利用上述理论模型方法,采用Matlab软件编写了短时间工作小型固体发动机内弹道计算程序^[7-8]。以某发动机试车结果为例,分析短时间工作小型固体发动机内弹道计算的精度。为了考核该内弹道计算方法的有效性,选择不同温度条件进行发动机试验,发动机试验过程的温度条件见表1。

表1 发动机试验温度条件

	试验编号
	#1	#2	#3	#4	#5	#6	#7
试验温度/℃	20	20	20	20	55	55	-20

图2给出了通过内弹道反算得到的喷管效率随已燃烧肉厚的变化曲线,可以看出:在发动机平衡段,喷管效率基本保持不变,但是在发动机工作初期及末段,喷管效率围绕某个值震荡。分析认为,在发动机平衡段,压强和推力较高、波动较小,故发动机在此段肉厚下计算喷管效率基本不变;在发动机工作初期及末段,压强和推力较低、波动较大,同时压强和推力测试存在不同步,导致发动机在此段肉厚下计算喷管效率存在震荡现象。

显示原图|下载原图ZIP|生成PPT

图2 喷管效率随已燃烧肉厚的变化曲线

图3给出了内弹道反算得到的理论推力系数随已燃烧肉厚的变化曲线,可以看出:在发动机平衡段,理论推力系数基本为恒定值,当压强低于P_L时,理论推力系数开始不断下降。分析认为,在发动机平衡段,工作压强较高,燃气处于完全膨胀或欠膨胀状态,理论推力系数基本为恒定值;当压强低于P_L时,燃气处于过膨胀状态,理论推力系数不断下降。

显示原图|下载原图ZIP|生成PPT

图3 理论推力系数随已燃烧肉厚的变化曲线

图4和图5给出了20 ℃下压强、推力预示结果与实测结果对比曲线,可以看出,预示结果与实测结果吻合。

显示原图|下载原图ZIP|生成PPT

图4 20 ℃下预示压强与实测压强的对比曲线

显示原图|下载原图ZIP|生成PPT

图5 20 ℃下预示推力与实测推力的对比曲线

表2给出了不同温度下试验发动机的内弹道性能预示精度,可以看出:该内弹道计算方法的预示精度在±5%以内,可以为短时间工作发动机内弹道预示提供依据。

表2 内弹道性能预示精度

序号	温度/℃	压强预示精度/%	推力预示精度/%
1	20	0.23~0.4	0.1~2.0
2	55	-1.6~0.7	-1.9~3
3	-20	-2.38	-2.16

3 结论

通过研究,得出以下结论:

1)建立了一种适用于短时间工作小型固体发动机的高精度内弹道计算方法,针对文中的算例,内弹道性能的预示精度在±5%以内,可以为短时间工作固体发动机的内弹道性能预示提供依据。

2)在发动机平衡段,喷管效率基本保持不变;在发动机工作初期及末段,喷管效率围绕某个值震荡,这是由于较低压强下压强和推力测试不同步引起的。

3)在发动机平衡段,理论推力系数基本为恒定值;当工作时间增加到一定值后,即压强低于燃气完全膨胀时工作压强P_L后,理论推力系数开始不断下降。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

[1]	郭凤美, 余梦伦. 导弹分离设计技术研究[J]. 导弹与航天运载技术, 2014(1):5-10.

[2]	陈汝训. 固体火箭发动机设计与研究:上册[M]. 北京: 中国宇航出版社, 2009:244-251.

[3]	杨月诚. 火箭发动机理论基础[M]. 西安: 西北工业大学出版社, 2010:115-123.

[4]	唐金兰, 刘佩进. 固体火箭发动机原理[M]. 北京: 国防工业出版社, 2013:266-285.

[5]	鲍福廷, 侯晓. 固体火箭发动机设计[M]. 北京: 中国宇航出版社, 2016:29-38.

[6]	韩旭里. 数值分析[M]. 北京: 高等教育出版社, 2011:65-70.

[7]	刘慧颖. MATLAB R2007基础教程[M]. 北京: 清华大学出版社, 2008:21-40.

[8]	MOORE Holly. MATLAB实用教程[M].高会生,刘童娜,李聪聪,译. 北京: 电子工业出版社, 2016:89-122.

Options

文章导航

模态框（Modal）标题

摘要

本文引用格式

Abstract

0 引言

1 理论模型

图1 短时间工作小型固体发动机内弹道计算方法流程图

1.1 发动机燃烧时间的确定

1.2 推进剂燃面计算

1.3 喷管效率和实际推力系数计算

1.4 发动机工作压强预示

1.5 发动机推力预示

1.6 工作时间预示

2 计算结果及分析

表1 发动机试验温度条件

图2 喷管效率随已燃烧肉厚的变化曲线

图3 理论推力系数随已燃烧肉厚的变化曲线

图4 20 ℃下预示压强与实测压强的对比曲线

图5 20 ℃下预示推力与实测推力的对比曲线

表2 内弹道性能预示精度

3 结论

参考文献