孔芳芳，王國(guó)新（大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116024）

磚混結(jié)構(gòu)地震安全性評(píng)價(jià)

孔芳芳，王國(guó)新
（大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116024）

磚混結(jié)構(gòu)的延性較差和變形能力較小，在歷次大地震中均發(fā)生了嚴(yán)重的破壞，因此，在傳統(tǒng)定性分析磚混結(jié)構(gòu)典型破壞特征基礎(chǔ)上開(kāi)展定量研究，對(duì)提高其抗震性能具有重要的現(xiàn)實(shí)意義。以三層的磚混結(jié)構(gòu)房屋為研究對(duì)象，利用ABAQUS對(duì)房屋進(jìn)行了不同地震作用下的彈塑性分析。通過(guò)對(duì)房屋受拉損傷云圖以及受拉損傷參數(shù)值的對(duì)比分析，明確了最大層間位移角與受拉損傷面積百分比均可以作為衡量結(jié)構(gòu)破壞程度的定量指標(biāo)。結(jié)果表明：在6、7、8度小震，6、7度中震以及6度大震作用下，房屋結(jié)構(gòu)處于基本完好狀態(tài)，在8度中震和7度大震作用下處于中等破壞狀態(tài)，而在8度大震作用下趨于嚴(yán)重破壞狀態(tài)。

磚混結(jié)構(gòu)；時(shí)程分析；破壞

0 引言

磚混結(jié)構(gòu)是我國(guó)主要的建筑結(jié)構(gòu)形式之一，由于砌體材料抗彎、抗剪和抗拉強(qiáng)度都較低，加之其連接措施不夠完善和牢固，整體性較差，導(dǎo)致磚混結(jié)構(gòu)抗震能力較低，容易發(fā)生脆性破壞，歷次震害均表明，此類(lèi)結(jié)構(gòu)在地震作用下倒塌和嚴(yán)重破壞的比例較大，造成大量的生命和財(cái)產(chǎn)損失，因而研究磚混結(jié)構(gòu)的破壞特征，以定量的方式體現(xiàn)結(jié)構(gòu)不同破壞程度和破壞部位，對(duì)改進(jìn)和提高磚混結(jié)構(gòu)的抗震設(shè)防能力，以及震后損失評(píng)估具有重要意義。本文應(yīng)用有限元分析軟件 ABAQUS開(kāi)展磚混結(jié)構(gòu)動(dòng)力彈塑性時(shí)程分析，分析不同強(qiáng)度的地震對(duì)磚混結(jié)構(gòu)縱橫墻、結(jié)構(gòu)整體的損傷，以及對(duì)層間位移角的影響。通過(guò)對(duì)拉裂損傷進(jìn)行對(duì)比、層間位移角的大小與位移角限值對(duì)比，以及受拉損傷面積百分比三個(gè)特征參數(shù)，以定量的方式綜合分析和展示不同強(qiáng)度地震作用下結(jié)構(gòu)的破壞狀態(tài)，并給出其安全性評(píng)價(jià)。

1 ABAQUS 中模型及參數(shù)設(shè)定

1.1 數(shù)值模擬時(shí)模型選取

鄭妮娜[1]提出了一整套完整的運(yùn)用ABAQUS分析砌體結(jié)構(gòu)的思路，詳細(xì)地介紹了ABAQUS中砌體結(jié)構(gòu)、混凝土結(jié)構(gòu)以及鋼筋材料的各本構(gòu)關(guān)系和參數(shù)的設(shè)置，并證明了參數(shù)的可靠性；殷園園[2]應(yīng)用ABAQUS對(duì)文獻(xiàn)[1]振動(dòng)臺(tái)試驗(yàn)進(jìn)行數(shù)值模擬，取得了較好的效果；劉露[3]應(yīng)用ABAQUS建立了三維有限元模型研究地震作用下砌體結(jié)構(gòu)房屋發(fā)生破壞的主要影響因素；唐春美[4]運(yùn)用ABAQUS研究了砂漿強(qiáng)度以及橫墻開(kāi)洞對(duì)砌體結(jié)構(gòu)彈塑性抗震性能的影響，提出用受拉損傷面積比作為界定結(jié)構(gòu)破壞程度的衡量指標(biāo)；牟道燦[5]應(yīng)用ABAQUS對(duì)砌體結(jié)構(gòu)進(jìn)行有限元分析，總結(jié)出提高砌體結(jié)構(gòu)抗震性能的措施；劉杰[6]以文獻(xiàn)[1]的試驗(yàn)為研究基礎(chǔ)，運(yùn)用ABAQUS整體式建?？偨Y(jié)出砌體結(jié)構(gòu)發(fā)生地震破壞時(shí)的薄弱部位，以及影響其抗震性能的主要因素；Bocca[7]和Pande[8]分別提出了離散裂縫模型和分布裂縫模型，指出在彈性階段兩種模型與試驗(yàn)值均有不錯(cuò)的吻合度；Bennett等[9]研究了有構(gòu)造柱約束的墻體的抗震性能，認(rèn)為當(dāng)墻體有約束構(gòu)件時(shí)，其開(kāi)裂會(huì)延緩，而承載力會(huì)有所提高。

