半無限大體中次表面裂紋間的相互影響

張言庫， 楊鴻達(dá)， 昝曉東， 刑帥兵， 江曉禹

(西南交通大學(xué) 力學(xué)與工程學(xué)院， 成都 610031)

半無限大體中次表面裂紋間的相互影響

張言庫， 楊鴻達(dá)， 昝曉東， 刑帥兵， 江曉禹

(西南交通大學(xué) 力學(xué)與工程學(xué)院， 成都 610031)

針對(duì)鋼軌次表面的微裂紋問題，采用含裂紋的半無限大體進(jìn)行簡化，研究了半無限大體中裂紋間的相互影響。利用彈性力學(xué)中半平面體在邊界上受法向分布力的應(yīng)力公式求出裂紋面上的應(yīng)力，并采用Kachanov簡化方法，得到了裂紋尖端的應(yīng)力強(qiáng)度因子。分析了裂紋與表面的距離、裂紋尺寸和裂紋之間的距離對(duì)裂紋尖端應(yīng)力強(qiáng)度因子的影響。結(jié)果表明：裂紋與表面的距離越大裂紋間的影響越小，裂紋尺寸越大裂紋間的影響越大，裂紋之間的距離越大裂紋間的影響越小。

應(yīng)力強(qiáng)度因子；半無限大體；裂紋

鋼軌在鐵路運(yùn)輸中起著支撐作用。近年來，鐵路運(yùn)輸向著高速、重載以及高密度使用的方向發(fā)展。隨著中國鐵路運(yùn)行速度的提高，鋼軌表面的破壞越來越嚴(yán)重[1]，這些破壞會(huì)對(duì)列車的行車安全造成危害。鋼軌的破壞很大程度上是因?yàn)殇撥壌伪砻娈a(chǎn)生裂紋，由于車輪對(duì)鋼軌的循環(huán)力作用使得鋼軌次表面產(chǎn)生微裂紋，隨著微裂紋擴(kuò)展到一定程度會(huì)導(dǎo)致鋼軌剝離甚至產(chǎn)生斷裂。研究鋼軌次表面的裂紋能更好地了解微裂紋的規(guī)律，更加清晰地認(rèn)識(shí)鋼軌的疲勞斷裂機(jī)理，對(duì)鐵路鋼軌維護(hù)起到強(qiáng)有力的助推作用。這也是裂紋問題一直得到廣大學(xué)者持續(xù)關(guān)注的根本原因所在。

針對(duì)裂紋問題很多專家學(xué)者進(jìn)行了研究[2-14]。閆相橋[2]利用一種邊界元方法來研究有限長主裂紋與微裂紋的相互作用。夏曉舟等[3]研究不同裂紋布置帶來的屏蔽和增強(qiáng)效應(yīng)。王恒等[4]用擴(kuò)展有限元的方法研究了微裂紋對(duì)于主裂紋的屏蔽和增強(qiáng)作用。Soh A K等[7]研究了多條微裂紋對(duì)主裂紋的影響。高山等[8]研究了界面主裂紋和微裂紋間的干涉。陳小翠等研究了基于擴(kuò)展有限元方法的主次裂紋之間的影響[9]。王慶豐等[10]研究了不同位置裂紋間的相互作用。卿海等[11]研究了多條表面裂紋相互作用的應(yīng)力。徐華等[12]利用有限元法直接求解了裂紋群裂尖應(yīng)力強(qiáng)度因子。高健等[13]研究了一維正方準(zhǔn)晶中半無限裂紋問題的解析解。 Chudnovsky A等[14]研究了微裂紋對(duì)主裂紋的影響。這些研究不同程度上都揭示了裂紋產(chǎn)生的原因，但由于真實(shí)的鋼軌裂紋是在鋼軌次表面上出現(xiàn)的，所以這些研究不能很好地解釋鋼軌次表面裂紋之間的影響。

本文利用彈性力學(xué)半平面體上的應(yīng)力公式和Kachanov方法的核心思想，研究鋼軌次表面裂紋之間的相互影響。采用簡化的鋼軌模型，以鋼軌上的裂紋作為研究對(duì)象，計(jì)算了次表面上裂紋之間的影響，得到了裂紋尖端應(yīng)力強(qiáng)度因子的解析解。

1 理論研究

1.1 研究模型

由于中國鐵路技術(shù)的飛速進(jìn)步，鐵路運(yùn)輸向著高速、重載和高密度方向發(fā)展，這就導(dǎo)致鋼軌的接觸疲勞損傷越來越嚴(yán)重。鋼軌表面、次表面的裂紋對(duì)鋼軌的剝落掉塊會(huì)有很大的影響。研究鋼軌的次表面裂紋對(duì)于保護(hù)鋼軌有很重要的意義。圖1是鋼軌實(shí)物[15],左圖是裂紋向左上擴(kuò)展，右圖是裂紋向右擴(kuò)展,橢圓圈標(biāo)出的是裂紋有分叉的位置。

