陳躍欣 任中俊

摘 要：該文利用有限元模擬有微裂紋的均質(zhì)彈性材料，在承受拉應(yīng)力時(shí)，微裂紋的裂紋面法向張開(kāi)位移，分析指定裂紋在植入不同數(shù)目裂紋情況下的裂紋面法向張開(kāi)位移變化，分別計(jì)算植入五種不同的裂紋數(shù)目時(shí)，指定裂紋所引起的平均應(yīng)變，并將其中4種裂紋數(shù)目下的平均應(yīng)變與裂紋數(shù)目為2時(shí)的平均應(yīng)變進(jìn)行了比較。計(jì)算結(jié)果表明，在僅有兩條微裂紋情況下，裂紋面法向的張開(kāi)位移隨離裂紋面中線距離呈現(xiàn)對(duì)稱的拱形變化，且變化平滑；裂紋數(shù)量逐漸增多，裂紋間發(fā)生相互作用，使得裂紋面法向的張開(kāi)位移增大，張開(kāi)位移曲線也呈拱形但是變化不平滑。由指定裂紋引起的平均應(yīng)變差值隨裂紋數(shù)量的增加而變大，且裂紋面法向與應(yīng)力成一定夾角的比裂紋面法向與應(yīng)力方向平行的要大。

關(guān)鍵詞：微裂紋相互作用 有限元分析 損傷 有限元建模

中圖分類號(hào)：TU452 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1674-098X（2015）10（a）-0211-02

材料的微觀缺陷是引起材料的損傷原因之一，而微裂紋又是微觀缺陷的最常見(jiàn)的形式。微觀缺陷存在會(huì)導(dǎo)致材料的力學(xué)性能的弱化，如強(qiáng)度、剛度、硬度、穩(wěn)定性以及壽命的改變，并最終引起材料的破壞。研究材料中微裂紋之間的相互作用對(duì)材料的性能影響具有十分重要的意義。

馮西橋[1]指出材料在復(fù)雜加載情況下的微裂紋損傷演化是一個(gè)復(fù)雜問(wèn)題，尤其在考慮微裂紋相互作用時(shí)是很困難的。Budiansky[2]將自洽方法應(yīng)用于計(jì)算微裂紋體的有效彈性模量，并采用能量等效原理，認(rèn)為微裂紋體的應(yīng)變能由無(wú)裂紋基體的應(yīng)變能與裂紋能組成，分析了含隨機(jī)分布橢圓形微裂紋的彈性固體的有效彈性模量。Huang[3]將廣義自洽方法用于描述微裂紋體。Mori-Tanaka[4]將單個(gè)微裂紋置于無(wú)損的彈性基體中，使微裂紋受到的載荷不再是遠(yuǎn)場(chǎng)應(yīng)力，是考慮了微裂紋損傷的有效應(yīng)力方法。Benvensite[4]將Mori-Tanaka方法用于描述微裂紋體，假設(shè)微裂紋置于彈性基體中且承受有效的應(yīng)力，總的應(yīng)變由基體的平均應(yīng)變與微裂紋引起的平均應(yīng)變組成，計(jì)算了微裂紋體的有效彈性模量。Hashin[5]將微裂紋體的有效彈性模量用基體的彈性模量與微裂紋密度參數(shù)，下一狀態(tài)的有效彈性模量則由當(dāng)前狀態(tài)的有效彈性模量與新增的微裂紋密度決定再利用初始條件通過(guò)微分方法計(jì)算含隨機(jī)分布橢圓形微裂紋的彈性固體的有效彈性模量。閆相橋[6]提出了平面彈性介質(zhì)中主裂紋與微裂紋相互作用問(wèn)題的有效數(shù)值計(jì)算方法。夏曉舟[7]等對(duì)裂紋間作用機(jī)制及微裂紋區(qū)對(duì)主裂紋的作用效應(yīng)進(jìn)行了探討。袁小平等[8]研究了特定巖石材料考慮微裂紋相互作用的細(xì)觀力學(xué)彈塑性損傷模型。楊衛(wèi)[9]研究了隨機(jī)分布共線微裂紋的強(qiáng)相互作用。

該文通過(guò)在均質(zhì)彈性正方形薄板材料中植入五種不同數(shù)目的微裂紋，運(yùn)用有限元分別模擬五種不同數(shù)目的微裂紋情況下承受一定拉應(yīng)力時(shí)薄板中指定裂紋的裂紋面法向張開(kāi)位移變化。第一種情況在薄板中植入一條方向與水平方向平行和一條與水平方向夾角為45°的位置相隔一定距離的兩條微裂紋，并確定這兩條裂紋為分析對(duì)象。其他情況的裂紋植入都以前一種情況為基準(zhǔn)，在原有裂紋數(shù)的基礎(chǔ)上增加相同數(shù)目的裂紋，裂紋的方向和位置隨機(jī)布置布置在薄板中，分析對(duì)象同樣為第一情況的確定的兩條裂紋。在垂直于水平方向的薄板兩邊施加均布拉應(yīng)力。輸出五種不同數(shù)目的微裂紋情況下，指定裂紋的裂紋面法向張開(kāi)位移，計(jì)算指定微裂紋引起的平均應(yīng)變。通過(guò)比較不同數(shù)目下原裂紋張開(kāi)位移和引起的平均應(yīng)變來(lái)分析微裂紋之間的相互影響。

