一种高原环境考虑防护与修复效果的混凝土寿命计算方法

1.本发明属于混凝土服役寿命计算领域，特别是涉及一种高原环境考虑防护与修复效果的混凝土寿命计算方法。


背景技术：

2.我国川藏铁路、青藏高速公路等国家重大交通工程正在往高原地区推进。该类工程区域具有最大28℃昼夜温差、30％-50％相对湿度、最低气温-32℃低温严寒、180天以上正负温交替高频冻融、盐浓度超标准上限数倍超盐渍土等极端环境特征。大温差、极端干燥、非饱和冻融与盐渍土的耦合作用，导致混凝土性能快速劣化、施工运维困难，对桥隧为主体的重大交通工程安全服役与耐久形成严峻挑战。表层防护、基体损伤修复与钢筋锈蚀修复是保障高频冻融、大温差、超盐渍土等高原环境下钢筋混凝土耐久性的重要途径。
3.目前，国内外学者开展了较多冻融、盐侵蚀、干湿交替环境下混凝土腐蚀劣化实验，收集氯离子、硫酸盐等离子浓度分布，从而预测混凝土服役寿命。然而，高原环境复杂，防护与修复材料种类多、性能参差不齐，现有研究方法对于单一高原环境下混凝土寿命计算不准确，无法适用于高原多因素耦合环境下的混凝土寿命计算，更无法计算考虑防护与修复作用的混凝土服役寿命。


技术实现要素：

