拉伸載荷作用下無限大板中三角孔-裂紋問題

劉寶良，閆相橋

(哈爾濱工業(yè)大學(xué)航天學(xué)院，150001哈爾濱，yanxiangqiao406@163.com)

由于孔邊應(yīng)力集中，孔邊很可能起源裂紋.因此，很多研究者［1］專注于孔邊裂紋問題.對(duì)單軸載荷作用下無限大板中源于圓孔的分支裂紋問題，其典型解答由Bowie［2］和 Newman［3］獲得.對(duì)單軸載荷作用下無限大板中源于橢圓孔的對(duì)稱分支裂紋問題，Newman利用邊界配置法、Nisitani等［4］利用體力法獲得其數(shù)值解答.對(duì)雙軸載荷作用下無限大板中源于圓孔的分支裂紋，Tweed等利用Mellin Transform Technique獲得其解答.對(duì)單軸載荷作用下無限大板中源于三角孔或方孔分支裂紋問題，Murakami［5］利用體力法獲得其數(shù)值解.對(duì)內(nèi)部壓力載荷作用下無限大板中源于圓孔的分支裂紋問題，Newman利用邊界配置法獲得其典型解答.

為此，通過把Bueckner原理普遍化，本文提出了一種求解遠(yuǎn)方載荷作用下無限大板中孔邊裂紋問題的數(shù)值方法，通過把適于單一裂紋的Bueckner原理［6］擴(kuò)充到遠(yuǎn)方載荷作用下無限大板中孔邊裂紋問題，將原問題分解為承受遠(yuǎn)方載荷不含孔邊裂紋的均勻問題，和在遠(yuǎn)方不承受載荷但在裂紋面上和孔的表面上承受面力的問題.于是，以應(yīng)力強(qiáng)度因子作為參量的問題可以通過考慮后者來解決，而利用文獻(xiàn)［7］最近提出的雜交位移不連續(xù)法，這種孔邊裂紋問題是很容易數(shù)值求解的.算例說明本數(shù)值方法對(duì)分析無限大板中孔邊裂紋問題既簡單又非常有效.進(jìn)而利用這種數(shù)值方法研究了拉伸載荷作用下無限大板中三角孔-裂紋問題.本文通過改變孔的幾何參數(shù)，通過與無限大板中心裂紋問題的應(yīng)力強(qiáng)度因子比較，揭示了孔的幾何參數(shù)對(duì)應(yīng)力強(qiáng)度因子的影響.發(fā)現(xiàn)孔對(duì)源于它的裂紋的應(yīng)力強(qiáng)度因子具有屏蔽效應(yīng)(shielding effect)和放大效應(yīng)(amplifying effect).這種屏蔽效應(yīng)和放大效應(yīng)隨著孔的幾何參數(shù)而變化.

1 Bueckner原理普遍化

Bueckner利用疊加原理得到一種重要結(jié)果，稱為Bueckner原理.由外部載荷引起的裂紋的應(yīng)力強(qiáng)度因子，可以通過考慮無外部載荷作用但裂紋面上承受面力的裂紋問題來確定，該面力與無裂紋幾何體在外部載荷作用下在裂紋所在位置處產(chǎn)生的應(yīng)力大小相等方向相反.

在這里，嘗試將適于單一裂紋的Bueckner原理擴(kuò)充到擴(kuò)充到遠(yuǎn)方載荷作用下無限大板中孔邊裂紋問題.為此，將原問題分解為承受遠(yuǎn)方載荷不含孔邊裂紋的均勻問題，和在遠(yuǎn)方不承受載荷但在裂紋面上和孔的表面上承受面力的問題.在裂紋面上承受的面力與承受遠(yuǎn)方載荷不含孔邊裂紋的均勻問題在該裂紋所在位置處產(chǎn)生的應(yīng)力大小相等方向相反，即

式中:φ為裂紋面相對(duì)于x軸的方位角;γ為把孔的表面上任意一點(diǎn)相對(duì)于x軸的方位角，則作用于改點(diǎn)處的面力為

當(dāng)然，由式(2)定義的面力隨孔的表面點(diǎn)不同而變化.

于是，以應(yīng)力強(qiáng)度因子作為參量的問題可以通過考慮后者來解決，而利用Yan［7］最近提出的雜交位移不連續(xù)法，這種孔邊裂紋問題是很容易數(shù)值求解的.

2 雜交位移不連續(xù)法

2.1 常位移不連續(xù)單元

無限大平面體中在位置|x|＜a，y=0處具有常位移不連續(xù)量Di=(Dx，Dy)的定義為

式中:|x|＜a，y=0.

此問題的解答是由Crouch等［8］獲得的，位移場和應(yīng)力場為

式中G和v分別為材料剪切模量和泊松比，Crouch等利用式(4)，(5)建立了常位移不連續(xù)邊界元法.

