大風干旱環境下開裂混凝土碳化過程數值分析

歐陽威，謝永江，朱長華

(1.中國鐵道科學研究院研究生部，北京100081;2.中國鐵道科學研究院鐵道建筑研究所，北京100081)

大風干旱環境下開裂混凝土碳化過程數值分析

歐陽威1，謝永江2，朱長華2

(1.中國鐵道科學研究院研究生部，北京100081;2.中國鐵道科學研究院鐵道建筑研究所，北京100081)

為了研究大風干旱地區開裂混凝土結構碳化侵蝕引起的耐久性問題，基于CO2在混凝土中擴散時變規律，通過修正低濕度下碳化規律系數和借助離散裂縫模型，并首次引入了碳化過程風作用影響系數，建立了大風干旱環境下開裂混凝土碳化過程的數值模型，給出了初始和邊界條件。根據熱傳導過程與碳化過程的相似性，利用有限元軟件ANSYS對該模型進行了數值計算，并用在役梁體現場檢測數據驗證了該模型的可靠性。結果表明，通過改變模型中裂縫寬度這一參數的大小，可利用該模型分析大風干旱環境下裂縫寬度對混凝土碳化速率的影響。

大風干旱環境 開裂混凝土 碳化 裂縫寬度 數值分析

處于西北大風干旱環境下的混凝土結構，在混凝土早期硬化時極易開裂。在開裂混凝土結構硬化后的服役期，大氣中的CO2是引起鋼筋脫鈍銹蝕，導致混凝土耐久性失效的主要原因。國內外學者主要采用試驗研究來揭示混凝土碳化規律，但混凝土耐久性試驗研究有著周期長、費用高、采集樣本數據有限的弊端，而基于數值分析方法的碳化過程仿真模擬能較好地彌補傳統方法的不足。

1 CO 2 在混凝土中擴散時變規律

1.1 CO2在混凝土中擴散傳質方程

混凝土碳化包括CO2由外部進入內部的物理擴散過程和CO2在混凝土內部與堿發生化學反應的過程，根據菲克擴散定律與質量守恒定律，可推導出混凝土中CO2擴散傳質方程［1］

式中:C為CO2濃度;t為時間;r為單位體積碳化反應CO2消耗速率;D為CO2在混凝土中擴散系數。

1.2 CO2擴散系數和消耗速率

標準環境為溫度20℃，相對濕度65%，無風。

1.2.1 D及其變化規律

考慮環境溫濕度、碳化程度和齡期等，CO2在混凝土中擴散系數D為

式中:Dc，0為標準環境擴散系數［2］，F1(T)為擴散系數溫度影響函數［3-4］，F2(H)為擴散系數相對濕度影響函數，F3(η)為擴散系數碳化程度影響函數［5-6］，F4(te)為擴散系數養護齡期影響函數［7］，表達式分別為

式(3)～(9)中:fcu，k為混凝土28 d抗壓強度標準值; Ea，c為擴散活化能，取35 kJ/mol;R為摩爾氣體常數; Tref為標準環境溫度;H為相對濕度;η為碳化程度; fcu，28為混凝土28 d抗壓強度值;fcu(te)為混凝土抗壓強度的時間函數;fcu，∞為混凝土完全水化后的抗壓強度;te為等效齡期，K為強度發展速率影響常數。

1.2.2 CO2消耗速率r及其變化規律

考慮環境溫濕度、CO2濃度、碳化程度等，碳化反應CO2消耗速率r為［8］

式中:r0為標準環境下碳化CO2消耗速率，r0=2.0×10－7kg/(s·m3)，f1(T)為消耗速率溫度影響函數［8］，f2(H)為消耗速率相對濕度影響函數，f3(C)為消耗速率濃度影響函數，f4(η)為消耗速率碳化程度影響函數，表達式分別為

式(11)～(16)中:U0為混凝土碳化活化能與氣體常數的比值，取為11.25 k J/mol;Cmax為單位混凝土中使碳化反應充分的CO2濃度最大值;Cr為碳化過程中單位混凝土累積消耗CO2的量，Cr，max為單位混凝土消耗CO2的最大值;Z為單位體積混凝土的水泥用量，可由水灰比求得;普通波特蘭水泥中［CaO］≈64%。

