标题

Size: px

Start display at page:

Download "标题"

砾刁
5 years ago
Views:

1 第 36 卷第 9 期哈尔滨工程大学学报 Vol 年 9 月 Journal of Harbin Engineering University Sep 摘非线性参数拟合的桥梁概率地震需求模型研究袁万城, 袁新哲, 庞于涛, 沈国煜 ( 同济大学土木工程防灾国家重点实验室, 上海 ) 要 : 针对云图法建立概率地震需求模型时部分工程需求参数 (EDP) 与地震动强度参数 ( IM) 间不满足对数线性回归关系假设的问题, 建立某三跨混凝土连续梁桥的 3D 有限元模型采用 IDA 方法建立桥梁概率地震需求模型, 并以 IDA 分析结果为基准, 在同一地震动参数下, 采用基于不同拟合方式的云图法建立概率地震需求模型并得到了结构在不同破坏状态下的易损性曲线研究表明, 当 EDP 与 IM 之间的数量关系不满足云图法对数线性关系基本假设时, 选择恰当的非线性参数拟合方式可以显著提高云图法建立概率地震需求模型的精度关键词 : 概率地震需求模型 ; 云图法 ;IDA 分析 ; 非线性拟合 ; 易损性分析 ; 连续桥梁 doi: / jheu 网络出版地址 :http: / / / kcms / detail / U html 中图分类号 :U442.5 文献标志码 :A 文章编号 : (2015) A probabilistic seismic demand model by nonlinear parameter fitting method YUAN Wancheng, YUAN Xinzhe, PANG Yutao, SHEN Guoyu ( State Key Laboratory for Disaster Reduction in Civil Engineering, Tongji University, Shanghai , China) Abstract:In order to solve the problem that when establishing a probabilistic seismic demand model by the cloud method, some Engineering Demand Parameters ( EDPs) and ground motion intensity ( IM) parameters do not meet the assumption of logarithmic linear regression relationship, a 3D finite element model of a three span concrete con tinuous beam bridge was established. The bridge probabilistic seismic demand model was first established by incre mental dynamic analysis( IDA). Taking the IDA analysis result as the benchmark, under the same seismic motion parameters, and using the cloud method based on various fitting modes to establish a probabilistic seismic demand model, the fragility curve under different breaking states was obtained. The results show that, when the quantitative relationship between EDP and IM does not meet the basic logarithmic linear relationship assumption of the cloud method, choosing appropriate nonlinear parameter fitting mode can significantly improve the precision of the proba bilistic seismic demand model established by the cloud method. Keywords:probabilistic seismic demand model; cloud method; IDA; nonlinear fitting; fragility analysis;continu ous beam bridge 美国太平洋地震工程研究中心 (PEER) 在 2005 年提出新一代基于性能的地震工程研究框架其中, 概率性的地震需求分析作为该框架的重要一环, 其核心在于建立概率地震需求模型概率地震需求模型通常被假设为对数正态分布模型 [1], 可以用 2 种方法来建立 : 云图法和 IDA 分析方法云图法假设工程需求参数 EDP 的均值与地震动参数 IM 之间满足对数线性关系, 而模型对数标准差在整个 IM 范围内保持不变 [2] IDA 分析方法通过将地震波调幅至不同强度, 统计 EDP 在不同 IM 时的分布特性来建立概率地震需求模型近年来, 国内外学者对云图法和 IDA 方法进行了广泛的研究与运用 :Kevin 收稿日期 : 网络出版时间 : 基金项目 : 土木工程防灾国家重点实验室基金资助项目 ( SLDRCE14 B 14); 国家自然科学基金资助项目 ( , ). 作者简介 : 袁万城 (1964 ), 男, 研究员, 博士生导师 ; 袁新哲 (1989 ), 男, 硕士研究生. 通信作者 : 袁新哲, E mail:yxz_hz@ 163.com. [3] Mackie 等利用云图法对结构进行非线性动力时程分析, 来获得结构的概率地震需求模型 ; 冯清 [4] 海运用 IDA 分析方法对特大桥梁的地震易损性与地震风险概率进行了分析 ;Payam Tehrani [5] 等运用快速 IDA 算法对不同地震类型作用下的连续四跨直梁桥地震反应进行了分析研究一般来说, 在使用云图法建立概率地震需求模型时, 为了保证其分析精度, 需要对所选取的 IM 参数的有效性进行评估, 确保其与 EDP 均值服从云图法对数线性关系 [6] 的假设 Luco Nicolas 等研究了速度脉冲型地震 IM 参数与剪切型多自由度体系顶层非线性变形需求的相关性, 分析了不同 IM 参数的有效性但这种 IM 参数有效性评估的方法只适用于特定的地震类型和结构类型, 由于地震动的强烈不确定性和结构类型的多样性,IM 参数的有效性也是不确定的因此, 通过改变 EDP 与 IM 参数之间的回归拟合方式解决 EDP 与 IM 参数之间的非线性对数关系问题, 可以提高云图法建立概率地震需求模型的精度

