网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

硫酸盐侵蚀混凝土的数值模拟  PDF

  • 蒋金洋 1
  • 郑皓睿 2
  • 孙国文 2
  • 王凤娟 1
  • 刘志勇 1
1. 东南大学 材料科学与工程学院,江苏 南京 211189; 2. 石家庄铁道大学 材料科学与工程学院,河北 石家庄 050043

中图分类号: TU528.01

最近更新:2023-10-31

DOI:10.3969/j.issn.1007-9629.2023.10.001

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

为了客观评价硫酸盐侵蚀下混凝土传输-反应-损伤的全过程,基于结晶压理论、体积膨胀理论及Fick第二定律,建立了考虑孔隙率、曲折度和临界损伤程度的传输模型,并通过交替隐式差分法实现了硫酸盐传输变系数求解.结果表明:预测模型结果与试验结果基本吻合,最大误差为26.7%,可较好地预测混凝土中SO42-的扩散规律;硫酸盐质量分数和临界损伤程度对混凝土剥落具有较大的影响,相同临界损伤程度下5%硫酸钠溶液中混凝土的剥落速率比其在3%硫酸钠溶液中快51.8%;相同硫酸盐质量分数下临界损伤程度从0.80增大至0.95时,混凝土的剥落厚度减小了42.9%.

盐渍土或滨海盐土地区的普通混凝土结构往往1~2 a就发生损

1,致使许多基础设施提前失效,其劣化的主要原因是硫酸盐通过混凝土孔/缝进入混凝土后的传输-反应-损伤所致.对混凝土的传输过程和损伤程度进行量化,是提高结构混凝土耐久性设计和寿命预测的关键所在.为此,国内外学者致力于研究硫酸盐在混凝土中的传输-反应-损伤劣化模型.

关博文

2、刘超3建立了混凝土内部受硫酸盐侵蚀的损伤劣化模型;左晓宝4将Fick第二定律与硫酸盐反应机理相结合,建立了硫酸盐在混凝土内部传输的化学-力学分析模型,由此来判断混凝土的开裂状态;Zuo5、Bary6还建立了不同的化学-力学模型.上述模型均能较好地模拟混凝土内离子传输、侵蚀产物生长和试件膨胀变形等过程,但并未反映出硫酸盐侵蚀损伤混凝土的全过程.

本文基于结晶压理论、体积膨胀理论及Fick第二定律,建立了考虑孔隙率、曲折度和临界损伤程度的传输模型,并通过剥落厚度来表征混凝土的损伤劣化,以期能更精准地对混凝土在硫酸盐中侵蚀性能进行评估.

1 试验方法

1.1 原材料与试件制备

选用P∙I 42.5基准水泥,成型3组水胶比mW/mB(质量比)为0.5的C30混凝土,其水泥、水、砂子、石子的用量分别为370、185、640、1 205 kg,试件尺寸为100 mm×100 mm×100 mm.试件成型24 h后,将其脱模并置于恒温养护箱养护28 d;将试件两侧蜡封以控制SO42-的扩散方向,再将其置于硫酸钠质量分数w(Na2SO4)为1%、3%、5%的溶液中浸泡.

1.2 测试方法

根据GB/T 2899—2017《工业沉淀硫酸钡》,采用硫酸钡沉淀法测试混凝土中SO42-的浓度c.

2 硫酸盐侵蚀混凝土模型的建立与求解

2.1 硫酸盐侵蚀过程及其机理

硫酸盐对混凝土的侵蚀损伤理论有体积膨胀理论、结晶压理论、吸水肿胀理论和拓扑反应理

7‑8.Ikumi9、Yin10通过大量试验和理论模拟证明,体积膨胀和盐结晶相结合,能更客观地反映硫酸盐传输-反应-损伤的全过程.

全侵蚀过程中硫酸盐侵蚀下,硫酸盐与其混凝土孔溶液中的铝酸盐反应形成钙钒石(AFt)或者石膏膨胀产

11.当这些膨胀产物达到过饱和时,形成晶体,并填充于不同尺度孔隙,主要是毛细孔4‑5,达到某临界点后产物继续膨胀,当混凝土自身的抗压强度不足以抵抗产物所产生的体积膨胀应力时,发生开裂、逐层破坏剥落.侵蚀反应过程中发生的化学反应12

2Ca2++2SO42-+C4AS¯H12+20H2OC6AS¯3H32 (1)
Ca2++SO42-+2H2OCS¯H2 (2)

2.2 侵蚀模型

硫酸盐侵蚀混凝土的初始和边界条件见图1.侵蚀模型的边界条件为:

图1  硫酸盐侵蚀混凝土的初始和边界条件

