截锥形弹体在液体介质中运动速度衰减规律分析

孔祥韶; 石干; 王旭阳; 周红昌; 吴卫国

doi:10.11883/bzycj-2020-0075

截锥形弹体在液体介质中运动速度衰减规律分析

doi: 10.11883/bzycj-2020-0075

孔祥韶^{1, 2,},
石干^{1, 2},
王旭阳^{1, 2, 3},
周红昌^{1, 2},
吴卫国^{1, 2}

1.
武汉理工大学绿色智能江海直达船舶与邮轮游艇研究中心，湖北武汉 430063
2.
武汉理工大学交通学院，湖北武汉 430063
3.
西安航天动力研究所液体火箭发动机技术重点实验室，陕西西安 710100

基金项目: 装备预研教育部联合基金（青年人才）（6141A02033108）

详细信息

作者简介:
孔祥韶（1983－　），男，博士，教授，博士生导师，kongxs@whut.edu.cn

中图分类号: O353.4
计量
- 文章访问数: 637
- HTML全文浏览量: 250
- PDF下载量: 78
- 被引次数: 0
出版历程
- 收稿日期: 2020-03-20
- 修回日期: 2020-06-24
- 刊出日期: 2021-01-05

On velocity attenuation of a truncated cone-shaped projectile vertically penetrating through liquid

KONG Xiangshao^{1, 2
,},
SHI Gan^{1, 2},
WANG Xuyang^{1, 2, 3},
ZHOU Hongchang^{1, 2},
WU Weiguo^{1, 2}

1.
Green & Smart River-Sea-Going Ship, Cruise and Yacht Research Center, Wuhan University of Technology, Wuhan 430063, Hubei, China
2.
School of Transportation, Wuhan University of Technology, Wuhan 430063, Hubei, China
3.
Science and Technology on Liquid Rocket Engine Laboratory, Xi’an Aerospace Propulsion Institute, Xi’an 710100, Shaanxi, China

摘要

摘要: 水面舰船被动防护体系中液舱的主要功能之一是阻止高速弹体（爆炸破片）对内部重要结构、设备和人员的威胁，高速弹体打击液舱的过程包含着复杂的能量传递与耗散。为了分析弹体形状对其在液体介质中运动速度衰减的影响, 开展了一系列不同头形因数的截锥形弹体在不同入水速度下弹体垂直侵彻液体介质过程的数值模拟，得到了垂直侵彻液体介质时弹体速度衰减特性，发现高速弹体在液体介质中运动的阻力因数与弹体形状和无量纲速度有关。基于对系列数值模拟计算结果的拟合分析，提出了计及头形因数的截锥形弹体垂直侵彻液体介质时的速度衰减经验公式，通过开展数值算例分析验证了公式计算结果的可靠性。本文中提出的经验公式可实现对高速弹体在液体介质中速度衰减的准确快速计算，为舰船防护液舱结构设计提供一定的参考。
- 截锥形弹体 /
- 头形因数 /
- 速度衰减 /
- 防护液舱
Abstract: One of the main functions of liquid tanks in the passive protective systems of surface warships is to prevent the damage by high-velocity projectiles (explosive fragments) to important internal structures, equipments and personnels. The process of high-velocity projectile penetrating through a liquid tank involves complex energy transfer and dissipation. To explore the influences of the head shape of a projectile and its water-entry velocity on the velocity attenuation of the projectile in fluid, a series of truncated cone-shaped projectiles with different head shape factors were developed to numerically simulate the processes of the truncated cone-shaped projectiles vertically penetrating through the fluid at different initial water-entry velocities. The velocity attenuation characteristics were obtained for the projectiles vertically penetrating through the fluid. The above results display that the resistance factor of a high-velocity projectile moving in the fluid is related with the projectile shape and the ratio of the instantaneous velocity of the projectile to its initial water-entry velocity. Based on the numerical simulations and the corresponding fitting results, an empirical formula was proposed by considering the projectile head factor to predict the velocity attenuation of the truncated cone-shaped projectiles vertically penetrating through the fluid. And a series of calculation examples were carried out to verify the formula. These calculation examples show that the formula is feasible and can be used to accurately calculate the velocity attenuation of high-velocity projectiles in fluid media and it is helpful for the structural design of the protective tanks of warships.
- truncated cone-shaped projectile /
- head shape factor /
- velocity attenuation /
- protective tank

