abaqus應變設置的案例
如何利用ABAQUS軟件在CAE界面中完成應變軟化子程序的設置？ ￥5
最近在ABAQUS中開展了CEL大變形分析，其中涉及到應變軟化子程序的嵌入，特此將最近的學習心得和各位分享一下，為大家避坑。
此文檔為VUSDFLD子程序如何在CAE中激活的步驟詳解，希望可以為有需要的朋友帶來幫助！如果有不正之處也請大家批評指正（新手小白的瑟瑟發抖）。
發現了一些問題，請查看最新版的文件??！
【Ls-dyna】Hypermesh&Ls-dyna聯合仿真時如何設置輸出單元應變和查看單元應變？
通常，使用Ls-dyna進行瞬態分析時，默認的輸出控制中是不包含單元的應變值的。但是，如果關心單元上的應變，或者需要對單元應變有更全面的了解，那么就需要在建模和分析中進行輸出控制，使計算結果中包含單元的應變。那么，問題來了，如何在Hypermesh中如何設置可以輸出單元的應變？如何在LS-POST中顯示單元的應變云圖？ 接下來，依次解答。
如何在Hypermesh中如何設置可以輸出單元的應變？
這里需要首先說明的是，本人使用Ls-dyna以來，一直使用的前處理不是Ls-prepost前處理，而是前處理軟件Hypermesh，后處理則使用Ls-prepost。所以，這里只說明Hypermesh&Ls-dyna聯合仿真時如何在Hypermesh中設置控制應變輸出。
前處理軟件工作界面
后處理軟件工作界面
首先啟動Hypermesh，點擊“user Profiles”，選擇Ls-dyna，其他默認，進入Ls-dyna分析模塊。
然后，在軟件的面板區域選擇“analysis”，并點擊“control card”。
點擊面板區域的Next，一直到出現“database-extent-binary”，點擊該按鈕
在關鍵字*Database-Extent-Binary的定義中，將第一行第四個參數【STRFLAG】的值設置為1，表示在二進制結果文件d3plot中輸出單元的應變。
在Hypermesh中按照上述方式就可以在d3plot中輸出單元的應變，由于d3plot是二進制文件，所以只能借助后處理LS-Prepost查看應變結果。
那么，如何在LS-POST查看應變，顯示應變云圖呢？
打開LS-PrePost，點擊【Fcomp】，選擇【strain】，就可以觀察計算后的應變云圖。
展開 應變壽命疲勞分析理論分析基礎及DesignLif參數設置 ￥6
? Strain-Life (EN) 應變疲勞分析理論基礎
? 討論循環應力-應變曲線和應變-壽命關系的關系
? 討論平均應力的影響
應變疲勞壽命分析理論基礎
? 應變壽命疲勞（EN）使用循環應變反轉和應變壽命關系方程評估疲勞損傷
–局部塑性應變導致疲勞
–適用于低周期和高周期應用
? 應力小于或大于屈服
–使用彈塑性應變
? 直接計算或根據彈性計算進行調整
? 相對較新的疲勞分析技術
–大約30年前開始使用
–難以手動計算
?僅限于CAE應用程序
展開 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?
展開 
LS-DYNA中的操作及設置（五）（應變率，質量縮放）
本文翻譯自官方文檔，原文鏈接：
https://www.dynasupport.com/howtos/general
一、應變率（Strain rate）
應變率效應不僅與加載速率有關，還與試件的尺寸、形狀有關。在單軸拉伸試驗中，假如試件發生均勻變形，也就是沒有頸縮等局部化現象，那么應變率在試件中的分布是均勻的，此時有：
長度的變化為deltaL = r * time
工程應變為deltaL/L = r * time/L
工程應變率為strain per time = r/L
真實應變為ln(1+ engineering strain) = ln(1+ r*time/L)
真實應變率隨時間的導數為d(true strain)/dt = [ln(1+r*time2/L) - ln(1+r*time1/L)]/(time2-time1)
其中，L為試件長度；r為加載速率。
當然,事實上試件中的應變率并不是均勻分布的，所以我們需要在分析中給定一個應變率的變化范圍。為了估算應變率，我們可以針對有代表性的單元進行高精度的預分析，并輸出應變率（set STRFLG=1 in *DATABASE_EXTENT_BINARY）。還可以使用*DATABASE_BINARY_D3THDT 和 *DATABASE_HISTORY_SHELL這兩個關鍵字來輔助完成這一目標。
展開 LS-DYNA中的操作及設置（二）（有效塑性應變，環境變量，狀態方程）
有效應變與有效塑性應變是不同的，它的張量表示是（ p. 461 of LS-DYNA Theory Manual 2006）：
sqrt(2/3(eps)ij*(eps)ij)
LS-DYNA中還可以顯示其他種類的應變，都是通過節點位移計算而來，例如：
FCOMP > Infin
FCOMP > Green
FCOMP > Almansi
二、環境變量（Environment variables）
**僅適用于UNIX系統**
940.1版本的LS-DYNA中引入了如下環境變量：
LSTC_FILE：用來定義許可證文件
默認文件目錄為：/usr/local/lstc/LSTC_FILE，可以使用setenv LSTC_FILE (license file name) 命令來指定文件名。
LSTC_SECURE：用來定義許可證文件的格式
目前支持三種文件格式：old，new和eta。若不設置這一環境變量，程序將全部檢查這三種格式。
展開 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提取單元平均應力/應變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據需要進行修改。
ABAQUS中對應力、應變的部分理解
對應力的部分理解
對應變的部分理解
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
Abaqus平均應力和應變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應力及應變進行自動提取并計算平均值
能夠得到每一幀的應力和應變平均值，并保存到CSV文件中
所得到的應力包括S11,S22,S33,S12,S13,S23以及Mises七個應力平均值，以及E11,E22,E33,E12,E13,E23六個應變平均值
ABAQUS變量解讀：教你讀懂應力/應變/損傷
主應變
與主應力類似，ABAQUS也提供主應變輸出：
Max/Mid/Min Principal Strain：第一、二、三主應變，分別對應最大、中間、最小主應變，在判斷第二強度理論時有奇效。
In-Plane Principal Strain：平面問題最大/最小主應變。
Max Principal(abs)：絕對值最大主應變。
3.應變張量
與應力張量方向類似，其中需要同學們注意的是：
E適用于幾何線性分析
LE為對數應變，適用于大變形分析（開啟幾何非線性）
PE為塑性應變張量，用于描述不可恢復的變形
三、損傷相關
損傷在ABAQUS中應用廣泛，尤其是材料失效分析中。
1. 混凝土損傷
這是大家喜聞樂見的損傷變量，有兩類：
DAMAGEC（dc）：壓縮損傷變量，從0到1，1表示完全損傷。主要用來判斷壓潰區域與剪壓開裂區域。
DAMAGET（dt）：拉伸損傷變量，同樣從0到1。主要用來判斷受拉開裂區域，如下圖。
2. 鋼材損傷
SDEG：剛度退化標量，也可用于混凝土。表示材料剛度的折減程度。
3. 內聚力模型損傷
CSDMG：描述cohesive單元進入軟化段后的損傷狀態。
4. 復合材料損傷
ABAQUS支持多種復合材料損傷變量：
DAMAGEFT/FC：用戶手冊中描述為：
Fiber tensile/ compressive damage variable.
這里毫無疑問表征了纖維縱向的拉伸/壓縮損傷，如下圖。
DAMAGEMT/MC：用戶手冊中描述為：
Matrix tensile/ compressive damage variable.
展開 ABAQUS學習筆記—對應力應變的部分理解
之前關于后處理的一些文章，由于一些原因全部刪除，故今天開始重新開始分享一些關于ABAQUS的一些知識，希望能夠對大家有所幫助，也希望大家能夠繼續支持筆者。
那么今天，我們再對‘’ABAQUS中應力應變的部分理解的‘’內容進行講解。
在ABAQUS中，一般是把X軸當成1軸，Y軸當成2軸，Z軸當成3軸，那么：
S11就是X軸向的應力，正值為拉應力，負值為壓應力；
S22就是Y軸向的應力，正值為拉應力，負值為壓應力；
S33就是Z軸向的應力，正值為拉應力，負值為壓應力；
S12就是在YZ平面上，沿Y向的剪力；
S13就是在YZ平面上，沿Z向的剪力；
S23就是在XZ平面上，沿Z向的剪力；
由于剪力的對稱性:S11=S21,S13=S31,S23=S32；
由以上可知，S11，S22，S33為主應力；S12，S13，S23為切應力；
主應力分別以σ1，σ2，σ3表示，按數值排序為：σ1≥σ2≥σ3。在ABAQUS中分別對應為：Max.principal；Mid.principal；min.principal。這三個量在任何坐標下都是不變量。
我們可利于最大應力判斷一些情況：比如最大主應力（拉應力）大于混凝土的抗拉強度，則認為混凝土開裂；通過顯示最大主應力的法線方向，則可大致表示出裂縫的發展影響。
應變中一些符號的含義
E——總應變
EP——主應變
EE——彈性形變
PE——塑性應變分量
Eij——應變分量
Ie——非彈性應變分量
PEEQ——等效塑性應變。若該值大于0，則認為已經屈服
注：在ABAQUS后處理中，盡量不要看Mises，其表示平均應力，更適合金屬材料；對于鋼筋混凝土結構，我們要看其單軸拉伸方向上的應力和對應的應變。
展開 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時輸出的為該方向位移）
*****************************************************************
展開 