王燕 徐善华 李安邦

摘 要

为研究寒冷及地震设防地区既有钢筋混凝土结构的抗震性能，对经历冻融循环作用的 8 根钢筋混凝土柱进行低周反复荷载试验研究。综合考虑冻融循环次数与轴压比对试件柱承载力和滞回性能的影响规律，通过理论分析与试验回归相结合的方法，提出了冻融损伤钢筋混凝土柱的三折线刚度退化骨架曲线，给出了试件柱屈服前后卸载刚度的计算式及具体的滞回规则，进而建立了冻融损伤钢筋混凝土柱“指向定点”的刚度退化三线型恢复力模型。通过试验曲线与计算曲线的对比分析表明:建议的恢复力模型预测结果较好，可为冻融环境下钢筋混凝土柱的地震反应分析提供理论依据。

引 言

在寒冷及地震多发地区，冻融循环的反复作用，导致混凝土材料的力学性能降低、混凝土与钢筋界面黏结性能退化，造成混凝土结构构件的几何截面损失、有效承载强度下降，严重威胁结构的安全运营 [1-3] 。地震作用下，受冻融损伤的钢筋混凝土结构或构件，其损伤积累效应继续演化，必然导致其滞回性能、耗能能力、延性、刚度等抗震性能退化，甚至引发结构整体破坏，从而造成不可挽回的人员伤亡和经济损失。因此，大力开展冻融环境下既有钢筋混凝土结构抗震性能的试验研究和理论分析，揭示冻融损伤混凝土结构的地震损伤累积行为及其抗震性能演变规律，不但符合客观实际，而且是寒冷地区混凝土结构安全评定与持久运营的迫切需求。对进一步完善钢筋混凝土结构抗震性能评估的理论和方法，对混凝土结构工程防灾减灾具有重大的理论意义和工程价值。

近年来，众多研究者在冻融损伤后结构或构件的整体性能层面展开了一系列的研究，在静力试验和理论方面积累了较多的研究成果 [4-7] ，在抗震性能方面的研究也已逐步开展。文献[8-9]通过 6 个不同轴压比的钢筋混凝土(RC)框架中节点梁柱组合体试件的冻融循环试验及低周反复加载试验发现:随着冻融循环次数的增加，贯穿节点的梁端纵筋与混凝土之间的黏结强度显著退化，导致低周反复荷载作用下节点核心区破坏模式发生转变，承载力降低，梁柱组合体耗能能力下降;基于试验结果，采用回归分析的方法建立了适应于冻融损伤的钢筋混凝土框架节点三折线型骨架曲线模型及相应的恢复力模型。文献[10-12]在对 8 片经历冻融循环的低矮 RC 剪力墙拟静力试验研究的基础上，综合考虑冻融循环次数、混凝土抗压强度和轴压比对剪力墙承载力和延性的影响，考虑累积损伤效应造成的承载力和刚度退化，引入第 i 次循环加载时的退化指数 β i ，提出了适用于遭受冻融循环的低矮钢筋混凝土剪力墙的滞回规则，进而建立了考虑冻融循环损伤作用的低矮钢筋混凝土剪力墙恢复力模型。上述研究均对冻融损伤后钢筋混凝土结构或构件的抗震演化规律进行了积极的探索。

然而，对于钢筋混凝土结构或构件而言，一方面，其冻融循环后的力学性能是混凝土强度、钢筋与混凝土黏结性能、混凝土保护层厚度、尺寸等多因素劣化后综合影响的结果;另一方面，现有的研究中由于冻融循环试验条件的差异，在破坏机理的认识上还存在一些争议 [13-15] ，试验研究及相应的理论模型仍然比较匮乏。鉴于此，本文拟通过 8 个钢筋混凝土柱的冻融循环试验及低周反复荷载试验，用理论推导和试验回归相结合的分析方法，建立冻融损伤钢筋混凝土柱的恢复力模型。

试验研究

1.1冻融循环试验