本文采用文獻(xiàn)[1]的本構(gòu)關(guān)系對(duì)一個(gè)三層磚混結(jié)構(gòu)用ABAQUS中的整體式模型研究結(jié)構(gòu)在地震作用下抗震性能，選用混凝土損傷塑性模型來(lái)建立砌體材料和混凝土材料的本構(gòu)關(guān)系。在定義強(qiáng)度時(shí)，參照混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范[10]中5.5.1條規(guī)定對(duì)材料強(qiáng)度統(tǒng)一取平均值；為了避免損傷因子過(guò)大而導(dǎo)致計(jì)算不收斂或收斂困難，設(shè)置材料最大損傷值為0.9。

1.2 材料本構(gòu)關(guān)系

本文的本構(gòu)關(guān)系參照文獻(xiàn)[1]選取，采用混凝土損傷塑性模型（CDP模型）描述混凝土的本構(gòu)關(guān)系，其單軸受拉（壓）應(yīng)力-應(yīng)變關(guān)系式依據(jù)混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范[10]C.2.4（C.2.3）條選取，鋼筋被認(rèn)為是一種理想的彈塑性材料。采用楊衛(wèi)忠[11]模型來(lái)模擬砌體材料的受壓本構(gòu)關(guān)系，強(qiáng)度指標(biāo)和基本力學(xué)參數(shù)按照《砌體結(jié)構(gòu)設(shè)計(jì)規(guī)范》（GB5003-2011）[12]選取，并按式（1）計(jì)算其受壓屈服應(yīng)變cmε。極限壓應(yīng)變?nèi)?0倍的屈服應(yīng)變[1]，極限拉應(yīng)變?nèi)?0倍的抗拉強(qiáng)度平均值對(duì)應(yīng)的應(yīng)變[1,6]。

2 磚混結(jié)構(gòu)模型模擬

2.1 磚混結(jié)構(gòu)模型參數(shù)

該模型共有三層，層高3.3m，墻厚240mm，預(yù)制樓板厚120mm，屋面板厚130mm，走廊處構(gòu)造柱的尺寸為240mm×240mm，箍筋為?8@200，未加密，除了走廊處構(gòu)造柱的縱筋采用4?16外，其余縱筋均采用?8；模型采用MU10燒結(jié)多孔磚，C25混凝土，M5水泥砂漿。圖1為模型平面圖，材料參數(shù)如表1所示。

圖1 計(jì)算模型平面圖Fig.1 Plan of calculation model

式中：fcm為抗壓強(qiáng)度平均值。

依據(jù)砌體規(guī)范，砌體受拉和受壓彈性模量取值相同，常用受壓彈性模量計(jì)算公式如下：

表1 模型材料參數(shù)

2.2 地震波的選取與調(diào)幅

本文計(jì)算分析中以EL CENTRO波作為地震動(dòng)輸入，地震動(dòng)持續(xù)時(shí)間為20s。加載雙向地震波時(shí)X方向（橫墻）最大幅值取表2中數(shù)值，Y方向（縱墻）幅值為X向0.85倍。

表2 地震加速度時(shí)程曲線(xiàn)的最大值（cm/s2）

2.3 模型的建立

首先在ANSYS中建模后再寫(xiě)成CDB文件導(dǎo)入到ABAQUS中，其中砌體和混凝土都采用C3D8R單元模擬，鋼筋采用T3D2單元，用Truss單元模擬鋼筋并嵌入(Embed)到構(gòu)造柱、樓板和屋面板實(shí)體中。

2.4 模態(tài)分析

