abaqus指定應(yīng)變的案例
批量提取Abaqus指定Step-x下 Set節(jié)點集的反力RF（Reaction force） ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發(fā)來批量提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中指定Step-x下的Set節(jié)點集的反力RF（Reaction force）。通過詳細(xì)的步驟說明、代碼示例和圖片展示，您將學(xué)會如何使用該腳本，自動化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，RF1（X），<span style="color: rgb(25, 27, 31);">RF2（Y），RF3（Z）</span>）。</p><p><br></p><p>如果還需要提取<strong><em><u>節(jié)點變形量 </u></em></strong>或者<strong><em><u>節(jié)點坐標(biāo) </u></em></strong>可以關(guān)注：<a href="https://www.yqgqt.org.cn/post/1939915" rel="noopener noreferrer" target="_blank">利用Python腳本 批量提取Abaqus的odb文件中Set節(jié)點集 初始節(jié)點坐標(biāo)、指定Step下的變形量、變形后節(jié)點坐標(biāo)_CAE Abaqus提取odb-技術(shù)鄰 (jishulink.com)</a></p><h2>1. 問題描述</h2><p>在工程仿真和分析領(lǐng)域，提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中的節(jié)點集的<span style="color: rgb(25, 27, 31);">反力RF（Reaction force）</span>是一項常見任務(wù)。然而，手動提取這些數(shù)據(jù)是一項繁瑣且容易出錯的工作。
展開 Abaqus DEM進階技術(shù)，指定顆粒級配
之前已有三期文章（
案例下載，
通過腳本生成顆粒，
粒子生成器-噴丸強化）不同程度的涉及到Abaqus DEM，今天我們來看一下粒子生成器-粒度級配的設(shè)定問題。
粒度級配（ASTM C136）
顆粒物的粒度級配性能將直接影響最終產(chǎn)品的質(zhì)量，比如混凝土的強度、抗開裂性能跟骨料尺寸分布有很大關(guān)系，所以，粒度級配在DEM分析中也是個十分重要的概念，它反映了不同尺寸顆粒所占的分量。在大量顆粒的填裝、混合、篩選等分析中，往往需要對不同尺寸的顆粒進行建模。
Abaqus的粒子生成器可以直接指定生成的粒子尺寸分布服從某一概率密度函數(shù)（Probability Density Function，以下簡稱PDF），因此可以解決DEM粒度級配的問題。
生成不同尺寸的粒子
Abaqus支持以下6種PDF：
Uniform/Normal/Log-normal
Piecewise linear/Discrete/Truncated
我們在下面的4個surface上，創(chuàng)建4個Particle generator（PG），每個PG分別引用不同的PDF，每一個面上指定生成2萬個顆粒。
4個粒子生成面
粒子尺寸分布信息
粒子生成器運行結(jié)果:
PG1正態(tài)分布-粒子半徑云圖
PG2對數(shù)正態(tài)分布-粒子半徑云圖
PG3均勻分布-粒子半徑云圖
PG4分段線性分布-粒子半徑云圖
從生成結(jié)果統(tǒng)計來看，顆粒尺寸分布整體上與指定的概率密度函數(shù)符合的很好，不過，對于微小顆粒而言（此問題中半徑小于0.001mm的顆粒），生成還是比較困難的，尤其是顆粒尺寸懸殊較大時。
展開 在指定文件夾中建立ABAQUS文件的快捷方法
在指定文件夾中建立ABAQUS文件有兩種方法：
1、一種是用通常的界面打開ABAQUS，建立模型后用另存為的方式放在指定文件夾里；
2、方法二介紹如下：
1）建立指定文件夾如ABAQUS exe
2)打開文件夾后，按住SHIFT鍵
3）右擊鼠標(biāo)，選擇“在此處打開命令窗口”
4）在文件夾名后輸入：abaqus cae
5) 即可打開ABAQUS CAE，且文件出現(xiàn)在指定文件夾里。
批量提取Abaqus的節(jié)點坐標(biāo)(初始坐標(biāo)、指定Step下的變形量、變形后節(jié)點坐標(biāo)) ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發(fā)來批量提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中指定Step下的Set節(jié)點集變形量。通過詳細(xì)的步驟說明、代碼示例和圖片展示，您將學(xué)會如何使用該腳本，自動化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，U1，U2）。</p><p>如果還需要按Increment提取每個增量下的變形后的節(jié)點坐標(biāo)的話，在提取變形量的基礎(chǔ)上，與初始坐標(biāo)進行簡單的計算就可以求得坐標(biāo)。 （備注：該代碼只提取了x，y方向的變形量）</p><h2>1. 問題描述</h2><p>在工程仿真和分析領(lǐng)域，提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中的節(jié)點集變形量是一項常見任務(wù)。然而，手動提取這些數(shù)據(jù)是一項繁瑣且容易出錯的工作。因此，需要一種自動化的方法來批量提取指定步驟下按節(jié)點集組織的變形量數(shù)據(jù)。</p><h2>2. 實例展示</h2><p>假設(shè)我們有一個名為`example.odb`的ODB文件，其中包含名為`Step-x`的步驟和名為`Set-x`的節(jié)點集。運行以上代碼后，腳本會自動將該步驟下節(jié)點集的變形量提取出來，并保存為`NodalDisplacement.csv`文件。
展開 
ABAQUS分析用戶手冊-指定條件、約束與相互作用卷
如題，這本書。
Abaqus中利用Python腳本獲取指定Set上的位移和應(yīng)力
Abaqus內(nèi)置了一個函數(shù)名為 getSubset 的函數(shù)可以獲取指定 集（Set）上的Field Data
比如要獲取如下圖片中名為SET-1的節(jié)點集的位移、名為SET-1的單元集的位移，可以通過圖3所示的代碼來實現(xiàn)。
abaqus子程序VUSDFLD——考慮應(yīng)變率與應(yīng)變軟化效應(yīng)的軟土模型 ￥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提出的西澳模型，用于求解軟黏土體劇烈變形后的強度變化，可應(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)力或平均應(yīng)變 ￥10
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)變適用于各種邊界條件，但需要對每個單元進行應(yīng)力（應(yīng)變）的輸出和計算。
ABAQUS用戶手冊：材料卷 、單元卷、分析卷、指定條件、約束與相互作用卷、介紹，空間建模，執(zhí)行與輸出 ￥1.2
ABAQUS用戶手冊及關(guān)鍵詞參考指南：初學(xué)者必備6件套
1材料卷
2單元卷
3分析卷
4指定條件、約束與相互作用卷
5介紹，空間建模，執(zhí)行與輸出
6工具包
7Abaqus關(guān)鍵詞參考指南
有限元理論基礎(chǔ)及Abaqus內(nèi)部實現(xiàn)方式研究系列27： Abaqus內(nèi)部計算和顯示的應(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)變后，進入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中采用自帶材料進行計算，材料參數(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ì)量矩陣研究。
展開 ABAQUS批量提交Job與Python讀取ODB結(jié)果應(yīng)力應(yīng)變數(shù)據(jù)
批量提交的核心是需要等待當(dāng)前計算任務(wù)的結(jié)束，上圖中若干個job的提交代碼如下：
# coding: utf-8
#微信公眾號：ABAQUS二次開發(fā)
#作者：阿信老師CAE
#email：axin_cae@163.com
#2022.03.17
from abaqus import *
from abaqusConstants import *
jobList = mdb.jobs.keys() #獲取所有計算任務(wù)的jobName
for jobName in jobList:
mdb.jobs[jobName].submit(consistencyChecking=OFF) #提交計算
mdb.jobs[jobName].waitForCompletion() #等待計算完成
print jobName , "is completed"
不過在實際的操作過程中，我們可能會需要避開一些job不提交，或者job窗口事先并沒有建立job，等等，總之實際問題永遠(yuǎn)比任何教程都復(fù)雜，不過只需要靈活面對就可以了，處理起來并不難。
展開 
Abaqus通過USDFLD子程序進行泥巖的應(yīng)變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應(yīng)變軟化特性的泥巖彈塑性本構(gòu)模型，使用Abaqus及其子程序?qū)δ鄮r的三軸壓縮試驗進行了數(shù)值模擬。泥巖在受壓過程中主要經(jīng)過了5個階段，即壓密階段、彈性變 形、應(yīng)變硬化、應(yīng)變軟化、殘余階段。泥巖應(yīng)變軟化模型如下所示。
式中，ξ為強度參數(shù)，ξp為峰值強度參數(shù)，ξr為殘余階段強度參數(shù)，η為應(yīng)變軟化參數(shù)，η*為殘余階段的應(yīng)變軟化參數(shù)初始值。對于三軸壓縮試驗，η用塑性剪切應(yīng)變來表示
塑性屈服準(zhǔn)則采用Mohr-Coulomb準(zhǔn)則，則粘聚力和內(nèi)摩擦角的參數(shù)演化可以用下式表示
在巖石的塑性變形過程中會產(chǎn)生比較明顯的剪脹現(xiàn)象，而用來描述這一現(xiàn)象的較常用的力學(xué)參數(shù)就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現(xiàn)，流程圖如下
有限元模型如下圖所示
計算得到的應(yīng)力云圖及不同圍壓下的載荷位移響應(yīng)如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學(xué)模型.中國科技論文,2016,11(21):2456-2461
有關(guān)于abaqus子程序開發(fā)的相關(guān)問題可以通過公眾號聯(lián)系我們.
公眾號: 320科技工作室
展開 Abaqus通過USDFLD子程序進行泥巖的應(yīng)變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應(yīng)變軟化特性的泥巖彈塑性本構(gòu)模型，使用Abaqus及其子程序?qū)δ鄮r的三軸壓縮試驗進行了數(shù)值模擬。泥巖在受壓過程中主要經(jīng)過了5個階段，即壓密階段、彈性變 形、應(yīng)變硬化、應(yīng)變軟化、殘余階段。泥巖應(yīng)變軟化模型如下所示。
式中，ξ為強度參數(shù)，ξp為峰值強度參數(shù)，ξr為殘余階段強度參數(shù)，η為應(yīng)變軟化參數(shù)，η*為殘余階段的應(yīng)變軟化參數(shù)初始值。對于三軸壓縮試驗，η用塑性剪切應(yīng)變來表示
塑性屈服準(zhǔn)則采用Mohr-Coulomb準(zhǔn)則，則粘聚力和內(nèi)摩擦角的參數(shù)演化可以用下式表示
在巖石的塑性變形過程中會產(chǎn)生比較明顯的剪脹現(xiàn)象，而用來描述這一現(xiàn)象的較常用的力學(xué)參數(shù)就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現(xiàn)，流程圖如下
有限元模型如下圖所示
計算得到的應(yīng)力云圖及不同圍壓下的載荷位移響應(yīng)如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學(xué)模型.中國科技論文,2016,11(21):2456-2461
有關(guān)于abaqus子程序開發(fā)的相關(guān)問題可以聯(lián)系扣扣1653004885或者關(guān)注cae320公眾號
展開 ABAQUS 小應(yīng)變分析(例3) 條形基礎(chǔ)或海洋淺基礎(chǔ)下壓模擬(Tresca 本構(gòu)) ￥67
ABAQUS 小應(yīng)變分析(例3) 條形基礎(chǔ)或海洋淺基礎(chǔ)下壓模擬(Tresca 本構(gòu))
條形基礎(chǔ)承載力是工程廣泛關(guān)注的問題，例如陸地條形基礎(chǔ)和海洋淺基礎(chǔ)。該模擬地基為飽和不排水的粘土，采用Tresca本構(gòu)，粘土強度su = 15 kPa。條形基礎(chǔ)處理成剛體。最終數(shù)模結(jié)果顯示，條形基礎(chǔ)的無量綱承載力Nc0 = F/Asu 近似于 pi + 2 = 5.14, 與傳統(tǒng)理論解極好的契合。
建模過程及結(jié)果：
荷載及位移邊界條件
網(wǎng)格劃分
局部網(wǎng)格劃分
條形基礎(chǔ)的力位移曲線(已達到極限承載力)
地基的土體應(yīng)力分布
地基的土體破壞模式
Abaqus平均應(yīng)力和應(yīng)變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應(yīng)力及應(yīng)變進行自動提取并計算平均值
能夠得到每一幀的應(yīng)力和應(yīng)變平均值，并保存到CSV文件中
所得到的應(yīng)力包括S11,S22,S33,S12,S13,S23以及Mises七個應(yīng)力平均值，以及E11,E22,E33,E12,E13,E23六個應(yīng)變平均值
