提拉缸式发射装置水下发射内弹道仿真研究

刁宏伟，李宗吉，李春来

(中国人民解放军91829部队, 辽宁 大连 116041； 2.海军工程大学 兵器工程系, 武汉 430033)

随着UUV技术的发展,未来海战中UUV[1]配备鱼雷武器[2]用于水下无人作战是必然趋势。但由于UUV一般排水量较小，传统潜艇带发射动力的鱼雷发射装置[3][4](如气动不平衡或水压平衡发射装置)因体积重量过大，无法直接应用在该发射平台上，目前的UUV一般采用自航发射方式发射鱼雷[5][6][7]。采用自航发射方式虽结构简单，但在发射鱼雷的过程中UUV自身的运动状态可能会因受到鱼雷自航运动的干扰而发生变化，影响平台的操作与控制。为此需要探索结构简单，适用于UUV这类小型平台的有发射动力的装置。提拉缸式发射装置是在参考借鉴防空导弹发射技术[8][9]基础上提出的，为了论证该方案在UUV上的可行性，本文综合运用MATLAB与FLUENT对UUV运用提拉缸发射装置发射鱼雷的过程进行了仿真研究。

REN Qian, ZENG Shu-xiong, TANG Shi-jie, DAI Li-he, XU Jin-shan, SUN Ying-hao, XU Chuan-liang

1 发射装置原理

本文提出一种以火药燃烧的化学能作为发射能源的鱼雷发射装置(见图1)。火药气缸长度约为发射管的一半，拉杆总行程为发射管的一半。发射前，发射管密封储存鱼雷，发射时发射管前后盖与发射管的连接由爆炸螺栓切断并与发射管脱离，此时发射管前后端连通海水。然后1号火药柱点火产生火药气体推动拉杆带动鱼雷向管口方向运动。当活塞运动1m时活塞越过2号火药柱并引燃2号火药柱，此时1、2号火药柱共同燃烧推动活塞运动(见图2)。活塞运动3m时拉杆与鱼雷脱离，鱼雷克服海水阻力继续运动直至出管。

2 仿真方法

鱼雷发射装置内弹道仿真通常采用数值求解内弹道微分方程组模型。这种方法较精细的模拟发射过程中发射管内部流场及鱼雷所受的流体阻力的变化情况非常困难，对内弹道仿真结果的误差较大。为此本文综合采用MATLAB与FLUENT[10]进行仿真[11]。在采用MATLAB仿真程序(程序数学模型见3.1小节)对发射过程中气缸中各物理量进行仿真的同时，运用FLUENT对鱼雷发射管内流体流动与鱼雷所受阻力进行仿真。MATLAB与FLUENT采用共享数据文件的方式进行数据交换，通过更新标志位的方法使MATLAB与FLUENT实时同步仿真。整个仿真流程见图3。

3 模型建立

3.1 数学模型

仿真数学模型由15个微分方程式组成，代入4阶定步长Runge-Kutta方法求解，其中鱼雷运动流体阻力由FLUENT模拟发射管流场计算，鱼雷速度、加速度由MATLAB解算，并传递给FLUENT。

鲁西的含氨基酸液体肥含多种游离态氨基酸，促进根系分化，加速根系发育，培育壮苗，促进插枝生根，提高存活率。而液钙、氨基钙叶面肥富含活性钙，防止裂果效果明显，在蔬菜和水果的结果期十分有效。提高作物的免疫力，抗外界不利影响能力，利用率100%。

(1)

式中：Uc为发射管内气体内能；Qe为火药燃烧热；Lt为鱼雷推动功。

(2)

式中：Tc为气缸内温度；Cv为气缸内气体定容比热容；Mc为气缸内气体总质量。

(3)

中国食品工业总产值从百亿元级到万亿元级的背后，是众多功勋人物、领军人物孜孜以求、攻坚克难的不懈努力，是众多食品企业争创百年老店、力求让消费者满意的永恒追求。

dMc=dMe1+dMe2

(4)

式中：Me1、Me2为两段火药的燃烧质量。

(5)

式中：Sc为气缸断面面积；Lp为活塞行程；Lc为气缸总长度。

