金属导体在强脉冲电流的作用下会因为能量的迅速积累而发生快速相变，转化为高温等离子体，进而出现爆炸现象。利用这种现象，军工科研人员研制了爆炸丝雷管，爆炸箔起爆器等，用于直列式起爆序列等应用场景。同时，电爆炸驱动飞片也常作为材料高压特性研究中的脉冲高压源。本报告基于爆炸箔驱动飞片过程进行仿真研究，以探究电爆炸驱动过程的特点。

       首先以南京理工大学曹始发的硕士论文中给出的尺寸作为我们算例中的结构尺寸，其中飞片材料为聚酰亚胺，厚度25微米。爆炸箔材质为铜，厚度5微米，长度和宽度均为0.3毫米。加速膛材质为蓝宝石玻璃，尺寸为Φ0.45毫米×0.40毫米。基板（底板）的材质为三氧化二铝陶瓷。仿真软件使用ANSY/LS-DYNA，采用ALE耦合方法建立二维1/2模型，其中的空气域和铜等离子体采用Euler算法，基板、加速膛和飞片采用Lagrange算法，选择爆炸力学数值模拟中常用的cm/g/μs单位制。仿真模型图如下图所示：

       考虑空气域，一共需要五种材料模型参数：三氧化二铝陶瓷、蓝宝石玻璃、聚酰亚胺、空气，和铜等离子体。其中除铜等离子体外其他四种均可以通过相关文献直接查到，在此不作赘述。而对于铜等离子，考虑到其在P-V膨胀做功关系上应该与γ律状态方程具有相似之处，使用LS-DYNA中的空白材料模型与线性多项式状态方程描述其本构关系及膨胀做功过程。

       其中，线性多项式状态方程的形式如（1）所示：

P=C₀+C₁μ+C₂μ²+C₃μ³+(C₄+C₅μ+C₆μ²)E

μ=(ρ/ρ₀)-1

当C₀=C₁=C₂=C₃=C₆=0，且C₄=C₅=γ-1时，此时该状态方程等效为γ律状态方程，不过要注意的是此时气体内能E的单位不再是焦耳，而应当与压力的单位相同。

       可以看到，要建立铜等离子材料模型，主要需要确定其初始密度ρ₀、等熵指数γ、内能E三种参数。对于ρ₀，考虑到铜箔在脉冲大电流作用下所发生的从固态到等离子体态之间的转变可视为瞬态过程，因此将ρ₀值设置为8.9g/cm³。对于γ，我们取其值为2.0。对于能量E，其值应当与输入电能有关，将沉积能量设置为0.5J，则可以算出E值：

E=(5×10^-6/(0.03×0.03×0.0005))/8.9≈1.25

就这样，我们得到了需要的参数，输入k文件中，然后带入ANSYS/LS-DYNA求解器进行求解，得到了仿真结果，如下所示：

从数值模拟结果的压力云图可以看到，在铜等离子体高压作用下聚酰亚胺被迅速剪切形成飞片，我们进一步导出飞片的速度，如下所示：

从飞片速度曲线中也可以看到飞片的速度在初始阶段迅速升高，然后逐渐趋于稳定。这也与采用PDV、VISAR等实验手段测得的飞片速度变化规律一致。

       不过需要注意的是，本报告所做的数值模拟并不能够反映爆炸箔初期在脉冲大电流作用下的迅速相变过程，这也是下一步需要仔细考虑和改进的地方。