4.本发明的目的是提供一种高原环境考虑防护与修复效果的混凝土寿命计算方法，以解决现有方法无法计算高原环境下考虑防护与修复效果的混凝土服役寿命问题。
5.为实现上述目的，本发明提供了一种高原环境考虑防护与修复效果的混凝土寿命计算方法，包括以下步骤：
6.s1.确定模拟高原环境的温度、湿度、盐浓度和内、外力荷载作用情况，建立高原环境、荷载参数与时间的映射关系；
7.s2.构建高原环境下考虑防护与修复作用的混凝土几何模型，并确定所述几何模型中各物相的材料属性；
8.s3.建立温度-湿度-盐浓度-荷载多因素耦合作用物理场方程，在所述混凝土几何模型的基础上设定内部、外部边界条件和初始条件，其中，所述外部边界条件为所述高原环境、荷载参数与时间的映射关系；
9.s4.求解高原环境和荷载作用下考虑防护与修复作用的混凝土内部温度、湿度、盐浓度、应力分布及其随时间的变化情况；
10.s5.计算考虑防护与修复作用的混凝土内部的钢筋达到临界锈蚀浓度发生脱钝状态的时间，确定混凝土服役寿命。
11.可选的，所述s1中，所述高原环境、荷载参数与时间的映射关系，为固定不变或任意形式变化的温度、湿度、盐浓度和内、外力荷载；盐为侵蚀介质，包括：氯化钠、氯化钙、氯化铯、硫酸钠、硫酸镁；内、外力荷载包括：重力、内部孔隙膨胀力、外部荷载。
12.可选的，所述s2中，所述混凝土几何模型的构建过程为：采用蒙特卡洛方法随机生成混凝土细观结构，在混凝土表层添加任意厚度的防护涂层几何模型，所述防护涂层几何模型在1维结构中为与混凝土相连的线段，在2维模型中为与混凝土贴合的四面体，在3维结构模型中为与混凝土相贴合的六面体，损伤修复材料是均匀分布在水泥浆体内部的点阵或空间域；
13.可选的，所述s2中，所述混凝土几何模型包含粗骨料、细骨料和水泥浆体，粗细骨料级配类型可自由设置，骨料体积分数在0％～80％之间，混凝土表面防护涂层是粘附在混凝土表面的连续涂层，混凝土修复材料是以场域形式均匀分布在水泥浆体内部，产生损伤后自动触发修复作用，降低侵蚀介质传输速率并恢复混凝土力学性能。
14.可选的，所述s2中，所述物相包括：骨料、浆体、涂层；
15.所述浆体的材料属性包括：弹性模量、泊松比、密度、抗拉强度、断裂能、相变温度、相变温度间隔、相变潜热、导热系数、恒压热容、离子扩散系数、孔隙率、孔径分布；
16.所述涂层的材料属性包括：导热系数、恒压热容、离子扩散系数；
17.所述骨料和所述涂层的相关实验参数通过实验测试获得；所述浆体由砂和水泥浆体组成，浆体的参数由细观力学等效介质理论计算得到。
18.可选的，所述s3中，所述多因素耦合作用物理场方程包括：温度传导方程、水分传输方程、盐传输方程和力作用方程；
19.可选的，所述温度传导方程为：
[0020][0021]
其中，ρ为密度，c为比热容，t为温度，t为时间，λ为导热系数，l为相变潜热，θ
l
为含水率；内部初始温度设定为与高原环境初始温度保持一致，边界条件选用dirichlet边界条件，与s1中确定的高原环境温度参数保持一致；
[0022]
所述水分传输方程考虑了气态水和液态水的两项传输作用，表达如下：
[0023][0024]
其中，pc为水分特征曲线，k
sl
为饱和状态下水在材料中的渗透系数，为孔隙率，η为水的粘度，k
rl
为液相的相对渗透系数，mw为水的摩尔质量，r为相对气体常数，ρ
l
为水的密度，d
g0
为水蒸气在空气中的自由扩散系数；p
gs
为饱和蒸气压；
[0025]
内部初始水分含量根据实际混凝土设定，范围在0～100rh％，边界条件选用dirichlet边界条件，与s1中确定的高原环境水分参数保持一致；
[0026]
盐离子传输仅考虑对流和扩散两个作用，不考虑化学反应时，所述盐传输方程如下：
[0027][0028]
其中，c为盐离子浓度，d0为浆体中的基准离子扩散系数，t0为基准离子扩散系数测定的时间，m为时间依赖系数，θc为孔隙临界饱和度，u为混凝土中的离子活化能，t
ref
为基准离子扩散系数测定时的温度，d为损伤因子，由损伤力学能量法求得，fr为修复材料对损伤的修复率，df为离子在水中的扩散系数，rb为混凝土对离子的结合影响因子；内部初始离子
浓度设定为0，边界条件选用dirichlet边界条件，与s1中确定的高原环境离子参数保持一致；
[0029]
所述力作用方程如下：
[0030][0031]
其中，σ为应力，b为biot系数，p
*
为环境作用等效荷载，所述环境作用等效荷载包括：冻融、硫酸盐因素引起的内膨胀等效荷载；
[0032][0033][0034][0035][0036]
其中，p
ft*
为冻融引起的内膨胀荷载力，r
max
和r
min
为基体内最大和最小孔径，k、k
l
、kc分别为浆体、水和冰的基体模量，∈为浆体的应变，b
l
、bc分别为水和冰的biot系数，αs为浆体的热膨胀系数，p
atm
为高原环境气压，p
atm0
为浇筑浆体时的环境气压，sc、s
l
、sg分别为对应温度下冰、水和气孔的体积分数，p
s*
为硫酸盐侵蚀引起的内膨胀等效荷载，v
aft
为膨胀钙矾石体积分数，ωi为对应环境作用影响因子，p
i*
分别代表冻融、硫酸盐、镁盐环境等效荷载。
[0037]
可选的，所述s5具体为：以钢筋开始锈蚀时的寿命为混凝土的最终服役寿命，通过步骤s1-s4求解不同深度处的离子浓度，对比临界浓度确定混凝土服役寿命，钢筋锈蚀修复材料会提高钢筋的临界离子浓度，因此混凝土服役寿命计算如下：
[0038][0039]
其中，tc为混凝土服役寿命，xd为钢筋埋置深度，f
rr
为钢筋锈蚀修复材料作用因子，ccr为钢筋固有的锈蚀临界浓度。
[0040]
本发明的技术效果为：
[0041]
本发明可以模拟不同高原环境下使用单一的表面防护涂层、裂缝修复材料和钢筋阻锈材料或同时使用几类防护与修复材料后混凝土内部温度、湿度、应力响应和服役寿命，从而为高原环境下混凝土结构设计与施工提供依据。
附图说明
[0042]
构成本技术的一部分的附图用来提供对本技术的进一步理解，本技术的示意性实施例及其说明用于解释本技术，并不构成对本技术的不当限定。在附图中：
[0043]
图1为本发明实施例中的为2维形式考虑防护与修复材料的几何模型示意图；
[0044]
图2为本发明实施例中的使用不同类型防护涂层后混凝土表面0-5mm处侵蚀介质
浓度分布情况示意图；
[0045]
图3为本发明实施例中的使用不同类型防护涂层后钢筋埋置深度处侵蚀介质浓度分布情况示意图；
[0046]
图4为本发明实施例中的使用不同厚度的防护涂层后混凝土内部的侵蚀介质浓度分布情况示意图；
[0047]
图5为本发明实施例中的使用不同寿命的防护涂层后混凝土内部的侵蚀介质浓度分布情况示意图；
[0048]
图6为本发明实施例中的使用修复率为50％的裂缝修复材料后混凝土内部的侵蚀介质浓度分布情况示意图；
[0049]
图7为本发明实施例中的使用修复率为70％的裂缝修复材料后混凝土内部的侵蚀介质浓度分布情况示意图；
[0050]
图8为本发明实施例中的使用修复率为90％的裂缝修复材料后混凝土内部的侵蚀介质浓度分布情况示意图；
[0051]
图9为本发明实施例中的使用防护涂层后钢筋深度处的混凝土侵蚀介质浓度和临界浓度比较情况示意图；
[0052]
图10为本发明实施例中的使用损伤修复材料和钢筋锈蚀修复后钢筋深度处的混凝土侵蚀介质浓度和临界浓度比较情况示意图。
具体实施方式
[0053]
需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本技术。
[0054]
本发明的目的是提供一种高原环境考虑防护与修复效果的混凝土寿命计算方法，解决现有方法无法计算高原环境下考虑防护与修复效果的混凝土服役寿命问题。
[0055]
一种高原环境考虑防护与修复效果的混凝土寿命计算方法，包括以下步骤：
[0056]
(1)确定模拟高原环境的温度、湿度、盐浓度和内、外力荷载作用情况，建立上述环境和荷载参数与时间的映射关系；
[0057]
(2)构建考虑防护与修复作用的混凝土细观结构几何模型，并确定几何模型中各物相的材料属性；
[0058]
(3)建立温度-湿度-盐浓度-荷载多因素耦合作用物理场方程，在步骤(2)建立的混凝土几何模型基础上设定内部、外部边界条件和初始条件，其中外部边界条件为步骤(1)中确定的高原环境和荷载参数与时间的映射关系；
[0059]
(4)求解高原环境和荷载作用下考虑防护与修复作用的混凝土内部温度、湿度、盐浓度、应力分布及其随时间的变化情况；
[0060]
(5)计算考虑防护与修复作用的混凝土内部的钢筋达到临界锈蚀浓度发生脱钝状态的时间，确定混凝土服役寿命。
[0061]
其中，所述步骤(1)中环境和荷载参数与时间的映射关系，可以是固定不变的或任意形式变化的温度、湿度和外荷载作用。
[0062]
所述步骤(2)中考虑防护与修复作用的混凝土几何模型，采用蒙特卡洛方法随机生成混凝土细观结构，在混凝土表层添加任意厚度的防护涂层几何模型，防护涂层几何模
型在1维结构中为与混凝土相连的线段，在2维模型中为与混凝土贴合的四面体，在3维结构模型中为与混凝土相贴合的6面体。损伤修复材料是均匀分布在水泥浆体内部的点阵或空间域。
[0063]
所述步骤(2)的考虑防护与修复作用的混凝土几何模型中存在骨料、浆体、涂层等物相，浆体需提供包括弹性模量、泊松比、密度、抗拉强度、断裂能、相变温度、相变温度间隔、相变潜热、导热系数、恒压热容、离子扩散系数、孔隙率、孔径分布等参数，涂层需提供导热系数、恒压热容、离子扩散系数。
[0064]
其中，骨料和涂层的相关实验参数通过实验测试获得。由于浆体由砂和水泥浆体组成，因此是一种复合材料，浆体的参数由细观力学等效介质理论计算得到。
[0065]
所述步骤(3)中的考虑防护与修复作用时，混凝土中的多因素耦合作用物理场方程包括温度传导方程、水分传输方程、盐传输方程和力作用方程。其中温度传导方程如下：
[0066][0067]
式中具体的ρ为密度，c为比热容，t为温度，t为时间，λ为导热系数，l为相变潜热，θ
l
为含水率。
[0068]
内部初始温度设定为与高原环境初始温度保持一致，边界条件选用dirichlet边界条件，与步骤(1)中确定的高原环境温度参数保持一致。
[0069]
其中水分传输方程考虑了气态水和液态水的两项传输作用，如下所示：
[0070][0071]
式中具体的pc为水分特征曲线，k
sl
为饱和状态下水在材料中的渗透系数。为孔隙率，η为水的粘度，k
rl
为液相的相对渗透系数，mw为水的摩尔质量，r为相对气体常数，ρ
l
为水的密度，d
g0
为水蒸气在空气中的自由扩散系数；p
gs
为饱和蒸气压。
[0072]
内部初始水分含量根据实际混凝土设定，范围在0～100rh％，边界条件选用dirichlet边界条件，与步骤(1)中确定的高原环境水分参数保持一致。
[0073]
其中盐离子传输仅考虑对流和扩散两个作用，不考虑化学反应时，控制方程如下：
[0074][0075]
式中具体的c为盐离子浓度，d0为浆体中的基准离子扩散系数，t0为基准离子扩散系数测定的时间，m为时间依赖系数，θc为孔隙临界饱和度，u为混凝土中的离子活化能，t
ref
为基准离子扩散系数测定时的温度，d为损伤因子，由损伤力学能量法求得，fr为修复材料对损伤的修复率，df为离子在水中的扩散系数，rb为混凝土对离子的结合影响因子。
[0076]
内部初始离子浓度设定为0，边界条件选用dirichlet边界条件，与步骤(1)中确定的高原环境离子参数保持一致。
[0077]
其中力作用方程如下：
[0078][0079]
式中具体的σ为应力，b为biot系数，p
*
为环境作用等效荷载，在高原环境中，主要包括冻融、硫酸盐等因素引起的内膨胀等效荷载。
[0080]
其中，高原环境下由冻融引起的内膨胀等效荷载按下列步骤求解：
①
先确定特定温度下的浆体孔隙中气态、液态和冰晶的含量；
②
依据气态、液态和冰晶含量求解计算内膨胀荷载。
[0081]
①
步骤计算按下式进行：
[0082][0083][0084][0085][0086]
式中具体的rc为对应温度下的临界冻结孔径，γ
cl
为水冰界面能，σm为水冰相变熵，δt为高原环境下实际温度与对应高原气压环境下溶液冰点的温度差，受高原气压、孔隙等影响，δ为孔壁边水膜厚度，厚度通过非饱和c-s-h凝胶孔中水分分布的分子动力学模拟获得，s
l0
(r)、sc(r)、s
l
(r)分别为对应孔径下初始液态水饱和度函数、冰晶饱和度函数和液态水饱和度函数，sg为浆体孔隙中气体饱和度，r
sc
为临界饱和孔隙同样由分子动力学模拟获取，为孔径小于r
sc
的孔隙率。
[0087]
②
步骤计算按下式进行：
[0088][0089][0090][0091]
式中具体的p
ft*
为冻融引起的内膨胀荷载力，r
max
和r
min
为基体内最大和最小孔径，k、k
l
、kc分别为浆体、水和冰的基体模量，∈为浆体的应变，b
l
、bc分别为水和冰的biot系数，
αs为浆体的热膨胀系数，p
atm
为高原环境气压，p
atm0
为浇筑浆体时的环境气压，m为计算系数，sc、s
l
、sg分别为对应温度下冰、水和气孔的体积分数。
[0092]
其中，高原环境下由硫酸盐等侵蚀性离子侵入并发生化学反应，在浆体孔隙中产生膨胀性侵蚀介质产生的内膨胀力按下列步骤求解：
①
考虑化学反应作用，求解不同离子在浆体中的浓度分布；
②
依据离子浓度分布，计算化学反应体积膨胀力。
[0093]
①
以硫酸盐侵蚀步骤计算按下式进行：
[0094][0095][0096][0097][0098]
式中具体的分别为硫酸根离子、钙离子、氢氧化钙、石膏和钙离子的浓度，k、k1、k2为化学反应系数，为化学反应系数，为硫酸根离子和钙离子扩散系数，q为混凝土腐蚀系数，为钙离子活度。
[0099]
②
步骤计算按下式进行：
[0100][0101][0102]
式中具体的p
s*
为硫酸盐侵蚀引起的内膨胀等效荷载，v
aft
为膨胀钙矾石体积分数，由步骤
①
中计算浓度结合化学反应方程确定。
[0103]
其他镁盐等腐蚀性产物在浆体中的内膨胀力可按其化学反应方程，结合步骤
①
进行计算。
[0104]
当高原环境同时存在冻融和硫酸盐等侵蚀介质作用时，浆体所受最终的内膨胀荷载计算如下：
[0105][0106]
式中具体的ωi为对应环境作用影响因子，p
i*
分别代表冻融、硫酸盐、镁盐等环境等效荷载。
[0107]
所述步骤(4)求解上述高原环境考虑防护与修复作用的混凝土数值模型，可采用有限元法进行求解。
[0108]
所述步骤(5)以钢筋开始锈蚀时的寿命为混凝土的最终服役寿命。通过步骤(1-4)求解不同深度处的离子浓度，对比临界浓度确定混凝土服役寿命。钢筋锈蚀修复材料会提高钢筋的临界离子浓度，因此混凝土服役寿命计算如下：
[0109][0110]
式中具体的tc为混凝土服役寿命，xd为钢筋埋置深度，f
rr
为钢筋锈蚀修复材料作用因子，ccr为钢筋固有的锈蚀临界浓度。
[0111]
实施例一
[0112]
如图1-10所示，本实施例中提供一种高原环境考虑防护与修复效果的混凝土寿命计算方法，具体过程如下：
[0113]
选择我国一处服役于高原环境的混凝土墩柱，墩柱接触水，温度范围为-8～20℃，氯盐浓度为20000mg/l，无硫酸盐，不承受额外的外荷载作用，使用不同类型的防护涂层(绿化聚乙烯涂层cpe/氯丁乳液涂层cr/氯偏乳液涂层pvdc/氯化聚氯乙烯涂层cpvc)，使用不同裂缝损伤修复材料(修复率50％/70％/90％)，使用不同的钢筋锈蚀修复材料(临界浓度提高10％/30％)。
[0114]
依据步骤(1-4)构建高原环境下考虑防护作用的数值模型(图1)，并求解获得各参数在混凝土中的分布情况，以最易引起钢筋锈蚀的氯离子浓度为例(图2-8)。
[0115]
由图2可知，在2年内，水分和离子均未穿过防护涂层，在第5年时，混凝土表层骨料较少的局部区域中出现了少量水分和氯离子，但在防护涂层生效的10年内，混凝土内部并未出现大量的水分和氯离子，这也表明了在防护涂层的使用寿命内防护涂层可有效隔绝水和离子的作用，直接提高混凝土抗侵蚀性能。在防护涂层的失效后，相当于水分和氯离子直接接触了混凝土表面，随即混凝土内部的水分和氯离子含量大幅增加。由于氯离子传输主要以水位载体，且混凝土内部存在吸附结合作用，水分传输速率远超过混凝土，在30年时，混凝土内部的整体饱和度接近于1，但氯离子仅在表层存在。使用服役寿命为10年的表层防护涂层，在15年、30年、50年、100年时，氯离子传输深度分别为8mm、37mm、34mm、40mm和48mm。
[0116]
由图3可知，混凝土氯离子浓度随时间的推移而不断增长，且增长趋势由急剧上升到缓慢提高。防护涂层的使用推迟了混凝土中氯离子的侵蚀时间点，延缓了氯离子在混凝土中的扩散速度。混凝土氯离子浓度与防护涂层寿命密切相关。随着防护涂层寿命的增加，混凝土中氯离子开始侵入的时间点逐渐被推迟，混凝土氯离子浓度也呈线性下降，表明防护涂层寿命的增加可有效降低混凝土使用过程全寿命周期中所含的氯离子浓度，提高防护涂层的寿命是阻隔延缓氯离子侵蚀和提高混凝土服役寿命的有效手段。
[0117]
由图4可知，防护涂层厚度的增加可以有效降低混凝土表面以及内部最终的氯离子浓度，且对混凝土内部氯离子浓度的降低比混凝土表面更加明显。当防护涂层厚度在1mm以内时，混凝土表面开始出现氯离子的时间点并没有被推迟；而当防护涂层厚度超过1mm时，混凝土表面氯离子浓度才开始逐渐上升。表明防护涂层厚度至少应在1mm以上时才可在涂层预期寿命内有效阻隔介质侵入混凝土中，但即使防护涂层厚度逐渐增加至3mm，依然会存在氯离子侵入混凝土表面。
[0118]
由图5可知，与未使用防护涂层相比，使用氯化聚乙烯涂层(cpe)、氯丁乳液涂层(cr)、氯偏乳液涂层(pvdc)和氯化聚氯乙烯涂层(cpvc)四种防护涂层对混凝土的最终氯离子浓度分别降低了86％、40％、70％和84％。其中，cpe、cr、pvdc和cpvc四种防护涂层均可有
效延缓混凝土氯离子侵入时间，降低混凝土内部氯离子浓度，提高混凝土服役寿命。
[0119]
由图6-8可知，损伤修复材料的修复率为70％时，对于100年时30mm处的氯离子含量为未修复时的73％，损伤修复材料的损伤修复率为90时，可将氯离子含量降低至未修复时的35％。但损伤修复率从70％～90％的变化，对于5mm内的侵蚀介质传输影响不大，因为并未完全修复损伤，由于溶液在表层裂纹中的抽吸作用，侵蚀介质随着通道早已快速侵入混凝土内部。修复材料的修复率对于10～30mm范围内的侵蚀介质浓度降低有较大提高。
[0120]
分别考虑防护涂层、损伤修复材料和钢筋锈蚀修复材料的作用效应，依据步骤(5)求解混凝土服役寿命(图9-10)。
[0121]
由图9可知，由环境引发损伤的混凝土服役寿命为50年，若使用0.3mm的防护涂层后，寿命可提高只56年，若使用3mm的防护涂层后，寿命可提高至70年。
[0122]
由图10可知，对于不使用防护材料的混凝土，考虑损伤修复效应和锈蚀修复效应，使用损伤修复率为80％的修复材料下，可达到78年服役寿命，若使用可使临界浓度提高10％的钢筋锈蚀修复材料时，服役寿命可达63年，若同时使用两类修复材料，可保障混凝土达到100年的服役寿命。
[0123]
因此，本方法可有效计算高原环境下考虑防护与修复作用的混凝土服役寿命，为工程实践提供重要的指导作用。
[0124]
以上所述，仅为本技术较佳的具体实施方式，但本技术的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本技术揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本技术的保护范围之内。因此，本技术的保护范围应该以权利要求的保护范围为准。

