考虑结构劣化和荷载历史的既有桥梁承载力更新
李全旺, 王草, 张龙
清华大学 土木工程系, 北京 100084

作者简介: 李全旺(1975—),男(汉),天津,副教授。Email:li_quanwang@mail.tsinghua.edu.cn

摘要

为提高桥梁承载力评估的准确性,该文基于结构时变可靠度理论,提出了利用桥梁已有信息(包括荷载历史、性能劣化规律等)对桥梁初评承载力进行更新的方法。将结构承载力的衰减描述成完全相关的随机过程,将荷载历史描述成不同时刻相互独立的随机过程,给出了同时考虑承载力衰减规律、荷载历史及其不确定性的既有桥梁承载力更新的显式公式,并通过具体算例验证了该方法的准确性。利用该文方法,对比了考虑与不考虑荷载历史时的桥梁初评与更新后的承载力,以及在更新时考虑与不考虑承载力衰减的承载力更新效果。结果表明: 历史荷载越大,对承载力更新的效果越明显; 但若未同时考虑承载力的衰减,会高估桥梁的实际承载力,造成评估结果的不安全。

关键词: 既有桥梁; 荷载历史; 结构劣化; 承载力; 评估
中图分类号:U445.7 文献标志码:A 文章编号:1000-0054(2015)01-0008-06
Updating for bearing capacity of existing bridges considering structural deterioration and loading history
Quanwang LI, Cao WANG, Long ZHANG
Department of Civil Engineering, Tsinghua University,Beijing 100084, China
Abstract

The load bearing capacity of existing bridges can be more realistically assessed using a method that assesses the capacity based on load history information using structural time-dependent reliability theory. The capacity degradation is described by a fully-correlated stochastic process and the loading history is described by a statistically independent stochastic process. An explicit formula is presented that takes into account the capacity deterioration and the uncertainties associated with the load history. Illustrative examples demonstrate the accuracy of this method. The method is then used to analyze the effects of loading history, with and without consideration of the resistance degradation, on the capacity assessment of existing bridges. Larger loading histories give a larger updated bridge bearing capacity. However, if the resistance degradation is ignored, the bearing capacity will be greatly overestimated, resulting in unsafe evaluation results.

Keyword: existing bridges; loading history; structural deterioration; bearing capacity; evaluation

近年来,中国现代化路网交通建设取得了长足进步。与此同时,大量的桥梁正承受着日益繁重的交通荷载[1],其安全问题已变得越来越突出。由于超重负载、材料劣化等因素,导致很多既有桥梁处于病害工作状态,这必然会对今后国民经济的发展埋下隐患[2]。为保证桥梁的使用安全状态,同时避免不必要的维护与更换,为既有桥梁的可靠性评估建立合理、准确的方法是非常必要的。

传统的评估方法主要有基于规范和基于概率理论2种途径[3,4]。基于规范的方法简便易行[5,6],但是由于不能考虑桥梁实际运营的情况,包括实际劣化的情况和实际荷载的情况,因而在业内广为诟病。基于概率的方法能够把与桥梁有关的信息合理地进行应用,从而得到更为合理准确的评估,因而基于概率的方法是目前这一领域的研究方向[7]。既有桥梁由于承载力衰减过程具有很大的不确定性[8],目前既有桥梁的承载力评估在准确性上还不够令人满意。

如何提高准确性,是当前研究的一个重点问题。李全旺等[9]指出,充分合理地利用已有信息是解决这一问题的重要途径。基于此,本文主要研究如何利用桥梁的已有信息(包括荷载历史、性能劣化规律等)对初评承载力进行更新。通过理论推导给出考虑结构劣化的承载力更新模型,以方便地在桥梁承载力更新时同时考虑结构承载力衰减和荷载历史信息,得到更为准确的结果。

1 结构时变可靠度

结构承载力极限状态方程可简单表示为[10]

Z=R-S. (1)

