• 中国科技核心期刊
  • JST收录期刊
  • Scopus收录期刊
  • DOAJ收录期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

水下爆炸冲击波数值仿真精度研究

敖启源 卢熹 姜智雅 康珀阁

敖启源, 卢熹, 姜智雅, 等. 水下爆炸冲击波数值仿真精度研究[J]. 水下无人系统学报, 2024, 32(1): 158-165 doi: 10.11993/j.issn.2096-3920.2023-0098
引用本文: 敖启源, 卢熹, 姜智雅, 等. 水下爆炸冲击波数值仿真精度研究[J]. 水下无人系统学报, 2024, 32(1): 158-165 doi: 10.11993/j.issn.2096-3920.2023-0098
AO Qiyuan, LU Xi, JIANG Zhiya, KANG Poge. Numerical Simulation Accuracy Study of Underwater Explosion Shock Waves[J]. Journal of Unmanned Undersea Systems, 2024, 32(1): 158-165. doi: 10.11993/j.issn.2096-3920.2023-0098
Citation: AO Qiyuan, LU Xi, JIANG Zhiya, KANG Poge. Numerical Simulation Accuracy Study of Underwater Explosion Shock Waves[J]. Journal of Unmanned Undersea Systems, 2024, 32(1): 158-165. doi: 10.11993/j.issn.2096-3920.2023-0098

水下爆炸冲击波数值仿真精度研究

doi: 10.11993/j.issn.2096-3920.2023-0098
详细信息
    作者简介:

    敖启源(1999-), 男, 在读硕士, 主要研究方向为水下爆炸

    通讯作者:

    卢 熹(1983-), 男, 博士, 副教授, 主要从事水下高效毁伤技术研究.

  • 中图分类号: TJ630.2; U674

