張昊宇，陳洪富，林旭川，陳 浩，毛晨曦

(中國地震局地震工程與工程振動重點實驗室(中國地震局工程力學研究所)，哈爾濱150080)

?

鋼筋混凝土框架柱端剪切破壞震害分析及對策

張昊宇，陳洪富，林旭川，陳 浩，毛晨曦

(中國地震局地震工程與工程振動重點實驗室(中國地震局工程力學研究所)，哈爾濱150080)

為避免鋼筋混凝土框架柱在地震中因填充墻附加剪力導致剪切破壞，本文總結了近年來國內(nèi)外破壞性地震中出現(xiàn)的典型柱端剪切破壞震害，采用ABAQUS對典型填充墻框架進行了水平推覆模擬，據(jù)此分析了框架柱端剪切破壞的受力機理，并以魯?shù)榈卣鹬械囊粭潎乐卣饟p的實際結構為例，采用ABAQUS對典型填充墻框架進行了水平推覆模擬和受剪承載力計算.分析表明，填充墻附加剪力是造成框架柱端剪切破壞的主要原因.依據(jù)分析結果和已有試驗成果給出了考慮填充墻附加剪力的柱端剪力荷載計算公式.震損結構案例分析表明，忽略填充墻附加剪力可能導致結構的抗倒塌能力嚴重下降，本文填充墻附加剪力公式計算結果與有限元分析結果及實際震害符合較好，可供工程設計參考.關鍵詞: 震害調(diào)查; 鋼筋混凝土框架; 框架柱; 剪切破壞; 填充墻

強柱弱梁、強剪弱彎、強節(jié)點弱構件是鋼筋混凝土(RC)框架結構的基本抗震設計原則，可確保結構在強震中具有較好的延性和抗倒塌能力.但近年來國內(nèi)外多次破壞性地震中很多RC框架結構出現(xiàn)了較為嚴重的薄弱層柱頂剪切破壞,嚴重威脅結構的“大震不倒”能力.本文作者歷經(jīng)2013年蘆山地震、2014年魯?shù)榈卣鸷?015年尼泊爾地震震后調(diào)查，總結了這3次地震中的RC框架結構柱剪切破壞，通過ABAQUS典型填充墻框架水平推覆分析，對該問題進行探討并初步給出防治措施.通過對魯?shù)榈卣鹬幸粭澱饟p結構的非線性分析，對柱端剪切破壞危害和防治措施有效性進行討論.

1 鋼筋混凝土框架柱的典型震害

圖1為2013年蘆山7.0級地震[1]及2015年尼泊爾8.1級地震中典型的RC框架柱剪切破壞；圖2為2014年云南魯?shù)榈卣鹬旋堫^山鎮(zhèn)典型RC框架的破壞[2].

上述RC框架結構震害呈現(xiàn)以下特點：1)絕大多數(shù)框架表現(xiàn)為強梁弱柱破壞；2)框架柱破壞模式與緊鄰填充墻相關：柱兩側無填充墻或毗鄰填充墻有較大洞口時，柱多為壓彎破壞；小開洞，無開洞填充墻，特別是較長的無洞口填充墻，相鄰框架柱多呈剪切破壞，且破壞程度明顯重于壓彎破壞柱.需要注意的是，與剪切破壞柱相鄰的填充墻，既有粘土磚墻(圖1(a)、(c))也有空心磚墻(圖1(b)).

剪切破壞顯著降低框架柱的延性，嚴重威脅結構的抗倒塌能力，應盡量避免.中國現(xiàn)行抗震規(guī)范[3]并未針對填充墻給框架柱的剪力影響給出明確的計算方法.89版抗震規(guī)范[4]曾在附錄三中提及填充墻附加剪力，但僅針對粘土磚填充墻框架，且計算方法有待細化.國內(nèi)外學者對填充墻RC框架開展了大量研究，但大部分側重于填充墻等效支撐模擬[5]、填充墻對整體結構或框架的影響[6]、有限元模擬方法[7-8]等.文獻[9]對填充墻導致框架柱的破壞進行了闡述，文獻[10]給出了填充墻框架整體的抗剪承載力計算方法，但如何從工程設計角度避免此類剪切破壞，仍有待討論.本文對此開展了研究.

