楊大鵬，潘海洋，劉邦先，張平萍，楊新華

(1.中鋼集團(tuán)鄭州金屬制品研究院有限公司，河南 鄭州，450001；2.華中科技大學(xué)土木工程與力學(xué)學(xué)院，湖北 武漢，430074；3.鄭州職業(yè)技術(shù)學(xué)院機(jī)械工程系，河南 鄭州，450121；4.鄭州職業(yè)技術(shù)學(xué)院城市軌道交通系，河南 鄭州，450121)

疲勞荷載作用下的三維彈塑性彎曲裂紋尖端張開位移

楊大鵬1,2,3，潘海洋4，劉邦先3，張平萍1，楊新華2

(1.中鋼集團(tuán)鄭州金屬制品研究院有限公司，河南 鄭州，450001；2.華中科技大學(xué)土木工程與力學(xué)學(xué)院，湖北 武漢，430074；3.鄭州職業(yè)技術(shù)學(xué)院機(jī)械工程系，河南 鄭州，450121；4.鄭州職業(yè)技術(shù)學(xué)院城市軌道交通系，河南 鄭州，450121)

綜合考慮疲勞作用應(yīng)力、三維塑性區(qū)域邊界上的交變正應(yīng)力與交變剪應(yīng)力，利用二階攝動(dòng)方法建立了計(jì)算疲勞載荷作用下三維彈塑性彎曲裂紋尖端張開位移的理論模型。用數(shù)值解法進(jìn)行求解，并作圖分析了三維彈塑性彎曲裂紋尖端張開位移的最大值和變化幅值與三維裂紋體幾何尺寸及外載荷之間的變化關(guān)系。結(jié)果表明，隨著裂紋體厚度的增大，三維彈塑性彎曲裂紋尖端張開位移的最大值與變化幅值不斷減??；當(dāng)裂紋體幾何尺寸相同時(shí)，彎曲裂紋尖端張開位移的最大值與變化幅值均隨外載荷的增加而逐漸增大。

三維裂紋體；彎曲裂紋；張開位移；疲勞荷載；二階攝動(dòng)方法

目前關(guān)于彎曲裂紋擴(kuò)展路徑的研究大多局限于二維線彈性斷裂和二維彈塑性斷裂問題[1-7]，而針對(duì)疲勞載荷作用下三維裂紋體彈塑性彎曲裂紋的路徑預(yù)測(cè)、三維裂紋體彈塑性彎曲裂紋尖端塑性區(qū)域交變應(yīng)力場(chǎng)的計(jì)算、張開位移最大值與變化幅值的確定等問題的研究相對(duì)缺乏。在工程實(shí)際中，裂紋體材料的厚度通常是不可忽略的，而且往往比較大。因此，研究三維疲勞彈塑性彎曲裂紋的斷裂特性是非常有必要的。本文將運(yùn)用二階攝動(dòng)方法計(jì)算出疲勞載荷作用下三維彈塑性彎曲裂紋尖端張開位移的最大值與變化幅值，從而更精確地服務(wù)于工程實(shí)際。

1 相關(guān)符號(hào)與說明

αZ、βZ、γZ為三維彎曲裂紋的形狀參數(shù)，kIZ、kIIZ、bIZ、bIIZ、TZ是與二維彎曲裂紋相對(duì)應(yīng)的三維彎曲裂紋斷裂特性參量；HZ為三維彎曲裂紋以及塑性裂紋的直線部分的長(zhǎng)度,RZ為三維彎曲裂紋尖端塑性區(qū)在裂紋直線部分延長(zhǎng)線上的射影長(zhǎng)度，RZ可以根據(jù)文獻(xiàn)[11-13]提供的思路、方法與力學(xué)模型求解，aZ為三維彎曲裂紋于裂紋直線部分延長(zhǎng)線上的投影長(zhǎng)度，令cZ=RZ+aZ，δZ為三維彎曲延伸裂紋尖端張開位移。彎曲裂紋尖端場(chǎng)的奇異性特征嚴(yán)重依賴于材料的本構(gòu)類型，本文中的材料本構(gòu)是理想塑性的。

2 三維疲勞彎曲裂紋尖端張開位移的確定

本文采用經(jīng)過拓展的彎曲裂紋Dugdale模型計(jì)算三維疲勞彎曲裂紋尖端張開位移，這個(gè)拓展模型的合理性已經(jīng)在文獻(xiàn)[13]中得到證明。根據(jù)文獻(xiàn)[11-13]，可得如下關(guān)系式：

