abaqus體積應(yīng)變的案例
代表性體積單元根據(jù)單元體積應(yīng)力應(yīng)變加權(quán)平均 ￥20
現(xiàn)如今，越來越多的人開始對復(fù)合材料性能進行研究，如何通過<a href="/major/abaqus">abaqus提取代表性體積單元是非常重要的，我提供了一種可以根據(jù)單元體積進行應(yīng)力應(yīng)變平均的代碼，希望對大家有用。
ABAQUS提取面積或體積
體積也是同樣的道理，在這就不演示啦！
附圖是站內(nèi)評論區(qū)看到的，再次感謝大佬！
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結(jié)果提取大于某值的區(qū)域體積（轉(zhuǎn)載）
經(jīng)常會有ABQer初學者會問我：“怎么獲得屈服區(qū)域的體積？”、“怎么獲得大于一定應(yīng)力的單元體積？”之類的一些問題。我們將帶大家探索ABAQUS后處理的兩種提取方法：CAE界面提取 和 Python后處理二次開發(fā)提取。本節(jié)先介紹第一種方法。
還是以最簡單的懸臂梁為例，提取加載歷程下大于100e6 Mises等效應(yīng)力下單元體積：
1）顯示應(yīng)力變形云圖，并通過云圖顯示設(shè)置，可將超過指定數(shù)值的云圖范圍顯式為灰色，便于觀察和對比：
2）顯示指定應(yīng)力范圍內(nèi)的單元：找到按鈕，或菜單欄Tools->DisplayGroup->Create...，點擊后激活Create Display Group對話框，選擇Item為Elements、method為Result value，并設(shè)置最小值（Min value）100e6和最大值（Max value）1e10，點擊Apply可顯示出范圍內(nèi)的單元，如下圖所示。
注意：同時點擊底部的Save As...按鈕，保存名字為DisplayGroup-2。
3）測量顯示單元體積：找到按鈕，或者菜單欄Tools->Query...命令，找到Mass properties選型，并點擊后激活查詢質(zhì)量流程，在提示信息欄將讓選擇單元，修改選擇方法為Display groups，并選擇剛剛生成的DisplayGroup-2；由于是測量質(zhì)量，針對材料中沒有定義密度的情況，ABAQUS友好的提供了Options選項，可以設(shè)置材料密度和平面shell的厚度，這里均設(shè)置為1。如下圖所示。
重新選擇顯式組后，在底部信息框中將顯示測量的體積、體積中心、質(zhì)量、質(zhì)量中心以及質(zhì)量距等數(shù)據(jù)。
展開 
ABAQUS結(jié)果提取大于某值的區(qū)域體積-CAE方法
以最簡單的懸臂梁為例，提取加載歷程下大于100e6 Mises等效應(yīng)力下單元體積：
1）顯示應(yīng)力變形云圖，并通過云圖顯示設(shè)置，可將超過指定數(shù)值的云圖范圍顯式為灰色，便于觀察和對比：
2）顯示指定應(yīng)力范圍內(nèi)的單元：找到按鈕，或菜單欄Tools->DisplayGroup->Create...，點擊后激活Create Display Group對話框，選擇Item為Elements、method為Result value，并設(shè)置最小值（Min value）100e6和最大值（Max value）1e10，點擊Apply可顯示出范圍內(nèi)的單元，如下圖所示。
注意：同時點擊底部的Save As...按鈕，保存名字為DisplayGroup-2。
3）測量顯示單元體積：找到按鈕，或者菜單欄Tools->Query...命令，找到Mass properties選型，并點擊后激活查詢質(zhì)量流程，在提示信息欄將讓選擇單元，修改選擇方法為Display groups，并選擇剛剛生成的DisplayGroup-2；由于是測量質(zhì)量，針對材料中沒有定義密度的情況，ABAQUS友好的提供了Options選項，可以設(shè)置材料密度和平面shell的厚度，這里均設(shè)置為1。如下圖所示。
重新選擇顯式組后，在底部信息框中將顯示測量的體積、體積中心、質(zhì)量、質(zhì)量中心以及質(zhì)量距等數(shù)據(jù)。
4）匯總不同時刻的數(shù)值，繪制時程曲線如下：
同時，該方法不僅適用于應(yīng)力數(shù)值范圍，基本所有的云圖范圍單元體積的測量都是可以滿足的。
作者：陳佳敏cn
來源：CAE技術(shù)資訊
展開 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)變）的輸出和計算。
有限元理論基礎(chǔ)及Abaqus內(nèi)部實現(xiàn)方式研究系列27： Abaqus內(nèi)部計算和顯示的應(yīng)變
（1）顯示應(yīng)變：Abaqus計算完畢后得到導入結(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)變。一般理論書都認為Abaqus是因為對數(shù)應(yīng)變計算復(fù)雜才采用別的應(yīng)變，但個人認為應(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ù)
批量提交的核心是需要等待當前計算任務(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，等等，總之實際問題永遠比任何教程都復(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)變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內(nèi)摩擦角的參數(shù)演化可以用下式表示
在巖石的塑性變形過程中會產(chǎn)生比較明顯的剪脹現(xiàn)象，而用來描述這一現(xiàn)象的較常用的力學參數(shù)就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現(xiàn)，流程圖如下
有限元模型如下圖所示
計算得到的應(yīng)力云圖及不同圍壓下的載荷位移響應(yīng)如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關(guān)于abaqus子程序開發(fā)的相關(guān)問題可以通過公眾號聯(lián)系我們.
公眾號: 320科技工作室
展開 星辰技文|ABAQUS結(jié)果提取大于某值的區(qū)域體積-Python方法 ￥10
上一篇技文《ABAQUS結(jié)果提取大于某值的區(qū)域體積-CAE方法》中帶大家使用ABAQUS CAE界面直接提取大于100e6 Mises等效應(yīng)力的區(qū)域體積，方法雖然比較好操作，但也存在明顯的缺點：
方法太繁瑣；統(tǒng)計歷程曲線時會顯得非常麻煩，因此我們找到了代步工具：Python。
精度較低；這是由于一個單元內(nèi)只需要一個積分點滿足數(shù)值要求，整個單元都會顯示出來，但一個單元可能存在多個積分點，比如演示案例中采用六面體二階減縮積分單元（C3D20R）存在8個積分點，單元內(nèi)部分積分點可能不滿足條件，也會被統(tǒng)計在內(nèi)，從而導致統(tǒng)計的體積偏大，我們對比下CAE方法和Python方法的結(jié)果差異如下圖，最少都差32%！
樣的誤差雖然可以通過提高網(wǎng)格密度來減緩，但并不能完全避免，因此也就有必要對這些單元內(nèi)所有積分點值進行判斷，獲得單元內(nèi)滿足條件的積分點所占比例，再對其體積進行加權(quán)相加。而這樣的操作過程顯然不適合手動統(tǒng)計，再一次將目光投向了我們的代步工具：Python！
沒有Python基礎(chǔ)的小伙伴，建議先看曹金鳳姐姐的《Python語言在Abaqus中的應(yīng)用》或江丙云哥哥的《ABAQUS Python二次開發(fā)攻略》，不然會有些吃力。雖然代碼備注的已經(jīng)非常詳細，但還是需要一些背景支撐的。
【編程需求與分解】
為了與CAE過程進行對比，我們還是將需求定義為：提取懸臂梁加載過程中Mises等效應(yīng)力大于100e6的區(qū)域體積。
展開 ABAQUS水力壓裂模擬|XFEM和Cohesive方法-多縫、交叉縫、體積縫、轉(zhuǎn)向縫、縫間干擾
ABAQUS水力壓裂模擬|XFEM和Cohesive方法
by 星辰北極星
關(guān)鍵字：單縫、多縫、交叉縫、體積縫、轉(zhuǎn)向縫、縫間干擾、儲隔層
水力壓裂，對于石油工程的朋友并不陌生，它是石油開采和增產(chǎn)的重要手段；也廣泛應(yīng)用于地熱開采、地基處理等領(lǐng)域。由于畢業(yè)于石油大學，所以有很多機會接觸這方面的問題，也關(guān)注著ABAQUS在壓裂領(lǐng)域的應(yīng)用。這個專題將分享自己在水力壓裂仿真中的一些積累，希望大家喜歡。
【主要內(nèi)容】
一、課程概述
二、仿真要點介紹
2.1 ABAQUS水力壓裂模擬常用仿真方法
2.2 地應(yīng)力平衡分析（Geostatic）
2.3 滲流-位移耦合分析（Soils）
2.4 材料與單位制講解
2.5 特殊的輸出需求與定義
2.6 交叉裂縫處理
三、實例講解
3.1 基于Cohesive單元的二維水力壓裂模擬
3.2 基于Cohesive單元的三維水力壓裂模擬
3.3 水力裂縫與天然裂縫相交模擬-Cohesive單元法
3.4 裂縫發(fā)育地層的水力壓裂模擬-Cohesive單元法
3.5 基于XFEM的水力裂縫轉(zhuǎn)向模擬
3.6 基于XFEM的水平井多段壓裂裂縫的縫間干擾問題研究
視頻地址：https://i.xue.taobao.com/detail.htm?spm=a2174.7765247.0.0.OHNzvF&courseId=89321
【二維水力壓裂模擬（Cohesive）】
通過這個簡單的案例講述采用Cohesive單元模擬水力壓裂的基本技巧，讓大家掌握注液、停泵憋壓等基本設(shè)置，以及前后處理的一些技巧。
展開 
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)變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內(nèi)摩擦角的參數(shù)演化可以用下式表示
在巖石的塑性變形過程中會產(chǎn)生比較明顯的剪脹現(xiàn)象，而用來描述這一現(xiàn)象的較常用的力學參數(shù)就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現(xiàn)，流程圖如下
有限元模型如下圖所示
計算得到的應(yīng)力云圖及不同圍壓下的載荷位移響應(yīng)如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關(guān)于abaqus子程序開發(fā)的相關(guān)問題可以聯(lián)系扣扣1653004885或者關(guān)注cae320公眾號
展開 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)變平均值
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)變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據(jù)需要進行修改。
