基于Weibull分布的鋼筋混凝土通電加速可靠性評(píng)估

路承功, 魏智強(qiáng), 喬宏霞,3, 李 刊, 喬國斌

(1.蘭州理工大學(xué) 土木工程學(xué)院, 甘肅 蘭州 730050; 2.蘭州理工大學(xué) 理學(xué)院, 甘肅 蘭州 730050; 3.蘭州理工大學(xué) 西部土木工程防災(zāi)減災(zāi)教育部工程研究中心, 甘肅 蘭州 730050)

服役于實(shí)際環(huán)境中的鋼筋混凝土結(jié)構(gòu)由于受到氣候和環(huán)境的綜合作用,使得其耐久性較早出現(xiàn)劣化,尤其對(duì)于東部沿海海工混凝土和中國西北內(nèi)陸鹽湖地區(qū)的建筑物,往往達(dá)不到設(shè)計(jì)使用壽命就過早失效,不僅造成巨大的經(jīng)濟(jì)損失,而且使其在使用過程中的可靠性降低[1-5].因此針對(duì)鋼筋混凝土在實(shí)際服役環(huán)境中受到的耐久性劣化因素,設(shè)計(jì)更貼近服役環(huán)境的耐久性試驗(yàn),并尋求合適的可靠性評(píng)價(jià)方法對(duì)于提升鋼筋混凝土結(jié)構(gòu)的安全性具有重要意義[6].

為了探究混凝土在海洋潮汐區(qū)受到氯鹽侵蝕作用的影響,韓學(xué)強(qiáng)等[7]設(shè)計(jì)室內(nèi)干濕循環(huán)試驗(yàn),通過對(duì)不同循環(huán)次數(shù)下混凝土孔隙結(jié)構(gòu)和內(nèi)部氯離子濃度檢測(cè),結(jié)合微觀手段,對(duì)腐蝕產(chǎn)物、表面形態(tài)進(jìn)行了分析,揭示了潮汐區(qū)混凝土腐蝕劣化機(jī)理.針對(duì)混凝土結(jié)構(gòu)在氯鹽侵蝕下的鋼筋銹蝕問題,彭建新等[8]研究者通過氯鹽噴霧試驗(yàn)?zāi)M氯鹽侵蝕環(huán)境,指出了氯離子侵蝕擴(kuò)散規(guī)律,建立了多因素氯離子擴(kuò)散模型,并在考慮空間變異性的基礎(chǔ)上得到了時(shí)空可靠度模型.施錦杰等[9-12]研究者則從鋼筋銹蝕模型、銹蝕速率、混凝土膨脹開裂、服役壽命等方面進(jìn)行了分析研究,使實(shí)際服役環(huán)境下的混凝土結(jié)構(gòu)劣化規(guī)律更加清晰.除了腐蝕性離子,自然界中的電流對(duì)混凝土結(jié)構(gòu)耐久性的影響也較大,當(dāng)電流流經(jīng)鋼筋混凝土結(jié)構(gòu)時(shí),鋼筋銹蝕速率和銹脹量明顯增大,Tang等[13-14]研究者主要從電流對(duì)混凝土結(jié)構(gòu)的損傷破壞方面開展研究,豐富了鋼筋混凝土耐久性的相關(guān)理論.

綜上所述,影響鋼筋混凝土耐久性的因素較多,所開展的耐久性試驗(yàn)也較多,但針對(duì)鋼筋混凝土在腐蝕離子和電流耦合作用下的耐久性研究較少,尤其對(duì)于西部鹽漬土地區(qū),實(shí)際服役環(huán)境中腐蝕離子種類較多,電流的存在極大加劇了鋼筋混凝土耐久性能的劣化.因此研究鹽漬土環(huán)境中鋼筋混凝土在電流作用下的耐久性變化規(guī)律,并對(duì)其可靠性進(jìn)行評(píng)估,對(duì)西部地區(qū)基礎(chǔ)建設(shè)具有重要作用.

1 試驗(yàn)過程及方案設(shè)計(jì)