HTML全文

韧性金属柱壳在高能炸药驱动下的膨胀断裂(简称外爆)是一个快速而复杂的动态过程。从炸药起爆到柱壳破坏只需数十微妙, 在此过程中柱壳经受了强加载、高应变率、大变形和复杂的加卸载。

R.W.Gurney^[1]给出了外爆下破片速度的经验公式并被广泛使用。G.I.Taylor^[2]根据高速摄影的实验图像, 提出了著名的Taylor断裂判据。C.R.Hoggatt等^[3]在回收的柱壳破片内壁面附近发现了绝热剪切带, 在Taylor思想的基础上, 提出外爆加载下在柱壳内壁面先生成剪切带并逐步向外扩展, 而后在拉应力控制下, 裂纹沿着剪切带失稳面由外向内传播的破坏机制。胡八一等^[4]对外爆中发生的剪切带进行了系统研究。胡海波等^[5]分析了柱壳破坏中出现的“单旋”剪切现象。汤铁钢等^[6]开展了一系列实验, 研究柱壳破坏的应变率效应。李永池等^[7]提出损伤演化方程, 模拟了中心线起爆方式下柱壳的破坏过程。

实验研究中, 炸药的起爆方式主要为单点起爆和平面起爆。但在理论分析和数值模拟中, 多数研究按照中心线起爆处理, 深入研究不同起爆方式对柱壳破裂影响的文章较少。本文中, 研究这3种起爆方式对柱壳破裂的影响。

1. 方法

外爆研究中, 炸药的起爆方式有3种:单点起爆、平面起爆和中心线起爆。单点起爆是指在炸药的一端插入雷管, 点燃雷管引爆炸药, 如图 1所示。平面起爆是指炸药一端各点同时起爆, 中心线起爆是指炸药柱中心线各点同时起爆。

图 1 实验装置示意图

Figure 1. Schematic of the experimental set-up

下载: 全尺寸图片幻灯片

对不同的起爆方式, 金属柱壳所经受的加载历程显然不同。但无论在理论研究还是在数值分析中, 为了便于理论推导和细化网格, 大多忽略结构的轴向影响, 把柱壳的破坏简化为二维平面应变问题, 并采用中心点起爆方式起爆炸药^[6-8]。这种简化对中心线起爆的柱壳中部是合理的, 但大部分实验采用的是单点起爆和平面起爆^[4-6]。此时这种简化是否也合理, 不同的起爆方式下柱壳所经受的加载历程又有何不同?本文中将对这些问题进行相应的数值模拟, 分析起爆方式对柱壳破裂的影响。

分别对3种不同的起爆方式, 使用Ls-Dyna软件模拟炸药爆轰驱动柱壳膨胀的过程。为了简化计算, 只计算了结构的四分之一, 计算模型如图 2。网格尺度为1 mm×1 mm×1 mm。

图 2 数值计算建模图

Figure 2. Geometry of the specimen used for simulation

下载: 全尺寸图片幻灯片

金属柱壳为4340钢, 强度采用Johnson-Cook模型:

$Y=\left(A+B \bar{\varepsilon}_{\mathrm{p}}^{n}\right)\left(1+C \ln \frac{\dot{\bar{\varepsilon}}_{\mathrm{p}}}{\dot{\bar{\varepsilon}}_{0}}\right)\left(1-\left(\frac{T-T_{0}}{T_{\mathrm{m}}-T_{0}}\right)^{m}\right)$

(1)