Fig.1  Initial and boundary conditions for sulfate attack concrete

cx,y,0=0 , x,yLcx,0,t=cs , 0,y,tcs , L,y,tcs , x,L,tcscs=φ0c0 (3)

式中:t为侵蚀时间;L为模型长度;s为损伤区域长度;cs为混凝土边界处SO42-的初始浓度;c0为溶液中SO42-的初始浓度;φ0为混凝土的初始孔隙率.

建立的硫酸盐传输模型为:

ct=xDx,ycx+yDx,ycy+crtDx,y=φtτtDc    s<x<Ls<y<L Dx,y=Ds=1-dφtτtDc    0xs0yscrt=-kcu1-αt+12βαt+βαte-16kcte-13kct (4)

式中:Dxy)为SO42-在不同扩散区域的扩散系数;crSO42-消耗浓度; φt)为混凝土侵蚀过程中可传输的时变孔隙率;τt)为混凝土侵蚀过程中的时变曲折度;Dc溶液SO42-扩散系数;Ds为损伤区域中SO42-扩散系数;d为混凝土的损伤程度;k为硫酸盐侵蚀混凝土的化学反应速率;u为混凝土中C3A的浓度;αt)为t时刻水泥的水化程度;β为水泥中石膏的掺量(质量分数).

2.2.1 可传输的时变孔隙率

侵蚀前期反应产物主要在混凝土毛细孔内生长,导致混凝土毛细孔被填充,引起混凝土孔隙率降低.当产物生长到一定程度时造成混凝土膨胀开裂,最终可反映在混凝土的体积应变ε

4.受侵蚀产物填充导致的孔隙率变化量φaε为:

φa=minε,pφ0 (5)
ε=mVC6AS¯3H32-1mV1-1+2mVCS¯H2-1-1f1+mVC6AS¯3H32-1mV2-1+3mVCS¯H2-1-1f2 (6)

式中:p为混凝土内毛细孔的有效填充分数,其取值范围为0.3~0.4,本文取值为0.4;mVC6AS¯3H32mVCS¯H2mV1mV2分别为C6AS¯3H32CS¯H2C4AS¯H12C3A的摩尔体积,其取值分别为0.14×104、1.349×104、0.314×104、1.126×104 mol/m3f1f2分别为C4AS¯H12C3A反应后的体积变化率,分别取0.55和1.31.

混凝土毛细孔隙率会因硫酸盐侵蚀产物不断填充和水泥水化作用而随着侵蚀时间的改变不断变化,毛细孔是硫酸盐传输的主要通

13‑15,其时变孔隙率φt)为:

φt=fcmW/mB-0.36αtmW/mB+0.32+φa (7)
fc=ρcmtρtmc+ρcmsρsmc+ρcρwmW/mB-1 (8)

式中:fc为水泥的体积分数;ρcρtρsρw分别为水泥、石、砂和水的密度,其值分别为1 300、2 600、2 650、1 000 kg/m3mtmsmc分别为石、砂和水泥的质量.

经系统研究,不同水灰比、不同龄期下的硅酸盐水泥水化程度αt)可表示

16

αt=0.697t0.097exp-0.133t0.026mW/mB (9)

2.2.2 时变曲折度

混凝土作为一种多孔介质材料,侵蚀介质在混凝土中传输的速率主要取决于曲折度的大小,在Zuo

15研究的基础上,经修正后的时变曲折度τt)为:

τt=ηsaηst1+0.35fsa1+0.35fstτcpt (10)
τcpt=ωmW/mBηrαt1-1-φt+141-φt1-1-φt (11)
ωmW/mB=1+7αtmW/mB-0.35 (12)

式中:ηsaηst分别为砂、石的形状系数,通过对不同粒径砂、石骨料形貌统计,ηsa取值范围为1.02~1.15,ηst取值范围为1.10~1.42,本文ηsaηst取值分别为1.06、1.28;fsafst分别为砂、石的体积分数,分别取26%和52%;τcpt)为硬化水泥浆体的曲折度;ωmW/mB为水胶比对硬化水泥浆体曲折度的影响修正系数;ηr为水泥颗粒水化后堆积体的形状系数,考虑到不同水灰比硬化浆体产物的形貌特征以及对传输路径的显著影响,对大量硬化水泥浆体的BSE和SEM图像结果进行统计,ηr取值范围一般为1.02~1.42,本文取值为1.20.

2.2.3 损伤程度函数

混凝土受硫酸盐侵蚀时,混凝土的损伤速率与传入混凝土内部的SO42-浓度成正比.设混凝土未损伤时损伤程度为0,完全损伤时的损伤程度为临界损伤程度dc,则损伤程度d为:

