0 引言
采用气体二次喷射推力矢量技术,能够提高控制效率、降低消极质量、提升工作稳定性,是未来极具发展潜力的推力矢量技术,因此吸引国内外学者开展大量理论与验证工作[1-8]。俄罗斯和乌克兰在燃气二次喷射推力矢量控制技术的研究较为成熟,已在地对地固体战略导弹“白杨/镰刀”、三级潜射洲际弹道导弹“P-39”等发动机上获得实际应用,美国NASA兰利研究中心开展了气体二次喷射推力矢量方法和试验研究。近年来国内开展了一些气体二次喷射试验和仿真工作,研究成果和结论具有重要意义,但未来仍需结合实际应用,对二次喷射的矢量控制机理和规律进行深入研究。文中采用优化方法对固体火箭发动机气体二次喷射复杂干扰喷管内流场进行了数值模拟,研究了二次喷射流与主流相互耦合作用机制以及二次喷射推力矢量流动机理,分析了不同喷射参数对气体二次喷射流场特征及侧向控制力的影响规律。
1 物理模型和计算方法
2 计算结果及交互作用分析
表1 因素水平 |
| 水平 | A喷射 位置/mm | B喷射 角度/(°) | C二次流 总压/MPa | D二次流 总温/K |
|---|---|---|---|---|
| 1 | 30.6 | 45 | 3 | 300 |
| 2 | 61.2 | 0 | 5 | 1 800 |
| 3 | 91.8 | -45 | 6 | 3 000 |
表2 表头设计 |
| 因素 | A | B | (A×B)1 | (A×B)2 | C | D | - | (B×C)1 | - | - | (B×C)2 | - | - |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 列号 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 |
矢量角定义为δ=arctan(Fn/Fa),Fn为发动机所受的侧向力,Fa为发动机所受轴向力。轴向力和侧向力均由对矢量喷管出口截面的压强和动量积分得到。
通过主效应影响分析(图2):二次射流喷射的位置、二次射流总压对矢量角影响较大,而二次射流的喷射角度、二次流总温对矢量角影响较小。采用通用线性回归进行方差分析发现:A、B、C、D 4个因素变化对矢量角总的贡献率为 66.10%;其中A的贡献率为24.2%,B为0.57%,C为36.94%,D为4.4%。如果不考虑因素间的交互作用,则选取好的因素水平组合为A3B1C3D1。AB的交互作用见表3,BC的交互作用见表4。从AB二元表分析看,A3B2或A3B3搭配组合获得的矢量角平均值较大;而从BC二元表分析看,B2C3搭配组合获得的矢量角平均值较大。所以因素B应取B2,此时满足A×B,B×C交互作用的最好搭配。总结出获得最大矢量角的组合方式:喷射位置为91.8 mm,喷射角度0°,二次流总压6 MPa,二次流总温300 K。
3 喷管内流场分析与试验验证
参考总结出的获得最大矢量角组合方式,确定气体二次喷射推力矢量控制方案:从附近的燃气发生器引入二次流,喷射位置在主喷管扩散中间,喷射角度为0°,二次流总压6 MPa,二次流总温1 500 K。
3.1 主次流干扰喷管内流场分析
图3为二次喷射矢量喷管内横向和纵向剖面马赫数分布图,分析流场结构发现射流出口压强高于附近主流,射流在主流区迅速膨胀,但是射流出口会受到超声速流压缩。当二次射流总压增大到一定程度,弓形激波就会越过中心线,逼近下壁面。继续增大二次流总压,弓形激波就与下壁面相交,产生反射激波,干扰到下壁面边界层,造成对应位置的边界层分离。改变喷射条件,不会改变主次流干扰流场结构,会改变分离区域的大小。二次喷射口区域存在高压区和低压区,这是由于存在流动边界层与激波的相互干扰形成高压区,气体剧烈膨胀和喷流引射作用导致喷口下游存在低压区,两者综合作用形成主次流干扰机理。
3.2 试验测试与分析
表5 主、副发动机尺寸 |
| 尺寸 | 主发动机 | 副发动机 |
|---|---|---|
| 燃烧室内径D/mm | 135 | 96 |
| 喷管喉径d/mm | 11.2 | 4.3 |
| 喷管扩张半角α/(°) | 12 | 12 |
| 喷管扩张比(de/dt) | 2.7 | 1.5 |
表6 低温双基药基本性能 |
| 性能参数 | 参数值 | 说明 |
|---|---|---|
| 稳定工作压强/MPa | 3~20 | |
| 燃速压强指数 | ≤0.2 | 3~10 MPa |
| 工作温度/℃ | -45~+60 | |
| 燃速温度系数/(%·K-1) | <0.2 | |
| 燃烧温度/℃ | ≤2 000 | |
| 静态燃速/(mm/s) | 4.5 | 常温7 MPa |
4 结论
对气体二次喷射推力矢量发动机的喷管内流场进行了模拟,得到二次射流与主流干扰的激波系结构,优化喷射参数总结出最大矢量角组合方式,通过试验验证了优化方案。
1)通过DOE设计与交互作用分析,发现影响气体二次喷射推力矢量角δ的主次顺序为:二次喷射的位置、二次流总压、二次流总温、喷射位置与喷射角度的交互作用,喷射角度与二次流总压的交互作用,喷射角度。
2)在相同的二次射流总压条件下,带一定角度的逆流喷射比垂直喷射、顺流喷射的侧向力大;二次流总温减少时,使得由小股喷流动量引起的侧向力升高,推力矢量偏角略微增加。
3)研究气体二次喷射复杂喷管内流场,分析详细流场结构,得到边界层与激波的相互干扰机理。根据总结出的最大矢量角组合方式及参数,开展试验研究,优化方案的矢量角为5.7°。