5 计算爆炸空气冲击波的传播,出现一些问题

使用 AUTODYN5.0 计算爆炸空气冲击波,遇到一些问题。我的计算思路是:

attachments-2017-06-UoNcsqce5956481d49c9f.png 在用 AUTODYN 软件进行数值计算时,应考虑将三个位置用“空白”空间进行填充,从而确定该位置处的冲击波压力(反射作用)。100g TNT 柱形装药在空气中爆炸,设初始 TNT 柱形药柱半径为 r,高为 h =2r。则,可以求出,在 AUTODYN 中用 2 维轴对称算法,应该用一半的药量进行计算。

传播空域和网格划分

在 X(I) 方向划分成 800 个网格,每个网格 7mm;

 在 Y(J) 方向划分成 600 个网格,每个网格 10mm;

100g TNT 装药的网格划分:I:1~4;J:121~125;

爆炸点的位置:X=0;Y=1200mm;

测量点的位置:

Gauges 1:1m=1000mm,对应位置 (I,J)=(144,1);

Gauges 2:3m=3000mm,对应位置(I,J)=(431,1);

Gauges 3:5m=5000mm,对应位置 (I,J)=(715,1);

100g TNT 柱形装药空气中爆炸冲击波峰值压力(超压 + 大气压力)

attachments-2017-06-9yv0OLOD5956495473737.png

结果相差太大,不知问题出在何处,难道是单位制的问题?

1)在该软件中,单位制该如何选用:“默认是 mm-ms-mg,但爆炸类问题常用 cm-g-s”,

该如何确定单位制,我觉得这和软件内其他参量默认值有关;

2) 如采用二维轴对称的算法,药量是否应用一半计算;

  • 1292
请先 登录 后评论

1 个回答

CSUA - 研究生
擅长:LS-DYNA

材料模型和状态方程得参数决定了你的单位制系统。
药量不能按一半算。
可能是你的单位制有问题,刚好差了近三个量级

请先 登录 后评论