李 妍，孟廣偉，尹新生

(1.吉林大學(xué)機(jī)械科學(xué)與工程學(xué)院，長春130022;2.吉林建筑工程學(xué)院寒地綠色建筑技術(shù)工程研究中心，長春130118，lessle@tom.com)

隨著我國墻體材料改革和節(jié)能政策的貫徹執(zhí)行，混凝土多孔磚成為我國大力發(fā)展的替代粘土磚的主導(dǎo)產(chǎn)品之一.但是混凝土多孔磚作為一種新型的墻體材料，其砌體結(jié)構(gòu)在施工和使用過程中普遍存在著裂縫問題，而裂縫的存在，不僅影響結(jié)構(gòu)的正常使用，而且對其抗震的性能有較大的影響［1－2］.所以如何對裂縫進(jìn)行控制，使其在安全使用的許可范圍內(nèi)，防止或減少裂縫帶來的危害就變得尤為重要.目前對砌體結(jié)構(gòu)中裂縫(群)相互作用及其破裂問題的研究，大多集中在結(jié)構(gòu)試驗(yàn)研究階段［3］，文獻(xiàn)［4－5］對不同荷載作用下砌體結(jié)構(gòu)的性能進(jìn)行了實(shí)驗(yàn)研究，并得出了滯回曲線;文獻(xiàn)［6－7］對砌體結(jié)構(gòu)在軸力、彎矩作用下的性能進(jìn)行了研究，提出了抗彎、抗剪承載力的計(jì)算公式.但總體上說，砌體破壞機(jī)理的研究無論是國外還是國內(nèi)主要是集中在假設(shè)模型的基礎(chǔ)上的靜力性能分析和利用有限元提出一些簡化模型對砌體的性能進(jìn)行研究［8－9］.而有關(guān)混凝土多孔磚墻破裂模式的研究，由于所處理破壞問題的復(fù)雜性，很難對其從變形、破裂全過程進(jìn)行模擬.

自Kaplan首次將斷裂力學(xué)應(yīng)用于混凝土以來，準(zhǔn)脆性材料的斷裂理論［10－11］分析日益成為國內(nèi)外專家學(xué)者的研究熱點(diǎn)之一.但是由于砌體材料復(fù)雜的開裂機(jī)理，尋找并建立適用于砌體裂縫擴(kuò)展的斷裂準(zhǔn)則、新的斷裂參數(shù)以及相應(yīng)的計(jì)算方法，長期以來在國內(nèi)外卻鮮有報(bào)道.本文基于固體力學(xué)和線彈性斷裂力學(xué)基本原理，運(yùn)用修正的剪滯理論［12－13］，考慮到建立模型的方便和計(jì)算結(jié)果的精確設(shè)置了三個(gè)變異層，結(jié)合能量法則推導(dǎo)出了混凝土多孔磚砌體墻的等效斷裂韌度的計(jì)算模型［14］:針對相關(guān)的試驗(yàn)數(shù)據(jù)，分析了對應(yīng)標(biāo)準(zhǔn)子層高為1，2，3，4，5 mm各組模型的等效斷裂韌度，通過對結(jié)果的比較分析，提出了適合計(jì)算的合理的標(biāo)準(zhǔn)子層高，同時(shí)驗(yàn)證了等效斷裂韌度是與試件尺寸無關(guān)的斷裂參數(shù)，從而為研究混凝土多孔磚砌體墻的開裂與破壞的預(yù)測預(yù)報(bào)提供了新的方法，在工程上具有較強(qiáng)的適用性.

1 剪滯分析模型

根據(jù)混凝土多孔磚砌體墻的實(shí)際受力情況，建立如圖1所示的墻體模型，假定上下表面受均布壓力q的作用，厚度為t，高度為h，含初始裂縫長度為a0.為研究裂縫擴(kuò)展中的應(yīng)力重分布，建立圖2所示分層剪滯模型，即將模型沿著其高度方向劃分為n個(gè)子層，其中無裂縫磚區(qū)、無裂縫砂漿區(qū)、裂縫區(qū)、無裂縫砂漿區(qū)、無裂縫磚區(qū)(詳見圖3)的子層數(shù)分別為p+1、q、r+1、s和np－q－r－s－2層，各子層的標(biāo)準(zhǔn)高度為為了計(jì)算方便，將每層磚的第1子層，每層砂漿的第1子層分別設(shè)為具有特殊高度的子層，其子層高度分別為db0，dm0，這樣就可以使得模型中不能整除的部分加到這幾個(gè)特殊高度的子層上，即設(shè)置變異層.為了更為合理地設(shè)置各區(qū)子層高度，將裂縫區(qū)中間層設(shè)置為變異層，其子層高度為d0.

假定每一子層在垂直y方向的截面上僅有正應(yīng)力，上下表面僅有切應(yīng)力，則各子層的位移ui應(yīng)該僅為x的函數(shù)(沿x軸正向?yàn)槲灰频恼?.取長為dx的第i子層的微段，其受力如圖4所示.

圖1 混凝土多孔磚砌體墻模型及受力分析圖

