因為要仿真混凝土破壞實驗，考慮用abaqus里面的CDP模型，查閱了相關資料進行了理論總結，并根據理論編寫計算程序。

ABAQUS中CDP 模型中采用的是混凝土在單軸受力狀態下的應力和非彈性應變，非彈性應變根據混凝土的單軸應力-應變曲線換算。

根據GB50010-2010混凝土結構設計規范，混凝土單軸應力應變關系如圖：

根據GB50010-2010混凝土結構設計規范，混凝土單軸受壓應力應變曲線可按下列公式確定：

混凝土單軸受拉應力應變曲線可按下列公式確定：

下面給出詳細步驟和程序代碼。