劉 暉，申韶丁，雷 電

(武漢理工大學(xué) 道路橋梁與結(jié)構(gòu)工程湖北省重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430070)

0 引 言

焊接結(jié)構(gòu)是工程中應(yīng)用非常廣泛的結(jié)構(gòu)形式之一，焊接結(jié)構(gòu)在焊接過(guò)程中不可避免地會(huì)存在初始裂紋缺陷，具有焊接初始裂紋缺陷的構(gòu)件或結(jié)構(gòu)疲勞裂紋擴(kuò)展非常迅速，進(jìn)而引起結(jié)構(gòu)破壞[1]。因此，研究焊接結(jié)構(gòu)存在初始裂紋擴(kuò)展行為，對(duì)預(yù)防焊接結(jié)構(gòu)破壞，保證結(jié)構(gòu)使用安全具有重要意義。

對(duì)于裂紋擴(kuò)展模擬分析，許多學(xué)者做了大量的研究，取得了一定的成果，并提出了很多模擬裂紋擴(kuò)展的方法[2-5]，其中有限元法是目前運(yùn)用最為廣泛并且發(fā)展最為迅速的數(shù)值計(jì)算方法[6-7]。解德等[8]采用ABAQUS有限元分析軟件實(shí)現(xiàn)了基于有限元的模擬裂縫擴(kuò)展方法。李長(zhǎng)安[9]將基于有限元的虛擬裂紋閉合法(VCCT)應(yīng)用于疲勞裂紋壽命預(yù)測(cè)，并降低了網(wǎng)格敏感性，使得在粗糙網(wǎng)格下也能獲得比較精確的結(jié)果。劉艷萍等[10]基于VCCT法，采用Zigzag模型近似模擬裂紋前沿，分析了具有不同裂紋形式的簡(jiǎn)單幾何體，獲得了比較理想的結(jié)果。由于有限元法及其衍生方法是一種基于網(wǎng)格的方法，在劃分單元時(shí)要防止單元畸變或纏結(jié)現(xiàn)象，當(dāng)結(jié)構(gòu)出現(xiàn)裂紋擴(kuò)展時(shí)，網(wǎng)格會(huì)發(fā)生扭曲甚至畸變，裂紋擴(kuò)展的每一步都要重新劃分網(wǎng)格，這要花費(fèi)大量的前處理時(shí)間。模擬平面問(wèn)題中任意放置的直線裂紋和彎曲裂紋的彎折擴(kuò)展時(shí)，單元網(wǎng)格會(huì)產(chǎn)生過(guò)度扭曲，或造成計(jì)算中斷，不能真正跟蹤裂紋擴(kuò)展路徑。

無(wú)網(wǎng)格法作為一種新興的數(shù)值計(jì)算方法，將求解域基于節(jié)點(diǎn)離散，無(wú)需劃分單元，并且在裂紋擴(kuò)展時(shí)僅需局部和少量人為介入[11-13]。

袁振等[14]提出將無(wú)網(wǎng)格伽遼金法用于復(fù)合型裂紋的擴(kuò)展模擬并估算其疲勞壽命。陳建等[15]通過(guò)無(wú)網(wǎng)格伽遼金法對(duì)含邊沿裂紋功能梯度材料板進(jìn)行應(yīng)力強(qiáng)度因子的計(jì)算。傳統(tǒng)無(wú)網(wǎng)格法處理不連續(xù)問(wèn)題時(shí)，對(duì)裂紋幾何形狀難以描述，用透視法等一些方法處理裂紋隔離作用時(shí)，人為增加了不連續(xù)性，計(jì)算過(guò)程復(fù)雜，不易推廣到三維問(wèn)題上。

