基于AUTODYN模拟破片高速冲击引爆炸药

本文以一个简单案例介绍破片冲击炸药导致炸药引爆，其中破片采用Cu，炸药采用COMP-BJJ1，材料参数均取自AUTODYN自带数据库，Cu飞片额外定义一个材料失效，防止计算过程中因网格变形太大使得时间步太小导致计算终止，其中破片起始速度3000m/s。

破片采用拉格朗日算法，空气域及炸药采用欧拉算法，空气域建立完毕后将炸药填充到空气PART中。节省计算资源，采用1/2模型。为防止炸药爆炸后在壁面发生反射造成结果不准确，在除对称面的边界上均施加流出边界；对称轴上每隔一段距离设置一个固定高斯监测点，用于后处理观察监测点处的压力变化。由于点火增长模型只支持cm-g-us单位制，故本仿真单位制采用cm-g-us。

炸药在飞片刚撞击到后就被引爆，各时刻炸药的压力云图如图2-4所示：

1us                      5us                        10us

不同时刻B炸药的压力云图

为观察炸药的反应程度，在计算前输出变量选择ALPHA，各时刻炸药的反应度云图如下：

1us                      5us                        10us

不同时刻B炸药的反应分数

AUTODYN自身的后处理功能同样可以绘出不同时刻的速度、位移、压力等曲线，并且可以直接用到文献中，当然也可以将其数据导出，自己用专业绘图软件绘制，各监测点的压力-时间曲线如下所示（软件自带曲线还挺漂亮的）：