2 第 9 期袁万城, 等 : 非线性参数拟合的桥梁概率地震需求模型研究概率地震需求模型与易损性分析 1.1 概率地震需求模型桥梁结构地震风险分析中最主要的步骤是建立桥梁结构的地震易损性曲线, 而桥梁结构的地震易损性曲线建立的关键步骤是进行概率地震需求分析, 得到概率地震需求模型分别可以采用云图法和 IDA 分析方法来得到桥梁结构的概率地震需求模型对于云图法而言, 重点在于运用能量法则进行统计回归分析得到桥梁结构构件的地震峰值响应与地 [7] 震动强度参数的关系其基本步骤为 : 选取可以反映结构所在场地区域地质情况的地震动纪录 ; 建立结构三维有限元模型 ; 进行非线性时程分析, 得到 IM 与 EDP 的样本 ; 通过能量法则, 对样本进行回归分析得到结构构件的概率地震需求模型根据能量法则,IM 与 EDP 的均值 μ D 一般服从对数线性关系 : μ D = aim b or ln( μ D ) = lna + bln( IM) (1) 式中 :a b 为回归参数, 两边取对数后方便运用最小二乘法进行拟合, 其拟合方差可由下式计算 : n i = 1 [ln(edp i ) - (lna + blnim i )] 2 σ lnd = n - 2 (2) 式中 : n 为非线性时程次数,EDP i 为结构的地震峰值响应得到模型均值和方差后, 即可得到相应的概率地震需求模型 : é P[EDP > d IM = x] = 1 - Φ lnd / μ D ù ë ê σ lnd û ú (3) IDA 分析方法与云图法的主要区别在于 IDA 通过统计结构需求在不同地震动强度下的分布特性, 获得模型均值和方差, 然后即可由式 (3) 建立结构的概率地震需求模型云图法基于对数线性关系假设建立的概率地震需求模型具有高度相似性, 而 IDA 分析方法统计在不同地震动强度下结构响应分布特性, 可以视为一种更加精确的概率地震需求模型建立方法因此, 以 IDA 分析结果为基准, 可以比较不同拟合方式下的云图法建立概率地震需求模型的精度, 研究非线性拟合方式下的概率地震需求模型 1.2 地震易损性分析概率地震需求模型中, 结构的地震需求与地震动输入强度息息相关, 而结构的能力通常用来定义不同的损伤状态函数 [8] 对于一个特定的损伤状态, 结构构件的失效概率 P f 可由桥梁结构的地震需求 S d 与结构的能力 S c 表示 : P f = P[ S d S c 1] (4) 由于地震需求 S d 与结构的能力 S c 均假设为对数正态分布, 因而失效概率 P f 可用标准正态分布函数表示为 é P f = Φ - ln( S c ê ù ) ú ê S d ú (5) ê ë β 2 + ú d β2 c û 式中 : β c 与 β d 分别为能力与需求的对数正态分布的方差由于 IDA 分析方法并不需要假定特定级别的地震动强度下的地震响应的概率分布函数的类型, 因而 IDA 方法在计算结构损伤概率的时候相对简单, 可由下式计算 : P[DI LS i IM] = n i (i = 1,2,3,4) (6) N 式中 : n i 为超过极限状态函数的 IDA 曲线的数目,N 为 IDA 曲线的总数,DI 为损伤指标,LS i 为给定损伤指标的极限状态值基于 IDA 易损性分析结果, 得到采用不同拟合方式的云图法下桥梁结构响应的概率地震需求模型, 并建立桥梁构件的易损性曲线, 研究不同拟合方式对易损性分析精度的影响 2 计算模型及地震波选取 2.1 模型简介以一座三跨混凝土连续梁桥为例, 墩截面采用正方形柱, 截面为 1.5 m 1.5 m, 墩高为 8 m, 跨径布置为 3 m 38 m 图 1 给出了混凝土桥梁的总体布置图主梁采用 C50 预应力混凝土, 桥墩混凝土为 C40, 钢筋采用 HRB335, 纵向配筋率为 0.8%, 配箍率 0.5% 使用 OpenSees 建立桥梁 3D 非线性有限元模型, 上部结构在地震作用下基本保持弹性, 使用线弹性梁柱单元模拟, 墩柱采用弹塑性纤维单元模拟, 用碰撞单元考虑地震作用下主梁与桥台的碰撞效应为简化起见, 本例中不考虑桩土共同作用, 假定墩底固结分别选取桥梁墩台位移 a D 固定支座位移 a Z 为需求参数, 建立其在同一地震动参数下的概率地震需求模型 Fig. 1 图 1 混凝土连续梁桥总体布置图 Overall layout of RC continuous beam bridge 2.2 地震波选取从 PEER 强震记录数据库中选取了 40 条地震波, 这些波的震级为 6.0 ~ 6.7, 并且各条地震波的震中距均大于 10 km, 地震持时均取为 30 s 图 2 给出了 40 条波阻尼为 5% 的伪加速度反应谱曲线以及反应谱概率分布考虑云图法模型精度受样本点数目

