塑性力學(xué)和損傷力學(xué)是常用的描述材料非線性破壞的理論。損傷力學(xué)可以模擬材料剛度下降，塑性力學(xué)可以模擬材料塑性變形，將兩者結(jié)合成塑性-損傷力學(xué)，就可以較為完整的描述材料的非線性破壞過程。許多商業(yè)軟件，例如ABAQUS，自帶實(shí)體單元的塑性損傷力學(xué)，但是不具備cohesive單元的塑性損傷模型。coheisve單元結(jié)合實(shí)體單元，可以仿真各種材料的開裂過程。本文根據(jù)以下參考文獻(xiàn)，將文獻(xiàn)中用于DEM的本構(gòu)模型進(jìn)行修改，實(shí)用vumat子程序，使得cohesive單元具備塑性-損傷破壞的性質(zhì)。

Nguyen, N. H., Bui, H. H., Nguyen, G. D., & Kodikara, J. (2017). A cohesive damage-plasticity model for DEM and its application for numerical investigation of soft rock fracture properties.International Journal of Plasticity,98, 175-196.

首先介紹本構(gòu)模型中的參數(shù)和模擬時(shí)相應(yīng)的數(shù)值：

具體建模結(jié)果和更多詳細(xì)內(nèi)容見知乎文章：https://zhuanlan.zhihu.com/p/359076668

子程序vumat和input文件為收費(fèi)內(nèi)容, 在附件中。