本文提出采用無(wú)網(wǎng)格和水平集耦合法來(lái)對(duì)焊接結(jié)構(gòu)存在初始裂紋幾何形狀進(jìn)行描述，并模擬裂紋擴(kuò)展。水平集法是一種追蹤界面移動(dòng)的數(shù)值模擬方法[16]，也是基于節(jié)點(diǎn)離散的方法，因此兩者可以很自然地實(shí)現(xiàn)耦合，水平集法采用2個(gè)相互正交的水平集函數(shù)對(duì)裂紋幾何形狀進(jìn)行描述及裂尖位置進(jìn)行定位。在無(wú)網(wǎng)格法的近似函數(shù)中增加了裂尖處的Westergard擴(kuò)展項(xiàng)和跳躍函數(shù)Heaviside[17]，使計(jì)算點(diǎn)不受裂紋隔離影響，裂紋擴(kuò)展時(shí)只需更新水平集函數(shù)，避免了無(wú)網(wǎng)格法計(jì)算每一步裂紋擴(kuò)展時(shí)都需要在裂尖位置進(jìn)行節(jié)點(diǎn)加密和人為增加的不連續(xù)項(xiàng)，保證了剛度矩陣的帶狀和稀疏，并且節(jié)點(diǎn)劃分只在裂尖部分的帶狀區(qū)域，既保證了計(jì)算精度又提高了計(jì)算效率。最后，本文以焊接節(jié)點(diǎn)為環(huán)狀形式截面且存在初始焊接裂紋為研究對(duì)象，編制了基于極坐標(biāo)的無(wú)網(wǎng)格和水平集耦合法程序，分析其裂紋擴(kuò)展，并跟蹤了裂紋擴(kuò)展路徑。

1 焊接結(jié)構(gòu)裂紋擴(kuò)展分析理論及流程

1.1 裂紋擴(kuò)展判別依據(jù)

最大周向拉應(yīng)力強(qiáng)度因子理論是判斷復(fù)合型裂紋斷裂應(yīng)用最為廣泛的準(zhǔn)則，該理論計(jì)算應(yīng)力強(qiáng)度因子簡(jiǎn)單，易于實(shí)現(xiàn)，被廣泛應(yīng)用于工程中，故本文采用此準(zhǔn)則作為裂紋擴(kuò)展判別依據(jù)。裂紋擴(kuò)展角θc應(yīng)滿足

KⅠsin(θc)+KⅡ[3cos(θc)-1]=0

(1)

式中：KⅠ，KⅡ分別為Ⅰ型和Ⅱ型裂紋的應(yīng)力強(qiáng)度因子(SIF)。

求解式(1)可得

(2)

式中：sign(·)為符號(hào)函數(shù)。

根據(jù)線彈性力學(xué)理論，裂紋尖端的應(yīng)力通常是趨近于無(wú)窮大的，不能作為判斷裂紋失效的準(zhǔn)則，因此采用應(yīng)力強(qiáng)度因子的概念來(lái)反映裂尖奇異應(yīng)力場(chǎng)的強(qiáng)弱。

復(fù)合型裂紋尖端附近極坐標(biāo)應(yīng)力場(chǎng)θ方向分量表達(dá)式為

(3)

當(dāng)θ=θc時(shí)

(4)

式中：r,θ為以裂紋尖為坐標(biāo)原點(diǎn)的極坐標(biāo)系；Kθmax為最大周向拉應(yīng)力強(qiáng)度因子，當(dāng)Kθmax大于材料斷裂韌度KIC時(shí)，裂紋開(kāi)始擴(kuò)展。

1.2 裂紋擴(kuò)展分析步驟及流程

裂尖處的SIF大小表示裂尖處應(yīng)力場(chǎng)的強(qiáng)弱，模擬裂紋擴(kuò)展就是計(jì)算裂尖處的SIF。因此，本文通過(guò)以下步驟計(jì)算SIF，最終實(shí)現(xiàn)對(duì)焊接節(jié)點(diǎn)的裂紋擴(kuò)展路徑跟蹤。

步驟1：采用無(wú)網(wǎng)格Garlerkin法(EFGM)和水平集法耦合建立模型，將模型求解域基于節(jié)點(diǎn)離散，通過(guò)計(jì)算節(jié)點(diǎn)水平集函數(shù)將節(jié)點(diǎn)劃分為常規(guī)節(jié)點(diǎn)、階躍擴(kuò)展節(jié)點(diǎn)和裂尖擴(kuò)展節(jié)點(diǎn)。