式中:A、B、n、C和m是材料参数; T₀为初始温度, T_m为熔点, 和分别为等效塑性应变和等效塑性应变率。4340钢的参数分别为:密度ρ=7 830 kg/m³, 体积模量K=164 GPa, 剪切模量

炸药为固态TNT, 爆轰产物采用JWL状态方程:

$p=A\left(1-\frac{\omega}{R_{1} V}\right) \mathrm{e}^{-R_{1} V}+B\left(1-\frac{\omega}{R_{2} V}\right) \mathrm{e}^{-R_{2} V}+\frac{\omega E}{V}$

(2)

式中:V为比体积, p为压力。TNT炸药的参数分别为:密度ρ=1 640 kg/m³, 爆速D=6 930 m/s, 爆压p_CJ=27 GPa, A=374 GPa, B=3.23 GPa, R₁=4.15, R₂=0.95, ω=0.3。

2. 结果与分析

柱壳轴向长度为12 cm, 分别提取柱壳内表轴向2、4、6、8和10 cm处的径向位移曲线和爆轰产物对柱壳的压力曲线, 处理数据得到炸药对柱壳内表初始单位面积的累积做功曲线。为了便于比较分析, 对累积做功曲线进行平移, 使爆轰产生的压力到达柱壳内表相应位置的时刻为零时刻。

2.1 平面起爆

平面起爆得到的曲线如图 3所示。柱壳内表所受压力载荷峰值在2 cm处由于边界效应的影响略低, 其他各处的压力峰值接近。虽然各处压力峰值接近, 但卸载速度不同。接近起爆端处, 直接接触空气, 压力到达峰值后卸载速度快, 炸药对柱壳的做功也小于远离起爆点处。在20 μs时, 8 cm处的做功是2 cm处的约1.73倍。

图 3 平面起爆的压力和累积做功曲线

Figure 3. The pressure and acting histories in planar detonation

下载: 全尺寸图片幻灯片

2.2 单点起爆

单点起爆得到的曲线如图 4所示。距离起爆点越远, 柱壳内表所受压力载荷峰值越小。因此在初期(约2 μs以前), 做功也随距离起爆点的位置越远而越小。但是由于距离起爆点近的位置, 压力卸载速度快, 因此在后期, 做功随距离起爆点的位置越远而越大。在20 μs时, 8 cm处的做功是2 cm处的约1.56倍。

图 4 单点起爆的压力和累积做功曲线

Figure 4. The pressure and acting histories in point detonation

下载: 全尺寸图片幻灯片

2.3 中心线起爆

中心线起爆得到的曲线如图 5所示。2和10 cm处的曲线几乎重合, 4、6和8 cm处的曲线几乎重合, 各处的内表压力峰值接近。在6.2 μs(压力曲线中为9.0 μs)时, 两侧稀疏波到达2和10 cm处, 此时该两处的压力, 特别是单位面积做功明显开始低于其他各处。炸药由中心线起爆后, 沿径向向外爆轰, 到达与柱壳的接触面, 冲击波在接触面反射, 向中心线会聚, 而后再向外反射传播, 在13.8 μs(压力曲线中为16.6 μs)时, 再次到达与柱壳的接触面, 二次加载柱壳。此时, 压力升高, 单位面积的做功也快速增加。这种重加载在平面起爆和单点起爆中没有出现。在20 μs时, 8 cm处的做功是2 cm处的约1.27倍。

图 5 中心线起爆的压力和累积做功曲线

Figure 5. The pressure and acting histories in centerline detonation

下载: 全尺寸图片幻灯片

2.4 不同起爆方式的差别

分别取平面起爆、单点起爆和中心线起爆下柱壳内表轴向8 cm处的数据, 对比3种不同的起爆方式下爆轰气体对柱壳作用的差异, 得到的曲线如图 6所示。

图 6 不同起爆方式的压力和累积做功曲线

Figure 6. The pressure and acting histories in different detonating methods

下载: 全尺寸图片幻灯片