自然界中雜散電流的存在,對(duì)鋼筋混凝土結(jié)構(gòu)造成了極大的破壞,為了更加客觀合理反映電流與腐蝕離子共同作用下鋼筋混凝土的劣化規(guī)律,本文設(shè)計(jì)如圖1所示的鹽漬土通電加速腐蝕系統(tǒng).鋼筋混凝土試件尺寸為100mm×100mm×400mm,采用直徑8mm 的HRB335螺紋鋼筋縱向貫穿于試件中間,且在試件縱向的底端鋼筋保護(hù)層約為50mm,縱向的頂端裸露約10mm.同時(shí)為了便于后期測(cè)試及避免裸露鋼筋養(yǎng)護(hù)時(shí)銹蝕,在鋼筋裸露處用銅線引出,并用環(huán)氧樹脂進(jìn)行包裹.試件周圍用從格爾木鹽漬土地區(qū)運(yùn)回的鹽漬土填埋,鹽漬土各類離子含量如表1所示.定期噴灑與表1中離子及其含量均相同的鹽溶液使其相對(duì)濕度保持在50%,以保證腐蝕離子具有較好的傳輸媒介.通過直流電源模擬電流輸出,鋼筋作為陽極,碳棒作為陰極.混凝土配合比與基本性能如表2所示,其中水泥為祁連山水泥廠生產(chǎn)的P·O 42.5水泥,粉煤灰為蘭州二熱廠提供的Ⅱ級(jí)粉煤灰,粗骨料及細(xì)骨料均由蘭州本地商砼站提供,級(jí)配良好,含泥量及相關(guān)性能指標(biāo)均符合要求.

圖1 鹽漬土通電加速腐蝕系統(tǒng)Fig.1 Electrically accelerated corrosion system of saline soil

表1 鹽漬土各類離子含量

表2 混凝土配合比與基本性能

參考文獻(xiàn)[15],采用恒電流加速法進(jìn)行通電加速,將電流密度設(shè)定為200μA/cm2,對(duì)應(yīng)的恒電流為20mA,其鋼筋理論質(zhì)量損失與通電時(shí)長的關(guān)系如式(1)所示,并在鋼筋理論質(zhì)量損失率分別為3%、6%、9%、12%及15%,即通電時(shí)長分別為220、440、660、880、1100h時(shí)進(jìn)行電化學(xué)測(cè)試,測(cè)試系統(tǒng)選用電化學(xué)經(jīng)典三電極體系.為避免電流損失,盛裝箱選用塑料制品.

(1)

式中：te為通電時(shí)長,單位s;Δm為理論鋼筋質(zhì)量損失,g;Z為反應(yīng)電極的化學(xué)價(jià),取值+2;F為法拉第常數(shù),取值96500C/mol;M為鐵的相對(duì)原子質(zhì)量,取值56;I為電流,A.

2 耐久性試驗(yàn)結(jié)果分析

不同通電時(shí)刻下鋼筋混凝土極化曲線如圖2所示.其中：E為腐蝕電位；i為腐蝕電流密度.由圖2可見：初始時(shí)刻,鋼筋混凝土的陰陽兩極極化曲線不對(duì)稱,表明陽極鋼筋的溶解過程阻力較大,鋼筋處于鈍化狀態(tài).通電開始后,電流持續(xù)通過,鋼筋表面鈍化膜在短時(shí)間內(nèi)出現(xiàn)破損,環(huán)境中的腐蝕性離子,特別是氯離子與鋼筋直接接觸,并在氧氣和水分同時(shí)存在時(shí)發(fā)生氧化還原反應(yīng),鋼筋作為陽極不斷被消耗,宏觀表現(xiàn)為陰陽兩極的極化曲線逐漸變得對(duì)稱,去極化過程存在的阻力大大減小.同時(shí)由于持續(xù)性電流的通過,在鋼筋附近也會(huì)形成磁場(chǎng),促使鹽漬土中的腐蝕性離子通過混凝土孔隙到達(dá)鋼筋表面,在增加混凝土自身導(dǎo)電性的同時(shí),也會(huì)使得混凝土內(nèi)部pH值大大降低,從而破壞鋼筋穩(wěn)定存在的堿性環(huán)境.

