abaqus計算塑性應變的案例
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>所開發的約束混凝土應力應變關系模型，以及損傷因子，其中受壓本構以及受拉本構以及其損傷因子均有，且附帶鋼材料的二次流塑模型，可直接輸入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混凝土損傷塑性模型-C30EXCEL計算表格（含計算結果）
這是根據GB50010-2010中混凝土結構設計規范中的混凝土本構模型，結合文獻所述的損傷因子定義，編制的計算C30混凝土非彈性應變和損傷因子的EXCEL表格。也是邊學變做，希望能和大家多交流。
C30砼本構（損傷塑性模型）.rar
ABAQUS混凝土損傷塑性損傷因子計算依據.rar
基于Abaqus的高純鋁不同應變率下單晶塑性變形的取向依賴性研究
文章題目：《Strain rate effect of high purity aluminum single
crystals: Experiments and simulations》
文章doi：10.1016/j.ijplas.2014.10.002
推薦理由：作者研究了高純鋁不同應變率下單晶塑性變形的取向依賴性，不同應變率下的流動應力情況通過Laue Back-Reflection 技術測量，并提出了兩類單晶本構模型用于預測單晶不同應變率的應力響應的能力，研究表明，相較于傳統的單晶冪律流動模型，所提出的另外的唯象和位錯密度模型很好捕捉了應變率效應，提出的唯象模型參數少，便于擬合，物理模型參數更多，但物理意義更明確，這在捕捉單晶多滑移系開動時提供了更準確的預測（更接近實驗結果）。
展開 有限元理論基礎及Abaqus內部實現方式研究系列27： Abaqus內部計算和顯示的應變
（1）顯示應變：Abaqus計算完畢后得到導入結果，在后處理中查看，應變E11=8.528e-1,E22=-5.173e-1如下：
（2）計算應變：Abaqus中采用UMAT子程序，利用我們的子程序調試插件DUS調試UMAT，在Visual Studio中查看dStran的值，發現在計算完應變后，進入UMAT時，E11=8.528e-1,E22=-5.173e-1，調試如下：
可以發現殼單元Abaqus的計算應變和顯示應變一樣，猜測都是對數應變。
1.5.3 iSolver的應變
iSolver中采用自帶材料進行計算，材料參數和UMAT的輸入完全一致。
為了計算和Abaqus完全一致，iSolver也采用對數應變計算方式，得到的應變顯示如下，可發現和Abaqus完全一致。
==總結==
由上可以看到，在實際計算中，對體單元，Abaqus和iSolver都采用變形率積分方式來計算應變，對殼單元，Abaqus和iSolver都采用對數應變。一般理論書都認為Abaqus是因為對數應變計算復雜才采用別的應變，但個人認為應該不是這個原因，因為Abaqus對體單元為了顯示對數應變，依然重新計算了一遍，說明Abaqus體單元采用變形率是有其它原因的，具體什么原因我也沒研究清楚，歡迎探討。
如果有任何其它疑問或者項目合作意向，也歡迎聯系我們：
snowwave02 From www.yqgqt.org.cn
email: snowwave02@qq.com
以往的系列文章：
1.7.1 ========第一階段========
第一篇：S4殼單元剛度矩陣研究。
http://www.yqgqt.org.cn/content/post/338859
第二篇：S4殼單元質量矩陣研究。
展開 
ABAQUS混凝土損傷塑性模型-C30EXCEL計算表格（含計算結果） ￥3.7
計算表格（如下）中標黃部分的參數可自行設定后，EXCEL程序會自動計算“抗拉強度、非彈性應變、受拉損傷因子”。
計算層合板的ABD陣/柔度陣/應力應變---abaqus插件 ￥29.9
Graphite/epoxy材料參數：
E1=181Gpa E2=10.3Gpa G12=7.17Gpa Nu12=0.28
例題中的剛度陣為：
插件輸出的剛度陣：
例題中的中面應變：
插件輸出的中面應變：
全局應變：
全局應力：
局部應變：
局部應力：
插件輸出的鋪層的全局和局部的應力應變：
對比上面的結果可以發現插件的計算結果和例題的計算結果是相吻合的。
線性強化彈塑性umat子程序-umat在abaqus計算流程中的意義及調試方法
abaqus子程序用fortran語言編寫，且子程序又通過abaqus調用，因而涉及到很多的力學及有限元理論知識，使得調試難度陡然增加。那么abaqus子程序常用的調試方法是什么呢？其實說出來也很簡單，就是打印變量法，即把程序中的變量打印出來（一般打印結果在當前工作目錄下的job-name.log的文件里），然后考察其合理性。
如在上述umat子程序中定義了一個save在內存的變量，用來記錄進入umat的次數。
直接打印變量的方法很有用，但是當單元數增加后，眾多的打印信息形成很多的干擾因素，這時候我們需要控制打印的頻率，常用的方法為針對某個單元的某個積分點打印，如下面的一段程序所示，紅色圓圈里內容的含義是當單元編號為1且積分點編號為1時才打印相關信息，這樣調試更具針對性。
我們通過打印信息發現，一次迭代二次進入umat，第一次進入umat僅是給abaqus返回雅可比矩陣，第二次進入umat的目的是為了更新應力等信息。
視頻教程有這個帖子的更詳細解釋，感興趣的點擊下面的鏈接觀看
http://www.yqgqt.org.cn/college/video/c14014
展開 ABAQUS混凝土損傷塑性模型-2010混凝土結構設計規范中C50混凝土-彈模34400Mpa-損傷因子計算及EXCEL
這是我自己計算的2010規范中ABAQUS混凝土損傷塑性模型-2010混凝土結構設計規范中C50混凝土-彈模34400Mpa-損傷因子計算及EXCEL
首先用自己的數據計算2010規范中規定的混凝土本構關系
然后借助文件夾中02版規范的方法，計算損傷因子。
以后還會有詳細計算方法，此數據僅供參考。
2010規范用C50混凝土損傷塑性本構關系數據-彈模34400MPa-帶損傷因子-自己數據計算得出.rar