2 大風干旱環境和裂縫對碳化規律的影響

西北地區氣候環境特殊，以蘭新二線為例，表現在3個方面:①大部分年平均相對濕度在20%～60%，濕度低;②途經安西風區等五大風區，風大且多;③冬季嚴寒，夏季日照強烈，季節性溫差明顯。在這種低濕度、大風、大溫差作用下，混凝土極易開裂。

2.1 低相對濕度對混凝土碳化規律的影響

2.1.1 擴散系數相對濕度影響函數

低濕度下仍可采用Saetta公式，見式(5)。

2.1.2 消耗速率相對濕度影響函數

濕度較高時，碳化反應速率較快，濕度較低時，碳化反應速率慢。大風干旱地區平均年相對濕度低于50%，故將消耗速率相對濕度影響函數修正為

2.2 風作用影響系數

風壓作用下，不僅存在由于濃度梯度引起的CO2在混凝土中的擴散作用，還有風壓引起的CO2在混凝土中的滲透作用［9］。引入表征風壓影響下碳化深度的函數——風作用影響系數模型

風作用影響系數f(W)與風壓作用下碳化滲透深度密切相關，滲透深度與風速v、風作用時間τ和風向φ有關。風壓僅對碳化過程中的擴散過程產生影響，在數值上

碳化理論模型［10］中碳化深度表達式為

式中，m0為單位體積混凝土吸收CO2的量。

由式(19)和式(20)，令F5(W)=f(W)，推理得到

令

根據文獻［11］工程經驗的研究成果:受風壓影響面與不受風壓影響面的碳化深度比值γ通常在1～3。本文根據大風干旱地區服役混凝土結構的碳化深度檢測結果，取大風干旱環境下γ=2，則大風干旱環境下CO2在混凝土中的擴散系數Df為

2.3 離散裂縫模型

為體現裂縫對混凝土碳化的影響，本文借助離散裂縫模型，將裂縫離散化。其基本思想是將開裂混凝土碳化分為完好混凝土中的碳化和在裂縫中CO2的傳輸兩部分，如圖1所示。

圖1 離散裂縫模型示意

CO2在裂縫中一維擴散控制方程為

式中:Ccr為t時刻裂縫中x深度的CO2濃度，Dcr為裂縫中CO2擴散系數，Dcr近似于空氣中以主擴散為主的擴散系數

式中:u—A是CO2分子均方根速度，λ為CO2分子的平均自由程;k為玻爾茲曼常數，m為氣體分子質量;d0為分子有效直徑，取為0.33 nm;p為氣體壓強，取為0.1 MPa。

經計算可得常溫常壓下Dz的數量級為10－5m2/s。而CO2在一般混凝土(C30以上)中的擴散系數為10－8m2/s。為簡化計算，本文裂縫擴散系數取為混凝土擴散系數的1 000倍。

3 大風干旱環境開裂混凝土碳化過程數值模型

3.1控制方程

引入裂隙巖體滲流工程和油藏工程中廣泛應用的雙重孔隙介質模型的處理方法，將混凝土和裂縫看作兩相介質，且CO2在該兩種介質中的擴散都符合菲克第二定律，則大風干旱環境下開裂混凝土CO2濃度分布的控制方程如下。

混凝土中CO2的擴散方程為

裂縫中CO2的擴散方程為

式中:C'為裂縫中CO2濃度值，Df，x，Df，y分別為大風干旱環境下混凝土x方向和y方向的CO2擴散系數。

CO2在開裂混凝土中的擴散模型如圖2所示?；炷量醋骶鶆虿牧蠒r，可令Df，x=Df，y=Df，大風干旱環境下裂縫中CO2擴散系數Dcr，f≈Dz。

圖2 CO2在開裂混凝土中的擴散模型

3.2 初始條件和邊界條件

3.2.1 初始條件

假定初始狀態下混凝土內部均勻且無碳化，則初始CO2濃度為

3.2.2 邊界條件

初始狀態及反應過程中，假定接觸邊界CO2濃度與環境中CO2質量濃度C0一致，即接觸邊界CO2濃度為

3.2.3 碳化鋒面深度臨界值