圖1 含裂紋鋼軌實(shí)物

從圖1中的裂紋尺寸可見：鋼軌上的裂紋相對(duì)于鋼軌尺寸(鋼軌高176 mm)是很小的，因此本文將鋼軌處理為半無限大體，且在距離半無限大體表面很近的位置有裂紋，以模擬鋼軌在次表面出現(xiàn)裂紋的情況，并通過施加局部荷載研究裂紋間相互影響作用。

1.2 問題求解

圖2是平面應(yīng)力狀態(tài)的理論分析模型，所示在半無限大體中在邊界ab上作用有分布力，它在各個(gè)點(diǎn)的集度為q=1 N/mm，在ab上距離坐標(biāo)原點(diǎn)ξ處，取微小長度dξ。將dF=qdξ看作一個(gè)微小的集中力，因此可以得出dF=qdξ在M點(diǎn)引起的應(yīng)力。

(1)

(2)

(3)

由式(1)～(3)可以得出：

(4)

(5)

(6)

根據(jù)式(4)～(6)可以求出裂紋上的應(yīng)力，然后采用Kachanov簡化方法[5-6]計(jì)算裂紋尖端的應(yīng)力強(qiáng)度因子：

(7)

圖2 理論示意圖

2條裂紋情況：2條裂紋長度為2L，距離為h，用KⅠ/KⅠ0表示裂紋間相互作用影響。KⅠ表示2條裂紋相互影響時(shí)A裂紋的應(yīng)力強(qiáng)度因子，KⅠ0表示在半無限大體與表面距離較大時(shí)單獨(dú)一條裂紋時(shí)的應(yīng)力強(qiáng)度因子。此時(shí)

(8)

p∞可以用σx在裂紋面上積分求出，

(9)

(10)

(11)

(12)

令圖3中a=1 mm,b=1 mm,L=0.4 mm,可以計(jì)算2條共線裂紋相互影響的結(jié)果。

圖3 2條裂紋的布置圖

3條裂紋情況：每條裂紋長度為2L，距離為h，用KⅠ/KⅠ0表示裂紋間相互作用影響。KⅠ表示3條裂紋相互影響時(shí)B裂紋的應(yīng)力強(qiáng)度因子，KⅠ0表示在半無限大體與表面距離較大時(shí)單獨(dú)一條裂紋時(shí)的應(yīng)力強(qiáng)度因子。此時(shí)：

(13)

p∞可以用σx在裂紋面上積分求出，

(14)

(15)

(16)

(17)

令a=1 mm,b=1 mm,可以計(jì)算3條共線裂紋相互影響的結(jié)果。

圖4 3條裂紋的布置圖

2 計(jì)算結(jié)果與討論

2.1 有限元模型

2.1.1 有限元模型的建立

為了驗(yàn)證第1部分的理論推導(dǎo)，建立了一個(gè)二維的有限元模型。模型的長度為40 mm，寬度為40 mm，裂紋長度為0.8 mm。材料的彈性模量為210 GPa，材料的泊松比為0.3，因此可以當(dāng)作一個(gè)半無限大體。單元類型是plane183。裂紋尖端用奇異三角形單元。在模型局部-1到1的表面施加q=1 N/mm的均布荷載。裂紋區(qū)域網(wǎng)格見圖5。

圖5 裂紋區(qū)域網(wǎng)格

2.1.2 結(jié)果修正

在理論計(jì)算時(shí)采取了半無限的模型，考慮無限域處理帶來的理論誤差，可以基于有限元的結(jié)果對(duì)模型做適當(dāng)?shù)男拚?圖6)。通過對(duì)數(shù)據(jù)的對(duì)比分析發(fā)現(xiàn)，隨著裂紋深度的不同誤差不同，所以采用了修正函數(shù)修正誤差。

(18)

圖6 誤差修正

2.2 2條共線裂紋間的相互影響

這一部分將分析半無限大體中，在局部載荷作用下，2條平行于表面的共線裂紋之間的相互影響。

圖7給出了裂紋尖端應(yīng)力強(qiáng)度因子隨著與表面距離變化的結(jié)果。模型的參數(shù)選取為：a=1 mm,b=-1 mm,L=0.4 mm。從圖7可以看出：隨著與表面距離的增大，裂紋尖端應(yīng)力強(qiáng)度因子越來越小，這是因?yàn)榕c表面距離越大，表面對(duì)裂紋的影響就越小。還可以看出：在與表面距離保持一定時(shí)，裂紋之間的距離越大，裂紋尖端的應(yīng)力強(qiáng)度因子越大，這是因?yàn)榱鸭y之間的距離越大，裂紋之間的影響就越小。

