舰船科学技术  2020, Vol. 42 Issue (11): 35-39    DOI: 10.3404/j.issn.1672-7649.2020.11.007   PDF    
舰载导弹在水下爆炸作用下的冲击环境数值仿真
谷鸿平1, 黄河2, 孟献策3, 张立建1, 畅博1, 吕永柱1, 张舵2     
1. 西安近代化学研究所,陕西 西安 710065;
2. 国防科技大学 文理学院,湖南 长沙 410073;
3. 火箭军装备部装备项目管理中心,北京 100085
摘要: 采用Ansys/Ls-dyna动力分析软件,针对水下爆炸对大型航空母舰的冲击响应过程及对其内舰载导弹的冲击环境问题进行数值模拟。引入冲击环境谱分析研究方法,建立航母的1∶1尺寸有限元模型以适用于冲击动力响应分析;研究水下爆炸对航母结构的冲击响应过程,分析航母导弹贮存典型位置处的冲击环境并给出了加速度谱。结果表明:水下爆炸会对航母结构造成较大冲击,舰体中间位置底部冲击响应最明显,航母整体结构表现出鞭状运动响应的特征;典型舰载导弹贮存位置处的冲击谱值较大,为了导弹的贮存安全,导弹装载平台应设置缓冲减振装置;冲击强度从底板往上开始衰减迅速,而后衰减渐趋于平缓。
关键词: 舰载导弹     水下爆炸     数值模拟     冲击环境     冲击谱    
Numerical simulation on dynamic response of aircraft carrier subjected to underwater blast landing
GU Hong-ping1, HUANG He2, MENG Xian-ce3, ZHANG Li-jian1, CHANG Bo1, LU Yong-zhu1, ZHANG Duo2     
1. Xi′an Modern Chemistry Research Institute, Xi′an 710065, China;
2. College of Liberal Arts and Sciences, National University of Defense Technology, Changsha 410073, China;
3. Project Management Center of Rocket Army Equipment Department, Beijing 100085, China
Abstract: Numerical simulation of the dynamic response process and the impact environment of shipborne missile of big aircraft carrier subjected to underwater blast landing are investigated with Ansys/Ls-dyna. The shock spectrum method is introduced and the finite element model of an aircraft carrier with real dimension is built for the analysis of dynamic response. The impact response of underwater explosion on aircraft carrier structure is studied. The impact environment of typical carrier missile storage location is analyzed and the acceleration spectrums are given. The results show that the underwater explosion will cause great impact on the aircraft carrier structure, the impact response of the bottom of the middle position is the most obvious, and the overall structure shows the characteristics of whip motion response. The shock spectrum value at the typical carrier missile storage location is relatively large. For the sake of the missile storage safety, the missile loading platform should be equipped with a shock absorber. The impact strength decay rapidly from the bottom at first and then gradually flatten out.
Key words: shipborne missile     underwater explosion     numerical simulation     impact environment     shock spectrum    
0 引 言

在海上作战中,作为“舰队核心”的航空母舰毫无疑问地会成为炸弹、导弹、鱼雷及水雷等武器的集群火力打击目标,舰载导弹及舰载设备等因此会面临更多的冲击考验,严重时会导致导弹失效甚至舰内引爆[1-3]。水下爆炸是一个极大能量瞬时转换的过程,反应场所位于水下极小面积上或体积内,通常划分为接触爆炸与非接触爆炸:接触爆炸会引起舰体破损,设备也会受到波及,但波及范围相对而言并不大;近距离非接触爆炸作用会导致舰体剧烈的振动及大程度的纵弯曲运动,使舰船受到中拱、中垂损伤甚至断裂沉底;非接触爆炸在中远距离发生时,气泡在多次脉动的过程中会诱导船体做低阶耦合运动,出现鞭状响应致使船身和众多设备运动变形,各类核心舰载装置因而受到程度严重且范围广泛的损伤,此为当下舰艇失效作战中一种关键的策略,备受各国海上军事力量瞩目[4-7]。因此,开展非接触水下爆炸作用下航空母舰舰载导弹的冲击环境研究,具有深远的现实影响及重大的军事应用价值。

本文基于Ansys/Ls-dyna软件数值模拟研究典型航母内导弹贮存平台位置处的冲击环境特性,并分析冲击强度在航母结构内的衰减过程,为提升舰载不敏感弹药研制及大型舰船生命力提供参考。

1 冲击环境研究方法

