实现低围压应变软化,高围压更接近理想塑性的煤的本构模型(用于分析图3),用连续损伤力学的热力学框架,采用隐式返回欧拉积分算法计算。要求:以论文1的框架为主,采用塑性和损伤耦合的修正方法,其中,屈服准则和屈服势函数采用图2论文中的模型;损伤准则和损伤势函数采用图1论文中的公式。或者主要复现论文1,在此基础上略加修改也行