塑性模型的案例
各版本亞塑性模型matlab及umat程序 ￥150
Wu和Kolymbas（1990）首先創建了亞塑性理論的基本方程,他用一個簡單的非線性張量函數來模擬非彈性材料的性能,其應力率由兩項組成：一項跟應變率呈線性關系,另一項跟應變率呈非線性關系,其非線性關系通過應變率的范數體現出來。Wu和Bauer（1994）結合無粘性土的變形特性,對亞塑性理論基本方程中關于應變率的線性部分和非線性部分的具體形式作了一些改進,于年提出了第一個較為實用的亞塑性模型—四參數亞塑性模型，又稱為Wu-Bauer亞塑性模型。為考慮密實度對土體響應的影響，Wu等（1996）進一步將臨界狀態理論引入到亞塑性模型中。隨后，為了更好的描述臨界應力狀態，Gudehhus（1996）、Bauer（1996）、Wolffersdorff（1996）將不同的屈服準則引入到模型中，大大提升了模型的模擬能力。為了更好的模擬循環荷載，Niemunis和Herle（1997）粒間應變張量的概念引入到亞塑性模型中。
若希望學習亞塑性模型的更多內容，可以觀看我的可從。以下是不同版本的亞塑性模型對排水三軸壓縮試驗、不排水三軸壓縮試驗以及循環三軸試驗的模擬結果展示。左側是論文中的實驗及模擬結果，右側是我利用編寫的matlab程序得到的模擬結果作為對比。付費內容是我編寫的亞塑性模型matlab程序，umat程序以及參考文獻，上述各版本的亞塑性模型均有。需要的同學可以購買（可以私戳我通過ZFB購買，可優惠20%）。
圖1.Wu-Bauer亞塑性模型對排水三軸壓縮實驗模擬結果（密砂）
圖2. Wu-Bauer亞塑性模型對排水三軸壓縮實驗模擬結果（松砂）
圖3. Wu-Bauer亞塑性模型對排水三軸壓縮實驗模擬結果（松砂）
圖4. 粒間應變張量亞塑性模型對排水三軸壓縮實驗模擬結果
展開 熱-彈-黏塑性晶體塑性模型文章推薦
傳統室溫本構模型通常需要依賴大量不同溫度、不同加載路徑下的實驗數據進行擬合，很難真正解釋“溫度如何影響晶體滑移和多晶塑性響應”。
Cyr 等人針對這一問題提出了一個三維熱-彈-黏塑性晶體塑性模型，即 TEV 模型，用于描述 FCC 多晶材料，特別是 AA5754 鋁合金在升溫條件下的力學行為。該模型的核心思想是：材料變形不僅包含彈性變形和晶體塑性滑移，還需要顯式考慮熱膨脹變形。因此，總變形梯度被分解為彈性/剛體轉動部分、熱變形部分和塑性變形部分。
在本構層面，作者保留了 FCC 晶體的 12 個 {111}<110> 滑移系，并采用冪律型滑移率方程描述率相關塑性流動。與常規晶體塑性模型不同的是，該模型把溫度效應系統地引入到多個關鍵物理量中：首先，單晶彈性常數 C11、C12、C44 隨溫度變化；其次，滑移阻力引入熱軟化函數，用來描述溫度升高后滑移更容易發生的現象；再次，單滑移硬化參數也被寫成溫度函數，包括參考臨界分切應力、初始硬化率和硬化指數。
這個模型的優勢在于，它不是簡單地給宏觀應力-應變曲線加一個溫度修正系數，而是從晶體滑移層面描述溫度對材料響應的影響。換句話說，它可以同時分析宏觀應力變化、微觀滑移活動、織構演化、局部應變集中和熱軟化機制。因此，它比普通經驗型熱塑性模型更適合用于多晶材料溫成形模擬。
作者首先利用 AA5754 鋁合金在 25 ℃、148 ℃、204 ℃ 和 232 ℃ 下的單軸拉伸實驗數據標定溫度相關硬化參數。隨后，又預測了 177 ℃ 和 260 ℃ 下的拉伸響應。
展開 【免費】ABAQUS中纖維增強水泥基復合材料/混凝土/SHCC/ECC/FRC的損傷塑性模型
本貼介紹纖維增強混凝土（SHCC/ECC/FRC）的損傷塑性模型，分為四個部分，首先介紹真實應力和真實應變的轉換，然后介紹SHCC/ECC/FRC的應力-應變關系，再介紹SHCC/ECC/FRC的損傷塑性模型，最后進行四點彎曲梁的累加循環仿真計算。
1. ABAQUS中真實應力與真實應變
ABAQUS中必須用真實應力和真實應變來定義塑性。而大多數實驗（單軸拉伸、單軸壓縮等）得到的是名義應力和名義應變。故必須將實驗得到的名義應力和名義應變轉換為真實應力和真實應變，從而得到ABAQUS中需要的材料參數。
考慮塑性變形的不可壓縮性，真實應力與名義應力間的關系:
得到：
令
可以得到真實應力與名義應力、名義應變的關系：
真實應變與名義應變的關系：
2. SHCC/ECC/FRC的應力-應變關系
SHCC/ECC在單軸拉伸時呈現應變硬化現象，其簡化的本構關系較為常見的有兩種：理想彈塑性模型和線性強化彈塑性模型。在本帖中，采用線性強化彈塑性模型。
SHCC/ECC簡化的拉伸應力-應變曲線，包括3個階段:
SHCC/ECC簡化的壓縮應力-應變曲線
3.
展開 混凝土塑性損傷模型 ￥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>的曲線。
展開 
黃永剛晶體塑性模型耦合相場方法模擬多晶斷裂
模擬包含200個晶粒的多晶模型，使用平面應變簡化，拉伸變形30%，模擬拉伸過程中裂紋的產生和發展，其中斷裂總能量包含彈性變形能和塑性耗散功兩部分，模擬結果如下
初始多晶模型：
網格劃分（CPE4網格）：
相場分布（0：材料完好，1：材料完全失效）：
退化程度分布：
歐拉角（phi）分布：
可以看到耦合相場的晶體塑性模型具有潛在的預測裂紋萌生和發展的能力，其準確程度取決于斷裂能參數的選擇，與更精細的實驗對比，如原位的ebsd拉伸將成為良好的校核手段。這可能成為介觀尺度下斷裂力學的應用提供良好的參考
雙線性彈塑性模型(四)
/f358.png', dpi = 300) #保存圖片
plt.show()
得到的迭代路徑
★★★★ 往期相關 ★★★★
雙線性彈塑性模型(三)
雙線性彈塑性模型(二)
雙線性彈塑性模型(一)
ABAQUS中UMAT中的經典循環塑性本構模型及相應子程序代碼 ￥20
ABAQUS中UMAT中的循環塑性模型，包含非線性各向同性強化彈塑性、線性各向同性強化彈塑性、線性隨動強化彈塑性模型，包含CAE文件、UMAT文件等。
一個好用的Abaqus晶體塑性模型生成插件-Voronoi模型V8.0
V7.0版本介紹：
一個好用的Abaqus晶體塑性模型生成插件-Voronoi模型
https://zhuanlan.zhihu.com/p/611427546
2. V8.0版本新增功能：
2.1 二維核殼晶體模塊
圖2.1 二維圓形核殼晶體模塊
圖2.2 二維多邊形核殼晶體模塊
2.2 三維核殼晶體模塊
圖2.3 三維球形核殼晶體模塊
圖2.4 三維多面體核殼晶體模塊
2.3 桁架模型模塊
圖2.5 桁架結構模型生成模塊
2.4 圓形和圓柱邊界加權晶體模塊
圖2.6 二維圓形邊界加權晶體模塊
圖2.7 三維圓柱邊界加權晶體模塊
2.5 二維梯度晶體模塊
圖2.8 二維梯度晶體模塊
2.6 三維圓柱邊界梯度模塊
圖2.9 三維圓柱邊界梯度晶體模塊
展開 ABAQUS系列講座——混凝土塑性損傷模型
【視頻地址】
ABAQUS混凝土塑性損傷模型-視頻課程
【PPT預覽】
如何理解ANSYS彈塑性分析中的強化模型
昨天在整理文檔的時候，發現很早以前有朋友和我探討ANSYS中強化模型的意義問題，當時我先把問題存在有道云筆記里，待有空的時候琢磨琢磨，結果后來竟然給忘記了，實在是不靠譜啊！那么既然如此，今天就把這個問題重新拿出來，聊一聊，不足的地方，還望各位同行補充。
先來回顧一些概念
什么時候才需要做彈塑性分析呢？線彈性分析階段就是應力和應變成正比唄，即應力=應變*彈性模量，卸載以后一切恢復原狀。一旦在達到材料的彈性極限后，繼續加載，使材料進入塑性階段，此時再卸載就無法恢復原狀。
那么在這個過程當中，構件產生的總應變就可以分為彈性應變和塑性應變兩部分，彈性應變依然和應力存在正比的關系，關鍵就是如何建立起來塑性應變與由此產生的應力之間的關系呢？這就需要引入塑性模型( Plasticity Models)了。
影響塑性應變的因素有很多，如加載歷史(這就是為什么彈塑性分析要涉及到荷載步了)、溫度、應力、應變率，以及一些內部因素，如材料的屈服強度、損傷等。
那么，塑性模型如何來描述塑性發展的過程呢？ANSYS用三個準則來解決這個問題：
屈服準則：加載過程中，一旦材料的等效應力超過屈服應力，程序判定進入塑性狀態，這是解決一個從彈性到塑性的過渡點問題；
流動準則：當構件發生塑性應變時，流動準則定義了應變方向，也就是說，流動準則可以描述在達到屈服后，在每一個荷載增量的作用下，塑性應變的各個分量是如何發展的；
強化準則：描述了初始屈服準則隨著塑性應變的增加是怎樣發展的。
關于“強化”，得多說幾句，當材料經過屈服階段的塑性變形后，卸載，再加載到屈服，新的屈服點要比原屈服點高一些。那第一次屈服點就對應著“初始屈服準則”，每一次的屈服都比上一次高一點，這個發展的過程就是強化。
展開 求助：混凝土塑性損傷模型
求助：請問混凝土塑性損傷模型的損傷因子什么情況下可以不定義，什么情況下要定義？比如我現在單調受彎加載該構件，那該構件（塑性損傷模型）需不需要定義損傷因子？
ABAQUS三維hill48彈塑性模型VUmat子程序（彈性為正交各向異性） ￥388
1.ABAQUS三維hill48彈塑性模型VUmat子程序
2.彈性階段為正交各項異性材料
3.hill48和正交各項異性材料參數參考ABAQUS靜力模塊自帶的模型參數
4.發貨方式為百度網盤鏈接，包含子程序及上面跑的兩個模型相關文件，包含Cae，inp文件，odb文件等
5.ABAQUS版本為2024，低版本可以利用導入inp文件的方式運行及修改
6.可以免費答疑三次，后續添加你自己的模型或者相關參數等輔導都可以優惠。
基于HUANG晶體塑性模型的三點彎曲案例------案例十
?基于HUANG晶體塑性模型的三點彎曲案例
案例教學如下
1，建立幾何模型：0.6*0.3*0.05（mm）的平板，半徑0.05，長度0.3的下壓棍和支撐并完成裝配
三點彎曲模型圖
2，分配材料屬性：分別賦予板材晶體屬性參數，棍子純彈性參數，并假設下壓輥和支撐的剛度遠遠大于平板，即變形過程中約束為剛體。
模型的材料屬性分配
3，建立接觸條件：板材與棍子支撐之間法向硬接觸，切向摩擦系數設置為0.2
建立接觸屬性
4，網格劃分對板材和棍子，支撐進行網格劃分，板材網格類型C3D4（保留晶界形狀），棍子支撐C3D8
5，載荷施加：固定下邊支撐，同時對上面的下壓輥施加向下的0.04的位移，其他自由度全部限制為0。
6，提交作業與后處理
三點彎曲位移分布情況
累計塑性應變分布情況
模型的應力分布情況
模型的變形流動方向分布情況
展開 考慮晶界影響的晶體塑性模型
這個思想非常適合和晶體塑性模型結合，因為晶體塑性本來就是逐滑移系計算的。
在作者的模型中，晶界通透性可以進一步轉化為晶界障礙應力。通透性越高，障礙應力越小；通透性越低，障礙應力越大。這樣一來，晶界對塑性滑移的影響就可以直接進入滑移率方程：只有當有效分切應力足夠克服晶界障礙時，晶界附近的滑移才能繼續發展。
這篇文章的另一個重要部分是位錯重分配。作者把有限元單元中的位錯內容等效成“超位錯”，再根據塑性滑移活動對可動位錯進行重新分布，并計算由這些位錯分布產生的背應力。這個處理很有啟發性：它不是直接追蹤每一根真實位錯，那樣計算量太大；但它也不是完全經驗化地加一個強化項，而是在連續體模擬和位錯物理之間做了一個折中。作者的模型概念圖：
積分流程圖：
從結果來看，作者的模型能夠再現單個位錯塞積問題中的位錯密度分布；在二維和三維規則晶粒陣列中，也能預測出與實驗同量級的 Hall–Petch 斜率。對于粗晶鐵多晶拉伸響應，這個兩尺度模型比傳統 CP-FEM 或 Taylor 類模型給出了更好的預測。此外，作者還比較了 Fe-3%Si 柱狀晶樣品中的晶格曲率，模型預測的晶界附近曲率峰值與實驗結果基本一致。
作者發現模型可以非常準確的預測晶粒尺寸效應：
我認為這篇文章的價值不只是“提出了一個更復雜的模型”，而是提供了一種很清楚的建模思路：晶界強化不一定只能通過經驗晶粒尺寸項來描述，也可以從滑移傳遞、位錯通量和局部障礙應力出發，逐步把晶界的物理作用放進晶體塑性框架中。
我們可以作者提出的模型完整的構建一個考慮晶界多尺度模型，演示如何計算每個滑移系對應的晶界通透系數，并將其轉化為晶界障礙應力引入晶體塑性本構中。通過對比是否考慮晶界障礙項，可以觀察晶界附近滑移活動、位錯密度分布以及應變局部化特征的變化。
展開 晶體塑性：構建Dream3D pipeline用于將EBSD模型制作成Abaqus可執行文件 ￥180
晶體塑性：構建Dream3D pipeline用于將EBSD模型制作成Abaqus可執行文件
案例實操
用于生成模型的Dream3D pipeline文件，只需要你設置EBSD數據的路徑和導出路徑即可，可以直接生成abaqus的晶體塑性模型，提供原始文件！
包含老版本Dream3D 6.5的管道文件，并且根據官方的使用說明文件已經成功移植到最新版Dream3D 7.4版本了。