圖7 KⅠ/KⅠ0隨著x的變化曲線

圖8給出了裂紋尖端應(yīng)力強(qiáng)度因子隨著裂紋之間的距離變化的結(jié)果。模型的參數(shù)選取為：a=1 mm,b=1 mm,L=0.4 mm。從圖8可以看出：隨著裂紋之間距離的增大，裂紋尖端應(yīng)力強(qiáng)度因子越來越小，這是因?yàn)榱鸭y之間距離越大，裂紋之間的相互影響就越小。還可以看出：在裂紋間的距離保持一定時(shí)，裂紋與表面距離越大，裂紋尖端的應(yīng)力強(qiáng)度因子越小，這是因?yàn)榕c表面距離越大，表面對(duì)裂紋的影響就越小。這種趨勢(shì)與圖7的趨勢(shì)是一樣的。

2.3 3條共線裂紋間的影響

這一部分將分析半無限大體中，在局部載荷作用下，3條平行于表面的共線裂紋之間的相互影響。

圖9給出的是在局部荷載作用下半無限大體中3條共線裂紋尖端應(yīng)力強(qiáng)度因子隨著裂紋間距離變化的結(jié)果。模型的參數(shù)選取為:a=1 mm,b=1 mm,L=0.4 mm。從圖9可以看出：隨著裂紋間距離的增大，裂紋尖端應(yīng)力強(qiáng)度因子越來越小，這是因?yàn)榱鸭y之間距離越大，裂紋之間的相互影響就越小。還可以看出：在裂紋間的距離保持一定時(shí)，裂紋與表面距離越大，裂紋尖端的應(yīng)力強(qiáng)度因子越小，這是因?yàn)榕c表面距離越大，表面對(duì)裂紋的影響就越小。這種趨勢(shì)與兩條裂紋時(shí)的趨勢(shì)是一樣的。

圖8 KⅠ/KⅠ0隨著h的變化曲線

圖9 隨著h變化的KⅠ/KⅠ0的曲線

圖10給出的是在局部荷載作用下半無限大體中3條共線平行裂紋隨著L變化的裂紋尖端應(yīng)力強(qiáng)度因子結(jié)果。模型的參數(shù)選取為：a=1 mm,b=1 mm,L=0.4 mm,h=2 mm。從圖10可以看出：隨著裂紋長度的增大，裂紋尖端應(yīng)力強(qiáng)度因子越來越大，這是因?yàn)榱鸭y長度越大，裂紋面上的應(yīng)力就越大，并且裂紋尖端會(huì)越接近。同時(shí)能夠看出，在裂紋長度保持一定時(shí)，裂紋之間的距離越大，裂紋尖端的應(yīng)力強(qiáng)度因子越小。

無論是3條共線裂紋還是2條共線裂紋，裂紋的應(yīng)力強(qiáng)度因子之所以會(huì)隨著裂紋長度、距離改變而變化，本質(zhì)上是裂紋上的應(yīng)力場(chǎng)受到其他裂紋的影響而發(fā)生變化。因此，研究這些影響有利于對(duì)裂紋的擴(kuò)展做出理論上的解釋，進(jìn)而解決工程上的問題。

圖10 隨著L變化的KⅠ/KⅠ0的曲線

3 結(jié)論

1) 裂紋與表面的距離對(duì)于裂紋之間影響很大，距離表面越近的裂紋彼此間影響會(huì)越大，容易造成裂紋在鋼軌中擴(kuò)展，使鋼軌表面產(chǎn)生微小剝離，所以應(yīng)該定期檢查鋼軌損壞情況。對(duì)于距離表面比較近的次表面裂紋，應(yīng)該采取預(yù)防性打磨的形式消除剝離的隱患。

2) 裂紋長度會(huì)對(duì)裂紋間的影響產(chǎn)生較大作用，對(duì)于長度比較大的次表面裂紋區(qū)域應(yīng)重點(diǎn)監(jiān)測(cè)檢查，這種長度的裂紋更加容易擴(kuò)展，會(huì)對(duì)鋼軌造成破壞，可以采用維修性打磨的形式防止危害。

3) 裂紋之間的距離對(duì)裂紋間的相互影響會(huì)有顯著作用，對(duì)于次表面裂紋密集的區(qū)域，因?yàn)榱鸭y之間距離很近，這些區(qū)域的裂紋極易相互影響擴(kuò)展形成較長的裂紋，對(duì)鋼軌造成比較嚴(yán)重的破壞，故應(yīng)根據(jù)實(shí)際情況采取維修性打磨形式或替換的措施保護(hù)鋼軌正常使用。

