LS-DYNA案例：用SPH方法模拟弹片高速冲击毁伤效果

弹片速度超过1km/s之后拉格朗日网格会严重畸变，删单元删到后面剩余质量都不对了，这种工况直接换SPH粒子法做，AUTODYN或者LS-DYNA都行，LS-DYNA的FEM-SPH自适应转换更省事，靶板主体用拉格朗日，接触区转粒子，碎云形状和碎片边界都能拿到。

建模的时候弹片用SPH填满，粒子间距取弹片最小特征尺寸的四分之一到五分之一，间距再小算不动，粒子数上百万之后单步时间太长。靶板靠近弹道的区域加密，远离的区域可以粗一点，SPH本身对粒子间距突变敏感，过渡区做一层渐变，不要一刀切。

材料模型，钢弹片用Johnson-Cook本构配Gruneisen状态方程，参数库里有现成的，铝靶也是一样，参数直接调用，别忘了把失效应变和侵蚀准则打开，不然穿透过程出不来碎云。陶瓷靶的话JH-2模型更合适，本构和EOS要配对，别混着用。

边界条件，靶板四周做成无反射边界，模拟无限大靶，不然应力波从边界反射回来会干扰毁伤区的判断。如果想看真实板件的边界效应，就按实际尺寸建，四边固支，计算量会大不少。

接触用automatic_nodes_to_surfaces配SPH，LS-DYNA里SPH粒子跟拉格朗日单元的耦合用/contact_automatic_nodes_to_surface，主面给靶板，从节点给SPH粒子组，接触刚度默认值偏软，穿透深度对不上的时候调一下slide参数。

求解时间步长由最小粒子间距和声速决定，SPH的步长比拉格朗日小很多，算几十微秒物理过程跑几个小时到一天是正常的，MPP版本开多核，不然等不起。输出间隔设密一点，碎云演化的中间帧对分析二次毁伤范围有用。

验证这块，拿Piekutowski的经典铝球撞铝板算例对一遍，文献里碎片云形态、弹坑直径、弹道极限速度都有实测数据，模型调到跟实验对得上再算自己的工况，直接上弹片打复合靶板，结果对不对都说不清。后处理里统计碎片云的空间分布和粒子速度，用MATLAB把粒子数据读出来做直方图，能量转化关系能看得比较清楚，弹体损失的能量大头变成破片动能，这个规律反复都能出来。

粒子间距做一次敏感性分析，间距减半算一遍，碎云轮廓和弹坑直径变化不大就可以认为收敛，变化还大就继续加密，SPH的粒子数敏感性比拉格朗日的网格敏感性更容易被忽略。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删