主动气动起旋离心抛撒技术具备结构简单、制造成本低、子母弹装填数量高、抛撒速度可控、子弹药之间不易产生相互干涉碰撞等优良特性,因此,得以在美国CBU-103型(见图1)、SUU-65/B型等航空子母炸弹系列产品中得到大规模的成熟应用[1],开展该项技术的研究对于我国的航空弹药和导弹工业而言具有重要的战略意义。

图1 美国CBU-103型制导航空子母炸弹

近年来,气动/运动一体化数值模拟方法在国内外得到广泛的应用。美国于2008年研制了飞行器计算模拟系统Kestrel,该系统将CFD求解器、结构动力学方程、刚体六自由度运动方程相耦合,且具备舵偏控制技术、动网格技术以及嵌套网格技术,能较准确地预测计算模型的操纵特性。国内相关领域学者通过对非定常N-S方程和六自由度运动方程进行同步耦合求解,并在求解过程中引入了Adams预估-校正法,采用动态网格技术对面对称飞行器加入舵偏控制的大机动过程进行数值模拟,验证了气动/运动紧耦合求解计算方法的可行性与准确性,计算结果和靶式数据较为一致,证明了该方法可以准确模拟飞行器“自由飞”过程[2]。

1 主动气动起旋过程数值模拟

母弹在空中的主动气动起旋是离心抛撒技术中的核心环节,它要求航空炸弹在空中实现飞行状态的平稳迅速切换,即从平飞状态切换到以高速滚转为主特征的飞行状态。掌握母弹的空中主动起旋技术是实现离心抛撒技术的必要前提,本节根据母弹的气动特性对弹体在空中主动气动起旋过程进行了数值模拟研究。

1.1 计算模型

本文以航空子母弹药为研究对象,计算模型采用轴对称气动布局,弹体头部为球缺型,并加装4片×型布局前翼(反安定面);弹身为圆截面体,制导控制尾舱以4片×型全动尾舵作为操纵面。全弹通过滚转舵偏实现“剪刀差”的主动气动起旋过程,滚转舵偏示意图如图2所示,图中,δx为滚转舵偏角。

图2 滚转舵偏δx=20°,30°,40°时计算模型示意图

1.2 数值计算方法

1.2.1 非定常流场计算

复合型人才的培养离不开高素质的师资队伍，培养复合型信息与计算科学人才，对师资队伍提出了很大挑战．需要制定师资队伍建设新方案，针对高素质复合型人才培养的要求，构建适应复合型人才培养的师资结构．一方面，要加强对骨干教师的专业培训，同时制定激励措施，鼓励专业教师到企业一线锻炼或有计划地选派有潜质的青年教师直接参与企业的项目实践；另一方面，聘请校外一些优秀的数据挖掘工程师、软件工程师、实践经验丰富的培训讲师等兼职任教、合作研究，扩大复合型人才培养的师资力量．

在曲线坐标系下,完全气体的三维无量纲化非定常N-S方程的守恒形式为

(1)

式中:Q为求解矢量;F,G,H分别为一般曲线坐标系下x,y,z 3个方向的无黏矢通量;Fv,Gv,Hv分别为3个方向的黏性矢通量[3]。

采用基于非结构网格的有限体积法进行空间离散,用近似通量差分格式和中心差分格式离散无黏通量和黏性通量。在非定常计算中,采用伪时间步为LU-SGS隐式计算格式的双时间步计算方法。在边界条件的设置中,采用局部一维Riemann不变量进行外域远场设置,以相对无滑动的固定壁面作为弹体表面边界条件,同时采用多重网格方法和MPI并行计算方法来提高收敛速度和计算效率[3]。

(1)分别于插管后5 min(T1)、建立气腹后5 min(T2)、体位建立后5 min(T3)、建立气腹后30 min(T4)、解除气腹后5 min(T5)时以旁流通气监测法检测2组呼吸力学指标，包括气道压(Paw)、气道峰压(PIP)，并计算肺泡动态顺应性(Cdyn);(2)于术前24 h及术后24 h时，采集2组患者清晨空腹静脉血5 mL，3 500 r/min离心10 min分离血清，采用ELISA法测定血清白细胞介素-6、10(IL-6、IL-10)及C反应蛋白(CRP)水平，免疫透射比浊法检测2组患者血清C3、C4补体水平;比较2组患者HR、SpO2及平均动脉压(MAP)。