dx,y,tt=qccx,y,t= ccΔtc-1cx,y,t , t0ttc0       ,     t<t0,t>tc                                                  (13)

式中:qc为混凝土腐蚀系数,取0.55

16t0tc分别为初始侵蚀时间、完全损伤时间;cc为完全损伤后SO42-恒定浓度.

2.3 求解方法

式(3)、(4)为变系数求解.取传输时间间隔为∆t,将传输总时间划分为kk=t/∆t)等份.各个节点坐标(xi,yj,tk)处的SO42-浓度用ci,jk表示,利用隐格式交替方向的有限差分

17求解偏微分方程,求解过程为:

(1)根据硫酸盐溶液的初始浓度c0,由式(4)计算SO42-在混凝土中的浓度分布;

(2)将SO42-的浓度分布代入式(3),得到由化学反应所消耗的SO42-浓度;

(3)根据化学反应所消耗的SO42-浓度,通过式(1)、(2)得到化学反应所消耗的C4ASH12体积和C3A体积;

(4)通过式(6)计算侵蚀产物生成引起的混凝土体积应变;

(5)通过式(5)、(7)计算混凝土的孔隙率变化量及时变孔隙率,并结合时变曲折度来计算硫酸盐溶液在混凝土中的有效扩散系数;

(6)通过式(4)、(13)计算不同时刻混凝土的损伤程度与损伤后混凝土的有效扩散系数,并通过损伤程度来判断混凝土是否完全损伤.

2.4 结果验证

模型参数为: u=673 mol/m3fc=0.3;α=0.74;β=4.8%;k=3.05×10-8 m3/(mol·s);dc=0.90;Dc=3.5×10-10 m2/s

4.将上述参数代入方程中,通过Matlab迭代求解.

将不同质量分数硫酸钠溶液中混凝土深度x=2 mm处SO42-浓度的模拟结果与试验测试结果进行对比,结果见图2.由图2可见,模拟结果和试验结果基本吻合,最大误差为20.6%.因混凝土组成复杂,其结构也是非均质的,误差在30.0%以内均可认为结果合

18.

图2  混凝土中SO42-的浓度

Fig.2  Concentration of SO42- in concrete

采用Qin

19的试验结果来进一步验证模型的可靠性.其试验中使用P·O 42.5R硅酸盐水泥制备了mW/mB=0.46的混凝土,试件尺寸为100 mm×100 mm×100 mm.在养护28 d后,将试件浸泡在w(Na2SO4)=5%的溶液中,获得了x=1 mm处SO42-在侵蚀时间为60、90、120、150 d的质量分数wSO42-).将试验参数代入到本文模型中,结果见图3(a).由图3(a)可见,模拟结果与试验结果基本吻合,最大误差为18.6%.

图3  硫酸根离子质量分数模拟与试验结果对比

Fig.3  Comparison between modeling and test results of mass fraction of SO42-

Min

20使用P·O 42.5普通硅酸盐水泥制备了mW/mB=0.57的混凝土,试件尺寸为100 mm×100 mm×100 mm,在w(Na2SO4)=5%的溶液浸泡300 d后对x=1 mm处的SO42-质量分数进行测定.将试验参数代入本文模型中,结果见图3(b).由图3(b)可见,试验数据趋势与本文的数值模拟结果较为接近,且相关参数的取值也符合客观规律,最大误差为26.7%,进一步证明了本模型预测结果可靠.

3 影响因素分析

3.1 侵蚀溶液质量分数与侵蚀时间的影响

模拟了mW/mB=0.50的C30混凝土在w(Na2SO4)为1%、3%、5%的溶液中侵蚀90、180、3600 d下SO42-的浓度.模型中相关参数为:试件标准养护28 d;L=100 mm;厚度增量∆x=0.5 mm;dc=0.90.模拟结果见图4.由图4可见:在相同龄期下,随着w(Na2SO4)的增加,混凝土中相同深度处的SO42-浓度增加;在相同浓度下,随着侵蚀龄期的增加,混凝土中相同深度处的SO42-浓度增加,且离子扩散深度增加的速率随着侵蚀时间的增加而变缓.侵蚀180 d,w(Na2SO4)由1%增加到5%时,SO42-扩散深度增加了11.1%;当w(Na2SO4)=5%时,侵蚀180 d比90 d的SO42-扩散深度增加了33.3%.

图4  不同侵蚀时间下SO42-浓度

Fig.4  Concentration of SO42- at different attack time

3.2 临界损伤程度的影响