其中: R表示结构抗力; S表示服役周期内最大荷载效应。二者均为随机变量。当 Z≥0时表示结构安全, Z<0表示结构失效。考虑到 R S的时效性,在评估期(0, T]内,结构的时变可靠度 Ps(0, T)为

Ps(0, T) =P{ Z( t) =R( t) -S( t) >0, t∈(0, T]} . (2)

其中: R( t)表示 t时刻结构的承载力; S( t)表示 t时刻结构承受的荷载效应。

在实际分析中,可将评估期(0, T]分为 n个区间: (0, t1], ( t1, t2],…, ( tn-1, tn],对应每个区间的最大荷载效应分别记为 Si, i=1,2, …, n。 如果每个区间的持时足够小,则可用 R( ti)近似代替区间( ti-1, ti]内结构的承载力。时间段(0, T]内的承载力过程和荷载随机过程可分别表示为离散的随机变量{ R1, R2, …, Rn }和{ S1, S2, …, Sn},则结构在该评估期内的时变可靠度改写为[11]

Ps(0, T) =P{ R( t1) >S1 R( t2) >S2∩…∩ R( tn) >Sn} .(3)

进一步地,考虑结构承载力随时间的衰减,衰减规律通常表示为[12]

R( t) =R0· g( t) . (4)

其中: R0表示桥梁刚建成时的承载力; g( t)是衰减函数。式(3)可进一步写为[3]

Ps(0, T) =P[ R0· g( t1) -S( t1) >0]· P[ R0· g( t2) -S( t2) >0 |R0· g( t1) -S( t1) >0]·…· P[ R0· g( tn) -S( tn) >0 |R0· g( tn-1) -S( tn-1) >0] .(5)

根据式(5),给定初始抗力、退化规律和荷载过程的概率模型,即可得到结构在评估使用年限内安全服役的概率。 Mori和Ellingwood给出了荷载过程为平稳Poisson过程时式(5)的显式表达[13]; Li等[14]在文[13]的基础上将平稳荷载过程推广至更加一般的非平稳过程; 秦权等[15]基于随机过程理论给出了式(5)的计算方法。

2 既有桥梁承载力更新模型

荷载试验是获得桥梁承载力信息的最直接的手段。如果桥梁成功通过了试验荷载q*,则其承载力概率密度即可从试验前的fR(r)更新到 fR'(r),即

fR'(r)=P(r|r>q*)= fR(r)1-FR(q*),r≥q*. (6)

其中FR ()表示承载力的累积密度函数。

桥梁的承载历史,等同于反复多次的荷载试验。如果将整个荷载历史中最大荷载效应的概率密度函数和累积概率密度分别表示为fS(s)和FS(s),则该桥成功服役到现在的概率为

Ps= -+FS(r)fR(r) dr. (7)

考虑到q*的随机性,则式(6)可改写为[16

] fR'(r)=FS(r)fR(r)-+FS(r)fR(r)dr. (8)

式(8)即是学界广为接受的基于荷载历史的桥梁承载力更新模型[16]。但这个模型没有考虑到桥梁的承载力是随着时间而衰减的,即在历史的某一时刻桥梁成功承受了某一荷载,不代表当前该桥梁依然能够承受该荷载,因此会高估当前桥梁的承载能力。 Stewart[11]曾注意到这个问题,并利用 Monte Carlo模拟的方法研究了考虑承载力劣化的车载历史对结构承载力更新的效果,但该法计算代价很大,难以推广应用。

3 考虑结构劣化和荷载历史的承载力更新模型

将桥梁的服役时间(从建成通车到现在)划分为n个时间区间,当区间划分足够小时,假设在该区间内承载力大小保持不变,承载力衰减函数如图1所示。

本文将承载力衰减看作完全相关的随机过程,在给定衰减曲线总体形状的前提下(如线性衰减、指数衰减等),若已知当前时刻(tn)的衰减函数值 g(tn),便可以确定任意时刻(ti)的衰减函数值 g(ti)。