2.2 裂尖位移不連續(xù)單元

基于式(4)～(5)，提出了裂尖位移不連續(xù)單元(可分為左、右裂尖位移不連續(xù)單元)以模擬裂尖附近的應(yīng)力奇異場.

對(duì)于左裂尖位移不連續(xù)單元，其位移不連續(xù)函數(shù)可取為

式中:Hs和Hn分別為裂尖單元中點(diǎn)處的切相和法向位移不連續(xù)量;atip為裂尖單元長度的1/2.在這里注意到裂尖位移不連續(xù)單元與常位移不連續(xù)單元具有相同的未知量，即2個(gè)，但由式(6)定義的位移不連續(xù)函數(shù)可以模擬裂尖附近的位移場，從而可以模擬裂尖附近的應(yīng)力奇異性.

基于式(4)～(5)，根據(jù)微積分學(xué)的理論易于求得由式(6)定義的裂尖位移不連續(xù)函數(shù)引起的在點(diǎn)(x，y)處的位移場和應(yīng)力場為

2.3 雜交位移不連續(xù)法的實(shí)施及其說明

Crouch等利用式(4)～(5)建立了常位移不連續(xù)邊界元法.可類似地利用式(7)～(8)針對(duì)裂尖單元建立邊界元方程，進(jìn)而把常位移不連續(xù)單元和裂尖位移不連續(xù)單元有機(jī)地結(jié)合在一起以建立一種對(duì)平面裂紋問題很有效的邊界元法.在該邊界元方法的實(shí)施過程中，左、右裂尖位移不連續(xù)單元分別置于裂紋的左、右裂尖處，而常位移不連續(xù)單元?jiǎng)t分布于除了裂尖位移不連續(xù)單元占據(jù)的位置之外的整個(gè)裂紋面及其它邊界.稱這種邊界元法為雜交位移不連續(xù)法.

雜交位移不連續(xù)法不同于Scavia［9］提出的雜交邊界元法.當(dāng)利用Scavia提出的雜交邊界元法分析分支裂紋問題時(shí)，計(jì)算量是巨大的.而利用雜交位移不連續(xù)法，可以既簡單而又準(zhǔn)確地處理分支裂紋問題［10-11］.

Pan［12］指出“常位移不連續(xù)法只適于無限大板中裂紋問題”.本文發(fā)現(xiàn)雜交位移不連續(xù)法不僅適于分析無限大板裂紋問題［13-17］而且適于處理有限大板中復(fù)雜裂紋問題［18-19］.

孔洞裂紋問題、多裂紋問題一直是斷裂力學(xué)中的重要課題.利用雜交位移不連續(xù)法可以準(zhǔn)確地處理這些問題［20-22］.

2.4 應(yīng)力強(qiáng)度因子的計(jì)算公式

線彈性裂紋分析的主要目標(biāo)是確定裂紋尖端的應(yīng)力強(qiáng)度因子KI和 KII.基于裂尖附近的位移場，有

3 數(shù)值算例

下面通過2個(gè)數(shù)值算例說明本文提出的數(shù)值方法對(duì)分析無限大板中孔邊裂紋問題既簡單又非常有效.

例1 拉伸載荷作用下無限大板中橢圓孔分支裂紋問題，如圖1所示.對(duì)此裂紋問題，關(guān)于x軸對(duì)稱性條件是可以利用的.考慮下列情況

圖1 拉伸載荷作用下無限大板中橢圓孔分支裂紋問題

