1、第 卷 第期 年 月广西大学学报（自然科学版）（）收稿日期：；修订日期：基金资助：国家“十三五”重点研发计划重点专项基金（、）；辽宁省重点研发计划指导计划基金（）；沈阳市建委城乡系统科研项目（）通讯作者：于玲（），女，陕西泾阳人，沈阳建筑大学教授；：。引文格式：包龙生，赵家康，张攀，等 装配式混凝土桥墩节点的抗冻性能研究 广西大学学报（自然科学版），（）：装配式混凝土桥墩节点的抗冻性能研究包龙生，赵家康，张攀，于玲，杨宇豪，包宇扬（沈阳建筑大学 交通与测绘工程学院，辽宁 沈阳 ；沈阳建筑大学 土木工程学院，辽宁 沈阳 ；沈阳工业大学 建筑与土木工程学院，辽宁 沈阳 ；沈阳大学 应用技术学院，辽
2、宁 沈阳 ）摘要：为了研究北方寒冷地区装配式桥墩节点的冻融损伤规律，基于热力学理论，建立适用于数值模拟的热分析表达式，结合实际工程，采用 有限元软件建立的缩尺模型，从应力、内外温差、塑性应变和位移等方面对装配式桥墩的不同部位开展冻融循环作用下的温度场和冻融损伤分析。结果表明：一次冻融循环作用下，墩帽、墩柱和承台的外表面温度应力与中心点相比，分别增加、，不同部位呈现一致的变化趋势，即中心点的温度应力变化趋势滞后于外表面；承台中心的温度极值波动最小，说明温度对尺寸较大的构件内部影响较小；随着冻融循环次数的增加，桥墩损伤越来越严重，次冻融循环作用下，桥墩节点的塑性应变和极限位移均大于其他部位。综合结
3、果来看，装配式桥墩节点是下部结构的薄弱部位。关键词：装配式桥梁；桥墩；冻融循环；数值模拟；二次开发中图分类号：文献标识码：文章编号：（），（，；，；，；，）：，广西大学学报（自然科学版）第 卷 ，：，：，：；引言在我国东北的严寒区域，冬季气温非常低并且早晚的温差大，导致装配式混凝土桥梁出现冻融损伤。装配式技术具有工期短、施工方便、安装精密度高等优点而逐渐被建筑业广泛应用，等，包龙生等指出装配式混凝土桥梁逐渐成为新型工业化发展的重要组成部分。随着建筑行业的快速发展，越来越多的学者开展了混凝土的冻害研究。宋玉普等进行了冻融循环后混凝土的力学性能方面的研究，建立了混凝土损伤后的单轴拉压数学公式。等使
4、用了超声影像跟 扫描相结合的技术观察冻融循环对混凝土的影响，进一步验证了冻融损伤分布的不均匀性。等研究了混凝土遭受冻融循环时的破坏规律，发现混凝土的应变曲线逐渐增加，说明混凝土的内部出现残余应变，导致损伤逐渐累积。等研究了混凝土受到冻融循环、冲击载荷同时作用时的能量耗散规律，发现混凝土具有从微观破坏到宏观断裂的特性。潘婷通过设计 种试验工况，开展冻融循环 干湿交替双因素作用时混凝土结构的耐久性研究，发现经过干湿交替 冻融循环的混凝土试件表层氯离子浓度比单独经历干湿交替的试件增长了 ，比单独经历冻融循环的试件增长了 ，耦合作用下的试件损伤程度比只经受干湿交替、冻融循环严重得多。等开展了考虑冻融损
5、伤影响的混凝土力学性能和氯离子渗透性能试验研究，通过设计潮汐循环模拟装置进行了一系列室内试验，结果表明冻融侵蚀造成的损害会随着氯离子环境下干湿循环次数的增加而放大，在相同的干湿循环次数下，混凝土表观氯离子扩散系数和表面氯离子浓度随冻融循环次数的增加而明显增加。目前混凝土结构的耐久性研究主要是以现浇结构为主，然而装配式混凝土桥墩采用分段预制，不同节段之间由灌浆套筒连接，节点接缝处采用坐浆料进行封闭，存在物理不连续的连接界面区，使得装配式桥墩节点部位更容易受到外部环境侵蚀。随着计算机技术的发展，采用数值模拟的方法开展耐久性研究成为可能 。等 通过建立二维有限元模型，预测了混凝土结构在冻融循环作用时
6、的膨胀量。等 研究了混凝土在冻融作用下微观层面的本构关系，采用数值模拟的方法对混凝土的失第期包龙生，等：装配式混凝土桥墩节点的抗冻性能研究效原因进行验证。等 通过数值模拟的方式，研究了饱和砂浆遭受冻结时的位移、压力、温度等变化规律。等 通过建立分离式的数值模型，开展混凝土在冻融循环作用下的裂缝发展模式及力学性能的研究，发现试样的尺寸和形状对裂缝发展模式有重要影响。吴文燕 为了研究混凝土遭受冻融循环影响时力学性能的变化情况，建立了有限元模型，模拟结果表明加载速率相同时，结构的塑性应变损伤和冻融次数是正相关的。等 采用数值模拟的方法对胶结材料内部因冻融作用造成的损伤影响进行了定量评估，结果表明孔隙
7、溶液中除冰盐的影响取决于饱和度，而临界饱和度取决于最低温度。本文结合学者的研究成果，针对装配式桥墩节点部位的特殊性，基于热力学理论，建立有限元热分析矩阵表达式，对 软件进行二次开发，能够将节点温度转化为积分点温度，同时满足冻融循环作用下损伤量的叠加计算，更加真实地反映实际情况，分析冻融循环作用对装配式桥墩节点部位的影响规律，为北方寒冷地区装配式桥墩的大规模推广与应用提供科学支撑。热传导分析结构经历冻融循环作用，受到稳定温度应力和瞬态温度应力双重作用，结构内部的温度分布情况比较复杂，导致传统试验方式难以精确测量温度场的分布特征，采用数值模拟的方式能够解决以上问题。由热力学第一理论，如式（）所示。
8、，（）式中：为热量；为外力产生的功；为动能；为内部能量；为势能。将模拟的对象离散并分解成个单元，即，之后选择内插法对每一个单元的温度场计算分析，有 ，代入泛函数如式（）所示。（）（）（），（）式中：（）（）；。将单元矩阵重组，可以得到全局矩阵。选取函数的极限数值，即，进而得到数值模拟所需要的矩阵形式，如式（）所示。，（）式中：为比热容的矩阵形式；为温度函数对时间进行一阶求导；为热量进行传导时的矩阵形式；为节点温度的向量形式；为温度载荷的向量形式。装配式混凝土桥墩的温度场分析 试件尺寸模型以东新高架桥下部结构工程为依托，下部结构的尺寸如下：盖梁高度为 ；承台高度为 ；墩柱的最大高度为 ，桥墩的横
9、断面均采用 的形式，节点之间采用灌浆套筒连接。比例尺采用，桥墩的横截面为 ，墩柱 高。墩帽的尺寸是 ，承台固定端要满足加载要求，其缩尺后的尺寸是 。桥墩主筋直径是，共根，箍筋直径，箍筋间距为 ，装配式桥墩和盖梁的拼接部位、桥墩和承台的拼接部位需要对箍筋进行加密，间距为，主筋和箍筋均是广西大学学报（自然科学版）第 卷 级钢筋。下部结构均采用强度等级为 混凝土，节点部位使用 型半灌浆套筒进行连接。装配式桥墩构造如图所示。（）立面图（）侧面图（）墩柱截面（）承台俯视图图桥墩尺寸 材料参数钢筋、混凝土的主要参数如表、所示。表钢筋的力学性能 密度（）线胀系数比热（）热传导系数（）弹性模量 泊松比 第期包
10、龙生，等：装配式混凝土桥墩节点的抗冻性能研究表混凝土的力学性能 密度（）线胀系数比热（）热传导系数 （）弹性模量 泊松比 建模过程分离式模型的特点是能够对混凝土和钢筋等结构进行独立建模，进而能够模拟混凝土与钢筋受力时产生的粘结 滑移现象，因此本文选择分离式的模型，混凝土采用实体单元（），钢筋单元采用桁架单元（），套筒及灌浆料均为轴对称实体单元（）。由于灌浆料强度高，因此选择强度等级为 混凝土的本构关系替代灌浆料，混凝土的本构关系为损伤塑性模型（），钢筋的本构关系选择三折线模型，钢筋与混凝土的约束形式选择嵌入式，各个混凝土部件之间采用面面接触。设置边界种子以及网格控制参数后进行网格划分，以均匀划
11、分网格为目的进行规律切割，其中半灌浆套筒模型进行局部加密，共设置 个单元，其中 单元 个。模拟的温度范围主要是我国北方寒冷地区，对北方寒区的装配式桥墩冻融损伤的研究有重大意义。设置环境温度为 ，温度范围为 ，通过设置幅值实现冻融循环，单次冻融的时间是，设置、次冻融循环。分析步需选择热传递，荷载与时间的关系设置成瞬态，相互作用需在热传递相应的分析步内进行设置，选择表面热交换的条件类型。半灌浆套筒模型如图所示，桥墩试件模型如图所示。图半灌浆套筒模型 （）钢筋骨架（）网格划分图桥墩的有限元模型 温度场分析为观察并研究外界温度变化对结构温度的影响，通过后处理对模型上相同水平路径的外表面处、套广西大学学
12、报（自然科学版）第 卷筒部位、墩帽的中心点处、墩柱的中心点处、承台的中心点等位置的 个冻融循环数据进行拟合，如图所示。依据图可得，桥墩冻融循环经历到第，墩帽的外表面和墩帽中心位置产生了最大的温度差值 ，相应的外表面及中心位置的温度数值分别是 、；与此同时，墩帽的外表面和套筒部位产生最大温差值 ，相应的外表面及套筒部位的温度数值分别是 、。当墩柱的外表面处和墩柱的中心点处在冻融后，产生最大的温度差值 ，相应的外表面及中心点处的温度数值分别是 、；与此同时，墩帽的外表面处和墩帽的套筒部位产生最大的温度差值 ，相应的外表面处和套筒部位的温度数值分别是 、。当承台的外表面处和承台的中心点处在冻融后，产
13、生最大的温度差值 ，相应的外表面及中心点处的温度数值分别是 、。冻融循环作用时，装配式桥墩的表面温度数值随着时间变化，并且与外界温度值近似相等，然而装配式桥墩的内部温度随着冻融循环次数的增加，产生了与外部温度相异的循环特征，并且波动的幅度比桥墩的外表面小。图温度变化曲线 分析结果可知，温度在桥墩内部传递的时候会出现损耗现象，导致中心处的温度数值逐渐减小，并且桥墩的尺寸比较大，造成温度的传导耗时更长，内部的温度值不能够达到环境温度，因此出现了温度差值。约次冻融循环之后，结构逐渐产生了新的平衡状态。墩帽、墩柱和承台均呈现一致的变化趋势，但温度数值存在差异，如承台的中心点处温度的极值波动是最小的，表
14、明温度对大尺寸结构的内部作用最小。在热传导的基础上进行热应力分析，温度变化趋于稳定状态的时候，可得出一次冻融循环历程的温度应力图，见图。由图可知，每个分析步中的温度应力云图，装配式桥墩的外部表面温度应力都比桥墩内部大，表明温度变化对装配式桥墩的表面部位产生最大影响。依据图结果，进一步得出装配式桥墩墩帽的外表面处与墩帽的中心点处、墩柱的外表面处与墩柱的中心点、承台的外表面处与承台的中心点处等个部位温度应力随时间的变化特征，如图所示。由图可知，装配式桥墩的外部表面遭受温度作用的影响比中心点要大，当冻融循环时间进行到时，墩帽、墩柱、承台处的外表面应力与中心点相比，分别增加了、，冻融作用对桥墩外表面处
15、的破坏程度是最大的。对于桥墩的温度应力，其中心点处的变化略滞后外部表面处，一个冻融循环周期内存在个波峰，与温度曲线的变化特征基本上保持一致。结构的尺寸越大，其中心点部位的温度应力数值反而越小，充分表明结构的温度应力由外表面向中心点处递减，可第期包龙生，等：装配式混凝土桥墩节点的抗冻性能研究知装配式桥墩的外部表面处为温度应力造成损伤最大的位置，当桥墩遭受冻融循环的影响时，外表面就会出现较为严重的剥落现象。（）分析步一（）分析步三（）分析步五（）分析步八图温度应力云图 图温度应力曲线图 装配式混凝土桥墩冻融损伤分析 冻融损伤的 子程序开发 有限元软件内部默认的温度场采用节点的温度而不是积分点的温度
16、，故需要通过 子程序使节点的温度值赋予给有限元模型的积分点之中。有限元软件内部默认不广西大学学报（自然科学版）第 卷考虑冻融损伤量的叠加，因此还需要导入 子程序，使软件能够计算装配式桥墩模型在冻融循环了、次时的累积损伤量。本文通过 语言编写相关的程序，如图、所示。图 子程序 图 子程序 子程序内部的 （，）（）表示的是把节点的温度数值赋予至模型的积分点中。根据静水压理论和渗透压假说 ，混凝土的损伤由内部结构自由水和正负温交替循环作用引起的累积损伤，可以通过编写 子程序实现冻融循环下的累积损伤。其中 语句的第一行表示选择温度的变化值低于零时的积分范围，即提取模型由于冻融循环作用造成的损伤阶段，语句的第二行表示单次冻融循环时的总体温度变化，语句的第三行表示单次冻融循环时的损伤程度，其中 表示模型所需的冻融循环作用的次数，依次设置为、次。子程序运行完成时可以判断和获取冻融循环作用下的损伤量，然后叠加循环的次数，最终得出相应冻融循环次数时的累积损伤量。装配式桥墩的冻融损伤分析通过调用个子程序，可以获得模型在温度稳定之后的循环状态，依次对建立的桥墩模型开展、次的冻融循环。需要在模型的底部设置边界