一种考虑锈蚀作用的非完全粘结预应力梁模拟方法

本发明涉及土木工程，特别是涉及一种考虑锈蚀作用的非完全粘结预应力梁模拟方法。

背景技术：

1、预应力混凝土结构因其优良的抗裂性能和承载能力，被广泛应用于桥梁、建筑等工程领域。然而，在长期服役过程中，由于环境因素和材料劣化，预应力钢筋会发生锈蚀，导致结构性能衰减。

2、钢筋锈蚀主要从以下几个方面影响预应力梁的受力性能：

3、（1）有效截面损失：锈蚀产物（铁锈）的体积约为原钢筋体积的2~4倍，产生的膨胀应力导致混凝土保护层开裂剥落，同时钢筋有效截面减小，承载能力下降。

4、（2）粘结性能退化：锈蚀破坏了纵筋（钢筋/预应力筋）与混凝土之间的机械咬合作用，导致粘结强度显著降低，进而影响预应力的传递和梁的协同受力性能。

5、（3）材料性能劣化：锈蚀不仅影响纵筋（钢筋/预应力筋）的屈服强度和延性，也会通过裂缝扩展和混凝土损伤影响整体结构的刚度。

6、传统的预应力梁分析方法多基于完全粘结假设，无法准确反映锈蚀引起的粘结滑移效应。向量式有限元方法虽然能处理大变形问题，但现有研究中尚未系统考虑锈蚀损伤对预应力梁力学行为的影响。

7、专利文献cn121919951a公开了预应力梁桥徐变-锈蚀-水动力疲劳耦合的变形计算与寿命评估方法、设备、存储介质，该方法通过建立含时变徐变本构和预应力筋的梁桥水动力仿真模型，逐年更新混凝土徐变、预应力松弛、预应力筋腐蚀深度及水动力疲劳开裂深度等参数，基于锈蚀与疲劳竞争准则确定预应力筋失效面积，并计算界面粘结性能。随后进行水动力与长期效应耦合仿真，提取跨中挠度、桥墩相对位移及预应力幅值，进而评估预应力筋腐蚀水动力疲劳寿命。

8、专利文献cn122021196a公开一种基于非完全粘结理论与向量式有限元的预应力梁模拟方法，包括：基于非完全粘结理论模型计算预应力梁截面在不同荷载作用下的力学响应，以构建弯矩转角关系；采用向量式有限元方法将预应力梁截面所在预应力梁结构离散为多个质量质点，并以无质量梁单元将多个质量质点进行连接，以构建的全局结构分析模型；将弯矩转角关系作为全局结构分析模型中无质量梁单元的力变形关系，计算全局结构分析模型中各质量质点的位移响应，在计算过程中根据无质量梁单元的变形和弯矩转角关系以更新内力，输出预应力梁结构的模拟结果。

9、因此，建立一种能够同时考虑非完全粘结特性和锈蚀损伤效应的预应力梁模拟方法，对于准确评估既有预应力混凝土结构的剩余寿命和安全状态具有重要意义。

技术实现思路

1、本发明的目的在于克服现有技术的不足，提供一种考虑锈蚀作用的非完全粘结预应力梁模拟方法，能够准确模拟锈蚀对预应力梁力学性能的影响，为既有预应力混凝土结构的性能评估提供有效工具。

2、为实现上述目的，本发明采用如下技术方案：

3、一种考虑锈蚀作用的非完全粘结预应力梁模拟方法，包括以下步骤：

4、步骤s1、根据锈蚀深度，建立考虑锈蚀效应的纵筋-混凝土界面粘结滑移本构模型；

5、步骤s2、根据锈蚀深度，计算锈蚀率、各折减系数、、及开裂损伤因子；

6、基于锈蚀率、截面刚度折减系数、锈蚀粘结强度折减系数、质量折减系数以及开裂损伤因子，计算因锈蚀导致的弹性模量损失、截面刚度损失、粘结强度损失；

7、步骤s3、将预应力梁离散为有限个质点和向量式梁单元，结合步骤s1建立的纵筋-混凝土界面粘结滑移本构模型，以及步骤s2得到的弹性模量损失、截面刚度损失以及粘结强度损失，构建锈蚀预应力梁动力平衡方程；

