考虑支座滑移及挡块破坏的中小跨径梁桥地震易损性研究1

2019-07-28宋帅王帅吴刚

震灾防御技术 2019年4期

宋帅王帅吴刚

1）太原理工大学，土木工程学院，太原 030024

2）华东交通大学，土木建筑学院，南昌 330013

引言

采用板式橡胶支座的中小跨径梁桥因施工简便、质量可靠，在我国公路交通网中占有较大比例。通过总结汶川地震中梁桥的破坏，发现板式橡胶支座普遍发生了支座滑移，主梁位移偏大，主梁与挡块或桥台间的碰撞现象明显，而下部结构的破坏却相对较轻。这种典型震害表现与既有延性抗震设计方法的预期目标具有一定差异。传统的延性抗震设计假定支座在地震作用下表现为线弹性，桥梁通过桥墩潜在塑性铰的转动实现耗能。针对此问题，王克海等（2014）通过总结既有震害，提出考虑支座摩擦滑移的中小跨径梁桥抗震设计方法；王凯睿等（2017）通过振动台试验研究不同设计条件下板式橡胶支座的滑移效应，提出相应的抗震措施。针对采用板式橡胶支座的中小跨径梁桥，李建中等（2015）提出一种新型隔震系统；赵宇翔（2018）对不同支座形式的简支梁桥地震响应开展专门研究。由此可见，在中小跨径梁桥抗震设计中，考虑板式支座滑移效应的设计理念已逐渐被广大学者所认识（吴刚等，2018）。然而，在目前中小跨径梁桥的地震易损性的研究中，尚未充分考虑支座滑移效应及挡块破坏的影响（张磊鑫等，2018）。

作为基于性能地震评估的重要组成部分，地震易损性分析在桥梁抗震领域中得到广泛应用，如桥梁抗震设计优化、抗震加固方法优选及地震风险评估等（李宏男等，2018；Billah等，2015）。本文采用基于离散概率的地震易损性分析方法，综合考虑地震动的随机性及梁桥结构中参数的不确定性，以一座典型的3 跨预应力混凝土连续梁桥为例，对比分析包含板式橡胶支座及挡块的不同分析模型计算结果，建立桥梁构件及系统的地震易损性曲线，研究支座滑移及挡块破坏对中小跨径梁桥地震易损性的影响。

1 桥梁地震易损性分析方法

1.1 基本理论

地震易损性是指在给定地震作用下，结构地震需求超过指定极限状态的概率水平。用结构地震需求及抗震能力可表示为式（1）所示的条件概率形式：

1.2 分析流程

结合结构概率地震需求分析及极限状态，基于回归分析即可建立结构地震易损性曲线，分析流程如图1 所示。

2 工程背景及分析模型

2.1 工程背景

以一座连续梁桥为例（3×30m），主梁为预应力混凝土单箱单室箱梁，材料为C50 混凝土。桥墩为圆形单墩形式，直径1.8m，桥墩高9m，材料为C30 混凝土，纵向HRB335 钢筋配筋率1.1%。支座为板式橡胶形式，橡胶层厚76mm。墩顶横桥向混凝土挡块高0.5m，厚0.3m，挡块内设置倒U 形直径16mm 的HRB335 钢筋19 根，水平箍筋间距0.1m。桥址为II类场地，设防烈度为7 度。

图1 桥梁地震易损性分析流程 Fig. 1 Flow chart of bridge seismic vulnerability analysis

2.2 分析模型

为分析支座滑移及挡块破坏的影响，采用OpenSees 建立3 种全桥对比模型为：

（1）模型一：基于传统延性设计，支座采用线弹性模型模拟，挡块采用双折线模型模拟；

（2）模型二：考虑支座摩擦滑移，支座采用基于库仑摩擦的双折线模型模拟，挡块采用双折线模型模拟；

（3）模型三：支座采用基于库仑摩擦的双折线模型模拟，采用考虑损伤的模型模拟挡块（汤虎等，2016）。

桥梁有限元模型如图2 所示，主梁采用弹性梁柱单元模拟。桥墩采用纤维单元模拟，将截面划分为钢筋纤维、约束混凝土和非约束混凝土纤维，分别采用Concrete02 材料模型及Steel02 材料模型模拟混凝土及钢筋，无约束混凝土的峰值应力及峰值应变均值分别取28.03MPa 及0.002，约束引起的混凝土强度增加系数取1.06；钢筋屈服应力及弹性模量均值分别取381.65MPa 及2.06×105MPa。对于桥台，参考加州桥梁抗震设计规范（Caltrans，2013），采用理想弹塑性单元进行简化（Aviran 等，2008），将弹塑性材料属性赋予零长度单元，材料的屈服力及刚度分别为1.0×105kN 及1.0×107kN/m。支座弹性刚度Kz取为6280kN/m，在库仑摩擦模型中，其屈服力取为420kN。采用程序中内置的Gap 单元模拟碰撞效应，主要包括主梁与桥台的碰撞及主梁与横向挡块间的碰撞。在挡块弹塑性模型中，其屈服力取支座竖向力的15%，初始刚度为1.0×105kN/m；在考虑损伤的挡块模型中，A、B、C 3 个关键点对应的力分别为593.5kN、878.6kN 及555.8kN，相应的位移分别为3.0mm、24.8mm 及76.4mm。

图2 桥梁有限元模型 Fig.2 Finite element model of the bridge