为了研究冻融损伤钢筋混凝土柱的抗震性能并提出相应的恢复力模型，本研究共设计了 8 根钢筋混凝土柱，其几何尺寸和配筋如图 1 和表 1 所示。柱子采用平卧方法浇筑，在自然环境下养护 28 d 后，放入特制的不锈钢水槽中充分浸泡 4 d，然后置于 ZHT/W2300 型人工气候模拟实验室进行冻融循环试验。考虑到寒冷地区的大多数混凝土工程，均因暴露在雨淋或蒸汽等一般大气环境下发生冻融损伤的实际情况，采用“气冻气融”的冻融循环方式。在冻融循环试验结束后，再给钢筋混凝土柱浇筑刚性底座。

图 1 柱的尺寸及配筋

Fig.1 Sizes and of

1.2低周反复荷载试验

低周反复荷载试验采用悬臂柱式加载方式。试件柱顶部用量程为 500 kN 的千斤顶施加竖向轴向压力;试件柱离基座 1 000 mm 高度处用 500 kN 的MTS 电液伺服作动器施加水平反复荷载，试验全过程由 MTS 电液伺服微机系统控制，由 TDS－602 数据采集仪采集数据。冻融循环试验及低周反复荷载试验如图 2 所示。

a—冻融循环试验; b—低周反复荷载试验。

图 2 冻融循环及低周反复荷载试验

Fig.2 -thaw and quasi- tests

1.3低周反复荷载试验结果分析

各试件柱的试验滞回曲线如图 3 所示。由图可知，不同设计参数的试件，其试验滞回曲线有以下共同特征:试件屈服前，荷载－位移关系近似呈线性关系，滞回圈面积很小，加卸载路径基本重合;随着荷载的增加，荷载－位移曲线逐渐偏离直线，卸载时有一定的残余变形，滞回圈面积逐渐变大，荷载逐渐达到峰值;峰值荷载过后，随着位移的继续增加，荷载不断下降，残余变形不断变大，直至试件破坏。

a—N1(0); b—N1(100); c—N1(200); d—N1(300); e—N2(0); f—N2(100); g—N2(200); h—N2(300)。

试验曲线; ---计算曲线。

图 3 8 根试验柱试验滞回曲线与计算滞回曲线对比

Fig.3 test and of 8 test

不同设计参数对试验滞回曲线的影响规律如下:

1) 冻融循环作用的影响。在轴压比为 0.2 的 4 个试件中，随着冻融循环次数的依次增加，滞回曲线由稳定丰满的“梭形”，逐渐趋于具有“捏缩效应”的“弓形”。当冻融循环次数达到 300 次时，试件 N1(300)的滞回圈数剧减，峰值荷载后滞回曲线稳定性显著下降，破坏较突然。这是混凝土强度及钢筋与混凝土间黏结性能遭冻融循环作用破坏的结果。

2) 轴压比的影响。在相同的冻融循环次数下，如 N1(200) 和 N2(200)，轴压比分别为 0.2 和 0.32，N1(200)的峰值荷载明显低于 N2(200)的峰值荷载;但 N1(200)的滞回曲线更加饱满，峰值荷载过后的滞回曲线更加稳定，破坏前的滞回圈数量更多。这说明轴压比小的柱虽然承载力较低，但其吸收或消耗地震能量的能力更强，抗震性能更好。

冻融损伤钢筋混凝土柱恢复力模型

分析 8 个钢筋混凝土试件柱的骨架曲线及相关试验数据发现，8 个试验柱均以弯曲破坏为主。由于箍筋的约束作用，试件柱在出现裂缝、纵筋屈服时刚度并没有显著降低，因此将骨架曲线屈服点以前简化为刚度没有变化的弹性阶段;屈服点到峰值点之间为骨架曲线的强化阶段;峰值点到破坏点这段，骨架曲线为明显的负刚度段，即强度下降段。因此将未冻融钢筋混凝土柱骨架曲线简化为特征点为屈服点、峰值点和破坏点的刚度退化的三折线模型 [16] ，同时假定冻融损伤钢筋混凝土柱的骨架曲线与未冻融构件的骨架曲线相似。则三线型的骨架曲线如图 4 所示。

图 4 三线型恢复力骨架曲线

Fig.4 force curve

2.1未冻融钢筋混凝土柱骨架曲线特征点