本文模態(tài)分析選用ABAQUS中的隱式求解功能中的線(xiàn)性攝動(dòng)頻率分析，用Lanczos特征求解器和單精度計(jì)算方法提取了如表3所示的前八階自振頻率和相應(yīng)的自振周期。

表3 模型前8階自振頻率及對(duì)應(yīng)的固有周期（s）

3 模型地震反應(yīng)對(duì)比

墻體剪切破壞主要表現(xiàn)為在墻體產(chǎn)生斜裂縫或者交叉斜裂縫，其實(shí)質(zhì)在于破壞點(diǎn)的主拉應(yīng)力或主壓應(yīng)力超過(guò)材料的抗拉或抗壓強(qiáng)度極限值。由于砌體屬于脆性材料，易發(fā)生受拉損傷，結(jié)構(gòu)破壞以受拉破壞為主，而受壓破壞相對(duì)較少，故本文用受拉損傷程度來(lái)表征結(jié)構(gòu)的破壞程度[13]。在對(duì)砌體結(jié)構(gòu)進(jìn)行整體式建模時(shí)，由于沒(méi)有引入有關(guān)裂縫的模擬參數(shù)，因此不會(huì)出現(xiàn)真正的裂縫，但受拉損傷參數(shù)云圖（DAMAGET）可以反映塑性應(yīng)變的大小，通過(guò)對(duì)DAMAGET中“模擬裂縫”出現(xiàn)的形狀和位置進(jìn)行觀(guān)察來(lái)近似模擬結(jié)構(gòu)的開(kāi)裂程度。在受拉損傷參數(shù)云圖中，顏色從藍(lán)到紅體現(xiàn)出受拉損傷破壞程度逐步加大的分布特點(diǎn)。

本文分析了九個(gè)工況，分別模擬結(jié)構(gòu)在6、7、8度小震、中震以及大震作用下裂縫的開(kāi)展情況，計(jì)算結(jié)果如圖2所示。

從計(jì)算結(jié)果可知，結(jié)構(gòu)在6、7、8度小震下完好，故本文重點(diǎn)是分析結(jié)構(gòu)在6、7、8度中震和大震下的破壞狀況。經(jīng)過(guò)系統(tǒng)總結(jié)、分析后可以得出以下規(guī)律性結(jié)果：隨著地震動(dòng)強(qiáng)度的增加，結(jié)構(gòu)破壞越來(lái)越嚴(yán)重，且破壞區(qū)域逐步加大。破壞區(qū)域主要集中在：（1）凸出墻面的轉(zhuǎn)角部位或墻體截面變化部位較易形成裂縫；（2）房屋墻角出現(xiàn)裂縫破壞；（3）在門(mén)、窗等洞口處，墻面出現(xiàn)斜裂縫或X形交叉裂縫，這與實(shí)際震害現(xiàn)象是相符的。其中一層破壞最為嚴(yán)重，二層、三層依次減弱。強(qiáng)烈地震作用時(shí)可導(dǎo)致墻體發(fā)生水平剪切破壞，在一層墻體底部甚至出現(xiàn)水平裂縫。模擬結(jié)果體現(xiàn)出墻體一層B軸、C軸破壞較為嚴(yán)重，另外由于外墻未設(shè)置構(gòu)造柱，所以①軸、④軸墻體損壞也較為嚴(yán)重。

圖2 EL Centro波作用下結(jié)構(gòu)受拉損傷云圖Fig.2 Structural tensile damage cloud under EL centro seismic wave

圖3 墻體受拉損傷參數(shù)Fig.3 Parameter of wall tension damage

3.1 模型的受拉損傷參數(shù)對(duì)比

本項(xiàng)研究對(duì)墻體破壞較為嚴(yán)重的B軸、C軸以及①軸、④軸的受拉損傷參數(shù)進(jìn)行了較為詳細(xì)的分析。限于篇幅，本文僅列出B軸和①軸的受拉損傷參數(shù)平均值，用平均值對(duì)圖3墻體受拉損傷云圖進(jìn)行量化。經(jīng)分析得出受拉損傷破壞均是在加載地震波0.6s后產(chǎn)生，并且在0.6～3s急劇發(fā)展，達(dá)到峰值后基本保持不變，在12s后達(dá)到破壞狀態(tài)。從圖3可以看出，8度大震作用下各層墻體的受拉損傷平均值最大，破壞程度可依次劃分為：8度大震＞7度大震＞8度中震＞ 6度大震＞ 7度中震＞ 6度中震，一層＞二層＞三層，B軸＞ C軸，另外由于未設(shè)置構(gòu)造柱，所以①軸、④軸的破壞也很?chē)?yán)重。