2.3 结构参数不确定性

桥梁结构参数的不确定性复杂多样。根据已有研究（宋帅，2017；Nielson，2005），汇总得到梁桥中主要的不确定参数及其概率分布，如表1 所示。

表1 梁桥不确定参数及其概率分布 Table 1 Uncertainty parameters and distribution types of girder bridges

2.4 地震动选择

地震动自身具有强烈的不确定性。为接近实际，根据地震场地条件，基于太平洋地震工程中心新一代强震记录数据库，选取100 条实际地震波，选取准则如下（Bommer 等，2004）：①地震波均来源于自由场地；②震中距小于100km；③震级大于5.5 级；④平均剪切波速为260—510m/s。

基于选取的地震动，绘制其加速度反应谱，如图3 所示。

采用OpenSees 有限元平台对3 种全桥模型进行非线性动力时程分析，将选取的地震动沿纵、横桥向同时输入。

图3 地震动加速度反应谱 Fig. 3 Acceleration spectrum of ground motions

3 构件破坏极限状态

3.1 桥墩破坏状态

在基于性能的地震评估中，常用曲率延性定义桥墩的破坏极限状态（李立峰等，2012），本文以墩底截面曲率延性比μφ作为桥墩极限状态指标，定义轻微、中等、严重及完全破坏4 种状态。结合单墩截面分析，并与位移延性对应，得到不同破坏状态下的桥墩曲率延性比，如表2 所示。

表2 构件极限破坏状态及其限值 Table 2 Damage states and limit values of components

3.2 支座破坏状态

已有研究中，对于板式橡胶支座，以剪切应变达100%、150%、200%和250%所对应的剪切变形作为支座轻微、中等、严重及完全破坏的极限指标（李立峰等，2012），相应的支座位移分别为76mm、114mm、152mm 和190mm，如表2 所示。

4 桥梁易损性分析

4.1 概率地震需求

基于不确定参数的概率分布特征，采用拉丁超立方分层抽样方法建立3 种全桥模型的概率需求分析样本，结合选取的100 条地震波，进行非线性动力时程分析，限于篇幅，仅列出模型一中支座及桥墩概率地震需求，如图4 所示。

图4 构件概率地震需求 Fig. 4 Seismic demands of bridge components

由图4 可知，在对数空间下采用线性函数对构件地震需求与地震动强度进行描述时，其95%的置信区间范围均较小，因此采用式（3）所示的指数函数对梁桥地震需求进行描述，满足统计分析中置信水平的要求。

4.2 支座易损性曲线

结合支座概率地震需求及其极限状态，根据式（2）计算得到3 种模型中板式橡胶支座易损性曲线，如图5 所示。

图5 支座易损性曲线对比 Fig. 5 Comparison of seismic vulnerability curves of bearings

由图5 可知，在不同破坏状态下，模型三中支座破坏概率均最大，模型一中支座破坏概率最小。可见由于模型三同时考虑了支座滑移效应及挡块损伤退化，混凝土挡块在变形达到一定值时，其承载力下降，限位能力也相应退化，甚至完全失去限位能力，导致主梁位移显著增加，支座滑动位移也明显增加。而模型一仅考虑挡块发生屈服，未模拟其限位能力的退化，且支座模型为线弹性，对主梁位移具有一定限制作用，导致橡胶支座的位移也较小。

4.3 桥墩易损性曲线

结合桥墩概率地震需求及其极限状态，根据式（2）计算得到3 种模型中桥墩易损性曲线，如图6 所示。

图6 桥墩易损性曲线对比 Fig. 6 Comparison of seismic vulnerability curves of piers

由图6 可知，模型一中桥墩破坏概率最大，而模型三中桥墩破坏概率最小，桥墩与支座破坏规律正好相反，主要原因为：支座及挡块的力学特性直接影响了传至桥墩的地震作用大小，模型三中支座发生了滑移，起到较好的隔震作用，同时挡块在一定变形后发生退化，导致总体传至桥墩的地震作用较小，而模型一中未考虑板式橡胶支座的滑移效应，上部结构的地震作用通过支座及挡块传至桥墩。此外，4 种破坏状态下，桥墩破坏概率均小于支座破坏概率，说明支座是梁桥中较为薄弱的一环。模型三中桥墩严重破坏及完全破坏的概率均较小，即桥墩发生混凝土压碎及倒塌的概率较小，这与汶川地震中观察到的中小跨径梁桥桥墩破坏表现一致。

4.4 系统易损性曲线

为对比支座滑移效应及挡块退化对桥梁整体系统地震易损性水平的影响，基于一阶界限法的上界建立系统易损性曲线，如图7 所示。

图7 桥梁系统易损性曲线对比 Fig. 7 Comparison of seismic vulnerability curves of bridge system

由图7 可知，在4 种破坏状态下，模型三的系统易损性总体上最大。模型一与模型二在轻微及中等破坏状态下，系统破坏概率接近，在严重及完全破坏状态下，模型二的系统破坏概率大于模型一。结合支座及桥墩易损性，对于中小跨径梁桥，其整体破坏状态主要取决于支座的破坏状态。由此可知，支座滑移效应及挡块对中小跨径梁桥桥梁构件及系统地震易损性的影响较显著，合理利用及准确模拟支座的滑移效应及挡块的破坏对桥梁抗震设计及地震评估具有重要意义。