每个时间区间的荷载历史均可视为对桥梁的一次荷载试验。对应(ti-1,ti]区间,其最大荷载效应的累积概率函数表示为FS,i,i=1, 2,3,…, n。 若当前时刻(tn)桥梁承载力的先验概率密度函数为 fR(r),则利用n次荷载试验结果对承载力评估进行更新,后验抗力密度函数 fR*(r)

fR*(r)= fR(r)i=1nFS,ig(ti)g(tn)r-+fR(r)i=1nFS,igtig(tn)rdr. (9)

其中g(ti)表示承载力衰减函数在(ti-1,ti]区间的取值。

为考虑承载力衰减函数的随机性,引入当前时刻衰减函数值g(tn)的概率密度函数fG(g)。若新建桥梁初始抗力为一确定值r0, 则由于退化过程假设为完全相关,fG(g)的不确定性体现为当前先验承载力评估的不确定性,因此fG(g)可表示为

fG(g)=r0·fR(r0·g). (10)

衰减函数随机过程模型如图2所示。

图2 结构抗力衰减的随机性描述

由式(10),考虑衰减过程的随机性,式(9)可改为

fR*(r)=fR(r)·i=1nFS,i(r0·g(ti))-+fR(r)·i=1nFS,i(r0·g(ti))dr. (11)

其中: g(tn)= rr0,且g(t1), g(t2), …, g(tn-1)可由g(tn)确定。

根据式(11),便可充分利用既有桥梁的既有信息,包括: 结构初始抗力(r0)、当前时刻承载力初评结果(r0fG(g))、历史车辆荷载累积概率函数(FS),对承载力进行更新。

若假设承载力不随时间衰减,即g(ti) =1.0, 代入式(11),则可得到式(8)。因此,式(8)是本文公式的一个特例,本文公式具有更好的普适性。

另外需要指出的是,本文所提出的承载力更新方法忽略荷载反复作用造成的结构损伤累积,不适用于车辆荷载造成桥梁疲劳破坏的情况。

4 算例分析

下面将结合算例说明本文所提出的更新方法的正确性。

某钢筋混凝土梁式桥已使用20 a,其横断面如图3所示。主梁采用T型截面,计算跨径15.5 m, T梁截面如图4所示[17]

梁体采用C30混凝土,纵筋为HRB335。根据规范[18]计算得截面的抗弯承载力为2 804 kN·m, 看作结构的初始承载力r0(忽略r0的随机性)。依据《公路桥梁承载能力检测评定规程》[6],考虑桥梁的实际劣化状态,评出当前状态的承载能力为 2 300 kN·m。 又考虑到桥梁承载力的不确定性,例如: 国内外关于钢筋混凝土梁承载力的统计分析表明[19],承载力服从对数正态分布,变异系数接近 0.15。 因此,在本研究中,当前状态下T型截面梁跨中抗弯承载力的初评结果为: 服从均值为 2 300 kN·m, 变异系数为0.15的对数正态分布,概率密度函数fR图5

图5 当前时刻先验承载力和年最大车载效应分布

考虑到结构自重引起的恒荷载D, D=702 kN·m, 式(11)修正为

fR*(r)=fR(r)·i=1nFS,i(r0·g(ti)-D)-+fR(r)·i=1nFS,i(r0·g(ti)-D)dr. (12)

该桥所在道路上安装了动态称重系统。利用测得的60 d车载数据,模拟出每年车载引起的梁跨中截面最大弯矩均服从极值I型分布,均值和标准差分别为1 020 kN·m和270 kN·m, 其累积密度函数(FS,1)和概率密度函数(fS,1)可分别表示为

FS,1(x)= exp -exp-x-898211, (13)

fS,1(x)=FS,1(x)· 1211exp-x-898211. (14)

其中fs,1的形状如图5所示。

