鄒俏俏

（遼寧省水利水電勘測設(shè)計研究院，遼寧 沈陽 110006）

在土壩壩坡穩(wěn)定分析方面，目前還是采用傳統(tǒng)的安全系數(shù)法，而把可靠度分析應(yīng)用到工程設(shè)計和現(xiàn)有工程的安全評價上是今后的發(fā)展趨勢。在土壩壩坡穩(wěn)定性分析方面，安全系數(shù)法是將影響土壩壩坡狀態(tài)的諸因素看作確定性變量，通過定值安全系數(shù)作為衡量土壩壩坡穩(wěn)定狀態(tài)的指標(biāo)。而大量的試驗和工程實踐證明，影響土壩壩坡穩(wěn)定狀態(tài)的因素中有許多具有很大的隨機性，如土壩壩坡的強度參數(shù)、外界荷載、自然條件、蓄水水位及土體內(nèi)各種不連續(xù)結(jié)構(gòu)面等，由于安全系數(shù)方法回避了計算參數(shù)的隨機性，從而給土壩壩坡穩(wěn)定性評價帶來很大的誤差，不少土壩壩坡工程從安全系數(shù)角度分析是安全的，但實際運營中卻發(fā)生滑坡破壞，具有同一安全系數(shù)的土壩壩坡實際穩(wěn)定狀態(tài)常常不同。而基于隨機性分析的土壩壩坡穩(wěn)定性分析方法是通過建立概率極限狀態(tài)方程，充分考慮土體參數(shù)的均值及其變異性，從實際意義上來說更為合理[1～2]。

對于大型水利工程，在定值安全系數(shù)分析的基礎(chǔ)上進行可靠度復(fù)核，可以提高工程穩(wěn)定性評價的精度。而一次二階矩法（FOSM）是常用的方法，包括中心點法和驗算點法，現(xiàn)結(jié)合定值安全系數(shù)與可靠度結(jié)合起來分析評價土壩壩坡穩(wěn)定性的方法。

1 基本理論

1.1 可靠度分析的基本理論

結(jié)構(gòu)破壞的可能性稱之為失效概率，以Pf表示，成功的概率為可靠度，以Pa表示，結(jié)構(gòu)在規(guī)定的時間內(nèi)，在規(guī)定的條件下具有預(yù)定功能的概率稱之為可靠度分析，且有Pa=1-Pf，工程設(shè)計中，可靠度分析與定值安全系數(shù)結(jié)合可以提高穩(wěn)定性分析的精度[3]。以可靠度的角度講安全系數(shù)應(yīng)滿足：Pf·DL＜△C，式中：DL為破壞損失；△C為更可靠設(shè)計所需造價?？梢姡?dāng)前兩者乘積較小時，使用較小的安全系數(shù)可大幅降低造價；反之，當(dāng)前兩者乘積較大時，使用較大的安全系數(shù)可以避免失事后的損失和善后費用[4]。

有資料顯示，在既定的土壩壩坡幾何形狀下，任何確定的破壞面總存在一個與之相應(yīng)的安全系數(shù)Fs和一個與之相應(yīng)的Pf。而且，在所有破壞面中，存在一個最小Fs（min）和最大Pf（max）的最危險破壞面。由于基本狀態(tài)參數(shù)是隨機變量，具有離散性，所以相應(yīng)最危險位置并非確定。但在土壩壩坡穩(wěn)定性問題中把破壞面位置固定而產(chǎn)生的影響并不大，采用基本變量均值確定相應(yīng)最小安全系數(shù)的破壞面作為可靠度分析的破壞面，完全可以滿足工程精度要求。

1.2 一次二階矩法探討

一次二階矩法（FOSM）進行土壩壩坡可靠度分析，是在影響土壩壩坡穩(wěn)定的隨機變量的分布尚不清楚時，采用一階矩（均值）和二階矩（標(biāo)準(zhǔn)差）的數(shù)學(xué)模型對土壩壩坡可靠度進行評價。建立功能函數(shù)，采用泰勒級數(shù)在平均值處將極限狀態(tài)功能函數(shù)展開，使之線性化后求解。

由于土體是一種高度非線性材料，在不同的應(yīng)力下具有不同的變形特性，因而，如果仍按照結(jié)構(gòu)工程中隨機變量互相獨立，服從正態(tài)分布的可靠度計算，將功能函數(shù)在基本變量的均值點上線性化，將出現(xiàn)迭代不收斂現(xiàn)象，影響可靠度評價的成果。

為了提高可靠度評價精度，計算結(jié)構(gòu)可靠指標(biāo)時，針對土體的一個有效方法就是將隨機變量在驗算點處“當(dāng)量正態(tài)”化為正態(tài)分布。這樣就可以采用與一次二階矩的方法求解結(jié)構(gòu)的可靠性指標(biāo)，由于此法被國際結(jié)構(gòu)安全度聯(lián)合委員會（JCSS）所采用，故也叫“JC”法。

1.3 一次二階矩法的基本原理

