循環(huán)塑性本構模型的案例
ABAQUS中UMAT中的經典循環(huán)塑性本構模型及相應子程序代碼 ￥20
ABAQUS中UMAT中的循環(huán)塑性模型，包含非線性各向同性強化彈塑性、線性各向同性強化彈塑性、線性隨動強化彈塑性模型，包含CAE文件、UMAT文件等。
約束混凝土cdp塑性損傷本構，mander混凝土本構模型 ￥10
約束混凝土本構，mander混凝土本構，自己做的箍筋約束方柱和圓柱本構模型，表格只要輸入相關參數(shù)，自動生成ABAQUS塑性損傷本構關系。
ABAQUS umat 理想彈塑性本構模型 ￥99
<p class="ql-align-justify"><span style="color: rgb(15, 17, 21);">本資源包含一份 PDF 文檔和可直接編譯運行的 Fortran UMAT 代碼，具體內容為：</span></p><p class="ql-align-justify">理想彈塑性本構 + 隱式積分 + 徑向返回</p><p class="ql-align-justify">完整公式推導 + Fortran 源碼直接編譯</p><p class="ql-align-justify">von Mises 屈服+ 一致切線模量全實現(xiàn)</p><p class="ql-align-justify">PDF 包含規(guī)范化的本構方程、隱式積分、徑向返回與一致切線模量推導，可供初學者學習。配套 UMAT 代碼可直接在 ABAQUS 編譯運行，采用全隱式積分搭配一致切線模量，收斂速度極快、計算精度極高，<span style="background-color: rgba(0, 0, 0, 0);">適合初學者快速入門。</span></p><p class="ql-align-justify"><span style="background-color: rgba(0, 0, 0, 0);">下圖展示了</span><span style="color: rgb(25, 27, 31);">部分</span><span style="background-color: rgba(0, 0, 0, 0);">PDF內容，及umat計算結果與abaqus內置模型對比，可以發(fā)現(xiàn)umat收斂速度極快，與abaqus內置模型幾乎一致。
展開 運用ABAQUS軟件對冰材料彈塑性本構模型改進及驗證（附源文件） ￥1300
<p class="ql-align-justify"><strong>內容：</strong></p><p class="ql-align-justify">基于參考文獻通過ABAQUS建立了冰材料彈塑性本構模型；對比已有試驗，對比裂紋演化現(xiàn)象和沖擊載荷曲線，驗證了冰材料本構模型的有效性。</p><p class="ql-align-justify"><img src="https://img.jishulink.com/202507/attachment/7b0d26ab81f645dc98e8b15335447247.png" width="1027"></p><p class="ql-align-justify"><br></p><p class="ql-align-justify"><br></p><p class="ql-align-justify"><br></p><p class="ql-align-justify"><br></p><figure style="text-align: center;" class="ql-align-center"><img src="https://img.jishulink.com/202510/attachment/7cbe0c886d1d4de59fdee40d233200d8.png" style="" width="616" data-mobile-src="https://img.jishulink.com/202510/attachment/7cbe0c886d1d4de59fdee40d233200d8.png?
展開 
金屬材料塑性本構模型（結合workbench）
工程中的金屬結構一般都處于彈性工作狀態(tài)，所以工程金屬結構分析大多數(shù)都使用線彈性材料本構模型。不過，塑性本構也是應該掌握的。
workbench中常見的四種塑性本構模型
涉及三個方面：
01 雙線性/多線性（bilinear / multilinear）
02 強化（hardening）
03 等向和隨動（isotropic / kinematic）
如圖所示：
01 雙線性和多線性的區(qū)別是一目了然的，即應力應變曲線是兩條折線或兩條以上折線（三條及以上）。
02 強化是指材料在屈服后，應力隨應變還會增加，與此相對應的是理想彈塑性，材料屈服后，應力不隨應變增加。
03 拉伸屈服點對壓縮屈服點存在影響（初始屈服影響后繼屈服）。等向模型中壓縮屈服點等于上一次最大拉應力；隨動模型中壓縮屈服點等于兩倍屈服應力減去上一次最大拉應力。由此可知，隨動和等向模型定義的是材料屈服條件的變化，在材料加載后卸載再加載的情況下（多次屈服）才發(fā)揮作用。對于單調加載（不存在卸載過程），實際起作用的定義只是雙線性強化或者多線性強化。
另外，材料的屈服條件（屈服面）也有不同的描述模型。比如Tresca屈服準則，Mises屈服準則，D-P屈服準則等。例如，對于二維應力狀態(tài)，Mises屈服準則在主應力空間中是橢圓形；對于三維應力狀態(tài)，Mises屈服準則在主應力空間中是圓柱形。
展開 再生混凝土塑性損傷本構模型 ￥20
由于再生混凝土相較于普通的混凝土 峰值應變更高，所以不可以采用普通的混凝土本構，本文檔采取了楊有福的公式，基于韓林海的進行了改進，適合鋼管混凝土構件
非局部晶體塑性本構模型實現(xiàn)與案例演示
在FCC晶體中，有12種滑移系統(tǒng)可能在塑性變形過程中被激活
通常，樣品的晶體學和應力狀態(tài)是決定滑移系統(tǒng)是否活躍的主要因素。試驗過程中試樣所經歷的塑性來自于激活滑移系統(tǒng)的貢獻。臨界分辨剪切應力是確定晶體滑移開始的標準，而FCC金屬材料塑性變形主要由位錯滑移貢獻。以位錯為內變量的本構方程可以對多晶材料的塑性變形做出更加物理的描述和預測，并與微尺度的實驗進行對比分析。
Ma和Roters引入的基于位錯密度的本構模型（Ma和Roter，2004；Ma、Roters和Raabe，2006a，b）使用移動位錯ρmα，沿著滑移系統(tǒng)α滑動，以適應部分外部塑性變形，在基于位錯的模型中，Orowan方程通常代替唯象的冪律流動方程
其中ρm是統(tǒng)計儲存位錯密度，b是伯格斯矢量，v是可移動位錯密度平均速度，統(tǒng)計儲存位錯密度表示為初始統(tǒng)計位錯密度和變形過程中統(tǒng)計位錯密度增量之和，統(tǒng)計位錯密度演化表示為
其中dαβ是位錯增殖相互作用張量，kc和knc分別作為控制共面和非共面滑移系統(tǒng)相互作用系數(shù)大小的常數(shù)。量rαc是位錯湮滅的位錯捕獲半徑，并隨溫度和變形速率的變化（Kocks，1976），通常使用考慮統(tǒng)計位錯密度的本構模型，即從一個材料點的加載歷史可以充分描述本構行為。對于多晶體的應力-應變曲線和織構預測，溫度效應，局部位錯模型已被證明是強大和有效的。
然而，如果模擬規(guī)模變小，例如在專注于納米壓痕（Zaafarani et al.，20082006）和微柱壓縮（Raabe，Ma和Roters，2007a）的研究中，則局部模型可能由于無法描述尺寸效應而不足，較小晶粒尺寸的強化效應是由于晶界附近非均勻塑性變形的體積分數(shù)較高。
展開 UMAT子程序（晶體塑性力學）經典案例-單胞模型（inp+UMAT文件+子程序對應的本構模型文件） ￥10
這份資料是從事固體力學研究幾年經典推薦教程，助你在有限元仿真理論部分有更深入的理解和認識，同時對有限元材料本構模型的UMAT子程序的編寫、材料參數(shù)的設置、ABAQUS的前處理有更加深入的感悟。
附件中文件：inp單胞模型，UMAT晶體塑性經典子程序，子程序對應的經典說明，固體力學國際經典教程（Computational Methods for Plasticity），ABAQUS工程項目前處理經典教程。
看好再購買，售出不退，謹慎入手
Abaqus Umat (子程序4) 彈塑性本構,米塞斯模型(Mises Model) ￥10
1）米塞斯模型為經典的彈塑性本構，主要用來模擬金屬材料在外荷載作用下的彈塑性行為
2）具體為金屬在各向均勻受壓狀態(tài)下不會產生塑性變形，只有在剪切作用下會發(fā)生塑性變形
3該Fortran代碼為Abaqus的外接子程序(user subroutine)，可用于學習最簡單的彈塑性本構的編寫過程
米塞斯模型的適用范圍及屈服面形狀
所編寫米塞斯模型UMAT子程序
非局部晶體塑性本構模型實現(xiàn)與案例演示
在FCC晶體中，有12種滑移系統(tǒng)可能在塑性變形過程中被激活
通常，樣品的晶體學和應力狀態(tài)是決定滑移系統(tǒng)是否活躍的主要因素。試驗過程中試樣所經歷的塑性來自于激活滑移系統(tǒng)的貢獻。臨界分辨剪切應力是確定晶體滑移開始的標準，而FCC金屬材料塑性變形主要由位錯滑移貢獻。以位錯為內變量的本構方程可以對多晶材料的塑性變形做出更加物理的描述和預測，并與微尺度的實驗進行對比分析。
Ma和Roters引入的基于位錯密度的本構模型（Ma和Roter，2004；Ma、Roters和Raabe，2006a，b）使用移動位錯ρmα，沿著滑移系統(tǒng)α滑動，以適應部分外部塑性變形，在基于位錯的模型中，Orowan方程通常代替唯象的冪律流動方程
其中ρm是統(tǒng)計儲存位錯密度，b是伯格斯矢量，v是可移動位錯密度平均速度，統(tǒng)計儲存位錯密度表示為初始統(tǒng)計位錯密度和變形過程中統(tǒng)計位錯密度增量之和，統(tǒng)計位錯密度演化表示為
其中dαβ是位錯增殖相互作用張量，kc和knc分別作為控制共面和非共面滑移系統(tǒng)相互作用系數(shù)大小的常數(shù)。量rαc是位錯湮滅的位錯捕獲半徑，并隨溫度和變形速率的變化（Kocks，1976），通常使用考慮統(tǒng)計位錯密度的本構模型，即從一個材料點的加載歷史可以充分描述本構行為。對于多晶體的應力-應變曲線和織構預測，溫度效應，局部位錯模型已被證明是強大和有效的。
然而，如果模擬規(guī)模變小，例如在專注于納米壓痕（Zaafarani et al.，20082006）和微柱壓縮（Raabe，Ma和Roters，2007a）的研究中，則局部模型可能由于無法描述尺寸效應而不足，較小晶粒尺寸的強化效應是由于晶界附近非均勻塑性變形的體積分數(shù)較高。
展開 混凝土塑性損傷本構模型-Matlab代碼 ￥50
本塑性損傷本構模型參照規(guī)范GB50010-2010（2015版），根據(jù)代碼注釋修改參數(shù)即可直接計算應力，應變和損傷因子具體數(shù)值以及繪制應力-應變關系曲線和損傷因子曲線。損傷因子的計算參照Sidiroff的能量等價原理，主要應用于ABAQUS中混凝土材料屬性中塑性應力應變關系和損傷因子的取值。下載包中還自帶規(guī)范PDF完整版。
如有需求，本人可提供相關abaqus技術指導。
下圖為利用本代碼所計算得出的C40標號混凝土材料受壓應力-應變關系曲線以及損傷因子曲線：
通過UMAT實現(xiàn)基于DP屈服準則的改進西原模型的三維粘彈塑性（蠕變）本構模型
傳統(tǒng)西原模型是目前可以比較好地描述巖石蠕變過程曲線的元件模型，但是，西原模型使用的元件為黏彈、黏塑性元件（如圖1），難以描述巖石屈服破壞后進入加速階段的蠕變變形。滑坡預報，特別是臨滑預報在地質災害防治領域具有重要意義。
通過編寫abaqus UMAT子程序，可得到如下結果：
（1）應力狀態(tài)較小時，僅發(fā)生彈性應變和粘彈性應變，最后隨時間趨于穩(wěn)定值。
（2）單元屈服時，發(fā)生粘彈塑性應變，應變隨加載時長逐漸增加，但尚未達到觸發(fā)應變，曲線呈現(xiàn)兩階段特性。
（3）隨著加載時長的增加，應變進一步增加，超越觸發(fā)應變后，進入快速蠕變階段，應變快速增加，曲線呈現(xiàn)三階段蠕變特性。
參考文獻：
[1] 齊亞靜, 姜清輝, 王志儉, 等. 改進西原模型的三維蠕變本構方程及其參數(shù)辨識[J]. 巖石力學與工程學報, 2012, 31(2): 347-355.
[2] 沈才華, 張兵, 王媛, 等. 基于DP屈服準則的西原本構模型及其運用[J]. 地下空間與工程學報, 2016, 12(2): 402-407.
最后，有需要歡迎通過微信公眾號聯(lián)系我們。
微信公眾號：320科技工作室。
展開 混凝土塑性損傷CDP模型的幾個問題 附2010規(guī)范用C50混凝土損傷塑性本構關系數(shù)據(jù)下載
下載地址：2010規(guī)范用C50混凝土損傷塑性本構關系數(shù)據(jù)
Abaqus鋼管混凝土塑性受拉及受壓應力應變本構模型及損傷因子 ￥5
<p class="ql-align-justify">本內容基于韓林海的約束混凝土模型所制作的Excel，可用于將其輸入直接到ABAQUS中，用于建立鋼管約束混凝土型，具體如下：</p><p class="ql-align-justify">模型介紹：</p><p class="ql-align-justify">本模型基于<span style="color: rgb(25, 27, 31);">韓林海</span>所開發(fā)的約束混凝土應力應變關系模型，以及損傷因子，其中受壓本構以及受拉本構以及其損傷因子均有，且附帶鋼材料的二次流塑模型，可直接輸入abaqus進行分析，均具有完美下降段。</p><p class="ql-align-justify"><br></p><figure style="text-align: center;" class="ql-align-center"><figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202511/attachment/66d3ae0e7f464f3c8a0a386084e4e134.png" style="display: inline-block;" data-regular="true"><img src="https://img.jishulink.com/202511/attachment/66d3ae0e7f464f3c8a0a386084e4e134.png" data-mobile-src="https://img.jishulink.com/202511/attachment/66d3ae0e7f464f3c8a0a386084e4e134.png?
展開 abaqus三維復合材料彈塑性+漸進損傷本構模型-3D VUMAT ￥145
對于纖維增強復合材料的模擬，在<a href="/major/ABAQUS中，集成了二維Hashin失效準則與多種損傷演化準則，但缺少三維的復合材料本構模型。
參考一篇已發(fā)表的SCI文章，使用Fortran語言建立三維平紋織物復合材料彈塑性、漸進損傷本構模型-Vumat子程序。平紋織物復合材料在1方向和2方向絲束性能近似相同。
該程序是博士期間學習復材子程序的小部分總結，編程結構并不是非常漂亮及完美，但確保能順利運行，且單元驗證結果與理論公式一致，介意請勿拍。
程序中塑性迭代部分并非主流的牛頓-拉夫遜和梯度下降方法，但經過驗證能夠適用于該模型，介意請勿拍。
附件內容：1. inp算例模型(低速沖擊工況，1/4模型，層間使用cohesive element) 2. 子程序 3 .使用方法 4.參考論文名稱
首先介紹該子程序的使用方法與效果
1. 在ABAQUS中建立三維復合材料模型，這里建立一個簡單的方塊。賦給材料方向，1,2方向分別表示絲束的方向，3方向表示垂直于1，2的方向，也就是面外方向。
2. 建立材料屬性
3. 建立顯示Explicit計算時間步，時間0.005，在場輸出中勾選輸出 SDV和 STATUS.
4. 劃分網格，賦給Explicit 3D stress單元類型，邊界條件根據(jù)需要設定即可，此處不再贅述。此處劃分為一個單元，使用12方向往復加載卸載。建立Job，提交模型前在Job中選擇該子程序，設置雙精度計算。
5. 查看結果，等效塑性應變在卸載時沒有變化，再次加載時剪切應力按照原來的路徑返回，剪切損傷在卸載時也保持不變。
6. 將該子程序應用在低速沖擊模型中，可以順利運行。
接下來簡要介紹該子程序的相關理論，子程序、參考的論文名稱以及輸入材料參數(shù)的對應含義打包在附件中。
展開 