模拟了mW/mB=0.50的C30混凝土在不同临界损伤程度下的硫酸盐传输-反应-损伤过程.相关参数为:标准养护28 d; L=100 mm;w(Na2SO4)为1%、3%、5%; ∆x=0.5 mm;∆t=1 d;dc为0.80、0.90和0.95.不同质量分数Na2SO4侵蚀下混凝土的剥落厚度见图5.由图5可见:在临界损伤程度相同的条件下,混凝土的剥落厚度变化呈锯齿形状;当Na2SO4溶液质量分数较低时,侵蚀早期未出现损伤剥落,这些规律刚好反映了硫酸盐在混凝土中的传输-反应-填充过程;在临界损伤程度相同的条件下,w(Na2SO4)越高,剥落出现得越早,剥落速率也越快,如当dc=0.95时,混凝土受到w(Na2SO4)=5%的溶液侵蚀时出现损伤剥落的时间比其在w(Na2SO4)=3%的溶液中提早51.8%.

图5  不同质量分数硫酸钠侵蚀下混凝土的剥落厚度

Fig.5  Spalling thickness of concrete under the attack of different w(Na2SO4

侵蚀时间、剥落厚度与损伤程度的关系见图6.由图6可见,随着临界损伤程度的增大,相同深度处混凝土的损伤程度增大,而相同侵蚀时间的剥落厚度减小;当临界损伤程度从0.80增大至0.95时,混凝土3 600 d的剥落厚度减小了42.9%,这表明临界损伤程度的大小对评价混凝土的剥落厚度至关重要.

图6  侵蚀时间、剥落厚度与损伤程度的关系

Fig.6  Relationship between attack time, spalling thickness and damage degree

3.3 水胶比的影响

模拟了mW/mB=0.20~0.60的混凝土相应孔隙率和曲折度变化对硫酸盐传输的影响.相关参数为:混凝土标养28 d;L=100 mm;w(Na2SO4)=3%;∆x=0.5 mm;∆t=1 d;dc=0.90.混凝土试件水胶比与初始孔隙率、曲折度的关系见图7.由图7可见:混凝土的初始孔隙率随着水胶比的增大而增大,且孔隙率增长速率逐渐减小;曲折度与砂石的形状系数和体积率有关(见式(10)~(12)),当砂石的形状系数和体积分数一定时,曲折度随着水胶比的增大而减小,当mW/mB<0.30时,水胶比对曲折度的影响较为显著,反之,水胶比对曲折度的影响较小;砂石对曲折度影响显著,如水胶比由0.20增长至0.60时,硬化水泥浆体的曲折度由13.5降低至4.7,而混凝土的曲折度由25.4降低至9.2.

图7  混凝土水胶比与初始孔隙率、曲折度的关系

Fig. 7  Relationship between mW/mB and initial porosity,tortuosity evolution of concrete

混凝土初始孔隙率与剥落时间的关系见图8.由图8可见,初始孔隙率越大,剥落现象出现得越早,且在相同初始孔隙率条件下,w(Na2SO4)越高出现剥落的时间越早;当初始孔隙率为2.0%时,浸泡在w(Na2SO4)=5%溶液中的混凝土出现剥落现象的时间比在w(Na2SO4)=1%的溶液中提早了约30 a.

图8  混凝土初始孔隙率与剥落时间的关系

Fig.8  Relationship between original porosity and spalling time of concretes

4 结论

(1)基于结晶压理论、体积膨胀理论及Fick第二定律建立了硫酸盐侵蚀下混凝土的传输-反应-损伤剥落模型,在模型中充分考虑了混凝土的孔隙率、曲折度以及混凝土的临界损伤程度等因素,通过试验和文献结果,验证了模型的可靠性,其最大误差为26.7%.

(2)定量分析了不同侵蚀时间和Na2SO4侵蚀溶液质量分数下SO42-的浓度.随着侵蚀溶液质量分数和侵蚀时间的增加,SO42-的扩散深度增加.180 d时Na2SO4溶液质量分数由1%增加到5%时,SO42-的扩散深度增加了11.1%.

(3)定量分析了不同临界损伤程度下,混凝土的剥落厚度与Na2SO4溶液质量分数的关系.在临界损伤程度一定的条件下,Na2SO4溶液质量分数越大,混凝土出现剥落的时间越早,剥落厚度也越大,其在5%Na2SO4溶液中出现剥落现象的时间比在3%Na2SO4溶液中提早了51.8%.

(4)定量分析了水胶比、初始孔隙率、曲折度与混凝土剥落时间的关系.初始孔隙率随着水胶比的增大而增大,而曲折度随着水胶比的增大而减小.随着初始孔隙率的不断增大,剥落时间会出现得更早.

参考文献

1

LI J PXIE FZHAO G Wet al. Experimental and numerical investigation of castinsitu concrete under external sulfate attack and dryingwetting cycles[J]. Construction and Building Materials2020249118789. [百度学术] 

2

关博文陈拴发李华平. 疲劳荷载作用下混凝土硫酸盐腐蚀寿命预测[J]. 建筑材料学报2012153):395398. [百度学术] 