3 1214 哈尔滨工程大学学报第 36 卷影响, 对选取的 40 条地震波进行调幅处理, 得到 80 条地震波输入来增加样本点数目, 避免因样本点数目较少而带来的分析误差在进行 IDA 分析时, 对选取的 40 条地震波进行 0.1g-1.0 g 的调幅处理, 来统计结构响应在不同地震动强度下的分布特性 (a) 墩台位移 a D (a) 1# (b) 2# 图 2 地震波反应谱及概率分布图 Fig. 2 Response spectrum and probability distribution of 40 seismic waves 地震动可以通过 PGA PGV SA(T) 等强度指标 [9] 来描述 Padgett 等通过对地震动参数进行效率可适用性以及危害的可计算性评价表明 :PGA 是比较适合的地震动强度指标因此, 本文选用 PGA 作为地震动强度参数由于选取的研究模型为三跨混凝土连续梁桥, 桥型较简单, 地震波输入时只考虑了沿纵桥向的单方向输入 3 计算分析 3.1 云图法与 IDA 分析结果分别采用基于对数线性假设的云图法和 IDA 方法得到桥梁墩台位移 a D 响应均值与 PGA 的对数回归曲线, 如图 3 所示从图中对比结果可以发现, 墩台位移 a D 和 PGA 间具有显著的对数线性关系, 满足云图法的基本假设但支座位移 a Z 与 PGA 之间具有显著的非线性对数关系, 其概率地震需求模型不再适用基于线性对数关系假设的云图法来建立通常 EDP 均值与 IM 间的对数线性关系假设不满足时, 应该选择其他 IM 参数然而如果找到一种能够拟合 EDP 均值与 IM 参数之间非线性对数关系的方式, 则同样可以提高云图法概率地震需求模型的精度 (b) 墩台位移 a Z (c) 墩台截面曲率 a Q 图 3 云图法与 IDA 分析结果比较 Fig. 3 Comparison between cloud method and IDA method 为了拟合结构响应 a Z,a Q 与 IM 间的非线性对数关系, 基于最小二乘原理采用一次线性二次抛物线和三次抛物线拟合 EDP 均值与 PGA 参数间的对数回归关系, 得到的桥梁墩台位移 a D 与 PGA 的对数曲线关系如图 4 所示由图可见, 虽然墩台位移 a D 与 PGA 间具有显著的线性对数回归关系, 但依然给出其不同拟合方式下的对数回归关系曲线, 来对比说明回归拟合方式的改变对概率地震需求模型精度的影响改变拟合方式后墩台位移 a D 基于非线性拟合的回归分析精度并未有提高, 非线性拟合曲线与 IDA 曲线偏离而支座位移 a Z, 桥墩截面曲率 a Q 采用非线

4 第 9 期袁万城, 等 : 非线性参数拟合的桥梁概率地震需求模型研究 1215 性拟合方式可以得到与 PGA 间的非线性对数关系从图 4(b),( c) 可以看出, 三次抛物线回归拟合方式得到的 IM 对数关系曲线与 IDA 结果最为吻合这表明当所选地震动参数不能满足云图法对数线性关系假设时, 采用三次抛物线可以建立 EDP 均值与 PGA 参数间的非线性对数关系, 根据不同地震动强度下桥梁结构响应, 即可得到基于非线性拟合方式的云图法概率地震需求模型以采用非线性拟合在不改变地震动参数 IM 的条件下建立 EDP 均值与 IM 的关系曲线, 可以提高云图法概率地震需求模型的精度同时需指出, 采用最小二乘法的高次抛物线回归拟合并非唯一的非线性回归拟合方式, 也可采用其他非线性回归拟合手段来建立更加精确的 EDP 与 IM 参数之间的数量关系 3.2 不同拟合方式的易损性分析采用非线性拟合方式建立结构概率地震需求模型后, 即可对构件易损性进行分析为比较不同拟合方式对易损性分析精度的影响, 得到了不同拟合方式下桥墩的易损性曲线选择桥墩而不是相对来说更加容易破坏的支座进行易损性分析是因为桥墩截面曲率 a Q 与 PGA 的非线性更加明显进行易损性分析时, 选用变形作为破坏准则 [10], 结构性能划分标准分为 :1) 轻微破坏 ;2) 中等破坏 ;3) 严重破坏 ;4) 倒塌 (a) 墩台位移 a D (a) 轻微损伤 (b) 支座位移 a Z (b) 中等损伤 (c) 墩台截面曲率 a Q 图 4 不同拟合方法的云图法与 IDA 分析结果比较 Fig. 4 Comparison between cloud methods under different fitting methods and IDA 因此, 当云图法对数线性关系假设不成立时, 可 (c) 严重损伤

