火箭橇试验战斗部动态攻角交互式计算与验证

第40卷第2期2020年4月
Vol.40No.2
Apr2020弹箭与制导学报
Journal of Projectiles,Rockets,Missiles and Guidance
DOI:10.15892/ki.djzdxb.2020.02.022
火箭橇试验战斗部动态攻角交互式计算与验证
张晨辉，徐健
(中国兵器工业试验测试研究院，陕西华阴714200)
摘要：针对火箭橇试验战斗部攻角控制问题,提出了一种基于气动环境分析的耦合交换数值计算方法，采
用自编程与CFD数据交换耦合方式，解决了火箭橇战斗部动态攻角工程计算难题,在设计、计算与优化的基
础上，对3种方法的特点及适用范围进行火箭橇试验验证。

该研究结果可用于精确、高效战斗部火箭橇试验
攻角设计与控制，可以为导弹引战系统终点效应试验及相关技术研究提供一定参考和技术支持。

关键词：战斗部；攻角；火箭橇试验；CFD；龙格-库塔方法
中图分类号：TJ760.6文献标志码：A
Rocket Sled Test Warhead Dynamic Attack Angle Interactive
Calculation and Verification
ZHANG Chenhui,XU Jian
(Test and Measuring Academy of China Ordnance Industries,Shaanxi Huayin714200,China)
Abstract:Aiming at the problem of rocket sled test warhead attack angle control,we put forward a coupling exchange numerical method based on aerodynamic environment analysis,which exchange data between self-programming and CFD,this method solve the engineering calculation problem of rocket sled warhead attack angle,and on the basis of design,calculation and optimization,we test and verify the characteristics and application of the three methods.The research can be used in precise and efficient warhead rocket sled test attack angle design and control,and it provides a certain reference and technical support for the missile fuze-warhead system end-point effect test and the related technology research.
Keywords:warhead；attack angle；rocket sled test；CFD；Runge-Kutta method
o引言
导弹战斗部终点弹道的攻角大小直接影响其侵彻能力，例如，与0。

攻角侵彻状态相比,单一动能侵彻战斗部攻角为3。

时，侵彻效应降低约5%,当攻角为
9。

时，侵彻效应降低约&7%;对于随进侵彻战斗部，初始攻角对其终点效应影响相对更大，当攻角为3°
时，侵彻效应降低约7%,当攻角为9。

时，侵彻效应降
低约30%，理论上撞击靶标时刻攻角越小越有利于
侵彻效果
火箭橇试验是导弹战斗部终点效应考核的一种
重要地面动态试验与验证方法,该方法运用高精度专
用地面滑轨设施，通过固体火箭发动机实现系统加
速，当达到指标速度时，利用战斗部与火箭橇时空分
离方式,使战斗部按预定要求的速度和攻角单独侵彻
目标，但随着导弹武器系统可靠性指标的不断提高,
相关型号研制试验中暴露出过小的控制战斗部初始
攻角并不能完全模拟出导弹引战系统的使用极限条件，无法实现对被试导弹武器系统的全面考核,因此,提出了在战斗部着靶时刻必须带有一定攻角的要求。

与导弹空中自由飞行不同,火箭橇试验中战斗部分离后飞行气动力学环境综合了战斗部和火箭橇的马赫数、战斗部姿态、三维相对时空位置、压心突变、局部激波振荡等变化的耦合作用结果,近年来,随着相关研究的推进,在超音速翼型火箭橇气动特性⑵、超音速火箭橇流场气动数值计算分析、战斗部攻角对升力系数、阻力系数及升阻比的影响⑸等方面取得了突破性研究成果,但是，在火箭橇试验战斗部动态攻角计算及控制中,需要将复杂的多维空间气动力学插值与战斗部6D飞行姿态计算相结合，相关数学建模、仿真网格划分和耦合求解都十分困难。

针对火箭橇试验中战斗部动态攻角计算和控制难题,结合大量火箭橇试验测试数据,文中采用自编程与CFD数据交换耦合方式，提出了一种交互式的火箭橇战斗部动态攻角工程计算方法,并结合验证试
*收稿日期：2019-04-16
基金项目:国防技术基础科研项目资助
第一作者简介:张晨辉（1988-）,男，山西长治人，助理研究员，硕士，研究方向:火箭橇试验CFD计算。

・90・弹箭与制导学报第40卷验,对计算精度进行了相关验证。

1火箭橇运动分析
1.1控制方程
流体分析时，流体的运动要符合经典力学的3个
基本守恒定律，即质量守恒定律、动量守恒定律和能
量守恒定律。

