王積靜?李春寶

摘 要：鋼筋混凝土結(jié)構(gòu)在受到氯離子等侵蝕介質(zhì)的作用后，往往會發(fā)生不同程度的鋼筋銹蝕乃至保護(hù)層剝落，造成混凝土結(jié)構(gòu)的提早破壞。本文利用COMSOL多物理場耦合軟件建立耦合腐蝕模型，得到鋼筋非均勻腐蝕分布;基于 ABAQUS軟件中的擴(kuò)展有限元，采用預(yù)先設(shè)置多條裂縫后選擇主裂縫的方法，模擬混凝土保護(hù)層銹脹開裂過程?？紤]到混凝土結(jié)構(gòu)中存在多根鋼筋，著重探討不同鋼筋排布方式對混凝土開裂的影響。結(jié)果表明：腐蝕程度和d/s（鋼筋直徑與鋼筋間距的比值）越大，裂縫更容易相互貫通并貫穿保護(hù)層;當(dāng)腐蝕程度一定時，隨著d/s的減小，混凝土中最大裂縫寬度先增大后減小，其中d/s為0.32的情況下裂縫寬度最大。這表明鋼筋排布方式對裂縫擴(kuò)展和貫通情況造成的影響不可忽略，在混凝土結(jié)構(gòu)耐久性設(shè)計和相關(guān)施工中有必要加以考慮。

關(guān)鍵詞：鋼筋銹蝕;鋼筋排布;混凝土保護(hù)層;裂縫擴(kuò)展

中圖分類號：TU375文獻(xiàn)標(biāo)識碼：文章編號：1006-8023（2018）06-0096-06

Study on the Influence of Different Rebar Arrangements on Concrete Cracking Induced by Corrosion Expansion

WANG Jijing? ?Li Chunbao*

（College of Pipeline and Civil Engineering， China University of Petroleum （East China）， Qingdao 266580）

Abstract： After being subjected to chloride ingress， the reinforced concrete structure often has different degrees of rebar corrosion and even cover spalling， resulting in premature destruction of concrete structures.The coupling corrosion model is established by COMSOL multi-physics coupling software， and the non-uniform corrosion distribution of steel bar is obtained.Based on the extended finite element in ABAQUS software， the methods of selecting main cracks and setting multiple cracks in advance are adopted， to simulate corrosion-induced cracking process of concrete cover.Considering the existence of many rebars in concrete structure， the influence of different rebar arrangements on concrete cracking is discussed. The results show that： the larger the corrosion degree and the ratio of d/s （the ratio of rebar diameter to rebar spacing）， the cracks are easier to get through and penetrate the concrete cover. With a certain corrosion degree， the maximum crack width in concrete increases first and then decreases with the decrease of d/s ratio. And in the case of d/s=0.32， the crack width is the largest.This shows that the effect of reinforcement arrangement on crack propagation and penetration cannot be neglected， and it is necessary to consider it in durability design and construction of concrete structures.

Keywords： Rebar corrosion; rebar arrangement; concrete cover; crack propagation

0 引言

鋼筋銹蝕是引起混凝土保護(hù)層開裂和剝落的重要原因，近年來許多國內(nèi)外學(xué)者針對氯離子侵蝕下混凝土結(jié)構(gòu)的腐蝕問題展開分析和研究[1-5]。大多數(shù)研究都采用數(shù)值簡化模型，只針對單根鋼筋的腐蝕膨脹行為進(jìn)行探討。實(shí)際上，混凝土結(jié)構(gòu)中存在多根鋼筋，而且按照強(qiáng)度和構(gòu)造要求進(jìn)行排布，相同截面面積的鋼筋組合由于鋼筋直徑和間距的變化，很可能導(dǎo)致開裂情況有所差別，目前對這一方面的研究相對較少。本文以氯鹽環(huán)境中的鋼筋混凝土梁為研究對象，將氯離子侵蝕過程、鋼筋電化學(xué)反應(yīng)和后期的銹脹開裂進(jìn)行綜合考慮，分析和討論了不同鋼筋排布方式下的混凝土開裂情況，從而為混凝土結(jié)構(gòu)的耐久性設(shè)計和相關(guān)施工提供了有效建議。

1 數(shù)值模型

1.1 耦合腐蝕模型