exp[-B1(z/B)(r/B)B2(z/B)]

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

(11)

(12)

(13)

(14)

(15)

式(11)～式(15)中：kI、kII為應(yīng)力強(qiáng)度因子；bI、bII、T均為系數(shù)，取決于裂紋彎曲延伸擴(kuò)展之前的邊界值問題的解。

同時(shí)，可以確定疲勞載荷作用下某一時(shí)刻三維彎曲裂紋尖端I型張開位移δ1Z、II型張開位移δ2Z的大小：

(16)

(17)

在疲勞荷載作用下，裂紋尖端塑性區(qū)的大小呈交替性變化，三維彎曲裂紋張開位移的周期變化是由于kIZ、kIIZ、bIZ、bIIZ、TZ、RZ等參量的循環(huán)增減而引起的，于是可用下列簡(jiǎn)化函數(shù)式表達(dá):

δ1Z=δ1Z(αZ,βZ,γZ,aZ,t2,σs,kIZ,

kIIZ,bIZ,bIIZ,TZ,RZ)

(18)

(δ1Z)max=δ1Z[αZ,βZ,γZ,aZ,t2,σs,(kIZ)max,

(kIIZ)max,(bIZ)max,(bIIZ)max,(TZ)max,(RZ)max]

(19)

(δ1Z)min=δ1Z[αZ,βZ,γZ,aZ,t2,σs,(kIZ)min,

(kIIZ)min,(bIZ)min,(bIIZ)min,(TZ)min,(RZ)min]

(20)

δ2Z=δ2Z(αZ,βZ,γZ,aZ,t2,σs,kIZ,

kIIZ,bIZ,bIIZ,TZ,RZ)

(21)

(δ2Z)max=δ2Z[αZ,βZ,γZ,aZ,t2,σs,(kIZ)max,

(kIIZ)max,(bIZ)max,(bIIZ)max,(TZ)max,(RZ)max]

(22)

(δ2Z)min=δ2Z[αZ,βZ,γZ,aZ,t2,σs,(kIZ)min,

(kIIZ)min,(bIZ)min,(bIIZ)min,(TZ)min,(RZ)min]

(23)

(24)

(25)

ΔδZ=(δZ)max-(δZ)min

(26)

3 三維疲勞彎曲裂紋尖端張開位移最大值與變化幅值的計(jì)算與分析

下面以碳鋼為例，根據(jù)上述公式用數(shù)值解法計(jì)算材料三維彈塑性彎曲裂紋尖端張開位移，研究張開位移的最大值和變化幅值與三維裂紋體幾何尺寸及外載荷之間的變化關(guān)系。碳鋼的泊松比ν=0.25，屈服極限σs=1725 MPa, 彈性模量E=2.15×105MPa。

當(dāng)垂直于初始直線裂紋方向的外應(yīng)力σaZ=0.5σs、平行于初始裂紋方向的外應(yīng)力σrZ=0.1σs、σZ=50 mm時(shí)，三維彎曲裂紋尖端張開位移的最大值隨裂紋形參αZ的變化關(guān)系如圖1所示。當(dāng)σaZ=0.33σs、σrZ=0.03σs、aZ=50 mm時(shí)，三維彎曲裂紋尖端張開位移的變化幅值隨裂紋形參αZ的變化關(guān)系如圖2所示。

(a)B=10 mm (b)B=20 mm (c)B=30 mm

(d)B=50 mm (e)B=100 mm

圖1三維彎曲裂紋尖端張開位移的最大值隨裂紋形參αZ的變化(σaZ=0.5σs,σrZ=0.1σs，aZ=50mm)

Fig.1Variationofthemaximumof3DcurvedcracktipopeningdisplacementwiththecrackshapeparameterαZ(σaZ=0.5σs,σrZ=0.1σs，aZ=50mm)

從圖1、圖2可以看出，在同樣的外載荷作用下，相同的三維裂紋直線部分長(zhǎng)度對(duì)應(yīng)的彎曲裂紋尖端張開位移的最大值與變化幅值隨著三維裂紋形狀參數(shù)αZ的不斷增大而減小。另外，在外載荷、三維裂紋直線部分長(zhǎng)度HZ和三維彎曲裂紋形狀參數(shù)αZ相同的條件下，當(dāng)三維裂紋體厚度逐漸增加時(shí)，裂紋尖端張開位移的最大值與變化幅值均不斷減小。