圖1 蘆山地震及尼泊爾地震中典型的RC框架柱剪切破壞[1]

2 框架柱剪切破壞機理分析

為探究柱端剪切破壞機理，對圖3的6個鋼筋混凝土框架建立有限元模型并進行水平推覆分析.

2.1 模型驗證及建模計算

圖3框架的混凝土強度設計等級C30，縱筋HRB335，箍筋HPB235，填充墻為Mu10實心磚，M5砂漿砌筑，厚度200 mm.材料強度均取用平均值.采用ABAQUS建模分析，混凝土及砌體均采用混凝土塑性損傷模型，混凝土應力應變關系選取文獻[11]偏壓偏拉本構(圖4(a))；砌體抗壓強度按文獻[12]附錄D計算，采用文獻[13]應力應變關系(圖4(b)).鋼筋選取三折線模型(圖4(c))，embed(嵌入)到構件中.采用無摩擦接觸模擬填充墻與框架間的相互作用.

對文獻[7]空框架、空心磚及實心磚填充墻框架試件(圖5(a))進行水平推覆模擬(空心砌塊抗壓強度按砌塊橫截面孔洞率進行折減).由圖5(b)可知，盡管在力-位移曲線的下降段存在一定誤差，但模擬承載力與試驗大體相符，可用于填充墻框架承載力分析.

對6個框架分別進行水平推覆分析：首先按設計軸壓比0.4在柱頂施加豎向力，之后在梁端施加水平位移，直至框架層間側移角達1/50或計算不再收斂.

圖3 計算框架(mm)

圖4 材料本構關系曲線

圖5 對文獻[7]試驗的模擬驗證

2.2 結果分析

圖3(a)框架加載至極限位移時的misses云圖如圖6(a)：填充墻所受剪力大部分傳至左側框架柱頂端，此即為填充墻附加剪力.圖6(b)左側柱頂箍筋的應變明顯大于右側，且已接近屈服應變.填充墻中剪力的傳遞路徑如圖6(c)：其一由右下角至左上角，其二由底邊均勻傳至左側柱體中部.當柱抗剪承載力較小時(圖3(d))，破壞模式變?yōu)閴w僅局部壓陷，柱頂剪切破壞.這與圖1、圖2震害現(xiàn)象吻合.圖7給出了框架a、b、c的總剪力及各框架柱剪力與層間位移的關系曲線.可見填充墻框架承受的水平剪力大于空框架，強柱框架a的層間剪力大于弱柱框架c.這是由于當柱有較大抗剪承載力時，墻體發(fā)生剪切破壞，其抗剪承載力充分發(fā)揮；當柱的抗剪承載力較弱時，墻體的抗剪承載力大于柱，造成柱端剪切破壞，墻體抗剪承載力未充分發(fā)揮.

圖8為框架a、b在1/120層間側移角時的柱剪力圖：框架a左側柱頂剪力約為空框架柱的2.8倍，可知在設計階段不考慮填充墻附加剪力會明顯低估框架柱的剪力荷載；圖9為框架a不同側移時填充墻作用于左柱的線荷載.最大荷載位置隨側移角的增大明顯下移，緣于墻角局部壓碎，造成墻體對柱的水平力合力點下移，在框架柱端形成短柱；圖10為不同層間側移角時框架a的柱剪力分布，可見絕大部分填充墻附加剪力在結構彈性階段(1/550)已施加，因此無論小震設計或大震驗算都應考慮填充墻附加剪力.

圖6 框架a、b、c模擬結果分析

圖7 框架a、b、c加載點力-位移曲線

圖8 框架a、b層間側移角1/120時柱剪力對比(kN)

圖9 不同側移角下填充墻作用于柱的線荷載(框架a)

圖10 框架a柱剪力圖(kN)