引入雷诺平均法后控制方程⑹如下：
连续方程舞+盏伽）二0
运动方程:
d（pu旳）__dp dt+dXj=一dX「
艇淞警+嚮二2（r d0
dxy込
式中“J二1,2,3；&和S分别是运动方程和输运方程的广义源项亦为广义扩散系数。

雷诺应力-p丽j属于新的未知量，对应3个正应力与3个切应力，共6个不同的雷诺应力项，因此，加入雷诺应力后,方程中共有10个未知量。

方程组不封闭，需要引入其他湍流模型联系时的均值与脉动值,使方程组封闭。

1.2气动仿真
三维复杂外流场气动仿真是以流体力学与固体力学理论为基础,采用大型计算流体动力学分析软件（CFD）进行模拟,仿真数据为结构气动特性分析提供基础。

现阶段国内外物体外流场的数值仿真研究中,通常应用N-S方程耦合双方程湍流模型对湍流现象进行模拟，而在双方程湍流模型中基于模型的SST （剪切应力输运）湍流模型⑺以其优势得到广泛应用:①该模型能适应压力梯度变化的各种物理现象;
②该模型中可应用粘性内层，可精确模拟边界层的现象，无需使用较容易失真的粘性衰减函数。

研究表明，采用SST模型进行外流场湍流模拟结果与风洞试验结果基本吻合,故采用SST模型对中文火箭橇复杂外流场问题进行研究。

1.3运动计算基本原理
龙格-库塔（Rrnige-Kutta）方法屈是用于非线性常微分方程解的重要的一种隐式或显式迭代法，该方法有如下优点:1）具有同一型式的计算流程图;2）不需要进行“起始”计算;3）在计算过程中很容易改变步长;4）能很好的适应方程右端为间断函数的情况。

四阶龙格-库塔法的截断误差为0（丽）。

程序框图如图1所示。

|开始］
|读入数据|
2战斗部动态攻角交互式计算方法
2.1气动力学环境分析
根据弹橇分离后战斗部的速度、姿态、时空位置和局部激波等因素，分析战斗部飞行过程的气动力学环境变化,得出战斗部气动阻力、升力和力矩曲线。

利用流体软件ANSYS CFX进行火箭橇流场计算,合理设置边界条件。

火箭橇模型表面采用无滑移的固壁边界条件，滑轨及地面采用滑移边界条件，采用SST（剪切应力输运）湍流模型，采用非结构网格求解可压缩雷诺平均N-S方程对火箭橇复杂外流场问题进行计算。

2.2耦合交换数值计算
结合流体软件ANSYS CFX和基于四阶龙格-库塔插值与牛顿运动定律的VB自编程序，建立战斗部攻角计算方法，形成了更为科学的计算理论。

1）结合导弹引战系统终点效应火箭橇试验气动力学环境分析与计算,进行数学建模,弹橇分离点与目标距离段的弹道数学模型的坐标轴为:选取弹橇分离时刻弹的质心位置为坐标原点,水平轴0%轴与滑轨平行，选取火箭橇飞行方向为正;沿直轴oy在弹纵向对称面内与水平轴0%相垂直，向上为正;轴垂直于。

矽面，翻转力矩以使被试品系统抬头的力矩为正,以使被试品低头的力矩为负，该坐标轴的建立基于实际的火箭橇试验要求。

初始求解程序数学模型如式（1）所示。

第2期张晨辉等:火箭橇试验战斗部动态攻角交互式计算与验证•91•
=F
Av
m不=S-mg
(1)
Ut M
式中-.m为战斗部质量（kg）；%为航向速度（m/s）-,v y 为纵向速度（m/s）茫为战斗部攻角角速度（rad/s）；%为航向位移（m）-y为纵向位移（m）；a为攻角（rad）; J z为战斗部赤道转动惯量（kg•m2）疥为航向空气阻力合力（N）；S为纵向空气阻力合力（N）；M为翻转力矩（N•m）。

2）为了仿真计算火箭橇在滑轨上的运动过程，得到火箭橇系统各部件运动轨迹,根据火箭橇结构参数建立火箭橇运动模型，用VB语言编程，运用四阶龙格-库塔对控制方程求解，进行初始化求解程序编制。

3）在初始化计算的基础上,与气动分析前处理软件ICEM进行战斗部分离后飞行气动力变换数据。

4）同时进行网格再划分CFD分析及后处理等,通过多次迭代计算，得出火箭橇试验战斗部弹橇分离后飞行姿态与攻角变化。

此方法已多次应用于火箭橇试验战斗部攻角控制中，并得到充分验证。