海洋環(huán)境中的混凝土結(jié)構(gòu)往往都處在非飽和狀態(tài)，而氯離子的侵蝕是擴(kuò)散和對流等多種復(fù)雜機(jī)制綜合作用的結(jié)果。綜合考慮溫度傳遞、水分運(yùn)輸、氯離子運(yùn)輸及結(jié)合作用[6，12-13]，相應(yīng)的偏微分方程如下：

式中：ρc為混凝土密度;cp為混凝土比熱容;DT為混凝土導(dǎo)熱系數(shù);we為孔隙水體積分?jǐn)?shù);Dh為相對水分?jǐn)U散系數(shù);h為孔隙相對濕度;為自由氯離子濃度;Cfc為自由氯離子濃度;Dac為氯離子表觀擴(kuò)散系數(shù)。

當(dāng)鋼筋表面的氯離子濃度達(dá)到臨界濃度時，鋼筋發(fā)生去鈍化，即銹蝕開始出現(xiàn)。本文采用國外學(xué)者C. Alonso[7]通過極化法測定的臨界氯離子濃度值0.364（氯離子質(zhì)量占混凝土質(zhì)量百分比）。

將混凝土視為導(dǎo)體來模擬其中的電位分布，根據(jù)歐姆定律通過實(shí)驗(yàn)可直接獲得有效電阻率ρ （Ω·m）的數(shù)值大小，本文取為200 Ω·m[14]。在混凝土中電位分布的控制平衡方程為[15-17]：

式中：E1為電位，V;ρ為有效電阻率，Ω·m。

鋼筋陽極主要發(fā)生活化極化（氯離子控制），陰極同時發(fā)生活化極化（氯離子控制）和濃差極化（氧氣控制）。活化區(qū)的陽極電流密度為（不考慮鈍化區(qū)的陽極電流密度）：

活化區(qū)和鈍化區(qū)的陰極電流密度為：

式中：i0Fe、i0O為陽極、陰極交換電流密度，分別為3×10-4 A/m2和1×10-5m A/m2;E0Fe、i0O為陽極、陰極平衡電位，分別為-0.78 V和0.16 V;βFe、βO為陽極、陰極Tafel斜率，分別為0.09V/dec和-0.18 V/dec;iL為極限電流密度。

由法拉第定律可得沿鋼筋周向的鋼筋半徑減小值r1表達(dá)式為[15]：

式中：itotal為總體腐蝕電流密度，A/m2;AFe為鐵原子量，55.85 g/mol;t1為鋼筋腐蝕開始的時間，s;ZFe為陽極反應(yīng)化合價，2;Ps為鋼筋密度，7 800 kg/m3;F1為法拉第常數(shù)。

不同的腐-蝕條件下產(chǎn)生的腐蝕產(chǎn)物差異很大，對應(yīng)的體積膨脹率β取值也各不相同，本文根據(jù)Kumiko Suda[8]的試驗(yàn)研究，取體積膨脹率β = 3。不考慮腐蝕產(chǎn)物和鋼筋之間的壓縮變形，腐蝕層厚度r2（m）表達(dá)式為：

r2（θ，t） = βr1（θ，t）? ?。? ? ? ? ? ? ? ? ? ?（8）

因此沿鋼筋周向腐蝕產(chǎn)物銹脹層厚度ur（m）表達(dá)式為：

ur（θ，t）=（β-1）r1? ? 。? ? ? ? ? ? ? ? ? ? ?（9）

1.2 混凝土開裂模型

為了求解的準(zhǔn)確性并便于動態(tài)追蹤裂縫的發(fā)展過程，本章運(yùn)用ABAQUS軟件中特有的擴(kuò)展有限元程序，通過采用預(yù)先設(shè)置多條裂縫后選擇主裂縫的方法，將鋼筋表面各個點(diǎn)的銹蝕層厚度作為位移荷載離散到ABAQUS軟件中的擴(kuò)展有限元模塊中，進(jìn)行開裂求解。將最大主應(yīng)力準(zhǔn)則作為初始開裂依據(jù)，并采用基于能量的線性開裂演化準(zhǔn)則?；炷翑嗔涯躘9]的計算表達(dá)式為：

