abaqus收斂規則的案例
Abaqus單元的命名規則
概述
Abaqus/Standard殼單元庫包括：
? 用于三維殼體幾何的單元；
? 具有軸對稱變形的軸對稱幾何的單元；
? 具有關于一個平面對稱的一般變形的軸對稱幾何的單元；
? 用于應力/位移、傳熱和完全耦合的溫度位移分析的單元；
? 通用單元，以及專門適用于分析“厚”或“薄”殼的單元；
? 使用簡化或減縮積分的通用、三維、一階單元；
? 考慮有限膜應變的單元；
? 在可能的情況下，每個節點使用五個自由度的單元，以及每個節點始終使用六個自由度的單元；
? 連續體殼單元。
Abaqus/Explicit 殼單元庫包括：
? 通用三維單元，用于模擬考慮有限膜應變的“厚”或“薄”殼；
? 小應變單元；
? 完全耦合的溫度位移分析單元；
? 用于具有軸對稱變形的軸對稱幾何形狀的單元；
? 連續體殼單元。
命名約定
殼單元的命名約定取決于單元的維度。
三維殼單元
Abaqus中的三維殼單元命名如下：
第一項：
常規應力/位移單元（S）
連續應力/位移單元（SC）
三角形應力/位移薄殼（STRI）
熱傳輸殼（DS）
第二項：節點數
第三項：減縮積分（可選）
第四項：
可選五自由度
耦合溫度-位移（T）
ABAQUS/Explicit中的小應變公式（S）
第五項：
在ABAQUS/Explicit中的小應變公式中考慮翹曲（可選）
例如，S4R是一個四節點四邊形應力/位移殼單元，具有減縮積分和大應變公式；SC8R是一個八節點、四邊形、一階插值、應力/位移連續體殼單元，具有減縮積分。
展開 ABAQUS 單元類型和選型規則
為了方便ABAQUS 的初學者，下面對ABAQUS單元類型選型規則匯總為如下表格
ABAQUS單元類型和選擇規則.png
ABAQUS顯式子程序調用規則及nblock變量解釋
在使用ABAQUS進行有限元仿真時，調用用戶子程序可以使用戶解決一些問題時具有很大的靈活性。在實現更加復雜的計算過程時，ABAQUS軟件與子程序之間的交互關系是需要了解的。Vxxxxxx(VUMAT/VUHARD/VUSDFLD等)子程序內需要將用戶自定義變形部分寫入(do k=1,nblock → end do)語句內，理解nblock的含義將對子程序的改進大有幫助。本文以VUHARD子程序為例，測試了ABAQUS顯式算法在調用用戶子程序時的基本規則。
首先建立了三維Part，為其賦予熱變形本構參數（用戶可參考Johnson-Cook模型），單元采用C3D8RT單元，具有8個節點，1個積分點。固定Z0面，施加沿Z+方向的位移載荷。
對模型劃分成不同網格數量進行測試，測試結果下表所示。
單元數
核數
nblock
k
8
1
8
1-8
64
1
64
1-64
125
1
125
1-125
216
1
136；80
1-136；1-80
343
1
136×2；71
(1-136)×2；1-71
1000
1
136×7；48
(1-136)×7；1-48
統計結果表明，nblock的最大值為136。
展開 abaqus用fortran寫子程序的語法規則
abaqus用fortran寫子程序的語法規則
abaqus用fortran寫子程序時個人覺得語法上偏FORTRAN77，因此，收集了下的FORTRAN 77 語法概述。
以下為轉載內容：http://blog.sina.com.cn/s/blog_6e0aa6420100pu9p.html
Abaqus中多相材料不規則表面輪廓提取
1 前言
在Abaqus切削仿真中，目前多采用二維正交模型來轉化和代替各種加工形式。目前對于Abaqus切削仿真的可查資料中，多是模型建立和一些注意事項，對于其后處理過程較少提及。加工表面的粗糙度是表面質量評價的一項重要指標，仿真得到的微觀結構的細觀變化也是切削仿真的一大研究重點。因此，對切削表面的輪廓提取是有必要的。
在Abaqus中，在后處理過程中，輪廓提取可以采用多種方式，例如建立路徑和導出連續節點坐標。但是多相材料建模通常采用不同的Part最后Assembly得到，建立路徑只能在獨立的Part中進行。除此之外，網格劃分與一般的規則形狀得到的網格也不相同，多相材料劃分得到的網格往往并不規則，因此導出連續節點也是不現實的。因此，要想導出多相復合材料的表面輪廓需要尋求一種別的方式。
本篇小節只要針對多相復合材料的切削表面輪廓進行講述，所使用的軟件包括Abaqus、AutoCAD、Excel和Origin。除了Abaqus切削仿真表面輪廓提取，也可對一些其它復雜形狀和結構的輪廓提取做出指導意義。
2 提取過程
如圖1所示為所選擇的一個案例的最終切削結果，其中顆粒和基體是單獨的Part，切削表面并不平整。提取目標既最上面的一條輪廓曲線。要求：其中某一條線的實際長度。
圖1 樣件的切削結果
第一步：在Abaqus Visualization界面下選擇命令Plot Deformed Shape，再點擊命令Render Model: Wireframe，得到如圖2所示的線框圖。
展開 abaqus應用之收斂篇 ￥1.66
<h1><strong>一、收斂的定義和重要性</strong></h1><h2><strong style="color: rgb(51, 51, 51);">1．收斂的多種含義</strong></h2><p>在有限元分析中，收斂具有多重意義。它包括網格收斂、時間積分精度和非線性程序收斂。</p><p><strong>l 網格收斂</strong>是指增加模型單元數量會使仿真解趨于解析解。對于線性和非線性問題都適用，AbaqUS 中使用 H 網格自適應技術來輔助實現網格收斂。當進一步加密網格時，結果變化很小或不變時，可認為網格達到收斂。但也存在一些例外情況，如網格奇異解或材料損傷累積在模型特定區域的局部問題。</p><p><strong>l 時間積分精度</strong>則是針對具有物理時間尺度的瞬態問題，AbaqUS 提供用戶定義參數，以控制對相關方程的積分精度。</p><p><strong>l 非線性程序收斂</strong>是本文重點討論的內容，要獲得精確解需要滿足網格收斂、瞬態問題的精確時間積分以及非線性求解過程收斂等條件。</p><h2><strong style="color: rgb(51, 51, 51);">2．收斂對分析結果的影響</strong></h2><p>收斂性直接關系到分析結果的準確性。如果模型不收斂，得到的結果可能毫無意義，甚至會誤導工程決策和學術研究。因此，理解和掌握 ABAQUS 中的收斂問題是正確使用該軟件進行有效分析的基礎。
展開 Abaqus有限元分析不收斂該怎么辦？ 附Abaqus 分析用戶手冊材料卷文檔下載
ABAQUS提供式和隱式兩種求解類型，其中顯式計算方法是“有條件收斂的”，只需滿足增量步小于限值，大多數情況均能順利計算完成；而隱式計算方法，在非線性情況下極易出現不收斂的情況，比如：欠約束、接觸、材料塑性或失效、斷裂、屈曲失穩等，都可能導致多次迭代不收斂，增量步大小一降再降，直到滿足終止條件而退出計算。作為老司機，使用了這么多年的軟件總有點心得吧，總結了五條經驗，分享給大家：
一、ABAQUS的任務提交流程
了解ABAQUS的任務提交流程，也就是讓我們學會找錯！當我們點擊Submit后會有兩個處理階段：1）預處理；2）任務計算。
結合ABAQUS Job Monitor窗口進行講解，兩個階段的分界點位于Data File子頁面的內容是否完成；也就是說，當出現Error，而Message File和Status File未激活（生成）時，表明還處于預處理階段，我們定義的模型一開始就存在問題，Errors子頁面都會一一列出，通常會有：信息不完整、材料參數不符合本構模型、特殊定義之間沖突、關鍵字輸入問題等，我們只需逐個修改即可。
當順利進入任務計算階段后，窗口上方的表格將實時更新為Status File（jobName.sta）中的內容，提示計算的進度，當后續再出現Error時，才可能是由于計算不收斂導致的錯誤。
二、收斂的基本條件
模型收斂是什么？很多初學者估計都不太清楚，從而提出一些奇奇怪怪的問題，比如：“我用彈性材料可以計算，換成復雜材料模型就計算不了，為什么？”。
展開 abaqus四點彎曲不收斂
Abaqus/Standard Analysis exited with an error - Please see the message file for possible error messages if the file exists.
Path based tracking is defined in contact pair (assembly__pickedsurf43,assembly_part-2-1_rigidsurface_). Path based tracking cannot be used with analytical rigid master surfaces, the state based tracking algorithm will be used instead.
Path based tracking is defined in contact pair (assembly__pickedsurf45,assembly_part-2-2_rigidsurface_). Path based tracking cannot be used with analytical rigid master surfaces, the state based tracking algorithm will be used instead.
Solver problem. Zero pivot when processing D.O.F. 1 of 1 nodes. The nodes have been identified in node set WarnNodeSolvProbZeroPiv_1_1_5_5_1.
展開 Abaqus不收斂怎么辦？
這里推薦大家有能力的可以深入閱讀ABAQUS幫助文檔分析手冊（Analysis Solution and Control）的內容,會對解決不收斂的問題有很大幫助。在2021版幫助文檔中入口如下圖所示
Abaqus中的求解類型分為顯示與隱士。顯示求解基于動力學方程，當前求解至于前一時刻的速度與位移有關，求解過程無需迭代，是有條件穩定，無條件收斂；隱式求解則基于虛功原理，一般需要迭代計算，無條件穩定，有條件收斂。
所以對于你提到的不收斂問題這里我就認為是針對隱式計算。
1.基礎不收斂問題
在模型計算報錯時，我們首先查看Job Monitor中的報錯信息，對于簡單的模型錯誤，例如材料、邊界、載荷定義錯誤、網格問題、關鍵字定義錯誤等都能在Job Monitor中直接看出。針對不同的問題針對修改就可以。這部分相信稍微有些經驗的CAEer都能自行解決。
2.不收斂的本質與進階解決方法
在接著講之前，我希望大家能夠了解模型收斂的本質是什么？
對于線性系統來說，一般不存在收斂問題。模型不收斂一般都是由于幾何非線性、材料非線性、邊界非線性。
有限元求解的過程是根據外力與內力平衡求解出各個節點的位移，根據位移再求解應力、應變等。對于非線性系統，載荷與位移的關系也通常是非線性的，如下圖所示。
我們的首要目的就是求解不同載荷下的位移。做法就是將一個完整求解過程細分為許多個小的過程。
這里就不得不提abaqus中的Step、increment與iterations。
這里舉一個例子，假如我們從冰箱里拿面包吃。
展開 abaqus 復合材料接觸不收斂
abaqus 復合材料接觸不收斂
第二章 abaqus分析收斂準則
abaqus分析收斂準則（外文書籍翻譯）.pdf
原文書籍：Troubleshooting Finite-Element Modeling with Abaqus With Application in Structural Engineering Analysis by Raphael Jean Boulbes (z-lib.org)
第二章 abaqus分析收斂準則
2.1 收斂問題的癥狀
收斂問題是與工程設計相關的一個典型的分析問題，涉及撓度、位移、應力、固有頻率、溫度分布等的預測。這些參數用于迭代材質參數和/或幾何體以優化其行為。傳統的方法，如手工計算，涉及理想化的物理模型使用簡單的方程來獲得解決方案。然而，這些近似使問題過于簡單化，而解析解只能提供保守估計。或者，有限元法和其他數值方法旨在提供一個考慮到更多細節的工程分析，這對于手工計算是不實際的。有限元法將物體分割成小塊，使位移沿這些單元邊界連續。對于那些使用有限元分析的人，通常使用收斂項。大多數線性問題不需要迭代求解過程。網格收斂是一個重要問題另外，還有在迭代過程中也需要考慮收斂性。在本節中，將調查收斂性問題并解決與此術語相關的問題。首先，要識別大多數收斂問題的癥狀，可以在消息文件（.msg）擴展名中找到。此外，.dat和（.sta）文件也可能包含問題的癥狀。有一些共同的信息可能表明收斂問題在求解有限元模型時造成數值困難。
展開 
Abaqus有限元分析不收斂該怎么辦？ 附ABAQUS非線性有限元分析實例下載
ABAQUS提供式和隱式兩種求解類型，其中顯式計算方法是“有條件收斂的”，只需滿足增量步小于限值，大多數情況均能順利計算完成；而隱式計算方法，在非線性情況下極易出現不收斂的情況，比如：欠約束、接觸、材料塑性或失效、斷裂、屈曲失穩等，都可能導致多次迭代不收斂，增量步大小一降再降，直到滿足終止條件而退出計算。作為老司機，使用了這么多年的軟件總有點心得吧，總結了五條經驗，分享給大家：
一、ABAQUS的任務提交流程
了解ABAQUS的任務提交流程，也就是讓我們學會找錯！當我們點擊Submit后會有兩個處理階段：1）預處理；2）任務計算。
結合ABAQUS Job Monitor窗口進行講解，兩個階段的分界點位于Data File子頁面的內容是否完成；也就是說，當出現Error，而Message File和Status File未激活（生成）時，表明還處于預處理階段，我們定義的模型一開始就存在問題，Errors子頁面都會一一列出，通常會有：信息不完整、材料參數不符合本構模型、特殊定義之間沖突、關鍵字輸入問題等，我們只需逐個修改即可。
當順利進入任務計算階段后，窗口上方的表格將實時更新為Status File（jobName.sta）中的內容，提示計算的進度，當后續再出現Error時，才可能是由于計算不收斂導致的錯誤。
二、收斂的基本條件
模型收斂是什么？很多初學者估計都不太清楚，從而提出一些奇奇怪怪的問題，比如：“我用彈性材料可以計算，換成復雜材料模型就計算不了，為什么？”。
展開 ABAQUS-接觸分析中收斂問題的解決方法【轉載】
5、平穩地建立接觸關系
如果在第一個分析步中就把全部載荷施加到模型上，使接觸狀態的發生劇烈改變，會有可能造成收斂的困難。因此一般應首先定義一個只有很小載荷的分析步，讓接觸關系平穩地建立起來，然后在下一個分析步中再施加真實的載荷。盡管這樣可能需要更多的分析步，但這減小了收斂的困難，會提高求解的效率。
6、細化網格
細化從面和主面的網格是解決收斂問題的一個重要方法。過于粗糙的網格會使ABAQUS難以確定接觸狀態，例如，如果在接觸面的寬度方向上只有一個單元，則常常會出現收斂問題。一般來說，如果從面上有90°的圓角，建議在此圓角處至少劃分10個單元。
7、使用一階單元
如果接觸屬性為默認的“硬”接觸，則不能使用六面體二次單元（ C3D20和C3D20R)，以及四面體二次單元（C3D10)，而應盡可能使用六面體一階單元。如果無法劃分六面體單元網格，可以使用修正的四面體二次單元（C3D10M)。
8、正確定義主面和從面
1）選擇剛度較大、網格較粗的面作為主面。
2）主面在發生接觸的部位不要有尖角或大的凹角。
3）如果是有限滑移，則在整個分析過程中，都盡量不要讓從面節點落到主面之外。
4）如果主面和從面在幾何位置上沒有發生重疊，則一個面的法線應指向另一個面所在的那一側(對于三維實體，法線應該指向外側)。
9、避免過約束
如果在節點的某個自由度上同時定義了兩個以上的約束條件，就會發生所謂“過約束”( overconstraint )。可能造成過約束的有以下主要因素。
1）接觸:從面節點會受到沿主面法線方向的約束。
2）邊界條件。
展開 ABAQUS彈塑性收斂問題
ABAQUS三維盾構隧道，莫爾庫侖本構。地應力平衡分析出現以后結果不收斂：The plasticity/creep/connector friction algorithm did not converge at 129415 points 去掉塑性參數后就可以收斂。請教大神這是什么問題
ABAQUS不收斂的原因一
一般說來,Mohr-Coulomb相對難收斂些,因為它在主應力空間的屈服面上存在尖角——因為在計算過程中需要對屈服函數求導,而如果曲線不光滑,尖角處也就沒法求導,向后歐拉算法也就沒法進行了。理論上基本是不能算的,但Abaqus內部對尖角處進行了一定的處理,替代為光滑的函數,這樣收斂性就得到了一定的保證,雖然還不夠好。而Drucker-Prager準則就較好地克服了這一點。因此,相對來說,建議考慮使用DP模型。尤其是ABAQUS中有非線性的DP模型,可以在一定程度上克服線性DP在剪拉區面積過大的情況,因此可以在一定程度上減小開挖卸載時土體的回彈。
一般材料相關的不收斂提示基本就是類似以下:
The plasticity/creep/connector friction algorithm did not converge at 364 points
碰到這個錯誤,一般直接原因就是材料應變太大了,塑性計算迭代不收斂。但這并不意味著材料參數給得不合適或材料強度太弱,很有可能是你的接觸、約束、荷載或邊界出了問題,導致計算中出現的特別大的位移。或者是初始條件(如初始地應力)出了問題,程序沒有計算,直接就報了這個錯誤。可以在下面位置看到計算不收斂的單元位置:
后處理 -> Tools->Job Diagnostics
不收斂需要單一調每個參數確定哪里有問題。
另外，在使用Soils分析步進行孔壓—應力分析時:
a. 墻裂建議各位注意單位,應力單位最好選擇kPa或MPa,不要用Pa,否則可能會遇到各種無腦錯誤提示。
b. 墻裂建議給定孔壓邊界,否則會遇到DOF. 8極大的情況(孔壓就是第8自由度)。
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
展開 