對圖3(d)、(e)、(f)模型進行水平推覆分析，可得圖11結果.圖11(a)窗臺上部左邊墻體以受彎為主，提供的附加剪力較??；圖11(b)框架的附加剪力主要為門口之間墻體提供；圖11(c)填充墻破壞模式與圖3框架a類似.圖11(d)表明附加剪力大小受墻體“有效水平截面”面積的影響.

3 填充墻附加剪力的計算建議

框架柱和填充墻的破壞互相制約：當墻體抗剪承載力大于框架柱抗剪承載力時，框架柱首先破壞，之后喪失承載力并卸載，墻體破壞即被推遲(圖6(d))；當墻體抗剪承載力小于框架柱時，柱子更可能發(fā)生壓彎破壞，而墻體發(fā)生整體剪切破壞(圖6(a)).相比之下，后者實現(xiàn)了強剪弱彎，保證了框架柱的延性，確保填充墻作為第二道抗震防線，為主體結構提供支撐和耗能，是希望出現(xiàn)的破壞模式.為實現(xiàn)該受力模式，宜在工程設計中定量考慮填充墻附加剪力.參考規(guī)范對底框結構嵌砌于框架之間的填充磚墻框架抗震受剪承載力驗算公式(7.2.9-3)[12]，建議毗鄰填充墻的框架柱按式(1)計算剪力：

(1)

(2)

式中:Vc為框架柱自身受彎產(chǎn)生的剪力，可按抗震規(guī)范[4]6.2.5、6.2.6條計算；Vf為填充墻附加剪力；fvE為填充墻砌體沿階梯形截面破壞的抗震抗剪強度設計值；Aw0為填充墻有效附加剪力墻段水平截面計算面積，無洞口時取實際截面的1.6倍(文獻[10]給出Aw0在無洞口時宜取1.25倍實際截面面積，但參考其他文獻和本文分析結果，無洞口時取實際截面1.6倍更安全)，有洞口時取截面凈面積，但不計入寬度小于洞口高度1/4的墻肢截面面積.

圖11 框架d、e、f內(nèi)力分析

偏于保守地認為附加剪力完全傳遞至一個框架柱上.當一根柱兩側均有填充墻時，附加剪力可取柱兩側填充墻附加剪力中的較大值.建議柱頂及柱底均按照該剪力進行抗剪承載力驗算.

為檢驗式(2)，采用該公式計算了文獻[7]、[9]、[14]試件及圖3框架 a、d、e、f的柱端剪力Vf，并將Vf與實際柱端剪力(試驗結果或有限元分析結果)之比繪于圖12.比值越接近1說明式(2)越合理.可見式(2)計算值基本符合試驗結果和分析結果，且大于大部分試驗及分析結果，略偏于安全.

圖12 附加剪力試驗值和有限元計算結果與式(2)計算值對比

Fig.12 Additional shear forces obtained by equation (2) vs. test & FEM analysis results

4 附加剪力計算公式的驗證

為驗證式(1)、(2)的有效性和填充墻附加剪力對整體結構抗震性能的影響，本節(jié)分別以圖3(a)、(c)框架，以及圖2龍頭山幼兒園為例進行對比分析.

4.1 單榀框架

分別計算圖3(a)、(c)框架單個柱子的抗剪承載力設計值(偏于保守地取剪跨比λ=1計算抗剪承載力)和剪力荷載設計值(按式(1)計算，其中Vc按抗規(guī)[3]式6.2.5-2，一級框架計算)，二者比值及填充墻附加剪力在總剪力荷載值中所占比重見表1.

表1 單榀框架抗剪承載力計算結果

對比表1及圖6(a)、(d)可知，填充墻附加剪力在柱端剪力荷載中占比較大(大于50%).框架a的柱子通過了抗剪承載力驗算，有限元模擬中實現(xiàn)了“強柱弱填充墻”破壞；框架c柱子抗剪承載力無法通過驗算，在模擬中則發(fā)生“弱柱強填充墻”破壞.因此公式(1)、(2)較好地預測了填充墻框架受力情況.

4.2 龍頭山鎮(zhèn)中心幼兒園

