張?jiān)戳?，王龍?/p>

(1.廣西北部灣投資集團(tuán)有限公司，廣西 南寧 530029；2.廣西交科集團(tuán)有限公司，廣西 南寧 530007)

0 引言

大跨斜拉橋由于外觀優(yōu)美、經(jīng)濟(jì)性能優(yōu)越以及跨越能力大等優(yōu)點(diǎn)，而逐漸成為大跨橋梁的首選橋型[1-3]。其中斜拉索作為斜拉橋的重要受力構(gòu)件，常發(fā)生如拉索腐蝕、錨頭銹蝕、索力損失等問(wèn)題，嚴(yán)重影響橋梁的安全性和可靠性[4-5]。因此，為保證斜拉橋的運(yùn)營(yíng)安全，對(duì)斜拉索進(jìn)行可靠度評(píng)估具有重要的意義。

斜拉索可靠度研究方面，劉綱等[6]通過(guò)分析斜拉索的疲勞損傷，基于概率密度演化理論研究了斜拉索疲勞可靠度；譚冬梅等[7]考慮風(fēng)荷載、車(chē)輛荷載以及覆冰荷載綜合作用，對(duì)斜拉索進(jìn)行了疲勞可靠度分析；劉發(fā)等[8]提出了腐蝕斜拉索承載力的退化模型，并研究了腐蝕程度和服役時(shí)間對(duì)斜拉索可靠度的影響；殷志祥等[9]基于Miner線(xiàn)性累積損傷理論評(píng)估了車(chē)輛荷載作用下斜拉索的疲勞可靠度水平；杜鵬剛等[10]通過(guò)響應(yīng)面法研究了斜拉橋靜力可靠度，分析表明，大跨斜拉橋主要失效路徑是由于外側(cè)斜拉索失效進(jìn)而導(dǎo)致主梁懸臂跨中發(fā)生彎曲破壞。上述研究可以看出，斜拉索可靠度的研究已取得較大進(jìn)展，眾多可靠度分析方法雖能較為準(zhǔn)確的進(jìn)行斜拉索可靠度分析，但大部分計(jì)算方法計(jì)算復(fù)雜，不便于工程應(yīng)用。欲將斜拉索可靠度分析應(yīng)用于工程實(shí)際，需首先選取高效可行的可靠度分析方法。目前，Monte Carlo法是可靠度分析領(lǐng)域公認(rèn)最為精確的分析方法，然而其計(jì)算成本昂貴，難以應(yīng)用于實(shí)際工程可靠度評(píng)估[11]。近年來(lái)，選點(diǎn)策略方法得到不斷發(fā)展，如重要性抽樣法[12]、子集模擬法[13]、線(xiàn)性抽樣法[14]以及Sobol序列抽樣法[15]。其中，Sobol序列抽樣法由于抽樣效率高，而廣泛應(yīng)用于結(jié)構(gòu)可靠度研究領(lǐng)域[16-18]。

鑒于此，本文基于Sobol序列和Monte Carlo法建立了斜拉索可靠度分析框架，以典型實(shí)際工程斜拉橋?yàn)檠芯繉?duì)象，通過(guò)Sobol序列選取隨機(jī)樣本點(diǎn)，并結(jié)合有限元分析方法和Monte Carlo法計(jì)算了該橋的斜拉索可靠指標(biāo)，以期為斜拉索可靠度分析提供參考。

1 基于Sobol序列的連續(xù)鋼構(gòu)橋可靠度分析方法

1.1 Sobol序列基本原理

Sobol序列是一種低差異序列，其分布均勻并且收斂較快，因此在各個(gè)領(lǐng)域得到了廣泛的應(yīng)用。假設(shè)ki為小于2i的正奇數(shù)，表達(dá)式為[19]：

(1)

式中，ki根據(jù)功能多項(xiàng)式產(chǎn)生，如式(2)所示。

f(x)=xr+c1xr-1+…+cr-1x+cr

(2)

當(dāng)i>r時(shí)，ki和mi存在以下遞歸公式：

ki=c1ki-1⊕c2vi-2⊕…⊕crvi-r⊕vi-r/2r

(3)

mi=2c1mi-1⊕22c2mi-2⊕…⊕2rcrmi-r⊕mi-r

(4)