5 1216 哈尔滨工程大学学报第 36 卷图 5 Fig. 5 (d) 完全破坏不同拟合方式下桥墩在不同损伤状态的易损性曲线 Fragility curves under different damage states by different fitting ways 图 5 分别对应给出了 4 种破坏状下基于 IDA 方法和不同拟合方式云图法得到的易损性曲线从图 5 可以看出, 采用三次抛物线拟合得到的易损性曲线与 IDA 分析结果最为逼近因此, 采用非线性拟合方法建立 EDP 与 IM 参数之间的关系, 可以对结构易损性进行更加精确的计算, 提高桥梁结构可靠度评估的准确性 4 结束语本文以某三跨 RC 连续梁桥为例, 选择 PGA 为地震动指标, 桥梁墩台位移 a D 为工程需求参数, 以 IDA 分析结果为基准, 采用非线性拟合方式的云图法建立了桥梁结构的概率地震需求模型, 并得到了桥墩构件的易损性曲线, 得出以下结论 : 对结构进行概率地震需求分析时, 如果所选 EDP 与 IM 参数之间不满足云图法对数线性回归关系的基本假设, 可以通过非线性拟合方式建立 EDP 与 IM 参数之间的数量关系, 建立结构的概率地震需求模型基于非线性拟合的概率地震需求分析与易损性分析结果与 IDA 分析结果吻合良好, 计算精度较高参考文献 : [1] 张菊辉, 胡世德. 桥梁地震易损性分析的研究现状 [ J]. 结构工程师, 2005, 21(5): ZHANG Juhui, HU Shide. State of the art of bridge seismic vulnerability analysis research [ J ]. Structural Engineers, 2005, 21(5): [2] 于晓辉, 吕大刚, 王光远. 关于概率地震需求模型的讨论 [J]. 工程力学, 2013, 30(8): YU Xiaohui, LYU Dagang, WANG Guangyuan. Discussions on probabilistic seismic demand models [ J ]. Engineering Mechanics, 2013, 30(8): [3] MACKIE K, STOJADONOVIC B. Probabilistic seismic de mand model for California highway bridges [ J]. Journal of Bridge Engineering, 2001, 6(6): [4] 冯清海. 特大桥梁地震易损性与风险概率分析 [ D]. 上海 : 同济大学, FENG Qinghai. Study on seismic vulnerability and risk prob ability analysis of super Large bridge[ D]. Shanghai: Tongji University, [5] TEHTANI P, MITCHELL D. Seismic response of bridges subjected to different earthquake types using IDA[ J]. Jour nal of Earthquake Engineering, 2013, 17(3): [ 6] LUCO N, CORNELL C A. Structure specific scalar intensity measures for near source and ordinary earthquake ground motions[ J]. Earthquake Spectra, 2007, 23(2): [7]PANG Yutao, WU Xun, SHEN Guoyu, et al. Seismic fra gility analysis of cable stayed bridges considering different sources of uncertainties[ J]. Journal of Bridge Engineering, 2014, 19(4): [8] 钟剑, 庞于涛, 沈国煜, 等. 采用拉索减震支座的斜拉桥地震易损性分析 [ J]. 同济大学学报 : 自然科学版, 2014, 42(3): ZHONG Jian, PANG Yutao, SHEN Guoyu, et al. Seismic fragility analysis of cable stayed bridges retrofitted with ca ble sliding friction aseismic bearing [ J]. Journal of Tongji University:Natural Science, 2014, 42(3): [9]PADGETT J E, NIELSON B G, DESROCHES R. Selection of optimal intensity measures in probabilistic seismic demand models of highway bridge portfolios [ J]. Earthquake Engi neering & Structure Dynamics, 2008, 37(5): [10] 沈国煜, 袁万城, 庞于涛. 基于 Nataf 变换的桥梁结构地震易损性分析 [J]. 工程力学, 2014, 31(6): SHEN Guoyu, YUAN Wancheng, PANG Yutao. Bridge seismic fragility analysis based on NATAF transformation [ J]. Engineering Mechanics, 2014, 31(6): [11] 庞于涛, 袁万城, 沈国煜, 等. 基于结构易损性分析的公路桥梁可靠度研究 [J]. 哈尔滨工程大学学报, 2012, 33(9): PANG Yutao, YUAN Wancheng, SHEN Guoyu, et al. Re liability analysis of highway bridges based on structural vul nerability analysis[ J]. Journal of Harbin Engineering Uni versity, 2012, 33(9):

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva 39 6 2011 12 Journal of Fuzhou University Natural Science Edition Vol 39 No 6 Dec 2011 DOI CNKI 35-1117 /N 20111220 0901 002 1000-2243 2011 06-0923 - 07 350108 105 m 14 69% TU311 3 A Seismic analysis of