2.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) 如采用二维轴对称的算法,药量是否应用一半计算; 




请先 登录 后评论

1 个回答

CSUA - 研究生
擅长:LS-DYNA

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

请先 登录 后评论
  • 1 关注
  • 0 收藏,366 浏览
  • 匿名 提出于 2017-06-30 20:52