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

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

這里有意思的是，我用前歐拉法推導出的tangent stiffness和Lourenco（1997）用后歐拉法推導出的tangent stiffness是一樣的。

tensile cut-off的返回算法：

模擬的視頻結果詳見：

https://zhuanlan.zhihu.com/p/113143055

模擬結果：

mises應力圖

ABAQUS建模：這是一個二維的單元，兩個彈性材料物體被一層cohesive單元（厚度為10）連接著。面外厚度我定位100（一般磚塊大約厚100mm）

步驟1：part

步驟2：material

步驟3：assembly

步驟4：step

步驟5：interaction

步驟6：load

步驟7：mesh