abaqus出現(xiàn)Too many attempts made for this increment怎么解決?
2017-06-14 by:CAE仿真在線 來(lái)源:互聯(lián)網(wǎng)
abaqus會(huì)進(jìn)程碰到Too many attempts made for this increment,如下圖
Threre may be many factors that you should check. Some of them you might want to check predominantly are,1_Material properties and equivalency of units 2_ Mesh size and type 3_Boundary conditions 4_constraints such as rigid body motion 5_Step size and no of increments Also look in Abaqus documentation for Common problems in convergence of solution.
模擬計(jì)算的加載過(guò)程包含單個(gè)或多個(gè)步驟,所以要定義分析步。它一般包含分析過(guò)程選擇,載荷選擇,和輸出要求選擇。而且每個(gè)分析步都可以采用不同的載荷、邊界條件、分析過(guò)程和輸出要求。例如:步驟一:將板材夾于剛性?shī)A具上。步驟二:加載使板材變形。步驟三:確定變形板材的自然頻率。增量步是分析步的一部分。在非線性分析中,一個(gè)分析步中施加的總載荷被分解為許多小的增量,這樣就可以按照非線性求解步驟來(lái)進(jìn)行計(jì)算。當(dāng)提出初始增量的大小后,ABAQUS會(huì)自動(dòng)選擇后繼的增量大小。每個(gè)增量步結(jié)束時(shí),結(jié)構(gòu)處于(近似)平衡狀態(tài),結(jié)果可以寫入輸出數(shù)據(jù)庫(kù)文件、重啟動(dòng)文件、數(shù)據(jù)文件或結(jié)果文件中。選擇某一增量步的計(jì)算結(jié)果寫入輸出數(shù)據(jù)庫(kù)文件的數(shù)據(jù)稱為幀。迭代步是在一增量步中找到平衡解的一種嘗試。如果模型在迭代結(jié)束時(shí)不是處于平衡狀態(tài),ABAQUS將進(jìn)行另一輪迭代。隨著每一次迭代,ABAQUS得到的解將更接近平衡狀態(tài);有時(shí)ABAQUS需要進(jìn)行許多次迭代才能得到一平衡解。當(dāng)平衡解得到以后一個(gè)增量步才完成,即結(jié)果只能在一個(gè)增量步的末尾才能獲得。
step,increment,attempt,iteration,的關(guān)系
1)step 分析步
2)increment 時(shí)間增量步
3)attempt 減小增量步的嘗試,即“cutback”
4)iteration 迭代
在一個(gè)計(jì)算中有可能用到多步分析,比如建一個(gè)土石壩,每激活(add)一個(gè)填筑層就是一個(gè)分析步step;
在每個(gè)step中,如果考慮非線性,step就會(huì)分成幾個(gè)增量步(increment)進(jìn)行計(jì)算;
在每個(gè)increment中,會(huì)有減小增量步的嘗試(attempt),在每個(gè)attemp中,要進(jìn)行迭代計(jì)算(iteration)。
如果迭代收斂,則在下一個(gè)increment中會(huì)增大時(shí)間增量步(比如第一個(gè)increment=0.2,則下一個(gè)會(huì)增大為0.3)
如果迭代無(wú)法達(dá)到收斂,則ABAQUS會(huì)自動(dòng)減小時(shí)間增量步(減小increment),即所謂的“cutback”,如果仍然不能收斂,則會(huì)繼續(xù)減小時(shí)間增量步,默認(rèn)的cutback最大次數(shù)為5次,也就是attempt最大=5,如果5次之后仍不能收斂則ABAQUS會(huì)停止分析,顯示錯(cuò)誤:too many attempts made for this increment:analysis terminated.
increment時(shí)間增量步有最小值,默認(rèn)的是1e-5,如果increment減小到比這還小,ABAQUS就會(huì)停止分析,出現(xiàn)錯(cuò)誤:time increment required is less than the minimum specified.
increment的值可以在關(guān)鍵字*static中修改:
*static 1., 1., 1e-05, 1.
分別為初始增量步,分析時(shí)間步,最小增量步,最大增量步
可以用關(guān)鍵字*Step設(shè)定一個(gè)分析步中increment的最大步數(shù),如:
*Step,INC=600 (the maximum number of increments in a step,默認(rèn)的是100 )
*static和*Step中的increment是相同的,*Step,INC默認(rèn)為100,而*static中默認(rèn)為1e-5,并不是100*(1e-5)=1,這兩個(gè)數(shù)都是限值,即number of increments最大為100,而increment最小為1e-5。
這種問(wèn)題怎么解決呢?
問(wèn)題:怎樣修改這個(gè)ABAQUS默認(rèn)的Cutback最大次數(shù)為5次的限制,因?yàn)槲业膍inmun increment size是1e-07,最后計(jì)算狀態(tài)是:
2 2775 1U 0 9 9 8.67 7.67 0.08653
2 2775 2U 0 5 5 8.67 7.67 0.02163
2 2775 3U 0 6 6 8.67 7.67 0.005408
2 2775 4U 0 6 6 8.67 7.67 0.001352
2 2775 5U 0 4 4 8.67 7.67 0.0003380
THE ANALYSIS HAS NOT BEEN COMPLETED
還遠(yuǎn)沒(méi)有達(dá)到我的minmun increment size=1e-07的限制,如果修改默認(rèn)的Cutback最大次數(shù),可能可以收斂。
** CONTROLS
**
*Controls, reset
*Controls, parameters=time incrementation
, , , , , , , 10, , ,
**
** OUTPUT REQUESTS
..........
上面inp文件中參數(shù)10就是我把原來(lái)缺省的Cutback最大5次限制設(shè)置到了10次
希望對(duì)各位有幫助,我試過(guò)了,可以收斂了。
** BOUNDARYCONDITIONS ** ......... ** CONTROLS ** *Controls, reset *Controls, parameters=time incrementation , , , , , , , 10, , , ** ** OUTPUT REQUESTS ** *Restart, write, frequency=0 ** **FIELDOUTPUT: F-Output-1 ** *Output, field, frequency=3 *Node Output
那要怎么修改這個(gè)值呢? 可以參考下圖,先導(dǎo)出inp,修改inp,在通過(guò)inp創(chuàng)建job,然后再執(zhí)行即可
這樣可以定位了在** BOUNDARY CONDITIONS和*Restart, write, frequency=0 之間 |
一般aba的分析精度到E-5,若是小于這個(gè)數(shù)量級(jí)還是不收斂,那么你幾時(shí)改小了初始時(shí)間增量,能計(jì)算出來(lái),結(jié)果的可靠性還是大打折扣。 個(gè)人認(rèn)為,出現(xiàn)這種問(wèn)題,屬于嚴(yán)重的不收斂問(wèn)題,最好還是從模型上找出原因,肯定是模型上有一些影響收斂的錯(cuò)誤。 還有另一種原因,可能是發(fā)生了計(jì)算過(guò)程中結(jié)果材料損傷嚴(yán)重塑性變形進(jìn)行高度非線性,這樣的話那么abaqus/Standard里面的隱式分析,速度很很受影響,也會(huì)出現(xiàn)收斂困難或者不收斂的問(wèn)題,對(duì)于動(dòng)力分析來(lái)講,可以考慮采用abaqus/Explicit的顯式積分算法,倒是一個(gè)不錯(cuò)的選擇。 |
I am using the macro tool to record myself importing geometry and then setting up contact between parts automatically.
Then I change the 2d iges files I am importing using another cadpackage(catia) and am trying to get a contatc analysis to run with the specific geometry I have imported.
Anyway the problem I am having is that sometimes when I import the new gemetry and submit ajobit works and other times I get the error 'too many attempts made for this increment error'
I am quite confused and hoping someone may know what this means;
The way of solving problems is to check all infos in your .msgfile:
first by running an analysis with *PREPRINT,CONTACT=YES in the model part of your .inp (between *HEADING and the first STEP) and adding *CONTACT PRINT in each step.
your analysis will give you very big .dat and .msg files with full info about all the contacts in the model.
Then many options :
- you read the contact details and detect "chattering" : a second order node opens, closes, opens, etc ... through the increments. Then use a more appropriate element type (linear element or 2nd order modified element - for example C3D10M instead of a C3D10)
- if you see that nearly all your contacts are realized and few are still opened or closed, you can allow code to make more attempts in step by using the *CONTROLS card (be careful this area is full of parameters, check it twice before re-running analysis) the default must be about 5 attempts you can define a bit more
- if you noticed a large occurence of the word "overclosure" followed by tiny values (i-e 1E-7) AND if your loading/boundary condition is not too large, you can switch the automatic tolerances by *CONTACT CONTROLS, AUTOMATIC TOLERANCES in the step part. ((reset to *CONTACT CONTROLS, RESET) in the next step .
To limit the problems with contact always try to use displacement approach instead of FOrce approach if possible.
I am a MSc student from BUET . i am modelling a Rcc beam using ABAQUS.The dimension of the beam is 200mmX300mmX3000mm. i choosed elemnt size 200mmX30mmX50mm . Used the general static step.
for concrete i euse
*elastic
33000, 0.2
*concrete
17.4,0? (i used 17.4=0.45*.85*f'c=fy)
34.5,0.0029 ? (34.5= 0.85*fic)
*tension stiffening
1,0
0,0.193
*FAILURE RATIOS
1.16 ,0.09,1.28,0.333
for steel i am using
*elastic
200000,0.31
*PLASTIC
439.964,0
753.13,0.1785
*embedded element,host elset=all_conc_elements
reinforcement
the problem is , i am getting much stiffer result. the initial portion of the load displacement curve is at a much higher (load almost double) position. i am also not getting the displacement correct. i get result up to 10 mm of displacement but the experimental curve has 80 mm displacement. In the msg file i find the same message everytime Too many attempts made for this increment. also there are negetive eigen values.
If any one can suggest me about this stiffer result, it would be a great help
相關(guān)標(biāo)簽搜索:abaqus出現(xiàn)Too many attempts made for this increment怎么解決? abaqus分析培訓(xùn) abaqus技術(shù)教程 abaqus巖土分析 鋼筋混凝土仿真 abaqus分析理論 abaqus軟件下載 abaqus umat用戶子程序編程 Abaqus代做 Abaqus基礎(chǔ)知識(shí) Fluent、CFX流體分析 HFSS電磁分析 Ansys培訓(xùn)