显式动力学侵彻笔记

*ERROR: ELEMENT 12345 HAS NEGATIVE VOLUME 弹出来侵彻计算基本废了，网格畸变太严重，单元被压到负体积，求解器直接终止。先看哪个部件，靶板的话通常是网格太粗或材料失效参数没给对。JC本构D1到D5，D1初始失效应变，D2指数，D3应力三轴度，D4应变率，D5温度，不设单元可以一直压下去不删，删了又质量不守恒。D1给0.05到0.1，D2给0.5到1.0。

侵彻网格尺寸比静力学敏感得多，靶板撞击区0.2mm以下，弹体0.5mm，接触面附近细化。六面体优先，四面体刚度偏大。ALE缓解畸变但数值扩散，弹速降得比实验快。SPH适合超大变形，边界处理不好。LS-DYNA在侵彻比Abaqus成熟，材料库和接触算法多，Autodyn某些军用工况有现成参数。

去年做钨合金弹打装甲钢，网格0.15mm，沙漏控制默认粘性系数，算到一半沙漏能飙到内能18%，弹体像橡皮弹回来。那天下午空调坏了，办公室三十四度，旁边同事吃瓜子咔咔响，我盯着能量曲线心里发毛，后来改刚度形式，系数从0.1降到0.03，加单元删除，压到5%以内。温度一高脑子慢，参数改错三次才发现。

接触用 *CONTACT_ERODING_SURFACE_TO_SURFACE 或Abaqus通用接触，摩擦系数0.1到0.3，法向刚度不能太大否则接触力振荡。侵蚀接触必须开，单元删了新表面继续接触。弹体网格细做主面，靶板做从面。初始速度给弹体，靶板固定边界离撞击区足够远，否则应力波反射影响结果。

显式时间步由最小单元尺寸和材料波速决定，钢波速5000m/s，0.2mm单元步长4e-8s，100微秒要250万步，算一天正常。质量缩放放大时间步，因子超2%总质量改变惯性，弹速偏低。固定时间步缩放目标1e-7，监控总质量增加不超过1%。并行用8核以上，domain decomposition别跨接触面。

能量平衡判断结果可信度。总能量守恒，动能转内能，沙漏能低于5%，伪应变能低于2%，有单元删除看删除带走的能量。历史输出ALLAE沙漏能，ALLIE内能，ALLKE动能，ETOTAL接近零。弹体剩余速度跟实验差10%以内就算不错。

重启文件别设太密，侵彻算例几十G，每50微秒写一次够了。场输出每1微秒，历史输出每0.1微秒，不然结果文件巨大。双精度更稳，单精度累积误差让弹体轨迹偏。GPU加速LS-DYNA能到3倍但只对特定求解器，Abaqus支持有限。

材料参数别直接用文献的，应变率效应占主导，JC模型A、B、n、C、m根据实验曲线拟合，C通常0.01到0.05。温度软化不能忽略，侵彻界面上千度。复合材料用Hashin加连续损伤，显式收敛困难，粘性正则化。算之前先跑低分辨率验证接触和边界，再上细网格。

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