式中，⊕表示二進(jìn)制中異或運(yùn)算。

存在任意整數(shù)O，唯一表示與數(shù)基b=2的表達(dá)式：

(5)

式中，p表示大于等于lbO的最小整數(shù)；aj取值為0或1。

對(duì)此，Sobol序列第z個(gè)元素可計(jì)算如下：

Θz=a1v1⊕a2v2⊕…⊕apvp

(6)

限于文章篇幅，Sobol序列采樣詳細(xì)構(gòu)造過(guò)程可參考文獻(xiàn)。以某9維隨機(jī)變量為例，采用Sobol序列抽取300個(gè)樣本點(diǎn)如下頁(yè)圖1所示。由圖1可以看出，Sobol序列采樣結(jié)果較為均勻，可高效地選取隨機(jī)樣本點(diǎn)，進(jìn)而應(yīng)用于結(jié)構(gòu)可靠度評(píng)估。

圖1 Sobol序列采樣結(jié)果示意圖

1.2 失效概率及可靠指標(biāo)計(jì)算

通過(guò)Sobol技術(shù)獲取樣本點(diǎn)后，可進(jìn)一步結(jié)合Monte Carlo法計(jì)算結(jié)構(gòu)的失效概率，計(jì)算公式如下：

(7)

已知結(jié)構(gòu)失效概率后，其可靠指標(biāo)可相應(yīng)計(jì)算得出：

β=-Φ-1(pf)

(8)

式中，Ф為標(biāo)準(zhǔn)正態(tài)分布函數(shù)。

2 基于Sobol序列的斜拉索可靠度分析流程及步驟

結(jié)合Sobol序列采樣以及Monte Carlo法，可對(duì)斜拉索開(kāi)展高效的可靠度分析工作，其具體流程如圖2所示，詳細(xì)步驟描述如下：

圖2 基于Sobol序列的斜拉索可靠度分析步驟示意圖

步驟1：基于通用有限元軟件ANSYS，建立斜拉橋非線(xiàn)性有限元分析模型；

步驟2：參考現(xiàn)有研究，確定斜拉索常見(jiàn)隨機(jī)變量并選取對(duì)應(yīng)的分布類(lèi)型；

步驟3：確定隨機(jī)模擬次數(shù)N，基于Sobol序列在步驟2生成N個(gè)確定的隨機(jī)變量樣本點(diǎn)；

步驟4：編制Matlab程序，實(shí)現(xiàn)樣本點(diǎn)批量寫(xiě)入有限元計(jì)算文件(Mac文件)，最終形成N個(gè)有限元分析樣本；

步驟5：基于N個(gè)有限元分析樣本，通過(guò)Matlab調(diào)用ANSYS并將計(jì)算結(jié)果存儲(chǔ)在指定文件夾；

步驟6：根據(jù)有限元分析結(jié)果，計(jì)算對(duì)應(yīng)斜拉索的功能函數(shù)值；

步驟7：基于Monte Carlo法，統(tǒng)計(jì)總樣本中失效次數(shù)，計(jì)算斜拉索失效概率和可靠指標(biāo)；

步驟8：根據(jù)斜拉索失效概率和可靠指標(biāo)，評(píng)價(jià)其安全性能。

3 基于Sobol序列抽樣的斜拉索可靠度分析

3.1 斜拉索有限元模型

結(jié)合上述分析方法，對(duì)某典型工程開(kāi)展斜拉索可靠度分析。該橋橋型為低塔斜拉橋，跨度布置為(110+220+110)m，斜拉索采用扇形布置，每個(gè)索面由10對(duì)斜拉索組成，其計(jì)算簡(jiǎn)圖以及斜拉索編號(hào)如圖3所示。

圖3 斜拉橋斜拉索編號(hào)示意圖

采用通用有限元軟件ANSYS建立斜拉橋有限元分析模型，斜拉索采用LINK10單元模擬，橋塔、主梁和橫向系梁采用BEAM4單元模擬。此外，橋塔底部以及兩側(cè)橫梁支座進(jìn)行節(jié)點(diǎn)約束，在橋塔橫梁和支座處設(shè)置彈簧約束。考慮到斜拉橋具有較高的幾何非線(xiàn)性，采用形狀迭代法確定斜拉索初始索力，最終斜拉橋有限元模型如圖4所示。