需要指出,实际车载效应往往是多峰分布[1],难以用单峰模型准确描述。最新的研究指出应着眼于荷载效应尾部分布特征,并建议用GEV(generalized extreme value)或者GPD(generalized Paretodistribution)模型描述[20,21],但过程复杂。本文采用的极值I型应用简单,虽然不够准确,但不妨碍本文方法的展示,且不影响定性结论。假定桥梁承载力在过去的20 a线性衰减,即式(4)中的 g( t)取关于时间 t的一次函数。

利用式(12)对当前T型截面梁的承载力进行更新。更新后的承载力分布如图6所示,均值为2 420 kN·m, 标准差为314 kN·m, 形状接近对数正态分布。为了进行比较,承载力初评结果也表示在了图6中。图中可以看出,荷载历史对桥梁承载力有着明显的更新效果,均值由2 300 kN·m增加到2 420 kN·m, 提高了5.2%; 标准差也略有减小,由345 kN·m减少到314 kN·m。

为了验证公式(12)的正确性,本文采用Monte Carlo数值模拟的方法考虑荷载历史和结构承载力衰减,对桥梁初评承载力进行了更新,具体步骤为:

Step 1 按均值为2 300 kN·m, 标准差 345 kN·m 的对数正态分布随机生成一个承载力样本 r;

Step 2 按均值为1 020 kN·m, 标准差 270 kN·m 的极值I型分布随机生成相互独立的荷载系列 s1, s2,…, s20;

Step 3 由样本 r和初始承载力值 r0( r0=2 804 kN·m),按照线性衰减规律得到服役的20 a中每一年的承载力值 r1, r2,…, r20;

Step 4 对于 i=1, 2, …,20, 若 ri-D<si, 表示结构在第 i年失效,舍弃此次模拟,重新回到Step 1; 若对于所有的 i, 均有 ri-D si, 则表示结构成功服役20年,记录下 r20;

Step 5 多次重复Step 1至Step 4,得到 r20的统计规律,即为更新后的承载力概率分布。

通过100万次的Monte Carlo数值模拟,得到T型截面梁承载力的概率分布,统计直方图也表示在图6中。统计直方图与公式(12)计算得到的曲线完全重合,证明了本文所提出的公式的正确性。

图6 公式计算与模拟结果的比较

为了研究荷载历史对承载力更新的影响,本文假设出另外2个年最大荷载效应的分布: 荷载1为均值为1 100 kN·m, 标准差为291 kN·m的极值I型; 荷载2为均值为900 kN·m, 标准差为 238 kN·m 的极值I型。再次利用式(12),分别得到根据荷载1和荷载2更新后的T型梁承载力的概率分布,并与实测荷载的更新结果进行了对比,如图7所示。

图7 不同车辆荷载对桥梁承载力更新的影响

由图中可以看出: 荷载历史对桥梁承载力更新有重要的影响,考虑3个荷载历史都对当前承载力有显著的更新效果。荷载越大,更新效果越明显(均值变大,标准差变小),例如: 3个荷载模型,由大到小分别为荷载1、 实测荷载、荷载2; 对应的更新后承载力的均值分别为: 2 468 kN·m, 2 420 kN·m、 2 364 kN·m; 标准差分别为 310 kN·m、 314 kN·m、 324 kN·m。

下面讨论在利用荷载历史对承载力更新时,是否同时考虑结构承载力衰减对更新结果的影响。对于只考虑荷载历史、不考虑结构承载力的衰减(即假设 g( ti)=1.0)的情况,Hall等[16]在1988年已经给出了承载力更新公式,式(8)。在式(8)中, fR表示承载力的初评概率密度函数,并假设在评估使用年限内承载力保持不变,对应本文算例, fR服从均值为2 300 kN·m, 变异系数为0 .15的对数正态分布; FS对应的不是每年最大荷载效应的最大值,而是整个荷载历史中荷载效应最大值的累积概率函数,对应本文算例,假设在过去的20 a内,每年的荷载效应最大值独立同分布,则

FS( x) = [FS,1(x)]20 . (15)

代入式(8),则可得到忽略结构承载力衰减的更新结果,并与考虑结构承载力衰减的更新结果进行比较,如图8所示。

