基于擴(kuò)展有限元法及虛擬裂縫模型的混凝土斷裂過程區(qū)分析

1、    基于擴(kuò)展有限元法及虛擬裂縫模型的混凝土斷裂過程區(qū)分析    摘要：提出一種利用擴(kuò)展有限元法（extended finite element method，xfem）和虛擬裂縫模型對混凝土斷裂過程區(qū)（fracture process zone， fpz）進(jìn)行研究的方法.利用該方法可以求出裂縫擴(kuò)展過程中混凝土fpz的長度及位移和應(yīng)力分布.利用該方法對一個三點彎曲混凝土梁進(jìn)行研究，考察骨料粒徑、不同軟化律和不同初始裂縫長度對fpz的影響.關(guān)鍵詞：混凝土； 擴(kuò)展有限元法； 虛擬裂縫模型； 斷裂過程區(qū)； 裂縫擴(kuò)展； 軟化律： tv313 文獻(xiàn)標(biāo)志碼：bab

2、stract：a method is proposed to study the fracture process zone（fpz） of concrete based on extended finite element method（xfem） and a fictitious crack model. the length， displacement and stress distribution of fpz during the crack growth process can be obtained by the method. a threepoints bending con

3、crete beam is studied and the effect of the aggregate size， different softening laws and different initial crack length on fpz is studied.key words：concrete； extended finite element method； fictitious crack model； fracture process zone； crack growth； softening law0 引 言混凝土材料是一種準(zhǔn)脆性材料，其斷裂行為的重要特征是真實裂紋前方

4、存在斷裂過程區(qū)（fracture process zone，fpz）.在這一區(qū)域存在微裂縫、集料互鎖、粗糙表面的接觸和摩擦等非線性現(xiàn)象，從而使得混凝土的斷裂行為呈現(xiàn)非線性特性.13由于fpz的尺寸通常與混凝土構(gòu)件尺寸在同一數(shù)量級，采用傳統(tǒng)的線彈性斷裂力學(xué)分析混凝土結(jié)構(gòu)的斷裂問題已不適合.fpz是混凝土斷裂力學(xué)中一個極為重要的概念.混凝土fpz的研究與混凝土裂縫模型的研究緊密相關(guān)，按采用裂縫模型的不同，對混凝土fpz的研究方法可大致分為2類.1）基于黏聚裂縫模型的方法.這類模型適合于數(shù)值分析，其中適用于混凝土材料的黏聚裂縫模型又以虛擬裂縫模型4和裂縫帶模型5為代表.2）基于彈性等效裂縫模型的方法

5、.這類模型一般用于解析計算混凝土結(jié)構(gòu)斷裂問題，按所使用的彈性等效原則的不同，這類裂縫模型又包括兩參數(shù)斷裂模型6、尺寸效應(yīng)模型7、等效裂縫模型8和雙k斷裂模型9.由于解析裂縫模型的適用范圍有限，本文重點關(guān)注黏聚裂縫模型，特別是虛擬裂縫模型.虛擬裂縫模型將fpz簡化為一條虛擬的集中裂縫并將fpz的非線性行為集中到該虛擬裂縫上考慮.采用虛擬裂縫模型結(jié)合傳統(tǒng)有限元法，可以得到較為準(zhǔn)確的裂縫擴(kuò)展路徑及虛擬裂縫面的張開位移和應(yīng)力等物理量，但由于在傳統(tǒng)有限元法中采用單元邊界模擬裂縫，因此隨著裂縫的擴(kuò)展需不斷進(jìn)行網(wǎng)格重剖以保持單元邊界與裂縫擴(kuò)展路徑一致10，這極大地限制該方法的使用.為克服傳統(tǒng)有限元法分析斷裂

6、問題時遇到的網(wǎng)格重剖難題，belytschko等11提出擴(kuò)展有限元法（extended finite element method，xfem），基于單位分解定理，在傳統(tǒng)有限元連續(xù)位移場中引入能描述裂紋兩側(cè)位移間斷特性的非連續(xù)位移項，使裂紋的描述獨立于計算網(wǎng)格，因而無須隨著裂紋的擴(kuò)展不斷進(jìn)行網(wǎng)格重剖.擴(kuò)展有限元法結(jié)合虛擬裂縫模型被很多學(xué)者應(yīng)用于混凝土結(jié)構(gòu)的斷裂問題研究，成功再現(xiàn)許多試驗結(jié)果，顯示出很好的應(yīng)用前景.1214本文提出一種基于xfem和虛擬裂縫模型的混凝土fpz計算方法，利用該方法可以準(zhǔn)確計算混凝土結(jié)構(gòu)裂縫擴(kuò)展過程中虛擬裂縫的范圍、張開位移和應(yīng)力分布等.1 虛擬裂縫模型的擴(kuò)展有限元法1