圖2 通電環(huán)境下鋼筋混凝土極化曲線Fig.2 Polarization curves of reinforced concrete in electrified environment

3 Weibull分布可靠性模型

可靠性是“產(chǎn)品”的重要屬性之一,用以衡量其在規(guī)定服役環(huán)境和時(shí)間內(nèi)完成預(yù)定功能的能力.Weibull分布是可靠性工程中一種較為常見的分布,該分布模型靈活性好,能夠通過小樣本得到比較精確的可靠度,常用以分析產(chǎn)品在加速壽命試驗(yàn)中性能在時(shí)間尺度上的變化規(guī)律[16-17].考慮到模型參數(shù)估計(jì)的難度,本文采用兩參數(shù)Weibull分布模型進(jìn)行可靠性分析.Weibull分布模型有2個(gè)參數(shù),即形狀參數(shù)U與尺度參數(shù)V,其中形狀參數(shù)主要影響分布曲線的形狀,而尺度參數(shù)主要影響曲線在橫坐標(biāo)軸上的離散程度.假設(shè)鋼筋混凝土在通電加速環(huán)境中的壽命即失效時(shí)間t服從Weibull分布,則其可靠度R(t)、密度f(t)及失效率λ(t)分別通過式2～4計(jì)算：

(2)

(3)

(4)

3.1 Weibull退化分布過程建立

對(duì)表征鋼筋混凝土腐蝕狀況的極化曲線進(jìn)行活化處理,得到鋼筋混凝土在220、440、660、880、1100h通電時(shí)長下的腐蝕電流密度,并通過指數(shù)函數(shù)、冪函數(shù)及二次函數(shù)建立腐蝕電流密度與通電時(shí)長的關(guān)系.參考文獻(xiàn)[18],當(dāng)腐蝕電流密度達(dá)到10μA/cm2時(shí)反推得到鋼筋混凝土在鹽漬土通電環(huán)境下的失效時(shí)間t,以此作為可靠性分析的失效數(shù)據(jù)進(jìn)行Weibull參數(shù)估計(jì),再將得到的參數(shù)估計(jì)值代入式(2)～(4)得到不同強(qiáng)度鋼筋混凝土通電加速環(huán)境下的可靠度R(t)、密度f(t)及失效率λ(t).

3.2 參數(shù)估計(jì)

Weibull分布參數(shù)估計(jì)方法較多,根據(jù)失效數(shù)據(jù)和Weibull分布函數(shù)的特點(diǎn),并對(duì)比不同參數(shù)方法對(duì)鋼筋混凝土結(jié)果可靠性的影響,本文分別采用中位秩法(MRM)、查表法(LTM)和極大似然法(MLM)進(jìn)行參數(shù)估計(jì)[19].

3.2.1中位秩法

假定有n個(gè)鋼筋混凝土失效時(shí)間,并將任意失效時(shí)間記為ti,得到兩參數(shù)Weibull分布函數(shù)F(ti)為：

(5)

(6)

試驗(yàn)中每個(gè)配合比試件的個(gè)數(shù)為9,為更加準(zhǔn)確反映鋼筋混凝土試件在鹽漬土通電加速環(huán)境中的腐蝕劣化規(guī)律,擴(kuò)大樣本容量,退化數(shù)據(jù)分別采用指數(shù)函數(shù)、冪函數(shù)及二次函數(shù)進(jìn)行回歸,將失效閾值I代入回歸函數(shù)便得到鋼筋混凝土試件在鹽漬土通電加速環(huán)境下的失效時(shí)間t,如表3所示,表中a1、b1、a2、b2、c2、a3、b3、c3分別為回歸函數(shù)的系數(shù).

表3 鋼筋混凝土試件在鹽漬土通電加速環(huán)境下的失效時(shí)間

將表3中的失效時(shí)間進(jìn)行升序排序編號(hào),并利用公式xi=lnti得到xi值;而F(ti)近似等于中位秩MR(ti),因而可得到y(tǒng)i的值.在統(tǒng)計(jì)數(shù)學(xué)理論中,中位秩計(jì)算方法主要有4種,如下式所示:

(7)

