本文介紹了如何用后歐拉算法模擬2D cohesive單元的復雜破壞，采用的本構模型是hyperbolic曲線。這種曲線可以將受拉破壞和受剪破壞耦合在一起，很好地反應了膠合接觸面地特性。

主要參考文獻為：Caballero, A., Willam, K.J. and Carol, I., 2008. Consistent tangent formulation for 3D interface modeling of cracking/fracture in quasi-brittle materials. Computer Methods in Applied Mechanics and Engineering, 197(33-40), pp.2804-2822.

Willam, K.J 和 Carol, I. 長期致力于模擬巖石、混凝土和磚結構地開裂破壞。

本文的很多數學推導是通過軟件Wolfram Mathematica獲得的。

本文所用的算法是Caballero（2008）中的算法1：monolithic iteration strategy without substepping。

具體建模過程和建模結果見知乎文章：https://zhuanlan.zhihu.com/p/113538156

umat子程序和代碼對應的

詳細解答見附件。