由柱壳内表的压力曲线可以看出, 不同起爆方式下柱壳内表的压力载荷历程差异显著。中心线起爆对应的压力峰值最大, 其次是单点起爆, 平面起爆下的压力峰值最小。达到峰值压力后, 不同起爆方式下压力的卸载速度也存在明显差异, 峰值压力最高的中心线起爆对应的压力卸载速度最快, 其次是单点起爆, 平面起爆下柱壳内表压力卸载速度最慢。根据应力波理论分析可知, 该压力峰值的差别是由爆轰波对金属柱壳的入射角不同所造成的。中心线起爆方式下, 爆轰波对金属柱壳正入射, 爆轰波阵面平行于柱壳内表面。单点起爆和平面起爆下, 爆轰波对金属柱壳均为斜入射, 而平面起爆下入射角最大, 因此压力峰值也最小。

从做功曲线可以看出, 对轴向8 cm处的数据, 在20 μs时中心线起爆比平面起爆下的做功低3.21%。即到过程末期, 不同起爆方式引起的总做功差异不大。但是做功的累积过程却存在显著差异。在10 μs时, 中心线起爆比平面起爆下的做功低19.5%。

由以上分析可以看出, 不同起爆方式下炸药对柱壳内表的压力、做功历史差别明显。即使同一起爆方式下, 轴向不同位置的压力、做功历史也存在显著不同。因此, 在数值模拟和理论分析中, 必须考虑起爆方式和轴向位置的影响。

3. 半解耦数值方法及分析

对外爆的数值模拟中, 炸药网格占总网格的比例很大, 计算效率低下^[7-8]; 另外, 粗糙的网格划分不能描述柱壳典型破坏中出现的剪切带(实验中剪切带宽度约50 μm), 甚至不能区分两个相邻的剪切带(实验中剪切带间距约1 mm)^[5-6]; 此外起爆方式和轴向位置对柱壳破裂带来的影响也不容忽视。基于此, 本文中提出一种炸药和柱壳计算半解耦的数值计算方法。

3.1 半解耦方法

首先, 对三维模型(包括炸药和柱壳)划分较粗网格, 不引入金属柱壳的破坏模型, 得到爆轰产物对柱壳内表的压力曲线。然后, 去除炸药, 作平面应变假定进行二维计算, 划分精细网格, 将粗网格得到的内表压力曲线施加于柱壳内表面, 同时引入金属柱壳的破坏模型。这样, 一方面充分细化柱壳网格使之能模拟剪切带等破坏现象, 另一方面可以准确研究不同起爆方式下柱壳不同位置的破坏。

3.2 断裂能量耗散

设柱壳破坏过程中形成的剪切带条数为N, 产生单位面积剪切断裂面耗能为λ, 柱壳向外运动的速度为v, 柱壳的厚度为δ, 密度为ρ, 柱壳内径为r, 剪切带间距为则形成新断裂面的断裂耗能与柱壳的动能比为:

$\frac{E_{\mathrm{s}}}{E_{\mathrm{k}}}=\frac{\sqrt{2} \lambda \delta N}{0.5 \rho \delta \times 2 \pi r v^{2}}=\frac{2 \lambda}{\rho L v^{2}}$

(3)

汤铁钢等^[6]对图 2所示尺寸的装置进行实验, 使用高速摄影技术拍摄了整个膨胀过程, 使用的金属材料是钢, 密度ρ=7 830 kg/m³, 测得剪切带间距L≈1 mm, 末速度v≈1.8 km/s, 单位表面积的表面能^[8]约16 KJ/m²。计算得E_s/E_k≈0.126%, 可见断裂面的表面能所引起的能量重分配对柱壳速度的影响微乎其微, 完全可以忽略不计。

3.3 柱壳膨胀断裂过程中的塑性应变能

分析柱壳在膨胀过程中的塑性应变能, 设初始时刻柱壳内径为r₀, 外径为R₀, 末状态柱壳内径为r_f, 外径为R_f。设R₀/r₀=p, R_f/r₀=q。忽略弹性应变, 做平面应变假定, 即ε_r=-ε_θ。r方向总应变则总塑性功为:

$\begin{array}{*{20}{c}} {W_{\rm{t}}} = L\int_0^{2\pi } {\rm{d}} \theta \int_{{r_0}}^{{R_0}} 2 \sigma {\varepsilon _r}r{\rm{d}}r = \\ \frac{M \sigma\left(q^{2} \ln q^{2}-p^{2} \ln p^{2}-\left(q^{2}-p^{2}+1\right) \ln \left(q^{2}-p^{2}+1\right)\right)}{\rho_{\mathrm{m}}\left(p^{2}-1\right)} \end{array}$

(4)

式中:L为柱壳的轴向长度, ρ_m为金属的密度, σ为金属流动应力, 假设变形过程中流动应力保持恒值为柱壳总质量。

在图 2构型中, 在柱壳贯穿破坏时刻, 取q=2^[2-3], 则总塑性功W_t≈75 kJ。应该注意到, 在柱壳破坏开始之前, 已经消耗大量塑性功, 在破坏开始时刻, 取q=1.7, 则破坏开始前塑性功W_a=52 kJ。因此材料破坏过程中, 如不引入破坏模型, 带来的塑性功W_i=W_t-W_a=23 kJ。柱壳动能E_k=Mv²/2=842 kJ。总塑性功与总动能的比为8.9%, 但不引入破坏模型带来的塑性功与柱壳总动能的比仅为2.7%。所以, 破坏过程中的塑性功所引起的能量分配对柱壳速度的影响较小, 可忽略。

3.4 半解耦法计算结果

对断裂耗能和塑性应变能的分析表明, 破坏产生新断面的能量耗散和破坏过程中的塑性应变能对柱壳动能的影响较小, 可忽略。因此半解耦数值方法是合理的, 有理论依据的。采用此半解耦数值方法, 模拟图 2结构在单点起爆方式下轴向6 cm处壳体的破坏过程, 网格尺度为90 μm×90 μm。局部放大破坏图像如图 7所示, 它成功地再现了炸药驱动下柱壳膨胀断裂过程中出现的绝热剪切带。

图 7 半解耦法计算结果

Figure 7. Modeling result using semi-decoupling method

下载: 全尺寸图片幻灯片

4. 总结

采用数值模拟方法, 研究了3种不同起爆方式下爆轰气体产物对韧性金属柱壳内表的加载历程, 发现不同起爆方式下柱壳内表所经受的压力历程差别巨大, 中心线起爆峰值压力最大, 单点起爆其次, 平面起爆最小。理论分析认为, 造成这种差别的原因是爆轰波对柱壳的入射角不同。即使在同一起爆方式下, 柱壳内表轴向不同点的压力载荷历程也存在巨大差异。因此, 在理论分析和数值模拟中, 必须考虑起爆方式和轴向位置的影响。提出一种模拟炸药驱动下韧性金属柱壳膨胀破坏的半解耦数值方法, 该方法可以实现对柱壳精细网格划分条件下的破裂模拟, 并能研究不同起爆方式的影响。通过理论分析说明了这种半解耦方法的合理性。

图 1 截锥形弹体示意图

Figure 1. Schematic diagram of a truncated cone-shaped projectile

下载: 全尺寸图片幻灯片

图 2 不同头形因数截锥形弹体有限元模型

Figure 2. Finite element models for truncated cone-shaped projectiles with different head type coefficients

下载: 全尺寸图片幻灯片

图 3 针对弹体直径为12.65 mm，长度为25.4 mm，入水初速度为603 m/s时，试验与数值计算的空泡尺寸

Figure 3. Comparison of cavitation sizes obtained experimentally and numerically for the projectile with the diameter of 12.65 cm and the length of 25.4 cm, water entering at 603 m/s

下载: 全尺寸图片幻灯片

图 4 长度为25.4 mm的弹体在两种入水速度工况下位移和速度的变化

