• ISSN 1673-5722
  • CN 11-5429/P

基于可靠度与易损性的地震动输入修正方法

刘浪 苗思雨 夏永庆 刘霍义 闫薄全

刘浪,苗思雨,夏永庆,刘霍义,闫薄全,2024. 基于可靠度与易损性的地震动输入修正方法. 震灾防御技术,19(2):378−386. doi:10.11899/zzfy20240217. doi: 10.11899/zzfy20240217
引用本文: 刘浪,苗思雨,夏永庆,刘霍义,闫薄全,2024. 基于可靠度与易损性的地震动输入修正方法. 震灾防御技术,19(2):378−386. doi:10.11899/zzfy20240217. doi: 10.11899/zzfy20240217
Liu Lang, Miao Siyu, Xia Yongqing, Liu Huoyi, Yan Boquan. Calibration of Seismic Input Parameter Based on Reliability and Fragility[J]. Technology for Earthquake Disaster Prevention, 2024, 19(2): 378-386. doi: 10.11899/zzfy20240217
Citation: Liu Lang, Miao Siyu, Xia Yongqing, Liu Huoyi, Yan Boquan. Calibration of Seismic Input Parameter Based on Reliability and Fragility[J]. Technology for Earthquake Disaster Prevention, 2024, 19(2): 378-386. doi: 10.11899/zzfy20240217

基于可靠度与易损性的地震动输入修正方法

doi: 10.11899/zzfy20240217
基金项目: 国家自然科学基金资助项目(51708069);重庆市基础研究与前沿探索项目(cstc2018 jcyjA2535)
详细信息
    作者简介:

    刘浪,女,生于1985年。副教授。主要从事随机荷载与桥梁可靠度、桥梁抗震方面的研究。 E-mail:yilupaolai2008@163.com