碳化鋒面臨界條件

根據文獻［12］，將混凝土碳化程度pH=9.0作為碳化鋒面臨界條件來判斷混凝土是否碳化。如圖3所示，當pH≤9.0時，混凝土完全碳化，當pH＞9.0時混凝土未碳化。

圖3 碳化深度臨界值判定

文獻［12］還研究了混凝土pH值與碳化程度η的關系。

化學試劑法測試混凝土碳化深度的常用試劑是1%濃度的酚酞酒精溶液，它以pH=9.0為界限，已碳化區呈無色，未碳化區呈粉紅色。

4 算例分析

4.1 問題描述及分析

以西北蘭新二線的環境條件為背景，選擇T1碳化環境下的鐵路在役混凝土結構為研究對象，保護層厚度取35 mm。西北大風干旱地區環境條件選定為:風速10 m/s，風壓作用時間約為50%，氣溫20℃，相對濕度40%，CO2濃度為0.032%。

圖4 開裂混凝土幾何模型(單位:mm)

開裂混凝土幾何模型如圖4所示。模型大小為200 mm×100 mm，裂縫深度L=35 mm，裂縫寬度為w。材料基本參數值見表1?；炷羶瘸跏糃O2為0，僅上表面接觸CO2，其他三面密封。

文獻［13］在青藏、蘭新、包蘭等西北大風干旱環境下鐵路沿線選擇了有代表性的鋼筋混凝土結構進行了碳化深度的檢測(見表2)。本文進行了如下計算:①裂寬w=0時，計算與表2相同參數下的碳化深度值，并與表2中的現場檢測數據進行比較。②計算不同裂縫寬度時C30混凝土碳化到鋼筋表面所經歷的最短時間。

表1 材料基本參數值

表2 在役梁體保護層混凝土碳化深度檢測結果

根據碳化過程與熱傳導過程的相似性，將有限元軟件ANSYS熱分析部分的材料屬性中的密度與比熱參數均設置為1，將材料熱傳導系數替換為CO2擴散系數D;將溫度替換為CO2濃度C;將生熱率函數替換為碳化反應速率－r;模型初始溫度為0℃，上表面施加溫度荷載，其他三面設置為絕緣邊界。通過數值計算，可得CO2質量濃度分布，再通過碳化鋒面臨界值限定條件，即可求得碳化深度。

4.2 有限元模型

混凝土單元選擇二維4節點平面實體單元PLANE55。當w=0.2 mm時，有限元模型共劃分為2 365個單元，2 456個節點。加載步長設為0.5年。

4.3 計算結果分析

4.3.1 CO2質量濃度分布

計算得出的CO2質量濃度分布云圖如圖5所示。可以看出，裂縫區域CO2擴散程度顯著大于無裂縫區域;距離表面越深，碳化程度越小，并且裂縫附近碳化程度稍微大于遠離裂縫位置，距離表面越深，這種現象越明顯。

圖5 CO2質量濃度分布云圖

4.3.2 數值計算值與現場檢測值的比較

為檢驗有限元模型的準確性，將裂縫寬度為0時的在役梁體保護層混凝土碳化深度檢測結果與大風干旱環境下混凝土碳化模型的計算結果進行了對比，見表3。由表3可知，碳化深度數值模擬結果與在役梁體檢測結果吻合程度尚可。

表3 混凝土碳化深度實測值與計算值的比較

圖6 裂縫影響區碳化到鋼筋表面所需時間

4.3.3 裂縫寬度的影響

以0.05 mm為步長，在0～0.5 mm范圍內逐漸增加裂縫寬度，計算在裂縫深度35 mm處(即T2環境鋼筋保護層厚度)，裂縫影響區混凝土碳化到鋼筋表面所需的時間，如圖6所示?？梢钥闯?，隨著裂縫寬度的增加，碳化到鋼筋表面的時間縮短，裂縫寬度對CO2的擴散速度影響顯著，尤其是當裂縫寬度＞0.2 mm時，碳化到鋼筋表面的速度急劇增加。

5 結語