式中：GF為混凝土的斷裂能，N/m;fc為混凝土標(biāo)準(zhǔn)試件平均抗壓強(qiáng)度，N/m2;dmax為最大集料粒徑，m;w/c為水灰比，1。

根據(jù)ABAQUS軟件模型中的標(biāo)準(zhǔn)尺度，將整個混凝土模型的尺寸離散到每一個矩形網(wǎng)格之中，從而可以準(zhǔn)確地求得混凝土中的最大裂縫寬度。

2 數(shù)值實(shí)例

2.1 數(shù)值建模

結(jié)合規(guī)范的具體要求，混凝土強(qiáng)度等級取為C30，保護(hù)層c取50 mm，鋼筋直徑表示為d，截面尺寸為b×l，鋼筋間距表示為s。為了便于計算忽略箍筋的影響，鋼筋混凝土梁截面選取如圖1所示（以3根鋼筋的模型為例），其中保護(hù)層上表面暴露在氯離子環(huán)境中。

為了探討鋼筋排布對保護(hù)層開裂的影響，現(xiàn)將混凝土保護(hù)層厚度c、混凝土截面設(shè)為定值，鋼筋排布分4種情況討論，另外為了保證具有相同的受力性能，以3Ф32鋼筋排布為參考標(biāo)準(zhǔn)，其余3種鋼筋排布類型分別為5Ф25、4Ф28和2Ф40，鋼筋間距s分別為42.75、61.3 、216 mm。相應(yīng)地，隨著鋼筋直徑的增大，4種排布類型對應(yīng)的d/s分別為0.585、0.457、0.32、0.185。

2.2 初始條件和邊界條件

本文取截面的一半進(jìn)行建模b×l/2（靠近氯離子侵蝕邊界一側(cè)）。在COMSOL軟件中，混凝土內(nèi)部初始孔隙相對濕度h0? = 0.73，內(nèi)部初始溫度取環(huán)境平均溫度T0 = 285.35 K[6];內(nèi)部初始電位E = 0，內(nèi)部初始氯離子濃度Cc0? = 0[15]。不考慮環(huán)境相對濕度和環(huán)境溫度周期性變化，混凝土表面孔隙相對濕度取為環(huán)境平均相對濕度hb = 0.73，表面溫度取環(huán)境平均溫度Tb? = 285.35 K。根據(jù)Duracrete研究[10]，當(dāng)水灰比為0.5時，混凝土的表面氯離子濃度Ccb為0.675%（氯離子質(zhì)量占混凝土質(zhì)量百分比）。在ABAQUS軟件中，將鋼筋混凝土梁截面的左、下、右3條邊界的垂直位移設(shè)為零，另一邊則設(shè)為自由端。

3 結(jié)果分析和討論

3.1 鋼筋腐蝕分布

如圖2所示，沿鋼筋周向腐蝕產(chǎn)物的分布是不均勻的，與文獻(xiàn)[18]得到的結(jié)果相近。鐵銹首先在靠近保護(hù)層一側(cè)的鋼筋表面累積，隨著暴露時間的增加，侵蝕角度θ和銹蝕層厚度也不斷增加。本文認(rèn)為當(dāng)氯離子侵蝕角度為零時，即沒有發(fā)生去鈍化現(xiàn)象;當(dāng)侵蝕角度為45°時，鋼筋發(fā)生四分之一腐蝕;當(dāng)侵蝕角度為90°時，鋼筋發(fā)生二分之一腐蝕;以此類推，氯離子侵蝕至鋼筋末端時（θ = ±180°），鋼筋發(fā)生全面腐蝕。

3.2 鋼筋排布影響下的銹脹規(guī)律

膨脹的銹蝕產(chǎn)物會使周圍混凝土遭受壓力，一旦最大主應(yīng)力超過混凝土的抗拉強(qiáng)度，就會出現(xiàn)裂縫。4種鋼筋排布類型下混凝土節(jié)點(diǎn)變形情況如圖3 ～圖6所示。

3.2.1 混凝土開裂路徑

在四分之一腐蝕時，混凝土承受位移荷載相對較小，開裂路徑基本沿與水平方向夾角呈45°方向擴(kuò)展，應(yīng)力最大處正對侵蝕最嚴(yán)重方向，此時裂縫未貫通也未穿越混凝土保護(hù)層。