1.2.2 运动方程求解及气动/运动耦合

1.4.2 滚转舵偏对起旋过程的影响

(2)

式中:为模型在弹体坐标系中的滚转角速度,为模型所受的滚转力矩,γ为模型的滚转角,上标“d”表示模型体轴系[4]。在非定常流场计算过程中,采用三阶Adams预估校正法同步求解运动方程,在确保耦合计算精度的基础上,对气动/运动进行紧耦合求解[5]。

1.3 计算网格方法

为实现尾舵舵偏形成的“剪刀差”外形的准确模拟,采用四面体非结构网格作为计算网格。同时,对模型表面曲率较大之处以及边界层区域引入各向异性非结构网格[6],网格划分情况见图3。

车间在需求申请赋码阶段计量单位基于车间的实际定额需求确定，比如油漆m，计量单位可以是升，也可以是kg；比如绑扎带，计量单位可以是个，也可以是m；而后期采购环节对于计量单位的要求，往往不完全与定额管理计量单位一致。

图3 全弹计算网格图(δx=20°)

本文采用刚性运动网格技术[6]模拟全弹整体运动过程,令计算网格随弹体一起做整体刚性运动,根据运动参数直接给出计算网格,计算量小,并且可以保持初始网格质量,计算精度较高[6]。

1.4 仿真结果与分析

弹体尾部的4片尾舵在旋转过程中同时具有驱动作用和主要的阻尼作用,因此需要研究尾舵附近流态在旋转过程中的变化。图7～图9为不同舵偏下t=0 s和t=0.5 s时刻的流线分布。可以看到,在初始时刻,尾舵背风面有较大面积的分离区,舵偏角越大,分离程度越强。与之相比,t=0.5 s时,由于存在一定的自旋速度,尾舵的当地迎角降低,背风面的流动分离大大减弱,这与图6的力矩变化趋势所揭示的规律是一致的。

文中介绍的分析方法是基于静止轨道同步卫星，对于其他具有透明转发器的卫星此方法也适用。另外，文章仅考虑了多普勒频移基本项，而没有考虑相对论影响，如引力频移和二阶多普勒频移，以及空间环境对本振的影响，有待继续研究。

1.4.1 静态滚转舵效分析

定常流场计算是开始动态数值计算的基础。图4为不同舵偏角滚转力矩系数数值计算与风洞试验对比曲线图,图中,Cl为滚转力矩系数,δx为滚转舵偏角。如图4所示,δx≤20°时计算与试验结果吻合较好,相对偏差约1.5%;由于计算方法无法准确模拟大舵偏角时的舵面失速现象,滚转舵偏越大计算误差越大。同时,前置小翼对滚转舵效有一定影响,对δx=30°的情况,去掉前置小翼后舵偏产生的滚转力矩系数增大约5%。

在设计相关的电气产品时，充分利用智能化技术能够使电气产品得到更好的改进。所以，在设计电气产品的过程当中，对于相关的影响因素便可以通过应用相关的智能化技术加以消除，替换人工操作的方式，充分发挥计算机技术应有的作用，这样就能够有效节省人力，使电气工程的设计工作具有更高的效率与质量。就现阶段而言，使用较为普遍的智能化设计有专家系统和遗传算法两种，其中，遗传算法基本上就是对操作的相应对象加以直接的控制，能够在很大程度上增强产品在运行方面的能力。而专家系统则主要是此领域中的专家基于自身的相关经验而进行推断，其都可以使电气产品得到有效的改进。

图4 不同舵偏角滚转力矩系数CFD计算与试验对比

本文仅考虑在弹体坐标系下策略模型绕弹体纵轴的滚转运动,对应的运动学方程与动力学方程如下:

图11给出了不同开舱条件下X向射程随抛撒高度的变化曲线,图中,θ为弹道倾角,L为射程,h为抛撒高度。

图5 不同舵偏角时转速随时间变化历程