圖2為一棟2011年建造的3層RC框架.混凝土強度設計等級C25，填充墻為180 mm厚的粘土磚墻，在2013年魯?shù)榈卣鹬性庥鰪娬穑讓佣鄠€柱頂嚴重剪切破壞，震后結構底層明顯側移，震損程度達毀壞.距該結構不足100 m的龍頭山強震臺站(圖2(a))記錄了地面加速度時程，東西向水平加速度峰值達949 gal.選取該結構設置3個力學模型進行模擬分析：

1)僅考慮填充墻重力荷載，不考慮填充墻和框架間相互作用的原型結構空框架；

2)考慮填充墻與框架相互影響(即為原型結構)；

3)在2基礎上，所有框架柱均采用圖2(d)右下角加強型配箍(加強后框架柱均滿足抗剪承載力要求).

采用ABAQUS進行分析.圖13為原型結構有限元模型.框架梁柱及填充墻采用實體單元C3D8R，樓板采用殼單元，框架柱及填充墻使用混凝土塑性損傷模型，樓蓋無明顯震損，近似采用彈性本構；框架柱的鋼筋采用桁架單元，embed(嵌入)到柱中.材料強度依據(jù)規(guī)范[12，15]取用平均值，材料應力應變關系與前述一致.按照圖14步驟分別對3個模型進行分析.分析過程分為兩個階段：1)非線性時程分析(采用顯式分析以獲得較好的收斂性)，得到結構各方向的位移反應；2)以得到的最大位移為條件，采用隱式算法對結構一層進行靜力推覆，得到較為精確的柱端箍筋應變(顯式分析的鋼筋應變結果精度較差).

圖13 原型結構有限元分析模型

圖14 結構分析步驟示意

先進行時程分析.加速度時程選取圖2中龍頭山臺站的3向地震動記錄(圖15)，截取0～10 s加速度記錄作為輸入.圖16為結構各樓層平均位移(同一層各點位移平均值)時程、最大層間側移、最大層間側移角及一層頂不同位置的層間位移時程.可知結構底層為薄弱層，最大層間側移角為x向1/29，y向1/20，且x、y向均出現(xiàn)明顯殘余變形(x向平均28 mm；y向平均39 mm)，這與圖2震害現(xiàn)象大體相符.

圖15 龍頭山臺站加速度時程記錄及5%阻尼比反應譜

之后采用隱式分析方法，按照最大位移(x向125 mm；y向180 mm)給一層結構分別加載 (圖16d曲線表明，結構在最大位移反應區(qū)段主要表現(xiàn)為x、y向平動，扭轉分量較小，因此沿兩個方向?qū)Y構進行靜力推覆，不考慮結構扭轉是合理的)，之后提取柱頂外圍箍筋的最大應變平均值，繪成柱狀圖，見圖17(由于結構左右對稱，故僅選取了左半部分結構).

圖16 原型結構時程分析位移結果

圖17 柱頂剪力荷載、抗剪承載力計算值及外圍箍筋平均應變對比

按照4.1節(jié)方法計算每個柱的抗剪承載力設計值和剪力荷載設計值，將二者比值及填充墻附加剪力在總剪力荷載值中所占比重繪于圖17.由此可知：

1)荷載及抗力計算表明，部分毗鄰填充墻的框架柱，填充墻附加剪力在總剪力荷載中占比很大，最大可達68%；圖17的 21根框架柱中，10根柱不滿足抗剪承載力要求，而實際震害中共有6根框架柱發(fā)生明顯剪切破壞(圖17中已用虛線框出)，全部在這10根柱子范圍內(nèi)，表明本文剪力計算公式符合實際.

2)靜力推覆分析顯示，不考慮填充墻的影響，柱端箍筋應變明顯降低，柱端剪力荷載被顯著低估；考慮填充墻影響后，柱端箍筋應變明顯增大，部分應變接近或超過箍筋屈服應變，且與荷載及抗力計算結果基本一致，尤其①、②、⑥軸與C、D軸相交的6根柱子，其箍筋應變、荷載抗力分析與圖2實際震害現(xiàn)象基本吻合；③軸與E軸相交處的柱頂實際發(fā)生了剪切破壞，而分析顯示箍筋并未屈服，這可能是填充墻材料強度的離散性及有限元模型精度不足所致；加強配箍后，所有柱端箍筋應變均低于屈服應變，說明適當增加配箍是有效的防治措施.