水平集函數(shù)Ψ(x,t)采用符號(hào)距離函數(shù)表示

(5)

式中：x，t分別為計(jì)算點(diǎn)的坐標(biāo)向量和時(shí)間變量；Γc為裂紋面；xΓ為位于Γc上的計(jì)算點(diǎn)。

如果計(jì)算點(diǎn)x位于Γc的上方，則水平集函數(shù)Ψ(x，t)為正，否則為負(fù)。

移動(dòng)界面Ψ(x，t)的演化方程為

Ψt+F‖Ψ‖=0

(6)

式中：Ψt為t時(shí)刻的界面；F為界面上點(diǎn)x∈Γc(t)在界面外法線方向的移動(dòng)速度向量；為梯度算子。

初始裂紋波前水平集函數(shù)為

(7)

水平集法通常采用裂紋面水平集函數(shù)和波前水平集函數(shù)描述裂紋位置及形狀，裂紋面水平集函數(shù)Ψ(x,t)=0表示裂紋及其延長(zhǎng)部分，初始裂紋由Ψ(x,t)=0和波前水平集函數(shù)φ≤0確定，那么裂紋面通過(guò)計(jì)算上述2個(gè)函數(shù)就可對(duì)節(jié)點(diǎn)進(jìn)行劃分，階躍擴(kuò)展節(jié)點(diǎn)集合M的表達(dá)式為

M={n∈Q∶φ<0且Ψ-rd≤0}

(8)

式中：rd為節(jié)點(diǎn)影響域半徑；Q為求解區(qū)域內(nèi)離散的節(jié)點(diǎn)集合；n為任意節(jié)點(diǎn)。

裂尖擴(kuò)展節(jié)點(diǎn)集合K的表達(dá)式為

(9)

當(dāng)節(jié)點(diǎn)既為階躍擴(kuò)展節(jié)點(diǎn)又為裂尖擴(kuò)展節(jié)點(diǎn)時(shí)，優(yōu)先按裂尖擴(kuò)展節(jié)點(diǎn)處理。

步驟2：根據(jù)EFGM和水平集耦合法推導(dǎo)整體剛度矩陣、荷載和位移向量，求解節(jié)點(diǎn)的應(yīng)力場(chǎng)和位移場(chǎng)。

根據(jù)線彈性斷裂力學(xué)的Westergaard解及EFGM的單位分解特性，耦合法的不連續(xù)位移場(chǎng)函數(shù)uI(x,t)可表示為

(10)

(11)

(12)

無(wú)網(wǎng)格和水平集耦合法的離散方程為

(13)

整體剛度矩陣的元素kij表達(dá)式為

(14)

式(14)上標(biāo)表示對(duì)應(yīng)的節(jié)點(diǎn)自由度，具體形式表達(dá)為

(15)

式中：B為應(yīng)變矩陣；D為彈性矩陣；Ω為計(jì)算域。

邊界項(xiàng)剛度矩陣的元素gij表達(dá)式為

(16)

荷載向量f中元素表達(dá)式為

(17)

式(17)上標(biāo)表示對(duì)應(yīng)的節(jié)點(diǎn)自由度，具體形式為

(18)

(19)

(20)

步驟3：根據(jù)式(13)計(jì)算得到節(jié)點(diǎn)的位移場(chǎng)，進(jìn)而得到應(yīng)力場(chǎng)，然后采用相互作用積分法[18]計(jì)算裂紋尖端區(qū)域的SIF。通過(guò)最大周向應(yīng)力準(zhǔn)則確定開(kāi)裂步長(zhǎng)，計(jì)算裂紋開(kāi)裂角。

步驟4：當(dāng)裂紋在荷載作用下擴(kuò)展時(shí)，通過(guò)水平集更新算法對(duì)節(jié)點(diǎn)水平集函數(shù)進(jìn)行更新[19]，從而對(duì)裂紋幾何形狀和裂尖位置進(jìn)行追蹤，實(shí)現(xiàn)對(duì)裂紋擴(kuò)展路徑的精確模擬。