对于未冻融钢筋混凝土柱，通过对其正截面的受力全过程进行分析 [17] ，可得 M - ϕ 关系曲线上的三个特征曲率:屈服曲率 ϕ y 、峰值曲率 ϕ m 及破坏曲率 ϕ u 。

对于形成于柱的端部区域的塑性铰区，其长度 l p 采用考虑轴压力影响 [18-19] 的计算式:

式中: f y 、 f ' y 分别为受拉钢筋、受压钢筋的屈服强度; μ s 、μ' s 分别为受拉钢筋、受压钢筋的配筋率; N 为轴压力; f c 为混凝土的抗压强度; b 、 h 、 h 0 分别为柱截面的宽度、高度、有效高度。

假定柱截面曲率沿高度呈线性分布，则柱顶端特征水平位移 Δ 与特征曲率 ϕ 的计算式为:

Δ y0 = ϕ y L 2 /3          (2a)

Δ m0 = ϕ y L 2 /3+( ϕ m − ϕ y ) l p ( L − l p )     (2b)

Δ u0 = ϕ y L 2 /3+( ϕ u − ϕ y ) l p ( L − l p )      (2c)

式中: Δ y0 、 Δ m0 、 Δ u0 分别为未冻融钢筋混凝土柱的屈服位移、峰值位移、破坏位移; L 为构件柱的净高。

考虑轴压力引起的 P － Δ 二阶效应及弯矩平衡，未冻融试件柱骨架曲线各特征点水平荷载的计算式为:

P y0 =( M y0 − NΔ y0 )/ L           (3a)

P m0 =( M m0 − NΔ m0 )/ L           (3b)

P u0 =( M u0 − NΔ u0 )/ L           (3c)

2.2冻融损伤钢筋混凝土柱骨架曲线特征值

回归分析试件柱骨架曲线上各相对特征点参数与冻融循环作用次数关系的试验结果，得到冻融循环后试件柱骨架曲线特征点的荷载和位移计算式:

1) 屈服荷载 P y 与屈服位移 Δ y 。

P y =(−10 −5 n 2 +0.0031 n +1) P y 0       (4a)

Δ y =(−0.0011 n +1) Δ y 0           (4b)

2) 峰值荷载 P m 与峰值位移 Δ m 。

P m =(−10 −5 n 2 +0.0035 n +1) P m 0      (5a)

Δ m =(−0.0012 n +1) Δ m 0          (5b)

3) 破坏荷载 P u 与破坏位移 Δ u 。

P u =(−10 −5 n 2 +0.0034 n +1) P u 0       (6a)

Δ u =(−0.0013 n +1) Δ u 0          (6b)

2.3刚度退化规则

1) 屈服前的卸载刚度计算。

分析本次低周反复荷载试验的 8 个试件柱的骨架曲线发现:钢筋混凝土柱在屈服之前，骨架曲线上并没有明显的转折点，而且卸载后残余变形很小，因此可以认为屈服前试件柱仍处于弹性阶段。于是，屈服前的卸载刚度 K r1 与弹性刚度 K e 相同，计算式为:

K r1 = K e = P y / Δ y            (7)

2) 屈服后的卸载刚度计算。

分析 8 个试件柱低周反复荷载试验的滞回曲线发现，试件柱在屈服后的卸载刚度与同级位移下的循环次数无关，而随着位移幅值的增加而逐渐减小。对 8 个试件柱屈服后不同位移幅值下的卸载刚度进行归一化处理，见图 5。则其屈服后在某一级位移下卸载刚度 K r2 的计算式为:

K r2 = K e ( Δ m / Δ y ) − α           (8)

式中: Δ m 为屈服后某一级的位移幅值; α 为刚度退化指数，本文取 0.667。

图 5 卸载刚度的归一化处理

Fig.5 of

2.4滞回规则

通过对 8 个试件柱滞回特性的分析及计算，建议恢复力模型的滞回规则见图 6，具体描述如下:

图 6 滞回规则

Fig.6 rules

1) 正反向的骨架曲线对称，均为简化的三折线模型;特征点分别为屈服点、峰值点以及破坏点。