7、.1 虛擬裂縫模型虛擬裂縫模型是在黏聚裂縫模型的基礎(chǔ)上發(fā)展而來的，該模型把混凝土裂縫分解為2部分（見圖1）：真實物理裂縫區(qū)（完全開裂區(qū)）和虛擬裂縫區(qū)（微裂區(qū)）.1前者代表宏觀的自由表面裂縫，裂縫表面無應(yīng)力作用；后者將帶狀fpz簡化為一條分離裂縫，即虛擬裂縫，虛擬裂縫表面由于微裂縫作用、集料互鎖等仍然可以承受應(yīng)力作用，且承擔(dān)應(yīng)力的大小與虛擬裂縫的張開位移相關(guān).虛擬裂縫模型概念清楚、模型簡潔，被認(rèn)為是基本的混凝土裂縫模型.虛擬裂縫模型采用如下假定：1）材料非線性完全集中于虛擬裂縫區(qū)，虛擬裂縫以外的區(qū)域為線彈性材料；2）虛擬裂縫表面存在的黏聚應(yīng)力與虛擬裂縫的張開位移w有關(guān)，一般可用w張拉軟化律曲線表

8、示.常用的線性和雙線性軟化律見圖2，曲線下方的面積即為混凝土材料的斷裂能gf=w00dw，其意義是形成單位面積裂縫所吸收的能量.虛擬模型的引入使得對fpz的研究可以通過對虛擬裂縫的研究進(jìn)行，通過研究虛擬裂縫張開位移及其表面的應(yīng)力分布即可獲取fpz的變形和應(yīng)力情況.1.2 xfem由于虛擬裂縫面上的黏聚應(yīng)力通常為虛擬裂縫張開位移的非線性函數(shù)，由式（5d）可知fcoh一般也為節(jié)點位移未知量的非線性函數(shù)，因此有限元方程系統(tǒng)的式（5a）和式（6）必須采用迭代方法求解.相關(guān)求解思路及裂縫擴(kuò)展過程模擬方法可參見文獻(xiàn)14，這里不再詳述. 2 虛擬裂縫計算方法虛擬裂縫的范圍及其上分布的黏聚應(yīng)力和虛擬裂縫張開位

9、移可以在有限元方程系統(tǒng)的迭代求解過程中予以求解，思路如下.虛擬裂縫擴(kuò)展路徑被單元邊界分割為許多候選虛擬裂縫段，見圖5a.在計算過程中，對每個候選虛擬裂縫段計算其高斯積分點上的裂縫張開位移w.1）首先由裂縫擴(kuò)展路徑計算裂縫與單元邊界交點在整體笛卡爾坐標(biāo)系下的坐標(biāo)xi和xj，見圖5b.求得高斯點的裂縫張開位移后，判斷候選虛擬裂縫段屬于黏聚裂縫還是真實裂縫，僅當(dāng)候選虛擬裂縫段上所有高斯點的張開位移w均滿足0ww0時，才認(rèn)為該段裂縫位于fpz，屬于黏聚裂縫；若有至少一個高斯點張開位移滿足w>w0，則該裂縫段成為自由表面裂縫.在每一次平衡迭代計算過程中均需對所有的候選虛擬裂縫段進(jìn)行上述計算.同時需

10、注意，一旦候選虛擬裂縫段變成真實裂縫的一部分，則在后續(xù)計算中該裂縫段將不再作為候選的虛擬裂縫段.在計算過程中，對每一個裂縫段采用3個高斯點計算其張開位移.當(dāng)系統(tǒng)平衡方程得到滿足時，黏聚裂縫段及其高斯點的張開位移即被確定，同時還可計算其長度和黏聚應(yīng)力分布，從而確定fpz的長度、張開位移和應(yīng)力分布.3 數(shù)值算例3.1 算例1 雙懸臂梁對于該算例，文獻(xiàn)13給出用偽邊界積分（pseudo boundary integral，pbi）方法得到的解，本算例以此作為比較的標(biāo)準(zhǔn).由于在偽邊界積分方法中裂紋的擴(kuò)展路徑是預(yù)先指定的，為保持一致性，本算例指定裂縫沿水平方向擴(kuò)展，即沿圖6中虛線所示路徑擴(kuò)展.有限元模型