在二分之一腐蝕時，當(dāng)d/s為0.585和0.457時，裂縫能夠達(dá)到貫通狀態(tài)但未穿越保護(hù)層，貫通方向與水平方向夾角呈10°，說明混凝土間粘結(jié)應(yīng)力已經(jīng)不足以抵抗銹脹應(yīng)力造成的開裂;此后隨著鋼筋間距的增大，裂縫之間未貫通也未穿越保護(hù)層，且裂縫長度也明顯減小。

在四分之三腐蝕時，當(dāng)d/s為0.585、0.457和0.32時，裂縫能夠達(dá)到貫通狀態(tài)且可以穿越保護(hù)層，貫通方向與水平方向大致呈13°角，此時由于膨脹位移較大，產(chǎn)生的裂縫寬度和長度對比前面兩種腐蝕狀態(tài)皆有明顯增大，且 d/s為0.32時達(dá)到最大值;當(dāng)d/s為0.185時，裂縫不再貫通也未穿越保護(hù)層，破壞情況已經(jīng)有所減輕。

當(dāng)發(fā)生全面腐蝕時，由于腐蝕程度較大，裂縫都能夠穿越保護(hù)層，但是對于d/s為0.185的情況，由于間距過大，致使裂縫并未達(dá)到貫通狀態(tài)。其他幾種情況下，裂縫貫通方向?yàn)榕c水平方向的夾角呈10°，裂縫長度和寬度均隨著d/s的減小而變大。

3.2.2 混凝土最大裂縫寬度

不同鋼筋排布類型下混凝土中的最大裂縫寬度如圖7所示。由圖7可知，隨著腐蝕程度的增加，混凝土中最大裂縫寬度皆呈現(xiàn)遞增趨勢，這與已有的試驗(yàn)結(jié)果[11，19]相吻合。對于前面3種排布方式而言，由于鋼筋間距較小裂縫間較容易貫通，裂縫寬度受腐蝕程度影響較大，尤其是在二分之一腐蝕到全面腐蝕之間，裂縫寬度會發(fā)生急劇變化。同時當(dāng)腐蝕程度一定時，隨著d/s的減小，最大裂縫寬度呈現(xiàn)出先增加后減小的趨勢，最大裂縫寬度出現(xiàn)在d/s為0.32處;當(dāng)d/s為0.185時，裂縫之間難以達(dá)到貫通狀態(tài)，且整體裂縫寬度相對較小。因此，為了減小裂縫最大寬度，在實(shí)際混凝土設(shè)計與施工中應(yīng)該避免d/s為0.32的鋼筋排布方式。

4 結(jié)論

（1）腐蝕產(chǎn)物沿鋼筋周向分布是不均勻的，距離保護(hù)層表面越近，銹蝕層厚度越大;暴露時間越長，侵蝕角度越大。

（2）腐蝕程度和 d/s的值越大，裂縫更容易相互貫通并貫穿整個保護(hù)層，這也說明保證一定的鋼筋間距，能夠延緩裂縫的擴(kuò)展和貫通。

（3）腐蝕程度越大，混凝土中裂縫寬度越寬。當(dāng)腐蝕程度一定時，隨著d/s的減小，最大裂縫寬度呈現(xiàn)出先增加后減小的趨勢;當(dāng)d/s為0.32時，裂縫寬度最大。

【參 考 文 獻(xiàn)】

[1]VAL D， CHERNIN L. Serviceability reliability of reinforced concrete beams with corroded reinforcement[J]. Journal of Structure Engineering， 2009， 135（8）： 896-905.

[2]O?BOLT J， BALABANI? G， KU?TER M. 3D numerical modelling of steel corrosion in concrete structures[J]. Corrosion Science， 2011， 53（12）： 4166-4177.

[3]趙羽習(xí)， 金偉良. 鋼筋銹蝕導(dǎo)致混凝土構(gòu)件保護(hù)層脹裂的全過程分析[J]. 水利學(xué)報， 2005， 36（8）： 939-945.

ZHAO Y X， JIN W L. Full-range analysis for cover cracking of concrete element due to reinforcement corrosion[J]. Journal of Hydraulic Engineering， 2005， 36（8）： 939-945