Figure 4. Changes of displacement and velocity of the projectile with the length of 25.4 mm at two initial velocities of water entry

下载: 全尺寸图片幻灯片

图 5 长度为38.1 mm的弹体在两种入水速度工况下位移和速度的变化

Figure 5. Changes of displacement and velocity of the projectile with the length of 38.1 mm at two initial velocities of water entry

下载: 全尺寸图片幻灯片

图 6 不同头形因数的截锥形弹体在不同入水速度工况下速度随入水距离的衰减

Figure 6. Velocity attenuation of truncated cone-shaped projectiles with different head coefficients with water-entry distance at different initial velocities of water entry

下载: 全尺寸图片幻灯片

图 7 不同入水速度下， $\psi = 0$ 的截锥形弹体阻力因数与 ${v_{\rm{b}}}/{v_0}$ 的关系

Figure 7. Resistance factor varying with ${v_{\rm{b}}}/{v_0}$ for the truncated cone-shaped projectile with $\psi = 0$ at different water-entry velocities

下载: 全尺寸图片幻灯片

图 8 对于头形因数不同的截锥形弹体，参数a₀和a₁与马赫数 $Ma$ 的关系

Figure 8. Changes of parameters a₀ and a₁ with Mach number $Ma$ for the truncated cone-shaped projectiles with different head shape factors

下载: 全尺寸图片幻灯片

图 9 $a$ 系列参数与头形因数 $\psi$ 的关系

Figure 9. Series parameters of a varyied with head shape factor $\psi$

下载: 全尺寸图片幻灯片

图 10 入水速度为400 m/s，头形因数不同的弹丸，速度随入水距离衰减的模拟结果与公式计算结果的比较

Figure 10. Comparison of velocity attenuation with distance between numerical simulation and formula calculation for the projectiles with different head shape factors at the initial water-entry velocity of 400 m/s

下载: 全尺寸图片幻灯片

图 11 头形因数为1/6，入水速度不同的弹丸，速度随入水距离衰减的模拟结果与公式计算结果的比较

Figure 11. Comparison of velocity attenuation with distance between numerical simulation and formula calculation for the projectile with the head shape factor of 1/6 at different initial water-entry velocities

下载: 全尺寸图片幻灯片

表 1 水的状态方程各项参数^[14]

Table 1. Parameters of equation of state for water^[14]

${A_1}{\rm{/GPa}}$	${A_{\rm{2}}}{\rm{/GPa}}$	${A_3}{\rm{/GPa}}$	${B_0}$	${B_1}$	${T_1}{\rm{/GPa}}$	${T_2}{\rm{/GPa}}$	${\rho }_{0}\rm{/(kg}·{\rm{m}}^{\rm{-3}}\rm{)}$
2.2	9.54	14.57	0.28	0.28	2.2	0	1 000

下载: 导出CSV

表 2 数值计算工况

Table 2. Numerical calculation conditions

工况	$\psi$	v₀/(m·s⁻¹)	工况	$\psi$	v₀/(m·s⁻¹)	工况	$\psi$	v₀/(m·s⁻¹)	工况	$\psi$	v₀/(m·s⁻¹)
1-1	0	400	2-1	1/3	400	3-1	2/3	400	4-1	1	400
1-2	0	650	2-2	1/3	650	3-2	2/3	650	4-2	1	650
1-3	0	900	2-3	1/3	900	3-3	2/3	900	4-3	1	900
1-4	0	1 150	2-4	1/3	1 150	3-4	2/3	1 150	4-4	1	1 150
1-5	0	1 400	2-5	1/3	1 400	3-5	2/3	1 400	4-5	1	1 400

下载: 导出CSV

表 3 不同工况下的参数 ${a_0}$ 和 ${a_1}$ 的数值

Table 3. Values of parameters ${a_0}$ and ${a_1}$ under different working conditions

