浙江大学学报(工学版), 2019, 53(4): 692-701 doi: 10.3785/j.issn.1008-973X.2019.04.010

土木工程、海洋工程

考虑河谷地形影响的多层框架结构地震易损性

黄博,, 廖凯龙, 赵宇,, 蒋建群

Seismic fragility analysis of multi-storey frame structure considering effect of valley terrain

HUANG Bo,, LIAO Kai-long, ZHAO Yu,, JIANG Jian-qun

通讯作者: 赵宇,男,副教授. orcid.org/0000-0003-0453-1960. E-mail: zhao_yu@zju.edu.cn

收稿日期: 2018-06-30  

Received: 2018-06-30  

作者简介 About authors

黄博(1973—),女,副教授,从事土动力学和实验土力学研究.orcid.org/0000-0002-7293-8618.E-mail:cehuangbo@zju.edu.cn , E-mail:cehuangbo@zju.edu.cn

摘要

针对杭嘉湖地区冲淤河谷地形分布广泛的特殊情况,结合当地土体工程特性和抗震设防区划设置现状,基于动力时程分析以及概率统计方法,给出适用于杭嘉湖地区多层框架结构的地震易损性曲线,对比不考虑地基、均质地基和河谷地基条件下多层框架结构的动力分析结果. 研究土-结动力相互作用以及局部河谷地形对结构动力响应和易损性的影响,考察房屋结构位于河谷不同位置处的差异. 结果表明,土-结动力相互作用对结构的动力分析结果有显著影响,不考虑土-结动力相互作用将低估房屋结构变形,易损性性能水平超越概率偏低,偏差可达7.9%,结果偏于危险;局部河谷地形对地震波独特的聚焦效应会增大结构变形,地震易损性超越概率与均质地基相比增幅可达12.1%,且河谷不同位置处的聚焦效果不同,在选取的地震波输入下表现为房屋结构距河谷中心越近,结构变形越大、地震易损性超越概率越大.

关键词: 地震易损性 ; 多层框架结构 ; 土-结动力相互作用 ; 局部河谷地形 ; 聚焦效应

Abstract

The seismic fragility curves of multi-storey frame structures suitable for Hangjiahu area were given based on dynamic time history analysis and probability statistics method combining with the local soil engineering characteristics and the fact of seismic fortification zones in view of the special situation that the scouring and silting valley topography is widely distributed in Hangjiahu area. The dynamic analysis results of multi-storey frame structures without considering foundation, homogeneous foundation and River Valley Foundation were compared. Results show that soil-structures dynamic interaction has a significant influence on the dynamic analysis results of structures. The deformation and the exceeding probability of fragility performance level of buildings were underestimated without considering soil-structures dynamic interaction. The deviation can reach 7.9%, which will lead to a dangerous result. The unique focusing effect of local valley topography on seismic waves will increase the structural deformation. The exceeding probability of seismic vulnerability can increase by 12.1% compared with homogeneous foundation, and the focusing effect is different at different locations of the valley. Under the selected seismic wave input, the closer the building structure is to the center of the valley, the greater the structural deformation and the greater the probability of exceeding the seismic fragility are.

Keywords: seismic fragility ; multi-storey frame structure ; soil-structure dynamic interaction ; local valley terrain ; focusing effect