1 有限元建模

有限元幾何模型為一正方形薄板，板的邊長(zhǎng)為50 mm。在薄板上分別植入微裂紋條數(shù)為2、4、8、16和32，微裂紋長(zhǎng)度取為2 mm，選取2條不同方向的微裂紋a、b來(lái)分析，其余的微裂紋隨機(jī)分布在微裂紋a、b周圍。其中裂紋條數(shù)為32的情況如圖1所示，圖中紅色的為裂紋a，方向與水平方向平行；圖中黃色的為裂紋b，方向與水平方向成45°夾角。該模型為各向同性彈性模型，均質(zhì)薄板的彈性模量取為40 GPa、泊松比取為0.25。荷載和邊界條件，在方形上下兩端均勻施加100 MPa的拉應(yīng)力，為滿足有限元積分的收斂性要求，在將模型中心一個(gè)節(jié)點(diǎn)在水平和鉛垂兩個(gè)方向的位移進(jìn)行約束。將幾何模型一共劃分了86163個(gè)四邊形單元。

2 有限元分析

薄板受到均勻拉應(yīng)力時(shí)裂紋將產(chǎn)生張開(kāi)位移，而一條裂紋中張開(kāi)位移在裂紋中線兩邊呈對(duì)稱性，取裂紋的一半來(lái)分析。將裂紋a張開(kāi)位移與裂紋位置變化分別在裂紋條數(shù)為2、4、8、16、32的情況下輸出并繪圖如圖2所示。同樣裂紋b張開(kāi)位移與裂紋位置的變化如圖3。從圖2和圖3可以看出，薄板中只有裂紋a、b的情況下，裂紋a、b的裂紋面法向張開(kāi)位移呈現(xiàn)拱形對(duì)稱變化，曲線變化比較平滑；裂紋a張開(kāi)位移比裂紋的b稍大。從圖2、圖3中還可以看到裂紋數(shù)目為4、8、16、32下裂紋a、b的張開(kāi)位移隨裂紋位置曲線都在裂紋數(shù)目為2的曲線之上，且曲線為不平滑的拱形曲線，說(shuō)明隨著裂紋數(shù)目的增加，新增的裂紋對(duì)原有裂紋面法向張開(kāi)位移會(huì)產(chǎn)生相互作用的影響，裂紋數(shù)目的增加增大了裂紋面的法向張開(kāi)位移。圖2中裂紋a的張開(kāi)位移在裂紋條數(shù)為32比裂紋條數(shù)為4、8、16 h；而圖3中裂紋b張開(kāi)位移隨裂紋數(shù)的增加是增大的。

運(yùn)用公式[10]（1）計(jì)算單個(gè)裂紋引起的平均應(yīng)變：

（1）

式中ui、uj為裂紋在水平和鉛垂兩個(gè)方向的張開(kāi)位移；nj、ni為裂紋面法線的方向余弦；V為正方形板的面積；da為裂紋長(zhǎng)度的微分。在植入不同裂紋數(shù)目的情況下，裂紋a的鉛垂方向的張開(kāi)位移如圖2所示，裂紋b的鉛垂方向的張開(kāi)位移如圖3所示；裂紋a的裂紋面法線的方向余弦為（0，1），裂紋b的裂紋面法線的方向余弦為（-0.707，0，707）。該文分別計(jì)算了植入5種不同的裂紋數(shù)目時(shí)，裂紋a和裂紋b所引起的平均應(yīng)變，并將其中4種裂紋數(shù)目下的平均應(yīng)變與裂紋數(shù)目為2時(shí)的平均應(yīng)變進(jìn)行了比較，其計(jì)算結(jié)果如表1、表2所示。

