对于第 1 点中，在骨试件测试二次开发的本构，带失效，和状态方程，没有问题，说明本构、失效、状态方程都没有问题。

第 3 点中，使用 MAT22 本构时没有问题，指的是在第 2 点中，模拟弹丸撞击靶板时，使用 MAT22 不会出现错误吗？另外，你说仅仅使用本构和状态方程，模拟也不会出问题，但是加了失效之后，就出现了，于是把问题定在状态方程有问题，这个理由不太充分，按照排除法，应该是失效存在问题。

说说第 2 点中，模拟弹丸撞击靶板，弹丸和靶板还没有接触的情况下，全部单元就都出现失效的问题。这个问题主要表现在仿真模拟出现了错误，错误的原因可能是材料模型的问题，或者是数值算法的问题，同时，也可能是数值建模的问题。

如果说，你有用 MAT22 本构模型，模拟了弹丸撞击靶板的过程，且没有报错，似乎可以证明数值建模没有问题，这个可以将数值建模是否有问题作为次要因素，但这也不能完全确定模型没有问题，例如，侵彻碰撞问题的计算，对网格的质量要求较高，网格的粗密和形状对计算结果有很大的影响。

然后是数值算法的问题，模拟弹丸对靶板的撞击，应该重点在定义侵蚀接触上。

最后，是题主比较关心的材料本构中的失效部分，在靶板的变形中，失效参数是一个重要的影响因素，但根据目前的情况看，尚不足以判断本构中的材料失效参数有问题。

所以，还是建议先排查一下数值建模、接触算法是否有问题，或者在问题详情中补充相关部分。

您好，请问您子程序中的单元失效删除是使用逻辑判断 failel=.true. 么，我想向我的子程序中引入单元删除，设置了最大应变的判据并用了这个逻辑判断 failel=.true.，但是单元应变超过了我设置的阈值仍不删除

材料参数可能有问题

你好，请问如果需要对状态方程进行二次开发，材料库是不是需要包含 dyn21b.f？

我怀疑是你的子程序写得有问题。我也遇上类似的问题。当二次开发采用应力为主的失效准则时，不至于第一步就大量单元被删除，但是侵彻过程中单元的删除明显多于弹孔，且与实验结果严重不符；当采用应变为主的失效准则时就会发生第一步删除靶板的所有单元的问题。