ABAQUS計算的案例
PyQt平臺開發Abaqus計算程序
最終實現的功能如圖4所示，生成的py文件：
圖4 Py腳本文件
（2）開始計算
啟動Abaqus調用Python腳本文件在PyQt里面比較好的實現方法是利用Abaqus Command命令提示符窗口。因此，在PyQt里面只需要通過subprocess函數啟動計算機里面的Abaqus Command并且在窗口輸入調用語句即可。
最終即可實現參數化計算，程序界面輸入參數，程序自行啟動Abaqus執行計算。
當然，程序還可以結合Python對Abaqus的計算結果進行后處理，此次不作闡釋。
展開 有限元理論基礎及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殼單元質量矩陣研究。
展開 案例13 基于模態的振動響應（Abaqus計算模態）
之前在superxjw版主的第二課中介紹了如何利用VL計算基于模態的振動響應，但是有網友是采用Abaqus計算模態，然后用VL來計算后續的振動響應以及聲學響應，然后就詢問如何導入Abaqus的模態分析結果，因此，做了一個導入Abaqus的模態結果，然后進行振動響應計算的案例，給大家分享一下。
superxjw版主的視頻教程：
LMS Virtual.Lab 11聲學視頻教程 第二課 基于模態的振動響應計算
對于VL的接口方面：
VL11SL2和VL12都是支持到Abaqus 6.12
所以，喜歡追求新版本，使用Abaqus6.13的朋友們就得注意一下版本的問題了。
感謝阿偉在本人學習LMS Virtual.Lab過程中的幫助！
本例視頻及Abaqus模態計算結果文件下載地址：http://pan.baidu.com/share/link?shareid=4100661600&uk=1728334102
LMS Virtual.Lab Acoustics 交流群 238339600
展開 abaqus單元平均應力的計算
abaqus計算結果可得到節點應力，那位大牛知道如何得到單元（或二維計算模型）的平均應力？
ABAQUS計算錯報 A LARGE NUMBER OF INCREMENTS
ABAQUS計算錯報 Error in job Job-concrete: THE ANALYSIS MAY NEED A LARGE NUMBER OF INCREMENTS (MORE THAN 20,000,000), AND IT MAY BE AFFECTED BY ROUND-OFF ERRORS. FOR ACCURACY, RUNNING DOUBLE PRECISION EXECUTABLE IS REQUIRED.
調試方案：注意單位問題，統一模型尺寸的單位和材料屬性里面的單位即可。是m還是mm單位制。
分享一個通過ABAQUS計算二維穿透裂紋應力強度因子的實例 ￥10
通過ABAQUS計算二維穿透裂紋應力強度因子。理論值為396，有限元計算值為407，誤差小于3%。
Python提取Abaqus計算結果的History Output并用Matplotlib繪圖
采用Python可以提取Abaqus計算結果的歷史輸出數據，提取出來后保存到本地并采用Matplitlib繪制折線圖，記錄一下代碼：以下三張圖片分別是讀取etotal歷史數據并保存到本地的代碼、繪制折線圖的代碼和etotal歷史變化折線圖。
基于Abaqus的VB計算程序開發
基于Abaqus的VB計算程序開發
1 概述
Python腳本文件可以用來編寫Abaqus執行有限元分析過程，而VB則可以生成Python腳本文件，結合兩者即可實現Abaqus計算程序的開發。
本次以一個平面開孔模型為例，如圖1所示，設計一個程序能夠通過程序界面設置相關參數并執行計算。
圖1 平面開孔模型
2 程序界面
參數化的參數包括矩形長度w、矩形寬度h、圓的半徑r、材料的彈性模量E、泊松比以及壓力載荷p，并設置三個按鈕：
代碼轉換：將Python文件編譯成可以直接嵌入在VB里面；
生成腳本：利用VB生成Abaqus可以執行的腳本文件；
開始計算：啟動Abaqus，調用生成的腳本文件進行有限元分析。
同時還添加了一個PictureBox控件，用于顯示圖片，方便參數輸入。
最終得到的程序界面如圖2所示，給所有的Textbox命名，方便后臺代碼調用，并設置默認值。
圖2 程序界面
3 功能代碼
（1）代碼轉換
雙擊代碼轉換按鈕，進入代碼編輯頁面，該段代碼需要打開已經存在的input_1.txt文件，并且逐行讀取里面的內容，然后在每行字符串添加一些字符，輸出到output_1.txt文件中。
展開 abaqus簡單立方體胞元周期性邊界條件施加計算腳本源代碼 ￥39.9
</p><p class="ql-align-justify">在abaqus中通過file→run script選擇該腳本運行，腳本會輸出X，Y，Z三個方向的計算結果ODB文件。</p><p><br></p>
【iSolver案例分享59】 水下爆炸實驗常用結構-簡化船體梁的模態計算與對比（Abaqus、文獻）
本文參考了Zhou等人發表的論文[1]，利用Abaqus、iSolver軟件對其中的簡化船體梁結構進行了模態計算，主要對水下爆炸中備受關注的一階垂向模態結果（干、濕）進行了對比，以評估自主有限元軟件iSolver在計算精度、可靠性和便利性等方面的表現。
1 模型介紹
根據論文提供的信息，建立如下所示的簡化船體梁結構模型：長2.8米，寬0.3米，高0.08米，板厚0.003米。結構材料采用Q235。
2 干模態的計算與對比
干模態的計算中，在Abaqus和iSolver使用相同的設置。Q235的密度取7850 kg/m^3,楊氏模量取2.1e11 Pa，泊松比取0.3。結構有3700個S4R單元。具體如下圖所示。
結果對比如下所示：
3 濕模態的計算與對比
濕模態的計算中，在Abaqus使用聲學單元建立水域，在iSolver直接使用軟件內置的施加虛擬流體質量設置（用戶手冊第4.14節）。結果對比如下所示：
4 結論
綜合上述對比，iSolver軟件計算結果分別在干、濕模態方面均與文獻結果、Abaqus計算結果展現出高度的吻合性，具有精度高、可靠性好的優點。且內置了施加虛擬流體質量的功能，對于船舶濕模態的計算更具有便利性，在不需要對水域進行建模的情況下，取得了比Abaqus更貼近實驗的結果，十分適合用于船舶行業的模態分析。
展開 Abaqus,Ansys等模型計算完成后自動提醒 無需編程
前言
你在工作站上批量算模型的時候是否會頻頻去檢查計算進度？
你是否有過信心滿滿提交計算作業，結果過段時間回來看第一步就不收斂？
你在趕ddl時是不是有著“人可以休息，電腦不可以休息”的心態？
如果您曾遇到過以上的煩惱，
TaskReminder_v1.0或許可以幫助你更加高效地進行計算任務。
軟件用途
監測程序運行情況，在程序完成或中斷時本軟件會通過郵件發送提醒。
使用場景
1、在工作站或其它電腦上運行程序
2、程序運行時間長
現有功能
1、 監測Abaqus單任務運行情況：在當前運行的Job計算完成或中斷時發送提醒消息；
2、 監測Abaqus多任務運行情況：在批量模型均計算完成或中斷時發送提醒消息；
3、 監測文件夾大小：當使用其他軟件如Ansys、Python、Matlab工作時，本軟件可以通過工作目錄文件夾大小來對程序運行進度監測，當工作目錄文件大小不變時發送提醒消息。
使用方法
將本軟件安裝在被檢測程序運行的電腦，選擇需要被監測的工作目錄和任務模式，輸入接受提醒的郵箱號（微信可綁定QQ郵箱接受消息），憑使用碼進行使用。
開啟TaskReminder軟件，點擊選擇路徑按鈕，選擇Abaqus工作目錄，或需要監測狀態的文件夾路徑。
點選任務模式，選擇需要進行的任務類型。
輸入接收消息提醒的郵箱地址，推薦采用qq郵箱（在微信中可以幾乎無延遲地收到消息提醒），輸入軟件使用碼。
點擊測試按鈕，測試軟件與服務器的連接狀態。
點擊開始任務，即可實時監測軟件運行狀態。任務完成后收到提醒。
軟件使用界面
測試后，連接正常，可以啟動任務
軟件正在運行中，可以監測到Abaqus計算進度。
展開 
ABAQUS 二次開發（UEL+UMAT C3D8 單元動力隱式計算）
概述：采用UEL接口二次開發實現八節點單元，考慮BBAR修正，避免體積自鎖，對標ABAQUS自帶的C3D8單元，計算的剛度矩陣、質量矩陣和阻尼矩陣均與ABAQUS保持一致。并且采用UMAT子程序進行應力和應變數據的可視化，計算的應力應變數據同樣與ABAQUS保持一致，可視化效果同ABAQUS。以方塊的受動力簡諧荷載為例，采用上述程序，應用動力隱式計算分析步，最終計算的位移、應變等時程曲線均與ABAQUS保持一致。
（）模型信息
模型尺寸為10x10x10，彈性模量1e10，密度2000，泊松比0.25，荷載和邊界條件示意圖為：
一面的所有節點均固定。另一面所有的節點施加簡諧荷載，簡諧荷載曲線為（詳細的參數見附件）：
計算的網格圖為：
網格尺寸為1，共計10x10x=1000個單元。總的計算時長為1s，這只固定增量步長為0.01s，所以總增量步數為100。
（）計算結果
以上面網格圖中中間角點為例，提取加載向位移時程曲線如下圖：
0.89s時刻x向應力云圖為：
注：左一為ABAQUS計算結果，中間為umat實現線彈性各向同性本構計算結果，右一為UEL+UMAT計算結果，后面的云圖也是這樣排布。
0.89s時刻y向應力云圖為：
寫文字好麻煩，不寫了！放個動圖算了，有興趣的移步附件。
展開 abaqus計算受重力作用下的土體表面重應力為什么不是零
abaqus計算受重力作用下的土體表面重應力為什么不是零呢，這個是下面的個人一些想法，可以供小白參考。
舉一個最簡單例子，假設土體大小10X10X10米，材料密度2000kg/m3.彈性模量100Mpa，泊松比0.3，摩擦角30度，粘聚力30Kpa，只受重力作用，重力加速度取10。單元尺寸大小分別取0.5、1、2、5m。
計算地表豎向應力分別為0.5X104pa、1X104pa、2X104pa、5X104pa，可以看出，單元尺寸越小，地表單元的應力就越小，結果偏于更準確。這是因為重力是作用在每個單元的重心位置，該模型標準矩形，單元也規整，第一層每個單元的標高是單元網格尺寸的一半，第一層重心位置的應力就是密度X重力加速度X該層單元格重心深度，再通過有限元原理轉化到每個單元的節點上，可想而知，要想地表網格節點尺寸為0，必須是單元網格大小足夠小，接近于0，這就是為什么abaqus模擬巖土工程不準確的地方，不可能做到足夠小，一般巖土工程的模型都是比較大的（幾十米幾百米幾千米），模型越大網格尺寸會劃分的很大，精確度也越低。更多案例可以關注抖音abaquser。
展開 如何提高ABAQUS計算速度
如何提高ABAQUS計算速度
當問題的自由度和復雜程度比較小的時候，計算速度也許不是什么太大的問題，但當自由度很多，比如幾百萬個自由度，又有復雜的接觸搜索計算時，計算速度就很重要了。我不精通計算機原理，但在Windows和Linux上安裝使用過32位和64位ABAQUS，也在64位工作站上用Linux并行計算過，對計算速度有一點自己的心得，寫出來和大家討論下，對計算機比較精通的，希望能解釋下原理和表達下自己的看法。
我認為影響計算速度的因素主要有3個：
1.計算模型的大小和復雜程度；
2.ABAQUS中關于內存和硬盤使用的設置；
3.計算機的配置。
在ABAQUS幫助文件中說，一個有限元模型在分析計算中，會生成兩種臨時文件，第一種臨時文件是必須放在內存中，第二種臨時文件可以放在硬盤里也可以放在內存中。由于內存讀寫速度比硬盤讀寫速度快，所以如果計算機配置不變的情況下，要達到最快計算速度，就要在計算時讓兩種臨時文件都放在內存中。在計算前進行datacheck，然后在dat文件中會告訴你能使模型計算所需最小的內存量，也就是在計算中所生成的第一種文件的大小，還會告訴你最小I/O交換情況所需的內存量，也就是兩種文件都放在內存中時所需要的內存量，我稱這叫“全速計算”所需內存量。在ABAQUS中，有個內存使用上限設置參數，也就是你允許你的計算機中有多少內存讓ABAQUS使用，如果你允許ABAQUS使用的內存上限大于“全速計算”所需內存，則ABAQUS會將兩種臨時文件放入內存使用。如果允許內存大于最小所需內存，小于“全速計算”計算內存，則ABAQUS會把第二種文件一部分放入內存，一部分放入硬盤。如果允許內存小于最小所需內存，則無法計算。
展開 基于Abaqus的VB計算程序開發
1 概述
Python腳本文件可以用來編寫Abaqus執行有限元分析過程，而VB則可以生成Python腳本文件，結合兩者即可實現Abaqus計算程序的開發。
本次以一個平面開孔模型為例，如圖1所示，設計一個程序能夠通過程序界面設置相關參數并執行計算。
圖1 平面開孔模型
2 程序界面
參數化的參數包括矩形長度w、矩形寬度h、圓的半徑r、材料的彈性模量E、泊松比以及壓力載荷p，并設置三個按鈕：
代碼轉換：將Python文件編譯成可以直接嵌入在VB里面；
生成腳本：利用VB生成Abaqus可以執行的腳本文件；
開始計算：啟動Abaqus，調用生成的腳本文件進行有限元分析。
同時還添加了一個PictureBox控件，用于顯示圖片，方便參數輸入。
最終得到的程序界面如圖2所示，給所有的Textbox命名，方便后臺代碼調用，并設置默認值。
圖2 程序界面
3 功能代碼
（1）代碼轉換
雙擊代碼轉換按鈕，進入代碼編輯頁面，該段代碼需要打開已經存在的input_1.txt文件，并且逐行讀取里面的內容，然后在每行字符串添加一些字符，輸出到output_1.txt文件中。
展開 