從表1可以看出，由裂紋a引起平均應(yīng)變隨周圍分布的裂紋的數(shù)量增加而增大，且隨著裂紋數(shù)遞增的而遞增，當(dāng)裂紋條數(shù)為4時(shí)和裂紋條數(shù)為2時(shí)的平均應(yīng)變差值已達(dá)到了11.22%；當(dāng)裂紋數(shù)為32時(shí)裂紋a平均應(yīng)變差值已達(dá)到了18.48%。同樣從表2也可以看出，由裂紋b引起的平均應(yīng)變隨周圍分布的裂紋的數(shù)量增加而增大的，但未呈現(xiàn)遞增關(guān)系，裂紋條數(shù)為8比裂紋條數(shù)為4的裂紋b引起平均應(yīng)變差值要小1.38%。從表1、表2能看出裂紋數(shù)目的增加對(duì)裂紋b引起的平均應(yīng)變的影響比對(duì)裂紋a引起的平均應(yīng)變的影響要大，當(dāng)裂紋數(shù)為4時(shí)裂紋a的平均應(yīng)變差值為11.22%而裂紋b的平均應(yīng)變差值為13.79，當(dāng)裂紋數(shù)為32時(shí)，裂紋a的平均應(yīng)變差值為18.48%，裂紋b的平均應(yīng)變差值為20.69%，平均應(yīng)變差值裂紋b比裂紋a已經(jīng)大了2.21%。

3 結(jié)語(yǔ)

該文基于有限元軟件，模擬植入不同數(shù)量的微裂紋均質(zhì)彈性薄板材料施加一定拉應(yīng)力作用下，指定裂紋的裂紋面法向張開(kāi)位移。通過(guò)模擬計(jì)算發(fā)現(xiàn)裂紋數(shù)目的變化會(huì)影響裂紋面法向張開(kāi)位移變化。在僅有2條微裂紋的情況下，裂紋面法向的張開(kāi)位移隨離裂紋中線距離呈現(xiàn)對(duì)稱的拱形變化，在裂紋面中線張開(kāi)位移最大并向裂紋兩端逐漸變小。裂紋的方向與應(yīng)力方向的夾角也會(huì)影響裂紋面法向張開(kāi)位移的變化，當(dāng)裂紋面法向方向和應(yīng)力方向平行時(shí)裂紋面張開(kāi)位移最大，當(dāng)裂紋面法向方向和應(yīng)力方向有夾角時(shí)裂紋面張開(kāi)位移變小。當(dāng)裂紋數(shù)量逐漸增多的情況時(shí)，裂紋間產(chǎn)生相互作用，從而使得原裂紋面法向張開(kāi)位移增大，裂紋面法向的張開(kāi)位移隨離裂紋中線距離也是呈現(xiàn)的拱形變化，但曲線變得不平滑。裂紋數(shù)量的增加使得原有裂紋引起的平均應(yīng)變?cè)龃螅鸭y面法向方向與應(yīng)力方向成一定夾角的比裂紋面法向方向與應(yīng)力方向平行的影響更大。從這些可以看到材料中微裂紋之間會(huì)產(chǎn)生相互影響，使材料的力學(xué)性質(zhì)發(fā)生了改變，如強(qiáng)度的降低剛度的變小變形等。我們?cè)谑褂貌牧蠒r(shí)必須考慮材料中的裂紋數(shù)量的對(duì)正常使用的影響，減少由微裂紋相互作用引起的故障和事故。

參考文獻(xiàn)

[1] 馮西橋，余壽文.包含微裂紋相互作用效應(yīng)的準(zhǔn)細(xì)觀損傷力學(xué)模型[C]//第十屆全國(guó)疲勞與斷裂學(xué)術(shù)會(huì)議，中國(guó)廣州.2000.

[2] B.Budiansky，R.J.OConnell.Elastic moduli of a cracked solid[J].International Journal of Solids and Structures， 1976，12（1）：81-95.

[3] Y.Huang，K.Hu，A.Chandra.A generalized self-consistent mechanics method for microcracked solids[J].Journal of the Mechanics and Physics of Solids，1994（42）：1273-1291.

[4] Y.Benvensite.On the Mori-Tanakas method in cracked bodies[J].Mechanics Research Communications，1986，13（4）：193-201.

[5] Z.Hashin.The differential scheme and its application to cracked materials[J].Journal of the Mechanics and Physics of Solids，1988，36（6）：719-734.

[6] 閆相橋.主裂紋與微裂紋相互作用的有效算法[J].力學(xué)學(xué)報(bào)，2006（1）：119-125.

[7] 夏曉舟，章青，喬丕忠，等.裂紋間作用機(jī)制探討及微裂紋區(qū)對(duì)主裂紋的作用效應(yīng)研究[J].應(yīng)用數(shù)學(xué)和力學(xué)，2010（1）：61-70.

[8] 袁小平，劉紅巖，王志喬.考慮微裂紋相互作用的巖石細(xì)觀力學(xué)彈塑性損傷模型研究[J].固體力學(xué)學(xué)報(bào)，2012（6）：592-602.

[9] 楊衛(wèi)，張宿林，李騰.隨機(jī)分布共線微裂紋的強(qiáng)相互作用[J].力學(xué)季刊，1997（6）： 189-195.

[10] 任中俊，萬(wàn)玲.復(fù)雜應(yīng)力下脆性巖石材料的微裂紋損傷特性[J].應(yīng)用力學(xué)學(xué)報(bào)，2013（2）：7-13.