彈塑性本構開發
關注創建者：匿名 創建時間：2026-01-05
彈塑性本構開發的視頻教程
土體彈塑性本構理論（臨界狀態理論，劍橋模型，狀態相關本構，邊界面模型）
本套課程將由淺入深教大家一些土體本構方面的理論知識。本構就是用數學公式描述土體的復雜特性，因此掌握一些本構相關知識對了解土體性質，進而開展巖土工程問題的研究具有重要意義。同時本課程的課件及參考資料都會在附件里贈送給大家。
￥68 4小時25分鐘 32463播放
查看
ABAQUS VUMAT 一維彈塑性本構VUMAT
本課程主要圍繞ABAQUS VUMAT在一維彈塑性本構中的應用展開，具體內容模塊如下： 第一部分：理論基礎夯實 各向同性線彈性本構：回顧基礎彈性理論。 屈服準則與流動法則：理解材料何時屈服及塑性流動方向。 應變分解與更新：掌握彈性應變與塑性應變的分離處理。 徑向返回法 ：學習顯式積分中處理彈塑性問題的核心數值算法。 牛頓迭代法：掌握求解非線性方程組的數值方法。
￥99 2小時28分鐘 31播放
查看
鋼結構S1 ABAQUS經典金屬彈塑性本構及模擬應用
ABAQUS經典金屬彈塑性本構及模擬應用 課程簡介： 本課程通過低碳鋼拉伸試驗及工字鋼梁低周往復模擬，詳細講解ABAQUS中經典金屬彈塑性本構計算及應用方法。 適用對象： ABAQUS入門學習者、ABAQUS鋼結構模擬用戶、土木工程、結構工程在校學生及從業人員。 模擬結果預覽： 購課提示： 購課后請登錄網頁版技術鄰，在課程下方下載課程相關資料。
￥100 1小時31分鐘 7364播放
查看
彈塑性本構開發的實例教程
問題介紹
木材的本構模型是采用連續體單元建模模擬木材彈塑性響應的基礎，然而木材復雜的力學性質常常為其本構模型的建立帶來困難。木材力學性質的復雜性主要表現在：
不同方向的強度值和剛度值各不相同；
2. 同一方向的抗拉強度和抗壓強度之間存在差異；
3. 不同形式荷載作用下材料的響應不同，壓力作用下材料的表現以延性為主， 而拉力和剪力作用下材料的破壞呈脆性。
木材在復雜應力狀態下的彈塑性本構模型。以經典彈塑性力學為框架，該本構模型建立在如下四個基本假設的基礎之上：
木材在彈性階段是理想的橫觀各向同性材料；
2. 材料的屈服符合簡化的 Hashin 屈服準則；
3. 材料在受拉和受剪屈服之前是理想線彈性的，屈服之后進入塑性流動階段；
4. 材料受壓初始屈服之前是理想線彈性的，屈服之后進入應變硬化階段，隨 著屈服面的轉移到達最終屈服面后進入完全的塑性流動。
二。子程序編寫流程
本工作室在三維hashing模型的基礎上，利用Abaqus軟件平臺，開發了完整的木材的彈塑性本構umat子程序，包含木材完整的彈性、塑性、強化以及軟化階段。編寫子程序的流程如下：
三。結果驗證
通過如下圖的木材模型進行驗證：
該模型在受壓、受剪及受拉的工況下，應力應變曲線如下所示：
該子程序還有以下特征：
能計算靜力非線性
2. 收斂性好
3. 能計算復雜應力狀態
附件為本子程序參考的文獻，供大家學習探討~
2. 木材的力學性質試驗研究及數值模擬方法.pdf
最后，大家有關于Abaqus二次開發的相關需求可以添加管理員扣扣：3045552826，微信：CAE320，同時也歡迎大家關注“320科技工作室”的微信公眾號，掃一掃二維碼即可關注~~
展開 RVE代表體積單元
RVE模型由纖維和樹脂構成，一般假設纖維是橫觀各向同性線彈性材料，樹脂則為彈塑性材料。本文通過在屈服準則中引入拉壓非對稱參量，研究了樹脂的拉壓不對稱彈塑性損傷行為。
由于樹脂的屈服行為與靜水壓力相關，這里采用下式所示的拋物面屈服準則。
式中J2為偏應力的第二不變量，I1為應力第一不變量，σt和σc為拉壓屈服應力
采用非關聯塑性流動準則，如下所示。
式中σvm為mises等效應力，P為靜水壓力，α為材料參數
損傷萌生準則如下所示
式中J2和I1為無損應力下的不變量。
為了降低模型的網格依賴性，損傷演化采用特征長度相關的指數模型
式中，rm為損傷內變量，am為特征長度相關的材料參數。
Melro的文章中給出了通過Simpson積分和弦截法計算Am的方法，實際計算發現通過該方法計算的am效果不是太理想，因此本文未對am進行迭代，直接采用其初值進行仿真計算，如下所示。
計算流程如下
計算流程圖
根據上文的彈塑性損傷模型編寫了vumat子程序，并通過單胞模型進行了驗證，計算結果如下圖所示。
abaqus單胞模型
拉伸載荷下的應力應變曲線
壓縮載荷下的應力應變曲線
展開 在寫彈塑性本構之前，我對塑性流動是干嘛使的沒有直觀概念。寫的時候我才明白，由于只能先算出來等效塑性應變，沒有流動方向的話，就無法把它轉換到各個應變分量，不知道應變分量就無法計算應力。這玩意從數學上講，是一個轉換公式。
我們目前重工業上大部分的結構材料還是金屬，盡管ABAQUS中有自帶的JC模型，但是如果要模擬更復雜的情況，學會寫彈塑性本構就十分必要。
本期就給一個彈塑性VUMAT拉伸失效的案例，結合單元刪除技術，模擬結構破壞過程。
本構模型
采用經典老演員JC模型描述本案例的彈塑性本構：
為了模擬結構破壞，采用如下準則判斷單元完全失效，滿足其一即可：
（1）材料Mises應力達到極限值；
（2）材料極限應變達到極限值。
子程序結構
子程序的基本結構如下：
1.初始化準備工作
程序首先進行初始化準備工作，讀入材料的彈性參數、強度參數、硬化參數以及應變率相關參數，然后構建彈性剛度矩陣，為后續計算奠定基礎。
2.進入材料點循環
接下來進入材料點循環，對每個積分點逐一進行計算。對于每個材料點，程序首先讀取上一步的狀態變量，包括累積的等效塑性應變、應力狀態以及背應力等內部變量。
3.失效判斷
程序隨后進行失效判斷，檢查材料是否滿足失效準則。判斷依據包括兩個方面：一是等效塑性應變是否超過極限應變閾值，二是等效應力是否達到破壞強度。一旦滿足任一失效條件，程序將材料標記為失效狀態，并大幅降低其剛度以模擬材料的承載能力喪失。
4.本構響應計算階段
在本構響應計算階段，程序考慮了應變率效應和材料硬化特性，更新當前的屈服應力。同時計算應力偏量，得到米塞斯等效應力和塑性流動方向，這些是判斷材料是否屈服的關鍵參數。
5.彈塑性判別
然后進行彈塑性判別。
展開 <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 屈服+ 一致切線模量全實現</p><p class="ql-align-justify">PDF 包含規范化的本構方程、隱式積分、徑向返回與一致切線模量推導，可供初學者學習。配套 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內置模型對比，可以發現umat收斂速度極快，與abaqus內置模型幾乎一致。
展開 在連續介質力學中，所有問題(包括運動、應力、應變以及守恒定律等)既可用物體變形前的初始構形B為參照構形(取x1為自變量)來描述，又可用物體變形后的新構形，B'為參照構形(取x1*為自變量)來描述，前者稱為拉格朗日(LagrangeJ L)描述，后者稱為歐拉(Euler L)描述。
在固體力學中，我們常采用拉格朗日描述;在流體力學中采用歐拉描述更為方便;而對大變形問題及一般的物理定律，采用拉格朗日坐標來建立它的數學表達式更為方便，但在求解具體問題時，又常以歐拉描述更方便，所以兩種描述都要采用。
—End—
CAE仿真與數值模擬微信公眾號，主要介紹CAE仿真與數值模擬的知識與應用公眾號主要介紹CAE仿真與數值模擬的知識與應用。通過論壇，博客，論文，案例等為大家帶來知識食糧。仿真軟件:abaqus、ansys、flunet、comsol、hypermesh、moldflow等，涉及領域有機械材料土木物理等。
展開 
彈塑性本構開發的相關專題、標簽、搜索
彈塑性本構開發的最新內容
<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
并不簡單的彈塑性本構子程序6個月前
初學材料力學就知道最常見的金屬一般都是彈塑性的。
所謂彈塑性，就是把材料性能劃分成了兩個階段，前面的階段是彈性，比較好理解，載荷與變形線性變化。后面塑性，就是指材料繼續變形，但是載荷不往上走了，或者即便走也變慢了。而且即便完全卸載，第二個階段的變形仍然會保留。
材料如此，人亦如此，過度消耗是補不回來的。彈塑性材料有屈服強度這個概念，就是指進入塑性后，本來向上的曲線開始低頭了
1 vumat與umat的區別
從程序實現的角度，我們重點關注以下幾點區別：
? vumat不需要輸出一致性切線剛度矩陣
? vumat中應力應變存儲順序與umat不同
? vumat中存儲的應變值為張量應變值，而umat中為工程應變
? vumat的應力和狀態變量的更新方式不同，其分為old和new兩個數組
Abaqus/Explicit在啟動計算前，會進行數據檢查
1 本構理論
本文講解如何將三維的率無關彈塑性理論應用到平面應力問題中。對于平面應變和軸對稱問題，由于是相應的應變分量為0，因為可以直接使用三維的本構，只需將相應的應變分量設為0作為本構的輸入即可。然后，對于平面應力問題，是相應的應力分量為0，由于本構是由應變驅動求得對應的應力，相應應力分量為0相當于對系統施加了相應的約束，因此三維的本構理論不可直接應用于平面應力問題中，需要將相應的約束考慮其中進行求解
<p class="ql-align-justify"><strong>內容：</strong></p><p class="ql-align-justify">基于參考文獻通過ABAQUS建立了冰材料彈塑性本構模型；對比已有試驗，對比裂紋演化現象和沖擊載荷曲線，驗證了冰材料本構模型的有效性。</p><p class="ql-align-justify"><img src="https://img.jishulink.com
<p>1 本構理論</p><p>1.1 率形式</p><p>本構方程為：</p><p><img src="https://img.jishulink.com/msimage/202402/dd344f7d07bc21faf8fcc073781a5aa2.png"></p><p>單軸拉伸的應力應變的硬化曲線如下：</p><p><img src="https://img.jishulink.com/msimage
各向同性硬化von Mises率無關彈塑性本構理論以及umat源代碼
1 本構理論
1.1 率形式
對于各向同性線彈性材料，其本構方程為：
式中假設了應變張量可以分解為彈性應變和塑性應變兩部分：
因此塑性本構的關鍵在于計算塑性應變的演化。對于率無關彈塑性的本構理論，需要確定以下三個部分：
(1)：屈服條件
(2)：流動法則
(3)：硬化法則
在此采用的是
混凝土在外載荷作用下的非線性行為中同時包含微裂縫和塑性流動這兩種微觀機制的影響。在考慮混凝士等準脆性材料的非彈性力學行為方面，連續損傷力學模型可以通過不同的方式來描述材料剛度和強度的退化以及單邊效應。真正意義上的彈塑性損傷本構模型：不僅考慮卸載時不可恢復塑性變形的影響，而且還應該考慮損傷和塑性的雙向耦合效應。
彈性階段應力應變滿足如下關系
通過對應力進行譜分解，可得
式中，
對于纖維增強復合材料的模擬，在<a href="/major/ABAQUS中，集成了二維Hashin失效準則與多種損傷演化準則，但缺少三維的復合材料本構模型。
參考一篇已發表的SCI文章，使用Fortran語言建立三維平紋織物復合材料彈塑性、漸進損傷本構模型-Vumat子程序。平紋織物復合材料在1方向和2方向絲束性能近似相同。
該程序是博士期間學習復材子程序的小部分總結，編程結構并不是非常漂亮及完美
我在隧道二次襯砌中采用了cdp模型，并進行了開挖模擬，但是開挖完成后損傷參數為0，沒有變化，這是為什么呢，很它這里的警告有關嗎，我看別的文獻里面參數值雖小，但是是有值的呀