GUAN BowenCHEN ShuanfaLI Huapinget al. Sulfate corrosion life of cement concrete under fatigue load[J]. Journal of Building Materials2012153):395398. (in Chinese) [百度学术] 

3

刘超姚羿舟刘化威. 硫酸盐干湿循环下再生复合微粉混凝土的劣化机理[J]. 建筑材料学报20222511):11281135. [百度学术] 

LIU ChaoYAO YizhouLIU Huaweiet al. Deterioration mechanism of recycled composite powder concrete under drywet cycles of sulfate[J]. Journal of Building Materials20222511):11281135. (in Chinese) [百度学术] 

4

左晓宝孙伟. 硫酸盐侵蚀下的混凝土损伤破坏全过程[J]. 硅酸盐学报2009377):10631067. [百度学术] 

ZUO XiaobaoSUN Wei. Full process analysis of damage and failure of concrete subjected to external sulfate attack[J]. Journal of the Chinese Ceramic Society2009377):10631067. (in Chinese) [百度学术] 

5

ZUO X BSUN WYU C. Numerical investigation on expansive volume strain in concrete subjected to sulfate attack[J]. Construction and Building Materials201236404410. [百度学术] 

6

BARY BLETERRIER NDEVILLE Eet al. Coupled chemotransportmechanical modelling and numerical simulation of external sulfate attack in mortar[J]. Cement and Concrete Composites2014497083. [百度学术] 

7

THAULOW NSAHU S. Mechanism of concrete deterioration due to salt crystallization [J]. Materials Characterization2004532‑4):123127. [百度学术] 

8

FLATT R JSCHERER G W. Thermodynamics of crystallization stresses in DEF [J]. Cement and Concrete Research2008383):325336. [百度学术] 

9

IKUMI TCAVALARO S H PSEGURA Iet al. Alternative methodology to consider damage and expansions in external sulfate attack modeling[J]. Cement and Concrete Research201463105116. [百度学术] 

10

YIN G JZUO X BSUN X Het al. Macromicroscopically numerical analysis on expansion response of hardened cement paste under external sulfate attack[J]. Construction and Building Materials2019207600615. [百度学术] 

11

YIN G JZUO X BTANG Y Jet al. Numerical simulation on timedependent mechanical behavior of concrete under coupled axial loading and sulfate attack[J]. Ocean Engineering2017142115124. [百度学术] 

12

BELLMANN FMÖSER BSTARK J. Influence of sulfate solution concentration on the formation of gypsum in sulfate resistance test specimen[J]. Cement and Concrete Research2006362):358363. [百度学术] 

13

CLIFTON J RBENTZ D PPONNERSHEIM J M. Sulfate diffusion in concrete[M]. GaithersburyNational Institute of Standards and Technology1994913. [百度学术] 

14

HUANG X FZHENG J JZHOU X Z. Simple analytical solution for the chloride diffusivity of cement paste[J]. Science and Technology of Overseas Building Materials2010312):46. [百度学术] 

15

ZUO X BSUN WLIU Z Yet al. Numerical investigation on tortuosity of transport paths in cementbased materials[J]. Computers and Concrete2014133):309323. [百度学术] 

16

中国工程建设标准化协会. 严酷环境混凝土结构耐久性设计标准:T/CECS 1203—2022 [S].北京中国工程建设标准化协会2022. [百度学术] 

China Association for Engineering Construction Standardization. Standard for design of concrete structure durability in severe environments:T/CECS 1203—2022[S]. BeijingChina Association for Engineering Construction Standardization2022.(in Chinese) [百度学术] 

17

AMES W F. Numerical methods for partial differential equations[M]. New YorkAcademic Press2014227. [百度学术] 

18

SUN G WSUN WZHANG Y Set al. Multiscale prediction of the effective chloride diffusion coefficient of concrete [J]. Construction and Building Materials20112510):38203831. [百度学术] 

19

QIN SZOU DLIU Tet al. A chemotransportdamage model for concrete under external sulfate attack[J]. Cement and Concrete Research2020132106048. [百度学术] 

20

MIN H GSUI L LXING Fet al. An effective transport model of sulfate attack in concrete[J]. Construction and Building Materials2019216365378. [百度学术]