图8 考虑结构劣化对桥梁承载力更新的影响

图8表明,无论是否考虑结构承载力衰减,荷载历史都对桥梁承载力评估有明显的提高作用。在利用荷载历史对承载力进行更新时,若忽略结构承载力随时间的衰减,更新后的承载力水平会比实际情况有所高估(均值由2 420 kN·m提高到 2 547 kN·m, 提高5.2%)。这是因为,公式(8)假设结构在建成初期的承载力与当前时刻相同(即: 不考虑承载力的衰减),对结构早期的承载力明显低估,因而夸大了结构早期荷载历史对承载力的更新效果,造成了更新结果的高估。这说明: 在利用桥梁的荷载历史对其承载力进行更新时,有必要同时考虑结构承载力的衰减,否则会高估既有桥梁的实际承载能力,低估桥梁的服役风险,造成评估结果的不安全。

5 结论

本文提出了依据既有桥梁的荷载历史对其承载力初评结果进行更新的方法,在考虑桥梁荷载历史的同时考虑承载力的衰减。随后利用该方法讨论了荷载历史和是否考虑结构承载力衰减对桥梁承载力评估的影响,得到如下主要结论:

1) 提出了利用既有桥梁荷载历史,同时考虑结构承载力衰减的既有桥梁承载力更新方法,并利用算例证明了该方法的准确性和实用性;

2) 考虑荷载历史能提高既有桥梁承载力评估的结果,荷载效应越大,对评估结果的提高作用越明显;

3) 在利用荷载历史对桥梁承载力更新时,必须同时考虑结构承载力的衰减,否则会高估桥梁的实际的承载力,造成评估结果的不安全。

The authors have declared that no competing interests exist.