式中：Sy为火药表面积；S1为火药初始表面积；β为药型系数；ψ为火药燃烧百分比数。

dQe=EqQb(dMe1+dMe2)

(6)

式中：Qe为火药燃烧总热量；Eq为热损失系数通常取0.8；Qb为火药燃烧热。

(7)

(8)

式中：U10为火药在单位压力下的燃烧速度；ε为火药燃烧时的通气系数。

(9)

如图5所示，当故障发生在区段Z6，FI1、FI2、FI3、FI5均正常上报故障信息，故障指示器FI4漏报故障信息，其他故障指示器也无故障信息上报。仍设每个区段发生故障的概率一样，λ取0.95，δ取0.1。

(10)

式中：Sm为火药端面面积；Sm0为火药初始端面面积。

(11)

(12)

式中：Fp为通过FLUENT仿真并拟合的气缸阻力(活塞后端联通海水，孔径0.01 m)，与活塞速度Vp有关。

dVp=(PcSc-Fflu-Fp)/(mt+mp)

(13)

式中：Sc为气缸端面面积；Fflu为通过FLUENT仿真传递的鱼雷阻力；mt、mp为鱼雷质量与活塞质量。

运行FLUENT及MATLAB软件，读取仿真模型与仿真初始值，进行发射水下发射鱼雷过程仿真。

dlt=Vp

(14)

式中：lt为鱼雷行程。

dLt=ScPcVp

(15)

3.2 发射管CFD模型

仿真中的主要参数为：鱼雷质量2 000 kg、气缸直径0.125 m、气缸长度3.5 m、提拉缸活塞质量20 kg、发射火药质量0.09 kg(两段)。

在CFD模型中立动网格，以鱼雷作为刚体动网格边界，沿X轴运动。由于本文建模采用非结构网格，网格更新采用弹簧光顺模型(smoothing)以及局部网格重划模型(remeshing)两种。

CFD模型网格共有节点28 798个，最小网格面积5.6×10-3m2，最大网格面积3.8×10-2m2，最大畸变率为0.64。发射管长6.5 m、直径0.537 m，鱼雷长6 m、直径0.534 m。发射管外流场直径为3 m，长度为7 m。前段边界条件为压力出口(pressure-outlet)，尾端为压力入口(pressure-inlet)。

牛哞悠长，挽篮入园，随心所欲采摘时鲜。拾级而上，山坡有亲子厨房。厂房车间式开放厨房，石垒砖砌灶台，铁锅劈柴烟囱，炊具一应俱全。可一试身手卖弄厨艺，亦可请教村夫农妇烹饪窍门。可见反客为主非诳语，此处可学故人邀诗人，把盏话桑麻。

4 发射过程内弹道仿真

记忆中，在我4岁那年，由中国国际电视总公司出品了一部41集古装神话剧——改编自明代小说家吴承恩同名文学古典名著《西游记》。1986年春节一经播出便轰动全国，可谓老少皆宜，获得了极高评价。至今仍是寒暑假期间被重播最多的电视剧之一，百看不厌，成为一部公认的、无法超越的经典。该剧讲述的是孙悟空、猪八戒、沙僧辅保大唐高僧玄奘(唐僧)去西天取经的故事，师徒四人一路抢滩涉险，降妖伏怪，历经九九八十一难，取回真经，终修成正果的故事。

在仿真中，0～0.3 s时1号火药柱燃烧，0.3 s以后1、2号火药柱共同燃烧；0.5 s时拉杆位移达到3 m并与鱼雷脱离。发射过程的仿真结果如图5～图9所示。

根据发射方案建立发射管CFD仿真模型。由于发射管及流体域轴对称，本文采用二维网格轴对称方式建立流场网格以减小计算量(见图4)。

仿真开始后，CFD软件将流场计算获得的作用在鱼雷上流体阻力，通过共享文件将数据传递给依据3.1数学模型编制的MATLAB仿真程序，仿真程序通过解微分方程组，再将计算出的鱼雷位移参数传递给共享文件CFD软件，CFD根据位移参数进行重新进行网格划分，重新进行新一时刻的流场计算，如此循环直至结束，其仿真流程见图3。

