ansys非線性收斂準則的案例
提高ANSYS非線性求解收斂性能的一般方法總結
提高ANSYS非線性求解收斂性能的一般方法總結
在采用ANSYS進行幾何非線性的求解過程中，如果采用系統的默認設置，有可能會因為參數的不合適而導致收斂困難。針對比較常見的非線性求解收斂困難，本文總結幾種比較常見的調整方法，僅做參考。
一、打開自動時間步(autots，on)
ANSYS在所有靜態和瞬態分析中，使用時間作為跟蹤參數，而不論分析是否依賴于時間。當我們收斂困難時，一個非常重要的方法是打開自動時間步。打開自動時間步長往往需要一個比較小的最小時間步長(或者最大的步數，采用DELTIM或者NSUBST定義)。在進行非線性求解過程中，初始時間步長如果太小，自動時間分步算法可能使你的運行時間太長；相反地，使你的最小時間步長太大，可能導致不收斂。因此合理設置初始時間步長是非常重要的，一般可以根據試算確定。
值得說明的是，當采用自動時間步長后，二分法會被自動激活。如果在一個太大時間步內收斂失敗，該特性能減半時間步，使得計算繼續。
二、Newton-Raphson 選項調整與自適應下降
在非線性計算收斂困難時，可以調整程序默認的牛頓—拉普森選項。雖然一般情況下采用采用自動默認的選項會獲得最佳的收斂特性，但也不排除會遇到其他選擇會更有效的情況。合理使用自適應下降因子也能增加某些非線性單元的收斂情況。
三、線性搜索
線性搜索可以看為是對自適應下降的一個替代，兩者不應同時使用。線性搜索一般情況下能使得分析得到收斂，但隨之帶來的是求解速度的大幅度降低，特別是針對有材料非線性的情況，一般而言，當結構采用力加載或者剛度增長的薄膜等類似分析時，可以打開線性搜索，線性搜索打開命令為LNSRCH。
四、調整收斂準則
ANSYS非線性收斂準則主要有四種，分別為力、位移、彎矩和轉角。
展開 ansys非線性收斂總結
文章來源于網絡，講解很系統，可以經典收藏，由于無法查證出處，無意冒犯，如有不妥，請聯系我
ansys非線性收斂總結
ansys計算非線性時會繪出收斂圖，其中橫坐標是cumulative iteration number 縱坐標是absolute convergence norm。他們分別是累積迭代次數和絕對收斂范數，用來判斷非線性分析是否收斂。
ansys在每荷載步的迭代中計算非線性的收斂判別準則和計算殘差。其中計算殘差是所有單元內力的范數，只有當殘差小于準則時，非線性疊代才算收斂。
ansys的收斂是基于力的收斂的,以力為基礎的收斂提供了收斂量的絕對值，而以位移為基礎的收斂僅提供表現收斂的相對量度。一般不單獨使用位移收斂準則,否則會產生一定偏差,有些情況會造成假收斂.(ansys非線性分析指南--基本過程Page.6) 。因此ansys官方建議用戶盡量以力為基礎（或力矩）的收斂誤差，如果需要也可以增加以位移為基礎的收斂檢查。
ANSYS缺省是用L2范數控制收斂。其它還有L1范數和L0范數，可用CNVTOL命令設置。在計算中L2值不斷變化，若L2<criterion的時候判斷為收斂了。也即不平衡力的L2范數小于設置的criterion時判斷為收斂。
由于ANSYS缺省的criterion計算是全部變量的平方和開平方（SRSS）*valuse(你設置的值)，所以crition也有小小變化。如有需要，也可自己指定crition為某一常數， CNVTOL,F,10000,0.0001,0就指定力的收斂控制值為10000*0.0001=1。
展開 ANSYS Workbench非線性分析收斂曲線解讀
進行非線性分析時，收斂性是大家非常關心的一個問題。在Ansys workbench中，可以通過Details of “Solution Information”中選擇“Solution Output=Force Convergence”來查看收斂情況，其中，最直觀的莫過于力收斂曲線了。
Solution Output選項
力收斂曲線如下圖所示:
力收斂曲線圖
判斷收斂的方法很簡單，只要“計算的力收斂曲線”落在“力收斂準則”曲線之下，就表示該載荷步或子步收斂了。
該模型中有兩個載荷步，分析設置中時間步長設置為“Program Contrlled”.
除了看上述的力收斂曲線圖，我們可以設置“Solution Output= Solve Output”查看計算輸出信息，從其中可以更詳細地看到收斂情況。
可以將計算輸出的信息與力收斂曲線圖對比起來看，就更容易理解力收斂圖了。
第1個載荷步中，第1個分析子步經過了15次迭代收斂（圖中每個圓點代表一次迭代）。
經過4個分析子步，第1個載荷步完成加載并收斂。第2個載荷步程序自動設置的信息如下：
初始子步數量為5，載荷步的分析時間為1s，因此初始的時間步長為0.2s。
第2個載荷步的第1個分析子步，經過25次計算迭代后，還不收斂。程序進行自動二分，將時間步長除以2，變為0.1s。
自動二分是一種用于解決非線性分析過程中收斂困難的策略。當收斂失敗發生在某個子步中，程序會自動減小時間步長，通常是前一個步長的一半左右。然后，程序會從前一個成功收斂的時間子步繼續求解。如果再次遇到收斂失敗，程序會繼續減小時間步長并繼續求解，直到達到收斂或達到指定的最小時間步長值。這種方法有助于逐步逼近正確解，并確保分析的穩定性和準確性。
展開 ANSYS非線性計算的收斂和速度
A：調整優化非線性計算的收斂和速度可以說幾乎是一種藝術，即沒有固定的可循規則。
我的經驗是，你的結構的"非線性"越小，非線性的變化越規則，就越容易收斂。 想象一下如果你是手算這個非線性問題，對你來講較容易的，對ANSYS的相應算法也會容易些。
可以把你的地震時程分析拿出幾點，做一下靜態的非線性分析，同時調整模型看看分析出來的結果是否合理。如果這一步還沒有做，那花大量時間做出的時程分析是廢品的可能性十分之大。
一定要記住有限元分析是一個"簡化"問題的過程。建立一個模型一定要由淺到深，線性的模型沒有搞透不要貿然進攻非線性，靜態沒有搞透不要碰時程分析。
A：影響非線性收斂穩定性及其速度的因素很多，我們可以看看這幾點：
1、模型——主要是結構剛度的大小。對于某些結構，從概念的角度看，我們可以認為它是幾何不變的穩定體系。但如果結構相近的幾個主要構件剛度相差懸殊，或者懸索結構的索預應力過?。此膭偠炔粔虼螅?，在數值計算中就可能導致數值計算的較大誤差，嚴重的可能會導致結構的幾何可變性——忽略小剛度構件的剛度貢獻。
有一種通用的方法判斷結構的幾何可變性，即det(K)=0。
在數值計算中，要得到det(K)恒等于零是不可能的，也就只能讓它較小時即認為結構是幾何可變的。
對于上述的結構，他們的K值是很小的，故而也可判斷為幾何可變體系。事實上這類結構在實際工程中也的確是非常危險的。
為此，看模型有沒有問題。如出現上述的結構，要分析它，就得降低剛度很大的構件單元的剛度，可以加細網格劃分，或著改用高階單元（BEAM->SHELL,SHELL->SOLID)。構件的連接形式（2剛接或鉸接）等也可能影響到結構的剛度。
2、線性算法（求解器）。
展開 
ANSYS Mechanical 非線性結構分析的收斂性 ￥5
ANSYS Mechanical 作為ANSYS致力于結構分析的模塊，可以對線性以及非線性結構分析問題進行仿真。其中非線性問題對于 用戶都是一種挑戰，分析過程中頻繁蹦出的“errors”and“Warnings”挑戰著分析人員的耐心，結果收斂成為大家最期待的結果。如果想順利進行非線性結構分析，學會診斷不收斂問題，就顯得至關重要了。
Ansys影響非線性收斂穩定性及其速度的因素分析
3非線性逼近技術。
在ANSYS里還是牛頓－拉普森法和弧長法。牛頓－拉普森法是常用的方法，收斂速度較快，但也和結構特點和步長有關?；￠L法常被某些人推崇備至，它能算出力加載和位移加載下的響應峰值和下降響應曲線。但也發現：在峰值點，弧長法仍可能失效，甚至在非線性計算的線性階段，它也可能會無法收斂。
為此，盡量不要從開始即激活弧長法，還是讓程序自己激活為好（否則出現莫名其妙的問題）。子步（時間步）的步長還是應適當，自動時間步長也是很有必要的。
4加快計算速度
在大規模結構計算中，計算速度是一個非常重要的問題。下面就如何提高計算速度作一些建議：
充分利用ANSYS MAP分網和SWEEP分網技術，盡可能獲得六面體網格，這一方面減小解題規模，另一方面提高計算精度。
在生成四面體網格時，用四面體單元而不要用退化的四面體單元。比如95號單元有20節點，可以退化為10節點四面體單元，而92號單元為10節點單元，在此情況下用92號單元將優于95號單元。
選擇正確的求解器。對大規模問題，建議采用PCG法。此法比波前法計算速度要快10倍以上（前提是您的計算機內存較大）。對于工程問題，可將ANSYS缺省的求解精度從1E-8改為1E-4或1E-5即可。
5荷載步的設置直接影響到收斂。
展開 關于ANSYS Workbench非線性分析收斂的學習筆記 ￥5
我們通過用應力——應變關系來描述材料的受力響應特征，如果某種材料的應力——應變關系呈現線性關系，那么我們基本就可以通過彈性模型模量以及泊松比來定義其材料性能，但是我們涉及的很多材料其應力-應變關系并不是非線性的，比如超彈材料，彈塑性材料、蠕變和黏彈性材料等，其本構關系需要通過一些本構模型來擬合表征，往往變得比較復雜。與此同時，如果在有限元分析中應用到非線性材料，無疑會對材料性能的確定、計算的設置、計算的收斂以及保證結果的精度增加難度。</li><li>幾何非線性。我們日常生活中，某些結構在載荷變化的過程中會發生突變，結構發生大變形，這將直接導致結構的響應規律發生重大變化，其中主要分為兩種，分別是大撓度和大應變。（1）大撓度：比如向上翹曲的殼零件，在很小載荷作用下發生很小的應變和位移，但當載荷作用加大，殼會向下凹陷，也就變成了小應變，大位移。（2）大應變:比如橡膠件在壓力作用下發生的變形，橡膠幾乎是不可壓縮的，應變很小，但是在拉伸時，應變很大。</li><li>狀態非線性。絕大部分有限元分析都不是簡單的零件分析，而且復雜的裝配體分析，很多的零件之間會存在接觸或者分離的狀態變化。比如齒輪的嚙合，兩個齒輪會存在接觸和分離的狀態變化，這時候結構剛度就會因為狀態變化而變化。邊界條件中的接觸就是狀態非線性的一種。</li></ol><p> 非線性無疑會增加有限元分析的難度和成本，在我們實際的模型中，以上三種非線性類型往往交叉出現，不僅具有材料非線性、幾何非線性，還有狀態非線性。對于此類問題，新手往往難以完成計算設計以實現計算收斂，下面簡單介紹相關非線性計算的收斂技巧。</p>
展開 干貨 | 接觸非線性應用——解決ANSYS 接觸不收斂問題的方法
根據ANSYS的使用者反饋，針對非線性接觸問題上的求解，經常會有客戶出現不收斂的情況，在調試收斂性上花費大量的時間。本文主要針對ANSYS 接觸不收斂問題進行方法上的技巧總結，希望通過本文使大家在ANSYS軟件的使用上有更好的體驗。
ANSYS接觸不收斂的原因有非常多的原因，針對每一種不收斂問題，選擇正確的方法都能使不收斂問題解決變得容易起來。在使用軟件中，ANSYS接觸不收斂原因主要有下面這些原因：
1、接觸算法的不正確選擇；
2、遺漏了相關的接觸對；
3、物體之間接觸剛度過大；
4、求解的載荷步較少；
5、奇異；
6、結構發生了剛體位移；
7、結構發生振蕩現象；
下面針對這些原因的解決辦法進行詳細的講解：
1
接觸算法的選取原則
ANSYS內部大體上包括5種算法，Pure Penalty，Augmented Lagrange，MPC，Pure Lagrange，Beam。
展開 Ansys非線性不收斂10大對策：讓你有“跡”可循，有“法”可醫
因此，接觸是強非線性問題。仿真中，80%以上的非線性不收斂主要是由于接觸問題引起的。
圖5 梁發生接觸后，結構剛度變大
非線性分析方程求解
ANSYS程序的方程求解器通過計算一系列的聯立線性方程組來預測工程系統的響應。然而非線性結構的行為，不能直接由線性方程求得，一種近似的非線性求解是將載荷分成一系列的載荷增量，可以在幾個載荷步內或者在一個載荷步的幾個子步內施加載荷增量，每一個增量確定一個平衡條件，在每一個增量的求解完成后，程序調整剛度矩陣以反映結構剛度的非線性變化。這種方法我們把它叫做牛頓-拉普森平衡迭代法，或者直接叫做牛頓迭代法。
其中：KiT為切向剛度矩陣；ΔUi為位移增量；F a是施加的載荷矢量；Fi nr為內力矢量。
下圖是一個載荷增量的迭代求解過程：
圖6 牛頓-拉普森平衡迭代過程
第一次迭代施加總載荷F a，對應的位移結果為X1，根據位移X1，計算內力F 1，若是 Fa≠ F1，系統不收斂，將進行剛度矩陣的修正，然后進行第二次迭代求解，第三次迭代……直至收斂。其中的差值Fa-Fi即外力與內力的偏差，也叫殘差力，殘差力需要足夠?。‵a≈F1 ，即內外力平衡）才能夠收斂，ANSYS程序中有相關的收斂準則定義。
非線性不收斂原因及ANSYS解決方案
ANSYS Mechanical具有強大的非線性計算能力，能夠對幾何非線性、材料非線性、接觸非線性、混合非線性等計算問題進行非常好的模擬仿真，是目前最強大的非線性問題計算軟件之一。針對非線性計算無法收斂的問題，我們主要可從以下方面著手：
l 首先從solution information中尋找突破點，找出報錯原因。
展開 ABAQUS非線性分析的平衡迭代過程和收斂原則 附ABAQUS非線性有限元分析與實例下載
12、如果滿足限制要求，則迭代結束，該增量步計算收斂并結束。如果不滿足要求，則迭代繼續。同時，迭代步不是無限制地增加，當超過16次迭代時，迭代將停止，增量步將不收斂。
下載地址：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中有時也采用’軟”接觸代替硬接觸來解決模擬中的數值收斂性問題。
那我們的問題是，這個‘有時’是指什么時候呢？
展開 
解決非線性分析不收斂的技巧
影響非線性收斂穩定性及其速度的因素很多：
1、模型——主要是結構剛度的大小。對于某些結構，從概念的角度看，可以認為它是幾何不變的穩定體系。但如果結構相近的幾個主要構件剛度相差懸殊，在數值計算中就可能導致數值計算的較大誤差，嚴重的可能會導致結構的幾何可變性——忽略小剛度構件的剛度貢獻。 如出現上述的結構，要分析它，就得降低剛度很大的構件單元的剛度，可以加細網格劃分，或著改用高階單元（BEAM->SHELL,SHELL->SOLID)。構件的連接形式（剛接或鉸接）等也可能影響到結構的剛度。
2、線性算法（求解器）。ANSYS中的非線性算法主要有：稀疏矩陣法(SPARSE DIRECT SOLVER)、預共軛梯度法(PCG SOLVER)和波前法(FRONT DIRECT SLOVER)。稀疏矩陣法是性能很強大的算法，一般默認即為稀疏矩陣法（除了子結構計算默認波前法外）。預共軛梯度法對于3-D實體結構而言是最優的算法，但當結構剛度呈現病態時，迭代不易收斂。為此推薦以下算法：
1）、BEAM單元結構，SHELL單元結構，或以此為主的含3-D SOLID的結構，用稀疏矩陣法；
2）、3-D SOLID的結構，用預共軛梯度法；
3）、當你的結構可能出現病態時，用稀疏矩陣法；
4）、當你不知道用什么時，可用稀疏矩陣法。
3、非線性逼近技術。在ANSYS里還是牛頓－拉普森法和弧長法。牛頓－拉普森法是常用的方法，收斂速度較快，但也和結構特點和步長有關?；￠L法常被某些人推崇備至，它能算出力加載和位移加載下的響應峰值和下降響應曲線。但也發現：在峰值點，弧長法仍可能失效，甚至在非線性計算的線性階段，它也可能會無法收斂。
為此，盡量不要從開始即激活弧長法，還是讓程序自己激活為好（否則出現莫名其妙的問題）。子步（時間步）的步長還是應適當，自動時間步長也是很有必要的。
展開 非線性不收斂原因及workbench解決方案
因此，接觸是強非線性問題。仿真中，80%以上的非線性不收斂主要是由于接觸問題引起的。
圖5 梁發生接觸后，結構剛度變大
4、非線性分析方程求解
ANSYS程序的方程求解器通過計算一系列的聯立線性方程組來預測工程系統的響應。然而非線性結構的行為，不能直接由線性方程求得，一種近似的非線性求解是將載荷分成一系列的載荷增量，可以在幾個載荷步內或者在一個載荷步的幾個子步內施加載荷增量，每一個增量確定一個平衡條件，在每一個增量的求解完成后，程序調整剛度矩陣以反映結構剛度的非線性變化。這種方法我們把它叫做牛頓-拉普森平衡迭代法，或者直接叫做牛頓迭代法。
其中：
KiT為切向剛度矩陣；ΔUi為位移增量；F a是施加的載荷矢量；Fi nr為內力矢量。
下圖是一個載荷增量的迭代求解過程：
圖6 牛頓-拉普森平衡迭代過程
第一次迭代施加總載荷F a，對應的位移結果為X1，根據位移X1，計算內力F 1，若是 Fa≠ F1，系統不收斂，將進行剛度矩陣的修正，然后進行第二次迭代求解，第三次迭代……直至收斂。其中的差值Fa-Fi即外力與內力的偏差，也叫殘差力，殘差力需要足夠小（Fa≈F1 ，即內外力平衡）才能夠收斂，ANSYS程序中有相關的收斂準則定義。
二、非線性不收斂原因及ANSYS解決方案
ANSYS Mechanical具有強大的非線性計算能力，能夠對幾何非線性、材料非線性、接觸非線性、混合非線性等計算問題進行非常好的模擬仿真，是目前最強大的非線性問題計算軟件之一。針對非線性計算無法收斂的問題，我們主要可從以下方面著手：
1、首先從solution information中尋找突破點，找出報錯原因
。
展開 幾何非線性不收斂問題要怎么辦
做幾何非線性屈曲分析，一直不收斂，實在是頭疼
非線性仿真之如何解決接觸仿真收斂問題
遇到接觸仿真無法收斂怎么辦？是不是嘗試增加更多子步緩慢加載還是無法解決收斂性問題？這篇文章給大家介紹一些關于非線性分析和收斂的重要背景知識，并討論克服收斂問題的不同方法。具體涉及接觸分析時，可以嘗試以下幾種方法來幫助提升計算收斂性：
1.消除剛體運動：
a. 開始時讓裝配體中的所有部件都相互接觸。這可以通過移動部件、添加接觸偏移量或添加穩定阻尼來實現。
b. 在接觸表面添加摩擦，避免切向毫無阻力的運動。
2.克服不收斂：
a. 降低接觸單元的剛度。（我的經驗表明，緩慢增加載荷和降低接觸剛度可以解決90%的收斂問題）。
b. 在接觸區域細化網格，以減少反復進入和脫離接觸的單元百分比。
在本文中，我將使用一個具體的例子來演示上述一些方法，并描述其他幾種有助于克服頑固的與接觸相關的收斂問題的方法。需要注意的是，許多有限元分析程序（如ANSYS）都內置了接觸算法，試圖設置程序默認值以實現快速收斂和準確的解。然而，不可能設計一種萬能的接觸算法，使其在每種接觸條件下都能自動工作。它們是為解決常見情況而設計的，但在某些情況下可能需要手動干預。
在這個例子中，一個板彈簧被一個承受作用力的扁平剛性板壓縮，如圖1所示。這個分析使用了ANSYS Workbench有限元軟件。為了得到收斂解，需要進行幾次嘗試。
第1次計算嘗試：
圖1 計算例子
第一次嘗試求解沒有收斂，并給出以下錯誤：“內部解的大小限制被超過”。這種類型的錯誤，以及其他如“小的負方程求解器主元項”或僅僅是 “遇到求解器主元警告或錯誤”，表明存在剛體運動。
展開 