8、步骤s4、基于当前时刻各质点受到的内力、外力和阻尼力，采用显式时间积分方法，求解锈蚀预应力梁动力平衡方程，以获得当前各质点的位移、速度和加速度响应。

9、具体的，根据锈蚀深度，建立考虑锈蚀效应的纵筋-混凝土界面粘结滑移本构模型，该纵筋包含钢筋或/和预应力筋；

10、该模型中，设为粘结应力，

s为滑移量，粘结应力-滑移本构关系可表达为：

11、；式中，为纵筋（钢筋/预应力筋）锈蚀深度，其与纵筋（钢筋/预应力筋）截面锈蚀率的关系为：

12、；其中，为纵筋（钢筋/预应力筋）初始直径。

13、锈蚀效应引起的粘结强度损失通过引入锈蚀粘结强度折减系数考虑，其表达式为：

14、；式中，和为材料参数。

15、具体的，所述锈预应力梁离散是将锈蚀预应力梁沿长度方向离散为个质点和个向量式梁单元，每个质点具有位置、速度、加速度和质量；锈蚀效应引起的有效质量损失通过质量折减系数考虑，其表达式为：

16、；式中，，和为材料参数。

17、具体的，锈蚀效应引起的截面刚度损失通过截面刚度折减系数考虑，其表达式为：

18、；式中，为初始截面惯性矩，为刚度衰减参数，取值范围为。

19、具体的，锈蚀产物膨胀会在混凝土保护层产生径向压力，当压力超过混凝土抗拉强度时，保护层开裂。锈蚀膨胀压力可表示为：

20、；式中，为混凝土弹性模量，为锈蚀产物体积，为原纵筋（钢筋/预应力筋）体积，为周围混凝土体积；当锈蚀引起的径向应力达到混凝土抗拉强度时，保护层发生开裂，开裂条件为：；式中，为环向应力，为保护层外半径，为纵筋（钢筋/预应力筋）半径；保护层发生开裂后，锈蚀效应引起的弹性模量损失通过引入开裂损伤因子考虑，其表达式为：；式中，为开裂临界锈蚀深度，为损伤演化参数。

21、具体的，所述锈蚀预应力梁动力平衡方程表示为：；式中：为锈蚀单元的内力向量；为外荷载向量；为阻尼力向量；锈蚀单元内力通过考虑锈蚀效应计算得到：；式中， 与  分别为向量式梁单元在当前时间步下的纯拉压变形与纯弯曲变形，表示弹性模量损失，表示截面刚度损失，表示粘结强度损失，采用显式中心差分法，基于质点质量  以及各质点所受的内力、外力和阻尼力数据输入，对锈蚀预应力梁动力平衡方程进行时间积分求解。质点 i 在 t+δt 时刻的位移和速度更新公式为：；；加速度由平衡方程计算：；时间步长需满足稳定性条件：；式中，为系统最小周期，由锈蚀后的截面刚度决定。

22、与现有技术相比，本发明的有益效果：

23、（1）将锈蚀效应与非完全粘结理论相结合，建立了考虑锈蚀的预应力梁模拟方法，能够准确反映不同锈蚀阶段预应力梁的力学行为。

24、（2）通过折减系数定量描述锈蚀对粘结强度、截面刚度、质量和混凝土弹性模量的影响，使模型具有明确的物理意义。

25、（3）基于向量式有限元方法，能够处理锈蚀引起的大变形和几何非线性问题，适用于从弹性阶段到破坏阶段的全过程模拟。

26、（4）算法采用显式时间积分，计算效率高，易于编程实现，便于工程应用。

27、本发明方法可广泛应用于以下领域：

28、（1）既有预应力混凝土桥梁的性能评估：通过检测纵筋（钢筋/预应力筋）锈蚀程度，采用本方法模拟结构实际承载能力，为桥梁养护决策提供依据。

29、（2）预应力混凝土结构的寿命预测：结合环境因素和锈蚀速率，预测结构在不同服役年限后的力学性能，制定合理的维修计划。

30、（3）锈蚀预应力梁加固设计：模拟加固前后结构受力性能，优化加固方案，提高加固效果。