Calibration of Seismic Input Parameter Based on Reliability and Fragility

  • 摘要: 在基于性能的地震工程理念与经典可靠度理论框架下,分析了地震易损性与结构可靠度的共同特性,阐释了将两者相结合的理论基础与概率方法;以一座三跨连续梁桥为例,基于OpenSees平台建立桥梁有限元模型,通过非线性动力时程分析对结构进行概率地震需求及易损性分析,基于易损性曲线考虑3种抗震目标对桥墩进行抗震设计,并对设计墩柱进行恒载与车辆活载作用下的可靠度分析;考虑不同地震重现期,根据泊松随机过程计算各类地震作用下的目标可靠度,基于概率方法对结构的初始可靠度进行修正,进一步结合地震易损性,确定不同抗震目标下的地震动水平,提出基于可靠度与易损性的地震动输入参数确定方法。结果表明,1000年、475年和50年地震重现期下,修正PGA均值分别为0.25 g、0.24 g和0.22 g,高于仅考虑易损性时的0.21 g。本方法将地震易损性与可靠度相结合,为结构不同层次的抗震设计提供了思路与方法。
  • 近年来,基于性能的地震工程(Performance-Based Earthquake Engineering,PBEE)理念(Deierlein等,2003Porter,2003Moehle等,2004)受到了地震工程与工程抗震领域的广泛关注,其主要思想是将结构的抗震设计从简单的、确定性方法向多样化、多层次功能目标以及可修复能力等进行过渡(Pathak等,2020)。在PBEE框架下,国内外学者开展了大量地震危险性分析、结构地震响应分析、结构损伤分析(结构地震响应和结构损伤分析亦统称为地震易损性分析)及地震灾害损失评估工作(Dong等,2014袁万城等,2016冯莉等,2020于晓辉等,2022),相关成果为确定结构的薄弱环节、抗震加固、风险评估等提供了重要支撑(Muntasir Billah等,2015李宏男等,2018)。

    另一方面,基于可靠度理论进行工程结构设计与评估已成为主要趋势(李杰,2018),其将结构自身的抗力与承受的外界荷载视为随机变量,建立功能函数分析两者的相对关系,并要求最终满足某一给定的可靠度指标。在此地震作用被视为一种随机作用,未明确地震作用下结构应当达到的可靠度水平(Yoon等,2019)。从本质上看,易损性与可靠度均从概率意义上刻画了结构完成某种既定目标的能力,并基于概率方法定量评定了工程结构的某种性能,如果通过概率方法将两者结合,在一定可靠度水平下进行结构易损性分析,则可为基于性能的结构设计提供更加明确的目标。

    为此,本文首先阐释易损性与可靠度相结合的理论基础,然后通过地震危险性(重现期)设定地震作用下的目标可靠度,基于概率方法对结构原可靠度进行修正,再利用结构易损性分析确定不同抗震目标下的地震动水平,最终提出基于可靠度与易损性的地震动输入参数确定方法,并以一座三跨连续梁桥为例进行应用说明。

    地震易损性作为PBEE的重要组成部分,表达了在不同强度的地震荷载工况下,结构发生各种程度破坏的概率(陆新征等,2010),建立了地震动强度与破坏程度间的关系,实现了对结构抗震性能的量化估计。但另一方面,我国当前所采用的工程结构可靠度设计方法中,尚未明确提出地震作用下的结构可靠度要求。如上文所述,地震易损性从概率意义上刻画了结构在不同地震设防水准下结构完成预定抗震目标的能力,表达了经典可靠度理论中“荷载”与“抗力”的关系(吕大刚等,2006),即地震危险性相当于结构的“广义荷载”,地震易损性则相当于“广义抗力”。从这个意义上讲,我们可以通过地震重现期(代表地震危险性)为结构设定地震作用下的“目标可靠度”,利用该目标可靠度对结构所具备的初始可靠度进行修正,再将修正的参数与地震易损性结合,即可得到不同地震危险性下基于易损性与可靠度分析的地震动输入参数。该方法基本原理如图1所示。

    图 1  基于易损性与可靠度的地震动修正流程
    Figure 1.  Flowchart of the proposed calibration method

    通常认为结构易损性曲线服从双参数的对数正态累积分布(Karim等,2001),可表示为:

    $$ {G}_{{\mathrm{h}}}\left(Y\geqslant D|PGA=x\right)=\Phi \left[\frac{\mathrm{ln}\left(\dfrac{x}{a}\right)}{b}\right] $$ (1)

    式中,$ \mathrm{\Phi }\left[\cdot \right] $为标准正态累计分布函数;$ a $$ b $分别为易损性函数的均值与标准差。

    $ a $$ b $的值可由“云图法”得到(Shome等,1998),即结构的地震需求服从对数正态分布,且工程需求参数的中值$ \overline {{M}_{\mathrm{E}\mathrm{D}\mathrm{P}}} $与地震动强度参数$ {T}_{\mathrm{I}\mathrm{M}} $之间满足如下指数关系:

    $$ \overline {{M}_{\mathrm{E}\mathrm{D}\mathrm{P}}}=\alpha \cdot {T}_{\mathrm{I}\mathrm{M}}^{\beta } $$ (2)

    式中,$ \alpha 、\beta $为估计参数,可以通过最小二乘回归分析得到。将式(2)改写为式(3)所示线性回归函数,即结构概率性地震需求模型:

    $$ \mathrm{ln}\left(\overline {\overline {{M}_{\mathrm{E}\mathrm{D}\mathrm{P}}}}\right)=\beta \cdot \mathrm{ln}\left({T}_{\mathrm{I}\mathrm{M}}\right)+\mathrm{ln}\left(\alpha \right) $$ (3)

    由于在不同地震下地震需求模型的对数标准差$ b $不同,根据Mackie等(2003)可假定在全部地震动PGA范围内对数标准差为一个定值,采用式(4)进行计算:

    $$ {b}_{{\mathrm{EDP}}|{\mathrm{IM}}}=\sqrt{\frac{\displaystyle\sum _{i-1}^{N}{\left[\mathrm{ln}\left({D}_{i}\right)-\mathrm{ln}\left(\overline {{M}_{\mathrm{E}\mathrm{D}\mathrm{P}}}\right)\right]}^{2}}{N-2}}$$ (4)

    式中,$ {b}_{{\mathrm{EDP}}|{\mathrm{IM}}} $表示在给定PGA范围内地震响应的条件对数标准差;$ {D}_{i} $为结构在第i条地震波作用下的需求响应峰值;N为时程分析次数。因此,地震易损性函数可表示为:

    $$ {G}_{{\mathrm{h}}}=\Phi \left(\dfrac{\mathrm{ln}\dfrac{\overline {{M}_{\mathrm{E}\mathrm{D}\mathrm{P}}}}{{T}_{\mathrm{c}}}}{\sqrt{{b}_{{\mathrm{EDP}}|{\mathrm{IM}}}^{2}+{b}_{{\mathrm{c}}}^{2}}}\right)$$ (5)

    式中,$ {T}_{\mathrm{c}} $$ {b}_{\mathrm{c}} $表示指定损伤极限状态(抗震能力)的均值和对数标准差。

    根据经典可靠度理论,假设抗力与荷载随机变量RS服从标准正态分布且相互独立,功能函数Z=R-S亦服从标准正态分布,则Z的概率密度函数为表示为:

    $$ \varnothing \left(Z\right)=\frac{1}{\sqrt{2\text{π} }}{\mathrm{e}}^{-\tfrac{{Z}^{2}}{2}} $$ (6)

    结构的失效概率与初始可靠度β满足:

    $$ {P}_{{\mathrm{f}}}=1-\Phi \left(\beta \right)={\int }_{\beta }^{+\infty }\frac{1}{\sqrt{2\text{π} }}{\mathrm{e}}^{-\tfrac{{Z}^{2}}{2}}\mathrm{d}Z$$ (7)

    其中,抗力随机变量R主要考虑构件几何尺寸与材料力学性能的变异,荷载随机变量S则为外部随机荷载产生的效应,再采用JC法即可计算出构件所具备的初始可靠度水准β

    地震学中常用泊松随机过程来描述t年内某地区发生地震的次数(Ghosh等,2015),则发生n次地震的概率Pn)可表示为:

    $$P\left(n\right)=\frac{{\left(vt\right)}^{n}\cdot {\mathrm{e}}^{\left(-vt\right)}}{n!} $$ (8)

    那么,t年内该地区不发生地震的概率为:

    $$ P\left(0\right)=\frac{{\left(vt\right)}^{0}\cdot {\mathrm{e}}^{\left(-vt\right)}}{0!}={\mathrm{e}}^{\left(-vt\right)} $$ (9)

    该地区在t年内至少发生一次地震的概率,即超越概率为:

    $$ F\left(t\right)=1-P\left(0\right)=1-{\mathrm{e}}^{\left(-vt\right)} $$ (10)

    以上各式中$ v $为地震年平均发生概率,其与重现期$ {T}_{0} $互为倒数。因此重现期T0与超越概率满足:

    $$ {T}_{0}=\frac{1}{v}=\frac{-t}{\mathrm{ln}\left(1-F\left(t\right)\right)} $$ (11)

    考虑到地震发生的偶然性,可以适当降低其可靠度水准。为便于后续对比不同重现期下的目标可靠度,本文将结构失效概率设为重现期的倒数,因此地震作用下的目标可靠度$ {\beta }{'} $可表达为:

    $$ {\beta }{{'}}={\mathrm{\Phi }}^{-1}\left(1-{P}'_{{\mathrm{f}}}\right)={\mathrm{\Phi }}^{-1}\left(1-\frac{1}{{T}_{0}}\right) $$ (12)

    设地震的超越概率密度函数为$ f\left(t\right) $$ f\left(t\right) $服从泊松分布,则在重现期T0和设计基准期100年条件下,其与结构可靠度功能密度函数$ \varnothing \left(Z\right) $联合概率函数为:

    $$ F\left(t,Z\right)=P\left[\beta < z < +\infty ,0 < t < 100\right]\\ ={\int }_{\beta }^{+\infty }\frac{1}{\sqrt{2\text{π} }}{\mathrm{e}}^{-\tfrac{{Z}^{2}}{2}}\mathrm{d}Z{\int }_{0}^{{T}_{0}}f\left(t\right)\mathrm{d}t =P(f) \times (1-{{\mathrm{e}}}^{\left(-vt\right)}) $$ (13)

    由前文知,代入相关参数即可对公式(13)进行求解。由推导过程可以看出,联合概率函数$ F\left(t,Z\right) $正是利用地震重现期下的目标可靠度$ {\beta }{{'}} $对结构初始$ \beta $进行修正,由此计算出结构的最终可靠度水准$ {\beta }'' $

    定义修正系数$ \alpha $为考虑地震发生概率(重现期)的结构可靠度水准$\beta '' $与初始可靠度水准$ \beta $的比值,即$ {\alpha =\beta }''/\beta $。将$ \alpha $代入公式$ {u}'_{\mathrm{D}\mathrm{I}}=\alpha \cdot {u}_{\mathrm{D}\mathrm{I}} $中对结构原损伤指标进行修正,修正后的损伤指数$ {u}'_{\mathrm{D}\mathrm{I}} $所对应的地震动PGA即为一定地震危险性下基于地震概率需求与结构可靠度水准的地震动输入参数。

    现以山区公路常用的连续梁桥为例,对基于可靠度与易损性的地震动修正方法进行说明。本文依托实桥跨径为40 m+70 m+40 m,上部为单箱单室箱梁结构,主梁采用C50混凝土,下部矩形截面桥墩采用C30混凝土,墩高10 m,墩身纵筋和箍筋均采用HRB400钢筋,直径分别为Φ28 mm和Φ16 mm,间距0.5 m。场地设防烈度为VII度,桥址区场地条件为III类。

    基于OpenSees平台对主梁、桥墩、支座等重要构件进行模拟,考虑材料非线性、几何非线性和接触非线性,建立桥梁有限元模型,如图2所示。其中,主梁采用弹性梁柱单元Elastic Beam Column模拟,支座采用零长度单元Zero Length模拟,桥墩采用纤维截面Nonlinear Beam Column单元模拟,将截面离散为保护层混凝土、核心混凝土以及钢筋纤维。混凝土采用Conerete 02 Kent-Park本构模型,钢筋采用Steel 02本构模型,桥墩墩身纤维截面如图3所示。本例修正以结合墩柱构件的可靠度为主,故不考虑桩-土相互作用,将墩底固结。

    图 2  OpenSees桥梁有限元模型
    Figure 2.  Finite element model of the bridge in OpenSees
    图 3  墩柱纤维截面
    Figure 3.  Fiber section of the pier

    根据桥梁设计反应谱及场地条件,并适当考虑地震动峰值及频率成分涵盖范围,本文从美国太平洋地震中心PEER选取了6条基础地震动记录(表1),其地震加速度反应谱曲线、平均反应谱曲线与设计反应谱曲线对比如图4所示。从图中可以看出,所选地震动记录以高频成分为主,加速度反应谱与设计反应谱在短周期段吻合较好,且平均反应谱与设计反应谱基本吻合,能够反映本文连续梁桥(自振周期T=0.35 s)的地震响应特性,可采用增量动力IDA方法进行结构地震反应分析。通过IDA调幅获得120条PGA在0.03~0.6 g范围内的地震动,为提高加载效率并获得结构的最大地震反应,本文截取每条记录中PGA最大的30 s作为地震动持时。

    表 1  基础地震动记录
    Table 1.  Ground motions set as benchmark
    序号地震记录记录地点震中距/kmPGA/g
    1BorregoEl Centro Array #956.880.066
    2Kern CountyLA - Hollywood Stor FF114.620.042
    3San Fernando2516 Via Tejon PV55.20.026
    4San FernandoLB - Terminal Island58.990.029
    5San FernandoPort Hueneme68.840.027
    6San FernandoSan Juan Capistrano108.010.043
    下载: 导出CSV 
    | 显示表格
    图 4  地震反应谱
    Figure 4.  Seismic response spectrum

    采用Hwang等(2001)提出的桥墩顶部位移延性比作为地震损伤指标,即某地震水平下最大墩顶位移与受拉钢筋首次屈服时的墩顶位移之比。其破坏状态分为无破坏、轻微破坏、中等破坏、严重破坏和完全破坏5个等级,相应的损伤指标如表2所示。

    表 2  桥梁损伤状态描述及损伤指标
    Table 2.  Bridge damage states and their measurements by displacement ductility ratios
    损伤状态 具体描述 准则 位移延性比
    无破坏墩柱无明显裂缝, 钢筋无屈服$ 0 < \mu \leqslant {\mu }_{\mathrm{c}\mathrm{y}1} $$ 0 < \mu \leqslant 1 $
    轻微破坏墩柱表面出现明显裂缝,最外侧钢筋首次出现理论屈服$ {\mu }_{\mathrm{c}\mathrm{y}1} < \mu \leqslant {\mu }_{\mathrm{c}\mathrm{y}} $$ 1 < \mu \leqslant 1.12 $
    中等破坏表层混凝土部分脱落,墩柱产生非线性变形,墩底塑性铰开始形成$ {\mu }_{\mathrm{c}\mathrm{y}} < \mu \leqslant {\mu }_{\mathrm{c}4} $$ 1.12 < \mu \leqslant 1.97 $
    严重破坏塑性铰完全形成,保护层混凝土全部剥落,核心混凝土部分开裂,纵筋大量屈服$ {\mu }_{\mathrm{c}4} < \mu \leqslant {\mu }_{\mathrm{c}\mathrm{m}\mathrm{a}\mathrm{x}} $$ 1.97 < \mu \leqslant 4.97 $
    完全破坏核心混凝土压碎,箍筋断裂$ {\mu }_{\mathrm{c}\mathrm{m}\mathrm{a}\mathrm{x}}\leqslant \mu $$ 4.97\leqslant \mu $
    下载: 导出CSV 
    | 显示表格

    用Xtract软件对墩底截面进行弯矩曲率分析,计算得到首次屈服曲率为$ 0.791\times {10}^{-3} $、等效屈服曲率为$ 0.8837\times {10}^{-3} $、混凝土压应变为0.004时的曲率为$ 3.15\times {10}^{-3} $。再对结构进行地震动作用下的非线性时程分析,最终得到桥墩易损性曲线如图5所示。

    图 5  地震易损性曲线
    Figure 5.  Seismic fragility curve

    参考JTG/T 2231—01—2020《公路桥梁抗震设计规范》的设防目标及以往桥梁震害调查资料,本文分别将80%概率出现轻微损伤、60%概率出现中度损伤和40%概率出现重度损伤作为桥墩的3种抗震性能目标,相应的地震动PGA分别为0.21 g、0.18 g、0.23 g,取其平均值0.21 g作为基于桥梁易损性的地震动输入参数,对原桥墩截面重新进行抗震设计。以抗震规范E1地震作用为例,采用能力保护方法进行墩柱截面抗震设计。主墩墩顶地震作用为5029.34 kN,墩高10 m,墩底弯矩为50293.44 kN·m。纵向配筋细节如图6所示。与原设计截面相比,基于易损性的抗震设计考虑了地震随机性及不同抗震性能,增大了纵桥向的截面钢筋配筋率。在钢筋和混凝土材料类别不变的情况下,主筋由62根增加至74根,钢筋间距由12.5 cm减小至11 cm。

    图 6  墩柱横截面配筋细节
    Figure 6.  Details of reinforcements in the cross section of the pier

    将截面尺寸、纵向钢筋及混凝土强度等几何参数和材料参数视作随机变量(张明等,2015),采用蒙特卡罗法生成5000组随机样本数据,计算墩底抗力R的均值与标准差;然后考虑上部结构传递下来的恒载和活载,计算墩柱所承受的荷载效应S。其中,恒载考虑结构自重加二期恒载,活载按照JTG D60—2015《公路桥涵设计通用规范》车道荷载模型,考虑车道横向折减,进行3个车道布载计算,计算恒载、活载作用下墩底的弯矩标准值,再参考JTG 2120—2020《公路工程结构可靠性设计统一标准》荷载效应建议系数值,计算墩柱荷载效应S的均值和标准差。RS的统计特征参数值如表3所示。最后采用JC法计算得到墩柱的初始可靠度$ \beta =3.3039 $

    表 3  抗力和荷载效应的统计特征参数
    Table 3.  Statistics parameters for resistance and load effect
    随机变量类别 均值/(kN·m) 标准差
    抗力 $ {\mu }_{{M}_{{\mathrm{R}}}}=48\;863.72 $ $ {\sigma }_{{M}_{{\mathrm{R}}}}=2\;855.02 $
    恒载 $ {\mu }_{{M}_{{\mathrm{G}}}}=19\;722.05 $ $ {\sigma }_{{M}_{{\mathrm{G}}}}=8\;492.84 $
    活载 $ {\mu }_{{M}_{{\mathrm{Q}}}}=1\;105.61 $ $ {\sigma }_{{M}_{{\mathrm{Q}}}}=173.55 $
    下载: 导出CSV 
    | 显示表格

    按照前文提出的理论计算方法,本文通过设定不同的地震重现期$ {T}_{0} $来考虑不同地震设防水准(即地震危险性)。同时,桥梁活荷载的设计基准期取为100年。以$ {T}_{0}=75 $年为例阐述该重现期下墩柱可靠度的修正过程,其他重现期的修正可靠度与修正系数计算结果如表4所示。

    表 4  不同地震重现期下的修正可靠度与修正系数
    Table 4.  Calibrated reliability index and calibration coefficient for various return periods
    重现期$ {T}_{0} $设计基准期t修正后$ {\beta }'' $初始$ \beta $修正系数$ \alpha $
    10001003.91413.30391.1847
    4751003.74413.30391.1332
    1001003.43043.30391.0383
    751003.38883.30391.0257
    501003.34453.30391.0123
    251003.30913.30391.0016
    101003.30393.30391.0000
    51003.30393.30391.0000
    下载: 导出CSV 
    | 显示表格

    $ {T}_{0}=75 $时,由式(12)可计算得到$ {\beta }{{'}} $=2.217。因此,在桥梁设计基准期100年内发生地震且墩柱的失效概率为:

    $$ {P}''_{{\mathrm{f}}}=F\left(t,Z\right)=\left(1-{{\mathrm{e}}}^{\left(-vt\right)}\right)\cdot {P}_{{\mathrm{f}}}=\left(1-{{\mathrm{e}}}^{\left(-100\times 1.33\times {10}^{-2}\right)}\right)\cdot {P}_{{\mathrm{f}}} $$ (14)

    根据式(7)计算得到墩柱初始可靠度$ \beta =3.3039 $,失效概率$ {P}_{{\mathrm{f}}}=4.77\times {10}^{-4} $,代入式(14)即可得到$ {P}_{{\mathrm{f}}}''=3.51\times {10}^{-4} $$ \mathrm{对}\mathrm{应}\mathrm{的}{\beta }''=3.3888 $

    利用表4中不同重现期对应的修正系数α,对结构损伤指标进行修正,再结合结构的易损性曲线,得到不同地震危险性下考虑结构可靠度水准的地震动PGA。通过前文的易损性计算可得,修正前$ {u}_{\mathrm{D}\mathrm{I}} $与PGA满足ln($ {u}_{\mathrm{D}\mathrm{I}} $)=1.3013 ln(PGA)+2.4567,经本文方法修正后其关系满足ln($ {u}'_{\mathrm{D}\mathrm{I}} $)=1.2939 ln(PGA)+2.375。限于篇幅,本文仅展示重现期为1000年、475年和50年时各参数的修正结果(表5表7)。由表可知,在文中的3种抗震性能目标下,经本文方法修正后考虑地震重现期的PGA均值分别为0.25 g、0.24 g和0.22 g,略高于仅考虑易损性时的0.21 g。3种抗震性能与地震重现期代表了不同的抗震设防目标,与传统的抗震规范地震动参数取值方法相比,本文提出的基于地震易损性与结构可靠度进行地震动参数修正的方法,能够将地震危险性、结构的抗震性能目标与可靠度需求紧密结合,为结构不同层次的抗震设计提供了思路与参考。

    表 5  重现期为1000年时的地震动参数修正值(${\boldsymbol{ \alpha }} $ = 1.1847)
    Table 5.  Calibrated peak ground accelerations for the return period of 1000 a(${\boldsymbol{ \alpha }} $ = 1.1847)
    损伤状态 修正前 修正后
    PGA/g $ {u}_{\mathrm{D}\mathrm{I}} $ $ {u}'_{\mathrm{D}\mathrm{I}} $ PGA/g
    80%轻微损伤 0.210 1.531 $ \alpha \cdot DI=1.814 $ 0.253
    60%中度损伤 0.183 1.280 $ \alpha \cdot DI=1.516 $ 0.220
    40%重度损伤 0.229 1.714 $ \alpha \cdot DI=2.030 $ 0.276
    下载: 导出CSV 
    | 显示表格
    表 6  重现期为475年时的地震动参数修正值($ {\boldsymbol{\alpha}} $ = 1.1332)
    Table 6.  Calibrated peak ground accelerations for the return period of 475 a(${\boldsymbol{ \alpha}} $ = 1.1332)
    损伤状态 修正前 修正后
    PGA/g $ {u}_{\mathrm{D}\mathrm{I}} $ $ {u}'_{\mathrm{D}\mathrm{I}} $ PGA/g
    80%轻微损伤 0.210 1.531 $ \alpha \cdot DI=1.735 $ 0.244
    60%中度损伤 0.183 1.280 $ \alpha \cdot DI=1.450 $ 0.213
    40%重度损伤 0.229 1.714 $ \alpha \cdot DI=1.942 $ 0.266
    下载: 导出CSV 
    | 显示表格
    表 7  重现期为50年时的地震动参数修正值(${\boldsymbol{ \alpha}} $ =1.0123)
    Table 7.  Calibrated peak ground accelerations for the return period of 50 a(${\boldsymbol{ \alpha }} $ = 1.0123)
    损伤状态 修正前 修正后
    PGA/g $ {u}_{\mathrm{D}\mathrm{I}} $ $ {u}'_{\mathrm{D}\mathrm{I}} $ PGA/g
    80%轻微损伤 0.210 1.531 $ \alpha \cdot DI=1.550 $ 0.224
    60%中度损伤 0.183 1.280 $ \alpha \cdot DI=1.296 $ 0.195
    40%重度损伤 0.229 1.714 $ \alpha \cdot DI=1.735 $ 0.244
    下载: 导出CSV 
    | 显示表格

    本文在基于性能的地震工程PBEE和经典可靠度理论框架下,阐释了结构易损性与可靠度相结合的可能性及其理论基础,通过设定地震重现期下的目标可靠度,基于概率方法对结构原可靠度进行修正,再结合地震易损性分析,确定不同抗震目标下的地震动水平,提出了基于可靠度与易损性的地震动输入参数确定方法,最后以一座三跨连续梁桥为例,对该方法进行应用说明,主要得到以下结论:

    (1)地震易损性与结构可靠度均从概率意义上刻画了结构完成某种既定目标的能力,通过概率方法将两者结合,在一定可靠度水平下进行结构易损性分析,可为基于性能的结构设计提供更加明确的目标。

    (2)本文提出的基于地震易损性与结构可靠度的地震动参数修正方法,将地震危险性、结构抗震性能目标、地震作用及一般荷载作用紧密结合,为结构不同层次的抗震设计提供了思路与方法。

    (3)在文中的3种抗震性能目标(80%轻微损伤、60%中度损伤和40%重度损伤)下,考虑地震重现期修正后的PGA均值分别为0.25 g、0.24 g和0.22 g,高于仅考虑易损性时的0.21 g

    (4)结构的地震易损性分析是基于性能的抗震设计中的重要一环,结构可靠度理论与其具有较多共同特点,将三者相互结合,是产生基于可靠度和性能的全生命周期抗震设计理论的新思路。

  • 图  1  基于易损性与可靠度的地震动修正流程

    Figure  1.  Flowchart of the proposed calibration method

    图  2  OpenSees桥梁有限元模型

    Figure  2.  Finite element model of the bridge in OpenSees

    图  3  墩柱纤维截面

    Figure  3.  Fiber section of the pier

    图  4  地震反应谱

    Figure  4.  Seismic response spectrum

    图  5  地震易损性曲线

    Figure  5.  Seismic fragility curve

    图  6  墩柱横截面配筋细节

    Figure  6.  Details of reinforcements in the cross section of the pier

    表  1  基础地震动记录

    Table  1.   Ground motions set as benchmark

    序号地震记录记录地点震中距/kmPGA/g
    1BorregoEl Centro Array #956.880.066
    2Kern CountyLA - Hollywood Stor FF114.620.042
    3San Fernando2516 Via Tejon PV55.20.026
    4San FernandoLB - Terminal Island58.990.029
    5San FernandoPort Hueneme68.840.027
    6San FernandoSan Juan Capistrano108.010.043
    下载: 导出CSV

    表  2  桥梁损伤状态描述及损伤指标

    Table  2.   Bridge damage states and their measurements by displacement ductility ratios

    损伤状态 具体描述 准则 位移延性比
    无破坏墩柱无明显裂缝, 钢筋无屈服$ 0 < \mu \leqslant {\mu }_{\mathrm{c}\mathrm{y}1} $$ 0 < \mu \leqslant 1 $
    轻微破坏墩柱表面出现明显裂缝,最外侧钢筋首次出现理论屈服$ {\mu }_{\mathrm{c}\mathrm{y}1} < \mu \leqslant {\mu }_{\mathrm{c}\mathrm{y}} $$ 1 < \mu \leqslant 1.12 $
    中等破坏表层混凝土部分脱落,墩柱产生非线性变形,墩底塑性铰开始形成$ {\mu }_{\mathrm{c}\mathrm{y}} < \mu \leqslant {\mu }_{\mathrm{c}4} $$ 1.12 < \mu \leqslant 1.97 $
    严重破坏塑性铰完全形成,保护层混凝土全部剥落,核心混凝土部分开裂,纵筋大量屈服$ {\mu }_{\mathrm{c}4} < \mu \leqslant {\mu }_{\mathrm{c}\mathrm{m}\mathrm{a}\mathrm{x}} $$ 1.97 < \mu \leqslant 4.97 $
    完全破坏核心混凝土压碎,箍筋断裂$ {\mu }_{\mathrm{c}\mathrm{m}\mathrm{a}\mathrm{x}}\leqslant \mu $$ 4.97\leqslant \mu $
    下载: 导出CSV

    表  3  抗力和荷载效应的统计特征参数

    Table  3.   Statistics parameters for resistance and load effect

    随机变量类别 均值/(kN·m) 标准差
    抗力 $ {\mu }_{{M}_{{\mathrm{R}}}}=48\;863.72 $ $ {\sigma }_{{M}_{{\mathrm{R}}}}=2\;855.02 $
    恒载 $ {\mu }_{{M}_{{\mathrm{G}}}}=19\;722.05 $ $ {\sigma }_{{M}_{{\mathrm{G}}}}=8\;492.84 $
    活载 $ {\mu }_{{M}_{{\mathrm{Q}}}}=1\;105.61 $ $ {\sigma }_{{M}_{{\mathrm{Q}}}}=173.55 $
    下载: 导出CSV

    表  4  不同地震重现期下的修正可靠度与修正系数

    Table  4.   Calibrated reliability index and calibration coefficient for various return periods

    重现期$ {T}_{0} $设计基准期t修正后$ {\beta }'' $初始$ \beta $修正系数$ \alpha $
    10001003.91413.30391.1847
    4751003.74413.30391.1332
    1001003.43043.30391.0383
    751003.38883.30391.0257
    501003.34453.30391.0123
    251003.30913.30391.0016
    101003.30393.30391.0000
    51003.30393.30391.0000
    下载: 导出CSV

    表  5  重现期为1000年时的地震动参数修正值(${\boldsymbol{ \alpha }} $ = 1.1847)

    Table  5.   Calibrated peak ground accelerations for the return period of 1000 a(${\boldsymbol{ \alpha }} $ = 1.1847)

    损伤状态 修正前 修正后
    PGA/g $ {u}_{\mathrm{D}\mathrm{I}} $ $ {u}'_{\mathrm{D}\mathrm{I}} $ PGA/g
    80%轻微损伤 0.210 1.531 $ \alpha \cdot DI=1.814 $ 0.253
    60%中度损伤 0.183 1.280 $ \alpha \cdot DI=1.516 $ 0.220
    40%重度损伤 0.229 1.714 $ \alpha \cdot DI=2.030 $ 0.276
    下载: 导出CSV

    表  6  重现期为475年时的地震动参数修正值($ {\boldsymbol{\alpha}} $ = 1.1332)

    Table  6.   Calibrated peak ground accelerations for the return period of 475 a(${\boldsymbol{ \alpha}} $ = 1.1332)

    损伤状态 修正前 修正后
    PGA/g $ {u}_{\mathrm{D}\mathrm{I}} $ $ {u}'_{\mathrm{D}\mathrm{I}} $ PGA/g
    80%轻微损伤 0.210 1.531 $ \alpha \cdot DI=1.735 $ 0.244
    60%中度损伤 0.183 1.280 $ \alpha \cdot DI=1.450 $ 0.213
    40%重度损伤 0.229 1.714 $ \alpha \cdot DI=1.942 $ 0.266
    下载: 导出CSV

    表  7  重现期为50年时的地震动参数修正值(${\boldsymbol{ \alpha}} $ =1.0123)

    Table  7.   Calibrated peak ground accelerations for the return period of 50 a(${\boldsymbol{ \alpha }} $ = 1.0123)

    损伤状态 修正前 修正后
    PGA/g $ {u}_{\mathrm{D}\mathrm{I}} $ $ {u}'_{\mathrm{D}\mathrm{I}} $ PGA/g
    80%轻微损伤 0.210 1.531 $ \alpha \cdot DI=1.550 $ 0.224
    60%中度损伤 0.183 1.280 $ \alpha \cdot DI=1.296 $ 0.195
    40%重度损伤 0.229 1.714 $ \alpha \cdot DI=1.735 $ 0.244
    下载: 导出CSV
  • 冯莉,樊燕燕,李子奇等,2020. 高速铁路连续梁桥全寿命周期概率地震损失分析. 铁道科学与工程学报,17(4):815−822. doi: 10.19713/j.cnki.43−1423/u.T20190690

    Feng L., Fan Y. Y., Li Z. Q., et al., 2020. Probabilistic whole life cycle seismic loss analysis of high-speed railway continuous girder bridge. Journal of Railway Science and Engineering, 17(4): 815−822. (in Chinese) doi: 10.19713/j.cnki.43−1423/u.T20190690
    李宏男,成虎,王东升,2018. 桥梁结构地震易损性研究进展述评. 工程力学,35(9):1−16. doi: 10.6052/j.issn.1000-4750.2017.04.0280

    Li H. N., Cheng H., Wang D. S., 2018. A review of advances in seismic fragility research on bridge structures. Engineering Mechanics, 35(9): 1−16. (in Chinese) doi: 10.6052/j.issn.1000-4750.2017.04.0280
    李杰,2018. 工程结构整体可靠性分析研究进展. 土木工程学报,51(8):1−10.

    Li J., 2018. Advances in global reliability analysis of engineering structures. China Civil Engineering Journal, 51(8): 1−10. (in Chinese)
    陆新征,叶列平,2010. 基于IDA分析的结构抗地震倒塌能力研究. 工程抗震与加固改造,32(1):13−18. doi: 10.3969/j.issn.1002-8412.2010.01.003

    Lu X. Z., Ye L. P., 2010. Study on the seismic collapse resistance of structural system. Earthquake Resistant Engineering and Retrofitting, 32(1): 13−18. (in Chinese) doi: 10.3969/j.issn.1002-8412.2010.01.003
    吕大刚,李晓鹏,王光远,2006. 基于可靠度和性能的结构整体地震易损性分析. 自然灾害学报,15(2):107−114. doi: 10.3969/j.issn.1004-4574.2006.02.018

    Lü D. G., Li X. P., Wang G. Y., 2006. Global seismic fragility analysis of structures based on reliability and performance. Journal of Natural Disasters, 15(2): 107−114. (in Chinese) doi: 10.3969/j.issn.1004-4574.2006.02.018
    于晓辉,周洲,吕大刚,2022. 钢筋混凝土框架结构的损伤状态相关余震易损性分析. 建筑结构学报,43(4):8−16. doi: 10.14006/j.jzjgxb.2020.0323

    Yu X. H., Zhou Z., Lü D. G., 2022. Damage-state-dependent aftershock fragility analysis for reinforced concrete frame structures. Journal of Building Structures, 43(4): 8−16. (in Chinese) doi: 10.14006/j.jzjgxb.2020.0323
    袁万城,王征南,庞于涛等,2016. 连续梁桥在主震-余震序列波下的地震易损性分析. 哈尔滨工程大学学报,37(12):1671−1676. doi: 10.11990/jheu.201509064

    Yuan W. C., Wang Z. N., Pang Y. T., et al., 2016. Seismic fragility analysis of a continuous girder bridge subject to an earthquake mainshock-aftershock sequence. Journal of Harbin Engineering University, 37(12): 1671−1676. (in Chinese) doi: 10.11990/jheu.201509064
    张明,金峰,2015. 结构可靠度计算. 北京:科学出版社.

    Zhang M. , Jin F. , 2015. Structural reliability computations. Beijing: Science Press. (in Chinese)
    Deierlein G. G. , Krawinkler H. , Cornell C. A. , 2003. A framework for performance-based earthquake engineering. In: Proceedings of 2003 Pacific Conference on Earthquake Engineering. Vancouver.
    Dong Y., Frangopol D. M., Saydam D., 2014. Sustainability of highway bridge networks under seismic hazard. Journal of Earthquake Engineering, 18(1): 41−66. doi: 10.1080/13632469.2013.841600
    Ghosh J., Padgett J. E., Sánchez-Silva M., 2015. Seismic damage accumulation in highway bridges in earthquake-prone regions. Earthquake Spectra, 31(1): 115−135. doi: 10.1193/120812EQS347M
    Hwang H. , Liu J. B. , Chiu Y. H. , 2001. Seismic fragility analysis of highway bridges. Memphis: Center for Earthquake Research and Information the University of Memphis.
    Karim K. R., Yamazaki F., 2001. Effect of earthquake ground motions on fragility curves of highway bridge piers based on numerical simulation. Earthquake Engineering & Structural Dynamics, 30(12): 1839−1856.
    Mackie K. , Stojadinović B. , 2003. Seismic demands for performance-based design of bridges. Berkeley: Pacific Earthquake Engineering Research Center, University of California, 1−113.
    Moehle J. , Deierlein G. G. , 2004. A framework methodology for performance-based earthquake engineering. In: Proceedings of the 13th World Conference on Earthquake Engineering. Vancouver: WCEE.
    Muntasir Billah A. H. M., Shahria Alam M., 2015. Seismic fragility assessment of highway bridges: a state-of-the-art review. Structure and Infrastructure Engineering, 11(6): 804−832. doi: 10.1080/15732479.2014.912243
    Pathak S., Watt S., Khennane A., et al., 2020. Dynamical systems approach for the evaluation of seismic structural collapse and its integration into PBEE framework. Soil Dynamics and Earthquake Engineering, 135: 106184. doi: 10.1016/j.soildyn.2020.106184
    Porter K. A. , 2003. An overview of PEER’s performance-based earthquake engineering methodology. In: Proceedings of the Ninth International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP9). San Francisco.
    Shome N., Allin Cornell C., Bazzurro P., et al., 1998. Earthquakes, records, and nonlinear responses. Earthquake Spectra, 14(3): 469−500. doi: 10.1193/1.1586011
    Yoon Y. H., Ataya S., Mahan M., et al, 2019. Probabilistic damage control application: implementation of performance-based earthquake engineering in seismic design of highway bridge columns. Journal of Bridge Engineering, 24(7): 04019068. doi: 10.1061/(ASCE)BE.1943-5592.0001443
  • 加载中
图(6) / 表(7)
计量
  • 文章访问数:  93
  • HTML全文浏览量:  40
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-11-27
  • 刊出日期:  2024-06-30

目录

/

返回文章
返回