高超声速飞行器是我国正在研制的重要武器装备, 由冲压发动机驱动、高超声速巡航, 具有飞行速度快、机动飞行、难以拦截等优势。马赫数4.5是高超声速飞行器、战略战术导弹等飞行器在上升段和再入段的关键马赫数, 也是亚燃冲压向超燃冲压过渡马赫数, 马赫数4.5的试验数据对飞行器研制至关重要。
目前国内工业性质的0.5 m量级高超声速风洞均无马赫数4.5的风洞试验能力。在高超声速风洞研制马赫数4.5喷管, 可以填补国内空白, 具有重要意义。通过马赫数4.5喷管的气动设计、计算验证、结构设计、加工、流场校测、标模试验, 验证了研制的马赫数4.5喷管流场满足试验要求, 可以为飞行器研制提供可靠的数据支撑。
1 气动设计气动设计采用较为成熟、可靠的特征线法。主要分为以下3步:
1) 根据等熵理想流体运动, 计算喷管的理想型面;
2) 按照气体的黏性影响只在喷管壁面附近, 计算喷管壁面的附面层位移厚度;
3) 修正喷管理想型面, 给出喷管的物理型线。
马赫数4.5喷管采用轴对称型面喷管, 喷管型面气动设计采取文献[1-3]方法计算: 在亚声速段收缩段采用维托辛斯基方法计算; 在喉道段采用3次曲线方法计算; 在超声速段采用特征线法计算, 并采用Sivells-Payne方法进行边界层位移厚度修正。设计计算图如图 1。
![]() |
图 1 轴对称喷管气动设计计算图 |
依照以上方法, 表 1给出了M4.5喷管型面设计计算的主要结果, T, A, B, C, D, E分别为图 1中的位置点, X为横坐标, Y为纵坐标, M为马赫数, T0为温度, 图 2给出了马赫数4.5的型面曲线。
点号 | X/mm | Y/mm | M | T0/(°) |
T | 0.000 | 61.43 | 1.000 | 0.000 |
A | 926.56 | 142.75 | 3.256 | 7.500 |
B | 1 596.60 | 0.00 | 4.300 | 0.000 |
C | 1 989.51 | 0.00 | 4.500 | 0.000 |
E | 3 086.38 | 250.00 | 4.500 | 0.000 |
D | 1 142.56 | 170.34 | 3.604 | 6.823 |
![]() |
图 2 M4.5喷管型面曲线 |
本文对M4.5喷管流场分布特性计算采用商用Gambit软件来生成网格。Gambit的主要优点体现在其突出的非结构化网格生成能力和生成过程较强的自动化能力, 包含了数十种网格生成方法, 能够针对复杂的几何外形生成三维四面体、六面体的非结构化网格及混合网格, 同时Gambit比较注重软件的图形用户界面(GUI), 交互式特点直观易学。
M4.5喷管的计算网格见图 3, 网格内部间距10 mm。当网格生成完毕, 要定义网格使用范围, 并定义好边界条件, 以备FLUNET软件计算使用。
![]() |
图 3 M4.5喷管流场计算网格 |
FLUENT软件提供了从不可压到可压、层流、湍流等范围很广的模拟计算能力。软件中输运现象的数学模型与所模拟的几何图形是结合在一起的。湍流模型是FLUENT中很重要的一部分, 湍流会影响到其他的物理现象如浮力和可压缩性。湍流模型具有很大的适用范围, 而不需要对特定的应用做特殊调节, 而且它涵括了其他物理现象的影响, 如浮力和可压缩性。通过使用扩展壁面函数和区域模型, 对近壁面精度问题有很好的处理能力。
本文利用FLUNET软件来模拟计算前面设计的M4.5喷管的流场分布。所计算喷管流场的计算区域由壁面、压力入口条件、压力出口条件和对称轴所包围的区域组成。图 4为层流模型计算残差收敛过程图, 从图中可以看出, 计算收敛过程较快, 说明网格质量和计算所选用模型较好。
![]() |
图 4 计算残差收敛过程 |
对于高超声速风洞来说速度场一般指喷管内和试验段的马赫数分布规律, 本文主要针对速度场, 温度场, 方向场进行简要分析, 目的在于检验本文计算的喷管型面坐标是否符合设计要求。
M4.5喷管的马赫数分布见图 5, 从图中可以看出, 在喷管内和喷管出口处的马赫数分布趋势同设计分布一致, 能得到试验所需的流场参数。
![]() |
图 5 马赫数云图 |
在喷管设计方面, 一般最关心的是喷管内部和出口处的流场特性情况, 为了更容易比较, 表 2利用FLUENT计算结果给出了M4.5喷管出口速度场指标分析结果。从表 2可以看出: 在最大马赫数偏差小于1.0%的情况下, k-ε湍流模型计算出的均匀区直径D为0.37 m, 平均马赫数Mj-AVE为4.479, 均方差σMj为0.005 3;而层流模型计算出的均匀区直径为0.48 m, 平均马赫数为4.499, 均方差为0.009;这满足参考文献[5]中规定的指标要求, 同时也符合参考文献[6]中规定的指标要求。其中Mj-MAX为最大马赫数, Mj-MIN为最小马赫数, ΔMj-MAX为最大马赫数偏差。
计算模型 | 湍流k-ε | 层流1 | 层流2 |
D | 0.37 | 0.48 | 0.42 |
Mj_AVE | 4.479 4 | 4.498 7 | 4.500 8 |
Mj_MAX | 4.485 3 | 4.509 2 | 4.509 2 |
Mj_MIN | 4.443 6 | 4.475 3 | 4.489 3 |
σMj | 0.005 3 | 0.009 0 | 0.006 5 |
ΔMj_MAX | 0.035 7 | 0.023 4 | 0.011 5 |
而从高超声速风洞现有的喷管流场校测结果来分析, 实际加工出的M4.5喷管的流场均匀区直径在两者之间, 其为0.42 m是比较合理的, 层流模型计算出的结果取均匀区直径为0.42 m时, 平均马赫数为4.500, 均方差为0.006 5;这满足参考文献[5]中规定的先进指标要求, 同时也符合参考文献[6]中规定的指标要求。可见如果喷管加工按要求进行, 本套M4.5喷管型面数据是符合要求的。
3 结构设计马赫数4.5喷管共3段, 即喷管一段、喷管二段和喷管三段, 喷管与喷管之间轴孔配合, 外端法兰连接, 为保证喷管的气密性要求, 在喷管之间安装有紫铜垫圈, 通过螺栓压紧紫铜, 紫铜垫变形的方式实现喷管段的密封。喷管一段为气体入口段, 也是通过法兰与风洞设备实现连接, 密封方式也采用紫铜垫圈, 喷管三段与风洞设备连接采用轴孔配合的方式, 配合端面采用石棉垫, 端面外侧通过多个压板实现喷管段与风洞设备的连接固定。其剖面图见图 6, 三维图见图 7。
![]() |
图 6 喷管剖面图 |
![]() |
图 7 喷管三维图 |
考虑喷管实际情况, 对喷管进行强度校核, 分别对喷管一段、二段、三段添加载荷约束后进行计算, 最后将3段装配完成后进行计算, 计算应力云图如图 8所示。由图中可以看出, 喷管最大应力为28.285 MPa, 计算其安全系数为: η=660/28.285≈23.33, 满足技术要求。
![]() |
图 8 喷管应力云图 |
速度场校测排架呈“十”字形, 共计49个测点。测点探头用Φ4.0 mm×1.0 mm不锈钢管制成, 为平头型总压探头, 伸出支架前缘50 mm, 用Φ3.0 mm×0.5 mm紫铜管连接。
试验时总压排架通过专用测压支臂固定在5D机构中部支架上, 测压管路从风洞超扩段专用引出口引出后, 再用聚四氟乙烯软管与Netscanner-9016电子扫描阀系统相连。试验中通过控制5D机构的X轴来实现排架在喷管出口X=-0.15 m~0.40 m范围内运行, 从而测量出试验段内不同截面的速度场特性。
速度场校测时, 测点马赫数的得出是通过把测量出的前室总压P01和波后总压P02代入正激波关系式通过迭代计算得到的。正激波关系式为
![]() |
(1) |
式中,γ=1.4。
根据上面给出的方法, 速度场校测时先计算出单个截面上均匀区内马赫数的平均值, 然后利用每个单截面的计算结果就可以计算出整个流校空间均匀区内的参数。
4.2 校测结果本期速度场校测总压范围为0.3~1.5 MPa, 其中P0=0.6 MPa, T0=336 K的开车条件下, 排架在X=-0.15~0.40 m范围内移动时, 单独截面数据处理结果为: 均匀区直径D=0.40~0.28 m, 平均马赫数Mj_AVE=4.443~4.461, 马赫数均方差σMj=0.005~0.014, 马赫数最大相对偏差εMj_MAX=0.32%~0.92%, 满足国军标中规定的马赫数最大相对偏差≤1.0%指标。速度场分布见图 9。每一条曲线中间为图中左边具体马赫数, 相应跳到点则在给出值上加减偏差小格区域的值。
![]() |
图 9 速度场分布 |
取X=0.00~0.30 m、均匀区直径D=0.40~0.31 m时, 空间上的数据处理结果为: 轴向马赫数梯度dM/dX=-0.040 6 m-1, 平均马赫数M-AVE=4.456, 马赫数均方差σM=0.012 2, 马赫数最大相对偏差εM-MAX=0.90%, 满足国军标中规定的马赫数最大相对偏差≤1.0%指标, 可以用于型号试验。
取X=-0.15~0.40 m、均匀区直径D=0.40~0.28 m时, 空间上的数据处理结果为: 轴向马赫数梯度dM/dX=-0.037 7 m-1, 平均马赫数M-AVE=4.456, 马赫数均方差σM=0.013, 马赫数最大相对偏差εM-MAX=0.96%, 满足国军标中规定的马赫数最大相对偏差≤1.0%指标, 但是由于在X=-0.15 m截面, 偏差点偏大, 建议不适用此截面进行试验, 因此, 马赫数4.5喷管从X=-0.10~0.40 m截面的流场可用于型号试验。
5 模型试验 5.1 试验模型试验采用HB-2(Φ70)标模进行试验, 其外形尺寸见图 10, 模型是常用的高超声速弹道导弹式导弹标模, 由钝锥-柱-裙组成, 头部半径为Rn=0.3D;底部直径db=1.6D;长度lm=4.9D, 参考面积取模型中部截面积πD2/4, 力矩参考长度取模型中部直径D, 力矩参考点取模型轴线上距头部顶点的1.95D处, 压心系数的无因次参考长度取模型长度lm, 图中D=0.07 m。零攻角时模型在风洞中的堵塞度约为5%, 满足风洞流场建立要求。试验时, 模型表面气流流态采用自由转捩方式进行。
![]() |
图 10 模型外形简图 |
首先对不同马赫数的气动力系数进行分析, 在Ma=4.46, 4.95, 5.96, 6.97试验条件下的CN, 在攻角-2°~2°时, 法向力系数CN的斜率分别为0.079 6, 0.075 5, 0.069 9, 0.063 1, 符合常规条件下的马赫数规律。其结果图见图 11, CA, Cm, Xcp也有同样的数据规律。
![]() |
图 11 不同马赫数的CN分布图 |
然后进行了同期重复性试验, 表 3给出了气动力系数的均方根误差结果, 从表可以看出, 所有气动力系数的精度满足国军标指标要求。说明风洞各系统稳定、可靠, 设备运行正常, 可以开展型号、课题试验。其中a为攻角, σCN, σCA, σCm, σXcp分别为法向力系数、轴向力系数、俯仰力矩系数、压心的均方根误差。
a | σCN | σCA | σCm | σXcp |
-4 | 0.000 52 | 0.001 25 | 0.001 31 | 0.000 62 |
-2 | 0.000 48 | 0.000 90 | 0.00121 | 0.001 16 |
-1 | 0.000 33 | 0.000 72 | 0.000 49 | 0.001 02 |
0 | 0.000 40 | 0.000 50 | 0.000 39 | |
1 | 0.000 31 | 0.000 56 | 0.000 60 | 0.001 12 |
2 | 0.000 32 | 0.000 69 | 0.000 64 | 0.000 65 |
4 | 0.000 60 | 0.000 47 | 0.000 77 | 0.000 38 |
6 | 0.000 81 | 0.000 43 | 0.000 65 | 0.000 16 |
8 | 0.000 68 | 0.000 65 | 0.000 51 | 0.000 08 |
10 | 0.000 66 | 0.000 70 | 0.000 47 | 0.000 11 |
12 | 0.001 02 | 0.000 62 | 0.000 51 | 0.000 10 |
14 | 0.000 52 | 0.001 25 | 0.001 31 | 0.000 62 |
通过对高超声速风洞马赫数4.5喷管的内型面设计、流场计算、结构设计加工, 最后进行流场校测和模型试验, 可以得出以下结论:
1) 马赫数4.5喷管的流场均匀性满足GJB4399指标要求, X=-100~400 mm的流场均匀区可用;
2) 从HB-2标模数据的结果表明, 马赫数4.5与马赫数5, 6, 7的测力数据规律一致, 重复性精度达到GJB4399指标要求;
3) 综合上面结论, 本文采用的气动设计和结构设计方法合理, 研制的马赫数4.5喷管可投入试验应用。
[1] |
杨永能. FL-31风洞增设M4.5喷管总体方案[R]. 中国空气动力研究与发展中心, EJ-202000095 YANG Yongneng. Overall scheme of adding Mach number 4.5 nozzle of FL-31 wind tunnel[R]. CARDC-EJ-202000095(in Chinese) |
[2] |
黄飓. FL-31风洞马赫数4.5喷管研制实施方案[R]. 中国空气动力研究与发展中心, EJ-20210665 HUANG Ju. Mach number 4.5 nozzle development and implementation scheme of FL-31 wind tunnel[R]. CARDC-EJ-20210665(in Chinese) |
[3] |
杨永能, 黄飓. FL-31风洞M4.5喷管型面设计计算[R]. 中国空气动力研究与发展中心, EJ-202000096 YANG Yongneng, HUANG Ju. Mach number 4.5 nozzle profile design calculation of FL-31 wind tunnel[R]. CARDC-EJ-202000096 |
[4] |
张凯钊. FL-31风洞增建马赫数4.5喷管方案论证报告[R]. 中国空气动力研究与发展中心, EJ-20200096 ZHANG Kaizhao. Demonstration report of Mach number 4.5 nozzle of FL-31 wind tunnel[R]. CARDC-EJ-20200096(in Chinese) |
[5] |
黄飓. FL-31风洞M4.5喷管流场校测和标模试验报告[R]. 中国空气动力研究与发展中心, EJ-2021100141 HUANG Ju. Report of Mach number 4.5 nozzle flow calibration and standard model test in FL-31 wind tunnel[R]. CARDC-EJ-2021100141(in Chinese) |
[6] | HAGEMAN G, I MMICH H. Critical assessment of the linear plug nozzle concept[R]. AIAA-2001-3683 |
[7] | Maureen Rowbo tham. XR S2200 linear aerospike engine: Use of proengineer for determing mass properties[R]. AIAA-1999-2334 |
[8] | BOOTH T E, VILJA JO, CAP D P, et al. The design of linear aerospike thrust cells[R]. AIAA-1993-2562 |
[9] | MEENAH R, HOFFMAN J D, MURTHY S N B. Design and performance computations in complex 3D nozzles[R]. AIAA-1999-0882 |
[10] | VAN L B. Towards the ultimate conservative difference scheme V: a second order sequel to Godunov's method[J]. J of Comp Phys, 1979, 32: 101-136. |