舰船设备的冲击环境是指在水下爆炸作用下舰船设备的基础输入,一般体现为设备响应的冲击速度、加速度及位移[12]。舰船设备的冲击环境是其抗冲击设计的基本依据。

在冲击环境的研究刚起步时,设备被简化等效为仅含一个独立变量的单自由度系统,以Duhamel积分为工具,对等效质量块的基础输入进行求解。把求解得到的位移称为冲击位移,等效质块在单个重力加速度载荷作用下的位移称为自重位移,冲击设计数通过前者比后者求出。这种描述方法对于具备多阶模态的设备遭遇不同类型攻击的情况显然不能较好适用。构建一种普适性更佳且能更直观描述冲击环境的模型迫在眉睫,“冲击响应谱”应运而生。把一系列固有频率不等的单位质量线性振子安装在同一基础上,在同一基础运动输入条件下,不同振子遭受大小相同的冲击加速度与冲击速度,可以得到振子响应的最大幅值与其固有频率的关系,用曲线表示该关系即可得到冲击响应谱[13]。其卓越之处在于定量分析,能得到设备响应冲击的绝对加速度和以基础为参考的相对速度与位移。

假设基础有一瞬态冲击运动,在此基础上设置众多频率不等的单自由度振子,(见图1),在同一基础运动输入下各个振子的响应峰值即为冲击谱的谱值。谱值由振子的绝对加速度以及相对于基础的位移和速度三者组成。

图 1 冲击谱的定义图解 Fig. 1 Definition diagram of shock spectra

在基础运动是zt)的条件下,振子i的绝对位移记作yi,相对位移记作xi,振子的绝对运动方程表示为:

${\ddot y_i} + \omega _i^2(y - z) = 0,$ (1)

相对运动方程表示为:

${\ddot x_i} + \omega _i^2x = - \ddot z(t),$ (2)

假设初始速度和位移均为0,则xi的解可得:

${x_i}(t) = - \frac{1}{{{\omega _i}}}\int\nolimits_0^1 {\ddot z\left( \tau \right)} \sin {\omega _i}\left( {t - \tau } \right){\rm d}\tau \text{。}$ (3)

振子i的圆频率记为ωi,其位移谱值即为ωi(t)绝对值的最大值。对式(3)求微分化为:

${\dot x_i}(t) = - \int\nolimits_0^1 {\ddot z\left( \tau \right)} \sin {\omega _i}\left( {t - \tau } \right){\rm d}\tau \text{。}$ (4)

振子的速度谱值即为振子相对于基础的速度绝对值的最大值。从式(1)发现振子的绝对加速度与振子相对于基础的位移二者有一定的比例关系,比例系数为圆频率的二次方。因此速度、加速度谱值基于位移谱值而定,它们之间存在的关系为:

${V_i} = \omega {D_i}\;{A_i} = \omega _i^2{D_i}\text{。}$ (5)

其中:Di为位移谱;Vi为速度谱;Ai为加速度谱。连接所有振子的谱值得到的曲线即为所求的冲击谱曲线。

2 计算模型

选取美国海军“尼米兹”级大型航空母舰为研究对象,该舰长332.8 m,宽40.8 m,高76.0 m,飞行甲板最宽处76.8 m,最大吃水11.3 m。“尼米兹”级整体采用内外层结构设计,中间以X型钢梁连接,在弹药库、发动机舱等重要位置加装了150 mm以上的凯夫拉和钢制装甲。舰体由上至下依次设有7层甲板,整个舰体被33道水密和防火隔舱划分为2000多个水密舱室。基于以上数据,利用Ansys与Hypermesh软件进行简化的有限元建模与网格划分。在建模过程中,合理选取质量、薄壳、三维梁单元,具体为MASS166,SHELL163,BEAM161。考虑到计算规模及效率的优化,计算水域选取水面至距航母底面30 m处。在Ansys建模中,航母没于水面下的部分通过布尔运算挖空,水域采用SOLID164 8节点实体单元进行网格划分。舰体模型基于Johnson-Cook材料模型[14],遵循Gruneisen状态方程。参数的选取可以参见文献[15-16]。建立的有限元计算模型整体效果如图2所示。

图 2 有限元模型 Fig. 2 Finite element model

目前大多数反舰导弹及鱼雷等的TNT当量都在100~500 kg之间,能对舰体造成较大影响的爆距大约在50 m以内。基于此,本文研究300 kg TNT当量炸弹在航空母舰底面下20 m处爆炸时对航母的冲击。监测点位置如图3所示。

图 3 导弹贮存平台监测点位置示意图 Fig. 3 Map of monitoring point location of missile storage platform
3 冲击响应过程分析