3.2 層間位移角

層間位移角是衡量結(jié)構(gòu)破壞程度的一個(gè)重要指標(biāo)，表4為文獻(xiàn)[5]劃分的層間位移角破壞限值。據(jù)此本文以層間位移角而劃分的結(jié)構(gòu)破壞狀況如表5所示，從中可以看出最大層間位移角都發(fā)生在第一層的X方向。將各工況下模型的最大層間位移角與界限值進(jìn)行對(duì)比分析，可以得到：6度、7度中震和6度大震作用下結(jié)構(gòu)處于完好狀態(tài)，7度大震和8度中震作用下結(jié)構(gòu)從輕微向中等破壞發(fā)展；8度大震作用下結(jié)構(gòu)達(dá)到中等破壞且接近嚴(yán)重破壞狀態(tài)。

3.3 受拉損傷面積百分比

為了判斷砌體結(jié)構(gòu)在地震作用下的破壞程度，在前兩節(jié)的基礎(chǔ)上參考文獻(xiàn)[3]增加一個(gè)墻片的受拉損傷面積百分比作為衡量模型破壞程度的指標(biāo)，其定義為在地震動(dòng)加載作用結(jié)束時(shí)墻片中各個(gè)單元受拉損傷參數(shù)值大于0的單元數(shù)與整個(gè)墻片單元總數(shù)的比值。

表6為本文計(jì)算的各層墻片的受拉損傷面積百分比，其中一層的受拉損傷面積比最大，二層、三層依次減少，由此得出一層的破壞遠(yuǎn)遠(yuǎn)大于二層、三層的破壞。

表4 層間位移角限值

表5 層間位移角

表6 受拉面積百分比

文獻(xiàn)[14]中給出的結(jié)構(gòu)破壞程度界限值為：當(dāng)受拉損傷面積百分比小于5%時(shí)為基本完好；大于5%而小于15%時(shí)為輕微破壞；大于15%而小于30%時(shí)為中等破壞；大于30%而小于50%時(shí)為嚴(yán)重破壞；大于50%時(shí)為倒塌破壞。據(jù)此劃分的結(jié)果為：在6度、7度中震和6度大震作用下結(jié)構(gòu)處于基本完好狀態(tài)；在7度大震和8度中震作用下結(jié)構(gòu)處于中等破壞狀態(tài)；在8度大震作用下結(jié)構(gòu)處于嚴(yán)重破壞甚至倒塌狀態(tài)。

4 結(jié)論

本文運(yùn)用ABAQUS對(duì)磚混結(jié)構(gòu)有限元模型進(jìn)行非線(xiàn)性數(shù)值模擬分析，選取受拉損傷參數(shù)平均值、最大層間位移角以及受拉損傷面積百分比作為衡量指標(biāo)，從多角度、定量的方式對(duì)磚混結(jié)構(gòu)在各等級(jí)地震動(dòng)作用下的破壞程度進(jìn)行綜合評(píng)價(jià)，得出如下結(jié)論：

（1）加載各等級(jí)地震作用后，通過(guò)對(duì)受拉損傷云圖進(jìn)行分析，可以得到磚混結(jié)構(gòu)的薄弱部位為門(mén)窗洞口處、一層的最底層墻體、墻體截面變化處、房屋的墻角處以及墻角凸出墻面的轉(zhuǎn)角處。

（2）通過(guò)對(duì)最大層間位移角和受拉損傷面積百分比與相應(yīng)界限值進(jìn)行對(duì)比分析，兩者得出的結(jié)論基本一致，即在6、7、8度小震以及6、7度中震作用和6度大震作用下，結(jié)構(gòu)處于基本完好狀態(tài)；在8度中震、7度大震作用下，結(jié)構(gòu)處于中等破壞狀態(tài)；在8度大震作用下，結(jié)構(gòu)處于中等破壞但接近嚴(yán)重破壞，個(gè)別位置處于倒塌狀態(tài)。

