王繼林 夏葉飛 袁微微 宗周紅
(1.華設(shè)設(shè)計(jì)集團(tuán)股份有限公司 南京 210014; 2.東南大學(xué)土木工程學(xué)院 南京 210096)
為了保障橋梁交通的安全暢通,迫切要求對(duì)既有橋梁全面評(píng)價(jià)全壽命周期技術(shù)性能。剩余使用壽命預(yù)測(cè)研究是一個(gè)復(fù)雜的課題,盡管已有大量學(xué)者開(kāi)展理論實(shí)踐研究,但依然存在許多值得探討和研究的問(wèn)題。如在橋梁體系的可靠度分析中,構(gòu)件層次已經(jīng)作了嘗試研究,但在結(jié)構(gòu)體系可靠性研究方面還需進(jìn)一步深入。國(guó)內(nèi)外眾多學(xué)者對(duì)在役混凝土結(jié)構(gòu)抗力影響因素方面進(jìn)行了大量的理論和試驗(yàn)研究,目前建筑結(jié)構(gòu)方面的研究成果較多,而橋梁結(jié)構(gòu)方面研究成果很少?,F(xiàn)今研究者多從考慮橋梁全壽命期內(nèi)的耐久性能或安全性能來(lái)進(jìn)行橋梁分析研究,較少綜合考慮耐久性能和安全性能。
本文綜合考慮結(jié)構(gòu)服役壽命周期內(nèi)橋梁實(shí)際耐久性能和安全性能,進(jìn)行結(jié)構(gòu)服役全壽命周期性能預(yù)測(cè)評(píng)估[1](life cycle performance assessment,LCPA),先計(jì)算結(jié)構(gòu)功能函數(shù)時(shí)變可靠度,通過(guò)可接受可靠度界限進(jìn)行橋梁剩余使用壽命預(yù)測(cè)研究。橋梁LCPA剩余使用壽命預(yù)測(cè)是在考慮壽命期荷載和耐久性因素對(duì)結(jié)構(gòu)的極限承載力指標(biāo)和正常使用狀態(tài)功能可靠度影響的基礎(chǔ)上,實(shí)現(xiàn)壽命預(yù)測(cè)。
結(jié)構(gòu)全周期壽命是指橋梁從交竣工驗(yàn)收至結(jié)構(gòu)狀況退化到臨界危險(xiǎn)水平的服役時(shí)間,橋梁服役期根據(jù)實(shí)際作用和結(jié)構(gòu)損傷實(shí)時(shí)預(yù)測(cè)結(jié)構(gòu)可靠度可以預(yù)測(cè)橋梁剩余的使用壽命。預(yù)測(cè)的剩余使用壽命與眾多因素有關(guān),如材料性能、結(jié)構(gòu)形式、劣化機(jī)理、維修方式等。通?;炷翗蛄旱牧踊侵T多因素共同作用的結(jié)果,現(xiàn)今還是無(wú)法全面地量化這些影響因素。本文嘗試通過(guò)隨機(jī)過(guò)程的方法綜合考慮結(jié)構(gòu)安全性能與耐久性能,通過(guò)同時(shí)考慮安全性和耐久性來(lái)進(jìn)行橋梁剩余使用壽命的預(yù)測(cè)。
對(duì)結(jié)構(gòu)進(jìn)行LCPA剩余使用壽命預(yù)測(cè),需先進(jìn)行結(jié)構(gòu)各要素時(shí)變可靠度的研究。當(dāng)前各類(lèi)時(shí)變可靠度領(lǐng)域發(fā)展迅速,這些成果為預(yù)測(cè)結(jié)構(gòu)剩余使用壽命理論奠定了堅(jiān)實(shí)的基礎(chǔ)。2011年Frangopol、Tong Guo等針對(duì)現(xiàn)存的預(yù)應(yīng)力混凝土箱梁橋承受不斷增加的交通荷載和處于惡劣環(huán)境下兩方面問(wèn)題,提出了多峰車(chē)載模型及箱梁橋的時(shí)變腐蝕模型,得出針對(duì)預(yù)應(yīng)力混凝土箱梁橋的時(shí)變可靠度分析方法。
本項(xiàng)目以京滬高速(G2)新沂河大橋?yàn)楣こ瘫尘埃瑢?duì)連續(xù)梁橋進(jìn)行LCPA剩余使用壽命的預(yù)測(cè)。主要步驟:①根據(jù)動(dòng)態(tài)稱(chēng)重系統(tǒng)和統(tǒng)計(jì)分析,研究車(chē)輛荷載效應(yīng)的時(shí)變可靠度模型,采用概率分布函數(shù)卷積的計(jì)算方法,研究總荷載效應(yīng)最大值概率分布模型。②根據(jù)材料退化模型計(jì)算新沂河大橋的安全性和耐久性時(shí)變可靠度模型,依據(jù)LCPA對(duì)新沂河大橋進(jìn)行剩余使用壽命預(yù)測(cè)。
新沂河大橋位于京滬(G2)高速公路,跨越新沂河,橋梁上部結(jié)構(gòu)采用30 m跨徑先簡(jiǎn)支后連續(xù)的裝配式部分預(yù)應(yīng)力混凝土連續(xù)箱梁橋。2001年建成通車(chē),2012年橋梁健康監(jiān)測(cè)系統(tǒng)安裝完成開(kāi)始采集數(shù)據(jù)。本文以新沂河大橋狀態(tài)識(shí)別研究為工程背景,立足于長(zhǎng)期性能監(jiān)測(cè)系統(tǒng)及動(dòng)態(tài)稱(chēng)重系統(tǒng),對(duì)連續(xù)梁橋進(jìn)行LCPA剩余使用壽命的預(yù)測(cè)。本研究中主要是預(yù)測(cè)新沂河大橋的功能時(shí)變可靠度,結(jié)合目標(biāo)可靠度指標(biāo)對(duì)新沂河大橋的剩余使用壽命進(jìn)行評(píng)估。主要內(nèi)容如下:①依據(jù)動(dòng)態(tài)稱(chēng)重系統(tǒng)數(shù)據(jù)得出既有橋梁的隨機(jī)車(chē)載模型,結(jié)合恒載概率模型建立既有橋梁的荷載概率模型;②考慮構(gòu)件耐久性參數(shù)隨時(shí)間變化建立既有橋梁構(gòu)件的抗力衰減概率模型;③計(jì)算箱梁抗力與效應(yīng)包絡(luò)時(shí)變可靠度,根據(jù)既有橋梁目標(biāo)可靠度指標(biāo),預(yù)測(cè)評(píng)估既有橋梁的剩余使用壽命。
2.1.1恒載效應(yīng)概率模型
一般認(rèn)為橋梁結(jié)構(gòu)的恒載服從正態(tài)分布,恒載效應(yīng)與恒載之間一般按線(xiàn)性關(guān)系考慮,恒載效應(yīng)與恒載相同的概率模型,均滿(mǎn)足正態(tài)分布,恒載彎矩效應(yīng)概率分布函數(shù)FSG(x)和統(tǒng)計(jì)參數(shù)如式(1)、式(2)。
(1)
μSG=κSGSGK;σSG=μSGδSG
(2)
式中:SGK為單梁恒載彎矩效應(yīng)標(biāo)準(zhǔn)值,SGK=3 434 kN·m;κSG=1.014 8;μSG、σSG、δSG分別為構(gòu)件恒載效應(yīng)的平均值、標(biāo)準(zhǔn)差和變異系數(shù),μSG=3 485 kN·m,σSG=150 kN·m,δSG=0.043 1。
2.1.2車(chē)輛活載效應(yīng)概率模型
新沂河大橋?qū)嶋H運(yùn)行車(chē)輛荷載為一般運(yùn)行狀態(tài)與密集運(yùn)行狀態(tài)的綜合,稱(chēng)為混合運(yùn)行狀態(tài)。按照2種運(yùn)行狀態(tài)下各自車(chē)隊(duì)總數(shù)的比例,將一般與密集運(yùn)行狀態(tài)下的彎矩概率函數(shù)進(jìn)行組合,得到混合運(yùn)行狀態(tài)下的概率密度函數(shù),進(jìn)而得到荷載重現(xiàn)期分別為1,5,10,30,50,70,100年時(shí)的混合狀態(tài)下的彎矩概率分布函數(shù)。下面采用2013年4月份-2014年4月份動(dòng)態(tài)稱(chēng)重系統(tǒng)采集的活載數(shù)據(jù),取重現(xiàn)期1年,計(jì)算混合運(yùn)行狀態(tài)下1年梁體跨中截面活載彎矩效用最大值的概率分布函數(shù)F1年(M)[2],概率分布函數(shù)見(jiàn)圖1。
圖1 1年重現(xiàn)期車(chē)輛活載彎距效應(yīng)概率F1年(M)分布圖
通常假定構(gòu)件抗力概率分布服從對(duì)數(shù)正態(tài)分布,在役橋梁構(gòu)件的彎矩抗力衰減概率模型可按式(3)計(jì)算。
R(t)=KT(t)KP(t)RP[fcd(t),fsd(t),
δel(t),ky(t),kb(t)…]
(3)
式中:R(t)為構(gòu)件抗力隨機(jī)過(guò)程模型;KT(t)為考慮抗力參數(shù)測(cè)試及預(yù)測(cè)影響的隨機(jī)變量;Kp(t)為抗力計(jì)算模式不定性隨機(jī)變量;Rp(·)為規(guī)范中的抗力計(jì)算公式;fcd(t)、fsd(t)為材料強(qiáng)度退化隨機(jī)值;δel(t)、ky(t)、kb(t)為第i根鋼筋銹蝕深度屈服強(qiáng)度協(xié)同工作系數(shù)預(yù)測(cè)值。JTG 3362-2018規(guī)范中,箱梁正截面受彎抗力預(yù)測(cè)概率模型如式(4)。
(4)
具體在推導(dǎo)構(gòu)件抗力預(yù)測(cè)概率模型時(shí)需要綜合考慮混凝土的強(qiáng)度、鋼筋銹蝕深度、鋼筋黏結(jié)性能等耐久因素退化。本項(xiàng)目工程橋梁處于江蘇淮安沭陽(yáng)地區(qū),相應(yīng)參數(shù)時(shí)變模型選擇牛荻濤[3]模型,包括混凝土強(qiáng)度經(jīng)時(shí)變化模型、鋼筋銹蝕率模型等。
1) 一般大氣環(huán)境下混凝土強(qiáng)度劣化模型、混凝土強(qiáng)度經(jīng)時(shí)變化模型,混凝土強(qiáng)度平均值μf(t)和標(biāo)準(zhǔn)差σf(t)的經(jīng)時(shí)模型如式(5)。
(5)
式中:μf0為混凝土強(qiáng)度回彈實(shí)測(cè)推定值,本橋取21.7 MPa;η(t)為混凝土強(qiáng)度平均值隨時(shí)間變化的函數(shù);σf0為混凝土強(qiáng)度回彈實(shí)測(cè)換算值標(biāo)準(zhǔn)差,本橋取1.05 MPa。
2) 鋼筋銹蝕深度(保護(hù)層脹裂前)計(jì)算采用牛荻濤[4]銹蝕深度預(yù)測(cè)模型如式(6)。
(6)
式中:δel(t)為混凝土保護(hù)層的鋼筋銹蝕深度,mm;λel為銹脹開(kāi)裂前的鋼筋銹蝕速度,mm/年;kcr為鋼筋位置修正系數(shù);kce為環(huán)境條件修正系數(shù);RH為環(huán)境相對(duì)濕度,%;d為混凝土保護(hù)層厚度,mm;fcu為混凝土立方體抗壓強(qiáng)度,MPa。
3) 計(jì)算銹蝕鋼筋屈服強(qiáng)度降低系數(shù)ky(t),屈服強(qiáng)度降低系數(shù)的計(jì)算采用惠云玲[5]模型,如式(7)。
(7)
式中:ηs(t)為鋼筋銹損率;δe(t)為鋼筋銹蝕深度,mm;D為鋼筋直徑,mm。
4) 銹蝕鋼筋與混凝土協(xié)同工作系數(shù)kb(t)計(jì)算采用金偉良[6]模型如式(8)。
kp(t)=
(8)
橋梁結(jié)構(gòu)的目標(biāo)可靠指標(biāo)受破壞影響程度、評(píng)估基準(zhǔn)期、構(gòu)件對(duì)結(jié)構(gòu)整體安全的影響程度、結(jié)構(gòu)的重要性等因素影響,目標(biāo)可靠指標(biāo)見(jiàn)表1。根據(jù)實(shí)際結(jié)構(gòu)破壞類(lèi)型特征,本次參照文獻(xiàn)[7-8]將0.85%作為結(jié)構(gòu)達(dá)到壽命終點(diǎn)標(biāo)志,取新沂河大橋目標(biāo)可靠指標(biāo)3.995,對(duì)應(yīng)目標(biāo)失效概率為1.1×10-6。應(yīng)用新沂河大橋荷載效應(yīng)概率模型及抗力衰減概率模型,依據(jù)JC法及Monte-Carlo法計(jì)算新沂河大橋的彎矩時(shí)變可靠指標(biāo),依據(jù)目標(biāo)可靠指標(biāo)預(yù)測(cè)新沂河大橋的剩余使用壽命。
表1 結(jié)構(gòu)破壞類(lèi)型對(duì)應(yīng)的折減目標(biāo)可靠指標(biāo)0.85β0和目標(biāo)失效概率
2.3.1JC法計(jì)算結(jié)構(gòu)可靠指標(biāo)
設(shè)結(jié)構(gòu)功能函數(shù)Z為隨機(jī)自變量Xi(i=1,2…n)的函數(shù),Z=g(x1,x2,…,xn)=R(x1,x2,…xn)-S(x1,x2,…,xn)。結(jié)構(gòu)可靠度即為求函數(shù)Z可靠度。為避免中心點(diǎn)法中的誤差,本文提出了驗(yàn)算點(diǎn)法計(jì)算函數(shù)可靠度。驗(yàn)算點(diǎn)法是將線(xiàn)性化點(diǎn)選在功能函數(shù)Z=0且為最大失效概率對(duì)應(yīng)點(diǎn)(x1*,x2*,…,xn*)上,在P*點(diǎn)上用Taylor級(jí)數(shù)線(xiàn)性化展開(kāi)函數(shù)求解結(jié)構(gòu)的可靠指標(biāo)β值,該方法被國(guó)際安全度委員會(huì)(JCSS)所推薦,簡(jiǎn)稱(chēng)為JC法。根據(jù)上述原理將功能函數(shù)用Taylor級(jí)數(shù)在P*點(diǎn)上展開(kāi),取一階項(xiàng)得功能函數(shù)Z均值和標(biāo)準(zhǔn)差。
則函數(shù)Z可靠指標(biāo)β計(jì)算方法如式(9)。
(9)
根據(jù)新沂河大橋的監(jiān)測(cè)數(shù)據(jù)統(tǒng)計(jì)規(guī)律,交通量年增長(zhǎng)按1%預(yù)測(cè),可得出評(píng)估期內(nèi)每年的總荷(活)載效應(yīng)分布,評(píng)估期內(nèi)的結(jié)構(gòu)抗力概率分布由式(9)計(jì)算,得出新沂河大橋函數(shù)Z時(shí)變可靠指標(biāo)β。JC法計(jì)算結(jié)果見(jiàn)表2和圖2,由此可知2021年新沂河大橋的時(shí)變可靠指標(biāo)降至4,而2022年時(shí)可靠指標(biāo)降至3.9,小于目標(biāo)可靠指標(biāo)3.995。因此,從2012年JC法預(yù)測(cè)新沂河大橋剩余使用壽命為10年。
圖2 JC法預(yù)測(cè)結(jié)構(gòu)壽命可靠度曲線(xiàn)(2012年對(duì)應(yīng)x=11)
表2 JC法預(yù)測(cè)結(jié)構(gòu)安全可靠指標(biāo)
2.3.2Monte-Carlo法計(jì)算結(jié)構(gòu)失效概率
Monte-carlo法又稱(chēng)統(tǒng)計(jì)試驗(yàn)法,它是利用計(jì)算機(jī)進(jìn)行大量抽樣模擬隨機(jī)變量的分布,反復(fù)進(jìn)行試驗(yàn)觀測(cè)取得服從相應(yīng)分布的隨機(jī)變量的子樣,用子樣的統(tǒng)計(jì)特征近似代表隨機(jī)變量母體的統(tǒng)計(jì)特征,進(jìn)而估計(jì)結(jié)構(gòu)體系的失效概率pf。新沂河大橋目標(biāo)失效概率為1.1×10-6,Monte-Carlo法計(jì)算結(jié)果見(jiàn)表3。由表3可知,2019年新沂河大橋的時(shí)變失效概率值達(dá)到6.3×10-6,首次超越目標(biāo)失效概率。因此,通過(guò)Monte-Carlo法計(jì)算新沂河大橋剩余使用壽命為8年。
表3 Monte-Carlo法預(yù)測(cè)結(jié)構(gòu)失效概率pf數(shù)據(jù)
2.3.3預(yù)測(cè)剩余壽命數(shù)據(jù)驗(yàn)證
對(duì)于新沂河大橋來(lái)說(shuō),JC法計(jì)算結(jié)果可作為Monte-Carlo法預(yù)測(cè)的剩余使用壽命結(jié)果的參考。JC法驗(yàn)證Monte-Carlo法預(yù)測(cè)剩余使用壽命結(jié)果是否可靠。JC法預(yù)測(cè)得出的新沂河大橋剩余使用壽命為10年,Monte-Carlo法預(yù)測(cè)得出的剩余使用壽命為8年,兩者結(jié)果較為接近。剩余壽命預(yù)測(cè)數(shù)據(jù)表明新沂河大橋技術(shù)狀況較差。
2013年6月利用在線(xiàn)監(jiān)測(cè)系統(tǒng)進(jìn)行靜載試驗(yàn)數(shù)據(jù)分析。新沂河大橋右幅第九聯(lián)中跨箱梁當(dāng)加載至4級(jí),相當(dāng)于汽車(chē)超-20級(jí)荷載等級(jí)時(shí),測(cè)試截面箱梁位移和應(yīng)變均存在實(shí)測(cè)值大于理論值的現(xiàn)象,并出現(xiàn)箱梁裂縫大量增加和裂縫寬度大幅增大現(xiàn)象。試驗(yàn)表明,新沂河大橋相應(yīng)結(jié)構(gòu)強(qiáng)度、剛度和抗裂性能不能滿(mǎn)足汽車(chē)超-20級(jí)荷載等級(jí)的使用要求。
1) 本項(xiàng)目立足于安全性和耐久性角度計(jì)算LCPA時(shí)變可靠度,采用JC法預(yù)測(cè)新沂河大橋剩余使用壽命為10年,采用Monte-Carlo法預(yù)測(cè)剩余使用壽命為8年,2種方法壽命預(yù)測(cè)結(jié)果較為接近。
2) 本項(xiàng)目建立了混凝土梁橋LCPA剩余使用壽命的預(yù)測(cè)流程和公式,結(jié)合項(xiàng)目實(shí)際應(yīng)用了荷載效應(yīng)模型及結(jié)構(gòu)抗力模型進(jìn)行可靠度評(píng)估,可為同類(lèi)橋梁的剩余使用壽命預(yù)測(cè)提供參考。
3) 本次LCPA剩余使用壽命預(yù)測(cè)結(jié)果與實(shí)橋荷載試驗(yàn)結(jié)果一致,2種不同評(píng)估結(jié)果均表明新沂河大橋劣化程度加劇,不滿(mǎn)足現(xiàn)行荷載的需求,需要及時(shí)進(jìn)行維修加固。
目前橋梁LCPA剩余使用壽命預(yù)測(cè)養(yǎng)護(hù)需求強(qiáng)烈,但研究還在發(fā)展中,尚有一些問(wèn)題待進(jìn)一步研究深化。
1) 本次LCPA橋梁剩余使用壽命預(yù)測(cè)方法僅考慮了結(jié)構(gòu)抗彎性能,可能的情況下還應(yīng)考慮抗剪、抗扭等其他承載力性能指標(biāo)。
2) 橋梁耐久性能劣化模型和車(chē)輛荷載(活載)效應(yīng)模型可分別從長(zhǎng)期監(jiān)(檢)測(cè)和動(dòng)態(tài)稱(chēng)重系統(tǒng)來(lái)進(jìn)行優(yōu)化,進(jìn)一步提高綜合耐久性和安全性橋梁LCPA剩余使用壽命預(yù)測(cè)方法的精度。
3) 本次LCPA橋梁剩余使用壽命預(yù)測(cè)仍然基于規(guī)范靜態(tài)抗力和靜態(tài)作用效應(yīng)力演變規(guī)律。后期還可以進(jìn)一步研究橋梁LCPA剩余使用壽命的動(dòng)態(tài)預(yù)測(cè)方法,例如基于神經(jīng)網(wǎng)絡(luò)算法、灰色系統(tǒng)模型、人工智能系統(tǒng)等各種預(yù)測(cè)方法。