TNT当量为300 kg的炸弹在距航空母舰舰底20 m距离处爆炸时,航母的响应过程如图4所示(位移云图的单位为m)。当冲击波传至舰底,舰底中心地带的外板及龙骨首当其冲并产生竖直向上的位移,如图4(a)所示。响应迅速向舰顶及舰两端传播,如图4(b)所示,舰体中部产生的竖向位移继续增大,同时向舰体两端传播;当舰体中部位移达到最大,如图4(c)所示,位移响应开始下降,而两端的位移开始增大,如图4(d)所示,类似于鞭状响应[9]。随后,舰体经历几次鞭状响应过程,并趋于稳定。

图 4 水下爆炸作用下航母的响应过程位移云图 Fig. 4 Response displacement nephogram of aircraft carrier under underwater explosio
4 导弹贮存平台冲击环境分析 4.1 冲击响应谱分析

分析图3的导弹贮存平台的典型位置,这些测点的加速度时程曲线及相应的加速度冲击响应谱如图5所示。测点B位于舰中部甲板中心,距爆点最近,所受的冲击作用最为严重(见表1),其加速度最大值达到129 g,冲击响应谱中与951 Hz频谱相对应的加速度谱值为578.6 g。相比之下AC点的加速度谱最大值明显小得多,分别为262.7 g和216.2 g,如图5表1所示。

表 1 典型位置冲击环境结果统计表 Tab.1 Statistical table of typical location impact environment results

图 5 监测点加速度时程曲线与加速度冲击响应谱 Fig. 5 Acceleration time history curve and acceleration impact response spectrum of monitoring points

现以冲击环境标准GJB150中的轻、中型设备为例,轻型设备(m<200 kg)在3个方向上的冲击环境存在的差异较小,谱加速度的最大值为2 000 g。中型冲击机用于考核中重量设备(120 kg<m<2700 kg)的冲击试验,其谱加速度的最大值为200 g。冲击环境在垂向上是在横向上的2倍[17]。因而在本文研究爆炸载荷(300 kgTNT当量舰底下20 m爆炸)作用下,轻型设备的冲击环境能够满足所需,而中重型设备的冲击环境难以保证设备的贮存安全,需要增设缓冲减振装置。

4.2 冲击强度衰减分析

为了研究航母底面水下爆炸冲击在航母内不同甲板层上产生响应的强烈程度,以工况1为例进行分析。对底板中心正上方各层甲板上的点进行监测,底板中心点记为DB,各层甲板上的点从主甲板往下分别记为D1D7,如图6所示。

图 6 各层甲板上的监测点示意 Fig. 6 Indication of monitoring points on each deck

基于计算结果易得位于底板中心点处的最大冲击谱值达2741 g。分析上述水下爆炸工况冲击强度的衰减规律,基于仿真求解各个位置加速度谱的最大值,并对其规律进行拟合,如图7所示。加速度谱的最大值记作Amax,横轴记为X,各个监测点记为DDB对应x=1,D7对应x=2,……),拟合各点处的加速度谱值,可得:

图 7 冲击强度衰减分析图 Fig. 7 Analysis chart of impact strength attenuation
${A_{\max }} = 658.65 + 3\;267.6{{\rm{e}}^{ - 0.42X}}\text{。}$ (9)

图7可知,冲击强度从底板向上到第4层甲板衰减速度较快,衰减了约64%,从第4层甲板向上衰减速度逐渐趋于平缓,到主甲板时与初始冲击强度相比衰减了约71%。

5 结 语

1)用于动力分析计算而建立的航母1∶1有限元模型和水域的计算模型的适用结果较好;

2)水下爆炸会给航母舰体结构带来较大冲击,诱导鞭状响应的出现,导致舰体及大范围设备纵向变形,舰体中间位置底部发生的冲击响应尤为明显;

3)典型航母舰舰载导弹贮存位置处的冲击谱值较高,应设置缓冲减振装置以提高弹药的稳定性贮存水平;

4)从航母底板向上冲击强度的衰减速度起初时较大,随后逐渐趋向稳定。