11、采用119×59的均勻網(wǎng)格，按平面應(yīng)變問題進(jìn)行計算，裂縫每一步擴(kuò)展長度為4 mm.圖7和8分別給出裂縫擴(kuò)展至第9、34和58步時的虛擬裂縫張開位移和虛擬裂縫表面的應(yīng)力分布.圖中的裂縫張開位移和裂縫表面應(yīng)力分別用臨界裂縫張開位移w0和混凝土抗拉強(qiáng)度ft進(jìn)行無量綱化處理.計算結(jié)果顯示采用本文方法計算得到的虛擬裂縫張開位移和表面應(yīng)力分布與文獻(xiàn)13給出的結(jié)果完全一致，這說明本文方法和程序準(zhǔn)確.3.2 算例2 三點彎曲混凝土梁試件混凝土梁試件見圖9.梁的幾何尺寸為：厚度t=b=150 mm，l=600 mm.取文獻(xiàn)15給出的c3等級的混凝土進(jìn)行研究，其彈性模量和抗拉強(qiáng)度分別為e=34.65 g

12、pa，ft=3.5 mpa，計算時取=0.1.3.2.1 不同骨料粒徑及不同軟化律的影響考慮3種不同骨料粒徑對fpz的影響，對應(yīng)的最大骨料粒徑dmax分別為8，16和32 mm.不同的骨料粒徑主要影響混凝土材料的斷裂能，3種骨料粒徑對應(yīng)的斷裂能及采用線性和雙線性軟化律所需的參數(shù)見表1.對應(yīng)峰值載荷時虛擬裂縫表面的張開位移和應(yīng)力分布分別見圖10和11.不論是采用線性軟化律還是雙線性軟化律，不同骨料粒徑對應(yīng)的虛擬裂縫表面張開位移和表面應(yīng)力分布呈大致相同的規(guī)律，但fpz的長度卻不同.大致來說，骨料粒徑越大，對應(yīng)峰值載荷的fpz長度越大，特別是采用線性軟化曲線時這種現(xiàn)象非常明顯.需注意的是，當(dāng)載荷達(dá)到

13、峰值載荷時，混凝土fpz并未發(fā)展至飽和程度，這可以從圖10中x=0時w0看出.裂縫擴(kuò)展過程中fpz長度最大時虛擬裂縫表面的位移和應(yīng)力分布分別見圖12和13.不同骨料粒徑對應(yīng)的虛擬裂縫表面的位移分布遵循的基本規(guī)律相同，但采用線性軟化律和雙線性軟化律得到的虛擬裂縫表面應(yīng)力分布卻大不相同.這是因為此時fpz已充分發(fā)展.從圖12和13中可以看出：x=0時，虛擬裂縫張開位移已經(jīng)達(dá)到臨界值w0，而對應(yīng)的虛擬裂縫表面應(yīng)力為0；在整個fpz上位移從0至w0連續(xù)變化，由于雙線性軟化律曲線存在轉(zhuǎn)折點，因此此時的應(yīng)力分布曲線也必然出現(xiàn)轉(zhuǎn)折點.加載過程中fpz長度隨所加集中載荷的變化曲線見圖14.不論是采用線性軟化律

14、還是雙線性軟化律，當(dāng)載荷相同時，fpz長度隨骨料粒徑增加而增大；fpz長度在峰值載荷處并未達(dá)到最大值，在峰值載荷過后裂縫進(jìn)入非穩(wěn)定擴(kuò)展路徑并擴(kuò)展至一定階段時，fpz長度達(dá)到最大值.在骨料粒徑相同的情況下，采用雙線性軟化律得到的fpz長度大于采用線性軟化律得到的長度；采用雙線性軟化律得到的峰值載荷和fpz長度最大時對應(yīng)的載荷均小于采用線性軟化律時得到的結(jié)果.3.2.2 不同初始裂縫長度影響為考察不同初始裂縫長度對虛擬裂縫張開位移和表面應(yīng)力分布的影響，考慮線性軟化律情況，令初始裂縫長度從0.01b變化至0.25b.不同初始裂縫長度時虛擬裂縫張開位移和表面應(yīng)力分布分別見圖15和16.隨著初始裂縫長度

