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

在用 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 柱形装药空气中爆炸冲击波峰值压力（超压 + 大气压力）

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

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

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

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