式中：Pc为气缸内气体压力；k为火药燃气比热比；Vc为气缸内气体体积。

以关系到千家万户的粮食加工行业为例，在解决粮食供应问题后，人们的生活开始由温饱型向小康型逐步过渡。如小型砻谷机，碾米机和200型榨油机、“东方红”牌粮机和“双狮”牌砻谷胶辊等具有代表性的机械，开始被运用到米、面、油的深加工、精加工上来，以提高原料的利用率。

从仿真结果来看，1号火药柱点火后发射温度、压力迅速上升至1 500 K、16 MPa，活塞运动后火药气体对外做功，温度、压力下降至500 K、10 MPa时(见图6、图7)，活塞越过2号火药柱并将其引燃，使温度下降梯度减小。0.5 s后活塞运动至气缸尾，火药气体不再对外做功，火药柱继续燃烧，温度回升。

从鱼雷速度曲线来看，速度在0.5 s时达到峰值12 m/s，然后拉杆脱离鱼雷，在海水阻力作用下鱼雷速度降至10 m/s,整个发射过程大约耗时0.6 s(见图8),基本满足水下发射鱼雷的要求。

漏电保护器按极数分类，有单极2线 双极2线3极3线 3极4线和4极4线等多种形式，其在低压配电线路中的接线如图1所示：RCD1—单极2线，RCD2—双极2线，RCD3—3极3线，RCD4—3极4线，RCD5—4极4线，QF—断路器，YR—漏电脱扣器.

通过进行发射过程流场仿真，不仅获得发射过程的发射管内部及外部流场压力变化情况(见图9～图12)，还可获得发射过程鱼雷所受流体阻力变化情况(见图13)。从鱼雷阻力曲线可以看到鱼雷阻力在0.5 s以及0.65 s时存在两次跳变。在0.5 s时鱼雷与拉杆脱离，鱼雷加速度由正值变为负值，鱼雷附加阻力降低为0；在0.65 s时，鱼雷尾锥部位经过发射管口，鱼雷与发射管抽吸作用减小，鱼雷压差阻力减小。

5 结论

为了对提拉缸式鱼雷发射装置在UUV上应用可行性进行分析与研究，综合采用MATLAB与FLUENT对发射过程进行仿真。结果表明，鱼雷出管速度达到10 m/s，满足水下发射的技术要求。与其他潜艇鱼雷发射装置相比，这种方案结构简单、体积小、重量轻，可作为UUV一次性使用发射装置，发射后可整体抛弃，适用于UUV之类的小型水下平台。本文所建的模型可为后续的发射装置总体方案论证及技术设计提供参考。

：

[1] 李德远,吴汪洋,李晓晨.军用UUV的发展与应用前景展望[J],舰船电子工程，2014,32(4):22-24.

[2] 周杰,王树宗.无人水下航行器发射鱼雷初始弹道分析[J].弹道学报,2009(3):86-89.

[3] 练永庆，王树宗.鱼雷发射装置设计原理[M].北京：国防工业出版社 2012:125-134.

[4] 朱清浩,宋汝刚.美国潜艇鱼雷发射装置使用方式初探[J].鱼雷技术,2012,20(3):215-219.

[5] 孙玉松,李宗吉,练永庆.UUV自航发射鱼雷过程仿真[J].鱼雷技术，2014,22(5):396-400.

[6] 章杰，张云海，赵加鹏.空投UUV入水冲击仿真研究[J].兵器装备工程学报，2016(1):70-73.

[7] 李亚男,叶慧娟,张西勇.水压平衡式鱼雷发射装置发射小型运载器内弹道仿真[J].兵器装备工程学报，2016(2):41-44.

[8] 杨风波,马大为,杨帆.高压弹射装置内弹道建模与计算[J].兵工学报，2013,34(5):527-534.

[9] 白鹏英，乔军.双级气缸式弹射装置内弹道分析[J].现代防御技术，2007(8)：44-49.

[10]蔡荣泉.船舶计算流体力学的发展与应用[J].船舶,2002,8(4):8-13.

[11]谢海斌，张代兵，沈林成.基于MATLAB/SIMULINK与FLUENT的协同仿真方法研究[J].系统仿真学报，2007，19(8):1824-1827.