冻损梯度对承压混凝土应力-应变关系的影响

刘西光，雷永洁，周佳丽，牛荻涛，*，刘 俊

（1.西安建筑科技大学 土木工程学院，陕西 西安 710055；2.西安建筑科大工程技术有限公司，陕西 西安 710055）

寒冷地区桥梁会遭受冻融损伤，同时荷载的作用使桥梁结构混凝土处于承载状态［1-3］，此时混凝土的冻融损伤规律与无应力状态下显著不同，荷载裂缝导致混凝土的饱水程度增加、冻融损伤加剧［4］.

国内外学者对冻融损伤混凝土基本力学性能开展了众多的试验研究，发现冻融循环导致混凝土的表面剥落，力学性能下降［5-7］.随着冻融循环次数的增加［8-10］，混凝土的峰值应力逐渐降低，峰值应变增大，脆性破坏特征更加明显［11-13］.许多学者对冻融损伤混凝土本构模型进行了研究：龙广成等［14］基于应变等效假设和统计损伤理论，建立了一种冻融损伤混凝土单调受压本构模型；Duan 等［15-16］建立了冻融循环作用下混凝土的随机损伤本构模型；关虓等［17］基于Weibull 强度理论，建立了冻融循环作用下混凝土的细观统计损伤本构模型；Qiu 等［18-19］基于塑性损伤理论，建立了冻融循环作用下混凝土的塑性损伤本构模型.

冻融后混凝土可分为损伤层、损伤过渡层和未损伤层［20］.已有冻融损伤混凝土力学性能研究均未考虑损伤梯度的影响，应力作用下混凝土冻融损伤沿试件深度方向的影响规律尚不清楚.因此，需要开展不同冻融损伤程度下承压混凝土的应力-应变关系研究.

1 试验

1.1 试件设计和材料

为得到不同冻融损伤程度的混凝土试件，设计了不同应力水平作用下的混凝土圆柱体试件，试件尺寸均为φ200×1 000 mm（见图1）.构件所用混凝土均为C40强度等级的商品混凝土，配合比见表1.其中，水为普通自来水，水泥为普通硅酸盐水泥，细骨料为天然河砂，粗骨料为5～20 mm 连续级配石子.预应力钢筋选用PSB830 级φT18 的精轧螺纹钢筋，屈服强度为918 MPa，抗拉强度为1 087 MPa，伸长率为9.1%.

图1 试件尺寸及张拉装置Fig.1 Specimen size and tensioning device（size：mm）

冻融循环次数（N）分别为0、100、200、250次；混凝土的应力水平（μ）为0、0.2、0.3和0.4，即张拉应力为混凝土轴心抗压强度平均值的0%、20%、30%和40%.

1.2 冻融循环试验

张拉试验开始前对试件进行为期4 d 的浸泡，采用自行设计加工的反力架进行张拉，张拉装置如图1所示.冻融循环试验采用ZHT/W2300 型气候环境模拟实验室进行气冻气融循环试验，温度范围为（-19～25）℃，升降温速率均为0.7～1.0 ℃/min.

1.3 分层取芯

考虑冻融损伤由表及里发展与不同深度处混凝土力学性能的非均匀退化，对冻融后的承压混凝土圆柱体试件进行分层取芯，得到圆环体试件，以研究不同冻融损伤程度下承压混凝土力学性能的退化规律.

冻融结束后对试件进行切割.为了消除端部复杂应力的影响，每个混凝土圆柱体试件沿长度方向两端各截掉200 mm，形成3 个φ200×200 mm 的圆柱试件，分4层进行取芯，见图2.用x表示试件外边缘到每层试件中心的距离，r表示圆柱体试件的半径，定义试件相对冻融深度λ=x/r，得λ=0.1、0.3、0.5、0.8.

2 结果与讨论

2.1 冻融损伤混凝土的表观形态

图3 为典型混凝土圆柱体构件表面裂缝的开展情况.由图3 可见，在冻融循环作用下，试件表面出现了沿长度方向发展的微裂缝，且随着冻融循环次数的增加，裂缝的数量和宽度也在逐步增加，表明混凝土的冻融损伤逐渐加剧.

图3 典型混凝土圆柱体构件表面裂缝的开展情况Fig.3 Development of surface cracks of typical concrete cylinder members

2.2 冻融损伤混凝土轴压试验结果分析

2.2.1 混凝土单轴受压的破坏形态

图4 为冻融损伤混凝土的单调受压破坏形态.由图4 可见：混凝土试件在单调压缩加载过程中，多条纵短裂缝从试件两端同时出现，部分受压裂缝有斜向开展的趋势；冻融循环次数越多，裂缝数量越少、发展越集中，呈现的脆性破坏特征越明显.

图4 冻融损伤混凝土的单调受压破坏形态Fig.4 Monotonic compression failure modes of freeze-thaw damaged concrete（λ=0.1）

2.2.2 应力-应变曲线

图5 为冻融损伤承压混凝土圆环体试件的单调受压应力-应变（σ-ε）曲线和理论结果.由图5 可见：