参考文献
[1] 梅刚, 秦权, 林道锦. 公路桥梁车辆荷载的双峰分布概率模型[J]. 清华大学学报: 自然科学版, 2003, 43(10): 1394-1396.
MEI Gang, QIN Quan, LIN Daojin. Bi-modal probabilistic model of highway and bridge vehicle loads[J]. Journal of Tsinghua University: Science and Technology, 2003, 43(10): 1394-1396. (in Chinese) [本文引用:2] [CJCR: 0.609]
[2] 钟曙亮. 公路既有桥梁承载力评定方法研究 [D]. 成都: 西南交通大学, 2010.
ZHONG Shuliang. Rating Method for Resistance of Highway Existing Bridges [D]. Chengdu: Southwest Jiaotong University, 2010. (in Chinese) [本文引用:1]
[3] 张宇贻, 秦权. 钢筋混凝土桥梁构件的时变可靠度分析[J]. 清华大学学报: 自然科学版, 2001, 41(12): 65-67.
ZHANG Yuyi, QIN Quan. Time dependent reliability analysis of deteriorating RC bridge girders[J]. Journal of Tsinghua University: Science and Technology, 2001, 41(12): 65-67. (in Chinese) [本文引用:2] [CJCR: 0.609]
[4] Enright M P, Frangopol D M. Service-life prediction of deteriorating concrete bridges[J]. Journal of Structural Engineering, 1998, 124(3), 309-317. [本文引用:1]
[5] JTG H11-2004. 公路桥涵养护规范 [S]. 北京: 人民交通出版社, 2004.
JTG H11-004, Code for Maintenance of Highway Bridges and Culvers[S]. Beijing: China Communications Press, 2004. (in Chinese) [本文引用:1]
[6] JTG/T J21-2011. 公路桥梁承载能力检测评定规程 [S]. 北京: 人民交通出版社, 2011.
JTG/T J21-2011. Specification for Inspection and Evaluation of Load-bearing Capacity of Highway Bridges [S]. Beijing: China Communications Press, 2011. (in Chinese) [本文引用:2]
[7] Ellingwood B R. Reliability-based condition assessment and LRFD for existing structures[J]. Structural Safety, 1996, 18(2): 67-80. [本文引用:1] [JCR: 2.391]
[8] Ellingwood B R. Risk-informed condition assessment of civil infrastructure: State of practice and research issues[J]. Structure and Infrastructure Engineering, 2005, 1: 7-18. [本文引用:1] [JCR: 0.954]
[9] 李全旺, 李春前, 孙健康. 基于结构可靠度理论的既有桥梁承载能力评估[J]. 工程力学, 2010, 27(II): 142-151.
LI Quanwang, LI Chunqian, SUN Jiankang. Capacity evaluation of existing bridges based on structural reliability theory[J]. Engineering Mechanics, 2010, 27(II): 142-151(in Chinese) [本文引用:1] [CJCR: 0.645]
[10] Ang A H S, Tang W H. Probability Concepts in Engineering Planning and Design[M]. 1st ed. New York, USA: Wiley, 1975. [本文引用:1]
[11] Stewart M G, Val D V. Role of load history in reliability-based decision analysis of aging bridges[J]. Journal of Structural Engineering, 2004, 125(7): 776-783. [本文引用:2]
[12] Enright M P, Frangopol D M. Condition prediction of deteriorating concrete bridges using Bayesian updating[J]. Journal of Structural Engineering, 1999, 125(10), 1118-1125. [本文引用:1]
[13] Mori Y, Ellingwood B R. Reliability-based service-life assessment of aging concrete structures[J]. Journal of Structural Engineering, 1993, 119(5): 1600-1621. [本文引用:2]
[14] LI Quanwang, WANG Cao, Ellingwood B R. Time dependent reliability of aging structures in the presence of non-stationary loads and degradation[J]. Structural Safety, 52(2015): 132-141. [本文引用:1] [JCR: 2.391]
[15] 秦权, 杨小刚. 退化结构时变可靠度分析[J]. 清华大学学报: 自然科学版, 2005, 45(6): 733-736.
QIN Quan, YANG Xiaogang. Time dependent reliability analysis of deteriorating structures[J]. Journal of Tsinghua University: Science and Technology, 2005, 45(6): 733-736. [本文引用:1] [CJCR: 0.609]
[16] Hall W B. Reliability of service-proven structures[J]. Journal of Structural Engineering, 1988, 114(3): 608-624. [本文引用:2]
[17] 袁伦一, 鲍卫刚. 《公路钢筋混凝土及预应力混凝土桥涵设计规范》条文应用算例 [M]. 北京: 人民交通出版社. 2005.
YUAN Lunyi, BAO Weigang. Explanations for the Code for Design of Highway Reinforced Concrete and Prestressed Concrete Bridges and Culverts [M]. Beijing: China Communications Press, 2005. (in Chinese) [本文引用:1]
[18] JTG D62-2004 D62-2004. 公路钢筋混凝土及预应力混凝土桥涵设计规范 [S]. 北京: 人民交通出版社. 2004.
JTG D62-2004 D62-2004. Code for Design of Highway Reinforced Concrete and Prestressed Concrete Bridges and Culverts [S]. Beijing: China Communications Press, 2004. (in Chinese) [本文引用:1]
[19] Ellingwood B R. Development of a Probability based Load Criterion for American National Standard A58: Building Code Requirements for Minimum Design Loads in Buildings and Other Structures[M]. Washington DC, USA: US Department of Commerce, National Bureau of Standards, 1980. [本文引用:1]
[20] O'Brien E J, Bordallo-Ruiz A, Enright B. Lifetime maximum load effects on short-span bridges subject to growing traffic volumes[J]. Structural Safety, 2014, 50: 113-122. [本文引用:1] [JCR: 2.391]
[21] 李植淮, 李春前, 孙健康, . 基于 GPD 模型的车辆荷载效应极值估计[J]. 工程力学, 2012, 6(增刊I): 166-171.
LI Zhihuai, LI Chunqian, SUN Jiankang, et al. Estimation of extreme vehicle load effect based on GPD model[J]. Engineering Mechanics, 2012, 6(suppl I): 166-171. [本文引用:1] [CJCR: 0.645]