圖2 基于豎向灰縫的分層剪滯模型

圖3 基于豎向灰縫擴(kuò)展的有限子層剪滯模型

圖4 第i子層微段的受力圖

其中:qi(x)為子層軸向壓力，pi(x)為由橫向變形對縱向變形的影響引起的相應(yīng)子層軸向拉力，τi－1(x)和τi(x)分別為上下層間切應(yīng)力.令第i子層的水平位移為ui(x)，豎向位移為vi(x).磚和砂漿的應(yīng)力－應(yīng)變關(guān)系均采用雙線性關(guān)系，通過切應(yīng)力來考慮磚和砂漿在分界處變形的相互影響，可以分別寫出磚和砂漿的pi(x)及τi(x)的表達(dá)式.

磚的表達(dá)式為

砂漿的表達(dá)式為

式中:Eb和Gb、Em和Gm分別為混凝土多孔磚和砂漿的彈性模量和切變模量;A=td為子層截面面積，其中d為標(biāo)準(zhǔn)子層厚度.

利用橫向變形與縱向變形的關(guān)系式

可以列出所有子層微段的剪滯平衡方程.

可見，第1至p+1子層微段的剪滯平衡方程組為

其中:A1=tdb0，列出第p+2至p+q+1子層微段的剪滯平衡方程組為

其中:A2=tdm0，需要注意的是，為便于利用對稱性計(jì)算，在采用本計(jì)算模式時(shí)，中間裂縫區(qū)的子層數(shù)r應(yīng)該為偶數(shù).在一個(gè)磚厚的砂漿裂縫區(qū)中，由于斷裂處的砂漿應(yīng)力為零，所以該區(qū)剪滯方程中彈性模量和泊松比均取自磚，即第p+q+2至p+ q+r+2子層微段的剪滯平衡方程組為

其中:A0=td0，同理，可得到各子層微段的剪滯平衡方程.

2 剪滯方程的求解

為便于求解，考慮引入無量綱參數(shù)，并根據(jù)相應(yīng)的邊界條件來求解上式的n階非齊次線性方程組，從而得到各子層的位移分布，則各子層的正應(yīng)力為

在ξ=0處，裂尖最鄰近子層的拉伸應(yīng)力集中因子則由下式求得:

采用有限子層模型模擬裂縫擴(kuò)展，假設(shè)裂縫擴(kuò)展長度a'為標(biāo)準(zhǔn)子層高度，無裂縫區(qū)第p+q+ 1層在x=0處的正應(yīng)力為σp+q+1(0)，則該混凝土多孔磚墻應(yīng)變能釋放率(即裂縫擴(kuò)展單位長度系統(tǒng)應(yīng)變能的下降值)為

由線彈性斷裂力學(xué)可知，在平面應(yīng)力條件下，

3 等效斷裂韌度的確定

由式(1)—(3)得到該混凝土多孔磚砌體墻等效斷裂韌度的解析表達(dá)式為

4 試驗(yàn)結(jié)果與分析

該試驗(yàn)采用吉林省祝焰新型建材有限公司生產(chǎn)的KP1型混凝土多孔磚，其材料配比分別為:水泥:18%;粉煤灰:19%;煤矸石:11%;中砂: 17%;5～10 mm的碎石:35%，并按照 GB/T 2542—2003的測試條件及要求測得該批磚的強(qiáng)度等級的平均值為fbm=12.3 MPa，材料常數(shù)Eb=1.65×104MPa，μb=0.22;吉林亞泰集團(tuán)鼎鹿牌水泥配制的水泥砂漿，測得其強(qiáng)度等級的平均值為 fmm=9.6 MPa，材料常數(shù) Em=0.4× 104MPa，μm=0.25，11個(gè)試件的尺寸分別為:5個(gè)L1系列(長度為240 mm)，3個(gè)L2系列(長度為490 mm)，3個(gè)L3系列(長度為740 mm)，其基本參數(shù)及荷載值如表1所示.

根據(jù)本文所建立的基于豎向灰縫模型的有限子層剪滯模型，對該模型所劃分的層數(shù)和層高如表2所列.利用有限元軟件ANSYS，在表3中列出了按本文計(jì)算的對應(yīng)標(biāo)準(zhǔn)子層高d分別為1，2，3，4，5 mm的等效斷裂韌度KC的解析解及由試驗(yàn)所得到的數(shù)值解.

表1 試件的基本參數(shù)及荷載值

表2 有限子層剪滯模型的基本參數(shù)

表3 混凝土多孔磚砌體墻等效斷裂韌度的計(jì)算值 MPa·m1/2

由表3可以得出，d分別為1，2，3，4，5 mm時(shí)，KC的平均值分別為0.302 8，0.301 7、0.300 6、0.300 0、0.299 1、0.301 8 MPa·m1/2，KC的均方差分別為0.010 1、0.010 3、0.010 6、0.010 0、0.009 8，KC變異系數(shù)分別為0.033 3、0.034 1、0.035 3、0.033 3、0.032 8.而實(shí)驗(yàn)所得到斷裂韌度的平均值、均方差、變異系數(shù)分別為0.301 8 MPa·m1/2、0.011 5、0.038 1.