三維彎曲裂紋尖端張開位移的最大值和變化幅值隨外載荷(用垂直于初始直線裂紋方向外應(yīng)力的最大值(σaZ)max及其變化值ΔσaZ與屈服極限σs的比值來表征)的變化分別如圖3和圖4所示。

(a)B=10 mm (b)B=20 mm (c)B=30 mm

(d)B=50 mm (e)B=100 mm

圖2三維彎曲裂紋尖端張開位移的變化幅值隨裂紋形參αZ的變化(σaZ=0.33σs,σrZ=0.03σs，aZ=50mm)

Fig.2Variationoftheamplitudeof3DcurvedcracktipopeningdisplacementwiththecrackshapeparameterαZ(σaZ=0.33σs,σrZ=0.03σs，aZ=50mm)

(a)B=10 mm (b)B=20 mm

圖3三維彎曲裂紋尖端張開位移的最大值隨外載荷的變化(aZ=50mm,αZ=0.5°)

Fig.3Variationofthemaximumof3Dcurvedcracktipopeningdisplacementwiththeexternalload(aZ=50mm,αZ=0.5°)

(a)B=30 mm (b)B=50 mm (c)B=100 mm

圖4三維彎曲裂紋尖端張開位移的變化幅值隨外載荷的變化(aZ=50mm,αZ=0.5° )

Fig.4Variationoftheamplitudeof3Dcurvedcracktipopeningdisplacementwiththeexternalload(aZ=50mm,αZ=0.5° )

從圖3、圖4可以看出，當(dāng)三維裂紋體厚度和裂紋直線部分長(zhǎng)度相同時(shí)，三維彎曲裂紋尖端張開位移的最大值與疲勞載荷最大值正相關(guān)，其變化幅值與疲勞載荷變化幅值正相關(guān)。另外,隨著三維裂紋體厚度的增加，對(duì)于相同的裂紋直線部分長(zhǎng)度和彎曲程度，三維彎曲裂紋尖端張開位移的最大值隨疲勞載荷最大值的提高而不斷增大的速度逐漸降低，其變化幅值隨疲勞載荷變化幅值的提高而不斷增大的速度也逐漸降低。

4 三維疲勞彈塑性彎曲裂紋擴(kuò)展的判據(jù)

根據(jù)文獻(xiàn)[11-13]可以得到：

(1)疲勞載荷作用下，三維彎曲裂紋的擴(kuò)展速率為

(27)

式中：R為應(yīng)力比，即R=(σZ)min/(σZ)max；N為疲勞荷載的循環(huán)次數(shù)；C′、m′均為材料常數(shù)；δCZ為三維彎曲裂紋的斷裂韌性；ΔδthZ為三維彎曲裂紋張開位移變化幅度的門檻值。

(2)三維疲勞彈塑性彎曲裂紋在遠(yuǎn)場(chǎng)交變應(yīng)力作用下穩(wěn)定存在的條件為

ΔδZ<ΔδthZ

(28)

(3)三維疲勞彈塑性彎曲裂紋在遠(yuǎn)場(chǎng)交變應(yīng)力作用下穩(wěn)定擴(kuò)展的條件為

ΔδZ≥ΔδthZ，且(δZ)max<δCZ

(29)

(4)三維疲勞彈塑性彎曲裂紋在遠(yuǎn)場(chǎng)交變應(yīng)力作用下失穩(wěn)擴(kuò)展的條件為

ΔδZ≥ΔδthZ，且(δZ)max≥δCZ

(30)

5 結(jié)語(yǔ)

本文運(yùn)用二階攝動(dòng)方法得出三維彈塑性彎曲裂紋尖端張開位移的最大值與變化幅值計(jì)算公式，作圖分析了三維彈塑性彎曲裂紋尖端張開位移的最大值、變化幅值與三維裂紋體厚度、裂紋直線部分長(zhǎng)度、裂紋形狀參數(shù)以及外載荷之間的關(guān)系，并進(jìn)一步討論了三維疲勞彈塑性彎曲裂紋擴(kuò)展的判據(jù)。由于三維疲勞彈塑性彎曲裂紋在土木建筑、航空航天、艦船潛艇等工程結(jié)構(gòu)中普遍存在，因此本研究具有重要的理論意義和工程實(shí)用價(jià)值。

[1] 黃學(xué)偉，蔡力勛，包陳，等.基于低周疲勞損傷的裂紋擴(kuò)展行為數(shù)值模擬新方法[J].工程力學(xué), 2011，28(10):202-208.