Numerical Simulation Accuracy Study of Underwater Explosion Shock Waves

  • 摘要: 在水下爆炸数值仿真研究中, 网格尺寸和一次项人工粘性系数对冲击波峰值压力计算结果有较大影响。在预定计算精度条件下, 快速确定网格尺寸及人工粘性对数值计算意义重大。为此, 文中基于LS-DYNA有限元软件, 建立78 g三硝基甲苯(TNT)二维水下爆炸数值计算模型, 重点分析网格尺寸和一次项粘性系数对水下爆炸冲击波峰值压力和整体计算误差的影响规律。结果表明, 随着网格密度因子的增大, 计算峰值压力对网格的敏感性降低, 且网格密度较大时, 过小的一次项系数会导致计算峰值压力与经验公式值的相对误差增大。在此基础上获得20%范围内误差与网格尺寸、粘性系数之间的关系, 并构建出可用于快速确定网格尺寸和一次项人工粘性系数的误差预估模型, 通过0.2 ~5 000 kg范围内的TNT柱形装药(长径比为1)和球形装药的水下爆炸计算, 验证了预估模型的普适性, 可为二维中近场范围内的水下爆炸冲击波数值仿真计算研究提供参考。

     

  • 水下武器作为舰船生命力的主要威胁之一, 其爆炸冲击波及气泡载荷会对舰船造成严重的局部和总体破坏[1]。各国相继开展了大量水下爆炸试验, 但实弹实测试验安全风险较高、成本高昂、重复性低且观测范围有限[2]。随着计算机性能和仿真技术的发展, 数值仿真以其较高计算精度、低成本和高可重复性等优点成为继实验、理论研究后第3种水下爆炸研究方法[3]。因此开展水下爆炸数值仿真精度研究具有十分重要的意义。

    数值仿真研究水下爆炸冲击波载荷问题时, 为了处理冲击波的强间断面, 抑制波阵面前后的数值振荡, 引入了人工粘性。由于人工粘性的引入会在几个网格宽度上光滑冲击波波阵面, 因此水下爆炸冲击波载荷计算结果直接受人工粘性的影响。同时,人工粘性的引入要求网格尺寸不能过大, 否则计算过程中网格会忽略部分波阵面信息, 致使峰值过低影响计算精度。Huang等[4]通过典型三硝基甲苯(Trinitrotolue,TNT)炸药水下爆炸数值分析, 探讨了一次、二次粘性系数对数值计算结果的影响, 并给出了一定比例爆距范围内的建议值。Wang等[5]研究了不同装药质量下网格尺寸对冲击波峰值的影响, 并引入与装药半径和单元边长相关的无因次量表征网格尺寸。胡亮亮等[6]以常规TNT水下爆炸为例, 对水的状态方程、人工粘性系数和网格尺寸对于数值仿真结果的影响进行了研究。张社荣等[7]基于有限元软件AUTODYN建立了不同炸药量的水下爆炸数值模型, 对比分析了网格尺寸对不同爆距处冲击波峰值压力及比冲量的影响。此外还有其他学者[8-13]讨论了网格尺寸及粘性系数对水下爆炸计算结果的影响, 但现有研究无法在预定精度下快速确定网格尺寸和人工粘性系数。因此文中以TNT水下自由场爆炸数值计算为例, 探究网格尺寸和一次项人工粘性系数对水下爆炸冲击波峰值压力的影响, 综合分析不同工况下网格密度因子和一次项系数与峰值压力平均误差间的关系, 构建出普适性较高的水下爆炸数值误差预估模型, 为预定精度的仿真模型设计提供依据。

    基于文献[14]中开展的爆炸水井试验, 建立二维轴对称计算模型。试验水域尺寸4.5 m×9 m, 空气域尺寸4.5 m×0.1 m, 炸药为直径2 cm的等高药柱, 炸药质量78 g, 放置在水深4 m位置处, 采用中心起爆方式。网格的排列走向和过渡方式会影响计算结果, 因此采用均匀网格划分方式。在距装药中心0.4~2.8 m/kg1/3比例爆距内选定16个观测点, 采用关键字*INITIAL_HYDROSTATIC_ALE初始化静水压力, 设置水域压力梯度, 以模拟真实条件下的水下压强环境。炸药、空气和水介质均选用ALE(Arbitrary Lagrangian-Eulerian)算法, 数值计算模型如图1所示。

    图  1  数值计算模型
    Figure  1.  Numerical calculation model

    装药选用典型的单质TNT炸药, 采用JWL(Jones Wilkins Lee)状态方程描述其爆炸过程, 具体形式为

    $$ P = A\left( {1 - \frac{\omega }{{{R_1}V}}} \right){{\text{e}}^{ - {R_1}V}} + B\left( {1 - \frac{\omega }{{{R_2}V}}} \right){{\text{e}}^{ - {R_2}V}} + \frac{{\omega E}}{V} $$ (1)

    式中: V为相对体积; ABR1R2为常数, 取值如表1[15]所示; P为压力; ω为药量; E 为单位体积内能。

    表  1  TNT状态方程参数
    Table  1.  State equation parameters for TNT
    ${A / { {\text{GPa} }} }$${B / { {\text{GPa} }} }$$ {R_1} $$ {R_2} $$ \omega $$ {E / {{\text{GPa}}}} $
    371.203.234.150.950.307.00
    下载: 导出CSV 
    | 显示表格

    数值仿真中常用POLYNOMIAL和GRUNEISEN状态方程描述水介质, 但所使用的压力梯度关键字对POLYNOMIAL状态方程参数设置有一定要求, 因此水介质采用GRUNEISEN状态方程描述, 其形式根据水所处状态而定。

    当水介质处于压缩状态时, 其状态方程为

    $$\begin{align} P =&\; \dfrac{{{\rho _0}{C^2}\left[ {1 + \left( {\left. {1 - {v_0}/2} \right)} \right.\mu - \left( {\alpha /2} \right){\mu ^2}} \right]}}{{{{\left[ {1 - \left( {{S_1} - 1} \right)\mu - {S_2}{\mu ^2}/\left( {\mu + 1} \right) - {S_3}{\mu ^3}/{{\left( {\mu + 1} \right)}^2}} \right]}^2}}} + \\&\;\left( {{\gamma _0} + \alpha \mu } \right)E \end{align} $$ (2)

    当水介质处于膨胀状态时, 其状态方程为

    $$ P = {\rho _0}{C^2}\mu + \left( {{\gamma _0} + \alpha \mu } \right)E $$ (3)

    式中: µ为水的压缩比; C为水中声速; S1S2S3为常数; γ0为GRUNEISEN系数; α为体积修正系数; V0为初始相对体积。以上参数取值如表2所示。

    表  2  水状态方程参数
    Table  2.  State equation parameters for water
    C/(m/s)$ {S_1} $$ {S_2} $$ {S_3} $$ \alpha $$E$$ {V_0} $
    14832.106−0.1740.011.0
    下载: 导出CSV 
    | 显示表格

    空气使用POLYNOMIAL状态方程进行描述, 其形式为

    $$ P = {C_0} + {C_1}\mu + {C_2}{\mu ^2} + {C_3}{\mu ^3} + \left( {{C_4} + {C_5}\mu + {C_6}{\mu ^2}} \right) $$ (4)

    式中: C0C1C2C3C4C5C6为常数。取值如表3所示。

    表  3  空气状态方程参数
    Table  3.  State equation parameters for air
    E/(J/kg)C0C1C2C3C6$ {C_4} $$ {C_5} $
    2.5×106000000.40.4
    下载: 导出CSV 
    | 显示表格

    根据文献[10], 长径比为1∶1的柱形装药可以近似为球形装药, 对于球形装药水下自由场爆炸冲击波的传播, Cole[16]通过大量试验标定了水下爆炸冲击波相似律的公式系数, 获得了不考虑水深影响的TNT炸药水下爆炸冲击波峰值压力计算公式, Zamyshlyaev[17]在其基础上将经验公式修正为

    $$ {P_{\text{m}}} = \left\{ {\begin{array}{*{20}{c}} {44.1 \times {{\left( {\dfrac{{\sqrt[3]{W}}}{R}} \right)}^{1.5}}},&{6 \leqslant \dfrac{R}{{{R_e}}} \leqslant 12} \\ {52.4 \times {{\left( {\dfrac{{\sqrt[3]{W}}}{R}} \right)}^{1.13}}},&{12 \leqslant \dfrac{R}{{{R_e}}} \leqslant 240} \end{array}} \right. $$ (5)

    式中: Pm为冲击波峰值压力; R为爆距; W为装药质量; Re为装药半径。

    根据式(5)可计算出水中一定范围内的冲击波峰值压力, 即表4中的经验值。为获得水下爆炸冲击波传播演化规律, 文献[14]在爆炸水井中开展了水下爆炸冲击波试验, 获得了78 g TNT装药不同爆距Z处的冲击波峰值压力, 即表4中的试验值。基于爆炸水井试验, 利用LS-DYNA有限元软件对78 g TNT炸药水下爆炸过程进行仿真, 网格尺寸0.25 cm, 一次项系数取值0.06, 获得不同爆距处的冲击波峰压力, 即数值解。将试验值、数值解及经验值进行对比。

    表  4  不同爆距处峰值压力对比
    Table  4.  Comparison of peak pressure at different scaled blast distances MPa
    PmZ /(m/kg1/3)
    0.751.171.882.80
    试验值70.1039.6023.3015.30
    数值解75.5042.0022.7013.60
    经验值72.6043.8025.7016.40
    下载: 导出CSV 
    | 显示表格

    对比表4数据, 可以发现数值解与试验值最小误差仅2.6%, 平均误差6.9%; 与经验值最小误差4%, 平均误差9.2%。因此网格尺寸和一次项粘性系数取值合理时, 数值模型可获得较高的计算精度。经验值与试验值峰值压力平均误差7.2%, 表明经验公式可以较为准确地预估水下爆炸冲击波峰值压力。因此, 在探究网格尺寸及人工粘性对数值计算精度的影响时, 仅比较数值解与经验值。

    在有限元计算中为了反映冲击波波阵面的强间断, 引入了人工粘性来光滑冲击波, 这使得数值计算峰值压力低于真实值。一次项人工粘性系数对冲击波峰值压力影响较大, 二次项人工粘性系数主要用于抑制冲击波衰减过程中的虚假振荡, 但对水中爆炸数值仿真计算中冲击波峰值影响较小。参阅文献[8], 二次项人工粘性系数取定值1.0。LS-DYNA中人工粘性形式为

    $$ q = \left\{ {\begin{array}{*{20}{l}} {\rho l\left( {{Q_1}l{\varepsilon ^2} - {Q_2}C\varepsilon } \right)},&{\varepsilon \lt 0} \\ 0,&{\varepsilon \leqslant 0} \end{array}} \right. $$ (6)

    式中: ρ为材料密度; Q1为二次项人工粘性系数; Q2为一次项人工粘性系数; l为特征长度; C为当地声速; ε为体积变化率。

    在水下爆炸数值计算中, 网格密度和排列方式对计算结果影响很大, 过大的网格尺寸在计算过程中会忽略冲击波波阵面信息, 冲击波爬升至峰值所需时间变长; 而过小的网格尺寸对计算资源带来的压力也不可忽视。同时, 水域网格长宽比尽量接近于1, 尤其是炸药附近的网格, 否则爆炸冲击波易出现失真[3]。在探究网格尺寸对数值计算结果的影响时, 引入与装药半径 R0和网格尺寸 L0相关的无因次量, 因此文中为研究网格尺寸对数值计算结果的影响[5, 7-8], 其具体形式为

    $$ \lambda = \frac{{{R_0}}}{{{L_0}}} $$ (7)

    为探究水下爆炸数值仿真中网格尺寸和Q2对冲击波峰值压力的影响, 对78 g TNT水下自由场爆炸过程进行数值仿真计算。当Q2>0.1时, 计算峰值压力误差较大[8], 因此调整分别取Q2=0.02、0.03、0.04、0.05、0.06、0.07、0.08和0.10。对于二维模型而言, 网格密度因子λ=8时便满足大部分二维数值模型, 因此取λ分别取λ=1、2、3、4、5、6、7和8, 共计算64个工况。为了从全局意义上分析网格尺寸和Q2Pm的影响, 引入平均误差, 定义为比例爆距0.375~2.8 m/kg1/3范围内选定观测点峰值压力相对于经验公式误差的平均值。

    网格尺寸及Q2Pm结果如图2~图6所示。由图2图3可知, 当λ较小时, Pm较低, 与经验公式值偏差较大, 随着λ的增加, Pm显著增大且近场峰值压力与经验值吻合较好。但随着λ增加, Pm趋于稳定, 不会随着λ的变化出现显著变化, 这表明随着网格密度的增加, Pm对网格的敏感性越来越低。与此同时, 网格尺寸对冲击波超压爬升至峰值所需时间影响很大, λ越小, 冲击波超压爬升至峰值时间越长且峰值压力越低, 如图4所示。

    图  2  Q2=0.08时, 不同λPm随比例爆距变化曲线
    Figure  2.  Variation of Pm with scaled blast distance for different λ at Q2=0.08
    图  3  Q2=0.08时, 不同比例爆距处Pmλ变化曲线
    Figure  3.  Variation of Pm with λ for different scaled blast distances at Q2=0.08
    图  4  λPm上升速度影响
    Figure  4.  Effect of λ on the rate of increase of Pm
    图  5  λ=3时, 不同Q2Pm随比例爆距变化曲线
    Figure  5.  Variation of Pm with scaled blast distance for different Q2 at λ=3
    图  6  λ=3时, 不同比例爆距处PmQ2变化曲线
    Figure  6.  Variation of Pm with Q2 for different blast scaled distances at λ =3

    图5图6所见, 随着Q2逐渐变小, Pm逐渐增大, 且近场冲击波峰值压力受Q2影响较大。由于Q2取值范围较小, 不同比例爆距处PmQ2的变化曲线与文献[8]中有所不同。在实际计算中发现, 当网格密度较大时, 过小的Q2反而导致Pm数值解与经验值的偏差增大, 如表5所示。通过计算, Q2=0.02时不同爆距处的平均误差为11.2%, 而Q2=0.06时平均误差仅4.2%, 因此有必要探究λQ2对峰值压力平均误差的影响。

    表  5  λ=6时, 不同比例爆距处冲击Pm对比
    Table  5.  Different scaled blast distances peak pressure of shock waves at λ =6 MPa
    工况类型Z/(m/kg−1/3)
    0.3270.5000.7000.9001.1701.700
    Q2=0.020253.000138.00092.80067.50048.80030.700
    Q2=0.060235.000127.00081.80058.30041.50025.400
    经验值235.800124.70078.40059.00043.90028.800
    下载: 导出CSV 
    | 显示表格

    对比数值解和经验值在不同爆距处的冲击波峰值压力, 获得不同Q2下峰值压力平均误差EPλ的变化关系, 如图7所示。整体来看, 随着λ的增大和Q2的减小, EP逐渐降低。但以Q2=0.10为例, 随着λ的增大, EP先减小后增大, 这是因为网格密度较大时, 过小的Q2会加大伪振荡, 致使Pm数值解与经验的相对误差增大, 因此在数值计算中不能一味地减小网格尺寸和一次项系数。同时, 当λ=1时, EP很大, 即使调整Q2也未能使EP满足工程精度。因此在水下爆炸数值计算中应首先确定λ, 同时调整Q2方能得到精度较高的计算结果。

    图  7  峰值压力平均误差
    Figure  7.  Average errors of peak pressure

    为了便于应用, 实现在预定精度下快速确定网格尺寸和一次项系数, 需构建出关于λQ2的误差EP预估模型。考虑到实际工程需要, 误差应控制在20%以内。同时, 由前面分析可知, 较大的网格密度以及过小的一次项系数均可能导致计算误差增大。因此, 将误差预估模型的变量区间范围限定在0.03≤Q2≤0.10和3≤λ≤8, 则图7中的数据在该区间范围显示为图8中曲线。

    图  8  限定区间后峰值压力平均误差曲线
    Figure  8.  Average errors of peak pressure curves after limited interval

    图8中的EPλ取对数得到lgEP和lgλ的关系曲线如图9所示, k为斜率, b为截距。可以看到, 各曲线呈近似平行的线性关系, 因此可构造线性表达式为

    图  9  lgEP和lgλ关系曲线
    Figure  9.  The relationship curves between lgEP and lgλ
    $$ \lg {E_{\rm{P}}} = k\lg \lambda + b $$ (8)

    利用式(8)对图9中的数据进行拟合得到拟合系数如表6所示。其中, 复合相关系数R1均值0.979, 决定系数R2均值0.959, 可见数据线性拟合精度较高。表6中各曲线斜率k值比较接近, 为了得到归一化的误差预估模型, 可以取k值平均值, 并且将截距b看作关于Q2的函数。如图10所示为截距b与lgQ2的关系, 可以看到二者近似呈线性关系, 因此可构建线性关系式

    表  6  lgEP关于lgλ的拟合参数
    Table  6.  Fitting parameters for λ in EP
    拟合参数Q2
    0.030 00.040 00.050 00.060 00.070 00.080 00.1 000
    k−0.690 0−0.790 0−0.590 0−0.510 0−0.530 0−0.530 0−0.490 0
    b1.300 01.420 01.410 01.410 01.500 01.550 01.610 0
    R1 0.940 40.967 60.963 30.985 80.996 10.999 30.999 9
    R20.884 30.936 30.927 90.971 70.992 20.998 60.999 8
    下载: 导出CSV 
    | 显示表格
    图  10  截距-对数粘性系数线性拟合
    Figure  10.  Linear fitting of intercept-logarithm viscosity coefficient
    $$ b = {k_1}\lg {Q_2} + {b_1} $$ (9)

    利用式(9)对图10中数据进行线性拟合, 结果如表7所示, 可以看出数据具有较高拟合精度。将式(9)代入式(8), 整理得

    表  7  b关于Q2的拟合参数
    Table  7.  Fitting parameters for b in Q2
    k1b1RR2
    0.586 02.180 00.971 60.944 1
    下载: 导出CSV 
    | 显示表格
    $$ {E_{\rm{P}}} = 151.36\times{\left( {\frac{{{Q_2}^{1.895}}}{\lambda }} \right)^{0.38}} $$ (10)

    从式(10)可以看到, 误差EP是关于${{{Q_2}^{1.895}} \mathord{\left/ {\vphantom {{{Q_2}^{1.895}} \lambda }} \right.} \lambda } $的函数。利用图8中的数据绘制EP${{{Q_2}^{1.895}} \mathord{\left/ {\vphantom {{{Q_2}^{1.895}} \lambda }} \right.} \lambda } $的关系曲线如图11所示。可以看到, EP${{{Q_2}^{1.895}} \mathord{\left/ {\vphantom {{{Q_2}^{1.895}} \lambda }} \right.} \lambda } $近似呈线性关系。因此, 为了进一步提高误差预估模型的拟合精度, 对图11中数据重新进行线性拟合, 其中, 拟合模型的复合相关系数 R=0.998, 决定系数 R2=0.996, 说明拟合精度较高, 则式(11)为最终得到的误差预估模型

    图  11  峰值压力平均误差拟合结果
    Figure  11.  Fitting results of peak pressure average error
    $$ {E_{\rm{p}}} = 1\;413.57\times{\left( {\frac{{{Q_2}^{1.895}}}{\lambda }} \right)^{0.76}} + 4.06 $$ (11)

    为验证误差预估模型具有较高的普适性, 需要对不同药量和装药形状的炸药进行仿真计算, 同时结合研究背景, 对0.2、5、500、1 500 和5 000 kg TNT柱形装药(长径比为1)和球形装药水下爆炸进行计算。在距装药中心0.327~2.8 m/kg1/3比例爆距间设置16个观测点, 取Q2=0.06, λ=6。不同工况下的峰值压力平均误差及预估误差如表8所示。可以看到, 预估误差与实际计算误差相近, 因此误差预估模型对于不同装药量和装药形状也具有很高的适用性, 有助于建立水下爆炸模型对整体计算精度和网格数量进行综合分析, 为预定精度的仿真模型设计提供依据。

    表  8  不同工况下峰值压力误差
    Table  8.  Errors of peak pressure for different operating conditions
    装药质量/kg装药形状平均误差/%预估误差/%
    0.2柱形9.7510.36
    0.2球形10.79
    5.0柱形9.45
    5.0球形10.45
    50.0柱形10.83
    50.0球形10.30
    500.0柱形9.62
    500.0球形10.18
    1 500.0柱形9.57
    1 500.0球形9.89
    5 000.0柱形9.58
    5 000.0球形9.76
    下载: 导出CSV 
    | 显示表格

    文中基于LS-DYNA有限元软件, 分析了网格尺寸和一次项系数对中近场冲击波峰值压力和整体计算误差的影响, 主要得到如下结论:

    1) 在水下爆炸数值仿真中, 近场冲击波峰值压力受一次项系数和网格尺寸影响较大, 随λ的增大计算峰值压力对网格的敏感性降低, 且网格密度较大时, 过小的一次项系数会加大伪振荡, 致使计算峰值压力与经验公式值的相对误差增大;

    2) 通过研究网格尺寸及一次项人工粘性对冲击波峰值压力的影响, 获得了20%范围内的平均误差EP与网格密度因子λ和一次项系数Q2之间的关系, 并进一步拟合, 获得预定精度下快速确定网格尺寸和一次项人工粘性的预估模型为${E_{\rm{p}}} = 1\;413.57\times{\left( {{{{Q_2}^{1.895}} / \lambda }} \right)^{0.76}} + 4.06$

    3) 通过0.2~5 000 kg范围内的TNT柱形装药(长径比为1)和球形装药的水下爆炸计算, 验证了误差预估模型可适用于不同装药量的柱形(长径比为1)和球形装药二维中近场范围内的水下爆炸计算问题。

  • 图  1  数值计算模型

    Figure  1.  Numerical calculation model

    图  2  Q2=0.08时, 不同λPm随比例爆距变化曲线

    Figure  2.  Variation of Pm with scaled blast distance for different λ at Q2=0.08

    图  3  Q2=0.08时, 不同比例爆距处Pmλ变化曲线

    Figure  3.  Variation of Pm with λ for different scaled blast distances at Q2=0.08

    图  4  λPm上升速度影响

    Figure  4.  Effect of λ on the rate of increase of Pm

    图  5  λ=3时, 不同Q2Pm随比例爆距变化曲线

    Figure  5.  Variation of Pm with scaled blast distance for different Q2 at λ=3

    图  6  λ=3时, 不同比例爆距处PmQ2变化曲线

    Figure  6.  Variation of Pm with Q2 for different blast scaled distances at λ =3

    图  7  峰值压力平均误差

    Figure  7.  Average errors of peak pressure

    图  8  限定区间后峰值压力平均误差曲线

    Figure  8.  Average errors of peak pressure curves after limited interval

    图  9  lgEP和lgλ关系曲线

    Figure  9.  The relationship curves between lgEP and lgλ

    图  10  截距-对数粘性系数线性拟合

    Figure  10.  Linear fitting of intercept-logarithm viscosity coefficient

    图  11  峰值压力平均误差拟合结果

    Figure  11.  Fitting results of peak pressure average error

    表  1  TNT状态方程参数

    Table  1.   State equation parameters for TNT

    ${A / { {\text{GPa} }} }$${B / { {\text{GPa} }} }$$ {R_1} $$ {R_2} $$ \omega $$ {E / {{\text{GPa}}}} $
    371.203.234.150.950.307.00
    下载: 导出CSV

    表  2  水状态方程参数

    Table  2.   State equation parameters for water

    C/(m/s)$ {S_1} $$ {S_2} $$ {S_3} $$ \alpha $$E$$ {V_0} $
    14832.106−0.1740.011.0
    下载: 导出CSV

    表  3  空气状态方程参数

    Table  3.   State equation parameters for air

    E/(J/kg)C0C1C2C3C6$ {C_4} $$ {C_5} $
    2.5×106000000.40.4
    下载: 导出CSV

    表  4  不同爆距处峰值压力对比

    Table  4.   Comparison of peak pressure at different scaled blast distances MPa

    PmZ /(m/kg1/3)
    0.751.171.882.80
    试验值70.1039.6023.3015.30
    数值解75.5042.0022.7013.60
    经验值72.6043.8025.7016.40
    下载: 导出CSV

    表  5  λ=6时, 不同比例爆距处冲击Pm对比

    Table  5.   Different scaled blast distances peak pressure of shock waves at λ =6 MPa

    工况类型Z/(m/kg−1/3)
    0.3270.5000.7000.9001.1701.700
    Q2=0.020253.000138.00092.80067.50048.80030.700
    Q2=0.060235.000127.00081.80058.30041.50025.400
    经验值235.800124.70078.40059.00043.90028.800
    下载: 导出CSV

    表  6  lgEP关于lgλ的拟合参数

    Table  6.   Fitting parameters for λ in EP

    拟合参数Q2
    0.030 00.040 00.050 00.060 00.070 00.080 00.1 000
    k−0.690 0−0.790 0−0.590 0−0.510 0−0.530 0−0.530 0−0.490 0
    b1.300 01.420 01.410 01.410 01.500 01.550 01.610 0
    R1 0.940 40.967 60.963 30.985 80.996 10.999 30.999 9
    R20.884 30.936 30.927 90.971 70.992 20.998 60.999 8
    下载: 导出CSV

    表  7  b关于Q2的拟合参数

    Table  7.   Fitting parameters for b in Q2

    k1b1RR2
    0.586 02.180 00.971 60.944 1
    下载: 导出CSV

    表  8  不同工况下峰值压力误差

    Table  8.   Errors of peak pressure for different operating conditions

    装药质量/kg装药形状平均误差/%预估误差/%
    0.2柱形9.7510.36
    0.2球形10.79
    5.0柱形9.45
    5.0球形10.45
    50.0柱形10.83
    50.0球形10.30
    500.0柱形9.62
    500.0球形10.18
    1 500.0柱形9.57
    1 500.0球形9.89
    5 000.0柱形9.58
    5 000.0球形9.76
    下载: 导出CSV
  • [1] 蔡尚. 水下爆炸作用下舰船毁伤效能评估及水雷布阵策略优化研究[D]. 哈尔滨: 哈尔滨工程大学, 2018.
    [2] 田影. 不同边界条件下近场水下爆炸载荷特性研究[D]. 大连: 大连理工大学, 2022.
    [3] 辛春亮, 秦健, 刘科种, 等. 基于LS-DYNA软件的水下爆炸数值模拟研究[J]. 弹箭与制导学报, 2008, 28(3): 156-158.

    Xin Chunliang, Qin Jian, Liu Kezhong, et al. Research on UNDEX numerical simulation based on LS-DYNA[J]. Journal of Projectiles Rockets Missiles and Guidance, 2008, 28(3): 156-158.
    [4] Huang H, Jiao J Q, Nie X J, et al. Numerical modeling of underwater explosion by one-dimensional ANSYS-AUTODYN[J]. Journal of Energetic Materials, 2011, 29(4): 292-325. doi: 10.1080/07370652.2010.527898
    [5] Wang G, Wang Y, Lu W, et al. On the determination of the mesh size for numerical simulations of shock wave propagation in near field underwater explosion[J]. Applied Ocean Research, 2016, 59: 1-9. doi: 10.1016/j.apor.2016.05.011
    [6] 胡亮亮, 黄瑞源, 李世超, 等. 水下爆炸冲击波数值仿真研究[J]. 高压物理学报, 2020, 34(1): 102-114.

    Hu Liangliang, Huang Ruiyuan, Li Shichao, et al. Shock wave simulation of underwater explosion[J]. Chinese Journal of High Pressure Physics, 2020, 34(1): 102-114.
    [7] 张社荣, 李宏璧, 王高辉, 等. 水下爆炸冲击波数值模拟的网格尺寸确定方法[J]. 振动与冲击, 2015, 34(8): 93-100.

    Zhang Sherong, Li Hongbi, Wang Gaohui, et al. A method to determine mesh size in numerical simulation of shock wave of underwater explosion[J]. Journal of Vibration and Shock, 2015, 34(8): 93-100.
    [8] 杨坤, 陈朗, 伍俊英, 等. 计算网格与人工粘性系数对炸药水中爆炸数值模拟计算的影响分析[J]. 兵工学报, 2014, 35(S2): 237-243.

    Yang Kun, Chen Lang, Wu Junying, et al. The effects of computing grid and artificial viscosity coefficient on underwater explosion numerical simulation[J]. Acta Armamentarii, 2014, 35(S2): 237-243.
    [9] 闫秋实, 常松. 水下爆炸三维数值模拟特征参量敏感性分析[J]. 北京工业大学学报, 2023, 49(10): 1099-1108.

    Yan Qiushi, Chang Song. Underwater explosion 3D numerical simulation characteristic parameter sensitivity analysis[J]. Journal of Beijing University of Technology, 2023, 49(10): 1099-1108.
    [10] 马腾, 王金相, 刘亮涛, 等. 不同长径比柱形装药水下爆炸冲击波演化规律[J]. 振动与冲击, 2022, 41(8): 149-157, 222.

    Ma Teng, Wang Jinxiang, Liu Liangtao, et al. Shock wave evolution of cylindrical charge with different slender ratios[J]. Journal of Vibration and Shock, 2022, 41(8): 149-157, 222.
    [11] 孟龙, 黄瑞源, 王金相, 等. 小当量梯恩梯水下爆炸气泡脉动的数值模拟[J]. 兵工学报, 2020, 41(S1): 64-71.

    Meng Long, Huang Ruiyuan, Wang Jinxiang, et al. Numerical simulation of bubble pulsation of small scaled TNT in underwater explosion[J]. Acta Armamentarii, 2020, 41(S1): 64-71.
    [12] Huang C, Liu M, Wang B, et al. Underwater explosion of slender explosives: Directional effects of shockwaves and structure responses[J]. International Journal of Impact Engineering, 2019, 130(8): 266-280.
    [13] Wang J, Yang L, Xu Z, et al. Numerical simulation on underwater explosion in small-sized containers[J]. Mathematical Modelling of Engineering Problems, 2016, 3(3): 151-156. doi: 10.18280/mmep.030307
    [14] 高源. 炸药深水爆炸载荷研究[D]. 北京: 北京理工大学, 2021.
    [15] 辛春亮. 有限元分析常用材料参数手册[M]. 北京: 机械工业出版社, 2022.
    [16] Cole R H. Underwater explosions[M]. New York: Dover Publications, 1965.
    [17] Zamyshlyaev B V, Yakovlev Y S. Dynamic loads in underwater explosion[R]. Washington, D. C: Naval Intelligence Support Center, 1973.
  • 期刊类型引用(1)

    1. 张竞元,李海龙,冯海云,王龙侃,张之凡. 双发CL-20基含铝装药水下爆炸载荷特性. 火炸药学报. 2024(08): 747-757 . 百度学术

    其他类型引用(2)

  • 加载中
图(11) / 表(8)
计量
  • 文章访问数:  330
  • HTML全文浏览量:  137
  • PDF下载量:  30
  • 被引次数: 3
出版历程
  • 收稿日期:  2023-08-18
  • 修回日期:  2023-09-10
  • 录用日期:  2023-09-18
  • 网络出版日期:  2024-01-18

目录

/

返回文章
返回
服务号
订阅号