式中：Ni為排序后腐蝕天數(shù)編號(hào)；n為失效數(shù)據(jù)樣本容量,本文中n=27.

在xi和yi已知的情況下,通過式(6)可以得到Weibull分布模型的形狀和尺度參數(shù)U、V,如表4所示.由表4可見：對(duì)于3種強(qiáng)度等級(jí)的鋼筋混凝土,采用4種中位秩法計(jì)算參數(shù)時(shí)的相關(guān)性都較好,但是得到的參數(shù)值略有差異,其中第3類中位秩形狀參數(shù)估計(jì)值最大,第4類方法次之,第1類方法形狀參數(shù)估計(jì)值最小;且2類參數(shù)中,形狀參數(shù)的值變化較大,而尺度參數(shù)的值變化較小.

表4 中位秩法參數(shù)估計(jì)匯總表

3.2.2查表法

中位秩與F(ti)之間屬于近似代替,這在一定程度上對(duì)計(jì)算結(jié)果有所影響.為了避免這種誤差,通過查表法對(duì)中間參數(shù)μ和σ進(jìn)行估計(jì).由于樣本容量為27>25,因此采用查表法中的簡單線性無偏估計(jì)方法對(duì)中間參數(shù)μ和σ進(jìn)行估計(jì)[20],結(jié)果如表5所示.

表5 簡單線性無偏估計(jì)參數(shù)匯總表

3.2.3極大似然法

為了判斷中位秩法和查表法所得參數(shù)的準(zhǔn)確性,本文還從Weibull分布函數(shù)的密度函數(shù)(式(3))出發(fā),構(gòu)造似然函數(shù),完全從數(shù)值計(jì)算角度來求解參數(shù)值,即最大似然法.

對(duì)于樣本容量為n的Weibull分布,形狀參數(shù)U與尺度參數(shù)V的似然函數(shù)L為：

(8)

對(duì)式(8)兩邊取對(duì)數(shù),得到對(duì)數(shù)似然函數(shù),并對(duì)U和V分別求偏導(dǎo),則有：

(9)

(10)

令式(9)、(10)結(jié)果為零,得到2個(gè)似然方程：

(11)

(12)

式(11)、(12)是參數(shù)U的超越方程,直接求解較為困難,因此需要采用數(shù)值求解方法求得近似解,為對(duì)比不同數(shù)值求解方法的收斂速度及求解精度,本文選用對(duì)分法(bisection method)、簡單迭代法(simple iteration method)、牛頓-割線法(Newton-Secant method)及斯蒂芬森迭代法(Stephenson iteration method),并利用Matlab程序?qū)崿F(xiàn)對(duì)上述方程的參數(shù)估計(jì).

為充分說明各數(shù)值求解方法的優(yōu)劣,令迭代初始值和控制迭代結(jié)束的精度水平均相同,極大似然方法中不同數(shù)值求解方法的具體迭代次數(shù)、計(jì)算結(jié)束時(shí)的精度及參數(shù)估計(jì)值如表6所示.由表6可見,極大似然法中所得參數(shù)值非常相近,4種數(shù)值求解方法中以斯蒂芬森迭代法的迭代效率最高,精度水平最大,牛頓-割線法迭代效率緊隨其后,而簡單迭代法的迭代效率最差,尤其對(duì)于C40鋼筋混凝土試件,其斯蒂芬森迭代次數(shù)僅是簡單迭代法次數(shù)的8.89%.

表6 極大似然法參數(shù)估計(jì)匯總表

與其他3類迭代方法相比,對(duì)分法雖然迭代原理簡單,易于實(shí)現(xiàn)Matlab編程,但是效率較差;簡單迭代法則主要從方程本身出發(fā)求解變量,易于理解,是4種方法中迭代效率最低的方法;牛頓-割線法利用導(dǎo)數(shù)思想,使得參數(shù)估計(jì)值的收斂速度大于線性收斂速度,同時(shí)利用增量比來代替導(dǎo)數(shù),也避免了求解導(dǎo)數(shù)的困難,是一種簡單高效的數(shù)值計(jì)算方法.整體來看,對(duì)分法和簡單迭代法原理簡單,易于實(shí)現(xiàn)編程,但是迭代次數(shù)較多;而牛頓-割線法和斯蒂芬森迭代法的迭代次數(shù)明顯減少,但需引入導(dǎo)數(shù)及中間參量,提高了Matlab編程難度.在程序可實(shí)現(xiàn)的情況下,后2種方法大大提高了運(yùn)算效率.

