ansys斷裂破壞分析的案例
利用ABAQUS進(jìn)行巖土內(nèi)部斷裂破壞的cohesive單元分析研究
問題引出
我們知道，ABAQUS中經(jīng)常用于模擬裂紋擴(kuò)展或斷裂行為的是Cohesive單元，Cohesive單元可理解為一種準(zhǔn)二維單內(nèi)聚單元，可以將它看作被一個(gè)厚度隔開的兩個(gè)面，這兩個(gè)面分別和其他實(shí)體單元連接。這種單元在很多行業(yè)中的結(jié)構(gòu)關(guān)于損傷仿真研究方面經(jīng)常用到。為此，本文通過建立金屬板的膠合模型，用cohesive單元模擬結(jié)構(gòu)的損傷演化。
2. 模型說明
兩塊金屬板用膠結(jié)合在一起，在法向拉力作用作用下將兩塊板分開，分析在對(duì)金屬板加載過程中膠層的應(yīng)力應(yīng)力及失效過程。金屬板層尺寸為10×10×1mm，膠層厚度為0.1mm，有限元模型如圖1所示，左右兩層體單元為金屬，中間一層為厚度為0的Cohesive單元，此次仿真用的單位制系統(tǒng)是mm、N、MPa。
圖1金屬板膠接模型
3有限元分析
利用ABAQUS進(jìn)行本次Cohesive單元損傷演化分析步驟如圖2所示，主要包括：創(chuàng)建部件及劃分網(wǎng)格、創(chuàng)建材料并給部件賦予材料屬性、裝配、創(chuàng)建參考點(diǎn)和剛體約束、創(chuàng)建分析步、設(shè)置輸出變量、創(chuàng)建邊界條件及施加位移載荷、創(chuàng)建分析作業(yè)并提交分析、可視化處理，其中在兩個(gè)仿真環(huán)節(jié)中會(huì)設(shè)計(jì)到cohesive單元的設(shè)置及后處理操作，因此，本文著重分析這兩個(gè)方面的內(nèi)容。在創(chuàng)建材料賦予材料屬性這里，軟件中是通過單擊工具箱中assign section，單擊sets按鈕，在region assignment對(duì)話框中選擇Sect-cohesive，單擊ok完成，具體操作界面見圖3所示，其他部件操作依次類推。此外在最后的可視化處理操作在后續(xù)介紹。
圖2分析步驟
圖3 region section界面設(shè)置
3可視化后處理操作
具體在ABAQUS中的操作見附件cae格式文件，最終的結(jié)果如下各圖所示。
展開 ANSYS斷裂分析實(shí)例
/POST1
*GET,K,CINT,1,CTIP,1,,5,,K1
*STATUS,K
兩個(gè)應(yīng)力強(qiáng)度因子的計(jì)算結(jié)果基本一致，將斷裂韌性除以K，就可以得出安全系數(shù)，判斷裂紋是否擴(kuò)展。
例三：（交互積分法求應(yīng)力強(qiáng)度因子）
（整理自ANSYS的HELP）
例子位置索引:
有限元模型：
FINISH$/CLEAR
!
關(guān)于ANSYS斷裂力學(xué)分析清單
9) 材料構(gòu)型力:材料力主要用于分析材料的缺陷，如位錯(cuò)，空隙，界面和裂紋等。材料力也成為構(gòu)型力，可以考慮夾雜物中的彈性固體（基體材料）。
對(duì)于線性或非線性彈性材料中，材料力矢量與裂紋面相切的分量代表了裂紋尖端的能量釋放率。此外，裂紋擴(kuò)展方向，非均勻性，缺陷和失配網(wǎng)格也可以使用材料力進(jìn)行表征。在彈塑性力學(xué)問題中，材料力矢量與裂紋面相切（平行）分量代表了裂紋擴(kuò)展驅(qū)動(dòng)力（J積分）。材料力的計(jì)算不考慮作用在裂紋表面的載荷。
10) C*積分:對(duì)于高溫蠕變裂紋擴(kuò)展的研究，目前廣泛采用的控制參量之一是穩(wěn)態(tài)蠕變C*積分。
正如各向同性彈性材料中的J積分一樣，C*積分表征了各向同性材料經(jīng)歷蠕變變形第二階段的裂紋特征。C*積分的表達(dá)式如下：
在彈塑性階段，用以描述裂紋尖端區(qū)域應(yīng)力、應(yīng)變場(chǎng)強(qiáng)度的主要是，積分，因此，積分也就成為了彈塑性斷裂的基本準(zhǔn)則。但材料蠕變條件下，J積分不再適用，此時(shí)能有效地反映裂紋尖端的應(yīng)力應(yīng)變場(chǎng)的是蠕變斷裂參量C*。
來源：本文來自CAE技術(shù)聯(lián)盟公眾號(hào)，版權(quán)歸作者所有。
展開 ANSYS 斷裂力學(xué)新功能之SMART自適應(yīng)裂紋萌生分析
分別定義裂紋生成計(jì)算編號(hào)，斷裂參數(shù)，計(jì)算積分圍線，裂紋萌生的ADPCI編號(hào)，裂紋表面組件名稱。裂紋表面組件名稱（crksurf1和crksurf1）自定義，在裂紋萌生和局部網(wǎng)格發(fā)生變化時(shí)，程序會(huì)自動(dòng)填充節(jié)點(diǎn)列表。如果不明確給出裂紋表面的節(jié)點(diǎn)組件名稱，程序會(huì)自動(dòng)生成兩個(gè)內(nèi)部節(jié)點(diǎn)組件。
CINT,NEW,11
CINT,TYPE,SIFS
CINT,NCON,4
CINT,INIT,1
CINT,SURF,CRKSURF1,CRKSURF2
!! 分別定義裂紋擴(kuò)展分析編號(hào)，對(duì)應(yīng)的擴(kuò)展裂紋編號(hào)，以及使用smart分析方法進(jìn)行裂紋擴(kuò)展分析
CGROW,NEW,31
CGROW,CID,11
CGROW,METHO,SMART,REME
!! 該命令為非必須插入項(xiàng)，由于使用smart分析方法會(huì)自動(dòng)在裂紋處加密網(wǎng)格，為減小計(jì)算量，設(shè)置裂紋擴(kuò)展網(wǎng)格粗化選項(xiàng)，可能會(huì)影響計(jì)算精度
CGROW,RMCONT,coarse,aggr
!! 其中CONS – 使用保守的網(wǎng)格粗化 (default)
MODE – 適中網(wǎng)格粗化.
AGGR – 激進(jìn)的網(wǎng)格粗化策略
5、設(shè)定載荷分析子步，建議設(shè)定較多的初始載荷子步和最小載荷子步捕捉裂紋擴(kuò)展過程，本例設(shè)定了40個(gè)載荷步。
6、提交計(jì)算，計(jì)算過程中，在求解信息中會(huì)出現(xiàn)如圖中是否達(dá)到裂紋插入準(zhǔn)則信息提示。達(dá)到準(zhǔn)則之后程序會(huì)自動(dòng)插入橢圓形裂紋和計(jì)算插入的橢圓形裂紋坐標(biāo)位置和長(zhǎng)短軸長(zhǎng)度，以及輸出使用smart方法計(jì)算時(shí)重劃分網(wǎng)格的數(shù)目信息。
圖4求解過程信息提示
程序在確認(rèn)橢圓形裂紋坐標(biāo)位置時(shí)，會(huì)在每個(gè)子步標(biāo)記所有滿足裂紋萌生條件的節(jié)點(diǎn)，并將它們分組到節(jié)點(diǎn)云中。節(jié)點(diǎn)云的幾何中心是橢圓的中心。
展開 