3)抗力驗算及有限元分析表明，未考慮填充墻附加剪力，該結構48%的框架柱受剪切破壞的威脅；震害調(diào)查則顯示，該結構29%的框架柱在地震中出現(xiàn)了明顯剪切破壞.框架柱一旦發(fā)生嚴重剪切破壞，幾乎完全喪失水平及豎向承載力，因此忽略填充墻附加剪力可能會明顯降低RC框架的抗倒塌能力.

5 結 論

近年來多次地震震害調(diào)查顯示，毗鄰填充墻的RC框架柱剪切破壞是RC框架結構的一種典型震害，往往在結構薄弱層的小開洞、無開洞填充墻相鄰的框架柱頂端出現(xiàn)，且破壞程度明顯重于同層的壓彎破壞柱.本文對這一破壞現(xiàn)象進行分析并提出了避免此類破壞的設計方法.主要結論如下：

1)基于ABAQUS建模的典型填充墻框架水平推覆分析表明，框架發(fā)生水平位移時，填充墻的附加剪力幾乎全部傳遞至框架一側的柱頂端，該剪力在水平位移初期即達較大值，且合力作用點位于梁柱節(jié)點下方.基于這些分析，參考底框結構相關設計方法，給出了填充墻附加剪力估算公式，供工程設計參考.

2)魯?shù)榈卣鹬旋堫^山鎮(zhèn)幼兒園RC框架有限元分析及荷載與抗力計算表明，該結構底層近半數(shù)框架柱可能發(fā)生剪切破壞，實際震害中29%的框架柱發(fā)生嚴重剪切破壞，嚴重威脅結構的抗倒塌能力；本文填充墻附加剪力公式的計算結果與有限元分析結果及實際震害均符合較好；該公式能夠較為準確地估計框架柱端剪力，從而為調(diào)整柱端配箍，避免框架柱剪切破壞提供依據(jù).

[1] 孫柏濤, 閆培雷, 王明振,等. 四川省蘆山“4·20”7.0級強烈地震建筑物震害圖集[M]. 北京: 地震出版社, 2014. SUN Baitao, YAN Peilei, WANG Mingzhen, et al. Atlas of Sichuan Province Lushan ‘4.20’ 7.0 violent earthquake damage [M]. Beijing: Seismological Press, 2014.

[2] LIN Xuchuan, ZHANG Haoyu, CHEN Hongfu, et al. Field investigation on severely damaged aseismic buildings in 2014 Ludian earthquake [J]. Earthquake Engineering and Engineering Vibration, 2015, 14 (1): 169-176.

[3] 建筑抗震設計規(guī)范:GB50011—2010[S]. 北京: 中國建筑工業(yè)出版社, 2010. Code for seismic design ofbuildings:GB50011—2010[S]. Beijing: China building industry press, 2010.

[4] 建筑抗震設計規(guī)范:GBJ11—89[S]. 北京: 中國建筑工業(yè)出版社，1989. Code for seismic design of buildings:GBJ11—89[S]. Beijing: China building industry press, 1989.

[5] ASTERIS P G. Lateral stiffness of brick masonry infilled plane frames [J]. Journal of structural Engineering, 2003, 129: 1071-1079. [6] ANIL O, ALTIN S. An experimental study on reinforced concrete partially infilled frames [J]. Engineering Structures, 2007, 29: 449-460.

[7] MEHRABI A B, SHING P B. Finite Element modeling of masonry-infilled RC frames [J]. Journal of structural Engineering, 1997(3): 604-613.