此方法依据一次初算结果,构建多间隔点的数值模型，可并行计算，但战斗部飞行距离长的情况下角速度的累积偏差导致各插值点的攻角将出现较大偏差,故计算时间短、精度较差,适合于转动惯量较大或飞行距离较短的战斗部攻角估算，在工程计算中应用较为广泛。

2）分段逐步外推方法
同理,先计算出弹橇分离起始点初始参数,结合弹橇分离点与靶标距离段的弹道模型,将该段距离均分N段，利用VB自编程推算在该气动力环境下运行至第1段段末位置时战斗部攻角等参数，同样建模计算出第1段段末位置参数，利用四阶龙格-库塔插值进行外推,得出第2段段末位置参数，以此类推，得到各位置战斗部攻角等参数。

此方法需要分段计算，需多次建模、多次插值计算，但这种方法推导出的每位置点攻角等参数更准确，故计算时间较长、精度较高,适合于转动惯量较小且精度要求高的战斗部攻角计算,耗时长，不适合于工程计算。

3）区间插值方法
此方法是第一种方法的演变,针对飞行距离较长情况将其分为2个或3个短距离的区间，在区间内采用一次初算多点插值方法，将上一区间末的输出结果作为下一区间的输入结果，最终获取战斗部攻角的变化趋势。

此方法对第一种方法进行了优化,第二区间或第三区间战斗部的攻角模拟更加准确，故计算时间较短、精度较高，适合于转动惯量较小或飞行距离较长的战斗部攻角计算。

3计算方法分类
根据计算与建模的不同，对上述动态气动力学环境数值计算方法进行分类,并阐述其特点及适用范围。

1）一次初算多点插值方法
利用ANSYS CFX对起始点气动力学环境进行分析,得到起始点战斗部及橇体阻力/升力、战斗部翻转力矩等初始参数，通过VB攻角计算程序，推算出弹橇分离过程中每点战斗部速度、攻角、航向飞行距离、角速度、橇体倾斜角度、航-纵向拉开距离等参数，据此以一定的距离间隔构建多个点的力学模型,分别对各力学模型进行气动特性分析,得到各间隔点战斗部及橇体力学参数，进行VB自编程与CFD数据交换耦合计算求解得到战斗部攻角变化曲线。

4火箭橇试验验证
下面通过火箭橇试验案例对3种计算方法特点及适用性进行分析,战斗部由模拟件代替，橇体进行了局部简化处理，建立模拟件火箭橇仿真物理模型,假设模拟件转动惯量及飞行距离等参数。

案例一:火箭橇结构外形如图2所示，模拟件以预置攻角0。

安装在橇体上，假定弹橇分离速度300
m/s,模拟件质量500kg,火箭橇橇体质量600kg,模拟件质心1400mm（距后端），模拟件转动惯量200
kg•n/,空气密度1.067kg/n?,飞行距离设定为30m o
下面分别采用一次初算多点插值方法及区间插
值方法，获得模拟件的攻角变化。

特征位置模拟件及
橇体表面压力云图如图3所示。

/模拟件更快捷地得出攻角变化趋势。

4
—-一次初算多点插值方法
T-区间插值方法
a
图2火箭橇结构外形图0^-———1-----------------1--------1-------->-
0510********
飞拟件飞行距离/m
图4两种方法攻角变化趋势对比
案例二:上述参数不变的情况下,为对比一次初
算多点插值方法及分段逐步外推方法,假定模拟件转
动惯量为100kg・m2,获得两种情况下模拟件的攻角
变化趋势，如图5所示。

(a)弹橇分享时刻
(。

(b)过弯轨时刻
00,500l.O OO(m)
"™50^^50
(c)着靶时刻
图3模拟件及橇体表面压力云图
如图4所示,采用一次初算多点插值方法及区间插值方法计算得出着靶攻角分别为4.055。

、3.774°,根据大量的试验测试数据进行反演推算，实际试验中着靶攻角为3.5。

区间插值法将30m飞行距离分为12m和18m两个区间，在第一区间，攻角变化趋势与一次初算多点插值方法一致，在第二区间，由于区间插值法获取的模拟件各时刻的攻角较为准确，攻角发散程度较小，计算精度较高，计算偏差为0.281°,1+算时间增加一倍。

由于案例中模拟件转动惯量较大,两种方法同一时刻模拟件攻角偏差不会引起较大的翻转力矩差异,故区间插值方法优势不明显;若飞行距离缩短为20m,两者偏差会更小，仅为0.03。

因此,区间插值法适合于转动惯量较小或飞行距离较长的攻角计算，反之，采用一次初算多点插值方法可以
图5两种方法攻角变化趋势对比
如图5所示，由一次初算多点插值方法及区间插
值方法计算得出的着靶攻角分别为7.275。

和
5.647°,根据大量的试验测试数据进行反演推算，实 际试验中着靶攻角为4.9。

由于分段逐步外推方法
能够精确模拟出各时刻力学参数，计算精度提高了
22%，但由于计算中进行了5次分段计算，故计算时间增加5倍,不适用于工程计算。

对于模拟件飞行距
离短的情况，譬如假定飞行距离20m,计算方法引起
的偏差不会太大,计算偏差仅为0.64。

5结论
1)文中针对火箭橇试验中战斗部攻角控制，开发了一种基于气动环境分析的耦合交换数值计算手段,
解决了火箭橇战斗部动态攻角工程计算难题。

