
發(fā)布
注冊
/
登錄增量步的案例
abaqus增量步的理解
abaqus 對于計(jì)算分解為了step(分析步)-increment(增量步)-迭代步三部分
這其中最需要理解的就是增量步,比如增量步是如何對結(jié)果進(jìn)行影響的,增量步設(shè)置中的innital min和max 三個屬性值如何設(shè)置。
談下自己對于增量步的理解,增量步的增量是針對什么來的呢?很多人以為是時間的增量,我感覺這樣理解是不正確的,因?yàn)?em>增量步針對的是荷載的增量
我們看這個公式 初始增量步中的荷載大小= 總的荷載大小×初始增量步時間/分析步總時間
所以abaqus在進(jìn)行尤其是一些非線性問題的計(jì)算的時候,將step拆分成很小的增量步,就是通過還很小的荷載增量下,得到結(jié)果,然后再進(jìn)行下一步的分析的。
我們看下,abaqus在增量步中是如何計(jì)算的(取自abaqus有限元常見問題解答)
BA首先根據(jù)用戶指定的初始增量步大小進(jìn)行迭代;
缺省狀態(tài)下,若在某一增量步16次迭代中仍不收斂(或者出現(xiàn)發(fā)散),ABA會放棄當(dāng)前增量步,并將增量步大小減小為當(dāng)前增量步大小的0.25倍,重新開始迭代嘗試;
若在減小后的增量步16次迭代中仍不收斂(或者出現(xiàn)發(fā)散),ABA會繼續(xù)放棄當(dāng)前的增量步,并將繼續(xù)減小增量步的大小至當(dāng)前增量步大小的0.25倍,重新開始迭代嘗試;
ABA總是利用比較小的載荷增量來嘗試找到收斂解。若此增量步仍不收斂,ABA將再次減小增量步的大小;
以此類推,直至得到收斂解,或者達(dá)到ABA所允許的一個增量步中最多5次增量步減小的限制,而終止分析;
如果連續(xù)2個增量步只需少于5次的迭代就獲得收斂解,ABA會自動地將增量步的大小增大至當(dāng)前增量步的1.5倍,繼續(xù)進(jìn)行下一個增量步的計(jì)算。
展開 ABAQUS軟件中分析步增量步如何設(shè)置?
在ABAQUS軟件中的分析步(Step)設(shè)置界面中,增量步大小的初始值、最小值、最大值以及最大增量步數(shù)這4 個量之間的關(guān)系怎樣?又應(yīng)如何設(shè)置?
首先,我們需要清楚ABAQUS的計(jì)算迭代過程:ABAQUS軟件首先用增量步的初始值進(jìn)行迭代計(jì)算,如果計(jì)算結(jié)果收斂,則以該值代入下一步計(jì)算,若計(jì)算結(jié)果依然收斂,為了節(jié)約計(jì)算成本,ABAQUS軟件會自動嘗試增加增量步大小進(jìn)行迭代計(jì)算;如果計(jì)算結(jié)果出現(xiàn)不收斂現(xiàn)象(監(jiān)控器屬性欄出現(xiàn)字母“U”),則ABQUS軟件自動減小時間步長重新計(jì)算,直至計(jì)算結(jié)果收斂,然后再將該值代入下一步計(jì)算中,依此往復(fù)迭代。如果時間步長減小到增量步的最小值時計(jì)算結(jié)果仍不收斂,ABAQUS軟件將中止計(jì)算,判定計(jì)算結(jié)果不收斂。
搞清楚迭代原理之后,我們就知道如何設(shè)置這四個量的具體參數(shù)值了。對于容易收斂的問題且對相關(guān)變量的過程變化不做要求的仿真分析,為了節(jié)約計(jì)算成本,增量步初始值一般保持默認(rèn),設(shè)為1即可。但是,對于難于收斂的非線性問題或者我們比較關(guān)心模型加載的過程,增量步初始值可適當(dāng)設(shè)小。需要說明的是增量初始值如果設(shè)置太小,會增加我們的計(jì)算時間,如果設(shè)置過大,ABAQUS被迫進(jìn)行多次“折減”,甚至直接導(dǎo)致計(jì)算不收斂。
增量步的最小值一般使用默認(rèn)值,對于復(fù)雜非線性問題,可酌情再減少1~2個數(shù)量級,如果計(jì)算還不收斂,可考慮減少空間步長(網(wǎng)格尺寸)。
增量步的最大值對收斂沒有影響,一般采用默認(rèn)值(分析步時間)。
最大增量步數(shù)默認(rèn)值為100.對于一些復(fù)雜的問題,可以酌情將此參數(shù)設(shè)置大些。
展開 abaqus增量步的理解
abaqus 對于計(jì)算分解為了step(分析步)-increment(增量步)-迭代步三部分
這其中最需要理解的就是增量步,比如增量步是如何對結(jié)果進(jìn)行影響的,增量步設(shè)置中的
innital min和max 三個屬性值如何設(shè)置。
談下自己對于增量步的理解,增量步的增量是針對什么來的呢?很多人以為是時間的增量,我感覺這樣理解是不正確的,因?yàn)?em>增量步針對的是荷載的增量
我們看這個公式 初始增量步中的荷載大小= 總的荷載大小×初始增量步時間/分析步總時間
所以abaqus在進(jìn)行尤其是一些非線性問題的計(jì)算的時候,將step拆分成很小的增量步,就是通過還很小的荷載增量下,得到結(jié)果,然后再進(jìn)行下一步的分析的。
我們看下,abaqus在增量步中是如何計(jì)算的(取自abaqus有限元常見問題解答)
BA首先根據(jù)用戶指定的初始增量步大小進(jìn)行迭代;
缺省狀態(tài)下,若在某一增量步16次迭代中仍不收斂(或者出現(xiàn)發(fā)散),ABA會放棄當(dāng)前增量步,并將增量步大小減小為當(dāng)前增量步大小的0.25倍,重新開始迭代嘗試;
若在減小后的增量步16次迭代中仍不收斂(或者出現(xiàn)發(fā)散),ABA會繼續(xù)放棄當(dāng)前的增量步,并將繼續(xù)減小增量步的大小至當(dāng)前增量步大小的0.25倍,重新開始迭代嘗試;
ABA總是利用比較小的載荷增量來嘗試找到收斂解。若此增量步仍不收斂,ABA將再次減小增量步的大小;
以此類推,直至得到收斂解,或者達(dá)到ABA所允許的一個增量步中最多5次增量步減小的限制,而終止分析;
如果連續(xù)2個增量步只需少于5次的迭代就獲得收斂解,ABA會自動地將增量步的大小增大至當(dāng)前增量步的1.5倍,繼續(xù)進(jìn)行下一個增量步的計(jì)算。
展開 abaqus增量步的理解
abaqus 對于計(jì)算分解為了step(分析步)-increment(增量步)-迭代步三部分
這其中最需要理解的就是增量步,比如增量步是如何對結(jié)果進(jìn)行影響的,增量步設(shè)置中的
innital min和max 三個屬性值如何設(shè)置。
談下自己對于增量步的理解,增量步的增量是針對什么來的呢?很多人以為是時間的增量,我感覺這樣理解是不正確的,因?yàn)?em>增量步針對的是荷載的增量
我們看這個公式 初始增量步中的荷載大小= 總的荷載大小×初始增量步時間/分析步總時間
所以abaqus在進(jìn)行尤其是一些非線性問題的計(jì)算的時候,將step拆分成很小的增量步,就是通過還很小的荷載增量下,得到結(jié)果,然后再進(jìn)行下一步的分析的。
我們看下,abaqus在增量步中是如何計(jì)算的(取自abaqus有限元常見問題解答)
BA首先根據(jù)用戶指定的初始增量步大小進(jìn)行迭代;
缺省狀態(tài)下,若在某一增量步16次迭代中仍不收斂(或者出現(xiàn)發(fā)散),ABA會放棄當(dāng)前增量步,并將增量步大小減小為當(dāng)前增量步大小的0.25倍,重新開始迭代嘗試;
若在減小后的增量步16次迭代中仍不收斂(或者出現(xiàn)發(fā)散),ABA會繼續(xù)放棄當(dāng)前的增量步,并將繼續(xù)減小增量步的大小至當(dāng)前增量步大小的0.25倍,重新開始迭代嘗試;
ABA總是利用比較小的載荷增量來嘗試找到收斂解。若此增量步仍不收斂,ABA將再次減小增量步的大小;
以此類推,直至得到收斂解,或者達(dá)到ABA所允許的一個增量步中最多5次增量步減小的限制,而終止分析;
如果連續(xù)2個增量步只需少于5次的迭代就獲得收斂解,ABA會自動地將增量步的大小增大至當(dāng)前增量步的1.5倍,繼續(xù)進(jìn)行下一個增量步的計(jì)算。
展開 
abaqus某分析步初始增量步被自動縮減
請教各位,最近abaqus standard使用python腳本設(shè)置某分析步初始增量步0.025,inp文件也是記錄0.025,但是實(shí)際腳本命令mdb.jobs[].submit計(jì)算發(fā)現(xiàn)初始增量步取了0.008,改用bat文件也是被自動“修正”為0.008。更改初始增量步或把inp文件放到別的電腦上也會出現(xiàn)這樣的縮減。不明白為什么會出現(xiàn)這種情況
最近一批相同構(gòu)造不同尺寸的模型在計(jì)算,只有一個模型出現(xiàn)這種情況
以下是inp文件中某一分析步(第三個分析步)設(shè)定的增量步參數(shù)
......
** ----------------------------------------------------------------
**
** step: ml
**
*step, name=ml, nlgeom=yes, inc=10000
*static
0.025, 1., 1e-08, 0.025
**
......
展開 Abaqus結(jié)果云圖中無增量步顯示(ODB加載問題)
<p>我做的是掃頻分析</p><p>掃頻分析步<strong>已添加場輸出(默認(rèn)是沒有場輸出的)</strong>,提交計(jì)算完成后,監(jiān)視器中有增量步,結(jié)果云圖中無增量步顯示(ODB加載問題)解決方法:</p><p><br></p><ol><li>上方菜單欄File-Close ODB-全部清除</li><li>Job模塊-重新打開計(jì)算文件的Results</li><li>查看云圖,增量步結(jié)果可查看</li></ol><p><br></p><p>總結(jié):關(guān)閉ODB,重新打開</p><p><br></p>
展開 umat等子程序用open和write輸出數(shù)據(jù),出現(xiàn)一直顯示運(yùn)行中,增量步不計(jì)算的現(xiàn)象。 ¥10
但有時候會一直顯示運(yùn)行中,abaqus沒有增量步計(jì)算。我發(fā)現(xiàn)大量論壇上很多人出現(xiàn)了這個問題,但一直未解決。如何解決這個問題,我嘗試了多種方法,調(diào)試了大量程序,最終發(fā)現(xiàn)可以通過以下方法解決。以下內(nèi)容包括
產(chǎn)生此現(xiàn)象”一直顯示運(yùn)行中,abaqus沒有增量步計(jì)算“的原因分析;
解決方案;
關(guān)于ABAQUS分析步的問答
問題:在建模過程中Step模塊的分析步以及每個分析步下的初始增量步、最大增量步、最小增量步的具體含義,ABAQUS在求解一個非線性問題時是怎樣進(jìn)行迭代的,如何去判斷每個增量步迭代的平衡條件?等等。
答ABAQUS/Standard對于非線性問題的求解采用的是Newton-Raphson算法來實(shí)現(xiàn)。通過對每一個分析步下的增量步進(jìn)行多次迭代,來使每個增量步達(dá)到收斂,進(jìn)而得到該分析步下的收斂解。在迭代的過程中,ABAQUS會根據(jù)收斂情況,自動地對增量步進(jìn)行擴(kuò)大或折減。具體過程如下:如果一個增量步在16次迭代之內(nèi)獲得了收斂解,則成功結(jié)束當(dāng)前的增量步,并開始求解下一個增量步。如果兩個連續(xù)的增量步都在5次迭代之內(nèi)就獲得了收斂解,ABAQUS/Standard自動將下一個增量步增大為當(dāng)前增量步的150%。這個過程叫做增量步的“擴(kuò)大”。如果一個增量步經(jīng)過16次迭代仍沒有獲得收斂解,或者計(jì)算結(jié)果是發(fā)散的, 系統(tǒng)會將增量步減小為當(dāng)前增量步的25%,重新開始迭代嘗試,此過程稱為“折減”。當(dāng)折減次數(shù)超過5次,那么就會出現(xiàn)我們經(jīng)常遇到的錯誤信息:“ERROR:?TOO?MANY?ATTEMPTS?MADE?FOR?THIS?INCREMENT:?ANALYSIS?TERMINATED?”。造成這樣的問題往往是因?yàn)槟P偷谋旧碛袉栴},例如存在剛體位移、過約束、接觸或者塑性材料定義不當(dāng)、網(wǎng)格過于粗糙或過于細(xì)化等。在進(jìn)行非線性分析前,一般都要對最大增量步的數(shù)目、初始增量步、最大和最小增量步進(jìn)行適當(dāng)?shù)脑O(shè)定,來保證求解的順利進(jìn)行。
1.初始增量步:對于很容易收斂的問題,一般設(shè)定為1即可;對于難以收斂的非線性問題,需減小初始增量步,如將分析步時間乘以或(這個需根據(jù)問題的具體情況決定)。
2.?
展開 ABAQUS分析步的一些基本概念
模擬計(jì)算的加載過程包含單個或多個步驟,所以要定義分析步。它一般包含分析過程選擇、載荷選擇、輸出要求選擇。且每個分析步均可采用不同的載荷、邊界條件、分析過程和輸出要求。例如在模擬蠕變成型時,涉及到典型的三個過程,分別是壓縮、壓縮狀態(tài)保持一定的時間、卸載,針對該過程需要定義三個分析步,分別模擬上述三個過程。
增量步是分析步的一部分。在非線性分析中,一個分析步中施加的總載荷被分解為許多小的增量,這樣就可以按照非線性求解步驟來進(jìn)行計(jì)算。當(dāng)提出初始增量的大小后,ABAQUS
會自動選擇后繼的增量大小。每個增量步結(jié)束時,結(jié)構(gòu)處于(近似)平衡狀態(tài),結(jié)果可以寫入輸出數(shù)據(jù)庫文件、重啟動文件、數(shù)據(jù)文件或結(jié)果文件中。選擇某一增量步的計(jì)算結(jié)果寫入輸出數(shù)據(jù)庫文件的數(shù)據(jù)稱為幀。迭代步是在一增量步中找到平衡解的一種嘗試。如果模型在迭代結(jié)束時不是處于平衡狀態(tài),ABAQUS
將進(jìn)行另一輪迭代。隨著每一次迭代,ABAQUS
得到的解將更接近平衡狀態(tài);有時ABAQUS
需要進(jìn)行許多次迭代才能得到一平衡解。當(dāng)平衡解得到以后一個增量步才完成,即結(jié)果只能在一個增量步的末尾才能獲得。
ABAQUS/Standard對于非線性問題的求解采用的是Newton-Raphson算法來實(shí)現(xiàn)。通過對每一個分析步下的增量步進(jìn)行多次迭代,來使每個增量步達(dá)到收斂,進(jìn)而得到該分析步下的收斂解。在迭代的過程中,ABAQUS會根據(jù)收斂情況,自動地對增量步進(jìn)行擴(kuò)大或折減。具體過程如下:
如果一個增量步在16次迭代之內(nèi)獲得了收斂解,則成功結(jié)束當(dāng)前的增量步,并開始求解下一個增量步。如果兩個連續(xù)的增量步都在5次迭代之內(nèi)就獲得了收斂解,ABAQUS/Standard自動將下一個增量步增大為當(dāng)前增量步的150%。這個過程叫做增量步的“擴(kuò)大”。
展開 ABAQUS中的“CLOAD”、“BOUDARY”關(guān)鍵字使用思考
依據(jù)增量迭代法的實(shí)現(xiàn)流程,我猜測ABAQUS的計(jì)算流程是這樣的:abaqus依據(jù)cload給出的荷載數(shù)值,按照增量迭代的計(jì)算方法,在每一個增量步計(jì)算中,將一定量的荷載施加在結(jié)構(gòu)上,當(dāng)增量步滿足收斂判斷時,進(jìn)入下一個增量步,在下一個增量步中,abaqus再次將一定量的荷載施加在結(jié)構(gòu)上,當(dāng)下一個增量步滿足收斂判斷時,再次開始一個新的增量步,依次類推,一直將荷載疊加代cload關(guān)鍵字后面給出的荷載數(shù)值且收斂時,計(jì)算結(jié)束。
但是,事實(shí)上不是這樣的,關(guān)鍵字cload后面的數(shù)值并不是總的荷載數(shù)值,準(zhǔn)確來說,在動力隱式、固定增量步長計(jì)算中,這個荷載數(shù)值是abaqus在每一個增量步計(jì)算中都會施加在結(jié)構(gòu)上的荷載數(shù)值。
()下面是驗(yàn)證部分
設(shè)計(jì)平面板水平拉伸算例,尺寸35x10,彈性模量:2e15,密度2400,泊松比0.25,荷載以及邊界條件愛你如下圖:左端全固定,右端水平拉伸。
展開 ABAQUS中沖擊動力學(xué)問題的求解方法
應(yīng)用這個速度的變化值加上前一個增量步中點(diǎn)的速度來確定當(dāng)前增量步中點(diǎn)的速度:
(2-3)
速度對時間的積分加上在增量步開始時的位移可以確定增量步結(jié)束時的位移:
(2-4)
這樣,在增量步開始時提供了滿足動力學(xué)平衡條件的加速度。得到了加速度,在時間上“顯式地”前推速度和位移。所謂“顯式”是指在增量步結(jié)束時的狀態(tài)僅依賴于該增量步開始時的位移、速度和加速度。這種方法精確地積分常值的加速度。為了使該方法得到精確的結(jié)果,時間增量必須相非常小,這樣在增量步中加速度幾乎為常數(shù)。由于時間增量步必須很小,所以一個典型的分析需要成千上萬個增量步。因?yàn)椴槐赝瑫r求解聯(lián)立方程組,所以每一個增量步的計(jì)算成本很低。大部分計(jì)算成本消耗在單元計(jì)算上,以此確定作用在節(jié)點(diǎn)上的單元內(nèi)力。單元的計(jì)算包括確定單元應(yīng)變和應(yīng)用材料本構(gòu)關(guān)系(單元剛度)確定單元應(yīng)力,從而進(jìn)一步計(jì)算內(nèi)力。
展開 
ABAQUS非線性分析的平衡迭代過程和收斂原則 附ABAQUS非線性有限元分析與實(shí)例下載
每一個增量步計(jì)算結(jié)束得出的結(jié)果只是一個無限接近解析解的結(jié)果,并不能得到精確解。因此要判斷某一個增量步是否完成了計(jì)算,需要判斷其結(jié)果增量是否滿足某個限制,當(dāng)結(jié)果增量滿足限制要求時,該增量步即可收斂。要完成每一個增量步的計(jì)算,有時需要進(jìn)行多個步驟的迭代才能達(dá)到收斂限值,如果一個迭代步的結(jié)果無法達(dá)到收斂限值,則增加迭代步繼續(xù)計(jì)算,直到達(dá)到限值要求。具體過程做如下分析。
1、通過軟件賦予部件材料本構(gòu)關(guān)系,從而得到材料曲線。
圖1 一個增量步中的迭代過程
2、根據(jù)初始增量步大小,得到初始增量荷載P0,根據(jù)材料本構(gòu)關(guān)系可以得到材料初始剛度K0。從而可以得到初始位移u0。
P0=K0u0
3、進(jìn)入第一個增量步,增量步大小發(fā)生變化,此時荷載增加了△P,變?yōu)镻1,利用初始剛度可以得出一個新的位置值u1。
P
1=K
0u
1
4、在根據(jù)材料本構(gòu)關(guān)系,可以得到材料產(chǎn)生的內(nèi)力I1。
I
1=K
1u
1
5、根據(jù)解析解,此時P1應(yīng)該與I1相等。但采用迭代計(jì)算,此時R1=P1-I1。產(chǎn)生了一個差值R1。
6、此時,需要判斷這個結(jié)果是否滿足要求。根據(jù)ABAQUS默認(rèn)判斷準(zhǔn)則,如果R1小于該整個時間段內(nèi)作用在結(jié)構(gòu)上的平均力的0.5%,則此次迭代的結(jié)果滿足要求。
7、當(dāng)判斷了力的限值以后,還需要判斷結(jié)果位移的限值。新的位移u1與初始u0有一個差值,定義為△u1=u1-u0。如果△u1小于位移增量的1%,則接受結(jié)果。
8、如果上述要求不滿足任意一項(xiàng),則在該增量步中進(jìn)行第二次迭代。這里注意,P1是該增量步中荷載的最大值,P0是該增量步中荷載的最小值。因此無論進(jìn)行多少次迭代,產(chǎn)生的內(nèi)力都必須在這二者之間。
展開 Abaqus子程序系列:FRIC(定義接觸表面的摩擦行為)
TAULM(NFDIR)
如果先前的迭代中,LM設(shè)置為1,該值是增量步結(jié)束時約束應(yīng)力(等效剪切應(yīng)力)的當(dāng)前值。否則,將是0。通過與臨界剪應(yīng)力的比較,可以確定此時這個點(diǎn)是否由粘著變成了滑動。
有些量是對于增量步更新的,有些量是對于增量步中的每次迭代更新的。如上面的DGAM和TAULM是每次迭代中更新的,而DSLIP是增量步中更新的。
PRESS
增量步結(jié)束時的接觸壓力
DPRESS
接觸壓力的增量
DDPDDH
在軟接觸的情況下,此時的接觸剛度
SLIP(NFDIR)
在增量步開始時,總的不可恢復(fù)滑動。這個值是從前面的增量步中,DSLIP(NFDIR)的累計(jì)值。
KSTEP
分析步號
KINC
增量步號
TIME(1)
增量步結(jié)束時,分析步的時間值
TIME(2)
增量步結(jié)束時,總時間的值
DTIME
這個增量步的時間尺寸
NOEL
接觸單元的單元編號。如果定義的是接觸表面,傳遞進(jìn)來的是0
CINAME
與摩擦定義關(guān)聯(lián)的用戶定義的表面的名稱,左對齊。對于接觸單元,它是與摩擦定義相關(guān)聯(lián)的界面定義所給定的單元集名稱。如果給界面定義分配了一個可選名稱,那么CINAME將作為這個名稱傳入,左對齊。
SLNAME
從表面名稱。如果使用接觸單元,傳入空格。
NSNAME
主表面名稱。如果使用接觸單元,傳入空格。
NPT
接觸單元的積分點(diǎn)數(shù)量。如果接觸表面被定義,傳入0。
NODE
與這個接觸點(diǎn)相關(guān)的用戶定義的全局從節(jié)點(diǎn)號(或根據(jù)部件實(shí)例的組裝定義的模型的內(nèi)部節(jié)點(diǎn)編號)。如果使用面對面接觸公式,則對應(yīng)于約束的主要從節(jié)點(diǎn)。如果從接觸單元調(diào)用,則作為零傳遞。
展開 有限元理論基礎(chǔ)及Abaqus內(nèi)部實(shí)現(xiàn)方式研究系列4:非線性問題的求解 ¥1
Abaqus的Genera Static分析步采用增量迭代法來求解非線性分析問題,默認(rèn)采用結(jié)合增量法和Newton-Raphson法的混合法。
下面,筆者將參照這sta文件和msg文件來講此算例中Abaqus非線性問題的求解步驟。
2.2.2 第一增量步
2.2.2.1 第一迭代步
2.2.2.2 第二迭代步
1.2.3 第二增量步
1.2.3.1 第一迭代步
2.2.4 第三增量步
迭代求解過程同上。
2.2.5 最終結(jié)果
最終應(yīng)力和應(yīng)變的計(jì)算結(jié)果如下圖所示。
最終計(jì)算的增量和每個增量步的迭代次數(shù)如下圖所示。
由上圖可知,Abaqus在分別在進(jìn)行第一次計(jì)算第四個增量步和第一次計(jì)算第五個增量步時,因?yàn)榍蠼獠皇諗浚?em>增量步大小減為原來的0.25。
2.3 總結(jié)
本文簡單介紹了非線性問題的幾種常用求解方法,并給出了在結(jié)構(gòu)有限元非線性求解時遵循的收斂準(zhǔn)則。同時,給出了一個非線性彈簧靜力分析的例子(詳見附件SpringNonLinear.cae),并通過這個例子描述了Abaqus采用增量迭代法求解非線性問題的過程。
展開 基于XFEM的裂紋擴(kuò)展仿真過程詳解和仿真經(jīng)驗(yàn)交流(二)(包括直接循環(huán)載荷步疲勞裂紋擴(kuò)展分析) ¥20
然后在create interaction中設(shè)置初始載荷步允許裂紋擴(kuò)展
圖5.1 定義裂紋富集域
圖5.1 允許裂紋擴(kuò)展
定義其他接觸條件:一般定義硬接觸就可以了,在contact property > mechanical > normal behavior >Hard Contact,對于受壓縮載荷的情況這里會有所不同,要考慮到裂紋閉合效應(yīng),需要定義其他的接觸準(zhǔn)則。
(6) 定義載荷步:這里要做的是疲勞裂紋擴(kuò)展,在載荷步的定義問題上網(wǎng)上存在一些分歧,有人認(rèn)為裂紋擴(kuò)展是準(zhǔn)靜態(tài)過程,應(yīng)該定義通用靜態(tài)載荷步,然后在載荷模塊使用循環(huán)載荷;還有人認(rèn)為裂紋擴(kuò)展是受交變載荷的疲勞過程,應(yīng)該采用專門的direct cyclic分析步。筆者同意后者的觀點(diǎn),因?yàn)榍懊婺欠N我沒做出來。
圖6.1 direct cyclic載荷步及其參數(shù)設(shè)置
數(shù)據(jù)說明:
basic頁面定義的1為載荷步的總時間,但在隱式求解中它并不是真實(shí)意義上的時間,這個不用改;
incrementation 頁面定義的是增量步的相關(guān)信息,將一個載荷步離散為多個增量步進(jìn)行迭代求解,可以選用自動增量步或者固定增量步,最大總增量步(10000)和增量步大小(0.01)是兩種離散的方法,實(shí)際的增量步為min(總時間除以增量步大小,最大總增量步),根據(jù)你的要求進(jìn)行設(shè)置,看你是想獲得準(zhǔn)確的增量步還是增量步數(shù)。最大迭代次數(shù)(1000),顧名思義是迭代次數(shù)的上限,在解非線性方程組時采用的迭代求解方法,如果第m步迭代不收斂,第m+1步將上一步的增量步減半再次求解,否則乘以1.5再次求解,直到求解總時間達(dá)到1,在monitor可以清楚的看到這一過程。20,25,5是傅里葉級數(shù)的項(xiàng)數(shù),與求解器有關(guān),我也不是很清楚。
展開