《工程科学学报》录用稿，htps:/doi.org/10.13374/i,issn2095-9389.2021.06.29.008©北京科技大学2020 工程科学学报 DOI: 111公式章1节1纵筋锈蚀无腹筋混凝土梁抗剪 性能细观数值研究 金浏”，杨鸿申”，张仁波2)区，杜修力) 1)北京工业大学城市与工程安全减灾教育部重点实验室，北京1001242)清华大学土木系，北京 ☒通信作者，E-mail:zhangrenbo99@126.com 摘要钢筋锈蚀是影响钢筋混凝土结构耐久性的主要原因。锈蚀减小钢筋藏面积和性能，钢筋对混凝土的销栓作用 减弱。钢筋锈蚀严重时，混凝土保护层开裂，钢筋与混凝土之间的结性能降低。本文以钢筋混凝土梁为研究对象 考虑钢筋非均匀锈蚀膨胀效应，建立三维钢筋混凝土梁剪切破坏分析的数值分析模型。锈蚀对结构力学行为的影响 通过多阶段分析方法探索（钢筋锈蚀阶段，构件性能退化阶段 钢筋的非均匀锈蚀膨胀以施加非均匀径向位移的方 式模拟，获得保护层的破坏状态，并以此“最终状态”作为后泥凝土梁静载试验的“初始条件”输入，进而模拟 构件的力学行为。在验证了多阶段数值模型合理性的基础上，犷析了纵筋锈蚀、剪跨比对无腹筋混凝土梁抗剪性能 的影响规律。结果表明：纵筋锈蚀使混凝土梁产生明显的纵向裂缝。纵筋锈蚀率增大，保护层开裂区域增加，梁的 抗剪承载力下降严重。另外，剪跨比对梁的抗剪承载力产生影响，剪跨比对未锈蚀梁的影响明显大于对锈蚀梁的影 响程度。最后，基于模拟结果对相关设计规范中的抗剪承载力计算公式进行了讨论，发展建立了考虑锈蚀影响的无 腹筋混凝土梁抗剪承载力计算方法。 关键词钢筋混凝土梁：纵筋：锈蚀八细观数值模拟：抗剪行为 分类号TU375.1 Study on the corrosion effects of longitudinal reinforcement on shear behavior of concrete beams without web reinforcement JIN Liu,YANG Hong-shen,ZHANG Ren-bo2 DU Xiu-li) 1)The Key Laboratory of Urban Security and Disaster Engineering,Ministry of Education,Beijing University of Technology,Beijing 100124, China 2)Department of Civil Engineering,Tsinghua University,Beijing,100084,China Corresponding author,E-mail:E-mail:zhangrenbo99@126.com ABSTRACT Rebar corrosion is the main factor affecting the service performance of reinforced concrete (RC)structures On the one hand,corrosion reduces the area and the performance of reinforcements,weakens the pin bolt effect of rebar on 收精日期2021-06-Xx; 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： 111 公式章 1 节 1 纵筋锈蚀无腹筋混凝土梁抗剪 性能细观数值研究 金 浏 1)，杨鸿申 1)，张仁波 1,2) ，杜修力 1) 1) 北京工业大学 城市与工程安全减灾教育部重点实验室, 北京 100124 2) 清华大学 土木系, 北京 100084 通信作者，E-mail: zhangrenbo99@126.com 摘 要 钢筋锈蚀是影响钢筋混凝土结构耐久性的主要原因。锈蚀减小钢筋截面积和性能，钢筋对混凝土的销栓作用 减弱。钢筋锈蚀严重时，混凝土保护层开裂，钢筋与混凝土之间的粘结性能降低。本文以钢筋混凝土梁为研究对象 考虑钢筋非均匀锈蚀膨胀效应，建立三维钢筋混凝土梁剪切破坏分析的数值分析模型。锈蚀对结构力学行为的影响 通过多阶段分析方法探索(钢筋锈蚀阶段，构件性能退化阶段)。钢筋的非均匀锈蚀膨胀以施加非均匀径向位移的方 式模拟，获得保护层的破坏状态，并以此“最终状态”作为之后混凝土梁静载试验的“初始条件”输入，进而模拟 构件的力学行为。在验证了多阶段数值模型合理性的基础上，分析了纵筋锈蚀、剪跨比对无腹筋混凝土梁抗剪性能 的影响规律。结果表明：纵筋锈蚀使混凝土梁产生明显的纵向裂缝。纵筋锈蚀率增大，保护层开裂区域增加，梁的 抗剪承载力下降严重。另外，剪跨比对梁的抗剪承载力产生影响，剪跨比对未锈蚀梁的影响明显大于对锈蚀梁的影 响程度。最后，基于模拟结果对相关设计规范中的抗剪承载力计算公式进行了讨论，发展建立了考虑锈蚀影响的无 腹筋混凝土梁抗剪承载力计算方法。 关键词 钢筋混凝土梁；纵筋；锈蚀；细观数值模拟；抗剪行为 分类号 TU375.1 Study on the corrosion effects of longitudinal reinforcement on shear behavior of concrete beams without web reinforcement JIN Liu1), YANG Hong-shen1), ZHANG Ren-bo1,2) , DU Xiu-li1) 1) The Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education, Beijing University of Technology, Beijing 100124, China 2) Department of Civil Engineering, Tsinghua University, Beijing, 100084, China Corresponding author, E-mail: E-mail: zhangrenbo99@126.com ABSTRACT Rebar corrosion is the main factor affecting the service performance of reinforced concrete (RC) structures. On the one hand, corrosion reduces the area and the performance of reinforcements, weakens the pin bolt effect of rebar on 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 《工程科学学报》录用稿，https://doi.org/10.13374/j.issn2095-9389.2021.06.29.008 ©北京科技大学 2020 录用稿件，非最终出版稿
工程科学学报 DOI: concrete.On the other hand,when the rebar is seriously rusted,the concrete cover cracks,the bond behavior between reinforcement and concrete is deteriorated,which all affect the mechanical properties of the RC structure.In this paper. taking RC beam as the research object,a three-dimensional numerical model for shear analysis was established involving the non-uniform corrosion of reinforcement.The corrosion effects on the mechanical behavior of RC beam was explored through a multi-stage analysis method (namely,corrosion induced expansion stage and structural deterioration stage).The non- uniform radial displacement was applied to concrete around the rebar to simulate the expansion of the corrosion products. The cracking process and the damage pattern of concrete caused by corrosion can be obtained.Then,taking the corrosion state as the initial condition,the static load was applied to analysis the mechanical behavior of the RC beam.After verifying the rationality of the multi-stage numerical model,the influence of corrosion of tensile reinforcement and the shear-span ratio on the shear behavior of concrete beams without web reinforcement was analyzed.The modeling analysis results show that corrosion of longitudinal reinforcement causes obvious longitudinal corrosion cracks in concrete beams.And with the development of corrosion,the cracking area of protective layer increases,the shear capacitydecreases seriously. Besides,the shear span ratio has an impact on the shear capacity of beams,and the impaet on non-corroded beams is significantly greater than that of corroded beams.Finally,based on the simulation result eulation formulas of shear capacity in relevant design codes were discussed,and a methodology for the prediction o hear capacity of RC beams without web reinforcement was proposed. KEY WORDS RC beam;longitudinal reinforcement;Corrosion;Meso-scale numeri ulation;Shear behavior 钢筋混凝土结构在服役过程中常因遭受碳化、有害物质侵蚀等不良因素作用而造成内部钢筋锈 蚀山。钢筋锈蚀后自身截面积减小，钢筋的力学性能降低。 锈蚀产物体积膨胀对周围混凝土造成挤 压，导致混凝土保护层开裂、剥落，混凝土的性能及钢筋写混凝土之间的粘结强度降低。锈蚀严 重影响了混凝土结构的安全性和耐久性，需要花费大量资金用于锈蚀混凝土结构的维修与改造)。 因此，研究锈蚀对钢筋混凝土结构力学性能的影彩响对于锈蚀后混凝土结构的性能评估和加固设计具 有重要的理论价值与实际意义。 钢筋混凝土梁是混凝土结构中的重要受力构件，国内外学者对于其内部钢筋锈蚀引起的构件性 能退化开展了大量的研究工作，主要包括试验研究、理论分析和数值模拟三方面。Alaskar等进行 了9根不同锈蚀率的混凝土梁剪切试验，发现箍筋锈蚀导致梁表面产生了竖向裂缝。在箍筋锈蚀率 超过10%时，梁的承载力下降严重。Sfem等开展了剪跨比和箍筋锈蚀对梁承载力的影响试验， 结果显示箍筋锈蚀率相同时锈烛对较小剪跨比的梁的承载力影响更为明显。黄天立、戴明江等) 对仅纵筋锈蚀的混凝土梁进剪场试验，发现随着纵筋锈蚀率的增加，锈胀裂缝宽度逐渐增大，梁 的承载能力降低。X©等的试验结果表明纵筋锈蚀对小剪跨比梁(=1.5)的承载力影响不明显， 但对剪跨比为2、4的混凝土梁承载力产生了较大影响。 理论研究方面，徐善华等例根据极限平衡理论，建立了考虑箍筋锈蚀影响的混凝土简支梁斜截 面承载力计算模型。Khan等1通过截面分析方法和拉压杆模型对锈蚀混凝土深梁的抗剪承载力进 行了分析讨论，发视截面法计算结果较试验相对保守，拉压杆模型能够更准确计算锈蚀深梁的承载 力。EI-Sayed四采用锈蚀裂缝宽度来评估锈蚀损伤，提出了适用于锈蚀钢筋混凝土梁的剪切承载力 计算方法，并且与已有试验数据进行了验证。Azam等2,1分析了纵筋锈蚀前后混凝土梁受力机制 的变化，认为纵筋锈蚀使梁的荷载传递机制由梁作用转变为拱作用，纵筋锈蚀使梁的承载力提高。 纵筋锈蚀不仅导致其力学性能下降，锈蚀严重时，还会造成混凝土保护层产生纵筋裂缝，严重破坏 钢筋与混凝土之间的粘结行为，梁的承载机制和抗剪能力发生变化。然而，目前针对纵筋锈蚀对 混凝土梁剪切性能影响的认识仍然不够深入，需要开展进一步研究工作。 收离日期2021-06-xx 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： concrete. On the other hand, when the rebar is seriously rusted, the concrete cover cracks, the bond behavior between reinforcement and concrete is deteriorated, which all affect the mechanical properties of the RC structure. In this paper, taking RC beam as the research object, a three-dimensional numerical model for shear analysis was established involving the non-uniform corrosion of reinforcement. The corrosion effects on the mechanical behavior of RC beam was explored through a multi-stage analysis method (namely, corrosion induced expansion stage and structural deterioration stage). The nonuniform radial displacement was applied to concrete around the rebar to simulate the expansion of the corrosion products. The cracking process and the damage pattern of concrete caused by corrosion can be obtained. Then, taking the corrosion state as the initial condition, the static load was applied to analysis the mechanical behavior of the RC beam. After verifying the rationality of the multi-stage numerical model, the influence of corrosion of tensile reinforcement and the shear-span ratio on the shear behavior of concrete beams without web reinforcement was analyzed. The modeling analysis results show that corrosion of longitudinal reinforcement causes obvious longitudinal corrosion cracks in concrete beams. And with the development of corrosion, the cracking area of protective layer increases, the shear capacity of beam decreases seriously. Besides, the shear span ratio has an impact on the shear capacity of beams, and the impact on non-corroded beams is significantly greater than that of corroded beams. Finally, based on the simulation results, the calculation formulas of shear capacity in relevant design codes were discussed, and a methodology for the prediction of the shear capacity of RC beams without web reinforcement was proposed. KEY WORDS RC beam; longitudinal reinforcement; Corrosion; Meso-scale numerical simulation; Shear behavior 钢筋混凝土结构在服役过程中常因遭受碳化、有害物质侵蚀等不良因素作用而造成内部钢筋锈 蚀[1]。钢筋锈蚀后自身截面积减小，钢筋的力学性能降低。锈蚀产物体积膨胀对周围混凝土造成挤 压，导致混凝土保护层开裂、剥落，混凝土的性能及钢筋与混凝土之间的粘结强度降低[2]。锈蚀严 重影响了混凝土结构的安全性和耐久性，需要花费大量资金用于锈蚀混凝土结构的维修与改造[3]。 因此，研究锈蚀对钢筋混凝土结构力学性能的影响对于锈蚀后混凝土结构的性能评估和加固设计具 有重要的理论价值与实际意义。 钢筋混凝土梁是混凝土结构中的重要受力构件，国内外学者对于其内部钢筋锈蚀引起的构件性 能退化开展了大量的研究工作，主要包括试验研究、理论分析和数值模拟三方面。Alaskar 等[4]进行 了 9 根不同锈蚀率的混凝土梁剪切试验，发现箍筋锈蚀导致梁表面产生了竖向裂缝。在箍筋锈蚀率 超过 10%时，梁的承载力下降严重。Suffern 等[5]开展了剪跨比和箍筋锈蚀对梁承载力的影响试验， 结果显示箍筋锈蚀率相同时，锈蚀对较小剪跨比的梁的承载力影响更为明显。黄天立、戴明江等[6,7] 对仅纵筋锈蚀的混凝土梁进行剪切试验，发现随着纵筋锈蚀率的增加，锈胀裂缝宽度逐渐增大，梁 的承载能力降低。Xue 等[8]的试验结果表明纵筋锈蚀对小剪跨比梁（λ=1.5）的承载力影响不明显， 但对剪跨比为 2、4 的混凝土梁承载力产生了较大影响。 理论研究方面，徐善华等[9]根据极限平衡理论，建立了考虑箍筋锈蚀影响的混凝土简支梁斜截 面承载力计算模型。Khan 等[10]通过截面分析方法和拉压杆模型对锈蚀混凝土深梁的抗剪承载力进 行了分析讨论，发现截面法计算结果较试验相对保守，拉压杆模型能够更准确计算锈蚀深梁的承载 力。EI-Sayed[11]采用锈蚀裂缝宽度来评估锈蚀损伤，提出了适用于锈蚀钢筋混凝土梁的剪切承载力 计算方法，并且与已有试验数据进行了验证。Azam 等[12,13]分析了纵筋锈蚀前后混凝土梁受力机制 的变化，认为纵筋锈蚀使梁的荷载传递机制由梁作用转变为拱作用，纵筋锈蚀使梁的承载力提高。 纵筋锈蚀不仅导致其力学性能下降，锈蚀严重时，还会造成混凝土保护层产生纵筋裂缝，严重破坏 钢筋与混凝土之间的粘结行为，梁的承载机制和抗剪能力发生变化[14]。然而，目前针对纵筋锈蚀对 混凝土梁剪切性能影响的认识仍然不够深入，需要开展进一步研究工作。 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 录用稿件，非最终出版稿
工程科学学报 DOI: 在试验研究和理论分析之外，数值模拟也是常用的研究方法之一。目前有关锈蚀数值模拟的工 作主要集中于混凝土锈胀开裂行为以及锈蚀混凝土构件的力学性能两个方面。在混凝土保护层锈胀 开裂模拟过程中，通过对钢筋周围混凝土施加温度或位移荷载6来模拟钢筋锈蚀膨胀对混凝土所 造成的挤压。在锈蚀混凝土构件力学性能分析中则通过等效化折减截面积或力学强度等方式来考虑 钢筋锈蚀导致的钢筋截面积减小、混凝土和钢筋力学性能及二者之间相互作用的劣化等变化，1。 尽管上述工作可以较好地反映混凝土锈胀开裂行为以及锈蚀混凝土构件的整体力学性能，然而现有 数值模拟工作中一般将这两个过程分别进行考虑，并在构件整体性能分析中进行等效化处理，这与 锈蚀钢筋混凝土构件的力学行为演化过程有所差异。 鉴于此，为了更好地体现钢筋锈蚀膨胀效应，本文将锈蚀引起的保护层开裂及构件性能劣化过 程相结合，建立了锈蚀混凝土梁多阶段数值分析模型，研究纵筋锈蚀对混凝土梁剪切行为的影响规 律。在与己有试验结果吻合良好的基础上，模拟分析不同纵筋锈蚀率下混凝土保护层和混凝土梁的 剪切破坏模式与失效机制，并探讨了剪跨比对梁剪切性能的影响规律。最后，通过理论值与模拟值 的对比，对现有的承载力计算公式进行了修正。 1钢筋混凝土梁细观分析模型 1.1数值模拟方法 自然环境中，钢筋锈蚀引起混凝土结构性能的变化是一个随时间丽发展的过程。环境中的有害 物质侵入混凝土保护层之后破坏钢筋表面的钝化膜，钢筋开始锈蚀钢筋锈蚀后生成锈蚀产物的体 积是消耗钢筋体积的2~6倍，锈蚀产物体积膨胀挤压周围混凝土，导致混凝土开裂、分层，甚至 剥落。 考虑混凝土保护层的锈胀破坏在研究锈蚀对混凝：结构力学性能的影响过程中是必不可少的。 现有研究对于锈蚀导致的混凝土保护层开裂以及构件性能恶化多进行单独分别考虑[6，，较少有将 两者相联系起来。为了更加接近于实际情况，文对二者进行综合考虑，通过一个模型的多步分析 来体现，建立了多阶段锈蚀数值模型，如图1所示： Stage 1:Rebar Stage 2:Mechanical corrosion,expansion behavior of RC beams 录用稿 Concrete cracks 圆1多阶段锈蚀数值模拟方法 Framework of the numerical approach:Multi-stage analysis strategy 阶段1：混凝土保护层锈胀开裂行为分析。对钢筋周围的混凝土进行强制位移加载，以此表示 锈蚀产物膨胀对钢筋周围混凝土的力学作用。通过钢筋锈蚀膨胀模拟，可以再现混凝土保护层锈胀 开裂的破坏过程。同时能够获得较为真实的混凝土破坏状态，例如混凝土的开裂模式、锈胀压力等。 阶段2：锈蚀构件力学行为模拟。考虑钢筋锈蚀导致的混凝土保护层开裂，在阶段1钢筋锈胀 模拟结果的基础上，进行钢筋混凝士梁三点加载的数值试验。分析纵筋锈蚀对钢筋混凝士梁剪切行 为的影响规律。 1.2细观数值横型立 为了分析纵筋锈蚀无腹筋混凝土梁的力学行为，本文以Han等2o试验为参照，建立如图2所示 收离日期2021-06-xx 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： 在试验研究和理论分析之外，数值模拟也是常用的研究方法之一。目前有关锈蚀数值模拟的工 作主要集中于混凝土锈胀开裂行为以及锈蚀混凝土构件的力学性能两个方面。在混凝土保护层锈胀 开裂模拟过程中，通过对钢筋周围混凝土施加温度[15]或位移荷载[16]来模拟钢筋锈蚀膨胀对混凝土所 造成的挤压。在锈蚀混凝土构件力学性能分析中则通过等效化折减截面积或力学强度等方式来考虑 钢筋锈蚀导致的钢筋截面积减小、混凝土和钢筋力学性能及二者之间相互作用的劣化等变化[17,18]。 尽管上述工作可以较好地反映混凝土锈胀开裂行为以及锈蚀混凝土构件的整体力学性能，然而现有 数值模拟工作中一般将这两个过程分别进行考虑，并在构件整体性能分析中进行等效化处理，这与 锈蚀钢筋混凝土构件的力学行为演化过程有所差异。 鉴于此，为了更好地体现钢筋锈蚀膨胀效应，本文将锈蚀引起的保护层开裂及构件性能劣化过 程相结合，建立了锈蚀混凝土梁多阶段数值分析模型，研究纵筋锈蚀对混凝土梁剪切行为的影响规 律。在与已有试验结果吻合良好的基础上，模拟分析不同纵筋锈蚀率下混凝土保护层和混凝土梁的 剪切破坏模式与失效机制，并探讨了剪跨比对梁剪切性能的影响规律。最后，通过理论值与模拟值 的对比，对现有的承载力计算公式进行了修正。 1 钢筋混凝土梁细观分析模型 1.1 数值模拟方法 自然环境中，钢筋锈蚀引起混凝土结构性能的变化是一个随时间而发展的过程。环境中的有害 物质侵入混凝土保护层之后破坏钢筋表面的钝化膜，钢筋开始锈蚀。钢筋锈蚀后生成锈蚀产物的体 积是消耗钢筋体积的 2~6 倍，锈蚀产物体积膨胀挤压周围混凝土[19]，导致混凝土开裂、分层，甚至 剥落。 考虑混凝土保护层的锈胀破坏在研究锈蚀对混凝土结构力学性能的影响过程中是必不可少的。 现有研究对于锈蚀导致的混凝土保护层开裂以及构件性能恶化多进行单独分别考虑[16,17]，较少有将 两者相联系起来。为了更加接近于实际情况，本文对二者进行综合考虑，通过一个模型的多步分析 来体现，建立了多阶段锈蚀数值模型，如图 1 所示： Stage 1: Rebar corrosion, expansion Forced diaplacement Concrete cracks Stage 2: Mechanical behavior of RC beams 图 1 多阶段锈蚀数值模拟方法 Fig.1 Framework of the numerical approach: Multi-stage analysis strategy 阶段 1：混凝土保护层锈胀开裂行为分析。对钢筋周围的混凝土进行强制位移加载，以此表示 锈蚀产物膨胀对钢筋周围混凝土的力学作用。通过钢筋锈蚀膨胀模拟，可以再现混凝土保护层锈胀 开裂的破坏过程。同时能够获得较为真实的混凝土破坏状态，例如混凝土的开裂模式、锈胀压力等。 阶段 2：锈蚀构件力学行为模拟。考虑钢筋锈蚀导致的混凝土保护层开裂，在阶段 1 钢筋锈胀 模拟结果的基础上，进行钢筋混凝土梁三点加载的数值试验。分析纵筋锈蚀对钢筋混凝土梁剪切行 为的影响规律。 1.2 细观数值模型建立 为了分析纵筋锈蚀无腹筋混凝土梁的力学行为，本文以 Han 等[20]试验为参照，建立如图 2 所示 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 录用稿件，非最终出版稿
工程科学学报 DOI: 的三维细观数值模型。梁的截面为矩形（170mm×250mm),梁长1400mm,保护层厚度和截面有 效高度分别为30mm、203mm,剪跨比为2.96。纵向钢筋采用直径22mm的带肋钢筋，屈服强度 为400MPa。 Loading point Beam Base Mortar☐Aggregate ITZ 圆2钢筋混凝土梁三维细观模型 Fig.2 Meso-scale numerical model of the RC beam. 版稿 node Concrete Steel node Steel element Spring ☑3弹簧单元 Fig.3 Sketch of the spring element 在细观数值模型中，将混凝土看作由砂浆、骨料及两者之间的界面过渡区组成的三相复合材料。 采用八节点实体单元划分网格后，通过Fortran自编程序完成骨料，砂浆以及界面各相单元的区分， 各相交界面处的单元共用节点)人保迹位移和应力连续。其中，界面过渡区厚度设定为2mm,骨料 形状假定为圆形，体积分数为30%。混凝土为二级配，等效粒径分别为15mm和8mm1。纵向钢 筋同样采用八节点实体单元进行网格划分。与文献22]类似，钢筋与混凝士相邻单元节点之间通过 三个相互垂直的非线性弹簧单元进行连接。如图3所示，弹簧单元系统由三个相互垂直的弹簧组成， 平行钢筋长度方向的弹簧模拟钢筋与混凝土的粘结滑移，垂直钢筋方向的弹簧模拟钢筋与混凝土的 径向作用。 1.3本构横型和力学参数 已有研究表，混凝土中骨料强度较高，在静力作用下很难破坏，因此假定骨料为弹性体，即 认为骨料只产生弹性变形。混凝土中的砂浆和界面过渡区采用塑性损伤本构关系模型。该模型中 假定材料只发生拉伸和压缩破坏，并且拉伸和压缩引起的材料刚度退化通过损伤变量来描述。材料 在单轴拉伸和压缩情况下的应力-应变关具体表达式为： o,=(1-d)E(e,-e) 22\*MERGEFORMAT() e=(1-de)Eo(se-s) 33\*MERGEFORMAT O 收离日期2021-06-xx 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： 的三维细观数值模型。梁的截面为矩形（170 mm×250 mm）,梁长 1400 mm，保护层厚度和截面有 效高度分别为 30 mm、203 mm，剪跨比为 2.96。纵向钢筋采用直径 22 mm 的带肋钢筋，屈服强度 为 400 MPa。 Mortar Aggregate ITZ Beam Loading point Base 图 2 钢筋混凝土梁三维细观模型 Fig.2 Meso-scale numerical model of the RC beam. y x z1 x1 y1 z2 x2 y2 Concrete element Steel element Concrete node Steel node z z y x Spring Spring Spring 图 3 弹簧单元 Fig.3 Sketch of the spring element 在细观数值模型中，将混凝土看作由砂浆、骨料及两者之间的界面过渡区组成的三相复合材料。 采用八节点实体单元划分网格后，通过 Fortran 自编程序完成骨料，砂浆以及界面各相单元的区分， 各相交界面处的单元共用节点，保证位移和应力连续。其中，界面过渡区厚度设定为 2 mm，骨料 形状假定为圆形，体积分数为 30%。混凝土为二级配，等效粒径分别为 15 mm 和 8 mm [21]。纵向钢 筋同样采用八节点实体单元进行网格划分。与文献[22]类似，钢筋与混凝土相邻单元节点之间通过 三个相互垂直的非线性弹簧单元进行连接。如图 3 所示，弹簧单元系统由三个相互垂直的弹簧组成， 平行钢筋长度方向的弹簧模拟钢筋与混凝土的粘结滑移，垂直钢筋方向的弹簧模拟钢筋与混凝土的 径向作用。 1.3 本构模型和力学参数 已有研究表明，混凝土中骨料强度较高，在静力作用下很难破坏，因此假定骨料为弹性体，即 认为骨料只产生弹性变形。混凝土中的砂浆和界面过渡区采用塑性损伤本构关系模型[23]。该模型中 假定材料只发生拉伸和压缩破坏，并且拉伸和压缩引起的材料刚度退化通过损伤变量来描述 。材料 在单轴拉伸和压缩情况下的应力-应变关具体表达式为： pl t t 0 t t (1 ) ( ) d E 22\* MERGEFORMAT () pl c c 0 c c (1 ) ( ) d E 33\* MERGEFORMAT () 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 录用稿件，非最终出版稿
工程科学学报 DOI: 式中，E为材料初始弹性模量，4，、：分别为拉伸损伤变量和压缩损伤变量，表示拉伸和压缩引 起的材料弹性刚度的退化，其变化范围为0~1（分别对应无损伤至完全破坏），&、分别为单 轴拉伸、压缩条件下的等效塑性应变。关于塑性损伤模型的详细描述，可参考文献。 混凝土细观成分的力学参数如表1所示。其中，骨料和砂浆的力学参数为试验实测值，界面过 渡区的力学参数按照砂浆强度的70%~85%进行反复试算，最终选取最接近试验实测抗压强度的一 组数据作为界面过渡区的参数。采用表1细观力学参数所测得棱柱体的抗压强度为39MP,与试验 所测棱柱体抗压强度(39.2MPa)较为接近，说明了所选参数的合理性。 表1混凝土细观组分力学参数 Table 1 Mechanical parameters for meso-components of concrete Aggregate Mortar ITZ Compressive strength oe/MPa *39 31 Tensile strength /MPa 3. 2.8 Elastic modulus E/GPa 70 最终版裤 26 Poisson ratio v 0.2 0.2 注：*代表试验实测值m,为反复试算取值。 ■4本构模型 Fig.4 Constitutive models used in the numerical analysis 钢筋本构关系采用双拆线性-强化模型，具体参数按照文献取值。钢筋与混凝士节点间的 三个相互垂直的弹簧设置为不同刚度。由于钢筋与混凝土在钢筋法向方向上难以互相侵入，并且一 般不会产生相对位移因此通常假定钢筋与混凝土在垂直钢筋长度的方向上刚结，法向刚度设置为 较大值。平行于钢筋长度方向弹簧的刚度按照《混凝士结构设计规范》规定的钢筋混凝土之间的 粘结滑移关系 图4) 采用式(3)计算获得，图4中具体参数的含义及取值详见表2。 F=t(s).A 44\*MERGEFORMAT() 其中，t(⊙)为钢筋与混凝土之间的粘结滑移关系，A为弹簧连接节点所在单元的面积。 囊2粘结滑移本构参数 Table 2 Parameters utilized in the bond-slip model of steel and concrete Key point Crack(cr) Peak value(u) Remnant(r） Stress t/MPa to=2.5h t.=3.0f 五=f Slip s/mm S=0.025d 3=0.04d S=0.55d 收离日期2021-06-xx 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： 式中， E0 为材料初始弹性模量， t d 、 c d 分别为拉伸损伤变量和压缩损伤变量，表示拉伸和压缩引 起的材料弹性刚度的退化，其变化范围为 0 ~ 1（分别对应无损伤至完全破坏）， pl t 、 pl c 分别为单 轴拉伸、压缩条件下的等效塑性应变。关于塑性损伤模型的详细描述，可参考文献[23]。 混凝土细观成分的力学参数如表 1 所示。其中，骨料和砂浆的力学参数为试验实测值，界面过 渡区的力学参数按照砂浆强度的 70%~85%进行反复试算，最终选取最接近试验实测抗压强度的一 组数据作为界面过渡区的参数。采用表 1 细观力学参数所测得棱柱体的抗压强度为 39MPa，与试验 所测棱柱体抗压强度（39.2MPa）较为接近，说明了所选参数的合理性。 表 1 混凝土细观组分力学参数 Table 1 Mechanical parameters for meso-components of concrete Aggregate Mortar ITZ Compressive strength σc /MPa - *39 ^31 Tensile strength σt /MPa - *3.5 ^2.8 Elastic modulus E /GPa 70 32* ^26 Poisson ratio ν 0.2 0.2 0.2 注：*代表试验实测值[20],^为反复试算取值。 0 cr u r scr su sr τ τu τcr τr s Slip Stress 图 4 本构模型 Fig.4 Constitutive models used in the numerical analysis 钢筋本构关系采用双折线弹性-强化模型，具体参数按照文献[20]取值。钢筋与混凝土节点间的 三个相互垂直的弹簧设置为不同刚度。由于钢筋与混凝土在钢筋法向方向上难以互相侵入，并且一 般不会产生相对位移，因此通常假定钢筋与混凝土在垂直钢筋长度的方向上刚结，法向刚度设置为 较大值。平行于钢筋长度方向弹簧的刚度按照《混凝土结构设计规范》[24]规定的钢筋混凝土之间的 粘结滑移关系（图 4）采用式(3)计算获得，图 4 中具体参数的含义及取值详见表 2。 F s A ( ) 44\* MERGEFORMAT () 其中， ( )s 为钢筋与混凝土之间的粘结滑移关系，A 为弹簧连接节点所在单元的面积。 表 2 粘结-滑移本构参数 Table 2 Parameters utilized in the bond-slip model of steel and concrete. Key point Crack(cr) Peak value(u) Remnant(r) Stress τ/MPa τcr= 2.5ft τu = 3.0 ft τr = ft Slip s/mm Scr=0.025d Su = 0.04d Sr = 0.55d 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 录用稿件，非最终出版稿
工程科学学报 DOI: 1.4师做效应 1.4.1锈层分布 自然环境中，特别是氯盐侵蚀环境下，钢筋截面的锈蚀往往是不均匀的。往往靠近保护层一侧 钢筋锈蚀严重，远离保护层一侧钢筋轻微锈蚀甚至不锈。Yua等2发现钢筋截面内的锈层分布近 似成半椭圆型。考虑到中间位置钢筋仅受到一个方向有害物质的侵蚀，而角部钢筋受到两个方向有 害物质的侵蚀，因此角部钢筋的锈层分布采用将两个中部钢筋锈层分布叠加的方式获得，如图5所示。 角部钢筋锈层分布用极坐标方程表示如式(4)。 钢筋锈蚀膨胀导致混凝土保护层开裂通常包含三个阶段：锈蚀产物自由膨胀阶段、混凝土保护 层受拉阶段和保护层开裂阶段。在锈蚀产物半椭圆型分布的基础上，薛圣广分析单位长度内钢 筋锈蚀体积与钢筋初始体积的关系，推导出有关钢筋锈蚀率与截面最大锈胀位移的关系式，如式 (5). 元2 π32 版稿 Middle rebar Middle rebar Comer reba 圆5钢筋锈层分布图 Fig.5 Rust distribution of the steel bar R+)-点久 4(0≤0<)< -R(≤0) (R+u)cos20+(R+u sin20 ue (R+4)(R+42) R( L≤0<2π) (+4)2sinθ+(R+4,)2cos20 551* MERGEFORMAT ( 稿 = 280+4+42 (P-1)R 66\*MERGEFORMAT() 其中，R为钢筋的初始竿径： 。为对应0角的混凝土锈层位移：4和分别对应钢筋截面内的最大、 最小锈胀位移2②0~304网这里取4=30：？为钢筋的锈蚀率；心为钢筋与混凝土之间空隙 过渡区的厚度，其值为10~20m2网，本文采用Liu等的建议，取6=12.5m:p为锈蚀产物体积 膨胀系数，p=2~69,本文中p=2,与文献1s1一致。根据式(5)及各参数之间的关系，最大锈胀位移 与锈蚀率之间的关系如式(6)所示。 30R(p-1n-26) 41=一 31 77八*MERGEFORMAT(O 收精日期：2021-06-xx; 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： 1.4 锈蚀效应 1.4.1 锈层分布 自然环境中，特别是氯盐侵蚀环境下，钢筋截面的锈蚀往往是不均匀的。往往靠近保护层一侧 钢筋锈蚀严重，远离保护层一侧钢筋轻微锈蚀甚至不锈。Yuan 等[25]发现钢筋截面内的锈层分布近 似成半椭圆型。考虑到中间位置钢筋仅受到一个方向有害物质的侵蚀，而角部钢筋受到两个方向有 害物质的侵蚀，因此角部钢筋的锈层分布采用将两个中部钢筋锈层分布叠加的方式获得,如图 5 所示。 角部钢筋锈层分布用极坐标方程表示如式(4)。 钢筋锈蚀膨胀导致混凝土保护层开裂通常包含三个阶段：锈蚀产物自由膨胀阶段、混凝土保护 层受拉阶段和保护层开裂阶段[26]。在锈蚀产物半椭圆型分布的基础上，薛圣广[27]分析单位长度内钢 筋锈蚀体积与钢筋初始体积的关系，推导出有关钢筋锈蚀率与截面最大锈胀位移的关系式，如式 (5)。 Middle rebar Corner rebar R θ u1 u2 0 π/2 π π3/2 Middle rebar 图 5 钢筋锈层分布图 Fig.5 Rust distribution of the steel bar 1 1 2 2 2 2 2 1 2 2 1 2 2 2 2 2 1 2 π (0 ) 2 ( ) ( ) π ( θ<π) ( ) cos ( ) sin 2 3π ( ) 2 ( ) ( ) 3π ( θ<2π) ( ) sin ( ) cos 2 u R u R u R R u R u u u R u R u R R u R u 55\* MERGEFORMAT () 2 0 1 2 ( 1) u u R 66\* MERGEFORMAT () 其中，R 为钢筋的初始半径；u 为对应 θ 角的混凝土锈层位移； 1 u 和 2 u 分别对应钢筋截面内的最大、 最小锈胀位移, 1 2 u u (20 ~ 30) [25] , 这里取 1 2 u u 30 ；η 为钢筋的锈蚀率； 0 为钢筋与混凝土之间空隙 过渡区的厚度，其值为 10 ~ 20 μm[28] ,本文采用 Liu 等[29]的建议，取 0 =12.5 μm；ρ 为锈蚀产物体积 膨胀系数，ρ = 2~6[19] , 本文中 ρ = 2，与文献[18]一致。根据式(5)及各参数之间的关系，最大锈胀位移 与锈蚀率之间的关系如式(6)所示。 0 1 30 1 2 31 R u 77\* MERGEFORMAT () 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 录用稿件，非最终出版稿
工程科学学报 DOI: 1.4.2锈蚀影响 锈蚀不仅导致了钢筋截面积的降低，同时造成了钢筋性能的降低。Sun等0对锈蚀钢筋进行单 轴拉伸试验，发现锈蚀使钢筋的屈服与极限强度均有所下降，且强度的降低与锈蚀率增加近似呈线 性关系。钢筋的性能通常用截面积、弹性模量、屈服和极限强度等参数来描述。本文模拟分析中， 暂不考虑锈蚀对钢筋弹性模量的影响，锈蚀后钢筋的屈服强度参考文献31]取值，如式（⑦）。 fe=(0.986-1.038n)f0 88\*MERGEFORMAT() 其中，50、人分别为未锈蚀钢筋以及锈蚀钢筋的屈服强度：”为钢筋锈蚀率。 钢筋锈蚀生成的锈蚀产物在钢筋与混凝土界面之间累积并造成保护层开裂，钢筋与混凝土之间 的粘结性能发生退化。研究者对于锈蚀后钢筋与混凝土之间的粘结性能展开了大量工作，建立了一 系列有效的锈蚀钢筋混凝土粘结滑移关系B2,。Bhargaval3]等提出了适用于锈蚀钢筋与混凝土粘结 强度计算的经验模型，并与已有实验数据吻合良好。本文采用文献[33]建立的锈独粘结滑移关系来 考虑锈蚀对钢筋与混凝土之间粘结强度的影响，表达式如下： t(n)=t(0(m≤1.5%) t()=t(0)e-19.8-15%(>1.5%) MERGEFORMAT ( 其中，t()为钢筋锈蚀率为%时钢筋与混凝土之间的粘结强度 @)为未锈蚀钢筋与混凝土之间的 粘结强度。 1.5边界条件及荷戴 混凝土保护层锈胀开裂阶段：考虑钢筋截面锈蚀的不均匀性，如式(4)所示的不均匀径向位移荷 载施加于钢筋周围的混凝土，进而模拟钢筋的不均匀锈胀。考虑到自然环境下钢筋锈蚀受到多种 条件的影响，且锈蚀位置具有较大的随机性。为了简化补算，假定钢筋截面内的锈蚀在钢筋长度方 向上一致。根据锈蚀率与锈胀位移的关系式6。通过锈蚀率反算锈胀位移大小，实现不同锈蚀率 下的钢筋锈蚀膨胀模拟工作。 锈蚀构件力学行为模拟阶段：钢筋混凝土梁两端采用铰支座，与验证试验2一致。将混凝土锈 胀开裂的“最终结果”作为构件力学行为研究的“初始条件”，在已有保护层开裂基础上对锈蚀构 件进行三点加载模拟，直至构件破坏。 2数值分析横型险证 Han等Po开展了无腹筋混凝梁通电加速锈蚀以及三点弯加载试验，对不同纵筋锈蚀率下钢筋 混凝土梁开裂以及梁破坏后的抗剪承载力和破坏模式进行对比分析。为验证本文数值模拟方法的合 理性，本文对Han第试验中的试件NS-0、NS-8、NS-l5进行了模拟，相应的锈蚀率分别为 0%、4.35%、8%。模拟结果与试验结果的破坏模式如图6所示。试验梁在完成加速锈蚀试验后，梁 侧面和底面均产生了不同程度的纵向裂纹，由于部分箍筋锈蚀，梁表面产生少量横向裂纹。由图6 可知，通过钢筋锈张模拟，梁底面和侧面出现沿纵筋方向的锈蚀裂缝，且随锈蚀率增大，锈蚀损伤 更为严重。荷载作用下，未锈蚀梁在破坏时产生了明显的剪切裂缝，且在靠近支座位置观察到了水 平裂缝。对于锈蚀梁，由于纵筋锈蚀的影响，部分斜裂缝不再从梁底部产生而在纵向锈蚀裂缝基础 上逐渐向加载点方向延伸，最终构件破坏。模拟结果与试验结果的破坏模式基本吻合。另外，对比 图7荷载位移曲线可以发现，模拟与试验获得的曲线形式基本一致仅在下降段略有差距。上述结果 说明，本文采用的多阶段数值模拟方法能够较好再现试验现象，从而证明了该方法的合理性。 收离日期2021-06-xx 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： 1.4.2 锈蚀影响 锈蚀不仅导致了钢筋截面积的降低，同时造成了钢筋性能的降低。Sun 等[30]对锈蚀钢筋进行单 轴拉伸试验，发现锈蚀使钢筋的屈服与极限强度均有所下降，且强度的降低与锈蚀率增加近似呈线 性关系。钢筋的性能通常用截面积、弹性模量、屈服和极限强度等参数来描述。本文模拟分析中， 暂不考虑锈蚀对钢筋弹性模量的影响，锈蚀后钢筋的屈服强度参考文献[31]取值，如式(7)。 yc 0 (0.986 1.038 ) y f f 88\* MERGEFORMAT () 其中，fy0、fyc分别为未锈蚀钢筋以及锈蚀钢筋的屈服强度；η 为钢筋锈蚀率。 钢筋锈蚀生成的锈蚀产物在钢筋与混凝土界面之间累积并造成保护层开裂，钢筋与混凝土之间 的粘结性能发生退化。研究者对于锈蚀后钢筋与混凝土之间的粘结性能展开了大量工作，建立了一 系列有效的锈蚀钢筋混凝土粘结滑移关系[32,33]。Bhargava[33]等提出了适用于锈蚀钢筋与混凝土粘结 强度计算的经验模型，并与已有实验数据吻合良好。本文采用文献[33]建立的锈蚀粘结滑移关系来 考虑锈蚀对钢筋与混凝土之间粘结强度的影响，表达式如下： 19.8( 1.5%) ( ) (0)( 1.5%) ( ) (0)e ( 1.5%) 99\* MERGEFORMAT () 其中，τ(η)为钢筋锈蚀率为 η/%时钢筋与混凝土之间的粘结强度，τ(0)为未锈蚀钢筋与混凝土之间的 粘结强度。 1.5 边界条件及荷载 混凝土保护层锈胀开裂阶段：考虑钢筋截面锈蚀的不均匀性，如式(4)所示的不均匀径向位移荷 载施加于钢筋周围的混凝土，进而模拟钢筋的不均匀锈胀[16]。考虑到自然环境下钢筋锈蚀受到多种 条件的影响，且锈蚀位置具有较大的随机性。为了简化计算，假定钢筋截面内的锈蚀在钢筋长度方 向上一致[18]。根据锈蚀率与锈胀位移的关系式(6)，通过锈蚀率反算锈胀位移大小，实现不同锈蚀率 下的钢筋锈蚀膨胀模拟工作。 锈蚀构件力学行为模拟阶段：钢筋混凝土梁两端采用铰支座，与验证试验[20]一致。将混凝土锈 胀开裂的“最终结果”作为构件力学行为研究的“初始条件”，在已有保护层开裂基础上对锈蚀构 件进行三点加载模拟，直至构件破坏。 2 数值分析模型验证 Han 等[20]开展了无腹筋混凝土梁通电加速锈蚀以及三点弯加载试验，对不同纵筋锈蚀率下钢筋 混凝土梁开裂以及梁破坏后的抗剪承载力和破坏模式进行对比分析。为验证本文数值模拟方法的合 理性，本文对 Han 等[20]试验中的试件 NS-0、NS-8、NS-15 进行了模拟，相应的锈蚀率分别为 0%、4.35%、8%。模拟结果与试验结果的破坏模式如图 6 所示。试验梁在完成加速锈蚀试验后，梁 侧面和底面均产生了不同程度的纵向裂纹，由于部分箍筋锈蚀，梁表面产生少量横向裂纹。由图 6 可知，通过钢筋锈胀模拟，梁底面和侧面出现沿纵筋方向的锈蚀裂缝，且随锈蚀率增大，锈蚀损伤 更为严重。荷载作用下，未锈蚀梁在破坏时产生了明显的剪切裂缝，且在靠近支座位置观察到了水 平裂缝。对于锈蚀梁，由于纵筋锈蚀的影响，部分斜裂缝不再从梁底部产生而在纵向锈蚀裂缝基础 上逐渐向加载点方向延伸，最终构件破坏。模拟结果与试验结果的破坏模式基本吻合。另外，对比 图 7 荷载位移曲线可以发现，模拟与试验获得的曲线形式基本一致仅在下降段略有差距。上述结果 说明，本文采用的多阶段数值模拟方法能够较好再现试验现象，从而证明了该方法的合理性。 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 录用稿件，非最终出版稿
工程科学学报 DOI: Damage 0 0.99 Side surface Exp.4.35% Sim.4.35% P Exp.0% Sim.0% Bottom surface 18 Exp.4.35% P Sim.4.35% 1 Side surface Exp.8% Sim.8% 02 Ep.8% Sim.8% Bottom surface 35 ■6模拟与试验结果梁破坏模式对比 Fig.6 Comparison of distribution of typical cracks in experimental and simula 150 120 n=0% 90 60 0 圆7构件荷载位移曲线对比 Fig.7 Comparison between the simulated mid-span deflection with the experimental 3纵筋畅触的无腹筋混凝土梁抗剪性能 3.1开破坏过程 3.1.1保护层锈胀开裂 图8给出了纵向钢筋锈蚀4.35%时混凝土的锈胀压力以及保护层开裂破坏的发展过程。由图可 知， 钢筋锈蚀导致混凝生产生损伤。这是由于钢筋锈蚀产物体积大于所消耗的铁的体积，从而使得 钢筋周围的混凝土受到压应力。在压应力作用下，混凝土内部的界面过渡区首先发生破坏。在混凝 土内部开裂前，人混凝土受到的锈胀压力达到最大值。混凝土开裂使钢筋与混凝土之间产生应力释放， 混凝土受到的压应力逐渐减小。随着锈蚀的继续发展，混凝土的内部裂缝贯通，并且破坏斜向混凝 土外表面发展，混凝士保护层的破坏面积逐渐增大。 3.1.2梁的破坏过程 在混凝土保护层锈蚀损伤的基础上，钢筋混凝土梁在静载条件下的破坏过程如图9。由于纵筋 锈蚀的影响，混凝土梁表面产生了明显的纵筋裂缝。在混凝土梁加载之前，梁的内部存在初始损伤 (纵筋裂缝)。荷载作用下，一方面纵筋裂纹继续发展，另一方面，在纵筋裂缝处产生向加载点方 向延伸的斜裂缝。随荷载继续增加，梁最终发生剪切破坏。 收精日期2021-06-Xx; 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： P P P Exp. 0% Sim. 0% P Exp. 4.35% Sim. 4.35% P Exp.8% P Sim. 8% 0 0.99 Damage Side surface Bottom surface Side surface Bottom surface Exp. 4.35% Sim. 4.35% Exp.8% Sim. 8% 图 6 模拟与试验结果梁破坏模式对比 Fig.6 Comparison of distribution of typical cracks in experimental and simulated failure modes. 0 2 4 6 8 10 0 30 60 90 120 150 Load/kN Deflection/mm Exp. Sim. η = 0% η = 4.35% η = 8% 图 7 构件荷载位移曲线对比 Fig.7 Comparison between the simulated mid-span deflection with the experimental 3 纵筋锈蚀的无腹筋混凝土梁抗剪性能 3.1 开裂破坏过程 3.1.1 保护层锈胀开裂 图 8 给出了纵向钢筋锈蚀 4.35%时混凝土的锈胀压力以及保护层开裂破坏的发展过程。由图可 知，钢筋锈蚀导致混凝土产生损伤。这是由于钢筋锈蚀产物体积大于所消耗的铁的体积，从而使得 钢筋周围的混凝土受到压应力。在压应力作用下，混凝土内部的界面过渡区首先发生破坏。在混凝 土内部开裂前，混凝土受到的锈胀压力达到最大值。混凝土开裂使钢筋与混凝土之间产生应力释放 ， 混凝土受到的压应力逐渐减小。随着锈蚀的继续发展，混凝土的内部裂缝贯通，并且破坏斜向混凝 土外表面发展，混凝土保护层的破坏面积逐渐增大。 3.1.2 梁的破坏过程 在混凝土保护层锈蚀损伤的基础上，钢筋混凝土梁在静载条件下的破坏过程如图 9。由于纵筋 锈蚀的影响，混凝土梁表面产生了明显的纵筋裂缝。在混凝土梁加载之前，梁的内部存在初始损伤 （纵筋裂缝）。荷载作用下，一方面纵筋裂纹继续发展，另一方面，在纵筋裂缝处产生向加载点方 向延伸的斜裂缝。随荷载继续增加，梁最终发生剪切破坏。 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 录用稿件，非最终出版稿
工程科学学报 DOI: 2 Ponit.1P2P3P.4 Damage 0 10.99 10 B 6 14 A:Generating B:Penetrating C.:Completing 0 cracks internal corrosion 0.0 0.2 0.40.6 0.8 1.0 圈8纵筋锈蚀4.35%下混凝土保护层破坏过程 Intinal 版稿 Fig.8 The cracking process of concrete cover under the corrosion level e Damage 0 0.99 ■9纵筋锈蚀43%混凝土梁破坏过程 Fig.9 The cracking process of RC beam with the corrosion level of 4.35%. 3.2娇烛率形响 本节给出了不同锈蚀率下保护层和梁的最终破坏模式，进而从宏观上分析锈蚀率对梁力学行为 的影响。图10和图11给出了不同锈烛率下混凝土保护层的破坏模式以及锈胀压力的分布。图12和 图13分别为不同锈蚀率下无腹筋混凝土梁的破坏模式和荷载位移曲线。 从图10可以看出， プ货面驳随锈蚀率增大，保护层破坏面积逐渐增大。在锈蚀率达到12%时， 混凝土保护层破坏严重，这后即使钢筋继续锈蚀，保护层破坏面积不再产生较大变化。另一方面， 混凝土保护层的开裂模式基本不变，破坏位置多集中相邻钢筋之间以及靠近混凝土外表面的钢筋斜 上方向。对应图中混凝土的锈胀压力分布，由于钢筋中间位置以及靠近混凝土外表面的钢筋斜 上方向混凝士破坏严重，此处混凝土能够承受的锈胀压力要小于其他位置。混凝土受到的锈蚀膨胀 是一个复杂的学行为，压力的大小与钢筋锈蚀程度、混凝土内部结构以及外部约束条件等多个因 素有关。因此，混凝土截面内锈胀压力的分布表现出较大的不均匀性。 Damage 0 ■0.99 (a)=4.35% (b)n=8% (c)=12% 收精日期：2021-06-xx 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： A B C 1 2 4 3 A:Generating cracks B:Penetrating internal C:Completing corrosion Damage 0 0.99 0.0 0.2 0.4 0.6 0.8 1.0 0 2 4 6 8 10 12 Ponit.1 P.2 P.3 P.4 Pressure/MPa uθ / umax 图 8 纵筋锈蚀 4.35%下混凝土保护层破坏过程 Fig.8 The cracking process of concrete cover under the corrosion level of 4.35%. Shear damage Shear cracks Intinal damage 0 0.99 Damage 图 9 纵筋锈蚀 4.35%下混凝土梁破坏过程 Fig.9 The cracking process of RC beam with the corrosion level of 4.35%. 3.2 锈蚀率影响 本节给出了不同锈蚀率下保护层和梁的最终破坏模式，进而从宏观上分析锈蚀率对梁力学行为 的影响。图 10 和图 11 给出了不同锈蚀率下混凝土保护层的破坏模式以及锈胀压力的分布。图 12 和 图 13 分别为不同锈蚀率下无腹筋混凝土梁的破坏模式和荷载位移曲线。 从图 10 可以看出，一方面，随锈蚀率增大，保护层破坏面积逐渐增大。在锈蚀率达到 12%时， 混凝土保护层破坏严重，之后即使钢筋继续锈蚀，保护层破坏面积不再产生较大变化。另一方面， 混凝土保护层的开裂模式基本不变，破坏位置多集中相邻钢筋之间以及靠近混凝土外表面的钢筋斜 上方向。对应图 11 中混凝土的锈胀压力分布，由于钢筋中间位置以及靠近混凝土外表面的钢筋斜 上方向混凝土破坏严重，此处混凝土能够承受的锈胀压力要小于其他位置。混凝土受到的锈蚀膨胀 是一个复杂的力学行为，压力的大小与钢筋锈蚀程度、混凝土内部结构以及外部约束条件等多个因 素有关。因此，混凝土截面内锈胀压力的分布表现出较大的不均匀性。 (a) η = 4.35% (b) η = 8% (c) η = 12% Damage 0 0.99 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 录用稿件，非最终出版稿
工程科学学报 DOI: 圆10不同锈蚀率下保护层破坏模式 Fig.10 Failure patterns of the cross-section of RC beams with various corrosion levels. 12 60 6 18 1=4.35% edw/ainssald 8260 62 18 8% 18 12 EdW/ainssald 606 12 18 越版稿 .Generating internal Generating exteral crack -Penetrating protective laye Completing corrosion 圆11不同锈蚀率下保护层锈胀压力分布 Fig.11 Distribution of the rust pressures of concrete cover with various corrosion levels. 图12为不同锈蚀率下梁的破坏模式。随着纵筋锈蚀率增大，外力所做的功更多由于局部裂缝 宽度增大而消耗，裂缝数量变化不明显。因此，随着锈蚀的增加，梁破坏时产生的裂缝数量逐渐减 少。由图13可知，梁的刚度以及承载力随锈蚀率增大而降低。当纵筋锈蚀率从0%增大到12%时， 梁的承载力下降了约6？%。纵筋锈蚀严重时，由于钢筋强度与钢筋与混凝土之间的粘结性能产生较 大损失，结构的协同江作的机制被削弱，梁的承载力下降严重。 收精日期2021-06-Xx 善盒项目：国家自然科学基金(51822801)：国家重点研发计划项目(2019Y℉C1511003)
工程科学学报 DOI： 图 10 不同锈蚀率下保护层破坏模式 Fig.10 Failure patterns of the cross-section of RC beams with various corrosion levels. 0 6 12 18 6 12 18 0 6 12 18 6 12 18 0 6 12 18 6 12 18 Generating internal crack Generating external crack Penetrating protective layer Completing corrosion Pressure/MPa η = 4.35% η = 8% η = 12% Pressure/MPa Pressure/MPa 0° θ 图 11 不同锈蚀率下保护层锈胀压力分布 Fig.11 Distribution of the rust pressures of concrete cover with various corrosion levels. 图 12 为不同锈蚀率下梁的破坏模式。随着纵筋锈蚀率增大，外力所做的功更多由于局部裂缝 宽度增大而消耗，裂缝数量变化不明显。因此，随着锈蚀的增加，梁破坏时产生的裂缝数量逐渐减 少。由图 13 可知，梁的刚度以及承载力随锈蚀率增大而降低。当纵筋锈蚀率从 0%增大到 12%时， 梁的承载力下降了约 67%。纵筋锈蚀严重时，由于钢筋强度与钢筋与混凝土之间的粘结性能产生较 大损失，结构的协同工作的机制被削弱，梁的承载力下降严重。 收稿日期:2021-06-xx; 基金项目:国家自然科学基金(51822801)；国家重点研发计划项目(2019YFC1511003). 录用稿件，非最终出版稿