水利工程論文-防洪堤的可靠性設(shè)計方法探討.doc

水利工程論文-防洪堤的可靠性設(shè)計方法探討摘要：基于可靠性理論，提出了防洪堤的概率設(shè)計方法。按照堤防設(shè)計規(guī)范，建立典型斜墻式防洪堤漫頂、滲透穩(wěn)定、邊坡滑動穩(wěn)定的極限狀態(tài)方程，采用MonteCarlo法計算相應(yīng)的可靠性指標(biāo)。進而通過變動堤防幾何參數(shù)(堤頂高程、頂寬和坡比)和土工參數(shù)，探討了堤高與漫頂可靠性、堤防幾何參數(shù)的變化對堤基滲透穩(wěn)定的可靠性、坡率對邊坡失穩(wěn)可靠性的影響。結(jié)果表明：隨著堤高的增加，堤防漫頂?shù)目煽啃灾笜?biāo)逐漸增大；堤防坡率、堤頂寬度、堤基粘土層厚度、前灘地寬度、堤后壓蓋厚度和寬度等變量的增大，將使堤基滲透穩(wěn)定的可靠性指標(biāo)增大；堤基粘土層和砂土層滲透系數(shù)的增加，將使堤基滲透穩(wěn)定可靠性指標(biāo)減?。坏谭郎?、下游坡率的增加，堤坡穩(wěn)定可靠性指標(biāo)增加。關(guān)鍵詞：堤防可靠性變異性分析可靠性設(shè)計由于受到堤身斷面形狀、堤身填土性質(zhì)、堤基地質(zhì)、水文、地形和施工條件等諸多因素的影響，堤防在汛期往往出現(xiàn)漫頂、堤身(基)滲透破壞、堤身滑坡等險情。依據(jù)堤防的設(shè)計與使用要求，在使用年限內(nèi)應(yīng)充分考慮影響這些險情的因素，以滿足滲透穩(wěn)定和堤坡整體穩(wěn)定的要求。按照現(xiàn)有的堤防設(shè)計規(guī)范，設(shè)計方法是基于極限平衡分析的安全系數(shù)(或分項安全系數(shù))法，它無法給出工程可靠性的評價指標(biāo)。概率設(shè)計方法將結(jié)構(gòu)的多種參數(shù)作為隨機變量，可根據(jù)不同堤防結(jié)構(gòu)重要性程度采取相應(yīng)設(shè)計，無疑有其先進性和科學(xué)性。因此，以概率理論為基礎(chǔ)的可靠性設(shè)計方法在國內(nèi)外均得到較快的發(fā)展，這為巖土工程的設(shè)計開辟了新的領(lǐng)域。荷蘭的堤壩和護岸工程的一些設(shè)計導(dǎo)則中提出了概率設(shè)計和風(fēng)險分析方法1。國內(nèi)一些學(xué)者2-4對邊坡等工程的可靠性設(shè)計方法進行了系統(tǒng)研究。結(jié)合防洪堤的結(jié)構(gòu)特征與材料組成特點，綜合考慮多種失穩(wěn)因素，提出簡便實用的可靠性設(shè)計方法是值得深入探討的課題。作者采用MonteCarlo法，以二元結(jié)構(gòu)堤基的斜墻式典型堤防為算例，探討幾何設(shè)計參數(shù)和土工參數(shù)對堤防漫頂、滲透穩(wěn)定、岸坡穩(wěn)定可靠性指標(biāo)的影響，并提出堤防結(jié)構(gòu)可靠性設(shè)計方法。這將有助于完善現(xiàn)有堤防設(shè)計規(guī)范，最終與穿堤建筑物、防滲墻和防浪墻等可靠性研究共同實現(xiàn)堤防工程系統(tǒng)設(shè)計方法向概率型轉(zhuǎn)變。1可靠性的計算原理工程結(jié)構(gòu)可靠性常用的計算方法，主要有驗算點法和MonteCarlo法。驗算點法是在設(shè)計驗算點處泰勒級數(shù)二項式展開，把非線性函數(shù)簡化為線性函數(shù)，故計算精度不高，對非線性函數(shù)存在收斂性問題。MonteCarlo法是通過大量抽樣模擬各隨機變量以達到計算可靠性之目的，對于極限狀態(tài)函數(shù)形式為非線性的堤防可靠性分析較為適用，且具有方法簡單、易于編程的特點，其基本求解原理5如下述。1.1隨機抽樣方法先在(0，1)區(qū)間上產(chǎn)生均勻分布隨機數(shù)?？刹捎玫姆椒ㄓ谐送喾ā⒒旌贤喾ǖ?，其中混合同余法的遞推公式為xi=(xi-1+C)(modM)ri=xi/M(i=1,2,n)(1)式中：、x0、C和M為選定的常數(shù)。式(1)表示xi-1+C除以M的余數(shù)為xi，將其再除以M即得(0，1)上均勻分布的隨機數(shù)ri。將隨機數(shù)序列ri轉(zhuǎn)換為(a,b)區(qū)間上的均勻分布隨機數(shù)序列RiRi=a+(b-a)ri(2)采用反函數(shù)法將均勻分布隨機數(shù)ri轉(zhuǎn)換為符合某一指定概率分布的隨機數(shù)，其前提是經(jīng)驗分布的反函數(shù)存在，否則采用隨機變量函數(shù)法。設(shè)X為具有分布函數(shù)FX(x)且反函數(shù)F-1X(x)存在的連續(xù)隨機變量，r是均勻分布隨機變量R(分布函數(shù)為FR(r)的值。若給定累計概率FX(x)=r，則有x=F-1X(r)(3)若已知ri，則可得符合FX(x)的隨機數(shù)序列：xi=F-1X(ri)(i=1,2,n)(4)1.2MonteCarlo模擬法的基本步驟(1)統(tǒng)計確定與可靠性分析有關(guān)的各基本變量(如材料參數(shù)、荷載等)的概率分布模型及其分布參數(shù)；(2)對所有基本變量按統(tǒng)計特征進行第一次隨機采樣，將其代入極限狀態(tài)方程式，獲得第一次模擬結(jié)果；(3)重復(fù)n次獨立隨機采樣，進而估算失效概率。1.3MonteCarlo模擬法的結(jié)果和精度在工程結(jié)構(gòu)可靠性分析中，設(shè)極限狀態(tài)方程Z=g(x1,x2,xn)，其失效概率為Pf=P(g(x1,x2,xn)0)(5)借助隨機抽樣對基本變量賦值時，計算結(jié)果只有g(shù)()0和g()0兩種可能，故可定義指標(biāo)函數(shù)：(6)根據(jù)伯努利大數(shù)定理及正態(tài)分布隨機變量的特性得到的失效概率Pf為(7)式中：M為N次模擬計算中g(shù)()0的總次數(shù)。MonteCarlo模擬法不僅可按式(7)計算失效概率，還可根據(jù)模擬計算結(jié)果擬合出功能函數(shù)的分布形式，據(jù)以估計其一、二階矩x和2x。進而按下式近似計算可靠度指標(biāo)(8)MonteCarlo模擬法的誤差一般可用x表示。功能函數(shù)值Z越離散，誤差將越大。當(dāng)模擬次數(shù)充分大時，已經(jīng)證明由模擬結(jié)果樣本求得估計值的標(biāo)準(zhǔn)差與模擬次數(shù)的平方根成反比。因此，加大模擬次數(shù)可望提高模擬精度，特別是當(dāng)計算得到的失效概率較小時。一般要求Montecarlo法的樣本數(shù)必須大于引起一次Z0所需的平均樣本數(shù)的100倍，即N100/Pf次時，可滿足精度要求。由此編制了相應(yīng)的計算程序。2某斜墻式防洪堤的概率設(shè)計以某典型斜墻式防洪堤的概率設(shè)計為例，如圖1所示。假定其堤基為由表層弱透水粘性土、下臥強透水砂層組成的二元結(jié)構(gòu)。表1中給出計算中采用的土工強度參數(shù)以及洪水位等參數(shù)的取值，表中均值一列給出典型堤防的幾何參數(shù)取值。2.1漫頂防洪水位是確定堤頂高程的決定性因素。設(shè)防洪堤發(fā)生漫頂?shù)臉O限狀態(tài)方程z1為z1=h0-hw-hs-e(9)式中：hs為波浪爬高；e為風(fēng)壅高度。由于影響堤頂高程的因素很多，如堤基沉降?？紤]到施工精度及竣工后固結(jié)沉降等因素，為保持設(shè)計高程，在設(shè)計時需預(yù)留沉降量。假定堤頂高程h0服從正態(tài)分布，在一定堤段內(nèi)預(yù)留沉降量超過0.1m的概率小于2%，由正態(tài)分布的表達式可得hd-h0/h0=1.96(10)式中：hd為堤防的竣工后高程。表1堤防參數(shù)的取值符號名稱和單位單位類型均值標(biāo)準(zhǔn)差h0hbcL1LkmWtsb堤頂高程堤基粘土層厚度前灘地寬度有效滲徑寬度坡率堤頂寬度透水壓蓋的厚度mmmmmm設(shè)計變量設(shè)計變量設(shè)計變量設(shè)計變量設(shè)計變量設(shè)計變量設(shè)計變量11.03.510.068.752.55.00hbs堤基砂土層厚m常量2.5nkwsb度粘土的浮容重水的容重壓蓋土體的容重kN/m2kN/m2kN/m2常量常量常量9.010.018.0dkshwkcks渠底粘土層有效厚度洪水位粘土的滲透系數(shù)砂土的滲透系數(shù)mmm/sm/s正態(tài)指數(shù)正態(tài)正態(tài)3.58.3410-810-50.70.9210-80.7510-5由式(10)可得h0=0.1/1.96=0.051(11)根據(jù)通常的經(jīng)驗，最高洪水位hw假定服從指數(shù)分布，取hw=8.34m，hw=0.9。波浪爬高hs假定服從正態(tài)分布，并取hs=0.696。根據(jù)文獻7代入蒲田波浪要素公式計算波高的均值hgh/V2=0.13th0.7(gha/V2)0.7th0.0018(gF/V2)0.450.1