abaqus单元删除的一般方法

如果问单元失效模拟首先就应该了解它的功能与目的，所以简单地介绍一下吧。
单元删除功能本身是为了克服有限元本身的缺陷的一项方法，由于有限元本身就是基于连续介质力学的，而在连续介质理学中，所研究的物体需要是连续的，既物质 域在空间中连续。在这样的理论假设框架下，单元本身是不会消失的。然而在实际情况下，由于损伤断裂的存在，势必会使得一些单元消失或者完全的失效，所以为 了能够模拟这种情况，abaqus提供了单元失效功能。
我所知道的单元失效情况一般可以概括为三种情况
1）单元损伤失效，这种单元失效可以用来模拟材料由于损伤，或其他原因导致刚度减小的情况。
2）单元直接删除技术，这种技术可以用来模拟基坑，隧道开挖而导致的材料消失情况。
3）vumat，这种方法本质上讲与第一种相类似，但是它可以根据用户自己的情况来删除单元，但是属于很高级的操作，难度也较大
当然，cohesive也可以算作一种单元删除的方法，只是本身很少做断裂模拟，对于那方面不是很了解。
在这三种方法中，我想讲一讲前两种情况。
1首先来讲一下第二种方法。
这种方法很简单
他的使用方法就是在edit keyword中加一行关键字
*MODEL CHANGE, TYPE=ELEMENT, REMOVE
Set-1
Set-1为所定义的单元集合，这一集合需要事先定义，而这一关键字的位置一般位于step后，如果用回想要在哪一分析部删除单元，就在哪一分析部后添加。
但是使用这一方法是要注意的是由于这个单元是一下子从模型中删除，所以在standard下可能造成模型总体刚度奇异，导致计算不收敛，所以最好的方法就 是能够在单元删除前，事先削减它的刚度，这一做法也很简单，就是利用温度变量作为一个独立的状态变量然后通过状态变量与单元刚度的对应关系，较小弹性模 量。
2现在来谈一谈使用第一种方法来删除单元的过程；
这一方法使大家最为关心的。后来想到原因可能是大家对于这一部分的理论与实际联系不上，所以还是有很多人不太明白。
其实我也是初学，只不过是对这方面比较有兴趣，所以和大家一起探讨一下吧，也算是抛砖引玉。
这种方法首先只能在