参考文献
[1]
姚熊亮, 侯健, 王玉红, 等. 水下爆炸冲击载荷作用时船舶冲击环境仿真[J]. 中国造船, 2003(44): 71-74.
YAO Xiong-liang, HOU Jian, WANG Yu-hong, et al. Research on simulation of underwater shock environment of ship[J]. Ship Building of China, 2003(44): 71-74.
[2]
王国治, 支李峰, 邓文彬. 水下爆炸环境中舰船浮筏装置冲击响应研究[J]. 江苏科技大学学报(自然科学版), 2006(5): 4-9.
[3]
吴庭翱, 张弩, 侯海量, 等. 水下接触爆炸下多舱防护结构载荷特性及动响应研究进展[J]. 中国舰船研究, 2018, 13(3): 1-12.
WU Ting-ao, ZANG Nu, HOU Hai-liang, et al. Research progress on load characteristics and dynamic response of multicamerate defense structure subjected to contact under water explosion[J]. Chinese Journal of Ship Research, 2018, 13(3): 1-12.
[4]
李磊, 冯顺山. 水下爆炸对舰船结构毁伤效应的研究现状及展望[J]. 舰船科学技术, 2008, 30(8): 26-30.
LI Lei, FENG Shun-shan. Present state and perspectives of damage effect to ship construction subjected to UNDEX[J]. Ship Science and Technology, 2008, 30(8): 26-30.
[5]
陈辉, 李玉节, 潘建强, 等. 水下爆炸条件下不同装药对水面舰船冲击环境的影响试验研究[J]. 振动与冲击, 2011, 30(7): 16-20.
[6]
史晓骏. 水下爆炸冲击波载荷作用下舰船板架结构冲击可靠性研究[J]. 舰船科学技术, 2017(6): 14-16.
SHI Xiao-jun. Study of ship plate frame structure reliability during underwater explosion shock[J]. Ship Science and Technology, 2017(6): 14-16.
[7]
李海涛, 张振华, 牟金磊, 等. 水下爆炸作用下弹塑性船体梁整体运动模型及损伤特性[J]. 工程力学, 2019, 36(1): 239-242.
[8]
库尔. 水下爆炸[M]. 罗耀杰, 等译. 北京: 国防工业出版社, 1960.
[9]
董海, 刘建湖, 吴有生. 水下爆炸气泡脉动作用下细长加筋圆柱壳的鞭状响应分[J]. 船舶力学, 2007, 11(2): 250-258.
DONG Hai, LIU Jian-hu, WU You-sheng. Whipping response analysis of slender stiffened cylindrical shell subjected to underwater explosion with bubble pulse[J]. Journal of Ship Mechanics, 2007, 11(2): 250-258.
[10]
ZAMYSHLYAYEV BV, YAKOVLEV YS. Dynamic loads in underwater explosion[R]. AD0757183, 1973. Chinese Journal of Ship Research, 2016, 11(1): 46-56, 71.
[11]
刘建湖, 周心桃, 潘建强, 等. 舰艇抗爆抗冲击技术现 状和发展途径[J]. 中国舰船研究, 2016, 1l(1): 46-56.
LIU Jian-hu, ZhOU Xin-tao, PAN Jian-qiang, et al. The state analysis and technical development routes for the anti-explosion and shock technology of naval ships[J]. Chinese Journal of Ship Research, 2016, 1l(1): 46-56.
[12]
王珂, 尹群, 嵇春艳, 等. 典型舰船水下爆炸冲击环境仿真[J]. 计算机辅助工程, 2006(15): 95-98.
[13]
姜涛, 由文立, 李建林. 动态设计分析方法的舰船甲板设备垂向冲击环境分析[J]. 兵工学报, 2010, 31(1): 33-36.
[14]
Johnson GR, Cook WH. Fracture Characteristics of Three Metals Subjected to Various Strains, Strain rates, Temperatures and Pressures[J]. Eng Fract Mech, 1985, 21(1): 31-48. DOI:10.1016/0013-7944(85)90052-9
[15]
姚术健, 蒋志刚, 卢芳云, 等. 汽车炸弹钢箱梁内部爆炸局部破坏效应分析[J]. 振动与冲击, 2015, 34(7): 222-227.
YAO Shu-jian, JIANG Zhi-gang, LU Fang-yun, et al. Analysis on local damage of steel box girder under internal blast loading of vehicle bomb[J]. Journal of Vibration and Shock, 2015, 34(7): 222-227.
[16]
Yao S J, Zhang D, Lu F Y, et al. Experimental and numerical studies on the failure modes of steel cabin structure subjected to internal blast loading[J]. Int J Impact Eng, 2017, 110: 279-287. DOI:10.1016/j.ijimpeng.2017.03.006
[17]
刘建湖, 潘建强, 何斌. 各主要国家海军舰艇抗冲击标准之比较[C]//. 第三届全国舰艇抗冲击技术交流会(哈尔滨), 2009: 146-154.