$\psi$	拟合参数	v₀/(m·s⁻¹)
$\psi$	拟合参数	400	650	900	1150	1400
0	a₀	0.617	0.612	0.633	0.587	0.367
0	a₁	−0.217	−0.212	−0.223	−0.153	0.128
1/3	a₀	0.329	0.459	0.477	0.440	0.288
1/3	a₁	−0.019	−0.158	−0.168	−0.113	0.077
2/3	a₀	0.768	0.707	0.784	0.714	0.553
2/3	a₁	−0.387	−0.297	−0.394	−0.289	−0.082
1	a₀	1.216	1.204	1.527	1.382	1.416
1	a₁	−0.344	−0.345	−0.823	−0.580	−0.640

下载: 导出CSV

表 4 参数 ${a_0}$ 、 ${a_1}$ 的拟合结果

Table 4. Fitting results of parameters ${a_0}$ and ${a_1}$

$\psi$	a₀₁	a₀₂	a₀₃	a₀₄	a₀₅	a₁₁	a₁₂	a₁₃	a₁₄	a₁₅
0	0.847	−1.459	2.575	−0.707	−0.887	−0.376	0.918	−1.096	−1.059	1.743
1/3	−0.611	6.078	−13.025	12.901	5.055	1.060	−7.042	15.285	−15.196	5.971
2/3	2.810	−15.067	37.573	−38.342	13.579	−3.246	21.137	−52.932	54.463	−19.504
1	8.033	−53.017	141.081	−153.418	58.737	−10.547	79.685	−213.159	232.929	−89.548

下载: 导出CSV

表 5 ${a_{ij}}$ 系列参数的拟合结果

Table 5. Fitting results of ${a_{ij}}$ series parameters

${a_{ij}}$	${a_{ij}}_0$	${a_{ij}}_1$	${a_{ij}}_2$
${a_0}_1$	0.693	−7.536	15.031
${a_{02}}$	−0.864	49.585	−102.332
${a_{03}}$	1.908	−128.117	267.957
${a_{04}}$	−0.654	136.681	−289.498
${a_{05}}$	0.814	−32.000	88.221
${a_1}_1$	−0.239	9.212	−19.657
${a_1}_2$	0.628	−70.273	149.620
${a_1}_3$	−1.463	185.983	−397.312
${a_1}_4$	0.187	−201.803	433.298
${a_1}_5$	1.001	77.280	−167.088

下载: 导出CSV

参考文献(14)