以二維問(wèn)題為例來(lái)說(shuō)明具體的實(shí)現(xiàn)步驟：

(21)

Ψm+1=±|(x-xi)(F/‖F(xiàn)‖)|=±|(x-xi)|·

(Fy/‖F(xiàn)‖)+(y-yi)|(Fx/‖F(xiàn)‖)|

(22)

式中：正負(fù)號(hào)表示裂紋兩側(cè)更新部分的符號(hào)與未更新部分的符號(hào)保持一致。

(23)

式中：Δαm為第m步的擴(kuò)展步長(zhǎng)或者增量。

2 數(shù)值算例

為了驗(yàn)證本文方法的適用性和程序設(shè)計(jì)的正確性，選取應(yīng)力強(qiáng)度因子手冊(cè)中經(jīng)典邊緣受拉裂紋算例進(jìn)行計(jì)算。板模型如圖1所示，平板寬度L=1 m，高度D=2 m，板的上表面承受均勻拉應(yīng)力σ=1 MPa，板的左側(cè)中間部分有一邊緣裂紋長(zhǎng)度a=0.45 m，彈性模量E=1×103Pa，泊松比ν=0.3，應(yīng)力強(qiáng)度因子的精確解為

圖1 邊緣受拉裂紋模型

(24)

采用本文方法，首先將求解域離散為12×24個(gè)均勻布置的節(jié)點(diǎn)，則背景網(wǎng)格分布為11×23，每個(gè)背景網(wǎng)格內(nèi)設(shè)置6×6的高斯積分點(diǎn)，基函數(shù)選用二次基函數(shù)P=[1,x,y]T，權(quán)函數(shù)選用3次樣條權(quán)函數(shù)，影響域半徑選為1.7倍的最大相鄰節(jié)點(diǎn)間距rd，采用水平基函數(shù)將節(jié)點(diǎn)劃分為階躍擴(kuò)展節(jié)點(diǎn)和裂尖擴(kuò)展節(jié)點(diǎn)，節(jié)點(diǎn)分布如圖2所示。

圖2 節(jié)點(diǎn)分布

節(jié)點(diǎn)變形如圖3所示，為繞裂尖位置發(fā)生受拉變形。模型應(yīng)力云圖如圖4所示，在裂尖位置應(yīng)力最大，隨著遠(yuǎn)離裂尖位置應(yīng)力逐漸減小。采用本文方法計(jì)算結(jié)果與解析解對(duì)比如圖5所示，本文解與解析解基本一致，說(shuō)明本文方法具有較高精度，驗(yàn)證了本文方法的可行性和程序設(shè)計(jì)的正確性。

圖3 節(jié)點(diǎn)變形

圖4 應(yīng)力云圖(單位:MPa)

圖5 不同裂紋長(zhǎng)度對(duì)應(yīng)的SIF

焊接節(jié)點(diǎn)的一種形式是管球焊接節(jié)點(diǎn)或鋼管對(duì)焊節(jié)點(diǎn)，因此，本文以焊接結(jié)構(gòu)存在初始裂紋、具有環(huán)狀橫截面為研究對(duì)象分析其裂紋擴(kuò)展行為?；跇O坐標(biāo)系建立了如圖6所示的模型，材料為Q345鋼材，彈性模量E=2.06×1011Pa，泊松比ν=0.3，焊縫標(biāo)準(zhǔn)為二級(jí)焊縫，在外邊界上受切向均布荷載P=40 N。裂紋線為水平方向位置從A點(diǎn)(30，π)到B點(diǎn)(25，π)，裂紋長(zhǎng)度a=5 mm，裂尖坐標(biāo)為(25，π)。

圖6 焊接節(jié)點(diǎn)模型示意

