利用ANSYS進(jìn)行三維貫穿裂紋的斷裂參數(shù)計(jì)算【轉(zhuǎn)載】

2016-11-16  by:CAE仿真在線  來(lái)源:互聯(lián)網(wǎng)

據(jù)一些工業(yè)化國(guó)家統(tǒng)計(jì),因材料和結(jié)構(gòu)的破壞所造成的損失占國(guó)民經(jīng)濟(jì)生產(chǎn)總值的8%-12%多。破壞事故所造成的人員傷亡的損失更不可估量。我國(guó)作為 一個(gè)發(fā)展中國(guó)家,在這方面的情況比西方發(fā)達(dá)國(guó)家更嚴(yán)重。因此無(wú)論是為了減少破壞事故的損失還是研發(fā)滿足現(xiàn)代工業(yè)所需要的新材料,都要求對(duì)材料的破斷過(guò)程有 科學(xué)的、全面的、定量化的認(rèn)識(shí)。

  三維裂紋作為工程中常見(jiàn)的裂紋形式,早在六十年代初就有不少研究者開(kāi)始研究,到現(xiàn)在已有大量的文獻(xiàn)資料論及這一問(wèn)題,出現(xiàn)了一些有特點(diǎn)的分析方 法。工程上常見(jiàn)的表面裂紋的斷裂分析,由于其實(shí)質(zhì)是三維問(wèn)題,也幾乎同時(shí)開(kāi)始被人們所關(guān)注。三維裂紋問(wèn)題的危害極大,斷裂造成了大量的災(zāi)難性事故發(fā)生,這 使得斷裂力學(xué)在機(jī)械工程、海洋工程、核工程,特別是今天的航空航天工程中受到更廣泛的重視和深入研究。

  因此對(duì)含三維裂紋結(jié)構(gòu)斷裂特性尤其對(duì)三維裂紋體的應(yīng)力強(qiáng)度因子的研究有重要的現(xiàn)實(shí)意義。本文使用ANSYS成功的計(jì)算了三維貫穿裂紋的應(yīng)力強(qiáng)度因子,為計(jì)算三維裂紋提供了一種便捷方式。

 1.模型的建立


利用ANSYS進(jìn)行三維貫穿裂紋的斷裂參數(shù)計(jì)算【轉(zhuǎn)載】 - yliuchengzhuy - yliuchengzhuy的博客

  圖1 三維貫穿裂紋模型

  本文三維裂紋模型長(zhǎng)度為L(zhǎng),高度為H,寬度為W,裂紋半長(zhǎng)為a,裂紋位于模型的中心部位。幾何參數(shù)見(jiàn)表1。模型的為線彈性材料,其彈性模量為2.1E11Pa,泊松比為0.3。模型的邊界條件為:底端固定,頂端承受拉應(yīng)力σ為2E6Pa。

  表1 模型的幾何參數(shù)

  

利用ANSYS進(jìn)行三維貫穿裂紋的斷裂參數(shù)計(jì)算【轉(zhuǎn)載】 - yliuchengzhuy - yliuchengzhuy的博客

  本文采用二維奇異單元PLANE183建立二維的裂紋模型,然后通過(guò)拉伸并使用三維奇異單元SOLID186來(lái)建立三維貫穿裂紋模型。圖2-圖5給出了二維裂紋模型和三維裂紋模型。

  在13.0中對(duì)應(yīng)力強(qiáng)度因子的計(jì)算增加了一種計(jì)算方法即互動(dòng)積分法(Interaction Integrals ),這種方法與計(jì)算J積分的主域積分法類似。在二維問(wèn)題進(jìn)行面積分,在三維問(wèn)題中進(jìn)行體積分來(lái)獲得應(yīng)力強(qiáng)度因子。這種方法與傳統(tǒng)的位移擴(kuò)展法相比精度高, 需要的單元數(shù)少。


利用ANSYS進(jìn)行三維貫穿裂紋的斷裂參數(shù)計(jì)算【轉(zhuǎn)載】 - yliuchengzhuy - yliuchengzhuy的博客

  圖2 二維裂紋模型 圖3 二維裂紋模型裂紋尖端網(wǎng)格


利用ANSYS進(jìn)行三維貫穿裂紋的斷裂參數(shù)計(jì)算【轉(zhuǎn)載】 - yliuchengzhuy - yliuchengzhuy的博客

  圖4 三維貫穿裂紋有限元模型 圖5 三維貫穿裂紋局部圖

 2.結(jié)果與討論

  圖6到圖7給出了三維貫穿裂紋的應(yīng)力等效云圖,通過(guò)圖7可知,SOLID186單元可以很好的模擬出裂紋尖端的應(yīng)力場(chǎng)。

  對(duì)于應(yīng)力強(qiáng)度因子,本文計(jì)算圖5中的三個(gè)位置即位置1(0.159,0.25 ,0.00),位置2(0.159,0.25 ,0.04)和位置3(0.159,0.25 ,0.1),并取8條積分線,計(jì)算結(jié)果列于表2中,通過(guò)表2可知,位置1和位置3的各條應(yīng)力強(qiáng)度因子積分都相同,而位置2只有后三條的應(yīng)力強(qiáng)度因子積分值 與位置1和位置3相同。 出現(xiàn)這種差異的原因,圖8給出了解釋,主要是2位置的等效應(yīng)力與位置1和位置3存在差異,最終造成斷裂參數(shù)的計(jì)算差異。這也說(shuō)明了對(duì)于該模型,不能采用對(duì) 稱的計(jì)算。由表2還可以推出,位置2的I型應(yīng)力強(qiáng)度因子要大于位置1和位置3,當(dāng)外載開(kāi)始增加是,裂紋將從模型的內(nèi)部開(kāi)始擴(kuò)展,這對(duì)于實(shí)際工程有這指導(dǎo)意 義。


利用ANSYS進(jìn)行三維貫穿裂紋的斷裂參數(shù)計(jì)算【轉(zhuǎn)載】 - yliuchengzhuy - yliuchengzhuy的博客

  圖6 三維貫穿裂紋的等效應(yīng)力云圖 圖7 三維貫穿裂紋的局部放大圖

  表2 應(yīng)力強(qiáng)度因子

  

利用ANSYS進(jìn)行三維貫穿裂紋的斷裂參數(shù)計(jì)算【轉(zhuǎn)載】 - yliuchengzhuy - yliuchengzhuy的博客


利用ANSYS進(jìn)行三維貫穿裂紋的斷裂參數(shù)計(jì)算【轉(zhuǎn)載】 - yliuchengzhuy - yliuchengzhuy的博客

  圖8 裂紋前緣的等效應(yīng)力云圖與位置關(guān)系圖

 3.結(jié)論

  通過(guò)以上計(jì)算可以得出一下結(jié)論:

  (1)互動(dòng)積分法可以方便的計(jì)算三維裂紋問(wèn)題;

  (2)對(duì)于本文模型,內(nèi)部裂紋的I型應(yīng)力強(qiáng)度因子大于外表面的I型應(yīng)力強(qiáng)度因子,這對(duì)于實(shí)際工程有這指導(dǎo)意義。


開(kāi)放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:利用ANSYS進(jìn)行三維貫穿裂紋的斷裂參數(shù)計(jì)算【轉(zhuǎn)載】 Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識(shí) ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓(xùn) 

編輯
在線報(bào)名:
  • 客服在線請(qǐng)直接聯(lián)系我們的客服,您也可以通過(guò)下面的方式進(jìn)行在線報(bào)名,我們會(huì)及時(shí)給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國(guó)服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號(hào)金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項(xiàng)目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com