混凝土在外載荷作用下的非線性行為中同時包含微裂縫和塑性流動這兩種微觀機制的影響。在考慮混凝士等準脆性材料的非彈性力學行為方面，連續損傷力學模型可以通過不同的方式來描述材料剛度和強度的退化以及單邊效應。真正意義上的彈塑性損傷本構模型：不僅考慮卸載時不可恢復塑性變形的影響，而且還應該考慮損傷和塑性的雙向耦合效應。

彈性階段應力應變滿足如下關系

通過對應力進行譜分解，可得

式中，σ 為名義應力，d 為損傷，d=1-exp(-εp/ρ0)，公式右端σ為有效應力的正負分解。

拉壓屈服函數如下所示

屈服后，塑性流動由下式定義

按照彈性預測-塑性修正-損傷修正的流程，通過在主應力空間進行譜分解，結合徑向返回算法，本文編寫了混凝土彈塑性損傷的VUMAT子程序。

通過對單胞的單向拉壓模擬可以計算得到混凝土的應力應變響應如下圖所示。

拉伸損傷演化過程

壓縮損傷演化過程

不同圍壓下的應力應變曲線