[1]	孔祥韶. 爆炸载荷及复合多层防护结构响应特性研究[D]. 武汉: 武汉理工大学, 2013: 1-26.DOI: 10.7666/d.Y2364126.
[2]	LEE M, LONGORIA R G, WILSON D E. Cavity dynamics in high-speed water entry [J]. Physics of Fluids, 1997, 9(3): 540–550. DOI: 10.1063/1.869472.
[3]	李营, 张磊, 朱海清, 等. 爆炸破片在液舱中的速度衰减特性研究 [J]. 中国造船, 2016, 57(1): 127–137. DOI: 10.3969/j.issn.1000-4882.2016.01.014. LI Y, ZHANG L, ZHU H Q, et al. Velocity attenuation of blast fragments in water tank [J]. Shipbuilding of China, 2016, 57(1): 127–137. DOI: 10.3969/j.issn.1000-4882.2016.01.014.
[4]	沈晓乐, 朱锡, 侯海量, 等. 高速破片侵彻防护液舱试验研究 [J]. 中国舰船研究, 2011, 6(3): 12–15. DOI: 10.3969/j.issn.1673-3185.2011.03.003. SHEN X L, ZHU X, HOU H L, et al. Experimental study on penetration properties of high velocity fragment into safety liquid cabin [J]. Chinese Journal of Ship Research, 2011, 6(3): 12–15. DOI: 10.3969/j.issn.1673-3185.2011.03.003.
[5]	郭子涛. 弹体入水特性及不同介质中金属靶的抗侵彻性能研究[D]. 哈尔滨: 哈尔滨工业大学, 2012: 21-36.DOI: 10.7666/d.D241209.
[6]	ZHAO B L, ZHAO J G, CUI C Y, et al. Growth model of cavity generated by the projectile impacting liquid-filled tank [J]. Defence Technology, 2020, 16(3): 609–616. DOI: 10.1016/j.dt.2019.09.013.
[7]	ZHANG Y, LI X B, LI S Y. Research on the velocity attenuation characteristics of the fragments during high-speed water entry [C] // Proceedings of the 37th International Conference on Ocean, Offshore and Arctic Engineering. Madrid: ASME, 2018.DOI: 10.1115/OMAE2018-78665.
[8]	VARAS D, ZAERA R, LÓPEZ-PUENTE J. Numerical modelling of partially filled aircraft fuel tanks submitted to Hydrodynamic Ram [J]. Aerospace Science and Technology, 2012, 16(1): 19–28. DOI: 10.1016/j.ast.2011.02.003.
[9]	BIRKHOFF G, CAYWOOD T E. Fluid flow patterns [J]. Journal of Applied Physics, 1949, 20(7): 646–659. DOI: 10.1063/1.1698450.
[10]	孔祥韶, 吴卫国, 刘芳, 等. 舰船舷侧防护液舱对爆炸破片的防御作用研究 [J]. 船舶力学, 2014, 18(8): 996–1004. DOI: 10.3969/j.issn.1007-7294.2014.08.015. KONG X S, WU W G, LIU F, et al. Research on protective effect of guarding fluid cabin under attacking by explosion fragments [J]. Journal of Ship Mechanics, 2014, 18(8): 996–1004. DOI: 10.3969/j.issn.1007-7294.2014.08.015.
[11]	LECYSYN N, BONY-DANDRIEUX A, APRIN L, et al. Experimental study of hydraulic ram effects on a liquid storage tank: analysis of overpressure and cavitation induced by a high-speed projectile [J]. Journal of Hazardous Materials, 2010, 178(1/2/3): 635–643. DOI: 10.1016/j.jhazmat.2010.01.132.
[12]	JOHNSON G R, COOK W H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures [C] // Proceedings of the 7th International Symposium on Ballistics. Netherlands: The Hague, 1983.
[13]	陈刚, 陈忠富, 陶俊林, 等. 45钢动态塑性本构参量与验证 [J]. 爆炸与冲击, 2005, 25(5): 451–456. DOI: 10.11883/1001-1455(2005)05-0451-06. CHEN G, CHEN Z F, TAO J L, et al. Investigation and validation on plastic constitutive parameters of 45 steel [J]. Explosion and Shock Waves, 2005, 25(5): 451–456. DOI: 10.11883/1001-1455(2005)05-0451-06.
[14]	李晓杰, 张程娇, 王小红, 等. 水的状态方程对水下爆炸影响的研究 [J]. 工程力学, 2014, 31(8): 46–52. DOI: 10.6052/j.issn.1000-4750.2013.03.0180. LI X J, ZHANG C J, WANG X H, et al. Numerical study on the effect of equations of state of water on underwater explosions [J]. Engineering Mechanics, 2014, 31(8): 46–52. DOI: 10.6052/j.issn.1000-4750.2013.03.0180.

施引文献

资源附件(0)

访问统计

图(11) / 表(5)

计量

文章访问数: 637
HTML全文浏览量: 250
PDF下载量: 78
被引次数: 0

1. 方法
2. 结果与分析
2.1 平面起爆
2.2 单点起爆
2.3 中心线起爆
2.4 不同起爆方式的差别
3. 半解耦数值方法及分析
3.1 半解耦方法
3.2 断裂能量耗散
3.3 柱壳膨胀断裂过程中的塑性应变能
3.4 半解耦法计算结果
4. 总结

截锥形弹体在液体介质中运动速度衰减规律分析

doi: 10.11883/bzycj-2020-0075

作者简介: 孔祥韶（1983－ ），男，博士，教授，博士生导师，kongxs@whut.edu.cn

计量

出版历程