關(guān)于邊界單元的劃分，在1/2橢圓孔上劃分具有大致相同尺寸的400個(gè)單元，在其他邊界上按照所有單元具有大致相同尺寸的限制條件劃分.計(jì)算得到的歸一化應(yīng)力強(qiáng)度因子(SIFs nor-如表1所示.為了比較起見，表1列出了文獻(xiàn)［1］報(bào)道的數(shù)值計(jì)算結(jié)果.由表1可見本數(shù)值結(jié)果與文獻(xiàn)［1］報(bào)道的數(shù)值計(jì)算結(jié)果非常一致.

表1 拉伸載荷作用下無限大板中橢圓孔分支裂紋問題的歸一化應(yīng)力強(qiáng)度因子

例2 內(nèi)部壓力作用下無限大板圓孔分支裂紋問題，如圖2所示.對(duì)此裂紋問題，關(guān)于x軸和y軸對(duì)稱性條件是可以利用的.考慮下列情況為

圖2 內(nèi)部壓力載荷作用下無限大板中圓孔分支裂紋問題

關(guān)于邊界單元的劃分，在1/4圓孔上劃分具有大致相同尺寸的200個(gè)單元，在其他邊界上按照所有單元具有大致相同尺寸的限制條件劃分.計(jì)算得到的歸一化應(yīng)力強(qiáng)度因子(SIFs normalized如表2所示.為了便于比較，表2列出了文獻(xiàn)［1］報(bào)道的數(shù)值計(jì)算結(jié)果.由表2可見本數(shù)值結(jié)果與文獻(xiàn)［1］報(bào)道的數(shù)值計(jì)算結(jié)果非常一致.

表2 內(nèi)部壓力作用下無限大板圓孔分支裂紋問題的歸一化應(yīng)力強(qiáng)度因子

4 數(shù)值結(jié)果與討論

利用提出的求解無限大板中孔-裂紋問題的數(shù)值方法，詳細(xì)分析了拉伸載荷作用下無限大板三角孔-裂紋問題，如圖3所示.為了方便起見，首先考慮dy=dx=d情況.若此種三角孔-裂紋問題的應(yīng)力強(qiáng)度因子在數(shù)學(xué)上用KIth(a/d)表示，則它與中心裂紋問題(長度是2a)的應(yīng)力強(qiáng)度因子(KIcc)的比率(FIth)為

圖3 拉伸載荷作用下無限大板中三角孔－裂紋問題的示意圖

式(10)稱為歸一化應(yīng)力強(qiáng)度因子(normalized SIFs).對(duì)于此種三角孔-裂紋問題，計(jì)算得到的歸一化應(yīng)力強(qiáng)度因子如圖4所示.

圖4 拉伸載荷作用下無限大板中三角孔－裂紋問題的歸一化應(yīng)力強(qiáng)度因子

為了討論方便引入一個(gè)參數(shù)ad=2a/d.從圖4可見:

1)隨著ad的增加(亦即隨著三角孔尺寸d的減小)，F(xiàn)Ith快速單調(diào)增加，并且在某一值ad=adm處，F(xiàn)Ith達(dá)到其最大值FIthm.在這里，adm=1.08，F(xiàn)Ithm=1.219 6.

2)ad的繼續(xù)增加導(dǎo)致FIth開始快速減小，隨后慢慢減小.當(dāng)ad足夠大時(shí)(亦即三角孔尺寸d相對(duì)于裂紋長度a足夠小時(shí))，F(xiàn)Ith幾乎等于1(亦即FIth幾乎等于KIcc).

因此，可以看出三角孔對(duì)源于它的裂紋的應(yīng)力強(qiáng)度因子具有放大效應(yīng)(amplifying effect)，但在本研究范圍內(nèi)，未發(fā)現(xiàn)屏蔽效應(yīng)(shielding effect).在趨勢上，應(yīng)該也存在屏蔽效應(yīng).

通過變化三角孔的長細(xì)比dy/dx，進(jìn)而研究了三角孔的幾何參數(shù)對(duì)應(yīng)力強(qiáng)度因子的影響.為此，又考慮了3種情況:dy/dx=0.5，dy/dx=1.5和dy/dx=2.0.在這里，三角孔-裂紋問題的應(yīng)力強(qiáng)度因子在數(shù)學(xué)上可仍然用KIth表示，令其與中心裂紋問題(長度是2a)應(yīng)力強(qiáng)度因子的比率(FIth)，則

也稱為歸一化應(yīng)力強(qiáng)度因子.對(duì)于dy/dx=0.5、dy/dx=2.0及dy/dx=1.0情況，計(jì)算得到的歸一化應(yīng)力強(qiáng)度因子如圖5所示.發(fā)現(xiàn)對(duì)于by/bx= 2.0情況，屏蔽效應(yīng)的確存在.在這里，重點(diǎn)討論它的放大效應(yīng).上面引入的無量綱參數(shù)adm和FIthm如表3所示.鑒于在已研究的dy/dx范圍內(nèi)，F(xiàn)Ithm隨著dy/dx明顯地變化，又考慮了2種情況:dy/ dx=2.5和dy/dx=3.0，其adm和FIthm計(jì)算結(jié)果也如表3所示.FIthm隨著dy/dx的變化如圖6所示.從表3和圖6可作結(jié)論為:隨著dy/dx增大FIthm而增加，在dy/dx=2.5時(shí)，F(xiàn)Ithm基本上達(dá)到其最大值.

圖5 拉伸載荷作用下無限大板中三角孔－裂紋問題的歸一化應(yīng)力強(qiáng)度因子

圖6 FIthm隨著dy/dx的變化

表3 無量綱參數(shù)adm和FIthm隨著dy/dx變化

5 結(jié)論

1)提出了求解無限大板孔邊裂紋應(yīng)力強(qiáng)度因子的一種數(shù)值方法，算例說明本數(shù)值方法既簡單又非常有效.

