abaqus判斷收斂的案例
三十二、Fluent收斂判斷標準及方法
wx_fmt=png"> </p><p><br></p><p>并不是所有的工況都會收斂。當各物理量的值基本不變時，即殘差很小時，工況才可能收斂。</p><p><br></p><p>但是對于瞬態，各物理量的值總是變化，如何收斂？正因如此瞬態才有時間步的概念，瞬態問題在每個時間步上都認為是穩態，所以瞬態問題的殘差圖總是波浪線型。</p><p><br></p><p><br></p><p><strong>2. 收斂標準</strong></p><p> </p><p><strong>2.1 殘差標準</strong></p><p><br></p><p>對所有的工況，沒有統一的判斷標準。對于大多數問題，默認的判斷標準已經足夠（For most problems, the default convergence criterion in ANSYS Fluent is sufficient. ）</p><p><strong>建議殘差達到設定值后，多算50步，確定殘差之后都是減小的趨勢。</strong></p><p> <img src="https://mmbiz.qpic.cn/mmbiz_png/8tJMdLVYZyicRKfutZe5xPg3lTv8ADcu3WMbaibmnIkR3DeVcWx3Rmp2bmlIeuAPKZoH9IA2OwL80dj5U7qmjthQ/640?wx_fmt=png"> </p><p> </p><p> </p><p><strong>2.2 監測物理量</strong></p><p><br></p><p>為了使結果更加精確，可根據工況輔助檢測一些物理量，如速度，流量等。<strong>當所關心的物理量基本不變時，說明達到了收斂。
展開 仿真應用 | 一種更實用的應力收斂判斷方法
收斂功能
應力收斂曲線
網格細化結果
一種更實用的應力收斂判斷方法
劃分網格
插入Body Sizing，CurvatureNormal Angel設置為9度（直角90度的1/10），目的是讓倒圓角位置的網格更密；Num Cells Across Gap默認為3，本次設置為2，目的是減少模型小面位置的網格數量，因為已知應力集中不會發生在模型的小面位置，當然也可以使用默認設置。
網格設置
網格劃分結果如下：
網格狀態
等效應力結果
如下圖所示，最大應力出現倒角處，大小為204.7MPa。
等效應力
應力收斂分析
本文推薦的更實用的應力收斂判斷方法是，對比節點平均應力和節點非平均應力，如果它們之間的差距小于3%，可認為應力已經收斂，如下表所示。
這種應力收斂判斷方法更實用的理由
雖然Convergence能夠智能地定位到絕大數需要細化網格的位置，但如果出現Convergence定位錯了，那Convergence方法就會失效，如果還按照此時的求解結果作出評估，就會出現評估失誤。
如果分析者并不關心Convergence推薦的位置，而是關注自己指定的位置，那么Convergence方法也無法使用。分析者需要細化關注位置的網格，再對比節點平均應力和節點非平均應力的結果，以判斷應力的收斂解，也就是本文推薦的方法。
展開 【腳本】ABAQUS中判斷某個節點是否在某個單元集 ￥1.99
〇、問題來源
本題目來源于技術鄰ABAQUS主題下的問題
Abaqus Python語言判斷node是否在單元里？
if node1 在集合set1中，請問有這種語句嗎，小白請假不勝感激。
http://www.yqgqt.org.cn/answer/1197962
一、問題分析
這個屬于簡單的腳本判斷，該問題也得到了技術鄰大佬“藍牙”和“君莫”的回答，但大佬的回答并不能夠直接幫助小白解決問題，所以斗膽一試，還請指正。
本問題應該遍歷set-1中所有單元的節點，然后與特定節點編號進行對比，若等于則說明該節點在set-1中，否則則不在。
根據以上思路編寫腳本。
二、腳本測試
在ABAQUS中隨便創建一個模型，劃分網格后形成裝配體，在裝配體中創建一個名為“Set-1”的集合，集合內包含一部分單元，如下
運行編寫的腳本
會彈出如下窗口，要求輸入要查詢的節點ID
然后，程序開始運行，若該節點在Set-1中，則輸出如下信息
否則，輸出如下信息
以上。
最后，感謝關注與支持，如果覺得有用，不如點個關注和收藏？
最后，歡迎關注微信公眾號：ABAQUS在線，讓我們一起愉快的成長吧。
展開 abaqus系列技巧1：使用 Automatic stabilization進行分析及判斷合理性 ￥25
也歡迎加入abaqus交流群516073058進行討論研究
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 復合材料接觸不收斂
第二章 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不收斂怎么辦？
我們通常說的不收斂是發生在“嚼面包”（iteration）過程中，當軟件嘗試很多次迭代發現物體的內力與外力差值仍然超過規定誤差時就會報出“不收斂”。
為了深入了解，我們這里舉個例子。將一個物體所受外力設為P，物體內力（作用在節點上）設為I。作用在節點上的內力是由與節點相連的單元的應力產生的。
一個物體處于平衡，則內力與外力相平衡。
為了求解上式，首先施加一個小的載荷增量 。根據位移為 時結構的切線剛度 以及載荷增量 計算位移修正值 ，最后計算出此時的位移 。
此時重新計算物體內力 ，可以得到你內力與外力的差值
這里的 就是衡量每次iteration是否收斂的標準。abaqus默認 小于結構所受平均力的0.5%則認為內力與外力 相平衡，但是此時還未認定為收斂，還需要檢查 是不是小于 ( )的1%（abaqus默認值）,如果是則判定收斂，進入下一個迭代。
下一次迭代計算遇上個步驟相同，最后還是通過判斷內力與外力的殘差與來判斷該次迭代是否收斂。
那么對于復雜難以收斂的模型，如果你調整其他參數發現毫無作用時，可以適當調整判斷迭代收斂的這兩個參數：與
具體入口為：Step->Other->Genreal Solution Controls->Filed Equations
還可以調整 以及 ，但是這里作者也沒有搞明白是和含義，希望大神不吝賜教。
還可以調整增量步（increment），其中 和 是重要的兩個參數。因材料非線性、幾何非線性、摩擦產生的非單調收斂需要增加；因復雜的材料本構模型而導致雅可矩陣非對稱引起的不收斂需要增加。
展開 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大世界
旨在分享，若侵即刪.
展開 [非線性]ABAQUS收斂調整(1):接觸屬性
作者：羅元元 來源：CAETube講堂
大部分情況下我們所討論的非線性分析不收斂都來源于接觸問題導致的不收斂。
接觸這件事，建模設置上而言不外乎定義接觸對、定義接觸屬性等等（如下圖），初級選手容易輕視，不收斂的時候也不知該如何診斷，不知該如何去做有效調試。
Figure-0: 接觸屬性設置
實際上，接觸，在仿真分析中，絕對是個看似青銅實則王者級別的難題。
有一些通用的解決辦法，大家可以在幫助文件的Interaction → Contact Difficulties and Diagnostics中找到，例如初始接觸狀況、穿透、突然分離造成的局部不穩定等等，但是確實沒有一概而論的措施，更多的情況下準確的診斷以及有效的改善還是要依靠經驗的累積。
這里，有一些是筆者團隊在日常工作中所積累的一些小經驗，與大家分享。
P.S．基于經驗累積和理解，可能會有些偏差或錯誤，不足之處，還請大家指正
1. ‘軟’接觸
參考幫助文件Interaction → Contact Property Models → MechanicalContact properties Contactpressure-overclosure relationships → using “softened” contact relationship
在幫助文件中指出: 適用于接觸面有一方或者皆是單薄軟面的模擬，例如墊片、表面涂層等；在Abaqus/Standard中有時也采用’軟”接觸代替硬接觸來解決模擬中的數值收斂性問題。
那我們的問題是，這個‘有時’是指什么時候呢？
展開 