2) 在试件柱临近屈服或屈服后，正、反向加载均指向相应某一“定点”;分析试验数据，得出正反向定点即为骨架曲线上荷载为±0.35 P y 的点。

3) 试件柱屈服前，加载和卸载均沿骨架曲线的弹性段进行:若加载至正向屈服点 1 后卸载，则沿加载路径返回，卸载刚度为 K r1 ，继续反向加载指向反向屈服点 2。

4) 试件柱屈服但受力未到峰值荷载时，正向加载沿骨架曲线 1—3—9 进行;若加载到 3 点开始卸载并到零荷载点 4，此时的卸载刚度为 K r2 ; 反向加载则指向反向定点 5，直到最大位移点6;然后反向开始卸载到零荷载点 7，继续正向加载并指向正向定点 8，而后到达峰值点 9。整个正向卸载、反向加载、反向卸载、再加载路径为 3—4—5—6—7—8—9。

5) 峰值荷载过后，加载沿着骨架曲线 9—11—15—16 进行;若加载到 11 点开始卸载到零荷载点12，则此时卸载刚度仍为 K r2 ; 反向加载仍指向反向定点 5 并到最大点 13;反向卸载到零荷载点 14，并继续正向加载指向正向定点 8，而后达到最大位移点 15。整个正向卸载、反向加载、反向卸载，再加载路径为 11—12—5—13—14—8—15。如此依次进行，直到试件破坏为止。

模型验证

3.1骨架曲线验证

利用上述算式计算得到的骨架曲线特征点的荷载及位移，与试验所得的骨架曲线特征点的荷载及位移对比，见表 2 及表 3。

由表 2、表 3 可以看出，大多数计算值与试验值的相对误差都在 10%以内，满足精度要求。个别计算值的相对误差较大，如试件 N2(300) 的屈服荷载、破坏荷载的计算相对误差分别达到 -19%、-17%，试件 N2(0)的峰值位移的计算相对误差达到 17%，试件 N2(200)的屈服位移的计算相对误差达到 18%。产生的原因可能是:1) 试件柱在试验过程中随着荷载的增加，混凝土开裂或纵筋屈服后截面各纤维的应变已不符合平截面，且在计算中未考虑冻融作用对钢筋与混凝土之间黏结的劣化作用。2) 本文试验拟合的数据量小，为小样本，也会造成一定的误差。

3.2滞回曲线验证

将上述骨架曲线计算模型结合相应的刚度退化规律和滞回规则，可得到冻融损伤后钢筋混凝土柱的恢复力特性曲线。两种轴压比下的试验滞回曲线与计算滞回曲线对比见图 3。

从图 3 可以看到，对于两种轴压比下经历 0、100、200 冻融循环作用的试件柱而言，计算滞回曲线能够较好地模拟试验滞回曲线。而对于经历 300次冻融循环作用的试件柱而言，可能是未能进一步考虑钢筋与混凝土之间剪切滑移引起的捏拢效应，计算滞回曲线模拟精度较差。

结束语

在冻融循环后钢筋混凝土柱低周反复荷载试验研究的基础上，通过理论分析和试验回归分析相结合的方法，提出了适用于冻融损伤钢筋混凝土柱的恢复力模型，得到以下结论:

1) 以弯曲破坏为特征的冻融损伤钢筋混凝土柱的骨架曲线可简化为三折线。未冻融试件柱的各特征点根据平截面假定推导计算得出;冻融后试件柱各特征点通过试验数据回归分析得到。各特征点的计算值与试验值的相对误差满足工程精度的要求。

2) 所提的恢复力模型为“定点指向”型，正反向定点即为骨架曲线上荷载为±0.35 P y 的点。模型考虑了冻融钢筋混凝土柱屈服后卸载刚度的退化，给出了相应的滞回规则，能较准确地描述冻融后的压弯构件在低周反复荷载下的滞回特性。

3) 从计算曲线与试验曲线的对比看到，所提恢复力模型有一定的工程精度且简单实用，其计算结果可进一步为冻融损伤后钢筋混凝土压弯构件的非线性地震反应分析提供参考。