Study on Mechanical Properties of Corroded Reinforced Concrete Frame Columns Under Cyclic Load
-
摘要: 采用有限元软件ABAQUS,以锈蚀率(0%、5%、10%、15%和20%)为变量,对5根钢筋混凝土柱的力学性能进行了数值模拟,研究各试件的滞回性能、骨架曲线、延性及耗能能力,分析钢筋锈蚀率对承载力、延性、耗能和塑性铰转动能力的影响。研究结果表明:模拟分析得到的锈蚀钢筋混凝土柱的强度和变形与试验结果吻合较好,建立的有限元模型可用于锈蚀钢筋混凝土柱的力学性能分析;混凝土开裂前,锈蚀构件的力学性能基本与未锈蚀构件相同,混凝土开裂后,构件的承载力、屈服荷载、极限位移、延性等均随钢筋锈蚀率的增大而降低;轻度锈蚀构件的滞回性能和破坏形式与未锈蚀构件类似,随着钢筋锈蚀率逐渐增大,滞回环的饱满程度降低,“捏拢”现象严重,滞回曲线由“弓形”逐渐发展成“反S形”,耗能能力降低,破坏形式趋于脆性破坏,位移延性系数、平均耗能系数等指标逐渐下降。Abstract: The numerical simulation of the mechanical properties of five corroded reinforced concrete columns with different degree of corrosion was presented with the finite element software ABAQUS. Hysteretic behavior, skeleton curve, ductility and energy dissipation capacity of specimens was studied. The influence of corrosion amount on bearing capacity, ductility, energy dissipation capacity and plastic hinge rotation capability were analyzed. Our results show that the strength and deformation of the specimens obtained by simulation are in good agreement with the test results. The numerical model established in this paper can be used for the mechanical properties analysis of corroded reinforced concrete columns. Before the cracking of concrete, the corrosion of steel bars has no obviously effect on the mechanical properties of concrete bending members. After the cracking of concrete, bearing capacity, yield load, ultimate displacement and ductility of specimens decrease with the increase of the corrosion rate. The hysteresis behavior and the failure mode of mildly corroded specimens are similar to those of non-corroded specimens. As the increase of corrosion amount, the hysteresis loop gradually deflates, and its pinch phenomenon is more serious. The shape of hysteresis curve gradually develops from bow-shaped to anti-S-shaped. The failure mode of specimens tend to be brittle failure. The displacement ductility factor, coefficient of average energy dissipation and other indicators decrease too.
-
Key words:
- Corrosion /
- Reinforced concrete column /
- Mechanical properties /
- Numerical simulation
-
引言
据河北省地震台网测定,2018年2月12日18时31分在河北省廊坊市永清县(39.37°N,116.67°E)发生4.3级地震,震源深度20km。此次地震未造成人员伤亡,在震中区域造成个别房屋墙皮脱落,河北省廊坊市、固安市震感较强,北京市、天津市普遍有感。地震发生后,中国地震局工程力学研究所对永清4.3级地震强震动记录进行了远程和现场记录回收。强震动记录是地震能量在地表或建筑结构上的真实反映,是工程地震、地震学研究的基础数据资料(朱升初等,2018),本文对此次地震构造背景及强震动记录进行了分析,对永清4.3级地震烈度采用建立虚拟台站方法进行了校正,为河北地区烈度调查、地震区划等提供基础资料。
1. 发震构造背景
永清4.3级地震区域构造主要以NE、NW向断裂为主,地震发生在廊固凹陷内(图 1),廊固凹陷是冀中坳陷北部的一个古近系湖盆,整个凹陷大体呈NE走向,南北长约90km,东西宽20—40km,面积约2600km2。该凹陷是受区域伸展并伴有走滑作用形成的典型箕状凹陷,整个凹陷表现为北断南超、西断东超的构造格局。由于受NE向断裂构造控制,表现出整体不均匀下降的箕状断陷特征。该凹陷内曾于1621年发生廊坊永清$ 5\frac{1}{2} $级地震(国家地震局震害防御司,1995)。
图 1 永清4.3级地震构造背景Figure 1. The seismotectonic background map of the Yongqing MS4.3 earthquake(1)施庄断裂(2)南口山前断裂(3)南口-孙河断裂(4)小汤山-东北旺断裂(5)黄庄-高丽营断裂 (6)顺义-前门-良乡断裂(7)南苑-通县断裂(8)新夏垫断裂(9)夏垫断裂(10)大兴凸起东缘断裂 (11)蓟县山前断裂(12)宝坻断裂(13)桐柏断裂(14)河西务断裂(15)天津断裂(16)大城东断裂 (17)沧东断裂(18)海河断裂(19)牛驼镇凸起东缘断裂(20)徐水断裂(21)东垒子-涞水断裂 (22)永定河断裂(23)徐水南断裂本次地震发生在廊固凹陷东侧控制断裂河西务断裂带附近,该断裂北起河西务西北,向西南经大王务,延伸至永清东码头镇,全长近50km,走向NE,倾向SE,为上陡下缓的铲形正断层,断层产状与震源机制参数基本一致(王晓山等,2018),所以永清4.3级地震发震构造可能为河西务断裂。周月玲等(2018)通过浅层地震勘探和钻孔联合剖面探测,结合年代样品测试,对河西务断裂活动性进行了综合研究,揭示断裂上断点埋深约150m或以浅,第四系底界面垂直断错为20—45m,断裂最新活动时代为晚更新世早期。
2. 观测台站及记录获取
2.1 强震动台站分布
本次地震强震动数据主要由中国地震局工程力学研究所提供,共获取3个分量自由场记录74组,常规处理后得到222条3分向加速度记录,获取记录的强震动台站均为土层台,记录台站分布如图 2所示,由图 2可知,触发台站主要分布在震中北部,南部无台站触发,这可能是因为首都圈台站较密集,而河北省境内强震动台站稀疏,或与仪器型号不同且很多台站老旧未改造有关。触发台站的震中距为29.4—185.5km,其中震中距<50km的记录有5组,震中距为50—100km的有29组,震中距为100—150km的有37组,震中距为150—200km的有3组,距震中最近的台站为长子营台(11ZZY),震中距为29.4km。虽然此次地震获取的强震动记录较全面,但由于震级较小,未造成大的灾害损失,且震中附近台站相对稀疏,能供工程使用的强震动记录较少。
2.2 强震动记录获取
在强震动记录使用分析前,首先对波形数据进行消除基线漂移,即原始强震动记录通过减去事前记录平均值的方法进行零线调整(吕俊强等,2015),通过分析计算得到震中距R<50km的5个台站强震动记录基本信息及相关参数,见表 1。其中距震中最近的长子营台记录到的最大峰值加速度为7.9cm/s2,觅子店台记录到的东西向最大峰值加速度为36.1cm/s2、南北向最大峰值加速度为40.4cm/s2,红星中学台记录到的垂直方向最大峰值加速度为10.8cm/s2,5个台站加速度时程曲线如图 3。由表 1可知,离震中较近的长子营台(震中距29.4km)较离震中较远的觅子店台(震中距38.4km)和于家务台(震中距36.8km)各分量加速度小,影响地震动强度的因素除震中距外,还包括仪器型号、场地特征条件和地质构造、场地土层结构、覆盖层厚度、地形及断裂分布等(王海云等,2010;温瑞智等,2013;郭秋娜等,2018)。本文根据温瑞智等(2015)提出的适用于我国强震动台站场地划分经验曲线,采用H/V谱比法对本次地震触发台站进行分类,分类结果见图 4。由图 4可知,记录到的本次地震台站大部分为Ⅰ类场地(51.35%),其次为Ⅱ类场地(35.14%)和Ⅲ类场地(13.51%),其中,觅子店台和红星中学台为Ⅰ类场地,于家务台和张家湾中学台为Ⅱ类场地,长子营台为Ⅲ类场地。长子营台与觅子店台和于家务台记录差异较大,这可能与场地类型有关。
表 1 强震动记录(震中距<50km)及相关参数Table 1. Strong motion records(epicenter distance < 50km) and related parameters台站名称 台站代码 仪器型号 震中距/km PGA/cm·s-2 EW SN UD 觅子店 11MZD REFTEK/SLJ-100 38.4 36.1 40.4 8.5 红星中学 11YHZ REFTEK/SLJ-100 43.4 -5.4 -4.6 10.8 于家务 11YJW ETNA/ES-T 36.8 -21.7 -18.5 -8.1 张家湾中学 11ZJW REFTEK/SLJ-100 47.9 -5.8 4.2 -4.4 长子营 11ZZY REFTEK/SLJ-100 29.4 5.0 -7.9 -3.5 3. 谱分析
加速度反应谱是目前国内外工程抗震设计的重要依据,反应谱形状、幅值及卓越周期的差异体现了场地条件、震级对地面运动反应谱的影响(温瑞智等,2013;徐钦等,2017)。为避免地震时建筑结构发生共振,确保结构安全,有必要对加速度反应谱成分进行分析。本文选取震中距<50km的5个典型强震动台站进行频谱分析,阻尼比取0.05,得到本次地震EW、SN、UD向加速度反应谱(图 5)。
由图 5可知,觅子店台记录到了本次地震EW向和SN向加速度反应谱最大值,红星中学台记录到了本次地震SN向加速度反应谱最大值。觅子店台与长子营台在EW和SN向差别最大,觅子店台与张家湾中学台在SN向差别最大。5个典型台站记录到的加速度反应谱高频成分丰富,谱卓越周期主要为0.1—0.3s,0.3s后迅速下降。根据现场调查,震区主要房屋类型为土木结构、木房顶砌体结构、预制板砌体结构、钢结构简易房屋4种,本次地震除少数土木结构房屋原有裂缝扩大、墙皮开裂外,其他类型房屋未发现新开裂、屋顶掉瓦、墙皮掉灰现象,中小城市主要建筑物结构自振周期为0.3—1s(王文才等,2018),因此本次地震对该震区自振周期范围内的建筑物破坏影响较小。通过对比不同震中距台站各方向反应谱特征周期发现,震中距对特征周期的影响不明显。
傅氏谱能反映地震波在频域上的幅值特性,地震动记录中包含多种振动频率,不同地区的主要振动频率不同,结构自身具有不同频率、幅频响应,因而在不同频率振动作用下,结构振动不同。为避免地震作用下结构发生共振现象,确保结构安全,有必要对地震频率成分进行分析(冉志杰等,2012)。
本文计算了永清4.3级地震震中距50km以内的5个典型台站强震动记录傅氏谱(图 6),由图 6可知典型台站傅氏谱谱型主要为单峰值型,强震动记录EW分量傅氏谱最大峰值为3.5—8.5Hz,SN分量傅氏谱最大峰值为2.9—6.2Hz,UD分量傅氏谱最大峰值为2.1—9.8Hz。总体来看,UD分量频率范围分布较宽,EW分量傅氏谱最大峰值大于SN分量和UD分量傅氏谱峰值。一般而言,随着震中距的增大,地震动记录长周期(低频)越显著,但本次地震5个典型台站特征不明显,可能与仪器型号及场地特征条件有关。另外,EW分量傅氏谱最大峰值随着震中距的增大呈先减小后增大的趋势,而SN分量和UD分量无该特征。
本文以东西向为正方向,统计峰值加速度反应谱与永清4.3级地震触发台站同震中所处方位角和震中距的关系(图 7)。由图 7可知,峰值加速度反应谱在同一震中距下整体随着方位角的增大而增大,达70°左右后下降,这可能与附近河西务断裂构造走向有关;峰值加速度反应谱值随着震中距的增加而减小。
4. 地震烈度分析
地震烈度是衡量地震破坏和地震动强弱程度的重要指标。Wald等(1999)认为当烈度较小时,用峰值加速度计算烈度较合适,因为人对加速度较敏感。随着强震动观测技术的发展,由于以各观测点峰值加速度作为依据可真实反映地震发生时各点地面运动强度(田秀丰等,2013),基于地震动记录快速确定地震仪器烈度,越来越受到研究者的重视(崔建文等,2016)。目前国内外基于强震动记录资料的地震烈度计算方法较多,不同方法利用不同的地震动参数,包括时间、频度和强度等,建立地震烈度与地震动参数之间的关系。本文采用美国ShakeMap烈度计算方法(李俊等,2010),计算永清4.3级地震事件记录台站烈度,并采用克里格插值方法进行网格化插值,绘制等值线图(图 8)。Wald等(1999)的研究认为,当烈度低于Ⅴ度时峰值加速度PGA与烈度相关性明显高于峰值速度PGV,而当烈度为Ⅴ—Ⅶ度时,PGA和PGV与烈度相关程度相当,考虑本次地震烈度较低,故采用峰值加速度计算烈度。由图 9可知,本次地震仪器烈度等震线长轴沿NW走向,且极震区位于震中东北部,而不是震中附近,但由现场宏观烈度调查(图 8)与震源机制解(王晓山等,2018)可知该地震等震线长轴沿NE向展布,这种差异现象的出现可能与台站分布不均匀有关,震中南部无强震动记录,所以图 8中震中北部附近出现数据空白现象,不利于获取实际烈度长轴方向真实性。
可靠快速生成地震烈度可为人员伤亡估计、应急救援决策和工程抢险修复决策提供重要依据,为解决由于缺少地震动观测值出现空白区域现象及烈度长轴方向和实际方向不一致的问题,本文首先利用永清4.3级地震中获取的所有台站加速度峰值,按地震动峰值衰减模型(王国新,2001)确定该地区地震动衰减关系,然后在整个地震区内建立空间随机假设台站进行补点插值,得到新的地震烈度分布图。地震动峰值衰减模型(王国新,2001)如下:
$$ \lg Y = {{\rm{C}}_1} + {{\rm{C}}_2}M + {{\rm{C}}_3}M\lg \left({R + {R_0}} \right) + {{\rm{C}}_4}\lg \left({R + {R_0}} \right) \pm \varepsilon $$ (1) 式中,Y为PGA(峰值加速度),M为震级,R为震中距,C1、C2、C3、C4为待求常数,R0为近场饱和因子,ε为统计分析误差项。采用通用全局优化法拟合得到地震动峰值加速度随震级和震中距的变化规律,考虑本次地震震级较小,故选取震中距<100km的Ⅰ类场地强震动数据进行拟合,其中水平向按EW、NS向加速度峰值矢量拟合回归分析,拟合曲线如图 9,计算得到的衰减关系式为:
$$ \lg Y = 1.6882 + 0.3926M - 0.1546M\lg \left({R + 5} \right) + 0.7000\lg \left({R + 5} \right) \pm 0.2287 $$ (2) 峰值加速度随震中距的变化如图 10所示,由图 10可知峰值加速度随着震中距的增加而减小。同一震中距不同台站测得的峰值加速度离散性较大。通过拟合各台站实际PGA值,与俞言祥等(2013)地震动衰减关系进行对比,可知20km范围内2条曲线吻合较好,而20km外拟合曲线呈衰减相对缓慢状态,这可能与本次地震震级较小、记录台站均为土层台、远场有一定放大作用有关(王海云,2011)。由观测值得到的预测方程虽存在一定离散性,但整体趋势与实际结果保持一致。
在震区空白区域随机产生N个空间坐标的假设台站,按式(2)计算水平峰值加速度,按美国ShakeMap烈度计算方法转换成仪器烈度,绘制增加假设台站的烈度分布图(图 11),由图 11可知,增加假设台站的仪器烈度分布图与现场烈度调查分布图极震区方向基本一致,进一步验证了河西务断裂为本次地震发震构造的可能,解决了台站分布不均匀导致地震烈度影响场计算出现缺值现象的问题,为以后缺少台站记录震区提供准确快速制作烈度分布图的思路,为震害调查和地震应急救援提供重要依据。
5. 结论
本文介绍了河北永清4.3级地震发震构造背景,并对国家强震动台网中心在本次地震中获取的加速度记录特征进行了分析,为河北省地震动特征和震害调查提供了一定参考,结论如下:
(1) 河北永清4.3级地震在廊固凹陷东侧控制断裂河西务断裂带附近,地震烈度极震区走向与该断裂展布方向一致,结合震源机制解初步判定本次地震发震构造为河西务断裂。
(2) 74个强震动台站获取到强震动加速度记录,震中距位于29.4—185.5km,其中觅子店台记录到的本次地震东西向最大峰值加速度为36.1cm/s2,垂直方向最大加速度为40.4cm/s2,于家务台记录到的南北向最大峰值加速度为18.5cm/s2。采用H/V谱比法对触发台站进行分类,可知Ⅰ类场地居多。
(3) 通过对5个典型台站反应谱分析发现永清4.3级地震加速度反应谱高频成分丰富,谱卓越周期主要为0.1—0.3s。傅氏谱谱型主要为单峰值型,UD分量频率范围分布较宽,EW分量傅氏谱最大峰值大于SN分量和UD分量。
(4) 通过分析峰值加速度反应谱与震中距、方位角的关系,发现峰值加速度反应谱整体随着方位角的增大而增大,达70°左右后下降,这可能与附近河西务断裂构造走向有关,有待进一步研究。其次,峰值加速度反应谱值整体随着震中距的增加而减小。
(5) 由本次地震仪器烈度分布图可知,台站分布不均匀性直接影响烈度分布的真实性,为避免此类问题,拟合出该地区地震动衰减关系,随机产生空间假设台站进行弥补,为以后缺少台站记录震区提供准确快速制作烈度分布图的思路,为震害调查和地震应急救援提供重要依据。由于收集到的强震动数据有限,文中仅采用本次地震获取的峰值加速度拟合当地地震动衰减关系,拟合结果可能存在误差,今后将收集更多的强震动数据进行拟合,以校正该地区地震动衰减关系。
-
表 1 CDP模型塑性参数取值
Table 1. The plastic parameters of the CDP model
参数 Ψ/° φ fb0/fc0 Kc μ ωt ωc 数值 33 0.1 1.16 0.667 0.001 0 0.3 表 2 承载力、变形参数试验值与模拟值对比
Table 2. Comparison of bearing capacity and deformation parameters between numerical simulation and physical test
开裂荷载/kN 屈服荷载/kN 极限荷载/kN 延性比 侧移角/rad 试验实测 20.00 41.78 50.15 4.46 0.038 有限元计算 20.47 43.41 53.85 4.75 0.041 相对误差/% 2.35 3.90 7.38 6.56 7.89 表 3 滞回环能量试验值与模拟值对比
Table 3. Comparison of hysteresis loop energy between numerical simulation and physical test
循环次数 试验值/kN·mm 模拟值/kN·mm 相对误差/% 1st 201.45 178.19 -11.54 2nd 590.07 463.68 -21.42 3rd 730.58 710.88 -2.69 4th 828.98 951.38 14.76 5th 1187.56 1270.92 7.02 6th 1929.19 1936.33 0.37 7th 2258.10 2332.90 3.31 8th 2825.36 3000.31 6.19 表 4 锈蚀钢筋计算参数
Table 4. Calculation parameters of corroded steel bars
钢筋截面锈蚀率/% 钢筋直径/mm 屈服强度/Mpa 极限塑性变形率/% 粘结强度降低系数 0 14.00 415.60 10.44 1 5 13.64 410.55 9.80 0.664 10 13.28 404.93 9.16 0.369 15 12.91 398.66 8.52 0.206 20 12.52 391.60 7.87 0.115 -
程玲, 贡金鑫, 李颖, 2012.基于Pushover方法分析的受腐蚀钢筋混凝土柱抗震性能评价.振动与冲击, 31(10):19-23. http://d.old.wanfangdata.com.cn/Periodical/zdycj201210005 丁发兴, 余志武, 欧进萍, 2008.混凝土单轴受力损伤本构模型.长安大学学报(自然科学版), 28(4):70-73. doi: 10.3321/j.issn:1671-8879.2008.04.016 蒋连接, 袁迎曙, 2011a.反复荷载下锈蚀钢筋混凝土压弯构件刚度退化的试验研究.混凝土, (5):3-6. http://d.old.wanfangdata.com.cn/Periodical/hnt201105002 蒋连接, 袁迎曙, 2011b.锈蚀钢筋混凝土压弯构件的恢复力模型.混凝土, (6):29-32, 40. http://d.old.wanfangdata.com.cn/Periodical/xajzkjdx200502002 李小军, 于爱勤, 甘朋霞等, 2008.汶川8.0级地震北川县城区灾害调查与分析.震灾防御技术, 3(4):352-362. doi: 10.3969/j.issn.1673-5722.2008.04.004 李悦, 颜超, 辜中伟, 2013.钢筋混凝土锈蚀膨胀力的数值模拟研究.建筑技术, 44(3):198-201. http://d.old.wanfangdata.com.cn/Periodical/jzjs201303001 刘婕, 2016.考虑钢筋局部锈蚀影响的混凝土桥梁抗震性能研究.杭州:浙江大学. 龙渝川, 李正良, 2011.模拟混凝土滞回行为的各向异性损伤模型.工程力学, 28(8):62-69. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201101701170 马颖, 贡金鑫, 赵顺波, 2011.锈蚀钢筋混凝土圆柱抗震性能的试验研究.防灾减灾工程学报, 31(2):152-159. http://d.old.wanfangdata.com.cn/Periodical/dzxk201102005 曲哲, 叶列平, 2009."破坏-安全"结构抗震理念及其应用.震灾防御技术, 4(3):241-255. doi: 10.3969/j.issn.1673-5722.2009.03.001 曲哲, 2010.摇摆墙-框架结构抗震损伤机制控制及设计方法研究.北京:清华大学. 史庆轩, 牛荻涛, 颜桂云, 2000.反复荷载作用下锈蚀钢筋混凝土压弯构件恢复力性能的试验研究.地震工程与工程振动, 20(4):44-50. doi: 10.3969/j.issn.1000-1301.2000.04.007 邢国华, 牛狄涛, 2014.锈蚀钢筋混凝土梁的受弯分析模型.中南大学学报(自然科学版), 45(1):193-201. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=SciencePaper201404040000013073 闫培雷, 孙柏涛, 陈洪富, 2010.汶川地震震中某钢筋混凝土框架结构的非线性地震反应分析.震灾防御技术, 5(3):364-369. doi: 10.3969/j.issn.1673-5722.2010.03.013 阎红霞, 杨庆山, 张丽英, 2010.ABAQUS在超高层结构动力弹塑性分析中的应用.震灾防御技术, 5(1):108-115. doi: 10.3969/j.issn.1673-5722.2010.01.013 张誉, 蒋利学, 张伟平等, 2003.混凝土结构耐久性概论.上海:上海科学技术出版社. Berto L., Vitaliani R., Saetta A., et al., 2009. Seismic assessment of existing RC structures affected by degradation phenomena. Structural Safety, 31(4):284-297. doi: 10.1016/j.strusafe.2008.09.006 Bhargava K., Ghosh A. K., Mori Y., et al., 2008. Suggested empirical models for corrosion-induced bond degradation in reinforced concrete. Journal of Structural Engineering, 134(2):221-230. doi: 10.1061/(ASCE)0733-9445(2008)134:2(221) Kivell A., Palermo A., Scott A., 2015. Complete model of corrosion-degraded cyclic bond performance in reinforced concrete. Journal of Structural Engineering, 141(9):04014222. doi: 10.1061/(ASCE)ST.1943-541X.0001195 Mehta P K., 1991. Durability of concrete-fifty years of progress. In: Malhotra V. M., Proceedings of the Second International Conference on Durabilitay of Concrete. Detroit: American Concrete Institute, 795-802. Yalciner H., Sensoy S., Eren O., 2012. Time-dependent seismic performance assessment of a single-degree-of-freedom frame subject to corrosion. Engineering Failure Analysis, 19:109-122. doi: 10.1016/j.engfailanal.2011.09.010 -