2)文中对计算方法分类、特点及适用范围进行介绍,并通过火箭橇试验案例分析其计算精度及计算时
间,一次初算多点插值方法计算时间短、精度较差,
适合于转动惯量较大或飞行距离较短的战斗部攻角
估算,在工程计算中应用较为广泛;分段逐步外推方
法计算时间较长、精度较高，适合于转动惯量较小且
精度要求高的战斗部攻角计算，但耗时长，不适合于
(下转第96
页)
降低，提高了测量的精度，此时测量的平均绝对误差
图7进行畸变校正后的双目测距结果4结语
文中深入研究了单目测距原理和双目测距原理,为双目测距系统的目标距离测量提供了一种基于计算机视觉的计算方法，并进行了验证实验，实验环境下单目测距的平均相对误差为3%，双目视觉测距的平均绝对误差为0.766cm,获得了符合精度要求的结果。

参考文献：
[1]丁欢欢.双目立体视觉测距系统关键技术研究[D].成
都：电子科技大学,2015.
[2]ZHOU H L,WEI L,CREIGHTON D,et al.Orientation a-
ware vehicle detection in aerial images[J].Electronics
Letters,2017,53(21)：1406-1408.
[3]SATZODA R K,LEE S,LU F,et al.Vision-based front
and rear surround understanding using embedded processors
[J].IEEE Transactions on Intelligent Vehicles,2017,1
(4)：335-345.
[4]周士超.无人机电力巡检视觉避障技术研究[D].北京：
北京理工大学,2016.
[5]裴聪•基于计算机视觉中双目立体匹配技术研究[D].
镇江:江苏大学,2010.
[6]尚倩.基于双目立体视觉的目标识别与定位[D].北京：
北京交通大学,2011.
(上接第92页)
工程计算;区间插值方法计算时间较短、精度较高，适合于转动惯量较小或飞行距离较长的战斗部攻角计算。

该研究结果可用于精确、高效战斗部火箭橇试验攻角设计与控制。

参考文献：
[1]张学伦，皮爱国,刘宗伟，等•攻角影响战斗部侵彻效应数
值模拟分析[J]•弹箭与制导学报,2012,32(6)：71-73. [2]张传侠，吕水燕,谢波涛，等.强地效环境下有翼火箭橇
侧翼气动特性研究[J].兵器装备工程学报,2018,39
(9)：89-92.
[3]张立乾，邓宗才,陈向东.超音速单轨火箭滑橇气动特性
数值模拟[J]•弹道学报,2011,23(4)：100-104.
[4]LOFTHOUSE A J,HUGHSON M C,PALAZOTTC A ­
putational aerodynamic analysis of the flow field about a hyper­
velocity test sled：AIAA2003-981[R].Redston:AIAA
Inc.,2003.
[5]张红英，孙姝，程克明，等.不同飞行工况下的某高超飞
行器流道特征及气动力特性研究[J].航空动力学报, 2008,23(12)：2262-2267.
[6]黄银柳，陶如意，王骁，等.某异形弹丸的气动特性研究
[J]•兵器装备工程学报,2018,39(4):42-45.
[7]胡朋，李永乐,廖海黎.基于SST k-co湍流模型的平衡大
气边界层模拟[J]•空气动力学学报,2012,30(6):737-
743.
[8]游培寒，胡瑜，缪昕，等•超音速悬索火箭橇下拉幅度运
动仿真分析[J]•弹箭与制导学报,2017,37(5)：147-
150.
[9]沈艳，张丽玲,张琦智，等.基于三阶和四阶龙格库塔法
的GM(1,1)模型优化及应用[J].数学的实践与认识, 2016(7)：168
-173.。