吳國(guó)榮，戎瑞亞，劉俊梅

(浙江海洋學(xué)院船舶與建筑工程學(xué)院，舟山 316004)

裂紋問(wèn)題廣泛存在于實(shí)際工程結(jié)構(gòu)中，特別是巖體工程中的軟弱結(jié)構(gòu)面、結(jié)構(gòu)工程及大型水工建筑物中的施工縫或開(kāi)裂面、機(jī)械工程中的接觸面和航天航空工程中復(fù)合材料的層面等均可能出現(xiàn)位移不連續(xù)現(xiàn)象，其中有些裂紋危及結(jié)構(gòu)的安全，甚至導(dǎo)致突發(fā)的災(zāi)難性后果，因此，如何有效地識(shí)別和檢測(cè)裂紋就具有十分重要的工程意義[1]。常見(jiàn)的裂紋識(shí)別方法有很多，如超聲檢測(cè)法、熒光劑檢測(cè)、電渦流與聲發(fā)射檢測(cè)等,這些方法在裂紋檢測(cè)上已取得了很好的效果，然而，這些傳統(tǒng)的裂紋檢測(cè)方法都存在著工藝復(fù)雜，難以定量分析等缺陷,很容易出現(xiàn)漏檢和誤判?；谡駝?dòng)特性的裂紋檢測(cè)法作為一種在線(xiàn)的無(wú)損檢測(cè)方法近年來(lái)越來(lái)越受到重視[2-7]。如何建立裂紋的力學(xué)模型準(zhǔn)確計(jì)算裂紋對(duì)結(jié)構(gòu)特性的影響是這種無(wú)損檢測(cè)方法獲得好的效果的前提條件。許多學(xué)者對(duì)此進(jìn)行了大量的研究，從解析解、半解析解到數(shù)值計(jì)算方面提出一系列解決策略[8-10]。

本文研究了含常開(kāi)裂紋梁的彎曲變形問(wèn)題，建立起含裂紋的梁?jiǎn)卧Ｐ?，基于斷裂力學(xué)中的含裂紋梁的能量表達(dá)式，利用卡氏定理導(dǎo)出了含裂紋梁?jiǎn)卧娜岫染仃嚕换诹Φ钠胶夥匠?，利用柔度矩陣和剛度矩陣的物理含義導(dǎo)出了含裂紋梁?jiǎn)卧膭偠染仃??；贛atlab軟件進(jìn)行了數(shù)值模擬運(yùn)算，給出了裂紋對(duì)梁撓度的影響。

1 柔度矩陣

對(duì)于含有裂紋的梁?jiǎn)卧?，由于裂紋處截面的轉(zhuǎn)角不連續(xù)，要構(gòu)造一個(gè)恰當(dāng)?shù)男魏瘮?shù)是非常困難的。斷裂力學(xué)對(duì)由于裂紋尖端的應(yīng)力集中而引起的應(yīng)變能的增加作了深入的研究，給出了其計(jì)算表達(dá)式。從而可以方便地導(dǎo)出含裂紋梁?jiǎn)卧娜岫认禂?shù)。

考慮一含裂紋梁?jiǎn)卧?，將右邊的單元作為?duì)含裂紋梁?jiǎn)卧a(chǎn)生主動(dòng)載荷，左邊的單元?jiǎng)t對(duì)其起固定約束作用，建立圖1所示的懸臂梁模型。不考慮裂紋時(shí)，該梁?jiǎn)卧淖冃文転椋?/p>

圖1 含裂紋梁?jiǎn)卧Ｐ虵ig.1 The model of the cracked beam element

其中，EI為梁的抗彎剛度，Q和M分別為右端截面上的剪力和彎矩，Le為單元的長(zhǎng)度。對(duì)于矩形截面梁，由于裂紋而額外增加的變形能為：

式中，a和b分別為裂紋的深度和截面的寬度，ГⅠ和ГⅡ分別為與張開(kāi)型和滑移型裂紋相對(duì)應(yīng)的應(yīng)力集中因子，其表達(dá)式分別為：

這里，a 為裂紋深度，b 和 h 分別為截面的寬和高，系數(shù) ξ=a/h。函數(shù) FⅠ(ξ)和 FⅡ(ξ)分別為[11]：

記C*為該梁?jiǎn)卧娜岫染仃嚕煽ㄊ隙ɡ?，可得?/p>

式中，

2 剛度矩陣

圖 2 為一含裂紋梁?jiǎn)卧?F=[Q(1)M(1)Q(2)M(2)]T和X=[v(1)θ(1)v(2)θ(2)]T分別為其節(jié)點(diǎn)載荷列陣和位移列陣。由力的平衡方程，可得：

圖2 含裂紋梁?jiǎn)卧狥ig.2 The cracked beam element

式中，

這里上標(biāo)“T”表示矩陣的轉(zhuǎn)置。注意到圖1模型中對(duì)節(jié)點(diǎn)1處位移的約束，則有：