图6为母弹起旋过程中不同舵偏角时滚转力矩系数随转速变化历程。如图6所示,随着转速的增加,相同转速下不同舵偏产生的滚转力矩的差量增大,即大舵偏角的滚转舵效随着转速增大而增大,这是由于转速增大时尾舵的当地迎角减小,因此尾舵背风面的分离情况会有一定的改善。

图6 不同舵偏角时滚转力矩系数随转速变化历程

1.4.3 初、末时刻尾舵附近的流态变化

通过对全弹进行气动/运动一体化数值模拟,并结合风洞静态测力试验结果,分析不同滚转舵偏(δx=20°,30°,40°)对主动气动起旋过程的影响。

综上所述,尾舵偏转产生的滚转力矩是影响转速的主要因素,同时由于大舵偏时存在舵面失速,滚转舵偏角在30°以上继续增大时,转速增加效果不再明显。为加快模型的自旋过程,增加尾舵展长是最行之有效的办法。这是因为增大展长能够有效提升尾舵效率并增大滚转力矩的作用力臂,同时,应在保持纵向操纵性能的前提下,尽量减少弹体表面的其他气动阻尼面。

图7 δx=20°时初、末时刻尾舵附近流态

图8 δx=30°时初、末时刻尾舵附近流态

图9 δx=40°时初、末时刻尾舵附近流态

2 离心抛撒子弹药弹道计算

为了验证主动气动起旋离心抛撒子弹药的性能,获取影响子弹药离心抛撒效果的母弹开舱条件,本文以小型杀爆子弹药作为研究对象,通过对不同开舱条件下子弹药弹道计算结果的对比分析,获得影响子弹药散布效果的规律。

式中，C0(mg·L-1)和Ce(mg·L-1)分别为亚甲基蓝溶液的初始浓度和平衡浓度；m(mg)为吸附剂投加量；Qe(mg·g-1)为平衡吸附容量。

于是始胠其橐，得故先君子时所藏与熹后所增益者，凡数十种，虽不多，要皆奇古可玩，悉加标饰，因其刻石大小，施横轴，悬之壁间，坐对循行卧起，恒不去目前，不待披筐箧卷舒把玩而后为适也。盖汉魏以前刻石制度简朴，或出奇诡，皆有可观，存之足以佐嗜古之癖，良非小助。其近世刻石本制小者，或为横卷若书帙，亦以意所便也。盖欧阳子书一千卷，赵氏书多倍之，而予欲以此数十种者追而与之并，则诚若不可冀，然安知积之久，则不若是其富也耶！姑首是书以俟。

图10 子弹药在弹仓内布局示意图

2.1 计算模型

在各飞行阶段中子弹药的质量不变、仅受重力和空气阻力作用的质点受力系的简化条件下,子弹药三轴向弹道运动微分方程组如下:

母弹离心抛撒子弹药的开舱条件主要由母弹转速、开舱高度、飞行速度、弹道倾角和目标海拔高度等参数组成。子弹药的散布效果主要由子弹药散布中心射程、散布面积、散布图形均匀性等指标体现[7-9]。母弹的旋转方向对离心抛撒子弹药没有影响,由主动气动起旋过程计算可知,采用不同滚转舵偏角,母弹在主动气动起旋后0.5 s可达到的转速绝对值约1 000～2 000 r/min的转速,且弹体姿态平稳,此时子母战斗部舱开舱,子弹药受离心力作用将以不同的速度和角度抛撒出舱。设定母弹含有12个弹舱,每个弹舱装填18枚子弹药,共216枚子弹药。子弹药在弹仓内的布局如图10所示。

(3)

式中:x,y,z;vx,vy,vz;Sx,Sy,Sz;Cx,Cy,Cz分别代表子弹药在弹道坐标系OXYZ三轴向的坐标位置、飞行速度、有效阻力面积、阻力系数函数;v0为子弹药初始速度;m为子弹药质量;n为母弹转速;r为母弹弹体半径;ρ为子弹药所在高度上的空气密度。采用四阶龙格库塔法求解弹道方程[10-11]。

2.2 计算结果与分析

2.2.1 子弹药散布中心射程分析

