abaqus各種應變的案例
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中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
解決Abaqus Python二次開發各種問題！ 附ABAQUS PYTHON二次開發攻略文檔下載
網絡接口的環境和abaqus的python2環境是天然隔離的，安裝各種第三方庫自然不在話下了。
說到成果保護，作者可以精準對具體某個函數進行訪問控制，包括“哪個賬號可以訪問？”、“可以訪問幾次？”、“在哪個時間段可以訪問？”等等。
而且，這里提到的網絡可以是局域網。最小的局域網就是同一臺計算機訪問自身的不同端口，所以完全可以在你的計算機上配置好網絡服務器應用，然后在abaqus的python2環境中訪問、交互，實現上述的各種功能。
2、函數上云
對于那些飽含作者心血的寶貴函數，完全可以部署到各種云服務器上，在代碼“絕對安全”和用戶“絕對方便”的前提上讓有需要的人能夠方便快捷的使用到。
以上，都可以基于Abaqus python 漸進式二次開發的技術來實現。代碼的使用方法很大程度決定了代碼使用的場景和邊界。
所以，希望突破限制復制粘貼“長長”腳本的你，歡迎參加《Abaqus python 漸進式二次開發（二）》直播課，從基礎的開發“套路”開始，重新認識Abaqus python二次開發。
下載地址：ABAQUS PYTHON二次開發攻略文檔下載
展開 ABAQUS中的各種文件類型詳解（轉）
的分析結果要寫入此文件中則需要轉換，convert=select 或convert=all
10. abaqus.rpy
記錄一次操作中幾乎所有的ABAQUS/CAE命令
11. job_name.lck
阻止并發寫入輸出數據庫，關閉輸出數據庫則自行刪除
12. model_database_name.rec
包含用于恢復內存中模型數據庫的ABAQUS/CAE命令
13. job_name.ods
場輸出變量的臨時操作運算結果，自動刪除
14. job_name.ipm
內部過程信息文件：啟動ABAQUS/CAE分析時開始寫入，記錄了從 ABAQUS/Standard或ABAQUS/Explicit 到 ABAQUS/CAE的過程日志
15．job_name.log
日志文件：包含了ABAQUS執行過程的起止時間等
16．job_name.abq
ABAQUS/Explicit模塊才有的狀態文件，記錄分析、繼續和恢復命令。為restart所需的文件。
17．job_name.mdl
模型文件：在ABAQUS/Standard和ABAQUS/Explicit中運行數據檢查后產生的文件，.在 analysis和continue 指令下被讀入并重寫，為restart所需的文件。
18．job_name.pac
打包文件：包含了模型信息，僅用于ABAQUS/Explicit ，該文件在執行 analysis、datacheck命令時寫入，執行 analysis, continue, recover 指令時讀入，restart時需要的文件。
19．job_name.prt
零件信息文件：包含了零件與裝配信息.。
展開 
有限元理論基礎及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-MPC各種約束搞清楚了嗎？
很多時候，在處理復雜裝配體的時候，我們都用到MPC和Coupling。那么這么多種約束各適用于什么情況，對于承載又有什么區別呢？今天就用一個shell單元的簡單案列和大家一起探討一下。
圖1 三維模型
如圖一的T形連接件，通常實際工程結構是焊接件，但在有限元分析中，我們如何處理這種連接呢？通常我們會用共節點來處理，共節點把兩部分當成一個部件，可以考慮連接處的變形。但也有時候，要實現共節點，需要較長時間進行前處理，因此考慮用MPC來進行連接。今天對于這個薄壁件（壁厚1mm），我建立6個模型分別是：1 共節點2kinematic coupling 3MPC BEAM 4MPC Link 5 MPC Tie 6 MPC Pin ，分別用7種方式處理這個連接，每個模型建立三個分析步，第一個分析步施加拉伸載荷，第二個分析步施加豎向載荷產生彎矩，第三個分析步施加扭轉載荷，如圖二所示。
圖2 三種載荷工況
結果分析
求解過程中，使用MPC-Pin和MPC-Link的兩個模型算到第二步（彎曲載荷）時第一個增量步不收斂，其余模型順利計算完成。究其原因，我們知道桿單元是不能承受彎曲載荷的，這里的Link便是將兩節點使用剛性桿單元連起來，而Pin是將兩節點的三個平移自由度綁定到一起，也不能承受彎曲載荷。同樣扭矩載荷也是不收斂的。
再看看拉伸載荷工況結果：6種處理方式在拉伸載荷下結果幾乎是一致的，可見，它們均可傳遞拉力載荷。
展開 ABAQUS批量提交Job與Python讀取ODB結果應力應變數據
問題2
如何從ODB結果文件中讀取計算結果數據
第二個問題是大家經常要面對的，常規的數據提取采用ABAQUS自帶的后處理工具已經能滿足需求了，但是如果我們需要進行一些寫出文本，數據篩選，數據的二次計算處理。。。，我們就需要用程序語言的幫助了。
針對今天的問題——
提取結果并寫出到文本
，如果用自帶的工具，可以導出abaqus.rpt文件，不過這種方式并不能滿足我們較為苛刻的要求，用python程序會的自由度會更大一些。
另外需要說明的是，如果我們僅僅只是看看某單元或者節點的各種計算結果，或者簡要的數據處理，臨時用一用，可以采用如下方式，不必費心寫復雜的程序，不值當浪費時間，
青春易逝，頭發寶貴
。當然了，通過阿信這里給出的基本程序改一改，能省幾根頭發。
限于時間關系，這里只列出基本的應力提取程序，寫入文本比較簡單，隨意百度即可實現，不再給出。
展開 Abaqus通過USDFLD子程序進行泥巖的應變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應變軟化特性的泥巖彈塑性本構模型，使用Abaqus及其子程序對泥巖的三軸壓縮試驗進行了數值模擬。泥巖在受壓過程中主要經過了5個階段，即壓密階段、彈性變 形、應變硬化、應變軟化、殘余階段。泥巖應變軟化模型如下所示。
式中，ξ為強度參數，ξp為峰值強度參數，ξr為殘余階段強度參數，η為應變軟化參數，η*為殘余階段的應變軟化參數初始值。對于三軸壓縮試驗，η用塑性剪切應變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內摩擦角的參數演化可以用下式表示
在巖石的塑性變形過程中會產生比較明顯的剪脹現象，而用來描述這一現象的較常用的力學參數就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現，流程圖如下
有限元模型如下圖所示
計算得到的應力云圖及不同圍壓下的載荷位移響應如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關于abaqus子程序開發的相關問題可以通過公眾號聯系我們.
公眾號: 320科技工作室
展開 ABAQUS材料庫（常用各種金屬鋼筋混凝土等） ￥2
1606716696(1).png
流程.pdf
在安裝ABAQUS的盤里尋找abaqus_plugins文件夾
然后將附件里的文件復制
粘貼到制定位置即可
Abaqus/Explicit分析重要概念（2）：各種阻尼的功能及設置方法/橡膠阻尼
在 Abaqus/Explicit 分析中，為了避免數值振蕩，一般都需要定義模型的阻尼，
定義方法主要包括以下幾種：
1）體積粘性（bulk viscosity）
體積粘性用于引入由于體積應變引起的阻尼，在研究高速動力分析的高階性能時，體積粘性是尤其必要的。體積粘性只是作為一個數值效應被引入，因此，材料點上的應力并不考慮體積粘性壓力的影響。
Abaqus/Explicit 有兩種體積粘性參數：線性體積粘性和二次體積粘性，可以在 Step 功能模塊中進行設置（如圖1所示）。
一般情況下，采用 Abaqus 的默認設置即可。
圖1 設置體積粘性參數
2）材料阻尼
常用的材料阻尼是瑞利（Rayleigh）阻尼，在Property模塊的Mechanical菜單下定義（如圖2所示），它包含兩個阻尼參數：
質量比例阻尼是關于質量矩陣的比例系數，主要用于消除低階振蕩；剛度比例阻尼是關于剛度矩陣的比例系數，主要用于消除高階振蕩。
圖2 設置材料阻尼
關于材料阻尼的詳細介紹，請參見 Abaqus 幫助文檔《Abaqus Analysis User’s Manual》第20.1.1節“Material damping”和《Abaqus Keywords User’s Manual》中的關鍵詞
* DAMPING。
3）阻尼器（dashpot）單元
在 Property 功能模塊和 Interaction 功能模塊的Special菜單中都可以定義阻尼器單元（如圖3所示），其優點是可以僅在必要的節點上定義阻尼，其阻尼力與單元的兩個節點相對速度成正比。阻尼器單元必須與其他單元（如彈簧單元或桁架單元）同時使用，一般不會引起穩定極限值的顯著變化
。
展開 ABAQUS提取單元平均應力/應變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據需要進行修改。
Abaqus平均應力和應變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應力及應變進行自動提取并計算平均值
能夠得到每一幀的應力和應變平均值，并保存到CSV文件中
所得到的應力包括S11,S22,S33,S12,S13,S23以及Mises七個應力平均值，以及E11,E22,E33,E12,E13,E23六個應變平均值
Abaqus通過USDFLD子程序進行泥巖的應變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應變軟化特性的泥巖彈塑性本構模型，使用Abaqus及其子程序對泥巖的三軸壓縮試驗進行了數值模擬。泥巖在受壓過程中主要經過了5個階段，即壓密階段、彈性變 形、應變硬化、應變軟化、殘余階段。泥巖應變軟化模型如下所示。
式中，ξ為強度參數，ξp為峰值強度參數，ξr為殘余階段強度參數，η為應變軟化參數，η*為殘余階段的應變軟化參數初始值。對于三軸壓縮試驗，η用塑性剪切應變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內摩擦角的參數演化可以用下式表示
在巖石的塑性變形過程中會產生比較明顯的剪脹現象，而用來描述這一現象的較常用的力學參數就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現，流程圖如下
有限元模型如下圖所示
計算得到的應力云圖及不同圍壓下的載荷位移響應如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關于abaqus子程序開發的相關問題可以聯系扣扣1653004885或者關注cae320公眾號
展開 ABAQUS中對應力、應變的部分理解
對應力的部分理解
對應變的部分理解
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
ABAQUS 小應變分析(例3) 條形基礎或海洋淺基礎下壓模擬(Tresca 本構) ￥67
ABAQUS 小應變分析(例3) 條形基礎或海洋淺基礎下壓模擬(Tresca 本構)
條形基礎承載力是工程廣泛關注的問題，例如陸地條形基礎和海洋淺基礎。該模擬地基為飽和不排水的粘土，采用Tresca本構，粘土強度su = 15 kPa。條形基礎處理成剛體。最終數模結果顯示，條形基礎的無量綱承載力Nc0 = F/Asu 近似于 pi + 2 = 5.14, 與傳統理論解極好的契合。
建模過程及結果：
荷載及位移邊界條件
網格劃分
局部網格劃分
條形基礎的力位移曲線(已達到極限承載力)
地基的土體應力分布
地基的土體破壞模式