（1）当相对冻融深度小于0.3 时，随着冻融循环次数和应力水平的增加，混凝土的应力-应变曲线逐渐趋于平缓，整体右移.当相对冻融深度为0.1，冻融循环次数为0～250 次时，混凝土的峰值应力降低了32%～77%，峰值应变增加了17%～110%.

（2）当相对冻融深度分别为0.3 和0.5 时，应力-应变曲线在100 次冻融循环时发生了明显的变化.当冻融循环次数增加到200 次时，2 种相对冻融深度下混凝土应力-应变曲线的发展规律相似.

（3）当相对冻融深度为0.8 时，不同工况下混凝土的应力-应变曲线与未冻融时无明显差异.这表明混凝土的浅层损伤大于深层，相对冻融深度超过0.8时混凝土几乎不受冻融影响.

2.2.3 相对峰值应力和相对峰值应变

图6 为冻融损伤混凝土相对峰值应力的退化规律.图7 为冻融损伤混凝土相对峰值应变的退化规律.由图6、7 可见：

图7 冻融损伤混凝土相对峰值应变的退化规律Fig.7 Degradation of relative peak stain of freeze-thaw damaged concrete

（1）当相对冻融深度小于0.3时，随着冻融循环次数的增加，混凝土的峰值应力逐渐降低而峰值应变增加；当相对冻融深度增加到0.5 时，100 次冻融循环前混凝土的峰值应力较未冻融时相差不大，超过100 次循环后的峰值应力显著降低；当相对冻融深度增加到0.8 时，不同冻融循环次数下混凝土的重复受压峰值应力相差不大，表明该层不受冻融循环的影响.当相对冻融深度为0.1，应力水平为0时，冻融100、200、250次试件的峰值应力分别降低了23%、40%和50%.

（2）在不同的冻融循环次数下，混凝土的相对冻融深度存在临界值.当冻融循环次数为100 次时，相对冻融深度临界值为0.5；当冻融循环次数增加到200 次时，相对冻融深度临界值为0.8.当相对冻融深度小于临界值时，混凝土的峰值应力和峰值应变呈线性退化，超过临界值后保持稳定.

3 考虑冻融损伤梯度的承压混凝土应力-应变关系

3.1 数学模型

基于Lemaitre 应变等价性假说［21］，冻融循环作用下承压混凝土的受压本构关系可以表示为式（1）.其中，冻融后混凝土初始损伤值可按式（2）进行计算.

式中：σn为经历n次冻融后的混凝土应力，MPa；E0为未冻融混凝土的弹性模量，MPa；Dm为混凝土的初始损伤值；En为经历n次冻融后混凝土的弹性模量，MPa.

假设当混凝土受到荷载作用后，其微元强度符合Weibull 概率分布［22］，则荷载损伤微元（Dc）可用破坏微元（Nd）与总微元（Nt）的比值表示：

式中：a和b分别为Weibull分布参数.

根据应变等价原理，可以得到荷载作用下混凝土的本构关系：

结合式（1）、（3），得到混凝土总损伤变量（D）：

将式（5）代入式（4），可以得到冻融损伤承压混凝土的单调受压应力-应变关系：

3.2 应力-应变关系的验证

通过式（6）对试验应力-应变曲线进行回归拟合，得到不同相对冻融深度下a、b的取值如图8、9 所示.由图8、9 可见：当相对冻融深度小于0.5 时，随着冻融循环次数的增加，参数a和b呈现递增的趋势；当相对冻融深度为0.8 时，参数a和b随冻融循环次数与应力比的变化较小.由此可以看出，随着冻融循环次数的增加，混凝土的峰值应变增加，应力-应变曲线下降段变陡，脆性特征更加明显.随着应力比的增加，混凝土的峰值应力不断降低，表明应力作用加剧了混凝土的损伤.并且冻融损伤表现出逐层递减的趋势，当相对冻融深度为0.8 时，冻融损伤不再对混凝土应力-应变曲线的形状造成影响.

图8 参数a 随冻融循环次数的变化Fig.8 Variations of parameter a with freeze-thaw cycles

图9 参数b 随冻融循环次数的变化Fig.9 Variations of parameter b with freeze-thaw cycles

将不同深度处参数a和b的取值代入本构方程中，得到理论计算结果如图5 所示.由图5 可见，该理论模型计算结果与试验结果吻合较好.

4 结论

（1）冻融损伤混凝土的浅层损伤大于深层.当相对冻融深度为0.1 时，混凝土的应力-应变曲线随着冻融次数的增加逐渐趋于扁平，整体右移.当相对冻融深度分别为0.3 和0.5 时，应力-应变曲线在100 次冻融循环时发生了明显的变化.当相对冻融深度为0.8 时，各冻融循环次数下的应力-应变曲线与未冻融时基本吻合，表明混凝土未有较大的损伤劣化行为.

（2）不同冻融循环次数对应不同的相对冻融深度临界值.当冻融循环次数为100 次时，相对冻融深度临界值为0.5；当冻融循环次数增加到200 次时，相对冻融深度临界值增加到0.8.当相对冻融深度小于临界值时，混凝土的峰值应力、峰值应变呈线性退化，超过临界值后保持稳定.

（3）基于应变等价性假说和统计损伤理论，建立了不同冻融损伤程度下承压混凝土的受压应力-应变关系模型.