2)利用這種數(shù)值方法研究了拉伸載荷作用下無限大板中三角孔-裂紋問題.通過改變孔的幾何參數(shù)，與無限大板中心裂紋問題的應(yīng)力強(qiáng)度因子比較，揭示了孔的幾何參數(shù)對(duì)應(yīng)力強(qiáng)度因子的影響.發(fā)現(xiàn)孔對(duì)源于它的裂紋的應(yīng)力強(qiáng)度因子具有屏蔽效應(yīng)和放大效應(yīng).

3)隨著dy/dx增大，F(xiàn)Ithm而增加，在dy/dx= 2.5時(shí)，F(xiàn)Ithm基本上達(dá)到其最大值.

4)在 dy/dx=2.5，adm=1.4時(shí)，F(xiàn)Ithm= 1.341 9.

5)對(duì)于一個(gè)中心裂紋板，通過在裂紋板中割去一個(gè)三角孔，使得三角孔的尺寸滿足條件dy/dx= 2.5，adm=1.4，則可利用三角孔的放大效應(yīng)來降低裂紋板的破壞載荷高達(dá)34%.

［1］ MURAKAMI Y.Stress Intensity Factors Handbook［M］.New York:Pergamon Press，1987.

［2］BOWIE O L.Analysis of an infinite plate containing radial cracks originating at the boundary of an internal circular hole［J］.J Math Phys，1956(35):60-71.

［3］Jr.NEWMAN J C.An improved method of collocation for the stress analysis of cracked plates with various shaped boundaries［C］//National Acronactics and Space Administration.Washington，DC:IEEE，1971: 1-45.

［4］NISITANI H，ISIDA M.Simple procedure for calculating KIof a notch with a crack of arbitrary size and its application to non-propagating fatigue crack［C］//Proceedings of the Ioint JSME-SESA Conf On Experimental Mechanics.［S.l.］:［s.n.］，1982:150-155.

［5］MURAKAMI Y.A method of stress intensity factor calculation for the crack emanating from an arbitrarily shaped hole or the crack in the vicinity of an arbitrarily shape hole［J］.Trans Japan Soc Mech Engrs，1978，44(378):423-432.

［6］BUCKNER H F.The propagation of cracks and the energy of elastic deformation［J］.ASME Appl Mech，1958(80):1225-1230.

［7］YAN X.An efficient and accurate numerical method of SIFs calculation of a branched crack［J］.ASME J Appl Mech，2005，72(3):330-340.

［8］CROUCH S L，STARFIELD A M.Boundary Element Method in Solid Mechanics，with Application in Rock Mechanics and Geological Mechanics［M］.London: Geore Allon＆Unwin，1983:1-220.

［9］SCAVIA C.A numerical technique for the analysis of cracks subjected to normal compressive stresses［J］. Int J Num Methods Eng，1992(33):929-942.

［10］YAN X.Stress intensity factors for asymmetric branched cracks in plane extension by using crack tip displacement discontinuity elements［J］.Mechanics Research Communications，2005，32(4):375-384.

［11］YAN X.A numerical analysis of stress intensity factors at bifurcated cracks［J］.Engineering Failure Analysis，2006，13(4):629-637.

［12］PAN E.A general boundary element analysis of 2-D linear elastic fracture mechanics［J］.Int J Fract，1997(88):41-59.

［13］YAN X.Analysis of the interference effect of arbitrary multiple parabolic cracks in plane elasticity by using a new boundary element method［J］.Computer Methods in Applied Mechanics and Engineering， 2003，192(47/48):5099-5121.

［14］YAN X.A numerical analysis of perpendicular cracks under general in-plane loading with a hybrid displacement discontinuity method［J］.Mechanical Research Communications，2004(31):175-183.

［15］YAN X.Interaction of arbitrary multiple cracks in an infinite plate［J］.J Strain Analysis for Engineering Design，2004(39):237-244.

［16］YAN X.An effective boundary element method for analysis of crack problems in a plane elastic plate［J］. Appl Math Mech，2005(26):814-822.

［17］YAN X.A numerical analysis of cracks emanating from a square hole in a rectangular plate under biaxial loads［J］.Eng Fracture Mech，2004(71):1615-1623.

［18］YAN X.A numerical analysis of cracks emanating from an elliptical hole in a 2-D elasticity plate［J］.European J of Mechanics-A，2006，(25):142-153.

［19］YAN X.Cracks emanating from circular hole or square hole in rectangular plate in tension［J］.Eng Fracture Mech，2006，73(12):1743-1754.

［20］YAN X.Microdefect interacting with a finite main crack［J］.J Strain Analysis for Engineering Design，2005(40):421-430.

［21］YAN X.A brief evaluation of approximation methods for microcrack shielding problems［J］.ASME J Appl Mech，2006，73(4):694-696.

［22］YAN X.An effective numerical approach for multiple void-crack interaction［J］.ASME J Appl Mech，2006，73(4):525-535.