2.3 酒糟有机肥不同用量对烤后烟叶化学成分的影响 从表3可以看出，各处理中部烤后烟叶主要化学成分指标除钾外均无显著差异，含有酒糟有机肥2 250 kg/hm2配施20%化肥的处理钾含量显著低于其他处理，其糖和淀粉含量均高于其他2个处理；3个处理的氯含量较低，均低于对照，其他化学成分含量各处理间以及各处理与对照间均无显著差异。

图5为母弹起旋过程中不同舵偏角时母弹转速随时间变化历程,图中,t为时间;n为母弹转速,母弹的旋转方向由右手法则定义,在正滚转舵偏的情况下,母弹的旋转方向为负,转速对应为负值。由图5可知,在起始阶段转速迅速增加,同时滚转力矩逐渐降低;随着时间的推移,转速和滚转力矩的变化速度略有减缓。到t=0.5 s时,δx=20°状态时转速可达到-1 136.47 r/min;δx=30°时转速相对提高约45.2%,达到-1 650.67 r/min;δx=40°时转速为-1 985.22 r/min,相较δx=20°时转速相对增大约74.7%。

图11 v=300 m/s,n=2 000 r/min,不同弹道倾角时抛撒高度与X向射程关系曲线

图12为不同弹道倾角时子弹药散布面积Sd随抛撒高度的变化曲线。

2.2.2 子弹药散布面积分析

他又向后翻看，后面皆是蝇头小字，密密麻麻地记载了九州各地的风俗物产、舆地历史等信息。他下意识地找到巴州的那部分，果然，里边有关于唐门的介绍。

由图11可知,在母弹开舱离心抛撒子弹药时,弹道倾角越小、抛撒高度越高,子弹药抛撒中心射程越高。

对人类社会的发展规律缺乏深刻、正确的认识和把握。主要表现为：一是“三观”的迷茫和错误。即世界观、人生观和价值观的迷茫和错误。二是“三个自信”的缺乏。即不认同中国特色社会主义道路、理论、制度。三是“三个主义”的侵蚀。即受拜金主义、享乐主义、极端个人主义的腐蚀。四是“三个意识”的淡薄。即法律意识、纪律意识和道德意识的不强，导致政治底线、法律底线、道德底线“三个底线”失守。

图12 v=300 m/s,n=2 000 r/min,不同弹道倾角时抛撒高度与子弹药散布面积关系曲线

图13为不同母弹转速时子弹药散布面积随抛撒高度变化曲线。

根据以上数据可以看出，我国个人理财市场扩大的速度非常惊人。事实上，在西方发达国家，理财产品一直是一种重要金融工具，而我国的相关业务起步较晚[2]，但是增长速度惊人，这一方面体现了我国居民理财意识的增强，另一方面也体现了个人理财产品容易被更多的个人投资者所接受。

图13 v=300 m/s,θ=-30°,不同母弹转速时子弹药散布面积与抛撒高度关系曲线

由图12和图13可知,在母弹开舱离心抛撒子弹药时,母弹转速越高、弹道倾角越小、抛撒高度越高,子弹药的散布面积越大。

2.2.3 子弹药散布均匀性分析

图14为单位区域内不同弹道倾角子弹药散布落点图。

图14 n=2 000 r/min,h=500 m,v=300 m/s,不同弹道倾角子弹药散布落点

由图14可知,在母弹开舱离心抛撒子弹药时,弹道倾角越小,子弹药在单位区域内的散布越均匀;弹道倾角越大,散布边缘子弹药数越多,中心位置的子弹药数越少,子弹药散布落点的重复性越高。子弹药散布落点的图形均匀性受母弹弹道倾角的影响较大。

3 结论

本文主要研究结论如下:

①尾舵偏转产生的滚转力矩是影响转速的主要因素,在单位时间内,当来流条件和姿态角保持不变时,初始滚转舵偏所产生的滚转力矩越大,全弹达到的转速越高;在保持全弹纵向操纵性能的前提下,增加尾舵展长并减少弹体表面除尾舵外的其他气动阻尼面有利于加快全弹的主动起旋过程,提高转速。

②在离心抛撒子弹过程中,母弹的转速越高、开舱高度越高、弹道倾角越小,则子弹药的散布面积越大、散布图形越均匀;母弹的末端飞行速度越大,抛撒子弹药的中心射程越大;同时,通过调整全弹初始舵偏角和开舱条件,可实现子弹药散布面积与散布密度的可控调节。