式中，

由柔度系數(shù)的物理含義，有：

由(9)式、(10)式和(11)式，可得：

其中

為含裂紋梁?jiǎn)卧膭偠染仃?。?dāng)a=0時(shí)，則退化為：

這正是采用三次插值形函數(shù)的無(wú)裂紋梁?jiǎn)卧獎(jiǎng)偠染仃嚒?/p>

3 數(shù)值算例

分別就含常開(kāi)裂紋的簡(jiǎn)支梁和懸臂梁進(jìn)行了數(shù)值模擬運(yùn)算。梁的幾何尺寸及材料常數(shù)為：長(zhǎng)L＝1.0 m、高h(yuǎn)＝0.04 m及寬b＝0.02 m，彈性模量E=200 GPa。梁受一均布外載q=1.0 kN/m荷作用。表1給出的是無(wú)量綱最大撓度(最大撓度與無(wú)裂紋梁的比值)隨裂紋深度的變化。表2和表3分別給出的是ξ=0.2時(shí)簡(jiǎn)支梁和懸臂梁的無(wú)量綱最大撓度隨裂紋位置的變化。圖3給出是ξ=0.2和無(wú)裂紋梁的撓曲線(xiàn)。圖3和表1 ～3中的xc分別表示裂紋所在的位置對(duì)于簡(jiǎn)支梁為到支座的距離和對(duì)于懸臂梁為到固定端的距離。

表1 無(wú)量綱最大撓度隨裂紋深度的變化Tab.1 The variation of the dimensionless deflection with the crack depth

表2 簡(jiǎn)支梁的無(wú)量綱最大撓度隨裂紋位置的變化Tab.2 The variation of the dimensionless deflection with the crack position of the simple supported beam

表3 懸臂梁的無(wú)量綱最大撓度隨裂紋位置的變化Tab.3 The variation of the dimensionless deflection with the crack position of the cantilever beam

圖3 梁的撓曲線(xiàn)：實(shí)線(xiàn)為含裂紋梁，虛線(xiàn)為無(wú)裂紋梁Fig.3 The beam deflections:solid line for the cracked beam,dashed line for the uncracked beam

4 結(jié)論

本文提出了一種簡(jiǎn)便的含裂紋梁的有限元模型，基于斷裂力學(xué)理論和力的平衡方程，利用柔度矩陣和剛度矩陣的物理含義，導(dǎo)出了含裂紋梁?jiǎn)卧膭偠染仃嚒７謩e就簡(jiǎn)支梁和懸臂梁進(jìn)行了數(shù)值計(jì)算，從結(jié)果可以看出：隨著裂紋深度的增加，撓度也隨之增加；對(duì)于簡(jiǎn)支梁，裂紋在中間位置對(duì)撓度的影響最明顯，對(duì)于懸臂梁，裂紋在根部位置影響最明顯。

[1]范天佑.斷裂理論基礎(chǔ)[M].北京：科學(xué)出版社,2003.

[2]張敬芬,趙德有.工程結(jié)構(gòu)裂紋損傷振動(dòng)診斷的發(fā)展和展望[J].振動(dòng)與沖擊,2002,21(4):22-26.

[3]GEORGE D,GOUNARIS A,CHRIS A.Papadopoulos,Crack identi.cation in rotating shafts by coupled response measurements[J].Engineering Fracture Mechanics,2002(69):339-352.

[4]RIZOS P,ASPRAGATHOS N,DIMAROGONAS A.Identi.cation of crack location and magnitude in a cantilever beam from the vibration modes[J].Journal of Sound and Vibration,1990,138:381-388.

[5]李 兵,陳雪峰,胡 橋,等.基于小波有限元的懸臂梁裂紋識(shí)別[J].振動(dòng)工程學(xué)報(bào),2004,17(2):159-164.

[6]ZHENG D Y,FAN S C.Natural frequencies of a non-uniform beam with multiple cracks via modified Fourier series[J].Journal of Sound and Vibration,2001,242(4):701-717.

[7]ZHENG D Y,FAN S C,Natural frequency changes of a cracked Timoshenko beam by modified Fourier series[J].Journal of Sound and Vibration,2001,246(2):297–317.

[8]陳夢(mèng)成,余荷根,湯任基.三維裂紋問(wèn)題的高精度數(shù)值解法[J].固體力學(xué)學(xué)報(bào),2002,23(2):207-211.

[9]尹 峰.關(guān)于表面裂紋理論的研究與進(jìn)展(上)[J].沈陽(yáng)航空工業(yè)學(xué)院學(xué)報(bào),2003,20(3):88-91.

[10]BAMNIOS Y,DOUKA E,TROCHIDIS A,Crack identification in beam structures using mechanical impedance[J].Journal of Sound and Vibration,2002,256(2):287-297.

[11]NAHVI H,JABBARI M,Crack detection in beams using experimental modal data and finite element model[J].International Journal of Mechanical Sciences,2005,47(2):1 477-1 497.