ansys計算橡膠不收斂
關注創建者：王靖雯 創建時間：2023-03-07
ansys計算橡膠不收斂的視頻教程
Ansys workbench不收斂解決方案
我們在用ansys workbench進行仿真計算時，對于大型模型，尤其是非線性計算時，經常會出現不收斂的情況。 通過調整計算子歩（substep）也沒有起到良好的效果。 那么我們應該如何操作才能使計算收斂，以得到我們的最終解呢？ 讓這次課程來告訴你答案。
￥10 1小時46分鐘 17800播放
查看
ansys計算橡膠不收斂的實例教程
ansys計算之后出現這個錯誤，這是什么原因，怎么解決
A large negative pivot value ( -1.685395134E+09 ) has been encountered
in the global assembled matrix at the UZ degree of freedom of node
2028351. This may be caused by a bad temperature-dependent material
property used in the model.
眾所周知，橡膠作為不可壓縮材料，具有很強的材料非線性，在計算過程中，往往會因為產生過大的變形導致網格畸變嚴重，求解終止。不僅如此，橡膠材料計算過程中往往還會帶有復雜的幾何非線性和接觸非線性，使求解收斂難度進一步提高，本篇文章將介紹在Workbench中如何對橡膠計算的收斂進行調試。
本實例研究橡膠的擠壓問題的收斂方法，如圖1，上下實體是結構鋼材料，中間實體為橡膠，固定約束下端實體底面以及中間橡膠體的兩個端面，同時約束橡膠體的表面，使其僅有X方向自由度，上端實體上表面施加沿-x方向的6mm位移，接觸均為摩擦接觸，摩擦系數0.1，計算結構的變形情況。
圖1
首先設置橡膠材料，然后設置接觸，將橡膠設置為接觸面，結構鋼設置為目標面，兩對接觸均為摩擦接觸，摩擦系數0.1，均設為增廣拉格朗日算法，且在高斯積分點進行接觸檢測，其余設置保持程序控制，接著對將橡膠實體進行六面體網格劃分，施加邊界條件，最后進行求解設置，采用直接求解器，打開大變形。
圖2
圖3
圖4
圖5
圖6
進行求解，經過171次迭代后程序報錯，無法收斂，力的收斂曲線如圖7。
圖7
對于橡膠這種不可壓縮材料，為了提高其收斂性，可以在分析設置中的非線性控制里，將求解器改成非對稱求解器，同時打開線性搜索功能，如圖8。
圖8
進行求解，經過111步迭代后程序報錯，求解停止，報錯提示是單元產生了高度扭曲，如圖9，對于這種問題，可以嘗試通過減小時間步長解決，但減少時間步長仍會報錯，讀者可自行嘗試。
圖9
由于橡膠表現為不可壓縮性，泊松比接近0.5，易發生體積自鎖，因此可以使用U-P單元技術解除體積自鎖，只需在橡膠體下插入命令，如圖10.
展開 ********************************
* Summary of the simulation *
********************************
The computation failed.
***************************
* Expert tool diagnostics *
***************************
Stack error
***************************
* Expert tool diagnostics *
***************************
The problem F.E.M. Task 1 has not converged.
***************************
* Expert tool Suggestions *
***************************
A serious internal error has occured (in DPGETB) Please, contact
Polyflow s.a. and send us, if possible, the mesh and data files.
We apologize for the inconvenience.
展開 *SET,RAD_I,1250/2
*SET,RAD_PANI,705
*SET,RAD_O,2450/2
*SET,RAD_PANO,1830
*SET,RAD_B,180/2
*SET,M_BOLT,160
*SET,M_BH1,350+170+50
*SET,M_BH2,350+170+250
*SET,M_NUT,280/2
*SET,H_NUT,170
*SET,RAD_DRILL,2140/2
*SET,RAD_M,750+150
*SET,H_M,18
*SET,RAD_BI,1950/2
*SET,RAD_BO,2330/2
*SET,RAD_BH,5
*SET,F_RAD,45
*SET,FH_RAD,1855/2
*SET,R_FILLT,125
*SET,N,15
*SET,KEY_W,290
*SET,KEY_H,130
*SET,TH,200
*SET,FLANG_H1,350
*SET,FLANG_H2,350
*SET,H_SHAFT,350+490
*SET,H_SH,H_SHAFT+500
*SET,PI,ACOS(-1)
*SET,ELEMSIZE,60
/GRAPHICS,POWER
!*************************************
! 加載參數
!*************************************
*SET,DISP_B,0.75
*SET,ZMAX,H_SH
*SET,F_EXT,22.9215E+5*12
*SET,F_EEF,30.827E+5*12
*SET,R_OUT,RAD_PANO
*SET,T_W,375
*SET,U_R,R_OUT-T_W
!*************************************
! 準備建立模型
展開 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里還是牛頓－拉普森法和弧長法。牛頓－拉普森法是我們常用的方法，收斂速度較快，但也和結構特點和步長有關。弧長法常被某些人推崇備至，它能算出力加載和位移加載下的響應峰值和下降響應曲線。但也發現：在峰值點，弧長法仍可能失效，甚至在非線性計算的線性階段，它也可能會無法收斂。
為此，盡量不要從開始即激活弧長法，還是讓程序自己激活為好（否則出現莫名其妙的問題）。子步（時間步）的步長還是應適當，自動時間步長也是很有必要的。
A：如何加快計算速度
在大規模結構計算中，計算速度是一個非常重要的問題。下面就如何提高計算速度作一些建議：
充分利用ANSYS MAP分網和SWEEP分網技術，盡可能獲得六面體網格，這一方面減小解題規模，另一方面提高計算精度。
在生成四面體網格時，用四面體單元而不要用退化的四面體單元。比如95號單元有20節點，可以退化為10節點四面體單元，而92號單元為10節點單元，在此情況下用92號單元將優于95號單元。
選擇正確的求解器。對大規模問題，建議采用PCG法。此法比波前法計算速度要快10倍以上（前提是您的計算機內存較大）。
展開 
ansys計算橡膠不收斂的相關專題、標簽、搜索
ansys計算橡膠不收斂的最新內容
fluent模擬相變材料的問題，初始化的時候出現警告，后續計算一直不收斂
眾所周知，橡膠作為不可壓縮材料，具有很強的材料非線性，在計算過程中，往往會因為產生過大的變形導致網格畸變嚴重，求解終止。不僅如此，橡膠材料計算過程中往往還會帶有復雜的幾何非線性和接觸非線性，使求解收斂難度進一步提高，本篇文章將介紹在Workbench中如何對橡膠計算的收斂進行調試。
本實例研究橡膠的擠壓問題的收斂方法，如圖1，上下實體是結構鋼材料，中間實體為橡膠，固定約束下端實體底面以及中間橡膠體的兩個端面
如題我們在ansys workbench進行仿真計算的時候總會發現整個AWB文件占用了超大的內存，動輒5/6個GB，多的高達200/300GB，這樣給我們為客戶或者拷貝文件帶來了很大難題，如何復制和傳輸這種大型文件十分不利，加上主流聊天軟件的流量限制，我們總得借助某度網盤，而某度又對會員限速，總之問題多多。
在此給大家分享一個便捷傳輸workbench文件的方法。
我們平常使用的workbench
非線性問題是什么？
在日常生活中，經常會遇到結構非線性問題。例如，當用釘書針釘紙張時，金屬釘書釘將永久地彎曲成一個不同的形狀(圖 1a)；在一個木架上放置重物，隨著時間的推移木架將越來越下垂(圖 1b)；汽車或卡車上裝載貨物時，輪胎和下面路面間接觸面將隨貨物重量變化(圖 1c)。如果將上述例子的載荷變形曲線畫出來，我們將發現它們都顯示了結構非線性的基本特征—結構剛度改變。
ansys計算之后出現這個錯誤，這是什么原因，怎么解決
A large negative pivot value ( -1.685395134E+09 ) has been encountered
in the global assembled matrix at the UZ degree of freedom of node
2028351. This may
根據ANSYS的使用者反饋，針對非線性接觸問題上的求解，經常會有客戶出現不收斂的情況，在調試收斂性上花費大量的時間。本文主要針對ANSYS 接觸不收斂問題進行方法上的技巧總結，希望通過本文使大家在ANSYS軟件的使用上有更好的體驗。
ANSYS接觸不收斂的原因有非常多的原因，針對每一種不收斂問題，選擇正確的方法都能使不收斂問題解決變得容易起來。在使用軟件中，ANSYS接觸不收斂原因主要有下面這些原因
粘彈性本構方程是研究聚合物的流動性質，polyflow提供了積分型和微分型本構方程，對于初學者在運用這兩種方程的時候經常會遇到一些收斂的問題(滿足網格質量要求情況下)，下面我們簡單來分析一下這其中的原因，為了更好的說明這種現象，研究從KBKZ積分粘彈性方程來說明。
polyflow中KBKZ粘彈性方程
K-BKZ模型能夠很好地描述粘彈性流體剪切變稀，拉伸黏度，以及彈性方面的第一第二法向應力差
我在計算一個大型結構，地震荷載，BEAM188計算時間太長一個小時可能計算了1秒總共40秒，而且越來越慢，不小心早上還停了電如何能使計算加快？
A：調整優化非線性計算的收斂和速度可以說幾乎是一種藝術，即沒有固定的可循規則。
我的經驗是，你的結構的"非線性"越小，非線性的變化越規則，就越容易收斂。 想象一下如果你是手算這個非線性問題，對你來講較容易的，對ANSYS的相應算法也會容易些。
可以把你的地震時程分析拿出幾點
********************************
* Summary of the simulation *
********************************
The computation failed.
***************************
* Expert tool diagnostics
*SET,RAD_I,1250/2
*SET,RAD_PANI,705
*SET,RAD_O,2450/2
*SET,RAD_PANO,1830
*SET,RAD_B,180/2
*SET,M_BOLT,160
*SET,M_BH1,350+170+50
*SET,M_BH2,350+170+250
*SET,M_NUT,280/2
*SET,H_NUT,170
*SET,RAD_DRILL,2140