技术特征：
1.一种高原环境考虑防护与修复效果的混凝土寿命计算方法，其特征在于，包括以下步骤：s1.确定模拟高原环境的温度、湿度、盐浓度和内、外力荷载作用情况，建立高原环境、荷载参数与时间的映射关系；s2.构建高原环境下考虑防护与修复作用的混凝土几何模型，并确定所述几何模型中各物相的材料属性；s3.建立温度-湿度-盐浓度-荷载多因素耦合作用物理场方程，在所述混凝土几何模型的基础上设定内部、外部边界条件和初始条件，其中，所述外部边界条件为所述高原环境、荷载参数与时间的映射关系；s4.求解高原环境和荷载作用下考虑防护与修复作用的混凝土内部温度、湿度、盐浓度、应力分布及其随时间的变化情况；s5.计算考虑防护与修复作用的混凝土内部的钢筋达到临界锈蚀浓度发生脱钝状态的时间，确定混凝土服役寿命。2.根据权利要求1所述的高原环境考虑防护与修复效果的混凝土寿命计算方法，其特征在于，所述s1中，所述高原环境、荷载参数与时间的映射关系，为固定不变或任意形式变化的温度、湿度、盐浓度和内、外力荷载；盐为侵蚀介质，包括：氯化钠、氯化钙、氯化铯、硫酸钠、硫酸镁；内、外力荷载包括：重力、内部孔隙膨胀力、外部荷载。3.根据权利要求1所述的高原环境考虑防护与修复效果的混凝土寿命计算方法，其特征在于，所述s2中，所述混凝土几何模型的构建过程为：采用蒙特卡洛方法随机生成混凝土细观结构，在混凝土表层添加任意厚度的防护涂层几何模型，所述防护涂层几何模型在1维结构中为与混凝土相连的线段，在2维模型中为与混凝土贴合的四面体，在3维结构模型中为与混凝土相贴合的六面体，损伤修复材料是均匀分布在水泥浆体内部的点阵或空间域。4.根据权利要求1所述的高原环境考虑防护与修复效果的混凝土寿命计算方法，其特征在于，所述s2中，所述混凝土几何模型包含粗骨料、细骨料和水泥浆体，粗细骨料级配类型可自由设置，骨料体积分数在0％～80％之间，混凝土表面防护涂层是粘附在混凝土表面的连续涂层，混凝土修复材料以场域形式均匀分布在水泥浆体内部。5.根据权利要求1所述的高原环境考虑防护与修复效果的混凝土寿命计算方法，其特征在于，所述s2中，所述物相包括：骨料、浆体、涂层；所述浆体的材料属性包括：弹性模量、泊松比、密度、抗拉强度、断裂能、相变温度、相变温度间隔、相变潜热、导热系数、恒压热容、离子扩散系数、孔隙率、孔径分布；所述涂层的材料属性包括：导热系数、恒压热容、离子扩散系数；所述骨料和所述涂层的参数通过实验测试获得；所述浆体由砂和水泥浆体组成，浆体的参数由细观力学等效介质理论计算得到。6.根据权利要求1所述的高原环境考虑防护与修复效果的混凝土寿命计算方法，其特征在于，
所述s3中，所述多因素耦合作用物理场方程包括：温度传导方程、水分传输方程、盐传输方程和力作用方程。7.根据权利要求6所述的高原环境考虑防护与修复效果的混凝土寿命计算方法，其特征在于，所述温度传导方程为：其中，ρ为密度，c为比热容，t为温度，t为时间，λ为导热系数，l为相变潜热，θ
l
为含水率；内部初始温度设定为与高原环境初始温度保持一致，边界条件选用dirichlet边界条件，与s1中确定的高原环境温度参数保持一致；所述水分传输方程考虑了气态水和液态水的两项传输作用，表达如下：其中，p
c
为水分特征曲线，k
sl
为饱和状态下水在材料中的渗透系数，为孔隙率，η为水的粘度，k
rl
为液相的相对渗透系数，m
w
为水的摩尔质量，r为相对气体常数，ρ
l
为水的密度，d
g0
为水蒸气在空气中的自由扩散系数；p
gs
为饱和蒸气压；内部初始水分含量根据实际混凝土设定，范围在0～100rh％，边界条件选用dirichlet边界条件，与s1中确定的高原环境水分参数保持一致；盐离子传输仅考虑对流和扩散两个作用，不考虑化学反应时，所述盐传输方程如下：其中，c为盐离子浓度，d0为浆体中的基准离子扩散系数，t0为基准离子扩散系数测定的时间，m为时间依赖系数，θ
c
为孔隙临界饱和度，u为混凝土中的离子活化能，t
ref
为基准离子扩散系数测定时的温度，d为损伤因子，由损伤力学能量法求得，f
r
为修复材料对损伤的修复率，d
f
为离子在水中的扩散系数，r
b
为混凝土对离子的结合影响因子；内部初始离子浓度设定为0，边界条件选用dirichlet边界条件，与s1中确定的高原环境离子参数保持一致；所述力作用方程如下：其中，σ为应力，b为biot系数，p
*
为环境作用等效荷载，所述环境作用等效荷载包括：冻融、硫酸盐因素引起的内膨胀等效荷载；融、硫酸盐因素引起的内膨胀等效荷载；
其中，p
ft*
为冻融引起的内膨胀荷载力，r
max
和r
min
为基体内最大和最小孔径，k、k
l
、k
c
分别为浆体、水和冰的基体模量，∈为浆体的应变，b
l
、b
c
分别为水和冰的biot系数，α
s
为浆体的热膨胀系数，p
atm
为高原环境气压，p
atm0
为浇筑浆体时的环境气压，s
c
、s
l
、s
g
分别为对应温度下冰、水和气孔的体积分数，p
s*
为硫酸盐侵蚀引起的内膨胀等效荷载，v
aft
为膨胀钙矾石体积分数，ω
i
为对应环境作用影响因子，p
i*
分别代表冻融、硫酸盐、镁盐环境等效荷载。8.根据权利要求1所述的高原环境考虑防护与修复效果的混凝土寿命计算方法，其特征在于，所述s5具体为：以钢筋开始锈蚀时的寿命为混凝土的最终服役寿命，通过s1-s4求解不同深度处的离子浓度，对比临界浓度确定混凝土服役寿命，混凝土服役寿命计算如下：其中，t
c
为混凝土服役寿命，x
d
为钢筋埋置深度，f
rr
为钢筋锈蚀修复材料作用因子，ccr为钢筋固有的锈蚀临界浓度。

技术总结
本发明公开了一种高原环境考虑防护与修复效果的混凝土寿命计算方法，包括：建立高原环境、荷载参数与时间的映射关系；构建混凝土几何模型，确定几何模型中各物相的材料属性；建立多因素耦合作用物理场方程，在混凝土几何模型的基础上设定边界条件和初始条件；求解混凝土内部温度、湿度、盐浓度、应力分布及其随时间的变化情况；计算混凝土内部的钢筋达到临界锈蚀浓度发生脱钝状态的时间，确定混凝土服役寿命；本发明模拟不同高原环境下使用单一的表面防护涂层、裂缝修复材料和钢筋锈蚀修复材料或同时使用几类防护与修复材料后混凝土内部温度、湿度、应力响应和服役寿命，从而为高原环境下混凝土结构设计与施工提供依据。境下混凝土结构设计与施工提供依据。境下混凝土结构设计与施工提供依据。


技术研发人员：王凤娟 王赟程 蒋金洋 穆松 刘志勇 李映泽
受保护的技术使用者：东南大学
技术研发日：2023.05.25
技术公布日：2023/9/23