ABAQUS 計算振幅的案例
ABAQUS混凝土損傷塑性模型-C30EXCEL計算表格（含計算結果）
這是根據GB50010-2010中混凝土結構設計規范中的混凝土本構模型，結合文獻所述的損傷因子定義，編制的計算C30混凝土非彈性應變和損傷因子的EXCEL表格。也是邊學變做，希望能和大家多交流。
C30砼本構（損傷塑性模型）.rar
ABAQUS混凝土損傷塑性損傷因子計算依據.rar
ABAQUS 在計算到step3后計算中斷，變形過大
我在原本可以計算的模型的基礎上修改的，只是刪了幾個樓板與梁連接的栓釘
Mesh Free-眼鏡剛強度校核計算，附Abaqus計算結果對比
采用Mesh Free對某品牌眼鏡整體剛度、強度進行校核，如下圖所示，Mesh Free支持在不用做幾何清理的前提下進行計算分析，導入模型部件可以包含細節特征，比如螺釘上的倒角。
眼鏡定義了5種線彈性材料：鏡架主體采用鈦合金；眼鏡片采用樹脂；螺釘等連接件采用鋼；鼻托和鏡片扎線采用兩種不同的尼龍材料。
眼鏡腿一只固定，另一只向上掰，加力1N，模擬分析此種工況下眼鏡整體結構的剛度、強度。
Mesh Free所有接觸面定義為完全剛性連接，Abaqus作同樣處理，不考慮非線性因素，對比二者的線性計算結果。
Mesh Free給出的眼鏡最大變形為23.92mm，Abaqus的結果為23.46mm。
Mesh Free給出的眼鏡最大應力為303.4MPa，Abaqus的結果為308.3MPa。
談談Mesh Free使用感受：
雖然我常用ABQ，但是不得不說，對于包含細節幾何特征的復雜裝配結構建模分析，Mesh Free真的要比Abaqus高效的多。
據我了解Mesh Free的非線性也在大力的開發之中，目前已經支持經典塑性材料非線性、邊界條件非線性也可以設置滑動和一般的摩擦接觸。
對不熟悉常規有限元操作的結構設計人員來說，不用幾何清理、不用劃網格是極好的體驗。
關鍵是Mesh Free的結果也確實很準，目前的CAE無非是追求更準的基礎上算的更快，這兩點Mesh Free無疑是滿足的。
Mesh Free
Abaqus
展開 ABAQUS混凝土損傷塑性模型-C30EXCEL計算表格（含計算結果） ￥3.7
計算表格（如下）中標黃部分的參數可自行設定后，EXCEL程序會自動計算“抗拉強度、非彈性應變、受拉損傷因子”。
隨機振動分析-abaqus（附一個電池包計算案例） ￥20
四、如何將時域隨機振動曲線轉換得到功率譜密度曲線
五、 隨機振動分析理論
附.常見功率譜密度曲線給出形式
附.以dB/oct形式給出的功率譜密度曲線如何計算
附.國標中定義的PSD譜總均方根加速度值是如何計算的？
六. 隨機振動分析案例-abaqus
第一步：計算結構模態，輸出位移和應力。
第二步：隨機振動分析
2.1 定義輸出頻率上下限和模態阻尼
2.2 定義PSD載荷及加載
2.3 定義輸出
2.4 隨機振動計算頭文件設置
2.5 隨機振動分析結果
2.6 隨機振動σ應力結果評價
Abaqus接觸非線性在有限元計算分析中的應用 附莊茁ABAQUS非線性有限元分析與實例下載
來源：有限元在線
ABAQUS的非線性主要在有三種：幾何非線性，材料非線性以及接觸非線性。接觸非線性在ABAQUS的有限元計算分析中應用非常廣泛，特別是動態顯式的求解，只要模型中包含兩個以上相互接觸的部件，就要用到接觸非線性。
ABAQUS接觸非線性的設置主要在Interation模塊中完成，設置接觸的屬性時，可以設置摩擦系數，阻尼系數，損壞，失效準則等非線性參數，如圖1所示。
如圖2所示，在接觸定義界面，可以選擇通用接觸、面-面接觸、自接觸等各種非線性接觸方式。
在接觸編輯界面，可以選擇機械約束方式為運動學接觸算法，或是懲罰接觸方式，還可選擇滑移方式為有限滑移或小滑移，如圖3所示。
這是對模型定義非線性接觸后得到的分析結果，以供參考。
下載地址：莊茁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殼單元質量矩陣研究。
展開 Abaqus+PyQt+Python平面變形歐拉角計算
下面以簡單例子介紹平面變形、指向歐拉角的計算，包括絕對歐拉角、相對歐拉角。
1 簡化模型
下面的六面體為表面殼模型，下面由三段梁支持，三段梁分別沿X、Y、Z軸向。六個面的厚度不同，在上側3個面施加不同的壓力，如下左圖所示。位移云圖如下右圖所示。
2 計算要求
計算六面體上面3個面的變形歐拉角，包括3個面的絕對歐拉角，平面2、3相對與平面1的相對歐拉角。平面1、2、3如下圖所示。
3 數據處理
使用平面節點坐標、位移數據計算平面變形歐拉角。可以使用Python腳本輸出平面節點編號、節點坐標（X、Y、Z）、節點位移（U1、U2、U3），如下圖所示。下圖為平面1的10個工況的數據文件，打開的文本文件中7列數據為節點編號、坐標、位移。
三個平面10個工況的節點數據文件如下圖所示。每個文件中包含一個工況一個平面的節點編號、坐標、位移數據。
4 絕對歐拉角計算
使用PyQt+Python開發了一個簡單的小軟件，計算絕對歐拉角、相對歐拉角。
首先計算各平面的絕對歐拉角。
計算平面1的10個工況的絕對歐拉角。
平面1變形的絕對歐拉角計算結果如下圖所示。
伴隨絕對歐拉角計算結果，軟件同時寫出了平面變形前后的坐標系數據，如下圖。每行18個數據，每3個數據為一個坐標軸向量，變形前后2個坐標系，6個坐標軸，18個數據。
5 相對歐拉角計算
利用計算絕對歐拉角時得到的坐標系文件，計算平面變形相對歐拉角，如下圖所示，計算平面2相當于平面1、平面3相對與平面1的相對歐拉角。
計算結果如下圖所示。
6 小結
上述軟件用的算法申請了發明專利，軟件申請了軟著。
展開 Abaqus-原來顯式計算也可以這么快
本文詳細介紹了 Abaqus/explicit 中的準靜態分析及其提高求解速度的方法。
我如何知道我的模擬是否是準靜態的？
準靜態問題是通常可以使用 Abaqus/Standard 解決的問題之一，但由于觸或材料的復雜性，可能難以收斂，從而導致大量迭代。具有挑戰性的非線性準靜態問題通常涉及：
接觸條件非常復雜，Abaqus/Standard 可能會因接觸問題而無法收斂；
非常大的變形可能導致嚴重的網格變形；
例如，通常在金屬成形分析中，我們面臨這樣的困難：用 Abaqus/Standard 來模擬這樣的問題確實很困難。
示例：模擬拉深過程中的撕裂
Abaqus/Explicit顯式準靜態分析問題
Abaqus/Explicit 對于高度非線性靜態（準靜態分析）問題的建模更加有效。對于涉及接觸和金屬成型等非常大變形的三維問題尤其如此。
應用 Abaqus/Explicit 模擬準靜態事件需要特別考慮。在自然時間段內對過程進行建模在計算上是不切實際的。從字面上看，需要數百萬次的時間增量。因此，我們在模擬中人為地提高加載過程的速度以獲得經濟的解決方案。
使用 Abaqus/Explicit 獲得經濟的準靜態分析解決方案的兩種方法是：
1、人為提高加載速率
我們可以通過提高加載率來人為地減少該過程的時間尺度。增加的加載速率會減少模擬的時間。將加載速率提高 f 倍，分析速度提高 f 倍。
2、采用質量縮放
它增加了穩定時間增量的大小，因此完成作業所需的增量更少。人為地將材料密度（質量縮放）增加 f*f 倍，可以將分析速度提高 f 倍。
在本文中，我們的重點是提高加載速率的方法。
展開 Abaqus批量計算工具
功能：
每個inp任務放到單獨的文件夾中，不會導致計算結果文件都混在一起
支持 .for 文件的子程序
輸出計算成功或失敗信息
4個文件的求解結果示例
剛開發完成，沒怎么測試，可能有bug，歡迎大家反饋。
Abaqus批量計算.zip
abaqus并行計算中CPU超限的解決辦法
abaqus并行計算中CPU超限的解決辦法-Abaqus Error: The number of cpus (16) exceeds the number of cpus available(8）
# 1.說明，
CPU配置為i9-9900K，8核心16線程
# 2. 出現問題：
abaqus job=XX user=XX.for input=XX.inp cpus=16 int
出現報錯：Abaqus Error: The number of cpus (16) exceeds the number of cpus available (8).
# 3. 問題解析
這種情況主要是配置文件關于并行計算設置導致的，尤其是ABAQUS2023，
# 4.解決辦法：
使用everything找到路徑下：X:\SIMULIA\EstProducts\2023\win_b64\SMA\site\custom_v6.env的
custom_v6.env文件，其中X代表盤符，
修改前
在文件最后加一句話：
import socket
mp_host_list=[[socket.gethostname(),16],]
其中16代表自己電腦的線程數，一定要與自己電腦的線程一樣
修改后
展開 
基于ABAQUS的橡膠懸置膠合件剛度仿真計算
本文將以一個懸置膠合件仿真的實例講解一下如何利用ABAQUS來獲取其三個方向的靜態特性。所用膠合件的數模圖如圖1所示。其設計圖紙上標注的三向剛度如表1所示，膠料硬度是邵氏50±5度。
圖1 膠合件結構
表1設計要求
1、 網格劃分
采用HYPERMESH對圖一懸置進行網格劃分到的有限元模型如圖2所示。
2、材料設置
把劃分好的網格導入ABAQUS中，設置其材料參數，由于不同本構模型對橡膠懸置膠合件剛度計算結果有一定的影響。結合何小靜，上官文斌發表的《橡膠隔振器靜態力－ 位移關系計算方法》一文的研究結果表明，Mooney-Rivlin 模型的計算精度最高，其相對誤差均小于10%，所以本文采用M-R模型進行計算。50度膠料的M-R材料常數C10=0.2969，C01=0.0584。
3、剛度求解
3.1求解X方向剛度
按表 1要求，做如下設置：在Z方向先預載8mm，再在X向加載500N。取值0~5.6mm,對X向靜剛度進行求解。
求得的力和位移關系見表2所示，用表中數據進行畫圖差值可得到圖3所示的X向靜剛度為38N/mm,與設計值非常接近，其變形云圖見圖4所示
表2 X向力和位移關系表
圖3 X向剛度差值結果
圖4 X向云變形圖
3.2求解Y方向剛度
按表 1要求，做如下設置：在Z方向先預載8mm，再在X向加載1000N。取值2~4mm,對Y向靜剛度進行求解。
求得的力和位移關系見表3所示，用表中數據進行畫圖差值可得如圖5到Y向靜剛度為98N/mm,與設計值80N/mm有一定差異，見圖4。其變形云圖見圖6.
展開 Abaqus中陶瓷本構模型及其數值計算應用
2 數值計算軟件中本構模型
陶瓷由于其波速高、模量大，具有良好的抗侵徹性能，在各類型裝甲設計中被廣泛應用。而JH本構形式簡單，易于理解，已成為Abaqus、LS-DYNA和Autodyn等商用軟件的內嵌本構模型，可一定程度上滿足日常使用及工程計算要求。
對于陶瓷材料Abaqus幫助中給出了3種本構模型，Extended Drucker-Prager本構（以下簡稱DP本構）、JH-2和JHB本構模型。DP本構多用來模擬巖土材料（粒狀土壤和巖石），擴展DP本構給出的應力與壓力的關系也與JH本構中未損傷時應力與壓力的關系類似，其損傷段定義采用等效塑性應變與應力三軸度的對應關系進行定義，狀態方程采用Mie-Grüneisen形式（詳見Abaqus相應部分幫助）。
Abaqus官方幫助中給出的JHB本構模型參數如表1所示。其中標紅部分與Abaqus幫助（2021版本）不同，應為幫助原文疏漏。
表1 JHB本構模型參數
JHB本構模型的應力與壓力關系主要分為完整（Intact）和損傷（Failed）兩部分，表1中下標帶有 i 的即為完整部分相應參數，下標 f 即代表損傷部分參數；雖然JHB本構模型公式中考慮了脆性材料的相變特性，表1標藍部分參數應為對應的相變參數，但幫助中全部設置為0，推知官方幫助中給出的這組參數不能考慮陶瓷相變的影響。
Abaqus官方幫助中給出的JH-2本構模型參數如表2所示。
表2 JH-2本構模型參數
JH-2本構模型以無量綱形式描述了應力和壓力的關系，以Hugoniot極限下的壓力對壓力變量進行了無量綱化。
JHB本構模型狀態變量如表3所示。
展開 Abaqus中殘余模態的計算 ￥200
給出了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單位制。