[8] 楊偉, 歐進萍. 抗震結構中填充墻倒塌數(shù)值模擬與分析 [J]. 哈爾濱工業(yè)大學學報, 2011, 43(6): 16-19. YANG Wei, OU Jinping. Collapse simulation and analysis of infill walls for earthquake-resistant structures [J]. Journal of Harbin Institute of Technology, 2011, 43(6): 16-19.

[9] 黃群賢, 郭子雄, 朱雁茹 等. 混凝土空心砌塊填充墻RC 框架抗震性能試驗研究 [J]. 建筑結構學報, 2012, 33 (2): 110-118. HUANG Qunxian, GUO Zixiong, ZHU Yanru, et al. Experimental study on seismic behavior of RC frames infilled with concrete hollow blocks [J]. Journal of Building Structures, 2012, 33(2):110-118.

[10]童岳生, 錢國芳. 磚填充墻鋼筋混凝土框架的變形性能及承載能力[J]. 西安冶金建筑學院學報. 1985,42(2): 1-21. TONG Yuesheng, QIAN Guofang. Deformation and bearing capacity of reinforced concrete frame infilled with bricks [J]. Journal of Xi'an Institute of metallurgy and construction. 1985,42(2): 1-21.

[11]過鎮(zhèn)海. 混凝土的強度和變形-試驗基礎和本構關系[M]. 北京: 清華大學出版社, 1997. GUO Zhenhai. Strength and deformation of concrete-experiment and constitutive relations [M]. Beijing: Tsinghai university press, 1997.

[12]砌體結構設計規(guī)范: GB50003—2001[S]. 北京: 中國建筑工業(yè)出版社, 2002. Code for seismic design of masonry structures:GB50003—2001 [S]. Beijing: China building industry press, 2002.

[13]劉桂秋. 砌體結構基本受力性能的研究[D]. 長沙：湖南大學, 2005. LIU Guiqiu. The reseatch on the basic mechanical behavior of masonry structure [D]. Changsha:Doctoral Dissertation of Hunan University, 2005 .

[14]DEMETRIOS J K, CHRISTOS G K. Experimental Investigation of Infilled Reinforced Concrete Frames with Openings [J]. ACI Structural Journal, 2009, 106(2): 132-141.

[15]混凝土結構設計規(guī)范：GB50010—2010[S]. 北京: 中國建筑工業(yè)出版社, 2010. Code for design of concrete structures：GB50010—2010[S]. Beijing: China building industry press, 2010.

(編輯 趙麗瑩)

Mechanismanalysis and prevention of seismic shear failure of RC frame columns

ZHANG Haoyu, CHEN Hongfu, LIN Xuchuan, CHEN Hao, MAO Chenxi

(Key Laboratory of Earthquake Engineering and Engineering Vibration (Institute of Engineering Mechanics,China Earthquake Administration), Harbin 150080, China)

To prevent reinforced concrete (RC) frame columns from seismic shear failure caused by additional shear force of the infill wall, severe shear failures of RC frame columns investigated in recent earthquakes were presented. Finite element modeling (FEM) pushover analysis with ABAQUS and mechanical characteristics analysis on 6 infilled RC frames were conducted. Then a RC frame structure subjected to LuDian Ms6.5 earthquake was analyzed by ABAQUS. The analysis indicates that the additional shear force of the infill wall results in the shear failure of the column. To ensure that a structure operates in a high-ductility mode, an estimation formula was developed for the additional shear force which could be used in structure design. The case study shows that without appropriate considerations about the additional shear force of the infill wall, anti-collapse capability of the structure will be threatened seriously. Analysis also demonstrates that the additional shear force estimation formula is in well accordance with the FEM analysis result and the seismic damage investigations, it can be a reference for engineering design.

earthquake damage investigation; RC frame; frame column; shear failure; infill wall

10.11918/j.issn.0367-6234.2016.12.005

2015-11-29

地震行業(yè)科研專項經(jīng)費項目(201308003);國家自然科學基金青年基金(51108433，51308511，51408562)

張昊宇(1978—), 男, 博士，副研究員

毛晨曦，maochenxi2006@163.com

TU375.3

A

0367-6234(2016)12-0043-08