假定隨機變量已轉(zhuǎn)化為統(tǒng)計獨立情況，并進行了等效正態(tài)化。將功能函數(shù) g（Y）=（y1，y2，…，yn）在某點 Y*=（y*1，y*2，…，y*n）展開，使之線性化，然后計算功能函數(shù)在Y*=（y*1，y*2，…，y*n）的均值、標(biāo)準(zhǔn)差和可靠度。

初始值 y*i=（i=1，2，…，n）通常使功能函數(shù) g（Y）有一個非零的d值：

對于下一次迭代△y*i，認為y*i+△y*i將位于破壞面上，即有：

將上式在y*i展開，忽略△y*i二階以上小量：

將式（3）減式（1）得：

式中：λi為Cz的特征值，Cz為矩陣值；β為可靠度指標(biāo)，是一變量，β+△β是下一個變量值。

在上式中忽略△α*iY得：

1.4 一次二階矩法具體計算步驟[5]：

1）定義 g（X）=（x1，x2，…，xn），假定初始破壞點x*i=（i=1，2，…，n），將其標(biāo)準(zhǔn)正態(tài)化為z*i=（i=1，2，…，n），如果xi統(tǒng)計相關(guān)，計算得到獨立標(biāo)準(zhǔn)化變量y*i=（i=1，2，…，n），并確定β的初始迭代值。

2）將y*i代入式（1）得d值。

4）由式（12）得到△β，把β+△β 作為β 下一次的迭代值。

5）利用上面求得的β，由式（5）和式（10）重新計算z*i=（i=1，2，…，n）或y*i=（i=1，2，…，n）值，作為下步迭代的初始值。

6）重復(fù)以上步驟，直至滿足收斂條件△β＜εβ，d＞ε（dεβ，εd為允許誤差），這樣用一次二階矩法求解的可靠度收斂效果很好。

2 土壩壩坡穩(wěn)定可靠度分析計算

2.1 統(tǒng)計相關(guān)變量的獨立標(biāo)準(zhǔn)化

在可靠度分析中，對隨機變量Zi（i=1，2，…，n）的統(tǒng)計相關(guān)情況，通常把這些變量通過正交變換轉(zhuǎn)化為獨立的標(biāo)準(zhǔn)化變量Yi（i=1，2，…，n），它們的相關(guān)系數(shù)表示為：

通過正交矩陣T，將Zi變換為不相關(guān)變量：

Y=（y1，y2，…，yn）T，Y=TTT

T的列向量為Cz的特征向量，可靠度指標(biāo)β變?yōu)椋?/p>

2.2 非正態(tài)隨機變量的當(dāng)量正態(tài)化

非正態(tài)隨機變量按照驗算點處的累積概率分布函數(shù)CDF值和概率密度函數(shù)PDF值與原來分布相同的原則當(dāng)量正態(tài)化。很多研究表明，基本變量f服從正態(tài)分布，而基本變量c服從對數(shù)正態(tài)分布，我們需要將其當(dāng)量正態(tài)化[4]，具體方法如下：

式中：m′xi，σ′xi，Vxi，X*i，mxi*分別為等效正態(tài)分布的均值、標(biāo)準(zhǔn)差、變量xi的變異系數(shù)、驗算點及變量在驗算點處原分布的均值。

2.3 建立土壩壩坡極限穩(wěn)定狀態(tài)功能函數(shù)及其分析

在定值安全系數(shù)法計算的基礎(chǔ)上，選定某些參數(shù)為基本變量，并令按定值安全系數(shù)法得到的安全系數(shù)Fs=1，用瑞典圓弧法、Bishop法及Janbu法作為定值分析方法[6]，列出相應(yīng)于各種分析方法的極限狀態(tài)方程。

1）瑞典圓弧法的極限狀態(tài)功能函數(shù)為：

2）Bishop法的極限狀態(tài)功能函數(shù)為：

3）Janbu法的極限狀態(tài)功能函數(shù)為：

其中：mαi=cosαi+tanφ′isinαi：xi=△Eiti/bi-Eitanαi；△xi=xi+1-xi；Ei為土條兩側(cè)法向作用力；Xi為土條兩側(cè)切向作用力。

2.4 土壩壩坡極限穩(wěn)定狀態(tài)功能函數(shù)分析

采用定值安全系數(shù)法和一次二階距可靠度分析法對土壩壩坡進行最危險滑裂面搜索，確定最大失效概率的滑裂面，再反算對應(yīng)該面的安全系數(shù)，從而把定值安全系數(shù)法和可靠度分析法結(jié)合起來進行對比評價。

針對3個功能函數(shù)，對應(yīng)最大失效概率的滑裂面，不同計算方法得出的安全系數(shù)、可靠指標(biāo)、失效概率是不同的，同時c，φ的相關(guān)性也會影響計算結(jié)果。經(jīng)驗算比較，瑞典圓弧法得到的可靠指標(biāo)最小，失效概率最大，Bishop法得到的可靠指標(biāo)居中，而Janbu法得到的可靠指標(biāo)最大，失效概率最低。這可以從3種方法的原理上加以解釋，瑞典法完全沒有考慮土條間的相互作用，Bishop法考慮了土條間的水平作用力，而沒有考慮了土條間的垂直作用力。Janbu法同時考慮了土條間的水平和垂直作用力，并且其可靠指標(biāo)的大小與安全系數(shù)的計算結(jié)果也具有一致性，彼此結(jié)果符合得較好。

