CN113935220B - 基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法 - Google Patents
基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法 Download PDFInfo
- Publication number
- CN113935220B CN113935220B CN202111462206.2A CN202111462206A CN113935220B CN 113935220 B CN113935220 B CN 113935220B CN 202111462206 A CN202111462206 A CN 202111462206A CN 113935220 B CN113935220 B CN 113935220B
- Authority
- CN
- China
- Prior art keywords
- concrete
- diffusion
- tunnel
- submarine
- chloride ion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000004567 concrete Substances 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 30
- VEXZGXHMUGYJMC-UHFFFAOYSA-M Chloride anion Chemical compound [Cl-] VEXZGXHMUGYJMC-UHFFFAOYSA-M 0.000 claims abstract description 118
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 37
- 239000010410 layer Substances 0.000 claims abstract description 10
- 239000011241 protective layer Substances 0.000 claims abstract description 7
- 238000009792 diffusion process Methods 0.000 claims description 55
- 239000011148 porous material Substances 0.000 claims description 41
- 230000004907 flux Effects 0.000 claims description 30
- 229920006395 saturated elastomer Polymers 0.000 claims description 14
- 230000035699 permeability Effects 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 8
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- HPNSNYBUADCFDR-UHFFFAOYSA-N chromafenozide Chemical compound CC1=CC(C)=CC(C(=O)N(NC(=O)C=2C(=C3CCCOC3=CC=2)C)C(C)(C)C)=C1 HPNSNYBUADCFDR-UHFFFAOYSA-N 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 229910052759 nickel Inorganic materials 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000012423 maintenance Methods 0.000 abstract description 2
- 230000007613 environmental effect Effects 0.000 abstract 1
- 230000008859 change Effects 0.000 description 9
- 206010019233 Headaches Diseases 0.000 description 6
- 230000009471 action Effects 0.000 description 5
- 230000037427 ion transport Effects 0.000 description 4
- 239000013535 sea water Substances 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000003628 erosive effect Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 108091006671 Ion Transporter Proteins 0.000 description 1
- 102000037862 Ion Transporter Human genes 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 239000004574 high-performance concrete Substances 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000003334 potential effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/38—Concrete; Lime; Mortar; Gypsum; Bricks; Ceramics; Glass
- G01N33/383—Concrete or cement
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- General Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Mathematical Optimization (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Fluid Mechanics (AREA)
- Algebra (AREA)
- Ceramic Engineering (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于扩散‑对流模型的海底混凝土隧道耐久寿命预测方法，包括步骤：建立海底隧道非饱和混凝土中氯离子扩散‑对流模型；判别海底隧道保护层处的最危险点；进行海底隧道耐久寿命预测。本发明的有益效果是：通过分析任意时间节点任意位置处的自由氯离子浓度分布，来确定隧道结构保护层位置处自由氯离子浓度最大的危险点，以及该点浓度达到临界氯离子浓度的服役时长；本发明的预测方法准确性更高，为海底隧道的耐久性设计提供科学的依据，对工程的耐久性能维护管理具有指导作用，并且为今后的海底混凝土隧道耐久设计研究提供了理论基础；适用于多种工况，包括不同环境自由浓度的海域、以及海底隧道的不同水头和混凝土初始饱和度。
Description
技术领域
本发明属于海底混凝土隧道耐久寿命预测领域，尤其涉及一种基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法。
背景技术
海底隧道是不可逆的重要工程，设计使用年限一般要求在100年以上。海底隧道长期处于高压力海水侵蚀的恶劣环境，氯离子侵蚀是其耐久寿命劣化的重要因素之一。氯离子在混凝土中的传输方式和传输速度是影响海底隧道耐久寿命的关键。目前国内外学者提出了多种基于氯离子扩散理论的海底隧道耐久寿命预测理论模型，以及采用遗传算法等隐式模糊预测模型。这些模型通常对氯离子扩散系数进行了多因素修正，然后基于Fick第一和第二扩散定律进行确定性或考虑模型不确定性的可靠度研究。
混凝土是一种孔隙材料，氯离子在混凝土孔隙溶液中除了依靠浓度梯度驱动的扩散方式外，还会以孔隙水为载体通过对流作用一起运移，包括压力渗透作用、毛细作用等。目前对海底隧道中氯离子的扩散-对流研究，有学者根据Duracrete模型直接在Fick扩散模型中设置了对流层深度来实现，也有学者考虑压力渗透作用下的饱和混凝土扩散-对流模型，基于Fick扩散定律和Darcy定律建立了考虑海水压力渗透作用的扩散-对流模型，并采用Ogata&Banks的解析解进行分析。丹麦、瑞典等水下混凝土工程的取芯实验发现，即使服役了数十年的水下结构，也仅在距表面30mm深度内孔隙液态水含量完全饱和，超过该深度的混凝土处于非饱和状态。非饱和混凝土中水分运动的主要驱动势是反映毛细吸收作用的基质势，相应地非饱和混凝土中氯离子的扩散-对流模型应同时考虑扩散、外水压力和基质势作用。
综上所述，目前还没有基于非饱和混凝土中氯离子传输机制的扩散-对流模型应用于海底隧道耐久寿命预测的研究成果；目前通常采用氯离子扩散模型来预测海底隧道耐久寿命，忽略了海底隧道服役环境中外水压力，以及水下非饱和混凝土的毛细吸力作用对混凝土孔隙中水分的对流驱动，因此可能低估了海底隧道中自由氯离子的浓度，相应的可能高估了海底隧道的耐久寿命。
发明内容
本发明的目的是克服现有技术中的不足，提供一种基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法：
步骤1、建立海底隧道非饱和混凝土中氯离子扩散-对流模型；
步骤1.1、基于Fick第一和第二定律，将混凝土孔隙溶液中自由氯离子的扩散通量作为扩散模型；
上式中，Jd为自由氯离子的扩散通量，单位为kg/(m2s)；为Nabla算子；Cf为孔隙溶液中自由氯离子浓度，单位为kg/m3；t为扩散时间，单位为s；Deff为氯离子有效扩散系数，单位为m2/s；
根据Atkinson和Nickerson模型，氯离子的有效扩散系数为：
Deff＝Ds0φβ
上式中，Ds0为饱和混凝土的氯离子扩散系数，单位符号为m2/s；φ为孔隙率；β为孔结构参数β；
孔结构参数为：
β＝φ1/3s10/3
上式中：s为孔隙水饱和度；
步骤1.2、将混凝土孔隙溶液中自由氯离子的对流通量Jc作为对流模型；
Jc＝-Cfu
上式中，Jc为自由氯离子的对流通量，单位为kg/(m2s)；u为混凝土孔隙中水的流速，单位为m/s；
u按下式计算：
上式中：Ks为饱和渗透系数，单位为m/s；ρw为水的密度，单位为kg/m3；g为重力加速度，单位为m/s2；kr为相对渗透率，0≤kr≤1，0代表完全干燥，1代表完全饱和；θ为含水率；ψ为驱动势，ψ＝ψm+ψp，单位为Pa；其中ψP为渗流路径两侧压力差形成的压力势，单位为Pa；ψm为单位容积水分的基质势，单位为Pa；kr(θ)表示相对渗透率的大小与含水率的关系；
单位容积水分的基质势ψm的计算模型为：
ψm＝-α(s-1/m-1)1-m
上式中：α，m为混凝土材料的经验参数；
步骤1.3、将总的自由氯离子通量Jcl作为非饱和混凝土中氯离子扩散-对流模型：
上式中，Jcl为总的自由氯离子通量，单位为kg/(m2s)；
步骤2、判别海底隧道保护层处的最危险点；
步骤3、进行海底隧道耐久寿命预测；取海底隧道保护层处的最危险点氯离子浓度Cf(x，t)达到临界自由氯离子浓度Cr的最短时间，作为海底隧道耐久寿命T1；Cf(x，t)表示扩散时间t时刻距隧道外壁的距离x处的孔隙溶液中自由氯离子浓度Cf；Cf(x，t)由步骤1中的扩散-对流模型计算得到。
作为优选，步骤1.2中，当混凝土处于完全饱和状态时，混凝土孔隙中水的流速u由Darcy定律计算：
作为优选，步骤1.2中，当混凝土处于非饱和状态时，混凝土孔隙中水的流速u由Richards方程计算：
作为优选，步骤2具体包括以下步骤：
步骤2.1、首先根据海底隧道外壁初始时刻的水头压力P0，以及混凝土初始饱和度s，分别计算各点处的渗流路径两侧压力差形成的压力势ψP和单位容积水分的基质势ψm；t为扩散时间；
步骤2.2、将步骤2.1计算得到的ψP、ψm代入对流模型Jc＝-Cfu，采用积分有限差分法计算得到下一个时间节点的混凝土孔隙中水的流速u和混凝土孔隙水饱和度s；
步骤2.3、代入海底隧道外壁初始氯离子浓度和隧道内混凝土的初始氯离子浓度，计算得到第一个时间节点各位置处的对流通量Jc、内外氯离子浓度梯度确定的扩散通量Jd，以及总的自由氯离子通量Jcl；并将Jc、Jd和Jcl作为下一个时间节点计算的初始数据；
步骤2.4、返回执行步骤2.1至步骤2.3，获得任意时间节点任意位置处的自由氯离子浓度值；最后将海底隧道保护层处氯离子浓度Cf(x，t)最大的点作为整个海底隧道保护层处的最危险点；Cf(x，t)表示扩散时间t时刻距隧道外壁的距离x处的孔隙溶液中自由氯离子浓度Cf；Cf(x，t)由步骤1中的扩散-对流模型计算得到。
本发明的有益效果是：
本发明是一种从外水压力作用下非饱和混凝土中氯离子传输机制入手，采用扩散-对流模型预测海底隧道耐久寿命的方法，该预测方法准确性更高，为海底隧道的耐久性设计提供科学的依据，对工程的耐久性能维护管理具有指导作用，并且为今后的海底混凝土隧道耐久设计研究提供了理论基础。
本发明基于非饱和混凝土中氯离子传输机制建立了海底隧道非饱和混凝土中氯离子输运的扩散-对流模型，通过分析任意时间节点任意位置处的自由氯离子浓度分布，来确定隧道结构保护层位置处自由氯离子浓度最大的危险点，以及该点浓度达到临界氯离子浓度的服役时长，本发明可以适用于多种工况，包括不同环境自由浓度的海域、以及海底隧道的不同水头和混凝土初始饱和度。
附图说明
图1-1至图1-4分别为本发明隧道管段内第10年、第50年、第100年、第133年的氯离子分布云图和最危险点示意图；
图2为本发明不同水头压力下危险点处氯离子浓度变化图和隧道耐久寿命曲线图；
图3为本发明不同混凝土初始饱和度下危险点处氯离子浓度变化图和隧道耐久寿命曲线图；
图4为本发明不同环境自由氯离子浓度下危险点处氯离子浓度变化图和隧道耐久寿命曲线图；
图5为本发明扩散模型、饱和扩散-对流模型、非饱和扩散-对流模型下危险点处氯离子浓度变化图和隧道耐久寿命曲线图。
具体实施方式
下面结合实施例对本发明做进一步描述。下述实施例的说明只是用于帮助理解本发明。应当指出，对于本技术领域的普通人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
本申请实施例一提供了一种基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法：
步骤1、考虑海底隧道的服役环境特点和混凝土的非饱和特性，建立海底隧道非饱和混凝土中氯离子扩散-对流模型；
步骤1.1、基于Fick第一和第二定律，将混凝土孔隙溶液中自由氯离子的扩散通量作为扩散模型；
上式中，Jd为自由氯离子的扩散通量，单位为kg/(m2s)；为Nabla算子；Cf为孔隙溶液中自由氯离子浓度，单位为kg/m3；t为扩散时间，单位为s；Deff为氯离子有效扩散系数，单位为m2/s；
根据Atkinson和Nickerson模型，氯离子的有效扩散系数为：
Deff＝Ds0φβ
上式中，Ds0为饱和混凝土的氯离子扩散系数，单位符号为m2/s；φ为孔隙率；β为孔结构参数β；
孔结构参数为：
β＝φ1/3s10/3
上式中：s为孔隙水饱和度；
步骤1.2、将混凝土孔隙溶液中自由氯离子的对流通量Jc作为对流模型；
Jc＝-Cfu
上式中，Jc为自由氯离子的对流通量，单位为kg/(m2s)；u为混凝土孔隙中水的流速，单位为m/s；
u按下式计算：
上式中：Ks为饱和渗透系数，单位为m/s；ρw为水的密度，单位为kg/m3；g为重力加速度，单位为m/s2；kr为相对渗透率，0≤kr≤1，0代表完全干燥，1代表完全饱和；θ为含水率；ψ为驱动势，ψ＝ψm+ψp，单位为Pa；其中ψP为渗流路径两侧压力差形成的压力势，单位为Pa；ψm为单位容积水分的基质势，单位为Pa；kr(θ)表示相对渗透率的大小与含水率的关系；
单位容积水分的基质势ψm的计算模型为：
ψm＝-α(s-1/m-1)1-m
上式中：α，m为混凝土材料的经验参数；
当混凝土处于完全饱和状态时，混凝土孔隙中水的流速u由Darcy定律计算：
当混凝土处于非饱和状态时，混凝土孔隙中水的流速u由Richards方程计算：
步骤1.3、海底隧道外壁接触海水，内壁接触大气，内外氯离子浓度梯度会引起扩散通量Jd；海底隧道混凝土在初始服役时处于非饱和状态，随着海水入渗混凝土饱和度逐渐提高，内部逐渐形成饱和区和非饱和区，对流通量Jc的驱动势也会从基质势向压力势转换；因此海底隧道混凝土中氯离子传输受到扩散和对流共同驱动，将总的自由氯离子通量Jcl作为非饱和混凝土中氯离子扩散-对流模型：
上式中，Jcl为总的自由氯离子通量，单位为kg/(m2s)；
步骤2、判别海底隧道保护层处的最危险点；
步骤2.1、首先根据海底隧道外壁初始时刻(t＝0)的水头压力P0，以及混凝土初始饱和度s，分别计算各点处的渗流路径两侧压力差形成的压力势ψP和单位容积水分的基质势ψm；t为扩散时间；
步骤2.2、将步骤2.1计算得到的ψP、ψm代入对流模型Jc＝-Cfu，采用积分有限差分法计算得到下一个时间节点(t＝1)的混凝土孔隙中水的流速u和混凝土孔隙水饱和度s；
步骤2.3、代入海底隧道外壁初始氯离子浓度和隧道内混凝土的初始氯离子浓度，计算得到第一个时间节点各位置处的对流通量Jc、内外氯离子浓度梯度确定的扩散通量Jd，以及总的自由氯离子通量Jcl；并将Jc、Jd和Jcl作为下一个时间节点计算的初始数据；
步骤2.4、返回执行步骤2.1至步骤2.3，获得任意时间节点任意位置处的自由氯离子浓度值；最后将海底隧道保护层处氯离子浓度Cf(x，t)最大的点作为整个海底隧道保护层处的最危险点；Cf(x，t)表示扩散时间t时刻距隧道外壁的距离x处的孔隙溶液中自由氯离子浓度Cf；Cf(x，t)由步骤1中的扩散-对流模型计算得到。
步骤3、进行海底隧道耐久寿命预测；取海底隧道保护层处的最危险点氯离子浓度Cf(x，t)达到临界自由氯离子浓度Cr的最短时间，作为海底隧道耐久寿命T1。
本实施例可以根据海底隧道服役环境的水头压力P0和环境自由氯离子浓度Cs(x)，以及混凝土的材料参数，如初始饱和度s0、孔隙率φ、饱和渗透系数Ks、氯离子有效扩散系数Deff等分别建立氯离子输运的扩散模型和扩散-对流模型，计算得到任意时间节点任意位置处的自由氯离子浓度值Cf(x，t)；Cf(x，t)表示扩散时间t时刻距隧道外壁的距离x处的孔隙溶液中自由氯离子浓度Cf；Cf(x，t)由步骤1中的扩散-对流模型计算得到。
本实施例可以通过比较自由氯离子浓度，判别海底隧道保护层位置处的最危险点；通过比较扩散模型和扩散-对流模型中海底隧道最危险点处自由氯离子浓度达到临界自由氯离子浓度Cr的最短时间，来确定海底隧道的耐久寿命T1。
实施例二
在实施例一的基础上，本申请实施例二提供了实施例一中基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法在舟山沈家门海底隧道中的应用情况：
本实例以舟山沈家门海底隧道进行算例分析，该隧道水头压力15m，环境自由氯离子浓度为0.005，混凝土初始饱和度为0.70，饱和渗透系数为1.29e-12，氯离子扩散系数为1.69e-12，保护层厚度为0.05m，临界氯离子浓度为0.1。
图1-1至图1-4为根据实施例1提出的非饱和扩散-对流模型计算得到的海底隧道管段内不同时间节点的氯离子浓度分布云图。顶板与侧墙相交的腋角拐点部位氯离子浓度增长最为明显，因此将该点作为危险点。
图2为通过改变隧道外壁的水头压力取值，根据Darcy方程改变隧道外壁和内部的压力势▽ψP，根据实施例1的步骤1和步骤2分别计算出危险点处氯离子浓度随时间的变化规律，根据步骤3确定该浓度达到临界氯离子浓度的时间作为海底隧道的耐久寿命T1。
图3为通过改变隧道混凝土的初始饱和度取值，根据Richards方程改变隧道外壁和内部的压力势▽ψm，根据实施例1的步骤1和步骤2分别计算出危险点处氯离子浓度随时间的变化规律，根据步骤3确定该浓度达到临界氯离子浓度的时间作为海底隧道的耐久寿命T1。
图4为通过改变隧道外壁的环境自由氯离子浓度取值，根据Fick扩散定律改变隧道的扩散通量Jd，根据实施例1的步骤1和步骤2分别计算出危险点处氯离子浓度随时间的变化规律，根据步骤3确定该浓度达到临界氯离子浓度的时间作为海底隧道的耐久寿命T1。
图5为分别采用扩散模型、饱和扩散-对流模型和非饱和扩散-对流模型，根据实施例1的步骤1和步骤2分别计算出危险点处氯离子浓度随时间的变化规律，根据步骤3确定只考虑氯离子浓度差引起的扩散模型计算得到的耐久寿命为247年，按饱和混凝土的扩散-对流模型计算得到的耐久寿命为225年，按非饱和混凝土扩散-对流模型计算得到的耐久寿命为133年，故最终预测的隧道耐久寿命T1为133年。
本实施例根据海底隧道内外壁存在压力差的服役特点，以及高性能混凝土的非饱和性，建立了海底混凝土隧道氯离子扩散-对流运移模型和耐久寿命预测模型，按照扩散-对流模型预测得到的隧道耐久寿命均小于扩散模型，说明只按扩散模型进行耐久寿命预测可能偏于不安全。两种氯离子扩散-对流模型中，考虑了混凝土初始非饱和状态的预测寿命要低于饱和扩散-对流模型的预测寿命，说明毛细吸力对氯离子运动的驱动作用比仅考虑压力驱动的作用更为明显。因此建议采用非饱和混凝土的扩散-对流模型进行海底隧道耐久寿命预测，可以为海底隧道工程提供准确安全的耐久寿命预测。
Claims (3)
1.基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法，其特征在于，包括以下步骤：
步骤1、建立海底隧道非饱和混凝土中氯离子扩散-对流模型；
步骤1.1、基于Fick第一和第二定律，将混凝土孔隙溶液中自由氯离子的扩散通量作为扩散模型；
上式中，Jd为自由氯离子的扩散通量，单位为kg/(m2s)；▽为Nabla算子；Cf为孔隙溶液中自由氯离子浓度，单位为kg/m3；t为扩散时间，单位为s；Deff为氯离子有效扩散系数，单位为m2/s；
根据Atkinson和Nickerson模型，氯离子的有效扩散系数为：
Deff＝Ds0φβ
上式中，Ds0为饱和混凝土的氯离子扩散系数，单位符号为m2/s；φ为孔隙率；β为孔结构参数β；
孔结构参数为：
β＝φ1/3s10/3
上式中：s为孔隙水饱和度；
步骤1.2、将混凝土孔隙溶液中自由氯离子的对流通量Jc作为对流模型；
Jc＝-Cfu
上式中，Jc为自由氯离子的对流通量，单位为kg/(m2s)；u为混凝土孔隙中水的流速，单位为m/s；
u按下式计算：
上式中：Ks为饱和渗透系数，单位为m/s；ρw为水的密度，单位为kg/m3；g为重力加速度，单位为m/s2；kr为相对渗透率，0≤kr≤1，0代表完全干燥，1代表完全饱和；θ为含水率；ψ为驱动势，ψ＝ψm+ψp，单位为Pa；其中ψP为渗流路径两侧压力差形成的压力势，单位为Pa；ψm为单位容积水分的基质势，单位为Pa；kr(θ)表示相对渗透率的大小与含水率的关系；
单位容积水分的基质势ψm的计算模型为：
ψm＝-α(s-1/m-1)1-m
上式中：α，m为混凝土材料的经验参数；
步骤1.3、将总的自由氯离子通量Jcl作为非饱和混凝土中氯离子扩散-对流模型：
上式中，Jcl为总的自由氯离子通量，单位为kg/(m2s)；
步骤2、判别海底隧道保护层处的最危险点；
步骤2具体包括以下步骤：
步骤2.1、首先根据海底隧道外壁初始时刻的水头压力P0，以及混凝土初始饱和度s，分别计算各点处的渗流路径两侧压力差形成的压力势ψP和单位容积水分的基质势ψm；t为扩散时间；
步骤2.2、将步骤2.1计算得到的ψP、ψm代入对流模型Jc＝-Cfu，采用积分有限差分法计算得到下一个时间节点的混凝土孔隙中水的流速u和混凝土孔隙水饱和度s；
步骤2.3、代入海底隧道外壁初始氯离子浓度和隧道内混凝土的初始氯离子浓度，计算得到第一个时间节点各位置处的对流通量Jc、内外氯离子浓度梯度确定的扩散通量Jd，以及总的自由氯离子通量Jcl；并将Jc、Jd和Jcl作为下一个时间节点计算的初始数据；
步骤2.4、返回执行步骤2.1至步骤2.3，获得任意时间节点任意位置处的自由氯离子浓度值；最后将海底隧道保护层处氯离子浓度Cf(x，t)最大的点作为整个海底隧道保护层处的最危险点；其中t为扩散时间，x为距隧道外壁的距离；Cf(x，t)表示扩散时间t时刻距隧道外壁的距离x处的孔隙溶液中自由氯离子浓度Cf；Cf(x，t)由步骤1中的扩散-对流模型计算得到；
步骤3、进行海底隧道耐久寿命预测；取海底隧道保护层处的最危险点氯离子浓度Cf(x，t)达到临界自由氯离子浓度Cr的最短时间，作为海底隧道耐久寿命T1；Cf(x，t)表示扩散时间t时刻距隧道外壁的距离x处的孔隙溶液中自由氯离子浓度Cf；Cf(x，t)由步骤1中的扩散-对流模型计算得到。
2.根据权利要求1所述基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法，其特征在于，步骤1.2中，当混凝土处于完全饱和状态时，混凝土孔隙中水的流速u由Darcy定律计算：
3.根据权利要求1所述基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法，其特征在于，步骤1.2中，当混凝土处于非饱和状态时，混凝土孔隙中水的流速u由Richards方程计算：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111462206.2A CN113935220B (zh) | 2021-12-02 | 2021-12-02 | 基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111462206.2A CN113935220B (zh) | 2021-12-02 | 2021-12-02 | 基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113935220A CN113935220A (zh) | 2022-01-14 |
| CN113935220B true CN113935220B (zh) | 2024-04-09 |
Family
ID=79288665
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111462206.2A Active CN113935220B (zh) | 2021-12-02 | 2021-12-02 | 基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113935220B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011196160A (ja) * | 2010-03-24 | 2011-10-06 | Tokyo Electric Power Co Inc:The | 地中鉄筋コンクリート構造物のひび割れ補修方法 |
| KR20150066000A (ko) * | 2013-12-05 | 2015-06-16 | 서울대학교산학협력단 | 해양 콘크리트 구조물의 복합열화 모의시험 방법 |
| KR20170111213A (ko) * | 2016-03-25 | 2017-10-12 | 강원대학교산학협력단 | 경화 콘크리트의 염해 내구성 평가 방법 및 염해 내구성 평가 컴퓨터 프로그램 |
| CN108333082A (zh) * | 2018-01-23 | 2018-07-27 | 广西大学 | 非饱和混凝土氯离子扩散系数多尺度预测模型的构建方法 |
| CN111797519A (zh) * | 2020-06-24 | 2020-10-20 | 河海大学 | 基于gem方程的混凝土钢筋的脱钝化时间的计算方法 |
| CN112580212A (zh) * | 2020-12-23 | 2021-03-30 | 北部湾大学 | 混凝土结构服役寿命迭代测定方法、系统、计算机设备 |
-
2021
- 2021-12-02 CN CN202111462206.2A patent/CN113935220B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011196160A (ja) * | 2010-03-24 | 2011-10-06 | Tokyo Electric Power Co Inc:The | 地中鉄筋コンクリート構造物のひび割れ補修方法 |
| KR20150066000A (ko) * | 2013-12-05 | 2015-06-16 | 서울대학교산학협력단 | 해양 콘크리트 구조물의 복합열화 모의시험 방법 |
| KR20170111213A (ko) * | 2016-03-25 | 2017-10-12 | 강원대학교산학협력단 | 경화 콘크리트의 염해 내구성 평가 방법 및 염해 내구성 평가 컴퓨터 프로그램 |
| CN108333082A (zh) * | 2018-01-23 | 2018-07-27 | 广西大学 | 非饱和混凝土氯离子扩散系数多尺度预测模型的构建方法 |
| CN111797519A (zh) * | 2020-06-24 | 2020-10-20 | 河海大学 | 基于gem方程的混凝土钢筋的脱钝化时间的计算方法 |
| CN112580212A (zh) * | 2020-12-23 | 2021-03-30 | 北部湾大学 | 混凝土结构服役寿命迭代测定方法、系统、计算机设备 |
Non-Patent Citations (3)
| Title |
|---|
| 基于氯离子渗透的海工混凝土寿命预测模型进展;刘志勇, 孙伟, 杨鼎宜, 周新刚;工业建筑;20040618(06);61-64 * |
| 海底隧道钢筋混凝土基于氯盐腐蚀的耐久性参数设计研究;曲立清;金祖权;赵铁军;李秋义;;岩石力学与工程学报;20071115(11);177-184 * |
| 海底隧道非饱和混凝土中氯离子渗透研究;李晓珍;硅酸盐通报;20140815;1929-1934 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113935220A (zh) | 2022-01-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Martın-Pérez et al. | A study of the effect of chloride binding on service life predictions | |
| Šavija et al. | Lattice modeling of rapid chloride migration in concrete | |
| Muthulingam et al. | Non-uniform time-to-corrosion initiation in steel reinforced concrete under chloride environment | |
| Jones et al. | Foam flow in a model porous medium: I. The effect of foam coarsening | |
| CN107290503A (zh) | 冻土中水非饱和流运动和溶质迁移通量及等效参数监测的方法 | |
| US11314909B2 (en) | Methods and systems for simulating multiphase flow through porous media | |
| Cardenes et al. | Analysis of the correlations between freeze–thaw and salt crystallization tests | |
| CN108957542B (zh) | 一种建立地震波衰减岩石物理图板的方法 | |
| CN113935220B (zh) | 基于扩散-对流模型的海底混凝土隧道耐久寿命预测方法 | |
| KR102479756B1 (ko) | 고압의 침지 환경에 노출된 철근 콘크리트 구조물의 철근 부식 예측 방법 | |
| CN104360412B (zh) | 致密深盆气成藏预测方法和装置 | |
| Carmeliet et al. | Hysteresis and moisture buffering of wood | |
| CN110095584B (zh) | 一种储层油水饱和度校正方法 | |
| Zhang et al. | Chloride transport in undersea concrete tunnel | |
| Hussain et al. | Computer-aided oxygen transport model of mass and energy simulation for corrosion of reinforced steel | |
| AL‐Ameeri et al. | Modelling chloride ingress into in‐service cracked reinforced concrete structures exposed to de‐icing salt environment and climate change: Part 1 | |
| CN204789265U (zh) | 桥梁结构混凝土气体渗透系数现场无损检测装置 | |
| CN112392471B (zh) | 一种碳酸盐岩储层孔隙度计算方法及装置 | |
| Brachaczek | Study of the impact of microstructure and sorption properties of the renovation plasters on the wall drying rate | |
| Frederiksen | On the need for more precise threshold values for chloride initiated corrosion | |
| Oliveira et al. | A multi-physics modelling based on coupled diffusion equations to simulate the carbonation process | |
| Padala et al. | Characterization of the Severity of Tidal Environment on Reinforced Concrete Members Using Moisture and Chloride Penetration Indices | |
| Hong et al. | Numerical prediction of chloride penetration in concrete exposed to a marine environment at tide | |
| Nicolai et al. | Development of a combined heat, moisture, and salt transport model for unsaturated porous building materials | |
| Yokoyama et al. | Flow path selection during capillary rise in rock: Effects of pore branching and pore radius variation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |