一种模拟节理岩体破坏的新方法:单元劈裂法

在有限单元法中,三结点三角单元在几何上具有独特的优点,即当一条裂纹贯穿单元(单元劈裂)时,总有一结点位于裂纹的一侧,而另外两个结点位于裂纹的另一侧。位于一侧的两个结点与位于另外一侧的一个结点可以潜在地构成两个接触点对,从而可以利用这两个接触点对来推导劈裂单元的刚度矩阵,以再现裂纹面之间的接触和摩擦效应。利用该三角单元的几何特点,提出了单元劈裂法,用于模拟节理或其它裂纹的扩展问题。由于劈裂单元刚度矩阵与原三角单元享有共同的结点,因而在模拟裂纹扩展过程中无需为设置节理单元而改变原网格划分方案,这为节理扩展的数值模拟提供了很大方便,提高了计算效率。在劈裂单元刚度矩阵推导过程中,没有考虑单元劈裂后所形成两个块体自身的弹塑性变形,所以该方法只是一种近似方法。对相应试验进行了数值模拟,结果表明该方法是有效的。