可見，混凝土多孔磚砌體墻的等效斷裂韌度是一個(gè)不隨試件尺寸變化的材料常數(shù)，可以作為控制其斷裂的材料參數(shù).而且，用本文方法計(jì)算的同一試件的等效斷裂韌度存在隨層數(shù)增多而逐漸增大的趨勢，這是因?yàn)楸疚募僭O(shè)的裂縫擴(kuò)展長度為標(biāo)準(zhǔn)子層高，隨著層數(shù)的增加，標(biāo)準(zhǔn)子層高逐漸減小，使得應(yīng)變能釋放率逐漸趨于零.

5 結(jié)論

1)基于豎向灰縫模型推導(dǎo)出了求解混凝土多孔磚砌體墻等效斷裂韌度的解析計(jì)算模式，并通過相關(guān)實(shí)驗(yàn)進(jìn)行了數(shù)值分析，可知:同相關(guān)試驗(yàn)對應(yīng)的數(shù)值解相比，解析解的均方差和變異系數(shù)更小，因此本文所建立的求解混凝土多孔磚砌體墻等效斷裂韌度的解析計(jì)算模式具有很好的魯棒性，是求解其等效斷裂韌度的一種新的、實(shí)用的解析計(jì)算方法.

2)從理論上驗(yàn)證了基于豎向灰縫模型的混凝土多孔磚砌體墻等效斷裂韌度是一個(gè)與試件尺寸無關(guān)的斷裂參數(shù).

3)通過分區(qū)設(shè)置變異層可以使得所建立的剪滯分層模型更加方便和精確，從而為剪滯理論在砌體斷裂力學(xué)中的推廣應(yīng)用提供了理論依據(jù).

4)揭示了標(biāo)準(zhǔn)子層的層高是控制剪滯分層模型分析混凝土多孔磚砌體墻等效斷裂韌度的計(jì)算敏感性的一個(gè)重要因素，建議標(biāo)準(zhǔn)子層高度取為3～5 mm，這樣既能夠簡化計(jì)算，又能夠保證取得理想的計(jì)算精度.

［1］門進(jìn)杰，史慶軒，周琦.豎向不規(guī)則鋼筋混凝土框架結(jié)構(gòu)基于性能的抗震設(shè)計(jì)方法［J］.土木工程學(xué)報(bào)，2008，41(9):67－75.

［2］楊偉軍，林立，楊春俠.混凝土多孔磚砌體房屋抗震可靠度分析［J］.工程力學(xué)，2007，24(10):100－104.

［3］王鳳來，李新.配筋砌塊短肢砌體剪力墻偏壓承載力試驗(yàn)［J］.哈爾濱工業(yè)大學(xué)學(xué)報(bào)，2008，40(10):20－24.

［4］MIHA T，NIGEL C，SHIVAS S.Seismic behavior of masonry walls:experimental simulation［J］.Journal of Structural Engineering，ASCE，1996，122(9):1040－1047.

［5］ABRAMS D P，PAULSON T J.Modeling earthquake response of concrete masonry building structures［J］.ACI Structural Journal，1991，88(4):475－485.

［6］CHIOU Y J，TZENG J C，LIOU Y W.Experimental and analytical study of masonry infilled frames［J］.Journal of Structural Engineering，1999，125(6):1109－1125.

［7］ALSHEBANI M M，SINHA S N.Stress-strain characteristics of brick masonry under uniaxial cyclic loading［J］.Journal of Structural Engineering，ASCE，1999，125 (6):600－604.

［8］XING Jingtang，PRICE W G，CHEN Yigen.A numerical method to simulate nonlinear fluid-rigid structure interaction problems［J］.Acta Mechanica Solida Sinica，2005，26(4):373－383.

［9］王述紅，孫豁然.砌體受壓開裂過程數(shù)值模擬及其力學(xué)特性研究［J］.哈爾濱工業(yè)大學(xué)學(xué)報(bào)，2006，38(4): 644－648.

［10］吳智敏，王金來，徐世烺.基于虛擬裂縫模型的混凝土等效斷裂韌度［J］.工程力學(xué)，2000，17(1):99－104.

［11］徐世烺，周厚貴，高洪波.各種級配大壩混凝土雙K斷裂參數(shù)試驗(yàn)研究——兼對《水工混凝土斷裂試驗(yàn)規(guī)程》制定的建議［J］.土木工程學(xué)報(bào)，2006，39 (11):50－62.

［12］CARPINTERI A，CORNETTI P，BARPI F.Cohesive crack model description of ductile to brittle size-scale transition:dimensional analysis vs.renormalization group theory［J］.Engineering Fracture Mechanics，2003，70:1809－1839.

［13］曾慶敦.復(fù)合材料的細(xì)觀破壞機(jī)制與強(qiáng)度［M］.北京:科學(xué)出版社，2002.

［14］王鐸.斷裂力學(xué)［M］.哈爾濱:哈爾濱工業(yè)大學(xué)出版社，2006.