大应变FRP加固连续梁桥地震易损性分析_白玉磊.pdf

第 卷 第 期 年 月哈 尔 滨 工 业 大 学 学 报 .：大应变 加固连续梁桥地震易损性分析白玉磊，张玉峰，梅世杰，贾俊峰，杜修力（城市与工程安全减灾教育部重点实验室（北京工业大学），北京）摘 要：为评估地震作用下大应变纤维增强聚合物（，）加固桥梁的抗震性能，以一典型三跨连续梁桥为例进行了研究。目前关于 约束混凝土往复轴压模型的研究较少，限制了 加固桥梁的抗震分析。因此介绍了作者课题组提出的简化 约束混凝土往复轴压模型，根据此模型并结合地震易损性分析理论，基于 平台，采用增量动力分析（，）方法，对三跨连续梁桥在 加固前后的地震易损性进行了分析。结果表明：轻微损伤和中等损伤状态下桥梁墩柱基本处于未受损坏的线弹性阶段，加固对其地震易损性影响不大；严重破坏和完全倒塌状态下，未加固桥梁的损坏超越概率明显大于 加固的桥梁。提出的简化 约束混凝土往复轴压模型可为 加固桥梁的抗震分析提供参考。关键词：桥梁工程；地震易损性；增量动力分析；桥梁墩柱；中图分类号：文献标志码：文章编号：（），（），）：（），（），：；收稿日期：；录用日期：；网络首发日期：网络首发地址：基金项目：国家自然科学基金（，）作者简介：白玉磊（），男，教授，博士生导师通信作者：贾俊峰，为提高桥梁的抗震能力，在过去 年里，纤维增强复合材料（，）在桥梁加固领域得到了广泛应用。常见的 材料如碳纤维（，），玻璃纤维（，）等具有线性的拉伸应力 应变关系和较小的断裂应变（小于），一般被称为传统。近年来，一种新发展的断裂应变大于的大应变 （，）得到了越来越多的学者关注。是由废弃塑料中的聚萘二甲酸乙二醇酯（，）和 聚 对 苯 二 甲 酸 乙 二 醇 酯（，）制成的，是一种绿色可再生环保材料。与传统 相比，表现出双线性的拉伸应力 应变关系和较大的断裂应变特点。文献 表明 约束混凝土表现出优越的延性行为和超高的耗能能力，因此 更适合作为以增强延性为主要目标的抗震加固材料，有望为桥梁加固提供一种更经济、延性更好的解决方案。准确地预测 约束混凝土在往复轴压荷载下的应力 应变行为是 加固桥墩抗震分析的基础。可是，目前大多数往复轴压模型 是基于传统 约束混凝土柱的试验结果提出的，关于 约束混凝土往复轴压模型的研究较少。由于约束混凝土在往复轴压荷载下的应力 应变响应取决于其外部约束材料的特性，而 的材料特性与传统 明显不同，所以文献，指出如果将传统 约束混凝土模型直接应用于 加固桥墩的抗震分析会低估其变形能力。最近作者课题组提出了考虑塑性损伤的简化 约束混凝土往复轴压应力 应变模型。该简化模型采用类似 等的线性加卸载准则来模拟 约束混凝土在往复轴压下的加卸载曲线并可以对 加固墩柱拟静力试验结果做出准确地预测。鉴于此，本文在作者课题组提出的简化 约束混凝土往复轴压模型基础上，对一典型三跨连续梁桥在 加固前后的地震易损性进行了分析，研究了 加固桥梁地震损伤变化规律，以期为同类工程的抗震加固提供一定的技术支持。简化 约束混凝土往复轴压模型 单调荷载下的应力 应变模型作者课题组发展的简化 约束混凝土往复轴压模型的包络线采用了文献提出的基于刚度的单向轴压模型来表示。该模型是针对 约束混凝土提出的，并且同样适用于传统 约束混凝土。此模型由三部分组成：第一部分抛物线阶段、第二部分线性阶段以及第三部分线性阶段，见图。第一部分和第二部分相交于点（，），第二部分与第三部分相交于点（，），模型详细信息参考文献。包络线模型的数学表达式：当第一段直线的斜率即刚度 时（强约束）（），（），|（）当第一段直线的斜率即刚度 时（弱约束），（），（），|（）式中：与 分别为 约束混凝土的压应力和压应变；为约束混凝土的极限压应变，为其对应的极限压应力；为未约束混凝土的抗压强度，为未约束混凝土的抗压弹性模量，和分别为包络线模型的线性第二段和第三段的斜率。部分卸载/再加载完全卸载/再加载轴向应变轴向应力包络线:文献4加卸载准则:文献16c(t，t)Ec(re，re)(cu，cu)(un，un)(cu，cu)(p1，pseudo，0)E2E*2(*t，*t)cfco图 简化 约束混凝土往复轴压模型 加卸载准则和塑性应变卸载起点处的应变和应力分别称为卸载应变和卸载应力，再加载起点处的应变和应力分别称为再加载应变 和再加载应力，见图。文献 提出的简化模型采用线性加卸载准则来模拟 约束混凝土往复轴压曲线的加卸载路径，为使预测的线性加卸载路径能准确代表整个试验加卸载曲线的位置，提出了伪塑性应变，概念。该模型对应的线性加卸载准则为，（，）（）伪塑性应变为，（），（），（）桥梁易损性分析方法桥梁的地震易损性是指在任意给定的地震强度 下结构的需求 达到或超出其抗震能力 的条件概率，可表示为|（）第 期白玉磊，等：大应变 加固连续梁桥地震易损性分析式中（）为标准正态分布函数，为桥梁地震需求对数标准差，为桥梁抗震能力对数标准差。根据文献 建议，当易损性曲线以地震动峰值加速度为自变量时，取为 。参考文献，假设结构地震需求服从对数正态分布时，结构地震需求与地震动参数之间的关系为 （）将上式两边同时取对数可得（）（）式中 与 为回归系数，可通过 计算提取桥梁结构响应数据进行统计回归得到。将式（）代入式（）可得 ()（）地震易损性分析的核心是得出易损性曲线。其主要分析过程可分为四部分：建立桥梁有限元模型；地震动记录选取和调幅；桥梁损伤状态等级划分；通过 法对桥梁模型进行地震反应分析得出易损性曲线。梁桥有限元分析模型 工程概况选取一座典型公路连续梁桥作为研究对象，桥梁总长为 ，跨径布置为 ，主梁为高 ，宽 的单箱双室箱梁，桥墩采用无盖梁的单排圆形双柱墩，直径为 ，墩身高度为，两墩的中心距为 ，桥址为 类场地，桥梁示意见图。主梁和桥墩均采用 混凝土，一期和二期恒载合计为 ，纵向钢筋和箍筋都采用 钢筋，墩柱截面配有 根直径 的纵筋，箍筋直径为 ，纵筋配筋率为 ，配箍率为 。支座采用盆式橡胶支座，其平面布置方式和可移动方向见图。(a)桥梁纵断面9 000800(b)桥梁横断面(c)主梁横断面尺寸1 3251808001405501 3251 2255050242.5242.565657108030 803080308030603060306030603018020505025255025180505025图 桥梁示意（）（）固定铰支座单向盆式橡胶支座单向盆式橡胶支座双向盆式橡胶支座30 m30 m30 m51627384图 支座平面布置 桥梁数值模拟主梁在地震作用下通常不会发生损伤，因此采用弹性梁单元模拟，根据图（）计算其截面几何特性见表。墩柱采用 中的 哈 尔 滨 工 业 大 学 学 报 第 卷 来模拟，每个墩柱被划分为 个单元，每个单元被分为 个积分点。钢筋采用 模型，考虑箍筋对混凝土的约束作用，将截面划分为核心区混凝土、保护层混凝土及钢筋纤维截面，见图。此外，由于算例墩柱为悬臂墩，因此仅对其塑性铰区进行外包 加固，材料特性见表。通过二次开发的方法将上述简化 约束混凝土模型添加到 的单轴材料库中来模拟 加固区混凝土的特性。对于未加固区混凝土，其保护层采用 模型，核心区采用 模型。表 主梁截面几何特性 单元类型截面面积弹性模量剪切模量截面扭转惯性矩 绕 轴截面抗弯惯性矩 绕 轴截面抗弯惯性矩 表 材料参数 类型层数 名义厚度 第一段弹性模量 第二段弹性模量 断裂应变，直径 D=1.4 m纵筋保护层核心区图 墩柱纤维截面 文献指出主梁与桥台之间的碰撞对桥梁的动力响应有重要影响，在桥梁抗震分析中有必要准确模拟桥台与主梁的碰撞效应。本文采用接触单元方法模拟结构构件间碰撞效应，通过在梁体与桥台之间设置 个零长度单元并对零长度单元赋予 内置的 模型来实现。该材料是 等提出的双折线非线性碰撞模型，见图（）。其中，为初始刚度，为屈服后刚度，为有效刚度，为伸缩缝宽度。算例 取，其他参数计算为 ，。活动盆式支座的恢复力模型见图（），为临界滑动摩擦力，为临界滑动位移。参考 公路桥梁抗震设计规范，算例支座的初始刚度 。背墙台后填土伸缩缝主梁被动主动桥梁支座前墙承台桩基(b)支座模型(a)Impact 模型XXFFK2K1KKeffFmFy-FmaxFmaxxyxyxm-xyGapGap图 桥台和支座模型 根据桥梁结构的整体构造布置，采用 建立三维非线性有限元分析模型。由于场地地质条件较好，不考虑基础弹性变形的影响，仅在墩底施加固定约束。主梁顺桥向为 轴，横桥向为 轴，竖向为 轴。为方便叙述，图 给出了桥梁各个主要组成构件的编号。由图 支座布置情况可知，支座为固定铰支座，其和 支座共同约束了主梁 向自由度，、和支座约束了主梁 向自由度。因此，和桥墩在 向为固定墩，和桥墩为滑动墩。整个有限元模型包含 个框架单元及 个 单元。支座弹簧桥台弹簧刚性构件4#墩1#桥台8#支座2#桥台5#支座6#支座7#支座4#支座3#支座2#支座1#支座1#墩2#墩3#墩横向 y竖向 z纵向 x图 有限元模型 第 期白玉磊，等：大应变 加固连续梁桥地震易损性分析 地震动记录选取选择地震动输入是进行桥梁地震易损性分析的前提，文献表明为了体现地震动的随机性，选择 条地震动记录能够达到一定精度。因此，根据“谱相容性原则”从美国太平洋地震研究中心（，）强震记录数据库选取震中距大于 的 条 类场地地震波，以地震动峰值加速度 作为地震动参数，图 给出了调幅后阻尼比为 的各条地震波的平均谱加速度曲线与 公路桥梁抗震设计规范中 类场 反应谱曲线的拟合情况。由图 可看出，所选取地震记录的均值谱与 类场的 地震谱曲线吻合较好，可用于算例桥梁的地震时程分析。为简化起见，仅研究在纵向地震作用下的桥梁结构响应，因此只考虑纵向地震动输入。20 条地震动反应谱20 条地震动反应谱的平均值-标准差20 条地震动反应谱的平均值+标准差20 条地震动反应谱的平均值规范 E1 反应谱0 1 2 3 4 5 6周期/s0.350.300.250.200.150.100.05伪谱加速度/g图 加速度反应谱 损伤指标的确定与破坏状态的划分在地震易损性分析中，文献指出桥梁构件的抗震能力可定义为当桥梁构件丧失某一特定功能时桥梁所达到的破坏状态。采用 软件对墩底截面进行弯矩 曲率分析，采用文献提出的划分极限状态值，得到桥墩 种不同损伤状态：轻微破坏、中等破坏、严重破坏、完全破坏所对应的曲率延性作为评判依据。虽然 加固可以明显地提高混凝土的变形能力，但其提供的是一种被动约束，即当混凝土膨胀时才会受到 的约束作用。因此，加固对于墩柱轻微和中等破坏状态指标影响较小，而对于严重和完全破坏状态指标提高幅度较大。此外，参考文献 ，选择活动支座的相对位移作为其损伤指标。加固前后桥墩和支座的各个极限状态值见表。概率地震需求模型基于 平台，将上述 条地震动记录分别输入到桥梁模型中进行 分析，每条地震动从 逐步调幅到 ，提取各个地震动作用下各个桥墩的最大截面曲率和支座的最大相对位移，根据式（）对这些散点进行线性回归分析，可得到其回归系数 值和 值，统计结果见表。表 桥梁构件损伤指标取值 破坏状态墩柱损伤指标 取值未加固 加固活动支座损伤指标 取值 轻微破坏中等破坏严重破坏完全破坏注：墩柱曲率延性 ，其中为截面曲率，为未加固墩柱初始屈服曲率；表示活动支座的水平位移。表 桥梁构件概率地震需求模型 构件名称拟合函数未加固模型 桥墩 未加固模型 支座 加固模型 桥墩 加固模型 支座 值得注意的是，算例中的桥梁为规则桥梁，在考察桥梁系统的地震易损性分析时，主要是针对 个墩柱与 个支座进行各个极限状态的超越概率分析。算例桥梁主要由纵向地震动作用控制，主梁的纵向地震惯性力由固定墩承担，滑动墩处仅存在滑动摩擦力，该墩柱不会因活动支座传递的水平剪力而进入屈服，总在弹性阶段工作，因此，无需对柱顶设置活动支座的桥墩进行易损性分析。墩和 墩由于刚度完全相同，其分担的地震作用水平剪力也一致，故其地震易损性结果也相同，仅需选择其一进行分析即可。考虑到纵向地震作用下可将主梁的运动视为刚体运动，活动支座的位移变形基本一致，其易损性分析结果相同，仅需选择 个活动支座进行分析。综上所述，对算例梁桥的地震易损性研究最终可转化为对 个固定墩和 个活动支座在纵向地震动作用下的易损性研究。因此，表 仅给出了墩柱与 支座的概率地震需求模型。桥梁易损性分析 构件易损性曲线结合桥梁有限元模型分析结果（表 和表），再利用式（）即可求得 桥墩和 支座的地震易损性曲线，分别见图、。图中横坐标表示地震动的大小，纵坐标表示地震作用下结构反应超越不同破坏状态的概率。哈 尔 滨 工 业 大 学 学 报 第 卷(a)轻微破坏(b)中等破坏地震动峰值加速度 Apg/g地震动峰值加速度 Apg/g(c)严重破坏(d)完全破坏地震动峰值加速度 Apg/g地震动峰值加速度 Apg/g未加固模型LRS FRP 加固模型未加固模型LRS FRP 加固模型未加固模型LRS FRP 加固模型1.00.80.60.40.2未加固模型LRS FRP 加固模型0 0.2 0.4 0.6 0.8 1.0 1.20 0.2 0.4 0.6 0.8 1.0 1.20 0.2 0.4 0.6 0.8 1.0 1.20 0.2 0.4 0.6 0.8 1.0 1.2中等破坏超越概率1.00.80.60.40.2轻微破坏超越概率1.00.80.60.40.2完全破坏超越概率1.00.80.60.40.2严重破坏超越概率图 桥墩易损性曲线 (c)严重破坏(d)完全破坏地震动峰值加速度 Apg/g地震动峰值加速度 Apg/g(a)轻微破坏(b)中等破坏地震动峰值加速度 Apg/g地震动峰值加速度 Apg/g0 0.2 0.4 0.6 0.8 1.0 1.20 0.2 0.4 0.6 0.8 1.0 1.20 0.2 0.4 0.6 0.8 1.0 1.20 0.2 0.4 0.6 0.8 1.0 1.2未加固模型LRS FRP 加固模型未加固模型LRS FRP 加固模型未加固模型LRS FRP 加固模型未加固模型LRS FRP 加固模型1.00.80.60.40.2中等破坏超越概率1.00.80.60.40.2轻微破坏超越概率1.00.80.60.40.2完全破坏超越概率1.00.80.60.40.2严重破坏超越概率图 支座易损性曲线 由图 可知，当 小于 时，桥墩轻微破坏概率小于，其他破坏等级概率趋于零，基本处于安全状态。当 为 时，未加固桥墩完全破坏即倒塌的概率达到了，存在倒塌的危险，第 期白玉磊，等：大应变 加固连续梁桥地震易损性分析经加固后，其倒塌概率降至了，加固效果较为明显。由图 可知，加固结构的支座易损性概率降低在 以内，表明采用 加固对支座的破坏概率降低不明显。系统易损性曲线桥梁整体中任何单个构件的失效都有可能导致全桥发生破坏，因此有必要从整体上评价桥梁的系统易损性。同时考虑桥墩、支座对桥梁系统易损性的影响，参考文献，采用一阶界限法计算桥梁系统失效概率的上界和下界，来逼近桥梁系统的易损性。（）（）（）（）式中（）为第个构件发生失效的概率，（）为桥梁系统失效概率，为桥梁系统构件数量。图 分别给出了未加固和 加固桥梁系统在各种破坏状态下易损性曲线的上界和下界。可看出对于“轻微破坏”与“中等破坏”，加固前后上下界基本接近，说明 加固对这两种上下界影响不大；对于“严重破坏”及“完全破坏”，加固桥梁的上下界小于未加固桥梁，尤其当地震强度较大时，这种差距更为突出；加固得到的桥梁系统易损性曲线上下界与上述单个墩柱的影响规律基本一致。以上研究结果表明作者课题组提出的简化模型可有效评估算例桥梁的抗震性能。此外，该简化模型还适用于模拟其他各类 加固圆形构件约束区混凝土材料的特性，可为 加固结构的抗震分析提供一定的分析基础。(a)轻微破坏(b)中等破坏地震动峰值加速度 Apg/g地震动峰值加速度 Apg/g1.00.80.60.40.2中等破坏超越概率1.00.80.60.40.2轻微破坏超越概率未加固模型-上界未加固模型-下界LRS FRP未加固模型-上界LRS FRP未加固模型-下界未加固模型-上界未加固模型-下界LRS FRP未加固模型-上界LRS FRP未加固模型-下界0 0.2 0.4 0.6 0.8 1.0 1.20 0.2 0.4 0.6 0.8 1.0 1.2(c)严重破坏(d)完全破坏地震动峰值加速度 Apg/g地震动峰值加速度 Apg/g1.00.80.60.40.2完全破坏超越概率1.00.80.60.40.2严重破坏超越概率未加固模型-上界未加固模型-下界LRS FRP未加固模型-上界LRS FRP未加固模型-下界未加固模型-上界未加固模型-下界LRS FRP未加固模型-上界LRS FRP未加固模型-下界0 0.2 0.4 0.6 0.8 1.0 1.20 0.2 0.4 0.6 0.8 1.0 1.2图 桥梁系统易损性曲线 结 论根据作者课题组提出的简化 约束混凝土往复轴压应力 应变模型，基于 平台，采用增量动力分析方法，对一典型三跨连续梁桥在 加固前后的地震易损性进行了分析，研究了 加固对桥梁抗震性能的影响，主要结论如下：）轻微破坏和中等破坏状态下桥墩基本处于未受损伤的线弹性状态，比较 加固前后两种状态桥墩的易损性曲线，二者基本吻合，说明 加固对桥梁轻微破坏和中等破坏的地震易损性影响不大。）严重破坏和完全破坏状态下，对比 加固桥墩的易损性曲线，在同一损伤状态下，未加固桥梁的损伤超越概率明显大于 加固的；当哈 尔 滨 工 业 大 学 学 报 第 卷为 时，未加固桥梁完全损伤即倒塌的概率达到了，存在倒塌的危险，经加固后，其倒塌概率降至了，加固效果较为明显。）对于“轻微损伤”与“中等损伤”，加固前后上、下界基本接近，说明 加固对这两种上、下界影响不大；对于“严重损伤”及“完全损伤”，加固结构的上、下界小于未加固结构，尤其当地震强度较大时，这种差距更为突出。）作者课题组提出的简化 约束混凝土往复轴压模型可准确地评估 加固桥梁的抗震性能，可为 约束混凝土的工程实际应用提供参考。参考文献 ，（）：，（）：滕锦光 新材料组合结构 土木工程学报，（）：，（）：，（），：，（）：（），：，：，（），（）：，（）：（），（）：（），（）：，（）：，（）：（），（）：（），：，：，（），：，（）：，（）：（）（）：（），（）：，：公路桥梁抗震设计规范：北京：人民交通出版社，：，：，李立峰，吴文朋，黄佳梅，等 地震作用下中等跨径 连续梁桥系统易损性研究 土木工程学报，（）：，（）：吴文朋，李立峰，王连华，等 基于 的高墩大跨桥梁地震易损性分析 地震工程与工程振动，（）：，（）：（编辑 赵丽莹）第 期白玉磊，等：大应变 加固连续梁桥地震易损性分析