本文利用 cohesive element 模擬兩個粘接的物體受剪破壞，利用的本構模型是Coulomb模型，有限元實現算法是帶有返回功能的前歐拉法（Crisfield在1991出版的Non-linear finite element analysis of solid and structure中的第6章有提到），軟件是abaqus，子程序為UMAT。

這種有限元模型主要應用在在膠合破壞的預測，比如磚結構的水泥砂漿、纖維復合材料加固結構中復合材料的剝落、多層玻璃的脫膠以及夾層板的滑移等等。可以說，與膠有關的結構都可以試著利用cohesive單元模擬，所以開發cohesive單元的本構模型是很多領域的熱點，因為相比起固體單元，cohesive單元的發展相對較晚。

和tension cut-off那篇一樣，我用前歐拉法推導出的tangent stiffness和Lourenco（1997）用后歐拉法推導出的tangent stiffness是一樣的。

Coulomb的返回算法：

視頻結果詳見：https://zhuanlan.zhihu.com/p/113150354

模擬結果：

應力：

應變：

步驟1：part （2D 模型）

步驟2：material

步驟3：assembly

步驟4：step

步驟5：interaction

步驟6：load

步驟7：mesh

步驟8：結果