箍筋锈蚀约束混凝土力学性能退化数值模拟.pdf

四川建 筑科 学研究 S i c h u a n B u i l d i n g S c i e n c e 第 3 9卷第 5期 2 0 1 3年 1 O月 箍筋锈蚀约束混凝土力学性能退化数值模拟 张 雯， 牛荻涛， 李 强 ( 西安建筑科技大学土木工程学院， 陕西 西安7 1 0 0 5 5 ) 摘要： 为研究箍筋锈蚀对约束混凝土的影响 ， 运用 A B A Q U S有限元软件对不同锈蚀率箍筋约束混凝土轴心受压 构件进行分析 ， 并引入配箍特征值 A 表征箍筋不同锈蚀率 7 7 对构件的影响。分析结果表明： 随着箍筋锈蚀率的增 大， 约束混凝土的承载力和变形能力均有不同程度的下降。通过对数据的回归分析， 建立了箍筋锈蚀后约束混凝 土应 力 应变模型。 关键词： 箍筋锈蚀； 有限元； 配箍特征值 ； 应力 应变模型 中图分类号 ： T U 3 1 2 文献标 志码 ： A 文章编号 ： 1 0 0 81 9 3 3 ( 2 0 1 3 ) 0 5 0 9 4 0 4 S t i r r u p s c o r r o s i o n o f c o nc r e t e c o l u m n wi t h a x i a l c o mp r e s s i o n f o r n u me r i c a l s i mu l a t i o n Z H A N G We n ， N I U D i t a o ， L I Q i a n g ( X i a n U n i v e r s i t y o f A r c h i t e c t u r e T e c h n o l o g y ， X i a n 7 1 0 0 5 5 ， C h i n a ) A b s t r a c t ： C o n s i d e r i n g t h e e ff e c t o f s t i r r u p c o r r o s i o n o n t h e p e rf o r ma n c e o f c o n f i n e d c o n c r e t e ， A B A Q U S fi n i t e e l e me n t s o ft w a r e i s u s e d t o r e a r c h o n d ife ren t C O rrO s i o n r a t e o fa x i a l c o mp r e s s i o n c o l u mn I n t r o d u c e t h e s t i r r u p e i g e n v a l u e入 t o d i ff e r e n t c o r r o s i o n r a t e o f s t i r r u p s 田 。o n t h e c o mp o n e n t s T h e res m s h o w t h a t ， wi t h t h e s t i r r u p s c o r r o s i o n r a t e i n c rea s e s ， t h e b e a r i n g c a p a c i t y a n d d e f o rm a t i o n c a p a c i t y o f c o n f i n e d c o n c ret e h a v e d i ff e r e n t d e g r e e s of d e c l i n e B y t h e r e g r e s s i o n a n aly s i s O i l t h e s i mu l a t i o n d a t a ， s ti r r u p s c o r r o s i o n c o n fi n e d c o n c r e t e s t res s s t r a i n mo d e l h a s b e e n e s t a b l i s h e d Ke y wo r d s： s t i r r u p c o r r o s i o n ； fi n i t e e l e me n t ； s t i r r u p c h a r a c t e r i s t i c v Mu e ； s t r e s s s t r a i n mo d e l 0 引 言 钢筋锈蚀是导致钢筋混凝土构件性能退化的主 要原因之一。钢筋锈蚀后 ， 其力学性 能会发生改 变， 强度和延性均降低， 从而造成服役结构构件力学 性能退化 ， 严重时可能造成结构提前失效甚至倒塌。 近年来 ， 国内外学者对锈后钢筋混凝土构件 承载力 退化进行了研究 ， 研究方法主要采用有限元分析 和试验研究两种。运 用有 限元方 法分析锈蚀构 件力学性能退化 ， 可 以克服快 速腐蚀试验与实际工 程构件 自然腐蚀间差异及试验周期长等缺陷， 因此 ， 有限元法成为解决钢筋混凝 土结构耐久性分析的有 效手段 。 目前 ， 国内外学者对纵筋锈蚀后构件的力学性 能退化研究较多 加 J ， 对于箍筋锈蚀相对较少。然 而， 从钢筋混凝土耐久性试验研究和现场检测发现， 在 同等条件下 ， 箍筋锈蚀程度 比纵 向钢筋锈蚀程度 收稿 日期 ： 2 0 1 2 - 0 3 - 2 7 作者简介： 张雯( 1 9 8 7一 ) ， 女， 江苏苏州人， 硕士研究生， 研究方 向： 混凝土结构耐久性。 基金项目： 国家杰出青年科学基金( 5 0 7 2 5 8 2 4 ) ； 陕西省“ 1 3 1 1 5 ” 科技 创 新工程项 目( 2 O 1 O z D K G_ 5 5 ) E ma i l ： 5l 6 9 9 l4 6 7 6 q q c o rn 大得多 J ， 有时当纵筋锈蚀量仅为 5 1 0 时 ， 箍 筋的锈蚀 已经很严重， 甚至有锈断现象 J 。箍筋锈 蚀后： 截面面积减小， 其屈服强度、 弹性模量和延 性均降低； 对混凝土约束作用降低； 保护层发生 锈涨开裂； 箍筋与混凝土间粘绪陛能退化； 构件 破坏形式变为脆性破坏。箍筋锈蚀 对混凝土结构的 影响很大 ， 对箍筋锈蚀后约束混凝 土力学性能退化 研究尤为重要 。本文采用 A B A Q U S有 限元软件 ， 对 不同箍筋锈蚀 程度的轴 心受压混凝土短 柱进行分 析 ， 研究锈蚀箍筋约束混凝土的力学性能退化。研 究成果不仅可以为在役结构 的检 贝 4 鉴定提供依据 ， 而且可 以为新结构的设计做指导。 1 有限元模型 1 1 模型概况 1 1 1 构件 概 况 为考虑箍筋在约束混凝土中的作用 ， 模型建立 时不考虑保 护层 作用 ， 采 用 尺寸 为 1 5 0 m m 1 5 0 m m 4 5 0 m m的箍筋约束混凝土短柱， 混凝土强度 为 2 5 6 4 MP a ， 构件尺寸如图 1 所示。 1 1 2加载与边界条件 本次模拟将受压柱柱底设 为 固接 ， 柱顶 自由。 2 0 1 3 N o 5 张 雯， 等 ： 箍筋锈蚀约束混凝土力学性能退化数值模拟 9 5 6 O 图 1 构件配筋及 尺寸 Fi g 1 Re i n f o r c e me n t a n d s i z e o f c o mp o n e n t 由于模型直接加载时容易在局部节点上产生应力集 中， 导致求解失败 ， 为 了使 昆凝土构件均匀受力 ， 建 立 了刚性体垫块。垫块截面取与柱截面尺寸同， 厚 6 0 mm， 模拟过程 中将荷载施加在刚性体垫块上 ， 通 过垫块进一步传递到构件上 。 1 1 3实体模 型 及 网格 划分 网格划分 - 6 9 的疏密直接影 响到所求结果的准 确性 ， 网格划分要求具有足够多的单元 和节点 ， 以满 足求解结果的精确性。考虑本文模型尺寸， 模型中 混凝土所选用的网格划分尺寸定为 2 5 m m 2 5 m i l l 2 5 m m。构件具体网格划分情况见图 2 。 图 2 构 件网格划分 F i g 2 Gr i d o f c o m p o n e n t 1 1 4单元选择 本文采用 AB A Q U S有 限元软件模 拟分析锈蚀 箍筋约束混凝土性能的变化， 混凝土单元采用三维 实体单元 C 3 D 8 R 。钢筋单元采用桁架单元， 即空 间二节点单元 T 3 D 2 ， 用 e m b e d d e d 方式将普通钢筋 嵌人混凝土 中。 1 2 混凝土及钢筋本构 混凝土单轴受力时 ， 应力一 应变关系反映了混凝 土受力全过程的重要力学特征， 是进行混凝土构件 应力分析、 建立承载力和变形计算理论的必要依据， 也是利用计算机进行有限元非线性分析的基础。 1 2 1混凝土本构 有限元分析软件 A B A Q U S中共提供 了三种不 同的混凝 土本 构模型 6 】 ： 弥散裂 纹混凝 土模 型 、 混 凝土损伤塑性模型 、 混凝土裂纹模型 。本文在建模 过程中综合考虑各因素， 采用如下混凝土本构模型 ， 混凝土在单轴加载下拉伸和压缩应力一 应变 曲线分 别如图 3和图 4所示。 图3 馄凝土单轴受压应力一 应变曲线 F i g 3 Un i a x i a i c o mp r e s s i v e s t r e ss- s t r a i n c u r v e o f c o n c r e t e 图4 混凝土单轴受拉应力一 应变曲线 Fi g 4 Un i a x i a l t e n s i l e s t r e ss- s t r a i n c u r v e o f c o n c r e t e 1 2 2钢 筋本构 在研究钢筋混凝土结构性能时 ， 钢筋应力 应变 关 系常采用 的模 型有 以下 三种【 6 】 ： 理想 弹塑 性模 型 、 弹塑性模型 和刚塑性模型。本文在建立有 限元 模 型时 ， 钢筋采用各 向同性塑性模型即理想 弹塑性 本构模型 ， 其应力一 应变关系如图 5所示 。钢筋初始 弹性模量为 2 11 0 MP a ， 屈服强度取为 2 8 0 MP a ， 泊松 比为 0 3 。 图 5钢筋受拉应 力一 应 变 曲线 Fi g 5 Te n s i l e s t r e s s - s t r a i n c u r v e o f s t e e l 1 3 箍筋锈蚀后性能折减 文献 7 对比分析了箍筋参数中影响约束混凝 土性能的因素： 配箍特征值 A 对约束混凝土性能的 影响较大 ， 箍筋 间距 次之， 箍筋构造 和形式影 响不 9 6 四川建筑科学研究 第 3 9卷 大。配箍特征值 A 也在一定程度上反映箍筋间距 影响， 因此， 配箍特征值 A 可作为箍筋对约束混凝 土作用的综合参数。本文模型中， 构件其余条件不 变 ， 仅箍筋发生锈蚀 ， 因此 ， 以配箍特征值 A 作为主 要参数 。 1 3 1 锈 后箍 筋 强度折 减 许多学 者 进行 了大 量锈 蚀 钢 筋 力 学 性 能试 验 引， 结果表 明， 钢筋的屈服强度 、 极限抗拉强度及 延伸率都 随钢筋锈蚀程度增大而降低。文 献 1 0 对现有的锈蚀钢筋屈服强度和极限强度公式进行了 分析验证， 给出了锈蚀钢筋屈服拉力和极限拉力的 计算式 ， 本文 中的锈蚀 箍筋后 力学 性 能采用 文献 1 0 中方法折减。 1 3 2 锈后箍筋截面积折减 一 般大气环境下 ， 光圆钢筋表面锈蚀较均匀 ， 本 文视箍筋锈蚀为均匀锈蚀 ， 采用箍筋截面损失率 来表达箍筋锈蚀程度 ， 见下式 ： A =A 0 ( 1一叼 ) ( 1 ) 式中A 为箍筋锈蚀前截面面积 ； A为箍筋锈蚀后 截面面积。 2 数值模型计算结果和分析 数值模拟分析得到的锈蚀箍筋约束混凝土应力 一 应变曲线 ， 与以往完好构件约束混凝 土应力一 应变 曲线形状类似 ， 主要包括上升段和下降段两个部分。 约束混凝土应力一 应变 曲线形状受箍筋锈蚀量影响 较大： 随着箍筋锈蚀率增大 ， 约束混凝 土应力一 应变 关系曲线的峰值点呈下降趋 势 ， 上升段向右下方偏 移 ， 下降段变得陡直。 随着箍筋锈蚀率 叼 的增大， 约束混凝土峰值应 变和峰值应力值都随之下降 ， 其 中箍筋锈蚀对 约束 混凝土峰值应力的影响要略大于峰值应变。当箍筋 锈蚀 1 0 时 ， 约束混凝土峰值应变减小 0 7 3 ， 峰 值应力减小5 7 ； 箍筋锈蚀 2 0 时， 约束混凝土峰 值应变减小4 6 ， 峰值应力减小 5 7 ； 箍筋锈蚀 3 0 时 ， 约束混凝土峰值应变减小 2 7 ， 峰值应力 减小5 8 ； 箍筋锈蚀4 0 时， 约束混凝土峰值应变 减小 5 7 ， 峰值应 力减 小 5 7 ； 箍 筋锈 蚀 5 0 时， 约束混凝土峰值应变减小 0 9 ， 峰值应力减小 5 7 。 2 1 锈蚀箍筋约束混凝土承载力 体积配箍率P与配箍特征值 入 之间的换算关 系： A ：乎 ( 2 ) J P： ( 3 ) j 式 中P为构件体积 配箍率 ； f 为平 面内箍筋 净间 距； A 为核心区混凝土面积； n为平面内箍筋肢数 ； s 为箍筋 间距 ； A 为单肢箍筋 截面积。钢筋锈蚀后 ， 其强度折减 为： = ( 1一1 0 7 7 q ) ( 4 ) 式中 为锈蚀钢筋的名义屈服强度。 可推导配箍筋锈蚀后 的配箍特征值 A 与箍筋 锈蚀率 7 7 。 之间的关系 ： A ， ： 二 l_ ：( 1 一 1 0 一 一 叼 。 ) A v 一 L lI u 卵 B ，v ( 5 ) 峰值应力相对减小值( 。 一 。 ) 。 与配箍特征 值 A 的关系如图 6所示 ， 两者大致呈线性关系。随 着锈蚀率 7 7 。 的增大， 配箍特征值 A 的值随之降低， 约束混凝土的峰值应力 。 值也逐渐减小 ， 以上数据 均列于表 1 。 图 6( 一 ) 与 A 的关系 F i g 6 r e l a t i o n s h i p b e t w e e n( 一 ) 。 a n d A ， 表 1 不 同锈蚀 率 田 。 对应 的峰值应力 和配箍特征值 A T a b l e 1 P e a k s t res s | a n d s t i r r u p e i g e n v a l u e At und e r di ffe r e nt c o r r o s i o n r a t e s 根据表 1中数据 ， ( 一 。 ) 。 与 A 的关 系示 于图 6， 回归数据得到以下关系式 ： ( 。一 。 ) 。=3 3 1 5 A 一0 2 5 9 2 ( 6 ) 即 =( 0 7 4 0 8+3 3 1 5 A ) 。 ( 7 ) 式中 。 为素混凝土轴心抗压强度。 再将式 ( 7 ) 代人式 ( 9 ) ： ( o 7 4 0 8+ ( 8 ) 2 2 锈蚀箍筋约束混凝土变形 峰值应力相对减小值 。 一 厂 c。 ) 与 A 的关系 2 0 1 3 N 。 5 张雯， 等 ： 箍筋锈蚀约束混凝土力学性能退化数值模拟 9 7 如图 7所示 ， 两者大致呈线性关 系。在应力一 应变曲 线下降段中， 应 力值 为峰值应力 5 0 的点 ， 所对应 的应变 6 为极 限应 变， 8 0 5 6 。 为混凝土应变延性 系数 x ， 以上数值均列于表 2 。 0 图 7( s 一8 ) 6 与 A 的关 系 F i g 7 Re l a t i o n s h i p b e t we e n ( s c c s c o ) e 。 a n d A ， 表 2 特定点 应变 值及应变延性 系数 Ta b l e 2 Va l u e o f s t r a i n a n d i t s d u c t i l i t y f a c t o r o f s o me s p e c i fi c p o i n ts 由表2中数据可知， 随锈蚀率 增大， 约束混 凝土峰值应 变 g 。 减 小 ， 应 变延性 系数 也减小。 根据表 2中数据 ， 回归 ( 6 。 一 。 。 ) 6 与 A 的关 系 得到以下表达式 ： ( 占 。一 。 。 ) 6 。=2 2 1 5 6 A +0 1 3 8 4 ( 9 ) 即 。=( 1 1 3 8 4+2 2 1 5 7 6 A ) 占 。 。 ( 1 O ) 式中 。 。 为素混凝土峰值应变值。 o 5=( 0 9 5 3 3+6 6 4 7 1 A ) s 。 。 ( 1 1 ) 再将式( 7 ) 代人式( 1 2 ) 、 ( 1 3 ) 。 s 。 。=( 1 1 3 8 4+ ) s ( 1 2 ) 3 箍筋锈蚀约束混凝土本构模型 3 1 模型选取 箍筋锈蚀后， 约束混凝土应力一 应变本构关系曲 线仍然分为上升段下降段两部分。因此 ， 本文选用 有理分式形式 ， 通过拟合比较发现， 箍筋锈蚀约束混 凝土本构模型可用过镇海模 型 J 。 f Y= +( 32 a ) x +( 一2 ) l 【 ) ， 式中 和 分别为约束混凝土应力一 应变曲线上 升段和下降段形状控制参数； 横坐标 取为 = 6 ，纵坐标 Y取为 Y= ， 其 中， 。 为混凝 土的单 轴抗压强度 ， 。 为与 。 相对应 的 昆凝 土峰值压应 变。 3 2 模型参数取值 用 O r i g i n拟合数据 ， 得到锈蚀 箍筋约束 混凝 土 应力 应变全曲线形状控制参数及相关系数平方值 ， 详细数据见表 3 。 表 3 应力一 应变 曲线拟合结果 Ta b l e 3 F i t t i n g r e s u l t o f s t r e s s - s t r a i n c ur v e 由表 3数据看出， 本文建议本构模型拟合得 到 曲线与拟合结果得到曲线吻合较好。对方程参数 数据进行回归分析。 。 = 2 6 7 0 4 A +0 4 6 6 ( 1 3 ) 拟合相关系数为 R = 0 9 6 8 6 。 ： 3 0 1 4 3 e - 6 +1 2 0 7 ( 1 4 ) |f ) o 拟合相关系数为 R = 0 9 8 3 。 式中 。 为完好箍筋约束混凝土受压应力一 应变曲 线上升段形状参数； 。 为完好箍筋约束混凝土受压 应力- 应变曲线下降段形状参数。 将式 ( 7 ) 代人式( 1 5 ) 和式( 1 6 ) ： o L 。= +0 ( 1 54 0 0 1 3 ) + = 3 0 1 4 3 e 一 而丽丽+1 2 0 7( 1 6 ) po 4 结论 基于 A B A Q U S非线性有限元模拟， 分析箍筋锈 蚀后约束混凝土性能变化， 得到如下结论： 1 ) 箍筋锈蚀后 ， 其强度、 延性 、 截面积及对混凝 土约束作用均降低， 随着箍筋锈蚀率的增大， 约束混 凝土应力一 应变曲线中峰值应变和峰值应力均减小 ， 箍筋锈蚀对约束混凝土的峰值应力影响要大于峰值 应变 ； 2 ) 分析箍筋锈蚀对约束混凝土性能的影响， 建 立了箍筋锈蚀率与约束混凝土承载力及变形性能之 间的定量关系式； 3 ) 建立锈蚀箍筋约束混凝土应力一 应变全曲线 模型， 为建立箍筋和主筋共同锈蚀本构模型提供基 础。 ( 下转第 l 1 1页) 徐咏， 等： 通济桥倒塌原因调查与分析 分钟内， 3号 、 4号桥拱随之垮塌。 图 7后续垮塌 Fi g 7 S u b s e q u e n t c o l l a p s e 2 3垮塌的受力分析 拱桥区别于其他类型桥梁的一个特性是其具有 连拱效应， 其概念可以表述为单拱上的荷载传至拱 脚会产生竖向压力和水平推力， 连拱的中间各个拱 脚 因左右来的水平推力方 向相反 ， 大小相 同而平衡 稳定， 只剩下竖向压力传递给基础。该桥因为长期 的河水 冲刷和河砂开采 ， 基础已被一定程度的掏空。 5 号桥墩基础失去承载力后 ， 桥墩外侧条石脱落， 其 上拱圈拱脚处条石因为失去传力途径而相继连续脱 落， 致使5号桥墩左右水平力不等而失稳， 同时拱圈 处条石因为完全失去传力途径而垮塌， 其他拱圈亦 因连拱效应而相继垮塌 。 3该拱桥垮塌原因总结 由以上 的分析可知 ， 垮塌 的客观原因是垮塌当 天和前后几天 的极端 恶劣天气 以及 采砂 的持续 工 作； 力学原因是4号基础被掏空导致桥墩外侧条石 脱落， 拱圈条石失去传力途径后条石继续脱落并引 起桥墩水平受力不均而失稳 ， 最后 ， 因为拱桥的连拱 效应致使这座桥垮塌 _ 6 J 。 随着气候环境与地理条件 的不断变化， 现今世 界范 围内的很多江河出现了运载量增大洪峰过境水 量大等新情况。面对这一情况， 很多古桥面临着诸 多新 的挑战 ， 比如如今 的桥面荷载相 比于修建古 桥 时候考虑的荷载而言已经有了几何级数的增长。即 使桥梁修建 时已经很 富足地考 虑 了桥梁 的承载 能 力， 在现在看来也完全不能满足现代化交通的需要； 又如现今 的环境条件造成水质普遍变差 ， 对桥墩 和 基础提出了新 的挑 战； 再如近两百年来 的采砂工作 造成了河床松软， 基础承载力变低等。所以， 如今的 桥梁工作者在大力发展 中国现代桥梁事业 的同时 ， 万万不能忽略我国古桥的存在及其现实性难题， 从 爱护文物与保证 交通 两方面出发去维护升级古桥 梁 ， 特别应该注意桥墩与基础 的检查与维护 ， 以此桥 梁为前车之鉴 ， 争取通过努力 ， 不再让类似事件发 生 。 鸣谢： 感谢双流县气象局、 水务局、 交通运输局、 望江水文站提供珍贵资料 ； 感谢该桥 两岸周 围居 民 热情讲述实时情况并提供珍贵视频参考。 参 考 文 献 ： 1 四川省双流县志编纂委员会 双流县志 M 成都： 四川人民 出版社 ， 1 9 9 2 2 滇安 拱桥为什么会坍塌 J 劳动保护杂志， 1 9 9 9 ( 5 ) ： 3 2 3 黄力山 砌体结构裂缝的成因及控制措施 J 安徽建筑， 2 0 0 3 ( 4 )： 6 7 - 6 8 4 金吉寅， 冯郁芬 桥梁附属构造与支座 M 北京： 人民交通出 版社 ， 1 9 9 9 5 杨文渊， 徐彝 桥梁施工工程师手册 M 北京： 人民交通出 版社 ， 1 9 9 7 6 吉伯海， 傅中秋 近年国内桥梁倒塌事故原因分析 J 土木工 程学报 ， 2 0 1 0 ( 4 3 ) ： 4 9 6 - 4 9 8 ( 上接第9 7页) 参 考 文 献 ： 1 C a i ms J ， D u Y， L a w D I n fl u e n c e o f c o r r o s i o n o n t h e f r i c t i o n c h a r - a c t e r i s t i c s o f t h e s t e e l c o n c r e t e i n t e r f a c e J C o n s t r u c t i o n a n d B u i l d i n g MA T E RI A L S， 2 0 0 7( 2 1 ) ： 1 9 0 1 9 7 2 C h r i s t o p h e r H i g g i n s ， w 1 l i a m F a r r o w T e a t s o f re i n f o r c e d c o n c r e t e b e a m s w i t h c o r r o s i o n d a m a g e d s t i r r u p s J J A C I S t r u c t u ra l J o u ma l ， 2 0 0 6， J a n u a r y F e b rua r y： 1 3 3 - 1 41 3 徐善华 ， 牛荻涛 锈蚀钢筋混凝土简支梁斜截面抗剪性能研究 J 建筑结构学报， 2 0 0 4 ( 5 ) 4 赵羽习， 金伟良 锈蚀箍筋混凝土梁的抗剪承载力分析 J 】 浙 江大学 学报 ， 2 0 0 8 ， 4 2 ( 1 ) ： 1 9 - 2 4 5 刘西拉 重大土木与水利工程安全性及耐久性的基础研究 J 土木工程学报， 2 0 0 1 ， 3 4 ( 6 ) ： 1 -4 6 王金昌， 陈页开 A B A Q U S 在土木工程中的应用 M 杭州 ： 浙 江大学 出版社 。 2 o o 6 7 过镇海， 时旭东 钢筋混凝土原理与分析 M 北京： 清华大学 出版社 ， 2 0 0 3 8 L e e H S ， N o g u c h i T ， T o m o s a w a F F E M A n a l y s i s f o r S t r u c t u r e P e r f o r ma n e e o f De t e r i o r a t e d RC S t r u c t u r e d u e t o Re b a r C o r r o s i o n c P r o c e e d i n g of t h e I n te rn a t i o n a l C o n c rete u n d e r S e v e r a l Co n d i t i o n s Tr o ms a， 1 9 9 8 9 R o b in s o n J 有限单元法的综合理论 M 邹焕新， 等译 北京： 人民教育出版社， 1 9 8 2 1 0 牛荻涛 混凝土结构耐久性与寿命预测 M j E 京： 科学出版 社 2 o o 3