Numerical Simulation Accuracy Study of Underwater Explosion Shock Waves
-
摘要: 在水下爆炸数值仿真研究中, 网格尺寸和一次项人工粘性系数对冲击波峰值压力计算结果有较大影响。在预定计算精度条件下, 快速确定网格尺寸及人工粘性对数值计算意义重大。为此, 文中基于LS-DYNA有限元软件, 建立78 g三硝基甲苯(TNT)二维水下爆炸数值计算模型, 重点分析网格尺寸和一次项粘性系数对水下爆炸冲击波峰值压力和整体计算误差的影响规律。结果表明, 随着网格密度因子的增大, 计算峰值压力对网格的敏感性降低, 且网格密度较大时, 过小的一次项系数会导致计算峰值压力与经验公式值的相对误差增大。在此基础上获得20%范围内误差与网格尺寸、粘性系数之间的关系, 并构建出可用于快速确定网格尺寸和一次项人工粘性系数的误差预估模型, 通过0.2 ~5 000 kg范围内的TNT柱形装药(长径比为1)和球形装药的水下爆炸计算, 验证了预估模型的普适性, 可为二维中近场范围内的水下爆炸冲击波数值仿真计算研究提供参考。Abstract: In the numerical simulation study of the underwater explosion, the grid size and the artificial viscosity coefficient of the primary term have a large impact on the calculation results of the peak pressure of the shock wave. Under the condition of predetermined calculation accuracy, it is of great significance to quickly determine the grid size and artificial viscosity for numerical calculation. For this reason, based on LS-DYNA finite element software, a two-dimensional underwater explosion numerical calculation model of 78 g trinitrotoluene (TNT) was established to analyze the influence of the grid size and the viscosity coefficient of the primary term on the peak pressure of the underwater explosion shock wave and the overall calculation error. The results show that with the increase in the grid density factor, the sensitivity of calculated peak pressure to the grid decreases. When the grid density is larger, a small primary term coefficient will cause the relative error between the calculated peak pressure and the empirical formula value to increase. On this basis, the relationship among the error, grid size, and viscosity coefficient within 20% is obtained, and an error prediction model that can be used to quickly determine the grid size and the artificial viscosity coefficient of the primary term is constructed. Through the underwater explosion calculation of cylindrical TNT charge (aspect ratio of 1) and spherical TNTcharge in the range of 0.2–5 000 kg, the universality of the prediction model is verified, which can provide a reference for the numerical simulation of underwater explosion shock wave in the two-dimensional near-field range.
-
0. 引言
水下武器作为舰船生命力的主要威胁之一, 其爆炸冲击波及气泡载荷会对舰船造成严重的局部和总体破坏[1]。各国相继开展了大量水下爆炸试验, 但实弹实测试验安全风险较高、成本高昂、重复性低且观测范围有限[2]。随着计算机性能和仿真技术的发展, 数值仿真以其较高计算精度、低成本和高可重复性等优点成为继实验、理论研究后第3种水下爆炸研究方法[3]。因此开展水下爆炸数值仿真精度研究具有十分重要的意义。
数值仿真研究水下爆炸冲击波载荷问题时, 为了处理冲击波的强间断面, 抑制波阵面前后的数值振荡, 引入了人工粘性。由于人工粘性的引入会在几个网格宽度上光滑冲击波波阵面, 因此水下爆炸冲击波载荷计算结果直接受人工粘性的影响。同时,人工粘性的引入要求网格尺寸不能过大, 否则计算过程中网格会忽略部分波阵面信息, 致使峰值过低影响计算精度。Huang等[4]通过典型三硝基甲苯(Trinitrotolue,TNT)炸药水下爆炸数值分析, 探讨了一次、二次粘性系数对数值计算结果的影响, 并给出了一定比例爆距范围内的建议值。Wang等[5]研究了不同装药质量下网格尺寸对冲击波峰值的影响, 并引入与装药半径和单元边长相关的无因次量表征网格尺寸。胡亮亮等[6]以常规TNT水下爆炸为例, 对水的状态方程、人工粘性系数和网格尺寸对于数值仿真结果的影响进行了研究。张社荣等[7]基于有限元软件AUTODYN建立了不同炸药量的水下爆炸数值模型, 对比分析了网格尺寸对不同爆距处冲击波峰值压力及比冲量的影响。此外还有其他学者[8-13]讨论了网格尺寸及粘性系数对水下爆炸计算结果的影响, 但现有研究无法在预定精度下快速确定网格尺寸和人工粘性系数。因此文中以TNT水下自由场爆炸数值计算为例, 探究网格尺寸和一次项人工粘性系数对水下爆炸冲击波峰值压力的影响, 综合分析不同工况下网格密度因子和一次项系数与峰值压力平均误差间的关系, 构建出普适性较高的水下爆炸数值误差预估模型, 为预定精度的仿真模型设计提供依据。
1. 水下爆炸数值仿真
1.1 数值模型
基于文献[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.2 数值模型
装药选用典型的单质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为相对体积; A、B、R1、R2为常数, 取值如表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.20 3.23 4.15 0.95 0.30 7.00 数值仿真中常用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为水中声速; S1、S2和S3为常数; γ0为GRUNEISEN系数; α为体积修正系数; V0为初始相对体积。以上参数取值如表2所示。
表 2 水状态方程参数Table 2. State equation parameters for waterC/(m/s) $ {S_1} $ $ {S_2} $ $ {S_3} $ $ \alpha $ $E$ $ {V_0} $ 1483 2.106 −0.174 0.01 − − 1.0 空气使用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) 式中: C0、C1、C2、C3、C4、C5和C6为常数。取值如表3所示。
表 3 空气状态方程参数Table 3. State equation parameters for airE/(J/kg) C0 C1 C2 C3 C6 $ {C_4} $ $ {C_5} $ 2.5×106 0 0 0 0 0 0.4 0.4 1.3 冲击波峰值压力计算结果
根据文献[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 distancesMPa Pm Z /(m/kg1/3) 0.75 1.17 1.88 2.80 试验值 70.10 39.60 23.30 15.30 数值解 75.50 42.00 22.70 13.60 经验值 72.60 43.80 25.70 16.40 对比表4数据, 可以发现数值解与试验值最小误差仅2.6%, 平均误差6.9%; 与经验值最小误差4%, 平均误差9.2%。因此网格尺寸和一次项粘性系数取值合理时, 数值模型可获得较高的计算精度。经验值与试验值峰值压力平均误差7.2%, 表明经验公式可以较为准确地预估水下爆炸冲击波峰值压力。因此, 在探究网格尺寸及人工粘性对数值计算精度的影响时, 仅比较数值解与经验值。
2. 网格尺寸及人工粘性对Pm的影响
在有限元计算中为了反映冲击波波阵面的强间断, 引入了人工粘性来光滑冲击波, 这使得数值计算峰值压力低于真实值。一次项人工粘性系数对冲击波峰值压力影响较大, 二次项人工粘性系数主要用于抑制冲击波衰减过程中的虚假振荡, 但对水中爆炸数值仿真计算中冲击波峰值影响较小。参阅文献[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个工况。为了从全局意义上分析网格尺寸和Q2对Pm的影响, 引入平均误差, 定义为比例爆距0.375~2.8 m/kg1/3范围内选定观测点峰值压力相对于经验公式误差的平均值。
网格尺寸及Q2对Pm结果如图2~图6所示。由图2和图3可知, 当λ较小时, Pm较低, 与经验公式值偏差较大, 随着λ的增加, Pm显著增大且近场峰值压力与经验值吻合较好。但随着λ增加, Pm趋于稳定, 不会随着λ的变化出现显著变化, 这表明随着网格密度的增加, Pm对网格的敏感性越来越低。与此同时, 网格尺寸对冲击波超压爬升至峰值所需时间影响很大, λ越小, 冲击波超压爬升至峰值时间越长且峰值压力越低, 如图4所示。
由图5和图6所见, 随着Q2逐渐变小, Pm逐渐增大, 且近场冲击波峰值压力受Q2影响较大。由于Q2取值范围较小, 不同比例爆距处Pm随Q2的变化曲线与文献[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 λ =6MPa 工况类型 Z/(m/kg−1/3) 0.327 0.500 0.700 0.900 1.170 1.700 Q2=0.020 253.000 138.000 92.800 67.500 48.800 30.700 Q2=0.060 235.000 127.000 81.800 58.300 41.500 25.400 经验值 235.800 124.700 78.400 59.000 43.900 28.800 对比数值解和经验值在不同爆距处的冲击波峰值压力, 获得不同Q2下峰值压力平均误差EP随λ的变化关系, 如图7所示。整体来看, 随着λ的增大和Q2的减小, EP逐渐降低。但以Q2=0.10为例, 随着λ的增大, EP先减小后增大, 这是因为网格密度较大时, 过小的Q2会加大伪振荡, 致使Pm数值解与经验的相对误差增大, 因此在数值计算中不能一味地减小网格尺寸和一次项系数。同时, 当λ=1时, EP很大, 即使调整Q2也未能使EP满足工程精度。因此在水下爆炸数值计算中应首先确定λ, 同时调整Q2方能得到精度较高的计算结果。
3. 计算误差预估模型
为了便于应用, 实现在预定精度下快速确定网格尺寸和一次项系数, 需构建出关于λ和Q2的误差EP预估模型。考虑到实际工程需要, 误差应控制在20%以内。同时, 由前面分析可知, 较大的网格密度以及过小的一次项系数均可能导致计算误差增大。因此, 将误差预估模型的变量区间范围限定在0.03≤Q2≤0.10和3≤λ≤8, 则图7中的数据在该区间范围显示为图8中曲线。
将图8中的EP和λ取对数得到lgEP和lgλ的关系曲线如图9所示, k为斜率, b为截距。可以看到, 各曲线呈近似平行的线性关系, 因此可构造线性表达式为
$$ \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 0 0.040 0 0.050 0 0.060 0 0.070 0 0.080 0 0.1 000 k −0.690 0 −0.790 0 −0.590 0 −0.510 0 −0.530 0 −0.530 0 −0.490 0 b 1.300 0 1.420 0 1.410 0 1.410 0 1.500 0 1.550 0 1.610 0 R1 0.940 4 0.967 6 0.963 3 0.985 8 0.996 1 0.999 3 0.999 9 R2 0.884 3 0.936 3 0.927 9 0.971 7 0.992 2 0.998 6 0.999 8 $$ b = {k_1}\lg {Q_2} + {b_1} $$ (9) 利用式(9)对图10中数据进行线性拟合, 结果如表7所示, 可以看出数据具有较高拟合精度。将式(9)代入式(8), 整理得
表 7 b关于Q2的拟合参数Table 7. Fitting parameters for b in Q2k1 b1 R R2 0.586 0 2.180 0 0.971 6 0.944 1 $$ {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)为最终得到的误差预估模型$$ {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.75 10.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 4. 结论
文中基于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 TNT状态方程参数
Table 1. State equation parameters for TNT
${A / { {\text{GPa} }} }$ ${B / { {\text{GPa} }} }$ $ {R_1} $ $ {R_2} $ $ \omega $ $ {E / {{\text{GPa}}}} $ 371.20 3.23 4.15 0.95 0.30 7.00 表 2 水状态方程参数
Table 2. State equation parameters for water
C/(m/s) $ {S_1} $ $ {S_2} $ $ {S_3} $ $ \alpha $ $E$ $ {V_0} $ 1483 2.106 −0.174 0.01 − − 1.0 表 3 空气状态方程参数
Table 3. State equation parameters for air
E/(J/kg) C0 C1 C2 C3 C6 $ {C_4} $ $ {C_5} $ 2.5×106 0 0 0 0 0 0.4 0.4 表 4 不同爆距处峰值压力对比
Table 4. Comparison of peak pressure at different scaled blast distances
MPa Pm Z /(m/kg1/3) 0.75 1.17 1.88 2.80 试验值 70.10 39.60 23.30 15.30 数值解 75.50 42.00 22.70 13.60 经验值 72.60 43.80 25.70 16.40 表 5 λ=6时, 不同比例爆距处冲击Pm对比
Table 5. Different scaled blast distances peak pressure of shock waves at λ =6
MPa 工况类型 Z/(m/kg−1/3) 0.327 0.500 0.700 0.900 1.170 1.700 Q2=0.020 253.000 138.000 92.800 67.500 48.800 30.700 Q2=0.060 235.000 127.000 81.800 58.300 41.500 25.400 经验值 235.800 124.700 78.400 59.000 43.900 28.800 表 6 lgEP关于lgλ的拟合参数
Table 6. Fitting parameters for λ in EP
拟合参数 Q2 0.030 0 0.040 0 0.050 0 0.060 0 0.070 0 0.080 0 0.1 000 k −0.690 0 −0.790 0 −0.590 0 −0.510 0 −0.530 0 −0.530 0 −0.490 0 b 1.300 0 1.420 0 1.410 0 1.410 0 1.500 0 1.550 0 1.610 0 R1 0.940 4 0.967 6 0.963 3 0.985 8 0.996 1 0.999 3 0.999 9 R2 0.884 3 0.936 3 0.927 9 0.971 7 0.992 2 0.998 6 0.999 8 表 7 b关于Q2的拟合参数
Table 7. Fitting parameters for b in Q2
k1 b1 R R2 0.586 0 2.180 0 0.971 6 0.944 1 表 8 不同工况下峰值压力误差
Table 8. Errors of peak pressure for different operating conditions
装药质量/kg 装药形状 平均误差/% 预估误差/% 0.2 柱形 9.75 10.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 -
[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)
-