abaqus應變含義的案例
ABAQUS網格劃分顏色的含義
作為初學ABAQUS的小白，對于ABAQUS使用各個過程中的部件顯示的顏色變化感到好奇。
在創建部件部件中為灰色；劃分網格時，綠色代表可劃分結構性網格，黃色代表可劃分掃掠網格，粉色代表可劃分自由網格；裝配時部件變為藍色。
ABAQUS中關于幅值的含義
幅值*位移值=真實位移量幅值*應力值=真實應力值
換句話說，幅值就是一個比例系數
Abaqus后處理-云圖變量含義（部分）
結果中分量說明：
S11、S22、S33指各軸正應力；
S12指作用于XZ平面(與“2”，即Y軸垂直的平面)內，沿1方向剪應力；
S13指作用于XY平面內，沿1方向剪應力；
S23指作用于XY平面內，沿2方向剪應力。
若為柱坐標，S12、S13、S23分別指：由徑向向環向的剪應力、由徑向向軸向的剪應力、由環向向軸向的剪應力。
ABAQUS 二次開發VUMAT中的NBLOCK的含義
1、幫助文件的解釋
Number of material points to be processed in this call to VUMAT.
3、nblock的含義：
Abaqus一次調入vumat進行計算的單元數，是個常數，最大為136。假設模型單元總數為n，令n=a×136+b，其中a整數，b是小于136的整數。則前a×136個單元的nblock值為136，后b個單元的nblock值為b。
4、Abaqus調用vumat的過程
計算一步的過程中每nblock個單元調入一次vumat，即n（n=a×136+b）個單元的計算模型計算一步需要調入a+1次vumat。
理由如下：1250個單元的計算過程計算了1516步，一次調用vumat在data.txt文件中生成1250個隨機數，計算一步調用vumat的次數為Int（1250/136）+1=10。所以data.txt文件中生成的總個數應為1250×10×1516=1.895E7個。打開data.txt文件發現數據總個數為2.33E7，略大于1.895E7，這是因為還有一些計算步被終止，并未顯示。
原作者附帶的案例說明文檔一并轉載到此
最近的一點收獲以前總沒確切的理解VUMAT里面Nblock的含義，只是簡單理解為計算模型總的高斯點數。
nblock .rar
展開 
Abaqus數據文件Data File中最后Job Time Summary各項含義 ￥2
</p><p>所以在Abaqus中怎么去查看計算時長呢？</p><p>（這里是看計算完成后實際用時，而不是提交計算時預計用時）</p><p><br></p>
abaqus子程序VUSDFLD——考慮應變率與應變軟化效應的軟土模型 ￥25
<p><strong>【注意】本貼子只包含子程序文件</strong></p><p>基于<a href="https://www.yqgqt.org.cn/qa/6302" rel="noopener noreferrer" target="_blank">abaqus子程序</a>VUSDFLD編寫的由Einav與Randolph提出的西澳模型，用于求解軟黏土體劇烈變形后的強度變化，可應用于的大變形計算。</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 UEL中UMAT材料參數設置問題，等效迭代數（Equll Iter）的含義
問題一：在使用UEL的時候，為了顯示自定義單元，通常我們是設置一層實體單元賦予UMAT材料來實現。在UMAT設置中，一般來說是設置楊氏模量為1E-11，泊松比為0.3。請問一下為什么要這樣設置呢？既然只是借用實體單元的積分點來傳遞 UEL計算得到的SDV，通過一個common來傳遞，那這個參數為什么要設置這么小？
我最近做了一些工作，發現如果把楊氏模量設置為實際的材料參數，得到的結果會和實際情況相比偏小，設置為1E-11的話，結果就基本一致。
問題二：implicit算法是需要設置時間步長，一般有初始值，最大值，最小值；也可以設置為固定步長。我最近的工作設置為了固定步長，在作業監控器界面出現了下圖的情況，沒有不連續的迭代，但是等效迭代欄（Equll Iter）出現了大于1的情況，請問一下大家，這個數代表了什么？我的理解是第一次嘗試計算不收斂，然后嘗試第二次，但是我設置了固定步長，那么第二次為什么會出現收斂？
展開 有限元理論基礎及Abaqus內部實現方式研究系列17： 幾何非線性的物理含義
配合本章的視頻解說和操作演示可看下方：
http://www.yqgqt.org.cn/college/video/c12884 20理論系列文章17-幾何非線性的物理含義
2.1 幾何非線性的物理含義
2.1.1 從線性到幾何非線性
一個物體從初始狀態A由于受到外部載荷運動，如果現在已知了另一種狀態B的位移，那么其它的任意狀態C的位移怎么求？
如果能直接從B的位移乘以一個常量就得到C，那么這個系統就是線性系統。譬如下面的800mm的懸臂梁問題，在Abaqus中用線性計算，載荷F和位移u是直線關系。
載荷1N的時候Abaqus計算得到最大位移時1.177mm。
那么載荷是1000N時是多少？顯然，不用計算也知道就是1177mm。 但1177這個值明顯有問題，已經超過了梁的長度，按生活經驗判斷這個梁估計都斷了或者極端扭曲了，所以這種情況需要用幾何非線性來計算。
2.1.2 幾何非線性簡單物理含義
虛功原理如下：
在物理上可解釋能量守恒原理，即在某一個時刻點，假定在外力作用下有個虛擬的位移，那么外力在虛擬位移下做的虛功=內部應變能的變化相同。
為了更好的理解上面的物理解釋，如果我們把當成真實的位移，那么外層加上對時間的積分，可以理解為外力在虛擬位移下做的虛功=內部應變能的一段小時間內對應變能的積分：
那么就和我們高中時學過的小球受重力作用后勢能、動能相互轉換是一樣的原理。
舉個簡單的有限元例子，譬如下方一個減縮積分S4R單元的右端受X方向兩個力F：
得到的位移為U（相當于虛位移），那么力F做的功是W=2F*U。
增加的應變能為S*E*V在對所有時刻點t的積分，S、E、V都是當前時刻的值。
展開 ABAQUS中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
有限元理論基礎及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平均應力和應變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應力及應變進行自動提取并計算平均值
能夠得到每一幀的應力和應變平均值，并保存到CSV文件中
所得到的應力包括S11,S22,S33,S12,S13,S23以及Mises七個應力平均值，以及E11,E22,E33,E12,E13,E23六個應變平均值
ABAQUS提取單元平均應力/應變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據需要進行修改。
ABAQUS批量提交Job與Python讀取ODB結果應力應變數據
批量提交的核心是需要等待當前計算任務的結束，上圖中若干個job的提交代碼如下：
# coding: utf-8
#微信公眾號：ABAQUS二次開發
#作者：阿信老師CAE
#email：axin_cae@163.com
#2022.03.17
from abaqus import *
from abaqusConstants import *
jobList = mdb.jobs.keys() #獲取所有計算任務的jobName
for jobName in jobList:
mdb.jobs[jobName].submit(consistencyChecking=OFF) #提交計算
mdb.jobs[jobName].waitForCompletion() #等待計算完成
print jobName , "is completed"
不過在實際的操作過程中，我們可能會需要避開一些job不提交，或者job窗口事先并沒有建立job，等等，總之實際問題永遠比任何教程都復雜，不過只需要靈活面對就可以了，處理起來并不難。
展開 Abaqus插件——平均應力應變提取 ￥60
通過該插件可實現:
1）提取所有幀的任意單元集合的平均應力（事先定義單元集合，如圖中的SET-1）
2）提取所有幀的任意區域的x、y、z方向的平均應變（事先定義節點集合，如圖中的SET-2）
3）將以上數據保存至excel文件（excel文件名為odb文件名稱+_Stress_Strain.csv）
*************************注意事項******************************
1、插件使用過程中，如有任何問題請發郵件至shenz1hao@126.com
2、插件僅做學習交流使用，尊重原創者，切勿以營利目的傳播
*****************************************************************
********************插件安裝及使用*******************************
1、電腦路徑下輸入 %homepath%\abaqus_plugins并回車
2、將Stress-Strain文件夾解壓至當前目錄下
3、打開abaqus，菜單欄中點擊plug-ins，里面找出Stress-Strain
4、輸入相應參數（hx、hy、hz表示x、y、z方向模型長度，當以上三參數取1時輸出的為該方向位移）
*****************************************************************
展開 ABAQUS中對應力、應變的部分理解
對應力的部分理解
對應變的部分理解
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