采用均勻布點(diǎn)法在求解域內(nèi)進(jìn)行離散，外環(huán)半徑R為30 mm，內(nèi)環(huán)半徑r為20 mm，徑向方向等間距布置8個(gè)節(jié)點(diǎn)，弧度方向?yàn)槊扛籀?36等弧度布置72個(gè)節(jié)點(diǎn)。

影響域半徑過(guò)小，裂紋線附近的階躍擴(kuò)展節(jié)點(diǎn)不能完全覆蓋裂紋區(qū)域，對(duì)裂紋描述不夠精確；影響域半徑過(guò)大，裂尖擴(kuò)展節(jié)點(diǎn)過(guò)多，放大了裂尖區(qū)域的奇異性。因此，通過(guò)對(duì)比分析，本文選取2.2倍節(jié)點(diǎn)間距rd作為影響域半徑，裂尖部分的奇異場(chǎng)能夠得到足夠的節(jié)點(diǎn)位移近似函數(shù)加強(qiáng)，同時(shí)裂紋線附近的階躍擴(kuò)展節(jié)點(diǎn)也能夠完全覆蓋裂紋區(qū)域，使得計(jì)算點(diǎn)影響域不受裂紋隔離的影響，如圖7所示。

圖7 在2.2rd下的節(jié)點(diǎn)劃分

根據(jù)式(13)計(jì)算得到應(yīng)力場(chǎng)和位移場(chǎng)以后，采用相互作用積分法計(jì)算應(yīng)力強(qiáng)度因子，本文選取每一步裂紋擴(kuò)展增量da為a/5，da過(guò)大不適用于線彈性斷裂力學(xué)，會(huì)導(dǎo)致計(jì)算結(jié)果不精確；da過(guò)小，模擬裂紋擴(kuò)展時(shí)計(jì)算量大，并且計(jì)算所得SIF和θc不具有代表性。通過(guò)水平集更新算法更新裂尖位置及水平集函數(shù)，從而實(shí)現(xiàn)對(duì)裂紋擴(kuò)展的追蹤，得到最終的每一步擴(kuò)展增量下裂尖處的應(yīng)力強(qiáng)度因子，如圖7所示。所得裂紋擴(kuò)展角如表1所示，應(yīng)力強(qiáng)度因子隨裂紋擴(kuò)展增量的變化如圖8所示，裂紋擴(kuò)展路徑如圖9所示。

表1 裂紋擴(kuò)展對(duì)應(yīng)的SIF及擴(kuò)展角

圖8 SIF隨裂紋擴(kuò)展增量的變化

圖9 裂紋擴(kuò)展路徑

3 結(jié) 語(yǔ)

(1)本文方法避免了單純的無(wú)網(wǎng)格法處理不連續(xù)問(wèn)題時(shí)所增加的人為不連續(xù)項(xiàng)，能夠組裝成具有帶狀、稀疏性的整體剛度矩陣；對(duì)于裂紋擴(kuò)展區(qū)域加強(qiáng)點(diǎn)選取僅限制在裂紋附近的帶狀區(qū)域內(nèi)，簡(jiǎn)化了加強(qiáng)點(diǎn)選取和附加函數(shù)的建立，不會(huì)過(guò)多增加整體自由度，并且無(wú)需求解水平集更新算法的演化方程，因此計(jì)算過(guò)程簡(jiǎn)單，易于在程序上實(shí)現(xiàn)，提高了計(jì)算效率。

(2)本文方法的計(jì)算精度較高，結(jié)果光滑協(xié)調(diào)，且能很好地跟蹤裂紋擴(kuò)展路徑，能夠預(yù)測(cè)焊接結(jié)構(gòu)存在初始裂紋的裂紋擴(kuò)展及走向。

(3)本文方法計(jì)算量小、精度高，可以實(shí)現(xiàn)裂紋路徑跟蹤，易于推廣到三維節(jié)點(diǎn)模型，為預(yù)測(cè)復(fù)雜工程結(jié)構(gòu)在復(fù)雜荷載下的斷裂行為提供了一種可靠的計(jì)算途徑和模擬思路。