PDF (1422KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

黄博, 廖凯龙, 赵宇, 蒋建群. 考虑河谷地形影响的多层框架结构地震易损性. 浙江大学学报(工学版)[J], 2019, 53(4): 692-701 doi:10.3785/j.issn.1008-973X.2019.04.010

HUANG Bo, LIAO Kai-long, ZHAO Yu, JIANG Jian-qun. Seismic fragility analysis of multi-storey frame structure considering effect of valley terrain. Journal of Zhejiang University(Engineering Science)[J], 2019, 53(4): 692-701 doi:10.3785/j.issn.1008-973X.2019.04.010

房屋结构的地震易损性分析是采用概率统计方法预测结构在不同强度的地震作用下,达到或超过各级性能损伤水平的概率. 分析时,首先确定系统的随机变量并建立易损性分析样本,再输入不同的地震波进行动力响应分析,计算不同地震条件下结构破损状态的概率分布[1]. 目前,采用易损性分析方法评估地震引起的损伤,已被广泛地应用于各种房屋结构类型的抗震能力评估中. 地震易损性曲线是易损性的图形化表示.

结构工程地震灾变与岩土体特性密不可分,根据相关调查资料[2]可知,地震波通过岩土体传播和频谱放大,对上部结构的破坏更严重. 地基土体动力特性直接决定了结构基础的破坏形式,特别是当上部结构与地基刚度相差较大时,地基对结构动力响应的影响更突出[3]. 人们对地震波传播过程中地基的放大效应已有充分的认识,并进一步认识到局部地形作为不均匀散射体,会使地震波发生复杂的散射、波形转换及相干作用,使地震波产生局部放大或缩幅效应,从而对地震波特性产生关键性影响[4]. Tafazzoli等[5]对含矩形软弱土区域的地基进行动力分析发现,当矩形区域的宽度与深度之比小于4时,局部地形对地震波的放大作用十分显著. 金峰等[6]对含半椭圆形软弱土区域的地基动力分析发现,沉积层表面各点的响应显著放大. 尤红兵等[7]利用间接边界元方法,在频域内求解层状场地中任意形状沉积河谷对平面SH波的散射效应时发现,沉积河谷的散射效应使得地表位移幅值显著增大.

目前,对结构地震易损性的研究多侧重于结构类型的多样化、自身模型的精细化等方面[8-11],将地震波直接作用在结构底部进行分析,综合考虑地震波-土-结构三者耦合的结构地震易损性分析仍不普遍,该方法忽视了场地放大效应和局部特殊地形对结构易损性的影响.

在新一代《中国地震波参数区划图(GB 18306—2015)》中,浙江省全面消除了非设防区. 杭嘉湖地区经济发达、人口密集、灾害波及效应大,该区域大部分地区从6度设防调整至7度设防[12],这意味着对区域内结构抗震性能要求的提高[13]. 杭嘉湖平原为典型的冲积平原,河流、河谷众多,所形成的河流淤积土具有强度低、土颗粒均匀、孔隙率高和地下水位高等特点;再加上杭嘉湖地区的城镇化程度高,回填池塘、河道并在回填后的地基上构筑房屋的情况非常普遍,导致典型的局部河谷、盆地特征. 基于杭嘉湖地区的典型地层、地质特点,研究地基软、硬土分布不均匀造成的局部河谷地形对结构动力响应和易损性的影响,对评估杭嘉湖地区房屋结构的抗震性能具有重要意义.

本文结合杭嘉湖地区的地质条件、结构特点以及抗震设防现状,选取18条地震波输入,对比不考虑地基、均质地基和河谷地基条件下典型民用住宅结构的动力响应. 考察房屋结构位于河谷不同位置处的响应差异,基于动力时程分析结果和概率统计方法,提出适用于杭嘉湖地区多层框架结构的地震易损性曲线. 以上工作为杭嘉湖地区典型的多层框架结构抗震性能的评估和抗震加固、维修决策提供了参考依据.

1. 分析模型

1.1. 结构模型

杭嘉湖地区历来经济发达,一方面近年来兴建了大量的高层、超高层公建、民建结构,另一方面存在大量老旧小区. 这些民用住宅以砖混结构与框架结构为主,楼层一般为5~7层,地基处理简单,抗震性能较差.

将这类典型民用住宅简化为5层3跨框架结构(下文中简称房屋结构),建立对应的二维模型,层高3.6 m,总高18 m. 梁的截面尺寸为200 mm×400 mm,C30混凝土、HRB335钢筋. 柱的截面尺寸为600 mm×600 mm,所用材料与梁一致. 基础形式为条形基础,埋深为1.5 m,C30混凝土,具体尺寸如图1所示. 楼面和屋面荷载取7 kN/m2.

图 1

图 1   房屋结构尺寸与条基尺寸

Fig.1   Size of house structure and strip footing


结构的梁、柱构件均采用Timoshenko梁单元,混凝土和钢筋材料分别采用编辑材料属性和修改关键字的方式进行设置. 其中混凝土采用ABAQUS自带的混凝土损伤塑性(concrete damage plasticity)模型,该模型考虑损伤效应,能够较好地模拟地震作用下的混凝土力学行为[14]. 基本参数的取值[15]如下:混凝土弹性模量Ee取为2.2×104 MPa,泊松比μ取为0.2,混凝土单轴受压峰值应力fc,r和应变εc,r分别取为3.0 MPa和1.36×10−4,混凝土单轴受拉峰值应力ft,r和应变εt,r分别取为39.6 MPa和3.15×10−3. 为了考虑钢筋混凝土受拉硬化效应和箍筋的约束作用,根据高向玲等[15]提出的方法对混凝土的应力-应变关系进行修正. 钢筋通过编辑*rebar关键字在梁单元中实现插入,用来模拟钢筋混凝土构件中的纵筋,截面刚度直接根据钢筋纤维和混凝土纤维的单轴应力-应变关系,依据平截面假定积分得到[16]. 条基采用四节点缩减积分平面应变单元,混凝土材料的设置方式和材料参数与房屋结构构件一致. 考虑到实际施工中,柱-条基之间的连接方式多与柱-承台之间相同[17],因此参考庄海洋等[18]《土-桩-隔震结构动力相互作用》第8章中结构与承台的绑定连接方法,定义房屋结构与条基之间的连接.

为了使二维模型真实地反映三维结构的动力响应,对房屋结构的材料属性沿进深方向进行等效处理,即把相关参数除以布置间距换算成每延米的对应量,需要等效的参数包括弹性模量、单轴应力-应变曲线中的应力以及材料密度.

1.2. 地基模型

为了研究地基特性对房屋结构动力响应及易损性的影响,分析2种不同的地基:均质地基和河谷地基,同时与不考虑土-结动力相互作用的房屋结构动力响应进行比较. 不考虑地基的模型只建立房屋结构,约束柱底非激振方向上的自由度.

在均质地基中,将地基简化为单一均质粉质黏土,地基厚100 m,宽500 m,如图2(a)所示;考虑沉积河谷沿河流方向延伸,顺流向接近均质地基,因此取垂直流向的横剖面研究局部河谷地形的影响. 将河谷地基中的软土淤积区域简化为半椭圆型,顶面最宽处为150 m,深度为50 m,如图2(b)所示. 淤积软土周边为粉质黏土. 地基土体选用Mohr-Coulomb弹塑性模型,采用杭嘉湖地区常见的土层参数,参数取值来自文献[19]及其他工程资料,如表1所示. 地基土体采用四节点缩减积分平面应变单元,为了避免缩减积分导致的“沙漏效应”,在网格模块中勾选Enhance命令来解决沙漏问题. 为了研究房屋结构处于河谷不同位置处的动力响应特征,在河谷地基模型中分别考虑结构在河谷中心、河谷边沿以及处于两者之间的3种情况. 为了表述方便,将结构轴心离河谷中心的距离定义为x,用x=0、60和30 m分别表示以上3种情况.

表 1   地基土体参数

Tab.1  Foundation soil parameters

土层 ρ/(kg·m−3) E/MPa μ c/kPa φ/(°)
淤泥质土 1 820 149 0.49 5 20
粉质黏土 2 000 329 0.49 8 25

新窗口打开| 下载CSV


图 2

图 2   地基模型概化图

Fig.2   Generalized figure of foundation model


条基与地基土体的接触采用《ABAQUS在岩土工程中的应用》[20]第5章介绍的接触面对法. 选择刚度较大的条基表面作为接触面对的主控面,接触面之间的法向行为定义为“硬接触”,限制基础与地基土体在接触面上的法向行为,避免相互穿透,两者接触时才能产生法向接触应力,分离时应力为0;切向行为采用“罚函数”中的库仑摩擦模型,极限切应力可以表示为

${\tau _{{\rm{crit}}}} = \mu {{p}},$

$\mu = \tan\; (0.75\varphi ).$

式中:μ为条基和地基土体间的摩擦系数,p为接触面的法向应力,φ为对应土层的内摩擦角.

由于瑞利阻尼为正交阻尼,可以解耦复杂的动力方程,且具有形式简单、数学处理方便等优点,因此采用瑞利阻尼模拟地基和结构的阻尼. 阻尼矩阵 ${{C}}$可以表示为

${{C}} = \alpha {{M}} + \beta {{K}}.$

式中:α为质量比例阻尼,β为刚度比例阻尼.

αβ的计算方法较多,选用Hudson等[21]提出的QUAD4法进行计算:

$\alpha = 2\zeta \frac{{{\omega _{\rm{1}}}{\omega _{\rm{2}}}}}{{{\omega _1}{\rm{ + }}{\omega _{\rm{2}}}}},$

$\beta = 2\zeta \frac{1}{{{\omega _{\rm{1}}}{\rm{ + }}{\omega _{\rm{2}}}}}.$

式中:ω1为模型的基频,ω2=1,其中n为大于ωe/ω1的奇数,ωe为地震波的主频; $\zeta $为阻尼比,Nejati等[22]建议土体阻尼比取2.6%;黄宗明等[23]建议结构阻尼比取5%.

该方法综合考虑了模型的频率特性和地震波的频谱特性,被证实在计算地基土体和结构的动力响应时可以获得较精确的结果[24-25].

1.3. 地震波输入和边界条件

分析时,仅考虑垂直入射SV波的作用. 不考虑土-结动力相互作用时,将地表的地震波记录以加速度的形式直接施加在房屋结构柱底端,来实现外部激励. 考虑土-结动力相互作用时,多通过等效线性化法[26]或时域非线性逐步积分法[27]等,将地表地震波记录反演至基岩(模型底部)处作为激励输入,保证地震波传播到地表时与实测记录相当. 本文利用等效线性地震场地响应分析软件EERA的反演功能,将均质地基地表地震波反演至地基底部,再将得到的加速度时程积分为速度时程,根据式(6)将速度转化为结点力施加在地基的底部边界上[28],作为激励输入. 在分析局部河谷地形的影响时,地震波输入与均质地基相同,考察相同输入下局部地形对地震波的放大作用. 利用等效线性化方法实现地震波反演的优势是计算原理和迭代思路简单,土体动力参数容易获取,地震分析计算量少,已被证明在数值模拟中结果合理[27].

${\sigma _{{xz}}} = \rho {c_{\rm{s}}}\frac{\partial {{ u}_{{z}}}}{{\partial {{t}}}}.$

式中:uz为地基底部边界上节点的切向运动位移,ρ为质量密度,cs为S波的波速.

为了消除边界反射,模拟无限地基辐射阻尼效应,在地基有限域的左、右以及底部边界上,设置无限单元(简称无限元)吸收反射波. 无限元的概念最早由Ungless[29]提出,该方法不仅能够模拟远场地基对地震波能量的吸收作用,而且能够准确模拟无穷远处位移为零的边界条件,采用无限元能够减少单元数量,节约计算成本[30].

为了考察无限元边界能量吸收的效果,在图2(a)所示的均质地基中,将地基视为弹性材料,密度ρ为2 000 kg/m3,弹性模量E为170 MPa,泊松比为0.3,剪切波速为188 m/s2. 根据上述地震波输入方法,在模型底部输入峰值为1 m/s2、周期为0.25 s的单个半正弦波,对比在固定边界和无限元边界条件下地表中心处的加速度响应,如图3所示.

图 3

图 3   地表中心点加速度时程曲线

Fig.3   Acceleration time history curve of surface center point


当地震波首次传播至地表时,两条曲线几乎重叠,加速度信号经自由表面反射形成峰值为2 m/s2的半正弦波. 之后,2条曲线的差异较大,在固定边界下,由于能量无法透射,地震波在计算区域内来回振荡,地表多次接收到加速度信号,产生多个峰值,从而引起较大误差;在无限元边界下,地震波反射至底部边界后能够很好地被吸收,地表不再接收到地震波信号,加速度趋近于0. 可见,无限元边界吸收能量的效果较好,应用无限元边界分析土-结动力相互作用问题是可靠的.

1.4. 地震波选取

目前,在进行地震易损性分析时,Zhao等[31-32]多从地震信息的角度出发,依据地震波的震级Ms、震中距R等参数选取地震波. 本文选取的18条地震波的MsR变化范围主要集中在2个区域,如表2所示. 一是R约为10 km、Ms为4.9~7.6的近场地震波,这主要是考虑到杭嘉湖地区内布有肖山-球川和嘉兴-松江等断裂带[33];另一个是R为28.6~70.4 km、Ms为6.9~7.4的地震波,由于在输入时进行调幅处理,主要是考虑到距离较远的中远场地震波对杭嘉湖地区房屋结构的影响.

表 2   实际地震记录

Tab.2  List of earthquake records

地震名称 年份 Ms R/km
EL Centro 1979 6.6 2.5
Chichi 1999 7.6 7.3
Taft 1952 7.4 41.0
Kobe 1995 6.9 15.5
Loma Prieta1 1989 6.9 70.4
Loma Prieta2 1989 6.9 28.6
Loma Prieta3 1989 6.9 35.5
Northridge1 1994 6.7 7.1
Northridge2 1994 6.7 6.2
Whittier Narrows 1987 5.3 4.4
Coalinga 1983 6.0 10.0
Kocaeli 1999 7.4 67.5
Mammoth Lake 1980 4.9 7.5
San Fernando 1971 6.6 12.0
Coyote 1979 5.7 12.6
Nahanni 1985 6.8 6.0
Duzce 1999 7.2 8.2
Cape Mendocino 1992 7.1 9.5

新窗口打开| 下载CSV


图4给出所选18条地震波在阻尼比为5%时的加速度反应谱. 为了便于观察、比较,将各地震波的谱加速度在均质地基模型基本周期(T1=1.78 s)处按比例调整为10 m/s2. 由图4可见,调整后的加速度信号在0.01 s到T1的区间内表现出较大的离散性,这体现了输入地震波信号的随机性.

图 4

图 4   所选18条地震波反应谱

Fig.4   Response spectrum of 18 seismic waves


2. 动力响应

2.1. 结构自振周期

不考虑地基和均质地基的一阶自振周期分别为0.88和1.78 s,河谷地形3处的自振周期均为1.79 s. 可以看出,不考虑地基模型的房屋结构一阶自振周期约为其他模型的一半,均质地基模型与河谷地基模型之间的差异很小. 这是由于考虑地基后,整个模型质量大大增加,导致自振周期大幅上升,该现象在考虑土-结动力相互作用的动力分析中普遍存在[34]. 此外,局部河谷地形对模型自振周期的影响不明显;房屋结构在河谷地基中的所处位置对模型的自振周期没有影响.

2.2. 局部河谷地形下地基的放大效应

地震波的反演分析表明,均质地基对地震波具有滤波与放大作用. 以Northridge1波为例,图5给出均质地基模型地面处和反演分析得到的基底处的加速度曲线. 当地表加速度峰值为1.0 m/s2时,反演得到基底加速度峰值为0.7 m/s2,放大系数为1.4. 在不同的地震波输入下,均质地基模型的加速度放大系数略有不同,变化范围为1.07~1.51. 如图6所示为图5中加速度时程对应的傅里叶谱. 地基对加速度的不同频率成分作用不同,如图6所示,Northridge1波频谱变化特征为低频加震和高频减震,其中频率为0.3~5.0 Hz的低频分量加震现象比较明显,频率>7.5 Hz的高频分量减震现象比较明显. 此外,不同地震波的频谱在均质地基的作用下表现出类似的变化规律,但加震和减震的频率范围不同.

图 5

图 5   均质地基模型地面、基底处的加速度时程曲线

Fig.5   Acceleration time history curve of homogeneous foundation model at ground and base


图 6

图 6   均质地基模型地面、基底的加速度傅里叶谱

Fig.6   Acceleration fourier spectrum of homogeneous foundation model at ground and base


由于局部河谷地形的特殊构造,河谷地基对加速度的放大作用比均质地基更显著,且河谷不同位置处的地震波放大效果明显不同. 该现象被称为聚焦效应[35],王力等[36]通过大型振动台试验对比河谷地基与均质地基的加速度响应,证实了该效应的存在.

在Northridge1波的输入下,均质地基模型柱底加速度峰值为1.0 m/s2,河谷地基模型x=0、30、60 m处加速度峰值分别为12.9、12.5、10.3 m/s2,可见河谷地基结构柱底处的加速度峰值整体上高于均质地基,且不同位置的地震放大效果不同.

图7给出均质地基、河谷地基模型柱底相对于底部输入边界的加速度反应谱比. 图中,ξ为放大系数. 可以看出,河谷地基与均质地基对底部输入的地震波均都有放大效应,且河谷地基放大系数均大于均质地基,这一现象与袁晓铭等[37]的研究成果一致. 在周期小于0.1 s的短周期范围内,周期的改变对各模型的放大系数无显著影响;当周期为0.62和1.5 s时反应谱比曲线出现2个峰值,在曲线最高点处,河谷地基x=0、30、60 m处的放大系数分别为4.12、3.36、2.84,体现出河谷地基以河谷中心为聚焦点的聚焦效应.

图 7

图 7   均质地基、河谷地基模型柱底加速度反应谱比

Fig.7   Acceleration response spectral ratio of homogeneous foundation and valley foundation model at bottom of column


在不同的地震波输入下,均呈现出地震波附加放大的聚焦效应,但聚焦点的位置因地震波的主频不同而变化,这一现象与Semblat等[38]的研究成果一致. 针对本文研究的工况,各地震波的聚焦点多出现在河谷中心处,少数出现在河谷边沿处.

2.3. 结构变形

以Northridge1波输入为例,说明局部河谷地形对房屋结构变形的影响. 计算表明,房屋结构各层的最大侧向位移发生在结构端点处,因此将结构各层左侧端点作为参考点进行分析. 图8(a)给出不考虑地基、均质地基和河谷地基模型的结构侧向最大位移包络线. 图中,L为层数,0层指柱底位置. 从图8可以看出,侧向位移有一个共同特征,即1层处楼房侧移最大、2~5层依次减小. 由于没有约束各模型柱底水平方向的自由度,允许柱底发生侧向位移,导致0层的侧向最大位移受到较大影响. 为了消除这一影响,给出如图8(b)所示的层间最大位移包络线. 层间最大位移是地震作用过程中相邻两参考点发生的最大相对位移. 与不考虑地基模型相比,考虑土-结动力相互作用后,层间最大位移发生位置从4层变为1层,这一结构变形规律与汶川地震灾害调查中的框架结构破坏现象一致[39];层间最大位移从2.53 mm变为3.33 mm(均质地基),放大了1.3倍,河谷地基模型中层间最大位移增大更明显,x=0、30、60 m处的层间最大位移分别为不考虑地基模型的1.8、1.5、1.4倍,可见局部河谷地形对结构变形具有明显的放大作用,且房屋结构离河谷中心越近,结构变形越大. 在其他地震波的作用下,呈现出类似的特征,表现为:在河谷地基中,房屋结构距聚焦点越近,结构变形越大.

图 8

图 8   结构变形曲线

Fig.8   Structural deformation curve


3. 房屋结构的地震易损性

3.1. 地震易损性分析方法

杭嘉湖地区的设防烈度为6、7度,根据建筑抗震设计规范[40]中所规定的抗震设防烈度和设计基本地震加速度峰值(PGA)取值的对应关系,同时考虑到易损性曲线的完整性,将所选18条地震波的PGA分别调整为0.5、1、2、3和4 m/s2;根据1.3节方法,将调幅后的地震波分别输入到各组模型中,则每组模型地震易损性分析的样本容量均为90个.

我国规范《GB/T24335-2009建(构)建物地震破坏等级划分》中,将结构的破坏等级分为5个区间:基本完好、轻微破坏、中等破坏、严重破坏和倒塌. 划分这5个破坏等级区间的限值点称为性能指标,分别为轻微破坏LS1、中等破坏LS2、严重破坏LS3和倒塌LS4. 美国HAZUS99系列损伤评估准则[41]采用层间最大位移角θmax作为性能指标,这是国内学者在分析房屋的地震易损性时最常使用的性能指标. 我国建筑抗震设计规范[40]将钢筋混凝土框架结构的4个性能水平限值点对应的层间最大位移角分别取值为:LS1=1/550,LS2=1/250,LS3=1/120,LS4=1/60.

将5组模型共450个样本进行时程分析,得到各组模型在不同地震峰值、不同地震波输入下的层间最大位移角,李刚[42]认为在相同的PGA下层间最大位移角服从对数正态分布:

$\ln {\theta _{\max }} \sim {\rm{N}}({\mu _{\ln {\theta _{\max }}}},\sigma _{\ln {\theta _{\max }}}^2)$

式中: ${\mu _{\ln {\theta _{\max }}}}$为层间最大位移角的对数均值, ${\sigma _{\ln {\theta _{\max }}}}$为层间最大位移角的对数标准差.

结构的破坏概率函数,即在某一特定的PGA下层间最大位移角超过某一性能水平限值的概率 ${P}{(}{\theta _{\max }}|{{\rm PGA} > {\rm L}}{{{\rm S}}_{i}})$可以表示为

${P}{(}{\theta _{\max }}|{{\rm PGA} > {\rm L}}{{{\rm S}}_{i}}) = 1 - \varPhi \left( {\frac{{\ln ({{\rm L}}{{\rm S}_{i}}) - {\mu _{\ln {\theta _{\max }}|{{\rm PGA}}}}}}{{\sigma _{\ln {\theta _{\max }}|{{\rm PGA}}}^2}}} \right).$

式中: ${\mu _{\ln {\theta _{\max }}|{{\rm PGA}}}}$${\sigma _{\ln {\theta _{\max }}|{{\rm PGA}}}}$分别为加速度峰值为PGA时层间最大位移角的对数均值与对数标准差,由时程分析结果计算统计得到;LSi为对应结构的4个性能水平限值,i=1~4;Φ( )为标准正态分布函数.

通过以上计算,可得各模型在不同PGA下的超越概率,将其描绘在以PGA为横坐标、超越概率为纵坐标的坐标系中,并依次连接,绘得对应的地震易损性曲线.

3.2. 考虑河谷地形影响的地震易损性曲线

图9给出5组模型的地震易损性曲线. 图中,P为超越概率. 如图9(a)~(d)所示,分别对应轻微破坏、中等破坏、严重破坏和倒塌性能水平.

图 9

图 9   考虑局部河谷地形的结构地震易损性曲线

Fig.9   Seismic fragility curves considering local valley terrain


图9(a)可知,在轻微破坏性能水平下,5组模型的易损性曲线随PGA变化的趋势基本一致,均表现为在PGA较小时超越概率增长较迅速,在PGA=1 m/s2处易损性曲线出现拐点,随后超越概率增长变缓并逐渐趋近于1. 不考虑地基模型的超越概率整体上小于其他2种模型,说明不考虑土-结动力相互作用,会低估轻微破坏性能水平下的超越概率. 河谷地基的超越概率整体上大于均质地基,可见位于河谷地基上的房屋结构更易发生破坏,河谷地基x=0、30、60 m处的易损性相差不大.

图9(b)~(d)所示的3组曲线具有类似的特征:在中等破坏、严重破坏和倒塌性能水平下,超越概率在PGA较小时随着PGA的增加缓慢增长,当曲线出现第1个拐点(对应3种性能水平拐点位置分别为:PGA=0.5、1、2 m/s2处)后超越概率增长加快,以中等破坏性能水平为例,曲线在PGA=2 m/s2处出现第2个拐点,随后超越概率增长变慢并逐渐趋近于1. 对于严重破坏和倒塌性能水平,本文虽然未绘出完整的易损性曲线,但可以推测随着PGA的增大,超越概率的变化趋势与中等破坏性能水平类似. 除此之外,还能从图9(b)~(d)中得到与图9(a)类似的结论:在考虑土-结动力相互作用和局部河谷地形的影响后,房屋结构在相应性能水平下的超越概率有一定程度的增大. 房屋结构在性能水平较低的中等破坏、严重破坏和倒塌状态下,河谷地基不同位置处易损性的差异较轻微破坏明显增大.

为了定量分析土-结动力相互作用与局部河谷地形对结构易损性的影响,不考虑土-结相互作用、河谷地基x=0、30、60 m处结构的超越概率分别与位于均质地基上结构的超越概率作差比较,得到相应的超越概率增加量ΔP,如图10所示. 可以看出,在各性能水平下,不考虑地基影响时,超越概率减小,降幅最大为−7.9%,发生在轻微破坏性能水平下. 河谷地基3处超越概率均增大,最大增幅为12.1%,发生在严重破坏性能水平下的x=0处. 在相同的性能水平下,3条增幅曲线的变化趋势一致. 小震时,三者重合,表明位于河谷地基不同位置处的增加量一致,当PGA增大时,不同位置处的增加量不同,以轻微破坏性能水平为例,位于x=0、30、60 m处的超越概率最大可以增加9.0%、7.8%、5.9%,即结构离河谷中心越近,越容易破坏. 在其他性能水平下超越概率增加量表现出类似的规律. 这是由于在所研究的工况中,大多数地震波的聚焦点出现在河谷中心处,位于河谷中心附近的结构将承受更大的地震波作用,导致在相同条件下结构距河谷中心越近、损伤程度越高,发生破坏的概率越大. 在土-结动力相互作用与局部河谷地形的影响下,结构在不同性能水平下超越概率增加量峰值所对应的PGA从小到大分别为: 轻微损伤、中等损伤、严重损伤、倒塌.

图 10

图 10   河谷地形影响的超越概率增加量

Fig.10   Exceedance probability increase of local valley terrain


4. 结 语

杭嘉湖地区为冲积平原区,土体强度较低,地基具有较典型的局部河谷、盆地特征. 本文结合杭嘉湖地区的地形特点、地基条件及抗震设防区划要求,建立考虑土-结动力相互作用的框架结构数值模型,在所选取的18条地震波作用下,对比不考虑地基、均质地基、河谷地基条件下房屋结构的动力响应和地震易损性,分析结构位于河谷地基不同位置处的响应差异.

不考虑土-结动力的相互作用,多会低估房屋结构的结构变形,对应的易损性性能水平超越概率偏低,最大偏差为7.9%,结果偏于危险. 局部河谷地形独特的聚焦效应会增大房屋结构的结构变形,加剧房屋结构震害. 河谷地基中不同位置处房屋结构的结构变形明显不同. 本文所研究的工况下,房屋结构距河谷中心越近,结构变形越大、相应性能水平的超越概率的增量越大,最大增量为12.1%. 对位于河谷地基上的房屋结构进行易损性分析时,应根据相应地质勘探资料计算地震波聚焦点范围,考虑结构与聚焦点的相对位置对房屋结构易损性的影响.

参考文献

朱健. 结构动力学原理与地震易损性分析[M]. 北京: 科学出版社, 2013.

[本文引用: 1]

郭迅

汶川大地震震害特点与成因分析

[J]. 地震工程与工程振动, 2009, 29 (6): 74- 87

[本文引用: 1]

GUO Xun

Characteristics and mechanism analysis of the Great Wen Chuan Earthquake

[J]. Journal of Earthquake Engineering and Engineering Vibration, 2009, 29 (6): 74- 87

[本文引用: 1]

尚守平, 贺志文, 王海东, 等

上部结构与地基相对刚度比对土-结构体系基频影响试验研究

[J]. 地震工程与工程振动, 2008, 28 (05): 94- 101

[本文引用: 1]

SHANG Shou-ping, HE Zhi-wen, WANG Hai-dong, et al

Experimental investigation on the effect of the relative stiffness ratio between superstructure and ground soil on the fundamental frequency of soil-structure system

[J]. Journal of Earthquake Engineering and Engineering Vibration, 2008, 28 (05): 94- 101

[本文引用: 1]

PANZERA F, LOMBARDO G, MONACO C, et al

Seismic site effects observed on sediments and basaltic lavas outcropping in a test site of Catania, Italy

[J]. Natural Hazards, 2015, 79 (1): 1- 27

[本文引用: 1]

TAFAZZOLI N, BAZIAR M H. Evaluation of basin effect on ground motion [C] // 4th International Conference on Geotechnical Earthquake Engineering. Thessaloniki, Greece: Aristotle University of Thessaloniki Greece, 2007: 1323.

[本文引用: 1]

金峰, 张楚汉

半椭圆形河谷上沉积层地震响应研究

[J]. 清华大学学报: 自然科学版, 1993, 33 (05): 23- 30

[本文引用: 1]

JIN Feng, ZHANG Chu-han

Seismic response of sedimentary layers in semi-oval valley

[J]. Journal of Tsinghua University: Science and Technology, 1993, 33 (05): 23- 30

[本文引用: 1]

尤红兵, 赵凤新, 李方杰

层状场地中任意形状沉积河谷对平面SH波的散射

[J]. 防灾减灾工程学报, 2008, 28 (3): 313- 318

[本文引用: 1]

YOU Hong-bing, ZHAO Feng-xin, LI Fang-jie

Scattering of planar SH waves by arbitrary shaped sediments in a layered field

[J]. Journal of Disaster Prevention and Mitigation Engineering, 2008, 28 (3): 313- 318

[本文引用: 1]

刘晶波, 刘阳冰, 闫秋实, 等

基于性能的方钢管混凝土框架结构地震易损性分析

[J]. 土木工程学报, 2010, 43 (02): 39- 47

[本文引用: 1]

LIU Jing-bo, LIU Yang-bing, YAN Qiu-shi, et al

Performance-based seismic fragility analysis of CFST frame structures

[J]. China Civil Engineering Journal, 2010, 43 (02): 39- 47

[本文引用: 1]

刘阳冰, 刘晶波. 基于性能的组合结构的地震易损性分析[C] // 全国土木工程研究生学术论坛. 北京: 中国土木工程学会, 2008.

LIU Yang-bing, LIU Jing-bo. Seismic fragility analysis of composite structures based on performance [C] // National Civil Engineering Forum for Graduate Students. Beijing: CSCE, 2008.

CELIK O C, ELLINGWOOD B R

Seismic fragilities for non-ductile reinforced concrete frames: role of aleatoric and epistemic uncertainties

[J]. Structural Safety, 2010, 32 (1): 1- 12

PARK J, KIM J

Fragility analysis of steel moment frames with various seismic connections subjected to sudden loss of a column

[J]. Engineering Structures, 2010, 32 (6): 1547- 1555

[本文引用: 1]

葛立勇

新版《中国地震波参数区划图》简介及对浙江省的影响

[J]. 浙江建筑, 2017, 34 (03): 26- 28

[本文引用: 1]

GE Li-yong

Influence of the brief introduction of the new edition of《Zoning Map of the Seismic Parameters of China》

[J]. Zhejiang Construction, 2017, 34 (03): 26- 28

[本文引用: 1]

施炜, 叶列平, 陆新征, 等

不同抗震设防RC框架结构抗倒塌能力的研究

[J]. 工程力学, 2011, 28 (3): 41- 48

[本文引用: 1]

SHI Wei, YE Lie-ping, LU Xin-zheng, et al

Study on the collapse-resistant capacity of RC frames with different seismic fortification levels

[J]. Engineering Mechanics, 2011, 28 (3): 41- 48

[本文引用: 1]

陆新征, 叶列平, 缪志伟. 建筑抗震弹塑性分析[M]. 北京: 中国建筑工业出版社, 2015.

[本文引用: 1]

高向玲, 张业树, 李杰

基于ABAQUS梁单元的钢筋混凝土框架结构数值模拟

[J]. 结构工程师, 2013, 29 (06): 19- 26

[本文引用: 2]

GAO Xiang-ling, ZHANG Ye-shu, LI Jie

Numerical simulation of a RC frame based on abaqus beam elements

[J]. Structural Engineers, 2013, 29 (06): 19- 26

[本文引用: 2]

张亚旭, 王修信, 庄海洋

土-桩-框架结构非线性相互作用的精细数值模型及其验证

[J]. 防灾减灾工程学报, 2010, 30 (5): 558- 566

[本文引用: 1]

ZHANG Ya-xu, WANG Xiu-xin, ZHUANG Hai-yang

The fine numerical model of soil-pile-frame structure nonlinear interaction and its verification

[J]. Journal of Disaster Prevention and Mitigation Engineering, 2010, 30 (5): 558- 566

[本文引用: 1]

褚浩存

嵌固部位柱钢筋的施工工艺

[J]. 建筑工人, 2016, 37 (04): 4- 5

[本文引用: 1]

CHU Hao-cun

Construction process of reinforced steel bar embedded in site

[J]. Construction Workers, 2016, 37 (04): 4- 5

[本文引用: 1]

庄海洋, 于旭. 土-桩-隔震结构动力相互作用[M]. 北京: 中国建筑工业出版社, 2016.

[本文引用: 1]

蔡伟忠, 周金龙, 张远芳, 等

嘉兴市软弱地基工程特性及桩基持力层评价

[J]. 新疆农业大学学报, 2011, 34 (4): 353- 357

[本文引用: 1]

CAI Wei-zhong, ZHOU Jin-long, ZHANG Yuan-fang, et al

Evaluation on features of soft foundation engineering andbearing layer of pile in Jiaxing city

[J]. Journal of Xinjiang Agricultural University, 2011, 34 (4): 353- 357

[本文引用: 1]

费康, 张建伟. ABAQUS在岩土工程中的应用[M]. 北京: 中国水利水电出版社, 2010.

[本文引用: 1]

HUDSON M, IDRISS I M, BEIKAE M. User manual for QUAD4m: a computer program to evaluate the seismic response of soil structures using finite element procedures and incorporating a compliant base [D]. Berkeley: University of California, 1994.

[本文引用: 1]

NEJATI H R, AHMADI M, HASHEMOLHOSSEINI H

Numerical analysis of ground surface vibration induced by underground train movement

[J]. Tunnelling and Underground Space Technology incorporating Trenchless Technology Research, 2012, 29 (5): 1- 9

[本文引用: 1]

黄宗明, 白绍良, 赖明

结构地震反应时程分析中的阻尼问题评述

[J]. 地震工程与工程振动, 1996, (2): 95- 105

[本文引用: 1]

HUANG Zong-ming, BAI Shao-liang, LAI Ming

Comments on damping in structural seismic response time analysis

[J]. Journal of Earthquake Engineering and Engineering Vibration, 1996, (2): 95- 105

[本文引用: 1]

邹德高, 徐斌, 孔宪京

瑞利阻尼系数确定方法对高土石坝地震反应的影响研究

[J]. 岩土力学, 2011, 32 (03): 797- 803

[本文引用: 1]

ZOU De-gao, XU Bin, KONG Xian-jing

Study of influence of different methods for calculating Rayleigh damping coefficient on high earth-rock dam seismic response

[J]. Rock and Soil Mechanics, 2011, 32 (03): 797- 803

[本文引用: 1]

李哲, 王贡献, 胡勇, 等

改进的瑞利阻尼系数计算方法在岸桥结构地震反应分析中的应用

[J]. 华南理工大学学报: 自然科学版, 2015, 43 (06): 103- 109

[本文引用: 1]

LI Zhe, WANG Gong-xian, HU Yong, et al

Application of improved Rayleigh damping coefficient calculation method in seismic response analysis of shore bridge structure

[J]. Journal of South China University of Technology: Science and Technology, 2015, 43 (06): 103- 109

[本文引用: 1]

李平, 刘红帅, 薄景山, 等

汶川MS8. 0地震河谷地形对汉源县城高烈度异常的影响

[J]. 地球物理学报, 2016, 1 (1): 174- 184

[本文引用: 1]

LI Ping, LIU Hong-shuai, BO Jing-shan, et al

Effects of river valley topography on anomalously high intensity in the Hanyuan town during the Wenchuan

[J]. Chinese Journal of Geophysics, 2016, 1 (1): 174- 184

[本文引用: 1]

杨文采

波动方程的逆散射与迭代线性化反演

[J]. 石油物探, 1995, 34 (03): 114- 121

[本文引用: 2]

YANG Wen-cai

Inverse scattering of wave equations and iterative linearization inversion

[J]. Geophysical Prospecting for Petroleum, 1995, 34 (03): 114- 121

[本文引用: 2]

戚玉亮, 大塚久哲

ABAQUS动力无限元人工边界研究

[J]. 岩土力学, 2014, 35 (10): 3007- 3012

[本文引用: 1]

QI Yu-liang, HISANORI Otsuka

Study of ABAQUS dynamic infinite element artificial boundary

[J]. Rock and Soil Mechanics, 2014, 35 (10): 3007- 3012

[本文引用: 1]

UNGLESS R F. An infinite finite element [D]. Vancouver: University of British Columbia, 1973.

[本文引用: 1]

阳栋, 王志亮

基于无限元和波场分离法的地震响应数值分析

[J]. 同济大学学报: 自然科学版, 2012, 40 (8): 1129- 1134

[本文引用: 1]

YANG Dong, WANG Zhi-liang

Numerical analysis of seismic response based on infinite element and wave field separation method

[J]. Journal of Tongji University: Natural Science, 2012, 40 (8): 1129- 1134

[本文引用: 1]

ZHAO Guo-chen, XU Long-jun, XIE Li-li

On near-fault ground motion characteristics through multi-scale method

[J]. Chinese Journal of Geophysics, 2013, 56 (6): 742- 753

[本文引用: 1]

Applied Technology Council

Quantification of building seismic performance factors

[J]. California: Federal Emergency Management Agency, 2008,

[本文引用: 1]

顾斌, 丁政

苏浙沪及邻区深部构造及其与地震活动性的关系

[J]. 地震学刊, 1987, (3): 29- 37

[本文引用: 1]

GU Bin, DING Zheng

Deep structures and their relation to earthquake activity in Suzhizi and neighboring areas

[J]. Journal of Seismology, 1987, (3): 29- 37

[本文引用: 1]

田石柱, 张常明, 林元铮

考虑桩-土-结构相互作用下斜拉桥桥塔地震反应模拟

[J]. 工程抗震与加固改造, 2015, 37 (06): 71- 77

[本文引用: 1]

TIAN Shi-zhu, ZHANG Chang-ming, LIN Yuan-zheng

Seismic response simulation of cable-stayed bridge tower considering pile-soil-structure interaction

[J]. Earthquake Resistant Engineering and Retrofitting, 2015, 37 (06): 71- 77

[本文引用: 1]

CHOI Y, STEWART J P, GRAVES R W

Empirical model for basin effects accounts for basin depth and source location

[J]. Bulletin of the Seismological Society of America, 2005, 95 (4): 1412- 1427

[本文引用: 1]

王力, 周叮, 刘伟庆, 等

场地盆地效应的振动台试验

[J]. 南京工业大学学报: 自然科学版, 2014, 36 (02): 107- 111

[本文引用: 1]

WANG Li, ZHOU Ding, LIU Wei-qing, et al

Shaking table test of site basin effect

[J]. Journal of Nanjing Tech University: Natural Science Edition, 2014, 36 (02): 107- 111

[本文引用: 1]

袁晓铭, 李雨润, 孙锐

圆弧状沉积盆地与软土单覆盖层出平面地表运动对比

[J]. 地震工程与工程振动, 2002, 22 (4): 16- 21

[本文引用: 1]

YUAN Xiao-ming, LI Yu-run, SUN Rui

Comparison of out-of-plane surface ground motion between a circular-arc alluvial valley and a single overburden soft layer

[J]. Earthquake Engineering and Engineering Vibration, 2002, 22 (4): 16- 21

[本文引用: 1]

SEMBLAT J F, PARARA E, KHAM M, et al. Site effects: basin geometry vs soil layering [C] // Joint Conference of the 11th International Conference on Soil Dynamics and Earthquake Engineering and the 3rd International Conference on Earthquake Geotechnical Engineering. Berkeley: University of California, 2004.

[本文引用: 1]

孙柏涛, 张桂欣

汶川8. 0级地震中各类建筑结构地震易损性统计分析

[J]. 土木工程学报, 2012, 45 (02): 26- 30

[本文引用: 1]

SUN Bo-tao, ZHANG Gui-xin

Statistical analysis of the seismic fragility of various types of building structures in WenchuanM8. 0 earthquake

[J]. China Civil Engineering Journal, 2012, 45 (02): 26- 30

[本文引用: 1]

GB50011-2010, 中华人民共和国国家标准建筑抗震设计规范[S]. 北京: 中国建筑工业出版社, 2008.

[本文引用: 2]

HAZUS-MH MR1. Multi-hazard loss estimation methodology-earthquake model [S]. Washington D. C. : Department of Homeland Security, FEMA, 2003.

[本文引用: 1]

李刚. 基于性能的结构抗震设计[M]. 北京: 科学出版社, 2004.

[本文引用: 1]

/