[4]陳兆軍，吳偉國.連續(xù)配筋混凝土路面裂縫發(fā)展規(guī)律淺析[J].公路工程，2014，39（4）：289-291.

CHEN Z J， WU W G. Study of the law of development of the continuously reinforced concrete pavement crack[J]. Highway Engineering，2014，39（4）：289-291.

[5]包琪， 許純梅. 混凝土橋梁耐久性環(huán)境因素及其防護(hù)涂料分析[J]. 森林工程， 2015， 31（1）： 130-133.

BAO Q， XU C M.Analysis of durability environmental factors and protective coatings for concrete bridge[J].Forest Engineering， 2015， 31（1）： 130-133.

[6]程旭東， 孫連方， 曹志烽，等. 沿海鋼筋混凝土結(jié)構(gòu)Cl～-侵蝕數(shù)值模擬方法研究[J]. 中國腐蝕與防護(hù)學(xué)報， 2015， 35（2）：144-150.

CHENG X D， SUN L F， CAO Z F， et al. Numerical simulation of chloride ion induced corrosion of reinforced concrete structures in marine environment[J]. Journal of Chinese Society for Corrosion and Protection， 2015， 35（2）：144-150.

[7]ALONSO C， CASTELLOTE M， ANDRADE C. Chloride threshold dependence of pitting potential of reinforcements[J]. Electrochimica Acta， 2002， 47（21）： 3469-3481.

[8]SUDA K， MISRA S， MOTOHASHI K. Corrosion products of reinforcing bars embedded in concrete[J]. Corrosion Science， 1993， 35（5）： 1543-1549.

[9]黃海燕， 張子明. 混凝土斷裂能的統(tǒng)計分析[J]. 水利水電技術(shù)， 2004， 35（9）：95-97.

HUANG H Y， ZHANG Z M. Statistical analysis of fracture energy of concrete[J].Water Resources and Hydropower Engineering， 2004， 35（9）：95-97.

[10]LINDVALL A. Duracrete-probabilistic performance based durability design of concrete structures[C]. Budapest： 2nd International PhD Symposium in Civil Engineering， 1998.

[11]施志業(yè). 混凝土保護(hù)層銹脹開裂試驗(yàn)研究[D]. 杭州：浙江工業(yè)大學(xué)， 2010.

SHI Z Y. Experimental study on expensive cracking of concrete cover[D]. Hangzhou： Zhejiang University of Technology， 2010.

[12]MUTHULINGAM S， RAO B N.Non-uniform time-to-corrosion initiation in steel reinforced concrete under chloride environment[J].Corrosion Science， 2014， 82（5）： 304-315.

[13]程旭東， 徐立， 范燕平， 等. 氯離子侵蝕下異類鋼筋的腐蝕行為[J]. 腐蝕與防護(hù)， 2016， 37 （5）： 407-413.

CHENG X D， XU L， FAN Y P， et al. Corrosion behavior of different bars in chloride ion corrosion[J].Corrosion and Protection， 2016， 37 （5）： 407-413.

[14]O?BOLT J， OR?ANI? F， BALABANI? G. Modeling pull-out resistance of corroded reinforcement in concrete： coupled three-dimensional finite element model[J].Cement and Concrete Composites， 2014， 46（2）： 41-55.

[15]CAO C， CHEUNG M. Non-uniform rust expansion for chloride-induced pitting corrosion in RC structures[J]. Construction and Building Materials， 2014， 51（1）： 75-81.

[16]O?BOLT J， OR?ANI? F， KU?TER M， et al. Modeling damage of concrete caused by corrosion of steel reinforcement[C]. International Conference on Fracture Mechanics of Concrete and Concrete Structures， 2013.

[17]ZHU X， ZI G. A 2D mechano-chemical model for the simulation of reinforcement corrosion and concrete damage[J]. Construction and Building Materials， 2017， 137（4）： 330-344.

[18]CHEN E， LEUNG C K Y. Finite element modeling of concrete cover cracking due to non-uniform steel corrosion[J]. Engineering Fracture Mechanics， 2015， 134（1）： 61-78.

[19]PEDROSA F， ANDRADE C. Corrosion-induced cracking： Effect of different corrosion rates on crack width evolution[J]. Construction and Building Materials， 2017， 133（2）： 525-533.