綜上,對(duì)比Weibull分布的3種估計(jì)方法所得參數(shù)值可知：4種中位秩法所得形狀參數(shù)估計(jì)值差異較大,而極大似然法的4種數(shù)值計(jì)算方法所得形狀參數(shù)值非常相近,表明極大似然法在Weibull退化分布參數(shù)估計(jì)中更具穩(wěn)定性;并且中位秩法所得形狀參數(shù)值偏大,而查表法所得形狀參數(shù)值偏小,尤其對(duì)于較高強(qiáng)度的C40、C45試件,其偏差程度更大.這主要是因?yàn)橹形恢确安楸矸ㄖ饕獜臄?shù)理統(tǒng)計(jì)方面著手,中位秩及查表所得參數(shù)是在大量統(tǒng)計(jì)基礎(chǔ)上得到的,是一種近似解;而極大似然法完全從數(shù)值角度出發(fā),所需數(shù)據(jù)均是真實(shí)的樣本數(shù)據(jù).此外,不同參數(shù)估計(jì)方法中,形狀參數(shù)變化更為敏感,而尺度參數(shù)變化較小,這與參數(shù)估計(jì)過程及分布函數(shù)的移項(xiàng)等價(jià)有關(guān).在整個(gè)參數(shù)估計(jì)方法中,首先得到形狀參數(shù),而形狀參數(shù)的估計(jì)值與多項(xiàng)指標(biāo)及參數(shù)有關(guān);在得到形狀參數(shù)值后可直接根據(jù)相關(guān)公式計(jì)算出尺度參數(shù),尺度參數(shù)與其他參數(shù)取值無關(guān),因而尺度參數(shù)值波動(dòng)較小.

4 可靠度結(jié)果分析

為形象反映各參數(shù)估計(jì)方法及混凝土強(qiáng)度等級(jí)對(duì)鋼筋混凝土在鹽漬土通電加速環(huán)境中可靠性的影響,將各方法所得參數(shù)值取平均值后分別代入式(2)～(4),得到可靠度R(t)、密度f(t)及失效率λ(t),曲線如圖3～5所示.

圖3 不同參數(shù)估計(jì)方法的可靠度曲線Fig.3 Reliability curves of different parameter estimation methods

圖4 不同參數(shù)估計(jì)方法的密度曲線Fig.4 Density curves of different parameter estimation methods

圖5 不同參數(shù)估計(jì)方法的失效率曲線Fig.5 Failure ratio curves of different parameter estimation methods

由圖3可見：鹽漬土通電加速環(huán)境下不同強(qiáng)度等級(jí)鋼筋混凝土試件的可靠度曲線均表現(xiàn)出明顯的3個(gè)階段,即可靠度保持為1基本不變,可靠度加速下降與可靠度完全為零;且第1階段持續(xù)時(shí)間長,第2階段持續(xù)時(shí)間短.這表明鹽漬土通電加速環(huán)境下鋼筋混凝土試件在較長的時(shí)間內(nèi)可保持較高的可靠度,該過程中腐蝕離子不斷向混凝土內(nèi)部及鋼筋表面遷移,但受到鋼筋表面鈍化膜和混凝土保護(hù)層的影響,氧氣、水分等腐蝕介質(zhì)需要經(jīng)過一段時(shí)間才能抵達(dá)鋼筋表面發(fā)生氧化還原反應(yīng),因而具有較高的可靠度.當(dāng)鹽漬土中的腐蝕性離子濃度達(dá)到鋼筋腐蝕臨界值且腐蝕介質(zhì)充足時(shí),鋼筋銹蝕速度大大加快,鋼筋混凝土試件可靠度急速下降至零,試件完全破壞.由圖3還可知：與C45試件相比,C35、C40試件在第2階段的持續(xù)時(shí)間明顯減少,這是因?yàn)楦邚?qiáng)度混凝土具有更好的密實(shí)性,使腐蝕離子穿過混凝土到達(dá)鋼筋表面更加困難,大大減緩了鋼筋銹蝕劣化程度.

