abaqus計算應(yīng)變
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時間：2023-02-27
abaqus計算應(yīng)變的視頻教程
ABAQUS：TNT在鋼管內(nèi)部爆炸，考慮應(yīng)變率
采用abaqus軟件CEL方法模擬TNT在鋼管內(nèi)部爆炸，鋼材在material中采用柔性與剪切損傷并添加應(yīng)變率來考慮在極端大變形狀態(tài)下材料的損傷與破壞
￥100 25分鐘
查看
abaqus計算應(yīng)變的實例教程
（1）顯示應(yīng)變：Abaqus計算完畢后得到導(dǎo)入結(jié)果，在后處理中查看，應(yīng)變E11=8.528e-1,E22=-5.173e-1如下：
（2）計算應(yīng)變：Abaqus中采用UMAT子程序，利用我們的子程序調(diào)試插件DUS調(diào)試UMAT，在Visual Studio中查看dStran的值，發(fā)現(xiàn)在計算完應(yīng)變后，進(jìn)入UMAT時，E11=8.528e-1,E22=-5.173e-1，調(diào)試如下：
可以發(fā)現(xiàn)殼單元Abaqus的計算應(yīng)變和顯示應(yīng)變一樣，猜測都是對數(shù)應(yīng)變。
1.5.3 iSolver的應(yīng)變
iSolver中采用自帶材料進(jìn)行計算，材料參數(shù)和UMAT的輸入完全一致。
為了計算和Abaqus完全一致，iSolver也采用對數(shù)應(yīng)變計算方式，得到的應(yīng)變顯示如下，可發(fā)現(xiàn)和Abaqus完全一致。
==總結(jié)==
由上可以看到，在實際計算中，對體單元，Abaqus和iSolver都采用變形率積分方式來計算應(yīng)變，對殼單元，Abaqus和iSolver都采用對數(shù)應(yīng)變。一般理論書都認(rèn)為Abaqus是因為對數(shù)應(yīng)變計算復(fù)雜才采用別的應(yīng)變，但個人認(rèn)為應(yīng)該不是這個原因，因為Abaqus對體單元為了顯示對數(shù)應(yīng)變，依然重新計算了一遍，說明Abaqus體單元采用變形率是有其它原因的，具體什么原因我也沒研究清楚，歡迎探討。
如果有任何其它疑問或者項目合作意向，也歡迎聯(lián)系我們：
snowwave02 From www.yqgqt.org.cn
email: snowwave02@qq.com
以往的系列文章：
1.7.1 ========第一階段========
第一篇：S4殼單元剛度矩陣研究。
http://www.yqgqt.org.cn/content/post/338859
第二篇：S4殼單元質(zhì)量矩陣研究。
展開 Graphite/epoxy材料參數(shù)：
E1=181Gpa E2=10.3Gpa G12=7.17Gpa Nu12=0.28
例題中的剛度陣為：
插件輸出的剛度陣：
例題中的中面應(yīng)變：
插件輸出的中面應(yīng)變：
全局應(yīng)變：
全局應(yīng)力：
局部應(yīng)變：
局部應(yīng)力：
插件輸出的鋪層的全局和局部的應(yīng)力應(yīng)變：
對比上面的結(jié)果可以發(fā)現(xiàn)插件的計算結(jié)果和例題的計算結(jié)果是相吻合的。
設(shè)t時刻的作用于試樣標(biāo)距段加載方向的載荷為F，工程應(yīng)力為，工程應(yīng)變為，真實應(yīng)力為，真實應(yīng)變為。規(guī)定試樣受拉伸長時，載荷、應(yīng)力及應(yīng)變取正值，反之，受壓縮短時，載荷、應(yīng)力及應(yīng)變取負(fù)值，則
式1
式2
式中表示0時刻至t時刻試樣長度的增量。
式3
式4
真實應(yīng)力 的求解變換中利用了材料變形過程中體積不變的假設(shè)，即。
繪制應(yīng)力-應(yīng)變曲線時，往往不管拉伸或者壓縮，都將應(yīng)力和應(yīng)變繪制成正值。這樣，拉伸時，按照公式<1>至<4>計算出的應(yīng)力和應(yīng)變均為正值，不需要進(jìn)行變換；壓縮時，按照上述方法計算出的應(yīng)力和應(yīng)度均為負(fù)值，需要進(jìn)行變換。變換方式為：對上述公式中所有的應(yīng)力和應(yīng)變乘-1。按照這種規(guī)則，壓縮時應(yīng)力、應(yīng)變用<5>至公式<8>進(jìn)行汁算。其中公式<5>中F取負(fù)值。
式5
式6
式7
式8
?
?
?
展開 <p><strong>【注意】本貼子只包含子程序文件</strong></p><p>基于<a href="https://www.yqgqt.org.cn/qa/6302" rel="noopener noreferrer" target="_blank">abaqus子程序</a>VUSDFLD編寫的由Einav與Randolph提出的西澳模型，用于求解軟黏土體劇烈變形后的強(qiáng)度變化，可應(yīng)用于的大變形計算。</p><div contenteditable="false" width="100%"><img src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png" title="8$U(VZ82]O{OEMQB}[P(ZMB.png" alt="8$U(VZ82]O{OEMQB}[P(ZMB.png" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/upload/202306/f69f50d42a81489ea1cb5e7a03da5c14.png?
展開 
abaqus計算應(yīng)變的相關(guān)專題、標(biāo)簽、搜索
abaqus計算應(yīng)變的最新內(nèi)容
經(jīng)建模驗證過的，考慮混凝土應(yīng)變率效應(yīng)的混凝土本構(gòu) 想要交流可以?v：wangh2444
有限元后處理直接與數(shù)據(jù)圖片處理、論文撰寫相關(guān)，除了典型的應(yīng)力張量與應(yīng)變張量外，ABAQUS還提供了大量可供使用者讀取的其他應(yīng)力/應(yīng)變/損傷參數(shù)，這都有助于結(jié)果的分析。今天喵星人就教你讀懂其中的應(yīng)力、應(yīng)變及損傷的后處理細(xì)節(jié)。
一、應(yīng)力相關(guān)
根據(jù)用戶手冊及后處理分類，ABAQUS提供了三類典型的后處理變量：
1.不變量
不變量的定義是指張量在坐標(biāo)旋轉(zhuǎn)下保持不變的量。這些量反映了材料內(nèi)在的力學(xué)狀態(tài)
一套基于 MATLAB/Fortran 編寫的二維鍵基近場動力學(xué)（Bond-based Peridynamics）數(shù)值仿真代碼。程序采用經(jīng)典的動態(tài)松弛算法（Dynamic Relaxation），將動力學(xué)方程轉(zhuǎn)化為解決準(zhǔn)靜態(tài)問題的工具，模擬二維材料在單軸壓縮載荷下的響應(yīng)及裂紋擴(kuò)展過程。
準(zhǔn)靜態(tài)模擬方案：利用動態(tài)松弛代碼，通過人為阻尼迭代，穩(wěn)定求解準(zhǔn)靜態(tài)單軸壓縮過程。
<p>因為要仿真混凝土破壞實驗，考慮用abaqus里面的CDP模型，查閱了相關(guān)資料進(jìn)行了理論總結(jié)，并根據(jù)理論編寫計算程序。</p><p>ABAQUS中CDP 模型中采用的是混凝土在單軸受力狀態(tài)下的應(yīng)力和非彈性應(yīng)變，非彈性應(yīng)變根據(jù)混凝土的單軸應(yīng)力-應(yīng)變曲線換算。</p><p>根據(jù)GB50010-2010混凝土結(jié)構(gòu)設(shè)計規(guī)范，混凝土單軸應(yīng)力應(yīng)變關(guān)系如圖：</p><p><img src="https://img.jishulink.com
利用python腳本對ODB文件中單元集里所有積分點的應(yīng)力及應(yīng)變進(jìn)行自動提取并計算平均值
能夠得到每一幀的應(yīng)力和應(yīng)變平均值，并保存到CSV文件中
所得到的應(yīng)力包括S11,S22,S33,S12,S13,S23以及Mises七個應(yīng)力平均值，以及E11,E22,E33,E12,E13,E23六個應(yīng)變平均值
<p class="ql-align-justify">本內(nèi)容基于韓林海的約束混凝土模型所制作的Excel，可用于將其輸入直接到ABAQUS中，用于建立鋼管約束混凝土型，具體如下：</p><p class="ql-align-justify">模型介紹：</p><p class="ql-align-justify">本模型基于<span style="color: rgb(25, 27, 31);"
文章題目：《Strain rate effect of high purity aluminum single
crystals: Experiments and simulations》
文章doi：10.1016/j.ijplas.2014.10.002
推薦理由：作者研究了高純鋁不同應(yīng)變率下單晶塑性變形的取向依賴性，不同應(yīng)變率下的流動應(yīng)力情況通過Laue Back-Reflection
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應(yīng)力和平均應(yīng)變。
3、應(yīng)用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應(yīng)力和平均真應(yīng)變，可通過對 RVE 內(nèi)每一個單元的真應(yīng)力 （真應(yīng)變）取均值獲得。使用一階均勻化計算方法輸出的應(yīng)力和應(yīng)變適用于各種邊界條件，但需要對每個單元進(jìn)行應(yīng)力（應(yīng)變）的輸出和計算。
聲明：貼主目前正在學(xué)習(xí)ABAQUS，對UMAT有一點淺淺的了解，若有不對的地方，請理性留言討論。
貼主的ABAQUS模型即使使用工作站，一運行也好幾天，苦惱不已，因此萌生了探討影響計算速度的相關(guān)因素的想法。
首先影響ABAQUS運行速度的最主要因素是模型的復(fù)雜程度，但往往模型是不易更改的，因此本文不做討論，而著重討論容易更改的部分，進(jìn)而提高ABAQUS的運行效率。以下對計算效率的討論均使用了使用
最近在開展分析時遇到錯誤如下：MAXIMUM SIZE OF STATIC WORKSPACE HAS BEEN EXCEEDED. CURRENT WORKSPACE SIZE IS 16384.00 MB. THE SIZE OF THE WORKSPACE CAN BE INCREASED USING THE SYSTEM ENVIRONMENT VARIABLE ABA_SINT_CAP.
