探空火箭减阻杆气动特性分析常耀予刘帆张家齐摘要：为了研究减阻杆对探空火箭气动力特性，通过采用SST两方程湍流模型、有限体积法求解N-S方程，对探空火箭高速流场进行数值模拟。计算结果显示，减阻杆能有效减小火箭阻力。亚跨声速（Ma0.8～1.2）最大减阻25%；高超声速（≥Ma6）阶段，最大减阻量35%，减阻效果随迎角增大而降低，到12°迎角时减阻量为12%。压跨声速及高超声速全箭升阻比增量随马赫数增大均增加，高超声速阶段升阻比增大18%。同时采用工程方法结合数值预示结果，评估减阻杆带来的气动热影响，结果显示，气动支杆的存在使得端头的平均热流密度下降了51%，并与飞行试验结果进行对比分析。计算得到热流结果与飞行实测热流结果相当，热环境预示比较准确，对于高超声速阶段的飞行器被动热防护技术研究具有良好的指导价值。关键词：减阻杆；气动热；气动力；高超声速；探空火箭：V211.73：ADOI：10.19452/j.issn1007-5453.2020.11.010近半个世纪，高超声速流动机理及控制方面较低速流动领域取得的进展较少[1]，高超声速流动控制主要集中在减阻防热控制方面，控制方法为在钝头体头部安装减阻杆进行流场重构，减小气动阻力[2]。减阻杆属于被动流动控制，通过前伸支杆外推钝头体头部弓形激波，转换为强度较弱的锥形斜激波，在头部形成低压回流区，减小激波阻力，降低表面热流。减阻杆以其显著的减阻和降热效果，在国外导弹上已经实现工程应用。美国飞行速度Ma24的三叉戟Ⅱ型彈道导弹（UGM-133Trident-II，1990年服役）采用伸缩式尖针状减阻杆[3]。导弹出水发动机点火不久，加速度触发伸出减阻杆，到位后锁定，这种设计便于导弹头部在保证装填空间的前提下，同时兼顾高速飞行时的低阻特性，导弹在高超声速情况下减阻达52%，同时导弹射程增加了550km[4]。由于减阻杆外形为针状结构，尖针杆会带来防热和流场不稳定问题[5-6]，本文研究了“球形头+圆锥”减阻杆构型下[7]，探空火箭在高超声速下的气动力特性与气动热特性。通过仿真手段对比分析了气动力减阻效果，同时采用工程算法及飞行试验对减阻杆带来的降热效果进行评估，为后续工程应用提供设计依据。1计算方法1.1气动力计算方法2探空火箭外形火箭采用轴对称气动布局，任务飞行器采用“球头+圆锥”气动外形，球头前端布置减阻杆；控制舱、发动机及尾段为柱段形式；空气舵和燃气舵共轴布置，均采用X形布局。火箭气动布局如图1所示。根据目前对减阻杆的气动布局特性研究，减阻杆采用“端头帽+细长杆”布局[8]，设计参数主要有减阻杆长度、端头直径及减阻杆直径，如图2所示。根据图3减阻降热的计算结果，优化设计，选取减阻杆长度为锥段球头直径的两倍，端头帽为球头直径的0.3倍，减阻杆直径为球头直径的0.1倍[9]。根据火箭几何参数，减阻杆端头半径45mm；圆柱长度555mm，直径30mm；减阻杆后部的锥段前缘半径150mm，半锥角12°，如图4所示。3减阻杆气动特性结果分析3.1减阻杆气动力分析气动力计算网格采用混合网格，为准确模拟边界层内黏性对气动力的影响，采用Octree八叉树法构建空间四面体网格，在壁面沿法向生成9层棱柱形网格，头部及尾段采用密度盒加密，总网格量1145万。此种网格划分方法、网格量已经在OS-X0，OS-X1两个型号飞行试验中进行验证，数据准确可靠，如图5所示。仿真计算设置压力远场边界条件，选择理想气体介质，计算马赫数从Ma0.6开始到发动机关机点Ma6，典型马赫数状态流场的马赫数和压力云图[10-11]如图6和图7所示。支杆顶部聚集起高压区，有效降低了端头附近的压力，从而降低了端头阻力。对减阻杆的减阻效果进行评估，亚跨声速（Ma0.8～1.2）最大减阻25%；高超声速（Ma6）阶段，最大减阻量35%，减阻效果随迎角增大而降低，到12°迎角时减阻量为12%。压心位置飞行全程变化量在1%以内，马赫数6小迎角状态压心前移3%，减阻杆对纵向压心的影响如图8所示。亚跨声速及高超声速全箭升阻比增量随马赫数增大均增加，高超声速阶段升阻比增大18%[12-14]，如图9所示。3.2减阻杆气动热分析根据工程算法，采用减阻杆后，支杆及端头的热流密度随时间的变化曲线如图10所示。支杆顶点最大热流密度1738kW/m2，端头最大热流密度9...