钢筋混凝土结构锈蚀开裂模型研究.pdf

第卷第期 年 月铁 道 学 报 文章编号：（）钢筋混凝土结构锈蚀开裂模型研究陈华鹏，刘昌雨，江 钰，肖林发（华东交通大学 交通运输工程学院，江西 南昌）摘 要：为揭示钢筋混凝土结构锈蚀开裂的机理和裂缝相关参数的发展规律，对钢筋锈蚀作用下混凝土的开裂全过程进行研究，并建立混凝土锈蚀开裂模型。模型采用混凝土材料受拉指数型软化形式和黏结裂缝理论，将锈蚀开裂过程分为保护层未开裂、部分开裂和完全开裂三个阶段。推导每一阶段对应的应力状态、径向位移和裂缝状态等力学参数的表达式，获得计算混凝土保护层裂缝宽度的控制方程，并且给出求解该方程的数值计算方法。基于所建立的模型，研究裂缝从钢筋黏结表面扩展到混凝土保护层表面的全过程行为，讨论裂缝宽度、环向应力等参数的变化规律，预测混凝土保护层表面开裂时间与相应的临界锈蚀率。研究结果表明：钢筋表面和保护层表面的裂缝宽度的差值随锈蚀时间逐渐减小并最后趋于零；混凝土裂缝宽度和锈蚀率之间表现为正线性相关关系；混凝土保护层厚度与钢筋直径之比、混凝土抗拉强度和锈蚀膨胀系数等因素直接影响着保护层锈蚀开裂时间；最后，基于黏结裂缝理论建立的混凝土结构锈蚀开裂模型能够有效地预测试验值，可为钢筋混凝土结构锈蚀开裂机理研究提供依据。关键词：混凝土开裂；钢筋锈蚀；开裂时间；裂缝扩展；残余强度中图分类号：；文献标志码：收稿日期：；修回日期：基金项目：国家重点研发计划（）；国家自然科学基金（，）；江西省科技合作专项重点项目（）；江西省重点研发计划（）第一作者：陈华鹏（），男，浙江台州人，教授，博士。：通信作者：陈华鹏（），男，浙江台州人，教授，博士。：，（，）：，：，：；第 期陈华鹏等：钢筋混凝土结构锈蚀开裂模型研究 钢筋和混凝土材料良好的黏结特性使得两种不同性质的建筑材料可以有效的协同工作，同时也让结构充分发挥了材料的承载性能。但是，当钢筋发生锈蚀时，混凝土结构的耐久性和承载性能就会受到严重的影响。报告显示，美国每年因锈蚀问题而产生的桥梁维修费用高达 亿美元，其中有近 的损失是由混凝土桥的钢筋锈蚀问题造成的。在英国，部分桥梁因锈蚀破坏造成的维修费用达到了造价的 倍。我国在基础设施建设过程中，钢筋混凝土结构耐久性问题也逐渐暴露，结构未达到设计年限就提前退役的现象时有发生，例如我国青海省盐湖地区的桥梁结构，使用仅 年就因严重腐蚀致使结构过早失效，其原因是钢筋表面的钝化膜被破坏，导致钢筋暴露在易腐蚀的环境中，并引起钢筋的锈蚀，而锈蚀物向外膨胀的过程会使混凝土保护层产生环向拉应力，当混凝土环向拉应力达到抗拉强度时混凝土就会开裂，产生的裂缝将进一步加速钢筋的锈蚀，使保护层剥落、结构承载性能退化，直至结构性能失效。因此，钢筋锈蚀所致混凝土的开裂问题是影响混凝土结构耐久性和承载性能的关键，揭示锈蚀引起混凝土的开裂机理对于工程结构寿命评估和维修加固等工作具有实际指导意义。国内外学者针对混凝土开裂模型、钢筋锈胀力、混凝土性能退化等问题，开展了大量的试验、数值以及理论研究。在试验研究方面：等针对钢筋位置、锈蚀速率、钢筋直径和保护层厚度等变量进行了加速锈蚀试验，建立了锈蚀深度和时间以及锈蚀速率的经验公式；等通过试验研究得到混凝土开裂的主要影响因素包括保护层厚度、钢筋尺寸、混凝土材料属性以及锈蚀速率等；等认为可以利用应变来判断开裂情况，并且建立临界锈蚀率和保护层厚度的幂函数经验公式；等进行锈蚀加速试验，引入混凝土抗拉强度、保护层厚度与钢筋直径的比值作为变量，得到保护层表面裂缝宽度和锈蚀深度的经验公式。在数值仿真方面：等建立混凝土锈蚀损伤的通用数值模型，分析混凝土的锈裂过程及性能退化；张伟平等利用结构的温度膨胀特性模拟钢筋锈蚀物的膨胀作用，建立锈胀仿真计算模型，得到混凝土保护层与钢筋直径之比是影响混凝土开裂的重要因素；金伟良等建立非线性有限元模型，分析均匀锈蚀和非均匀锈蚀锈胀力作用下结构的开裂过程，计算开裂时刻的锈胀力大小，并提出临界锈蚀率的计算方法。在理论研究方面：提出海洋环境下锈蚀钢筋混凝土结构的数学物理模型，研究锈蚀过程中氧、氯离子的扩散和锈蚀物的组成，简化了保护层锈胀开裂的计算方法；等将混凝土开裂考虑为各向异性的轴对称问题，建立一种新的黏聚裂缝模型，分析裂缝的扩展全过程，预测了保护层表面的裂缝宽度变化；伦培元等考虑在混凝土存在初始缺陷的情况下，提出三维半椭球形的混凝土结构锈胀开裂理论模型，认为增加混凝土保护层厚度有助于延缓结构的锈胀开裂，从而提高结构的耐久性；王海龙等基于断裂力学理论建立预测开裂时刻和临界锈蚀率的混凝土锈蚀开裂模型。以上研究较好地分析了结构锈胀开裂过程中锈胀力的变化，预测了保护层开裂的时刻和临界锈蚀率，探讨了混凝土锈蚀开裂的主要影响因素，但是对混凝土开裂全过程以及裂缝扩展的理论分析较少，无法揭示裂缝的演变机理。为此，本文基于厚壁圆筒应力理论，假设开裂混凝土的软化曲线为指数形式，提出了一种新的由钢筋锈蚀引起的混凝土保护层开裂模型。该模型将混凝土锈蚀开裂过程分为保护层未裂、部分开裂和完全开裂三个阶段，推导了每一阶段对应的应力状态、径向位移和裂缝状态等力学参数的表达式，得到了求解裂缝宽度的控制方程。此外，本文还研究了裂缝前沿的扩展规律，预测了混凝土保护层开裂的时间和锈蚀深度，揭示了裂缝的演化过程以及混凝土的刚度、强度变化规律。最后将所得结果同已有的试验数据对比分析，验证本文提出的模型。锈蚀裂缝模型最早利用弹性力学中的厚壁圆筒理论建立了模拟混凝土保护层的模型，分析开裂混凝土与钢筋的残余黏结强度，该模型也成为了研究锈蚀混凝土开裂问题的主要方法之一。本文基于混凝土材料开裂后呈各向异性的特性以及钢筋均匀锈蚀的假设，建立了三阶段厚壁圆筒开裂混凝土模型。采取钢筋均匀锈蚀假设的原因主要有：钢筋混凝土在均匀的锈胀力作用下可以等效为受均匀内压作用的厚壁圆筒结构，将机理复杂的钢筋锈蚀导致的混凝土开裂问题考虑为轴对称问题，可以简化计算过程；目前国内外学者为了研究混凝土锈蚀开裂问题而建立的锈裂模型大多基于均匀锈蚀的假设，且同试验数据吻合较好；自然情况下的非均匀锈蚀分布形态受到环境、人为因素影响较大，目前仍没有一个准确的非均匀锈蚀分布预测模型，关于非均匀锈蚀导致的混凝土开裂的研究也较为缺乏，相关试验与实测数据不足。建立的混凝土保护层锈蚀开裂模型见图，其中，为钢筋的初始半径，为裂缝前端半径，为保护铁 道 学 报第 卷层厚度，为包含保护层厚度 的外部圆的半径，为裂缝的前端半径，为钢筋直径，包围着 的圆环代表钢筋被消耗掉一部分而产生的锈蚀氧化物。图 混凝土保护层锈蚀开裂模型 锈蚀过程钢筋锈蚀主要是因为钢筋表面钝化膜被破坏，氯离子侵入混凝土内部，使得钢筋生锈，而根据锈蚀程度的不同，锈蚀物也不同，生成的锈蚀产物导致质量发生变化，单位长度的锈蚀物质量 为（）（）式中：（）为单位长度钢筋表面的年平均锈蚀电流；为锈蚀时间；为经验系数取 。根据锈蚀物的质量，推算钢筋的体积变化 为 （）式中：为单位长度内被消耗的钢筋质量；、分别为锈蚀物、钢筋的密度；为锈蚀产物膨胀系数，研究表明，锈蚀物可能产生 倍的体积膨胀。锈蚀物会随着体积增大开始向外膨胀，导致钢筋表面产生径向位移，假设 且钢筋均匀锈蚀，可以表示为（）式中：为钢筋与混凝土表面间的孔隙厚度。工程实际中均匀锈蚀的情况出现较少，更多的以点蚀或坑蚀为主。但是当锈蚀分布系数不大于 时，锈蚀产物接近均匀分布，既有研究表明，在研究裂缝扩展和结构残余寿命的问题时，采用均匀锈蚀或非均匀锈蚀两种假设造成的差异较小。因此，均匀锈蚀的假设在本文研究中是合理的。混凝土锈蚀开裂模型基于钢筋均匀锈蚀的假设，模型可以简化为轴对称模型。同时混凝土在工作条件下往往产生弯曲裂纹，可以将该问题简化为平面应力问题。因此，本文模型中的环向应力为主拉应力，径向应力为主压应力。当混凝土开裂后，由于混凝土的软化特性，材料的抗拉强度会随着裂缝宽度的增大而减小，根据裂缝和抗拉强度的关系引入材料的拉伸软化曲线，其主要形式有直线、双线性及指数型曲线。直线型软化曲线由于较为理想化，与试验结果差异较大。目前大多数研究主要是基于双线性的假设，而针对指数型软化曲线的研究较少。由于指数型曲线比较符合混凝土抗拉特性，因此本文利用指数型软化曲线，建立钢筋锈蚀混凝土开裂模型。软化曲线是材料的断裂能、裂缝宽度 和抗拉强度 等参数的函数，这些参数可通过试验测得，再根据断裂能公式（），其中 为结构失效时极限裂缝宽度，最后拟合出拉伸软化曲线。大量试验研究表明，不同骨料尺寸及混合比对软化曲线的影响比较有限，而指数型软化曲线被认为是较为合理的近似曲线，其定义为（）式中：为混凝土残余抗拉强度；为软化参数。为规范化裂缝宽度，与实际裂缝宽度 的转换关系为（）不同软化参数对应的指数拉伸软化曲线与典型的双线性曲线的形状对比见图。图 不同软化参数 的指数软化曲线与双线性软化曲线的对比由于研究对象是结构沿纵向钢筋的截面，而且由于锈蚀深度相对钢筋的初始半径 非常小，因此单位锈蚀深度 的求解可以进一步简化为单位长度锈蚀产物的体积除以其表面积求得（）引入锈蚀率 的概念，并定义其为单位锈蚀深度与原钢筋半径 的比值，即第 期陈华鹏等：钢筋混凝土结构锈蚀开裂模型研究（）上述的钢筋锈蚀过程忽略了结构所承受外力产生的影响。裂缝扩展全过程分析钢筋锈蚀物向外膨胀，使得混凝土产生径向位移，当混凝土的环向拉应力达到混凝土抗拉强度时，混凝土开裂，产生径向裂缝。随着时间的推移，钢筋锈蚀程度加深，裂缝逐渐加宽且向外扩展，直至混凝土保护层表面。混凝土裂缝从黏结表面扩展到保护层表面的过程中，可以分为开裂部分和未裂部分。未裂部分为各向同性问题并运用弹性理论进行分析。开裂部分的径向刚度保持不变，切向刚度受开裂影响而降低，因此可以将开裂混凝土视为正交各向异性材料，通过混凝土材料拉伸软化曲线，引入残余抗拉强度，结合断裂力学计算断裂应变及应力，考虑钢筋锈蚀产生位移边界条件，联立得到各向异性弹性体的位移协调方程，最终得到裂缝宽度的控制方程。整个混凝土开裂过程可分为裂缝萌生、部分开裂及完全开裂三个阶段。当钢筋锈蚀物填充满混凝土孔隙区后，继续向外膨胀，当钢筋表面混凝土的环向应力达到混凝土抗拉强度时，混凝土开裂，进入部分开裂阶段。随着时间的推移，锈蚀物逐渐增加，混凝土保护层的环向拉应力越来越大，裂缝也随之向外扩展，引入裂缝前端半径。当裂缝扩展到保护层表面时（），结构表面开始出现裂缝，利用钢筋表面混凝土裂缝宽度 及其对径向坐标 的一阶导数值，求解裂缝控制方程式，从而得到保护层表面裂缝宽度变化。裂缝萌生此时的圆筒模型没有产生裂缝，可以采用弹性力学理论，因此各向同性弹性体的位移随径向坐标 的关系为 （）式中：为径向位移；、为由边界条件确定的常系数；径向应力 和环向应力 可表示为 （）（）式中：为混凝土材料的泊松比。钢筋表面处的位移和混凝土保护层的径向应力边界条件为 （）式中：为黏结面位移。将边界条件式（）分别代入式（）和式（）中，解得（）（）（）（）（）（）（）（）将解得的常系数代入式（）中，得到径向应力 和环向应力 表达式为（）（）（）（）（）（）部分开裂阶段此阶段，圆筒内表面已经开裂，裂缝前端向外延伸扩展，但是还未扩展至保护层表面，可以将开裂模型分为开裂部分的内筒（）和弹性部分的外筒（）。外筒部分的边界条件为 （）式中：为混凝土初始抗拉强度。利用边界条件式（）和式（）、式（），径向应力 和环向应力 分别为 （）（）裂缝前端的表达式可以通过裂缝前端处未裂区与开裂区的变形协调关系得到，已知弹性区与开裂区在裂缝尖端处（）的径向位移分别为 （）（）（）式中：为混凝土弹性模量。将式（）、式（）代入式（）整理得 （）（）（）（）（）联立式（）和式（），利用位移协调关系即 ，且裂缝前端处裂缝宽度为，求得裂缝前端半径 的表达式为铁 道 学 报第 卷（）（）（）（）（）完全开裂阶段将开裂后的混凝土视为正交各向异性材料，为了方便求解，引入黏聚裂缝的概念，将实际的径向裂缝等效为对称分布在模型内的 条裂缝，构成一个轴对称问题。基于 等提出的断裂带理论，混凝土保护层表面的裂缝总数 为（）式中：为最小允许裂缝带宽度 ，为混凝土的最大骨料尺寸。开裂模型的环向总应变 由断裂应变 及线弹性应变 两部分组成，断裂应变由 条裂缝产生，线弹性应变则与残余抗拉强度有关，分别表示为（）（）式中：，为 等定义的特征长度，其计算式为（）因此，环向总应变 为 （）根据几何方程可得径向位移 和径向应变 分别为（）将式（）、式（）、式（）和式（）代入式（）得径向应变 和径向位移 分别为（）（）（）（）开裂混凝土的环向刚度会随着裂缝宽度的增加而逐渐减小，引入刚度折减系数 为 （）根据 等提出的正交各向异性弹性柱体的位移协调方程，即 （）将式（）和式（）代入式（），得到裂缝宽度关于半径 的控制方程为（）（）式（）是一个非线性二阶微分方程，难以直接求得规范化裂缝宽度 的显式公式。本文利用数值方法求解该微分方程，首先将方程降阶，取得某一点的 值及其一阶导数值，运用这两个初始边界条件，可以得到方程的数值解，再通过式（）的转换就可以得到实际裂缝的宽度。在裂缝前端到达保护层表面时，模型进入完全开裂阶段，结构表面开始出现裂缝。此时保护层表面的裂缝宽度可利用钢筋表面处的裂缝宽度 及其一阶导数值 两个初始边界条件求解式（）得到。联立式（）和式（），得到钢筋表面裂缝宽度为 ，（）式中：函数（，）为 的反函数，常应用于求解带有指数项的微分方程。开裂后混凝土可以看作正交各向异性材料，其应力应变关系可表示为 ()（）()（）考虑开裂部分环向拉应力等于其残余抗拉强度，因此将式（）、式（）和式（）代入式（），再联立式（）和式（），解得 对于径向坐标 的一阶导数为 ()（）（）利用式（）、式（）作为二阶微分方程的两个初始条件求解裂缝宽度控制方程式（），得到裂缝宽度 与径向坐标 的关系，再令 得到混凝土保护层表面的裂缝宽度。数值算例 模型验证为验证本文理论模型的准确性，采用 等得第 期陈华鹏等：钢筋混凝土结构锈蚀开裂模型研究到的加速锈蚀试验数据进行对比验证。根据文献中试 验 提 供 的 材 料 参 数，混 凝 土 抗 压 强 度 ，抗拉强度 ，弹性模量 ，泊松比 ，混凝土蠕变系数为，钢筋密度 ，钢筋锈蚀物密度 ，取。假设最大骨料尺寸 ，软化参数，断裂能 ，裂缝总条数 ，极限裂缝宽度 ，体积膨胀系数 取 和，分别表示取值的不确定性与可能的取值范围，用以研究体积膨胀对裂缝扩展的影响。利用本文方法计算保护层完全开裂时间并与试验结果相比较，见表。由表 可知，保护层完全开裂的时间会随着保护层厚度的增大而明显增加，试验结果被包含在预测的时间区间之内或非常接近预测区间，表明预测结果与试验结果吻合较好，从而验证了本文模型的准确性和有效性。表 不同体积膨胀系数下保护层开裂时间钢筋直径 保护层厚度 锈蚀电流密度（）预测区间试验结果 参数分析不同膨胀系数下钢筋表面处的径向位移 随时间的变化关系见图。由图 可知，相同时间下锈蚀膨胀系数越大，径向位移也越大。在锈蚀初期，径向位移增加的速率较快，随着时间的推移曲线趋于平缓。锈蚀时间低于 时，不同膨胀系数对径向位移的影响不明显，超过该年限后，径向位移的差异显著增大。并且分别在 、处达到 ，对应混凝土表面裂缝宽度达到极限值。图 径向位移 随时间变化不同体积膨胀系数下钢筋表面处的裂缝宽度随时间的变化见图，由图 可知，钢筋表面处的裂缝宽度随时间的变化同径向位移相似，同样也是在 、时达到裂缝极限值。图 不同体积膨胀系数下钢筋表面处裂缝宽度随时间变化为了解钢筋表面和保护层表面的裂缝宽度的时间特征差异，绘制了 时两者随时间的变化曲线见图。由图 可知，保护层表面在初期并没有出现裂缝，当锈蚀时间达到 时，表面突然出现裂缝，且宽度瞬间升至 ，这是由于混凝土在开裂时存在应力瞬间释放现象。钢筋表面的裂缝则随时间的增长持续增大，且由于锈蚀物向外膨胀，其裂缝宽度要略大于保护层裂缝宽度。随着时间的推移，两者的差值会逐渐减小并趋于相同，最终在 时达到极限裂缝宽度值 。这主要是由于锈蚀时间达到一定时，锈蚀产物逐渐增多，锈胀力会持续增大，使得钢筋表面的裂缝逐渐向保护层裂缝贯通。图 裂缝宽度 随时间变化当 时，裂缝前端从钢筋表面扩展到保护层表面的全过程见图。由图 可知，裂缝前端在 左右的时间点到达了保护层表面。裂缝前端在接近钢筋表面和保护层表面时，其扩展速率明显快于中间部分，同样是由于在开裂初始阶段和最终阶段应力释放现象导致的。钢筋表面和保护层表面的环向应力随时间的变化曲线见图。钢筋表面由于开裂较早，所以在经历了短暂的初裂过程之后环向应力由混凝土抗拉强度值开始呈指数下降，直至接近为。保护层表面处环向应力由于受到膨胀作用，在 的时间内升高至 ，达到混凝土的抗拉强度，随后保护层表面开铁 道 学 报第 卷裂，并同样以指数形式急剧降低并逐渐趋于。图 裂缝前端 随时间变化图 环向应力随时间的变化（）为研究临界钢筋锈蚀深度（混凝土保护层开裂时锈蚀深度）同保护层厚度与钢筋直径之比（）的关系，根据文献提供的试验参数计算得到理论预测结果，见图。同相应试验结果相比，可以发现，本文模型的预测结果与试验值吻合较好，临界锈蚀深度与保护层和钢筋直径之比 的关系呈正线性相关规律。图 不同保护层厚度与钢筋直径之比的临界锈蚀深度等效裂缝宽度随锈蚀率的变化曲线见图。锈蚀水平较低时，无裂缝产生。当锈蚀率达到某一临界值时（约为 ），保护层表面完全开裂。裂缝达到 后，随着锈蚀率的增加裂缝宽度近似呈线性增长。进一步地，将本文的预测结果与大量试验结果进行对比，可以发现本文的模型预测结果与大部分试验结果吻合程度较好，也进一步地验证了本文所提出方法的准确性。混凝土表面裂缝宽度在不同锈蚀电流的情况下随时间的变化关系见图。由图 可知，锈蚀电流的图 裂缝宽度随锈蚀率的变化与已有试验结果比较增大会导致裂缝提前出现。随着锈蚀电流从 增加至 ，在相同的锈蚀时间下，裂缝宽度将会明显变大。此外，锈蚀电流越大，裂缝的增长速率也显著增大，尤其是锈蚀电流超过 以后。由此可见，锈蚀电流对混凝土开裂与裂缝扩展有着显著的影响。图 锈蚀电流对裂缝宽度变化的影响钢筋表面和保护层表面处刚度折减系数 随时间的变化见图。未开裂时，刚度折减系数被假定为。由图 可知，在钢筋表面，刚度折减系数 伴随着混凝土开裂急剧下降，在 时即保护层表面开裂的时刻就已经降到了 以下，而保护层表面的刚度从开裂后骤降至 左右，随后继续快速下降。因此，混凝土保护层在开裂过程中其刚度下降的速度要快于强度衰减的速度。在不同的锈蚀物膨胀系数下，混凝土结构的临界锈蚀深度随混凝土抗拉强度的变化见图。结果表明，锈蚀产物膨胀系数从 增大到 过程中，临界锈蚀深度有了明显的减小。因此，锈蚀产物的膨胀系数越大，结构的锈蚀破坏发生得越快。随着混凝土抗拉强度的增大，混凝土结构的临界锈蚀深度逐渐增大，但是随着锈蚀产物膨胀系数的增大，抗拉强度对临界锈蚀深度的影响会有所削弱，其曲线也相对比较平缓。第 期陈华鹏等：钢筋混凝土结构锈蚀开裂模型研究图 刚度折减系数随时间的变化图 不同膨胀系数下临界锈蚀深度随混凝土抗拉强度的变化 结论本文考虑混凝土开裂后的抗拉强度指数型软化形式，提出了基于时间及锈蚀率的混凝土锈蚀开裂模型，研究了混凝土锈胀开裂全过程，主要结论如下：（）推导了裂缝宽度沿着径向坐标 的控制方程并给出了数值求解方法，预测了混凝土保护层表面完全开裂的时间，并利用试验数据进行对比，验证了模型的有效性。（）研究了裂缝扩展全过程，推导了裂缝前端及各阶段相关力学参数的公式，利用数值算例研究了径向位移和裂缝宽度随时间的变化规律。（）在锈蚀初期，混凝土表面的环向应力会迅速上升，当裂缝到达保护层表面时，混凝土抗拉强度急剧下降，而切向刚度则将会以更快的速度衰减。（）混凝土保护层厚度与钢筋直径之比（）、混凝土抗拉强度和锈蚀膨胀系数等因素直接影响着保护层锈蚀开裂时间。锈蚀电流直接与锈蚀速率相关，工程应用中可以通过测定、控制锈蚀电流来解决相应的混凝土结构耐久性问题。（）通过引入锈蚀率的概念来描述钢筋锈蚀程度，发现锈蚀率和裂缝宽度之间呈正线性相关，并通过大量试验数据对比分析，发现试验结果与预测结果吻合程度较好，进一步验证了本文模型的有效性和准确性。（）利用保护层外表面的裂缝宽度，结合结构允许的最大裂缝宽度以及钢筋锈蚀率，可以预测结构剩余寿命；同时，结合时变结构失效风险分析与全寿命成本分析，可以指导混凝土结构的维修加固以及优化维护方案。参考文献：胡志坚，夏雷雷，程晨，等钢筋混凝土构件锈蚀开裂与锈胀力分析哈尔滨工业大学学报，（）：，（）：，（）：孙伟现代混凝土材料与结构服役特性的研究进展混凝土世界，（）：，（）：，：，（）：，（）：，：，：，（）：张伟平，张誉混凝土中钢筋锈胀过程的计算机仿真分析 同 济 大 学 学 报（自 然 科 学 版），（）：，（），（）：金伟良，赵羽习，鄢飞钢筋混凝土构件的均匀钢筋锈胀力的机理研究水利学报，（）：，（）：，（）：，铁 道 学 报第 卷，：伦培元，张小刚，张强，等基于初始缺陷形状的混凝土结构锈胀开裂理论模型东北大学学报（自然科学版），（）：，（），（）：王海龙，金伟良，孙晓燕基于断裂力学的钢筋混凝土保护层锈胀开裂模型水利学报，（）：，（）：赵羽习钢筋锈蚀引起混凝土结构锈裂综述东南大学学报（自然科学版），（）：（），（）：郑建军，周欣竹钢筋混凝土结构锈蚀损伤的解析解水利学报，（）：，（）：赵羽习，金伟良混凝土锈胀时刻钢筋锈蚀率的数值分析方法 浙 江 大 学 学 报（工 学 版），（）：，（），（）：，（）：，（）：，（）：，（）：，（）：吴乔长期荷载条件下混凝土断裂特性研究大连：大连理工大学，：，（）：，（）：，：，（）：，（）：，（）：，（）：，：，：，（）：陈华鹏，江钰，刘昌雨锈蚀钢筋混凝土柱残余承载力研究铁道学报，（）：，（）：曹国辉，曹一青，周超，等栓钉锈蚀后钢混凝土组合梁抗弯性能试验研究铁道学报，（）：，（）：李福海，靳贺松，胡丁涵，等锈蚀钢筋混凝土黏结滑移试验研究铁道学报，（）：，（）：（责任编辑 聂现会）