15、的增大，虛擬裂縫張開位移增大，而虛擬裂縫的表面應(yīng)力變小.需指出的是，在初始裂縫從0.01b變化至0.25b過程中，虛擬裂縫長度保持不變；若初始裂縫長度達(dá)到0.30b，則虛擬裂縫的長度會發(fā)生變化.4 結(jié) 論提出一種利用xfem結(jié)合虛擬裂縫模型對混凝土fpz進(jìn)行研究的方法，數(shù)值算例表明該方法準(zhǔn)確可靠.對一個三點彎曲混凝土梁試件在裂縫擴(kuò)展過程中的fpz特性進(jìn)行研究，結(jié)果如下.1）不同骨料粒徑對應(yīng)的混凝土fpz的位移和應(yīng)力分布遵循大致相同的規(guī)律.2）在外載荷達(dá)到峰值載荷時，混凝土fpz并未發(fā)展至飽和程度；在峰值載荷過后，當(dāng)結(jié)構(gòu)裂縫繼續(xù)非穩(wěn)定擴(kuò)展至某一階段時混凝土fpz方發(fā)展至飽和.3）在相同載荷水平下

16、，骨料粒徑越大，對應(yīng)的fpz長度越大.4）隨著結(jié)構(gòu)初始裂縫長度的增大，fpz的張開位移增大，而應(yīng)力卻減小.本文提出的方法原則上適用于任意復(fù)雜結(jié)構(gòu)裂縫擴(kuò)展過程中的fpz的研究，且可以考慮任意軟化律，但本文只對i型斷裂問題進(jìn)行驗證，對于ii型甚至復(fù)合型斷裂問題，本文方法是否有效尚需進(jìn)一步研究. 參考文獻(xiàn)：1徐世烺. 混凝土斷裂力學(xué)m. 北京： 科學(xué)出版社， 2011： 129169.2陳瑛， 姜弘道， 喬丕忠， 等. 混凝土黏聚開裂模型若干進(jìn)展j. 力學(xué)進(jìn)展， 2005， 35（3）： 377390.chen ying， jiang hongdao， qiao pizhong， et al. de

17、velopments on concrete cohesive modelj. advances in mechanics， 2005， 35（3）： 377390.3卿龍邦， 李慶斌， 管俊峰， 等. 基于虛擬裂縫模型的混凝土斷裂過程區(qū)研究j. 工程力學(xué)， 2012， 29（9）： 112116.qing longbang， li qingbin， guan junfeng， et al. study of concrete fracture process zone based on fictitious crack modelj. engineering mechanics， 2012，

18、 29（9）： 112116.4hillerborg a， moder m， peterson p e. analysis of crack propagation and crack growth in concrete by means of fracture mechanics and finite elementsj. cement and concrete research， 1976， 6（6）： 773782.5baant z p， oh b h. crack band theory for fracture of concretej. matériaux et con

19、struction， 1983， 16（3）： 155177.6jenq y s， shah s p. two parameter fracture model for concretej. journal of engineering mechanics， 1985， 111（10）： 12271241.7baant z p. size effect in blunt fracture： concrete， rock， metalj. journal of engineering mechanics， 1983， 110（4）： 518535.8nallathambi p， karihalo

20、o b l. determination of specimensize independent fracture toughness of plain concretej. magazine of concrete research， 1986， 38（135）： 6776.9xu s l， reinhardt h w. a simplified method for determining doublek fracture parameters for threepoint bending testsj. international journal of fracture， 2000， 1

21、04（2）： 181209.10楊慶生， 楊衛(wèi). 斷裂過程的有限元模擬j. 計算力學(xué)學(xué)報， 1997， 14（4）： 3338.yang qingsheng， yang wei. finite element simulation of fracture processj. chinese journal of computational mechanics， 1997， 14（4）： 3338.11belytschko t， gracie r， ventura g. a review of extended/generalized finite element methods for material modelingj. modelling and simulation in materials science a