③使用主动气动起旋离心抛撒技术可减少弹上抛撒机构设备,有效增加子弹药的装填数量,在提高作战使用效率的同时,降低了成本,提高了子弹药抛撒稳定性。该项技术的应用可为不同类型航空子母弹药的抛撒方案设计提供参考和借鉴。

参考文献

[1] 谢汉桥,刘述,陈振教. 大展弦比轴对称气动布局应用研究[J]. 弹道学报,2016,28(4):42-46.

XIE Hanqiao,LIU Shu,CHEN Zhenjiao. Study on application of high-aspect-ratio axial-symmetrical aerodynamic configuration[J]. Journal of Ballistics,2016,28(4):42-46. (in Chinese)

[2] 王晓冰,赵忠良. 窄条翼导弹俯仰机动中滚转失稳及其控制过程[J]. 航空学报,2016,37(8):2 517-2 524.

WANG Xiaobing,ZHAO Zhongliang. Roll instability and control during pitching maneuver for a missile with strake wings[J]. Acta Aeronautica et Astronautica Sinica,2016,37(8):2 517-2 524. (in Chinese)

[3] DA Xingya,TAO Yang,ZHAO Zhongliang. Virtual flight Navier-Stokes solver and its application[J]. Procedia Engineering,2012,31:75-79.

[4] 童继进,刘忠,毛超. 基于高速旋转飞行器的箔条抛撒运动特性[J]. 海军工程大学学报,2013,25(5):83-87.

TONG Jijin,LIU Zhong,MAO Chao. On kinematics characteristic of chaff release based on aerocraft flighting and revolving around an axis at high speed[J]. Journal of Naval University of Engineering,2013,25(5):83-87. (in Chinese)

[5] KHALIL M,ABDALLA H,KAMAL O. Dispersion analysis for spinning artillery projectile[C]//The 13th International Conference on Aerospace Sciences and Aviation Technology. Cairo,Egypt:Military Technical College,2009:21-29.

[6] 阎超. 计算流体力学方法及应用[M]. 北京:北京航空航天大学出版社,2006.

YAN Chao. Computational fluid dynamics method and its application[M]. Beijing:Beihang University Press,2006. (in Chinese)

[7] 韩子鹏. 弹箭外弹道学[M]. 北京:北京理工大学出版社,2008.

HAN Zipeng. Rocket exterior ballistics[M]. Beijing:Beijing Institute of Technology Press,2008. (in Chinese)

[8] 陶如意,吴艳滨,王浩. 时序抛撒子母弹多体干扰气动特性的数值模拟[J]. 弹道学报,2011,23(3):53-57.

TAO Ruyi,WU Yanbin,WANG Hao. Numerical simulation of multi-body interference aerodynamic performance of cluster munition dispensed by sequential[J]. Journal of Ballistics,2011,23(3):53-57. (in Chinese)

[9] 李东光,殷跃强,杨瑞伟,等. 高速旋转子母弹开舱高度自适应技术研究[J]. 北京理工大学学报,2015,35(2):167-169.

LI Dongguang,YIN Yueqiang,YANG Ruiwei. Research on self-adaption dispersal altitude selection of cluster munition with high rotation[J]. Transaction of Beijing Institute of Technology,2015,35(2):167-169. (in Chinese)

[10] 谭大成,袁曾凤. 母弹抛撒技术及其内弹道计算[J]. 北京理工大学学报,2000,20(3):273-277.

TAN Dacheng,YUAN Cengfeng. Cluster warhead dispering technique and its interior ballistic calculation[J]. Transaction of Beijing Insitute of Technology,2000,20(3):273-277. (in Chinese)

[11] 杨瑞伟,殷跃强,李东光. 基于母弹弹道诸元误差值的自适应开舱方法[J]. 兵工学报,2014,35(6):801-807.

YANG Reiwei,YIN Yueqiang,LI Dongguang. Self-adaptive cabin opening method based on firing data error of shrapnel[J]. Acta Armamentarii,2014,35(6):801-807. (in Chinese)