[2] 廖芳芳，王偉，陳以一.往復(fù)荷載下鋼結(jié)構(gòu)節(jié)點(diǎn)的超低周疲勞斷裂預(yù)測(cè)[J].同濟(jì)大學(xué)學(xué)報(bào)：自然科學(xué)版，2014，42(4)：539-547.

[3] 嵇醒.斷裂力學(xué)判據(jù)的評(píng)述[J].力學(xué)學(xué)報(bào), 2016，48(4)：741-753.

[4] 王清遠(yuǎn)，劉永杰.結(jié)構(gòu)金屬材料超高周疲勞破壞行為[J].固體力學(xué)學(xué)報(bào), 2010,31(5):496-503.

[5] 魏國(guó)前，岳旭東，黨章，等.結(jié)合S-N曲線和斷裂力學(xué)的焊接結(jié)構(gòu)疲勞壽命分析[J].焊接學(xué)報(bào),2017，38(2):23-27.

[6] Zhao Minghao, Dang Huayang, Xu Guangtao, et al. Dielectric breakdown model for an electrically semi-permeable penny-shaped crack in three-dimensional piezoelectric media[J]. Acta Mechanica Solida Sinica, 2016,29(5):536-546.

[7] Yang Shengqi, Huang Yanhua, Ranjith P G, et al. Discrete element modeling on the crack evolution behavior of brittle sandstone containing three fissures under uniaxial compression[J].Acta Mechanica Sinica, 2015,31(6):871-889.

[8] 丁遂棟,孫利民.斷裂力學(xué)[M].北京：機(jī)械工業(yè)出版社,1997:148-169.

[9] 楊衛(wèi).宏微觀斷裂力學(xué)[M].北京：國(guó)防工業(yè)出版社,1995：65-70.

[10] 楊大鵬,趙耀,白玲.準(zhǔn)靜載作用下彈塑性微彎裂紋尖端塑性區(qū)[J].應(yīng)用力學(xué)學(xué)報(bào),2010,27(2): 401-405.

[11] 郭萬(wàn)林，于培師.構(gòu)件三維斷裂與疲勞力學(xué)及其在航空工程中的應(yīng)用[J].固體力學(xué)學(xué)報(bào)，2010，31(5):553-571.

[12] 張斌.材料結(jié)構(gòu)宏觀三維斷裂和微觀破壞行為研究[D].南京：南京航空航天大學(xué)，2005.

[13] 楊大鵬.微彎延伸裂紋斷裂特性的研究[D].武漢：華中科技大學(xué)，2006.

Threedimensionalelastic-plasticcurvedcracktipopeningdisplacementunderfatigueload

YangDapeng1,2,3,PanHaiyang4,LiuBangxian3,ZhangPingping1,YangXinhua2

(1.Sinosteel Zhengzhou Research Institute of Steel Wire Products Co., Ltd., Zhengzhou 450001, China; 2.School of Civil Engineering and Mechanics, Huazhong University of Science and Technology, Wuhan 430074, China; 3.Department of Mechanical Engineering, Zhengzhou Technical College, Zhengzhou 450121, China; 4.Urban Rail Transport Department, Zhengzhou Technical College, Zhengzhou 450121, China)

Taking the effects of fatigue stress, cyclic normal and shear stresses on the boundaries of three dimensional plastic area into consideration, this paper has built the theoretical models to calculate three dimensional elastic-plastic curved crack tip opening displacement (CCTOD) under fatigue loads by using second order perturbation method. Numerical solutions were made，and the relation curves between the maximum or amplitude of CCTOD and the crack body geometrical dimensions as well as the external loads were analyzed. The results show that the maximum and amplitude of CCTOD decline with the increase of crack body thickness. When crack bodies have the same size, the maximum and amplitude of CCTOD increase with increasing external loads.

three dimensional crack body; curved crack; opening displacement; fatigue load; second order perturbation method

2017-04-10

國(guó)家自然科學(xué)基金重大研究計(jì)劃項(xiàng)目(91016026)；河南省博士后科研資助項(xiàng)目(博士后編號(hào)：166053).

楊大鵬(1976-)，男，中鋼集團(tuán)鄭州金屬制品研究院有限公司博士后工作站以及華中科技大學(xué)土木工程與力學(xué)學(xué)院力學(xué)系博士后流動(dòng)站研究人員，鄭州職業(yè)技術(shù)學(xué)院副教授.E-mail: ydpzpysh@163.com

10.3969/j.issn.1674-3644.2017.05.014

TB301

A

1674-3644(2017)05-0395-06

[責(zé)任編輯尚晶]