塑性損傷模型 Python，的案例
混凝土塑性損傷CDP模型的幾個問題 附2010規范用C50混凝土損傷塑性本構關系數據下載
作者技術鄰用戶名：星辰_北極星，2012年開始從事Abaqus仿真相關工作，服務大小課題逾百項; 主要仿真領域：石油工程、巖土工程和金屬加工工藝；重點研究方向：ABAQUS GUI二次開發、固體力學、斷裂以及損傷等。
如有項目咨詢可點擊閱讀原文聯系作者。
在以往的課程和技文中都曾描述過混凝土塑性損傷（CDP）模型，但由于描述不夠完整、清晰，還是給讀者和學員留下不少的疑問，在這里表示歉意。我們先看看經常會被問到的幾個問題：
單調荷載下，損傷定義是否有影響？
輸出的單元應力應變曲線為什么和輸入的不一樣？
單元應力為什么比屈服強度還高？
單元應力超過定義的最大屈服應力后的發展趨勢是怎樣的？為什么會出現應力增大情況？
混凝土輸入的是真實應力應變曲線還是名義應力應變曲線？
為了解決上面的這些問題，這里準備了一個非常簡單的模型，一個平面的四邊形單元，A點進行完全固定，B點約束Y方向位移，CD點施加-Y方向位移載荷，對單元形成單軸壓縮分析；右側的模型是對單元四個面額外施加P的圍壓。
后面幾幅圖是基于上面兩個模型，輸入不同參數條件獲得的，現在我們就來看圖說話。
【圖一】
單軸壓縮模型（左側）、材料模型CDP定義塑性、但沒有定義損傷參數。
結論：
當模型選擇平面應力單元時，得到的應力-應變曲線與輸入的材料一致；
當模型選擇平面應變單元時，應力-應變曲線各處均大于輸入的材料性質，這是由于平面應變單元存在平面外法向約束，因此并不是純粹的單軸壓縮模型；
輸入定義的末尾應變為0.37，后續的應力-應變曲線是水平延伸的，單軸條件下，不存在增大情況。
展開 基于塑性損傷模型（CDP）FRP約束混凝土ABAQUS有限元模型 ￥12.99
本模型為基于CDP的FRP約束混凝土ABAQUS有限元模型
1. 在部件的建立上，使用殼體模擬FRP，實體模擬混凝土
2. 在材料屬性上，混凝土采用CDP模型，基于混規。FRP材料的單層板模型，并且采用常規殼方式進行鋪層，自定義了“離散”坐標系。
3. 在分析部上，打開幾何非線性，輸出參考點RP-1的力和位移。
4. 在相互作用上，將加載力的平面耦合到參考點RP-1上，并將FRP與混凝土進行綁定
5. 在荷載上，對混凝土底端進行完全固定，限制上表面除了U3方向其他方向的位移。給予U3方向一定位移，采用位移加載。
6. 在網格部分，混凝土采用C3D8R，FRP采用S4R。
得到模型后，可以根據FRP層數、材料屬性進行修改，根據混凝土實際強度進行修改，輸出應力應變曲線或者其他需要的部分即可
以下為模型的CAE文件：
展開 ABAQUS UMAT - 混凝土塑性損傷模型的實現 ￥1500
混凝土塑性損傷模型在工程上應用較為廣泛，同類型的本構模型多內置于各類仿真軟件中，供用戶模擬混凝土結構的破壞和受力情況。本文根據Peter Grassl 和 Milan Jira′sek 2006年的文章《Damage-plastic model for concrete failure》進行本構模型代碼復現，并對文中的模型進行了一些簡化。
UMAT代碼和INPUT文件見付費內容
ABAQUS混凝土塑性損傷模型（CDP模型）excle簡便 ￥20
本excle簡捷易懂，只需在excle表中更改彈模以及軸心抗壓強度自動生成數據，表中列出了公式以及只需要輸入ABAQUS中的數據，十分容易上手
ABAQUS UMAT-混凝土受拉狀態下塑性損傷模型的簡單實現 ￥600
本構模型的實現算法摘抄自DeBorst的書籍《Nonlinear Finite Element Analysis of Solids and Structures》，基本如下：
為了簡化模型，筆者將書中損傷部分做了簡化，不再采用損傷屈服面進行判定。損傷影子w的計算直接由塑性等效應變確定。
在ABAQUS中建立100*100*100的立方體塊，試件的底部固定，頂部反復加載-卸載，通過UMAT得到的模擬結果如下：
混凝土塑性損傷模型 ￥5.99
<p>混凝土塑性損傷模型是基于拉、壓各向同性塑性的連續線性損傷模型，用于描述混凝土的非線性行為。采用通用有限元分析軟件ABAQUS/Standard分析，在此軟件中的混凝土塑性損傷模型具有以下特點：</p><p>1. 適于各種單元（梁、桿、殼、實體）的混凝土或其他類似的脆性材料的模擬，用于殼元時，沿厚度方向的積分點數達到9個通常可以保證計算的準確性；</p><p>2. 雖然它主要致力于鋼筋混凝土結構的分析，但可以用于素混凝土；</p><p>3. 可用“rebar”選項模擬混凝土中的鋼筋；</p><p>4. 適于低圍壓下混凝土單調、往復和動力荷載下的計算；</p><p>5. 是非相關多軸硬化塑性和各向同性線性損傷模型的綜合，用于描述由于混凝土斷裂引 起的不可恢復的損傷；</p><p>6. 允許循環加載過程中用戶對于剛度恢復進行控制；</p><p>7. 可定義與應變速率的相關性；</p><p>8. 應用粘性系數修正，可提高軟化階段的收斂效率；</p><p>9. 要求材料的彈性行為應為各向同性且為線性的。</p><h2>1 線性損傷模型與塑性模型</h2><p>本節簡要介紹構成混凝土塑性損傷模型的線性損傷模型與塑性模型(Hibbitt等，2003)。</p><h2>1.1 線性損傷模型</h2><p>混凝土塑性損傷模型包括混凝土受拉開裂和壓碎兩種破壞機制，分別由<span style="background-color: yellow;">等效拉壓塑性應變</span>決定。單軸應力-應變關系轉變為<span style="background-color: yellow;">應力與塑性應變</span>的曲線。
展開 huang晶體塑性umat耦合Johnson-cook 損傷模型，實現晶體材料彈-塑-損傷模擬分析
Johnson-cook 損傷起始準則是延性損傷準則模型的一個特例，用于預測延性金屬中孔洞的形核、生長和聚結導致的損傷起始。該模型假設損傷開始時的等效塑性應變是應力三軸性和應變率的函數。同時可以考慮溫度的影響。
包含的材料參數有：
失效相關參數：d1-d5。
ABAQUS混凝土損傷塑性模型損傷因子對本構關系影響 附c40~c45混凝土損傷因子ABAQUS輸入
但是ABAQUS塑性損傷模型除了能模擬單調加載的混凝土行為外，更重要的功能就是模擬循環、動態荷載下的混凝土反應，在結構的抗震性能分析能起到很好的作用。
在動荷載作用下，混凝土在受力過程中拉伸和壓縮都會產生損傷造成的裂縫開展，從而導致材料剛度退化。CDP 模型就假定混凝土材料主要因為拉伸開裂和壓縮破碎而破壞，拉伸和壓縮采用不同的損傷因子來描述這種剛度退化，詳見圖 1、圖 2。
圖中E0是材料初始未受損的彈性剛度。損傷變量dc和dt分別為壓縮和拉伸條件下的損傷因子，表示彈性剛度的退化。損傷后的彈性模量為（1-dc）E0，或（1-dt）E0。損傷因子dc或dt=0時表示沒有損傷，dc或dt=1時表示材料失去強度。
那么混凝土的塑性損傷本構模型中的損傷因子到底對混凝土的應力-應變曲線有什么影響呢？讓我們采用100mm*100mm*300mm的混凝土棱柱體模型來做個測試看一下。
依然采用C110級混凝土的本構關系，混凝土的屈服應力和非彈性應變表格如下。子選項中損傷參數和非彈性應變關系的表格也在圖中給出。
但是注意上圖中紅色框部分默認是不填的，即下圖中的混凝土壓縮損傷——拉伸恢復因子wt，混凝土拉伸損傷——壓縮復原因子wc，默認是不填的。
因為CDP模型假定混凝土從拉伸到壓縮時裂縫會閉合，剛度會恢復；從壓縮到拉伸時裂縫仍然存在，剛度不會恢復。因此在ABAQUS中不填的話默認wt(拉伸剛度恢復因子)=0，wc(壓縮剛性恢復因子)=1.
下圖為損傷因子和剛度恢復因子在混凝土載荷循環中對混凝土本構模型的影響。
展開 求助：混凝土塑性損傷模型
求助：請問混凝土塑性損傷模型的損傷因子什么情況下可以不定義，什么情況下要定義？比如我現在單調受彎加載該構件，那該構件（塑性損傷模型）需不需要定義損傷因子？
混凝土塑性損傷模型 ￥1
單軸壓縮下塑性損傷發展：
應力—應變曲線：
混凝土塑性損傷模型的材料參數：
混凝土塑性損傷模型 ￥15
如圖 紅色部分輸入即可 采用韓林海本構模型 藍色部分直接復制粘貼
受拉使用應力斷裂能 輸入簡單
整個模型收斂速度快 擬合程度較好
【JY】JYCDP插件：ABAQUS混凝土CDP模型插件分享 | 混凝土損傷塑性模型 ￥59.9
(推薦6.14、2016版本，仍保留固流分析耦合模塊，后版本取消該模塊)，
文后附 6.14-4 軟件下載鏈接及子程序相關下載，
【簡介】
為簡便鋼筋混凝土構件或者結構的本構模型設置，本期給大家推薦一款Abaqus混凝土CDP模型插件，供大家應用參考。這個插件無需繁瑣的Excel操作，僅需選擇混凝土等級即可在Abaqus前處理界面一鍵生成混凝土CDP本構曲線，且可任意調整本構曲線長度，并可對極限強度進行修正，且適用于不同的力、位移單位，可用于各類混凝土構件及結構的精細化分析。
對于鋼筋混凝土構件或者結構而言，正確合理的本構模型是對構件或結構進行非線性分析的關鍵。ABAQUS提供三種混凝土本構關系模型，分別為脆性開裂模型、彌散開裂模型及損傷塑性模型，其中，混凝土損傷塑性 (Concrete Damaged Plasticity，CDP)模型是通過將各向同性下損傷彈性與拉伸和壓縮塑性相結合的方式來對混凝土的非彈性行為進行描述的，適用于Standard和Explicit兩大求解模塊，可用于模擬混凝土在任意荷載作用下的受力情況，同時考慮了由于拉、壓塑性應變導致的彈性剛度的退化以及循環荷載作用下剛度的恢復，具有較好的收斂性。有關CDP模型的介紹及應用可見推文：
【JY】淺談混凝土損傷模型及Abaqus中CDP的應用
【程序可解決的問題】
采用ABAQUS模擬梁柱節點時，ABAQUS中CDP模型損傷系數計算到0.9和損傷系數計算到0.99所得的滯回曲線相差甚大，筆者建立了現澆梁柱節點模型對此進行了驗證。
CDP模型本構曲線末尾段的選取，對滯回曲線下降段的影響較大。
展開 再生混凝土塑性損傷本構模型 ￥20
由于再生混凝土相較于普通的混凝土 峰值應變更高，所以不可以采用普通的混凝土本構，本文檔采取了楊有福的公式，基于韓林海的進行了改進，適合鋼管混凝土構件
Abaqus中混凝土損傷塑性模型（CDP） ￥40
制作了一個計算混凝土受壓受拉本構和損傷因子的Excel表格，包含
一、4種混凝土本構計算方法
①10規范
②02規范（過鎮海模型）
③丁發興模型
④ConcreteD模型
二、3種損傷因子計算方法
①圖解法
②能量法
③規范法
三、鋼筋雙折線模型
四、Mander模型計算核心區混凝土
Excel部分截圖如下：
我做的Excel表格優點有：
1、附帶公式，清晰明了，簡單易懂；
2、內容全面，含各種混凝土及損傷因子的計算方法；
3、修改簡單，適用于各級混凝土，可操作性強。
需要的話可付費購買，下載Excel即可。
視頻地址如下：
http://www.yqgqt.org.cn/college/video/c14528
1、4種混凝土本構計算方法及對比；
2、3種損傷因子計算方法及對比；
3、真實應力、應變與名義應力、名義應變之間的轉換；
4、本構及損傷因子的選取與截斷。
注意，視頻和帖子價格不同。
如有不懂問題，加我微信YClarie問我即可。
若你不想購買，讓你少走彎路：給你提供學習的資料如下：
計算不難，編制Excel表格也不難！可自己嘗試
展開 混凝土損傷塑性模型在abaqus里的模擬
混凝土規范上把混凝土本構用分段函數來表達的，有上升段和水平段，極限壓縮應變為0.0033（日本0.0035，與中國規范上升段化簡后一樣）如下圖：
中國規范：
日本規范：
而實際上混凝土本構包含上升段和下降段，特別對于下降段，混凝土內裂縫逐漸發展，卸載時彈性軟化，非線性彈性和彈塑性理論很難描述這一特性。
損傷模型可以考慮由于損傷積累而導致裂縫逐漸擴展，從而導致混凝土剛度軟化。混凝土塑性損傷模型可以用來分析混凝土結構在低周往復荷載和動載作用下，混凝土微觀損傷導致剛度軟化的模擬.
abaqus里CDP-model的應力應變關系可以采用《混凝土結構設計規范》附錄C給出的本構關系（如下圖）：
本次模擬采用C50,具體參數如下：
與規范應力應變關系對比圖（模擬的不太理想，以后改進）：
混凝土受壓卸載時剛度軟化示意圖如下:
本次abaqus模擬的混凝土試件往復加載卸載受壓時的應力應變關系履歷如下圖所示，無損傷混凝土卸載剛度與初始彈性模量一致（綠色），也是我們通常的處理方式，而考慮損傷的混凝土卸載剛度變小，且隨著塑性應變的發展，剛度越來越小（黃綠色，圖中斜線方便觀察剛性斜率而設）。
來源：有限元仿真
展開 