（3）通過(guò)數(shù)值分析確定了磚混結(jié)構(gòu)易發(fā)生破壞的薄弱部位、破壞程度、定量的性能指標(biāo)和參數(shù)，建立了與各種破壞狀態(tài)相對(duì)應(yīng)的定量參數(shù)，可以有效的將以往主要依據(jù)宏觀(guān)破壞現(xiàn)象定性評(píng)定結(jié)構(gòu)破壞程度的方法與本文定量的評(píng)價(jià)方法相結(jié)合，為合理評(píng)價(jià)震后磚混結(jié)構(gòu)破壞程度提供科學(xué)的理論依據(jù)。

[1] 鄭妮娜. 裝配式構(gòu)造柱約束砌體結(jié)構(gòu)抗震性能研究[D].重慶:重慶大學(xué), 2010.

[2] 殷園園. 芯柱式構(gòu)造柱約束砌體結(jié)構(gòu)抗震性能評(píng)價(jià)[D].重慶:重慶大, 2011.

[3] 劉露. 砌體結(jié)構(gòu)強(qiáng)震破壞模式研究[D]. 重慶:重慶大學(xué), 2012.

[4] 唐春美. 砌體結(jié)構(gòu)的非線(xiàn)性性能模擬及大震不倒對(duì)策研究[D]. 重慶:重慶大學(xué), 2013.

[5] 牟道燦. 砌體結(jié)構(gòu)抗地震倒塌能力分析[D]. 重慶:重慶大學(xué), 2012.

[6] 劉杰. 基于ABAQUS整體式模型下砌體結(jié)構(gòu)抗震性能影響因素研究[D]. 長(zhǎng)沙:湖南大學(xué), 2014.

[7] Bocca P, Carpinteri A, Valente S. Fracture mechanics of brick masonry: size effects and snap-back analysis[J]. Materials and Structures, 1989，22（5）: 364-373.

[8] Pande G N, Middleton J, Lee J S, et al. Numerical simulation of cracking and collapse of masonry panels subject to lateral loading[C]//Proc. 10th International Brick and Block Masonry Conference, 1994.

[9] Banting B R, El-Dakhakhni W W. Force-and displacementbased seismic performance parameters for reinforced masonry structural walls with boundary elements[J]. Journal of Structural Engineering, 2012，138（12）: 1477-1491.

[10] GB 50010-2010混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范[S]. 北京:中國(guó)建筑工業(yè)出版社, 2010.

[11] 楊衛(wèi)忠. 砌體受壓本構(gòu)關(guān)系模型[J]. 建筑結(jié)構(gòu), 2008, 38（10）：80-82.

[12] GB 50003-2011 砌體結(jié)構(gòu)設(shè)計(jì)規(guī)范[S]. 北京:中國(guó)建筑工業(yè)出版社.

[13] 楊玉成, 楊柳, 高云學(xué), 等. 多層磚房的地震破壞和抗裂抗倒設(shè)計(jì)[M]. 北京:地震出版社, 1981.

[14] 李佳. 基于數(shù)值模擬的砌體結(jié)構(gòu)倒塌影響因素分析及抗倒塌影響因素分析及抗倒塌措施初探[D]. 重慶:重慶大學(xué), 2013.

The Evaluation of Seismic Safety for Brick-Concrete Structure

KONG Fang-fang, WANG Guo-xin
（State Key Laboratory of Coastal and Offshore Engineering， Faculty of Infrastructure Engineering， Dalian University of Technology， Liaoning Dalian 116024，China）

The brick-concrete structures were badly damaged usually in the previous earthquakes due to their poor ductility and small deformation capacity. So the quantitative research on seismic performance of such structures based on traditional qualitative one is more important for designers to improve the structural seismic capability. The elastic and plastic analysis of one three-story brick concrete structure under different earthquake intensity is carried out by ABAQUS. By comparing and analyzing the tensile damage parameter values and tensile damage clouds of the structure， it is known that the maximum story displacement angle and the percentage of damage area can be chosen as proper measurement index to represent destructive degree of such structure， and also such indexes are used to describe structural different damage statuses under earthquakes.

brick concrete structure; time history analysis; damage

P315.9

A

10.13693/j.cnki.cn21-1573.2017.02.009

1674-8565（2017）02-0052-07

國(guó)家自然科學(xué)基金項(xiàng)目（51378092/51578113）

2017-01-22

2017-03-05

孔芳芳（1991-），女，河南省商丘市人，大連理工大學(xué)防震減災(zāi)工程及防護(hù)工程專(zhuān)業(yè)，碩士研究生，現(xiàn)主要從事數(shù)值模擬、震后安全性評(píng)價(jià)方面的工作。 E-mail：1538004519@qq.com