鋼軌的破壞會(huì)影響列車的行車安全，而鋼軌的主要破壞形式是鋼軌踏面的剝離掉塊，因此，消除鋼軌次表面微裂紋的影響帶來的剝離可以有效地延長鋼軌使用壽命，提高鋼軌使用效率。

[1] 郭火明,王文健,劉騰飛,等.重載鐵路鋼軌損傷行為分析[J].中國機(jī)械工程,2014,25(2):267-272.

[2] 閆相橋.主裂紋與微裂紋相互作用的有效算法[J].力學(xué)學(xué)報(bào),2006,38(1):119-124.

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

[4] 王恒,柳占立,許丹丹，等.微裂紋群與主裂紋相互作用的擴(kuò)展有限元分析[C]//北京力學(xué)會(huì)學(xué)術(shù)年會(huì)論文集.北京：北京力學(xué)會(huì)，2015.

[5] KACHANOV M.Elastic solids with many cracks:A simple method of analysis[J].International Journal of Solids & Structures,1987,23(1):23-43.

[6] KACHANOV M.A simple technique of stress analysis in elastic solids with many cracks[J].International Journal of Fracture,1985,28(1):R11-R19.

[7] SOH A K,YANG C H.Numerical modeling of interactions between a macro-crack and a cluster of micro-defects[J].Engineering Fracture Mechanics,2004,71(2):193-217.

[8] 高山,王曉明.界面主裂紋和微裂紋之間的干涉[J].力學(xué)學(xué)報(bào),1997,29(3):287-295.

[9] 陳小翠,杜成斌,江守燕.基于XFEM的主次裂紋間應(yīng)力強(qiáng)度因子相互作用[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,42(4):327-331.

[10] 王慶豐,黃小平,崔維成.不同位置裂紋間的互相作用及其影響規(guī)律的有限元分析[J].江蘇科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2006,20(1):16-19.

[11] 卿海,楊衛(wèi),呂堅(jiān).多條表面裂紋相互作用的應(yīng)力分析[J].工程力學(xué),2009,26(1):1-6.

[12] 徐華，徐德峰，楊綠峰.裂紋群應(yīng)力強(qiáng)度因子分析的廣義參數(shù)有限元法[J].應(yīng)用數(shù)學(xué)和力學(xué)，2016，37(10)：1039-1049.

[13] 高健，劉官廳.一維正方準(zhǔn)晶中半無限裂紋問題的解析解[J].應(yīng)用數(shù)學(xué)和力學(xué)，2015,36(9):945-955.

[14] CHUDNOVSKY A,DOLGOPOLSKY A,KACHANOV M.Elastic interaction of a crack with microcracks[J].Fracture,1984(4):825-832.

[15] 周劍華，任安超，吉玉，等.U71Mn 鋼軌踏面剝離掉塊缺陷分析[J].中國鐵道科學(xué)，2013，34(2)：1-6.

(責(zé)任編輯劉 舸)

InteractionAmongSubsurfaceCracksinanElasticHalf-Plane

ZHANG Yanku, YANG Hongda, ZAN Xiaodong, XING Shuaibing, JIANG Xiaoyu

(School of Mechanics and Engineering, Southwest Jiaotong University, Chengdu 610031, China)

It is simplified as a half-plane contains micro-cracks for the problem of micro-cracks on the subsurface of the rail. It analyzed the interaction between the micro-cracks in the half-plane. The tractions along the cracks are obtained from the stress formula of a half-plane subjected to a normal distribution force at the boundary in the elastic mechanics. Based on Kachanov’s simplified method, the solution of a half-plane containing the micro-cracks is presented. The stress intensity factor at the crack tip is obtained and the interaction among collinear micro-cracks is analyzed. The results show that the interaction among cracks decreases with the distance among the cracks increasing, and the distance between the micro-crack and the surface increasing and the crack size decreasing.

stress intensity factor; semi-infinite; cracks

2017-06-28

國家自然科學(xué)基金資助項(xiàng)目(11472230)；國家自然科學(xué)基金重點(diǎn)資助項(xiàng)目(U1134202,E050303);四川省青年科技創(chuàng)新團(tuán)隊(duì)資助項(xiàng)目(2013TD0004)

張言庫(1990—)，男，山東聊城人，碩士，主要從事應(yīng)用力學(xué)研究，E-mail:1911778208@qq.com。

張言庫， 楊鴻達(dá)， 昝曉東， 等.半無限大體中次表面裂紋間的相互影響[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2017(10):90-95.

formatZHANG Yanku, YANG Hongda, ZAN Xiaodong, et al.Interaction Among Subsurface Cracks in an Elastic Half-Plane[J].Journal of Chongqing University of Technology(Natural Science)，2017(10):90-95.

10.3969/j.issn.1674-8425(z).2017.10.015

O346.1

A

1674-8425(2017)10-0090-06