考慮土層抗剪強度參數(shù)c，φ的相關(guān)性時，安全系數(shù)與可靠指標(biāo)及失效概率的計算結(jié)果與不考慮其相關(guān)性時相差不大，變量之間的相關(guān)性對可靠指標(biāo)的影響很小，可以認為是相互獨立的。

2.5 清河水庫除險加固工程背水坡壩坡穩(wěn)定可靠度計算

清河水庫是遼寧省省屬大型水庫，2010年12月水利部大壩安全管理中心對清河水庫大壩進行安全鑒定復(fù)核，認定為Ⅲ類壩，需要進行除險加固。在除險加固設(shè)計過程中，經(jīng)過水文計算的復(fù)核，清河水庫大壩需要加高1.15m，經(jīng)過反復(fù)方案對比分析，決定采取帶帽加高的方式。為了節(jié)省投資，大壩背水坡采取改變原坡比的辦法進行加固處理，采用溢洪道尾渠的開挖料填筑。為此，對清河水庫大壩背水坡壩坡穩(wěn)定計算是非常必要的，同時對穩(wěn)定可靠度進行了計算。出險加固前背水坡比自上而下分別為：1∶2.25，1∶2.25，1∶2.5，出險加固后壩坡坡比分別為 1∶2.0，1∶2.25，1∶2.5；其填筑料為3層，即壩體填筑砂料，溢洪道尾渠開挖料，背坡護坡碎石料。

采用一次二階矩法進行可靠度計算，設(shè)計安全系數(shù)1.40。c，φ，γ作為基本變量，概率分布分別采用對數(shù)正態(tài)分布、非正態(tài)分布和正態(tài)分布，根據(jù)文中所述方法對其數(shù)據(jù)進行相應(yīng)轉(zhuǎn)化處理。得出表1的數(shù)據(jù)。

表1 安全系數(shù)及可靠指標(biāo)計算結(jié)果

由表1計算結(jié)果看出，不同計算方法得出的安全系數(shù)、可靠指標(biāo)、失效概率是不同的，同時c，φ的相關(guān)性也影響計算結(jié)果。瑞典圓弧法得到的可靠指標(biāo)最小，失效率最大；Bishop法居中，Janbu法得出的可靠指標(biāo)最大，失效率最小。瑞典法沒有考慮土條間的相互作用，而Bishop法只考慮了土條的水平力，Janbu法即考慮水平力也考慮垂直力的作用。

3 結(jié)語

在土壩壩坡穩(wěn)定性分析中，考慮分條的作用力可以使其穩(wěn)定系數(shù)、可靠指標(biāo)有較大提高；影響土壩壩坡穩(wěn)定性的主要因素是土體抗剪強度指標(biāo)，內(nèi)摩擦角的變異性對可靠指標(biāo)的影響大于內(nèi)聚力的變異性對可靠指標(biāo)的影響，可靠指標(biāo)對內(nèi)摩擦角的變化更為敏感；變量之間的相關(guān)性對可靠指標(biāo)的影響很小，可認為是相互獨立的；在隨機變量一定之條件下，可靠指標(biāo)β隨著安全系數(shù)Fs的增加而增加。

由于土壩壩坡工程中天然土體材料有較大的變異性，采用安全系數(shù)大小來評價土壩壩坡的穩(wěn)定性，掩蓋了土壩壩坡的破壞風(fēng)險，因而在土壩壩坡穩(wěn)定分析中引入可靠度分析是很有必要的。雖然土體主要結(jié)構(gòu)面的主要力學(xué)參數(shù)的取值方法和試驗條件差異較大，計算可靠度指標(biāo)所依據(jù)的功能函數(shù)也不固定，但由于可靠度分析是以定性分析為基礎(chǔ)，可以克服傳統(tǒng)定值安全系數(shù)分析中的缺陷，因而在定值安全系數(shù)計算基礎(chǔ)上進行可靠度復(fù)核，將其結(jié)果作為定值安全系數(shù)法的修正和補充，可以更加科學(xué)準(zhǔn)確的分析土壩壩坡的穩(wěn)定性。

［1］章建軍，陳傳堯，揭敏.邊坡抗滑穩(wěn)定可靠性分析［J］.華中理工大學(xué)學(xué)報，1997（04）：90-92.

［2］李獻民，王永和.黏性土坡穩(wěn)定的安全系數(shù)和可靠度［J］.湘潭礦業(yè)學(xué)院學(xué)報，2002（02）：82-86.

［3］祝玉學(xué).邊坡可靠性分析［M］.冶金工業(yè)出版社，1993：139-143.

［4］胡向陽.邊坡穩(wěn)定計算中的可靠度分析［J］.西北水電，2002（03）：12-16

［5］陳祖煜.土質(zhì)邊坡穩(wěn)定分析—原理方法程序［M］.中國水利水電出版社，2003：294-300.

［6］錢家歡，殷宗澤.土工原理與計算（第二版）［M］.北京：中國水利水電出版社，1996.