湯志杰 廖開星 李 毅 孔祥龍 遆文新
(蘇州熱工研究院有限公司壽命中心,蘇州 215004)
Mehta教授[1]曾在《混凝土耐久性——五十年進展》主旨報告中指出:“當今世界混凝土破壞原因,按重要性遞減順序排列是:鋼筋腐蝕、凍害、物理化學(xué)作用?!倍谝痄摻钿P蝕的因素之中,氯離子侵蝕作用最為顯著。由于國內(nèi)大部分核電站處于高溫高濕的濱海地區(qū),海風和霧氣中含有的氯離子極易沉積在構(gòu)筑物鋼筋混凝土表面,并向其內(nèi)部滲透,引起鋼筋的銹蝕,極大地危害了構(gòu)筑物的可靠性能和安全性能。國內(nèi)外許多專家和學(xué)者對基于氯離子滲透的鋼筋混凝土壽命預(yù)測模型進行了大量的研究,但提出的模型大都為確定型[2],即求解時只要對其中的每個參數(shù)輸入一個數(shù)值,如平均值,即能求出具體的壽命值。但在氯離子侵蝕過程中,涉及許多材料與環(huán)境變量,這些因素都是隨機性的變量。因此,用概率的思想,用更成熟的可靠性理論來研究混凝土結(jié)構(gòu)的氯離子侵蝕問題是更為合理的。
本文通過對某核電站安全相關(guān)構(gòu)筑物現(xiàn)場檢測獲取的數(shù)據(jù),結(jié)合模擬試驗獲得的參數(shù),建立了核電廠安全相關(guān)構(gòu)筑物的壽命預(yù)測模型,并采用Monte-Carlo模擬,對氯離子侵蝕進行了評估。
氯離子在鋼筋混凝土中的擴散符合Fick’s第二擴散定律[3],其表達式為
(1)
初始條件可以表示為
C(x,0)=C,x>0
(2)
邊界條件可以表示為
C(0,t)=CS,t>0
(3)
C(∞,t)=C0,t>0
(4)
由于混凝土內(nèi)部不斷發(fā)生水化反應(yīng),氯離子擴散系數(shù)也會隨時間逐漸變小,因此氯離子擴散模型也應(yīng)考慮時間因素,且同時引入環(huán)境系數(shù)K,該值主要與環(huán)境溫濕度有關(guān),滿足下列關(guān)系式:
(5)
D=KDt
(6)
由式(1)-式(6)得出修正后的Fick’s第二擴散定律的解析解為
(7)
式(7)中,如果氯離子的臨界濃度為Ccr,混凝土保護層厚度為d,則可以推出氯離子侵蝕壽命預(yù)測公式為
(8)
考慮到氯離子侵蝕誘發(fā)的結(jié)構(gòu)耐久性失效的關(guān)鍵點是鋼筋開始銹蝕,本文就將其作為耐久性失效的主要模式。這樣,結(jié)構(gòu)“失效”的狀態(tài)可以定義為在混凝土保護層厚度處,氯離子濃度達到引起鋼筋銹蝕的臨界氯離子濃度Ccr,因此抵抗“破壞”的抗力R可以定義為鋼筋混凝土保護層厚度d,而環(huán)境荷載S定義為混凝土中氯離子的侵蝕深度dc。則氯離子侵蝕環(huán)境正常使用耐久極限狀態(tài)方程可寫為
Z=R-S=d-dc
(9)
(10)
則失效概率Pf為
Pf=P(Z=R-S≤0)=
(11)
用失效概率Pf度量結(jié)構(gòu)的可靠性具有明顯的物理意義,但這種方法最主要的缺點就是,當變量較多時,需要計算多維積分,比較困難,因而可以采納由Siemes和Rostam[4]提出的可靠性指標β來代替Pf度量結(jié)構(gòu)的可靠性??煽啃灾笜甩驴砂聪率蕉x:
(12)
μz=μR-μS
(13)
(14)
式中,P目標為由可靠性指標β所決定的目標失效概率,Pf達到P目標即可判定結(jié)構(gòu)失效;μz為結(jié)構(gòu)函數(shù)Z的均值;σz為結(jié)構(gòu)函數(shù)Z的標準差;μR為抗力R的均值;σR為抗力R的標準差;μs為環(huán)境載荷S的均值;σs為環(huán)境載荷S的標準差。
由于鋼筋混凝土保護層厚度d和氯離子擴散系數(shù)D0是影響鋼筋銹蝕的最重要參數(shù),對混凝土結(jié)構(gòu)的耐久性極大,因此本文假定這兩個參數(shù)為隨機量,其他參數(shù)根據(jù)實際情況取為定值。從而隨機變量dc的均值μdc和標準差σdc可由下式表示:
(15)
(16)
β值可由下式表示:
(17)
處于濱海大氣環(huán)境下的某核電站安全相關(guān)構(gòu)筑物,其混凝土強度等級為C50,表面無涂層保護,已服役20余年。經(jīng)現(xiàn)場檢測,對其鋼筋混凝土保護層厚度d進行統(tǒng)計分析,見圖1,經(jīng)K-S檢驗不否定其服從正態(tài)分布假設(shè),其分布特征為N(58,92),mm。同時對該構(gòu)筑物迎風面墻體取混凝土粉末,測得平均水溶性表面氯離子濃度Cs的值為0.144%(占膠凝材料的質(zhì)量百分比),根據(jù)ACI-318《混凝土結(jié)構(gòu)設(shè)計規(guī)范》[5]的規(guī)定,對于預(yù)應(yīng)力鋼筋混凝土,其臨界水溶性氯離子濃度Ccr應(yīng)取0.06%(占膠凝材料的質(zhì)量百分比)。
圖1 鋼筋保護層厚度頻率直方圖Fig.1 Relative frequency of cover
另外,對該安全相關(guān)構(gòu)筑物的模擬混凝土分別進行了RCM測試和鹽霧試驗[6],獲得臨期為28 d的氯離子擴散系數(shù)D0的分布特征為N(162.91,32.582)×10-12m2/s和擴散系數(shù)衰減指數(shù)m為0.56。結(jié)合該電站的實際環(huán)境,環(huán)境系數(shù)K為0.448。具體相關(guān)影響參數(shù)取值見表1。
表1模型參數(shù)的概率分布特征及取值
Table 1 The value of model parameters
利用表1中參數(shù)代入式(8)中進行Monte-Carlo模擬(模擬次數(shù)為50 000次),獲得氯離子侵蝕條件下構(gòu)筑物耐久性壽命T的分布特征(圖2)。從中可以看出結(jié)構(gòu)耐久性壽命呈現(xiàn)偏態(tài)分布特征,初銹時間峰值點T為6 700年左右。
圖2 結(jié)構(gòu)耐久性壽命的頻率直方圖Fig.2 Relative frequency of life of concrete structure
同樣利用表1中參數(shù)代入式(11)進行Monte-Carlo模擬(模擬次數(shù)為50 000次),得到結(jié)構(gòu)服役失效概率與服役時間的曲線圖,如圖3所示。從圖3中可以看出,失效概率為Pf隨著服役時間的增長而增大。在鋼筋混凝土耐久性設(shè)計中,一般鋼筋銹蝕概率控制在5%~10%是比較符合工程的實際情況[7-8]。核電站安全相關(guān)構(gòu)筑物一般設(shè)計壽命為60年以上,從圖中可以看出,服役時間為100年時,其失效概率Pf僅為2×10-3%。因此該電站安全相關(guān)構(gòu)筑物能夠確保安全服役100年以上。
氯離子侵蝕耐久性極限狀態(tài)的設(shè)置標準見表2。對于氯離子侵蝕導(dǎo)致的鋼筋開始銹蝕或者銹蝕導(dǎo)致的保護層開裂的狀態(tài)而言,從正常使用耐久性能的失效到結(jié)構(gòu)完全破壞(承載力失效)的時間較長,而荷載作用的周期也較長,因而可以將其所對應(yīng)的β值取得較小,失效概率Pf取得高些[9]。依據(jù)耐久性極限狀態(tài)的設(shè)置標準以及一些調(diào)查結(jié)果,本文將β值取為1.645,作為規(guī)定的耐久性極限最小可靠度指標。
圖3 失效概率與結(jié)構(gòu)服役時間關(guān)系Fig.3 The relationship of Pf and service time
表2氯離子侵蝕耐久性極限狀態(tài)的設(shè)置標準
Table 2 Standard of durability limit state of chloride ingress
采用式(17)可以計算可靠度指標與服役時間的關(guān)系,其曲線圖如圖4所示。從圖4中可以看出,當可靠度指標β取為1.645時,t為2 840年。
圖4 可靠度指標與結(jié)構(gòu)服役時間關(guān)系Fig.4 The relationship of β and service time
本文考慮了時間及環(huán)境因素,修正了Fick’s第二定律的誤差函數(shù)解,并以此建立鋼筋混凝土氯離子侵蝕隨機概率壽命預(yù)測模型;通過Monte-Carlo模擬發(fā)現(xiàn)氯離子侵蝕條件下某核電站安全相關(guān)構(gòu)筑物在服役100年后的失效概率Pf僅為2×10-3%,當可靠度指標β取為1.645時,t為2 840年,因此該核電站安全相關(guān)構(gòu)筑物能夠確保安全服役100年以上。實際上,影響氯離子侵蝕的因素還很多,因此在使用本文模型時還需依據(jù)大量實驗數(shù)據(jù)和相關(guān)參數(shù)現(xiàn)場檢測結(jié)果進行取值,這樣才能使預(yù)測結(jié)果更好地服務(wù)于實際工程。
[1] 洪定海.混凝土中鋼筋的腐蝕與保護[M].北京:中國鐵道出版社,1998.
Hong Dinghai.Corrosion and protection of reinforcement in concrete [M].Beijing:China Railway Publishing House,1998.(in Chinese)
[2] 練繼建,張杰,吳慕丹,等.海洋環(huán)境要素綜合影響下混凝土中氯離子的傳輸模型[J].天津大學(xué)學(xué)報(自然科學(xué)與工程技術(shù)版),2016,49(4):333-340.
Lian Jijian,Zhang Jie,Wu Mudan,et al.Chloride ion transport model in concrete under combined influence of marine environment factors[J].Journal of Tianjin University (Science and Technology),2016,49(4):333-340.(in Chinese)
[3] 彭建新,夏偉.預(yù)應(yīng)力混凝土中氯離子擴散模型分析[J].交通科學(xué)與工程,2015,31(3):73-77.
Peng Jianxin,Xia Wei.Diffusion model of chloride ion in prestressed concrete[J].Journal of transport science and engineering,2015,31(3):73-77.(in Chinese)
[4] 岳著文,李鏡培,楊博.考慮對流條件下氯離子在混凝土中的輸運計算[J].同濟大學(xué)學(xué)報(自然科學(xué)版),2015,43(1):60-66.
Yue Zhuwen,Li Jingpei,Yang Bo.Calculation of chloride ions transportation in concrete considering convection[J].Journal of Tongji University (Natural Science),2015,43(1):60-66.(in Chinese)
[5] 高遠,陸春華,袁思奇,等.海工混凝土氯離子分布概率模型分析與應(yīng)用[J].水利水運工程學(xué)報,2016(1):37-43.
Gao Yuan,Lu Chunhua,Yuan Siqi,et al.Application analysis of probability model for chloride ion erosion distribution in marine concrete structure[J].Hydro-Science and Engineering,2016( 1) :37-43.(in Chinese)
[6] 彭建新,王華.氯鹽侵蝕下鋼筋混凝土氯離子擴散系數(shù)試驗及其概率分析[J].中國公路學(xué)報,2014,27(6):77-81.
Peng Jianxin,Wang Hua.Experimental study on chloride ion diffusion coefficient in concrete and probabilistic analysis of its influence on deterioration performance of RC structures[J].China Journal of Highway and Transport,2014,27(6):77-81.(in Chinese)
[7] Collepardi M,Marcialis A,Turrizzani R.The kinetics of penetration of chloride ions into the concrete[J].Cem and Concr Res,1970 (4):157-164.
[8] Tang L,Gulikers J.On the mathematics of time-dependent apparent chloride diffusion coefficient in concrete[J].Cem Concr Res,2007,37:589.
[9] 余紅發(fā),孫偉.混凝土在多重因素作用下的氯離子擴散方程[J].建筑材料學(xué)報,2002,3(5):240-246.
Yu Hongfa,Sun Wei.Diffusion equations of chloride ion in concrete under the combined action of durability factors[J].Journal Of Building Materials,2002,3(5):240-246.(in Chinese)
[10] 趙國藩.工程結(jié)構(gòu)可靠性理論與應(yīng)用[M].大連:大連理工大學(xué)出版社,1996.
Zhao Guofan.Engineering structure reliability theory and application[M].Dalian:Dalian University of Technology Press,1996.(in Chinese)
[11] Siemes A J M,Rostam S.Durable safety and serviceability-A performance based design format[R].IABSE Report 74:Proceedings IABSE Colloquium ‘Basis of Design and Actions on Structures-Background and Application of Eurocode 1,Delft,1996.
[12] 元成方,牛荻濤.基于可靠度理論的海工混凝土結(jié)構(gòu)壽命預(yù)測[J].武漢理工大學(xué)學(xué)報,2013,35(4):73-78.
Yuan Chengfang,Niu Ditao.Life prediction of the marine concrete structure based on the reliability theory[J].Journal of Wuhan University of Technology,2013,35(4):73-78.(in Chinese)
[13] Building Code requirements for structural concrete (ACI 318M-11) and commentary[S].American Concrete Institute,USA,2011.
[14] 劉志勇,孫偉,周新剛.基于Monte Carlo隨機計算的海工混凝土使用壽命預(yù)測[J].東南大學(xué)學(xué)報(自然科學(xué)版),2006,36(增刊II):221-225.
Liu Zhiyong,Sun Wei,Zhou Xingang.Service life prediction of marine concrete based on Monte Carlo stochastic calculation[J].Journal of Southeast University (Natural Science Edition),2006,36(Sup II):221-225.(in Chinese)
[15] 張建強.混凝土氯離子侵蝕的概率分析[J].結(jié)構(gòu)工程師,2009,25(4):48-52.
Zhang Jianqiang.Probabilistic analysis for chloride penetration inside concrete[J].Structural Engineers,2009,25(4):48-52.(in Chinese)
[16] 湯志杰,廖開星.濱海環(huán)境安全殼混凝土氯離子侵蝕規(guī)律研究[J].混凝土,2015(5):66-69.
Tang Zhijie,Liao Kaixing.Diffusion of chloride ion into containment concrete under coastal environment[J].Concrete,2015(5):66-69.(in Chinese)
[17] 劉秉京.混凝土結(jié)構(gòu)耐久性設(shè)計[M].北京:人民交通出版社,2007.
Liu Binjing.Durability design of concrete structure[M].Beijing:China Communications Press,2007.(in Chinese)
[18] 中國土木工程學(xué)會.CCES 01—2004 混凝土結(jié)構(gòu)耐久性設(shè)計與施工指南[S].2004.
China Civil Engineering Society.CCES 01—2004 Guide to durability design and construction of concrete structures[S].China’s Academy of Civil Engineering,2004.(in Chinese)
[19] 馬亞麗.基于可靠性分析的鋼筋混凝土結(jié)構(gòu)耐久壽命預(yù)測[D].北京:北京工業(yè)大學(xué),2006.
Ma Yali.Reliability based durable life prediction of reinforced concrete structures[D].Beijing:Beijing University of Technology,2006.(in Chinese)