圖4 斜拉橋有限元力學(xué)模型圖

3.2 隨機(jī)變量及功能函數(shù)

隨機(jī)變量方面，本文主要考慮斜拉索彈性模量E1、容重γ1和主梁二期恒載q1三個(gè)隨機(jī)參數(shù)，其具體分布類(lèi)型、均值和標(biāo)準(zhǔn)差如表1所示。

表1 斜拉橋隨機(jī)變量及其分布類(lèi)型數(shù)值表

對(duì)于斜拉索而言，其在運(yùn)營(yíng)過(guò)程中常發(fā)生強(qiáng)度破壞，進(jìn)而影響斜拉橋整體承載性能。因此，本文可靠度分析主要考慮斜拉索的強(qiáng)度破壞，其功能函數(shù)可分別表示為：

(9)

3.3 可靠度分析

基于Sobol序列抽樣方法，即可對(duì)本文斜拉索進(jìn)行可靠度評(píng)估工作。首先，考慮到斜拉索失效破壞屬于小概率事件，因此根據(jù)表1中隨機(jī)變量參數(shù)及其分布類(lèi)型，采用Sobol序列抽樣方法抽取106個(gè)樣本點(diǎn)；隨后，調(diào)用斜拉橋有限元模型計(jì)算各個(gè)樣本點(diǎn)對(duì)應(yīng)的響應(yīng)信息；最后，結(jié)合Monte Carlo法統(tǒng)計(jì)拉索失效次數(shù)，計(jì)算斜拉索失效概率pf和可靠指標(biāo)β。此外，考慮到本文斜拉橋算例中斜拉索沿中跨跨中左右對(duì)稱(chēng)，因此本文僅計(jì)算了編號(hào)為A1～A10的斜拉索可靠度，其結(jié)果如表2所示。

表2 斜拉索可靠度計(jì)算結(jié)果表

由表2可知，斜拉索失效概率主要集中在10-4～10-5，對(duì)應(yīng)可靠指標(biāo)分布在3.334～4.226以?xún)?nèi)，屬于小概率失效事件。根據(jù)《公路鋼筋混凝土及預(yù)應(yīng)力混凝土橋涵設(shè)計(jì)規(guī)范》(JTG D62-2012)[20]要求：橋梁結(jié)構(gòu)一級(jí)結(jié)構(gòu)承載能力極限狀態(tài)下結(jié)構(gòu)發(fā)生延性破壞的可靠指標(biāo)β不應(yīng)低于4.7。根據(jù)表2計(jì)算結(jié)果，本文算例的斜拉索可靠指標(biāo)最大為4.226，低于規(guī)范中可靠度要求，建議對(duì)該橋采取相應(yīng)的加固措施。

圖5給出了斜拉索位置和可靠指標(biāo)的對(duì)應(yīng)關(guān)系，由圖5可知，隨著斜拉索編號(hào)的增加，即拉索長(zhǎng)度的減小，斜拉索的可靠指標(biāo)逐漸增大，即安全水平逐漸提高。而最外圍拉索A1由于所承擔(dān)的索力最大，因此可靠指標(biāo)最小，拉索最易失效，應(yīng)予以重視。

圖5 斜拉索可靠指標(biāo)曲線(xiàn)圖

4 結(jié)語(yǔ)

斜拉索的可靠度水平對(duì)斜拉橋的運(yùn)營(yíng)安全具有重要的影響，本文以某大跨斜拉橋?yàn)槔?，基于Sobol序列抽樣方法對(duì)斜拉索進(jìn)行高效的隨機(jī)抽樣，并結(jié)合Monte Carlo法計(jì)算了斜拉索可靠度水平。分析結(jié)果表明：該算例中斜拉索可靠指標(biāo)相較于規(guī)范要求偏低，建議對(duì)斜拉索采取一定的加固措施；隨著拉索長(zhǎng)度的增大，其承受的索力越大，可靠指標(biāo)隨著降低，應(yīng)重視外圍拉索的損傷狀態(tài)；此外，基于Sobol序列抽樣方法應(yīng)用簡(jiǎn)單，能切實(shí)可行的應(yīng)用于斜拉索可靠度分析，可為斜拉索概率評(píng)估提供參考。