本文針對大風干旱環境下開裂混凝土的碳化過程，綜合考慮溫度、濕度、大風、時間、混凝土強度、CO2濃度等因素對混凝土碳化的影響，借助離散裂縫模型，進行了低濕度環境下混凝土碳化規律修正，首次引入風作用影響系數，建立了大風干旱環境下開裂混凝土碳化過程的數值模型，并用有限元軟件ANSYS對該模型進行了仿真模擬。通過把數值結果與現場在役檢測結果相比較，驗證了該模擬方法的有效性。

［1］王路平.基于ANSYS的混凝土結構碳化過程數值模擬分析［D］.長沙:中南大學，2010.

［2］Federation International de la Precontrainte.CEB-FIP Model Code 1990［S］.London:Thomas Telford，1990.

［3］SAETTA A V，VITALIANI R V.Experimental Investigation and Numerical Modeling of Carbonation PartⅠ:Theoretical Formulation［J］.Cement and Concrete Research，2004，34 (4):571-579.

［4］SAETTA A V，VITALIANI R V.Experimental Investigation and Numerical Modeling of Carbonation PartⅡ:Practical Applications［J］.Cement and Concrete Research，2005，35 (5):958-967.

［5］SAEKI T，OHGA H，NAGATAKI S.Mechanism of Carbonation and Prediction of Carbon/Carbon Composites［J］.Carbon，2003(41):2625-2634.

［6］SONG H W，KWON S J，BYUN K J，et al.Predicting Carbonation in Early-aged Cracked Concrete［J］.Cement and Concrete Research，2006(36):979-989.

［7］SAETTA A V，SCHREFLER B A，VITALIANI.2-D Model for Carbonation and Moisture/heat Flow in Porus Materials［J］.Cement and Concrete Research，1995，25(8):1703-1702.

［8］STEFFENS A，DINKLER D，AHRENS H.Modeling Carbonation for Corrosion Risk Prediction of Concrete Structures［J］.Cement and Concrete Research，2002(32): 935-941.

［9］屈文俊，張譽，車惠民.既有混凝土結構的碳化預測［J］.建筑結構，1999，36(4):31-37.

［10］牛荻濤.混凝土結構耐久性與壽命預測［M］.北京:科技出版社，2003.

［11］郭猛，劉洪，白文靜.混凝土橋梁結構的耐久性設計［J］.建筑科學與工程學報，2010，27(1):117-121.

［12］PAPADAKIS V G，VAYENS C G.Fardis M N.A Reaction Engineering Approach to the Problem of Concrete Carbonation［J］.AIChEJournal，1989，35(10):1639-1950.

［13］王春芬.旱區鐵路混凝土橋梁耐久性及安全性評估［D］.西安:西安建筑科技大學，2011.

Numerical analysis of carbonation process in cracked concrete under w indy and arid environment

OUYANG Wei1，XIE Yongjiang2，ZHU Changhua2
(1.Postgraduate School，China Academy of Railway Sciences，Beijing 100081，China; 2.Railway Engineering Research Institute，China Academy of Railway Sciences，Beijing 100081，China)

T he paper looks in to the carbonation and durability of cracked concrete in w indy ad arid areas，and incorporates in to its analysis CO2proliferation in concrete m aterial，revised carbonation factors under sligh tly hum id environment and discrete crack m odel.It proposes for the first time in the field that w ind-effect coefficients shall be considered.On this basis，the paper builds the num erical m odel for the carbonation of cracked concrete in w indy and arid regions and clarifies the precond itions and boundary cond itions accordingly.G iven that the carbonation process bears sim ilarities w ith heat conduction，it conducts num erical calcu lation w ith the help of finite elem ent softw are ANSYS，whose validity has been verified in the in-situ testing of in-service beam.T he results indicate that crack w id th stands as a relevant facto r and the m odel can be used for the analysis on relation betw een w id th change and carbonation speed.

W indy and arid environm en t;C racked concrete;Carbonation;C rack depth;Num erical analysis

TU528.01

A

10.3969/j.issn.1003-1995.2015.06.42

1003-1995(2015)06-0163-05

(責任審編周彥彥)

2014-12-04;

2015-01-07

中國鐵路總公司科技研究開發計劃項目(J2013G004)

歐陽威(1987—)，男，湖北鄂州人，研實員，碩士研究生。