密度函數(shù)能夠表征產(chǎn)品在t時(shí)刻的單位時(shí)間里發(fā)生失效的概率,從密度函數(shù)曲線(圖4)可以看到：3類試件的失效時(shí)間明顯不同,其中C35試件最先發(fā)生失效,失效時(shí)間為820～1000d,其次是C40試件,失效時(shí)間為1000～1300d,C45試件失效時(shí)間最遲,為1100～2700d;從密度極值來看,C35試件極值最大,C40試件次之,C45試件最小,表明鹽漬土通電加速環(huán)境下C35試件失效速率最大,而C45試件最小,這從3類試件可靠度曲線的第2階段也可以得到證實(shí).整體來看,強(qiáng)度最大的C45試件失效時(shí)間最長,密度極值最低,失效程度較為緩慢.對(duì)比中位秩法和極大似然法,其失效時(shí)間基本相同,但中位秩法的密度極值略高于極大似然法,表明極大似然法可靠度失效程度較中位秩法稍大.

失效率曲線不僅能夠反映失效機(jī)理,而且可以定量描述產(chǎn)品失效率的快慢.由失效率曲線(圖5)可見：不同強(qiáng)度鋼筋混凝土試件的失效率曲線變化相似,均呈現(xiàn)出逐漸增大的趨勢(shì),且初始較長的時(shí)間段內(nèi)失效率增加極為緩慢,當(dāng)接近于可靠度曲線第2階段(即失效開始發(fā)生)時(shí),失效率急劇增加.這表明3類試件的失效機(jī)理相同,均為耗損型失效,在通電初始時(shí)刻失效率低,可靠度降低程度小,隨著通電腐蝕時(shí)間增加,腐蝕性離子不斷向鋼筋表面聚集,鋼筋銹蝕加劇,試件失效率逐漸增大,直至試件破壞.

綜上所述,強(qiáng)度等級(jí)對(duì)鋼筋混凝土在鹽漬土環(huán)境中的抗腐蝕影響顯著,強(qiáng)度等級(jí)越高,相同腐蝕時(shí)間下的可靠性越高.參數(shù)估計(jì)方法對(duì)可靠性曲線影響略有差異,其中中位秩法和極大似然法所得鋼筋混凝土試件在鹽漬土通電加速環(huán)境下的可靠度較查表法要高.

5 結(jié)論

(1)Weibull分布的3種參數(shù)估計(jì)法方法中：中位秩法所得形狀參數(shù)估計(jì)值差異較大,極大似然法參數(shù)估計(jì)值非常相近;且中位秩法所得形狀參數(shù)值偏大,而查表法所得形狀參數(shù)值偏小,極大似然法所得形狀參數(shù)值則介于兩者之間且穩(wěn)定性較高;形狀參數(shù)值受各估計(jì)方法的影響較大,而尺度參數(shù)值波動(dòng)較小.

(2)極大似然法利用分布函數(shù)的密度函數(shù),所得參數(shù)估計(jì)值差異性較小,其中對(duì)分法和簡單迭代法原理簡單,易于實(shí)現(xiàn)Matlab編程,但迭代次數(shù)較多,牛頓-割線法和斯蒂芬森迭代法迭代效率較高,但提高了編程難度.

(3)鹽漬土通電加速環(huán)境下強(qiáng)度對(duì)鋼筋混凝土可靠性影響顯著,C35和C40試件的可靠度曲線類似,第1階段的持續(xù)時(shí)間較C45試件略短,且一旦可靠度開始下降,就迅速降至0,C45試件可靠度第2階段持續(xù)時(shí)間顯著增加.3類試件中,C35試件最先發(fā)生失效,且其密度極值最大,失效時(shí)間最短.隨著強(qiáng)度等級(jí)的增加,鋼筋混凝土密度極值變小,失效時(shí)間延長,失效率減小.