洪水作用下丁壩可靠度分析及剩余壽命預(yù)測(cè)

韓林峰，王平義，劉懷漢，杜 飛

(1．重慶交通大學(xué) 河海學(xué)院，重慶 400074;2．長(zhǎng)江航道規(guī)劃設(shè)計(jì)研究院，湖北 武漢 430011)

航道整治成敗的關(guān)鍵在于治理效果，而整治建筑物的穩(wěn)定性是確保治理效果的重要基礎(chǔ)．在所有結(jié)構(gòu)工程領(lǐng)域，運(yùn)用可靠度理論和方法進(jìn)行分析和設(shè)計(jì)已成為一種趨勢(shì)．與巖土、土木、水利等建筑物破壞形式不同，航道整治建筑物破壞因素極為復(fù)雜［1－3］:一方面，整治建筑物密實(shí)度較差、結(jié)構(gòu)強(qiáng)度低、基礎(chǔ)可動(dòng)性強(qiáng);另一方面，水流、河勢(shì)等外部因素變化大，加之其工作環(huán)境一般位于水下，條件復(fù)雜且難以進(jìn)行實(shí)際觀測(cè)．正是由于影響整治建筑物破壞的因素尚不明確，對(duì)其水毀機(jī)理的認(rèn)識(shí)不夠深入，因此尚未形成公認(rèn)的整治建筑物失穩(wěn)破壞計(jì)算模型和極限狀態(tài)方程，對(duì)于局部穩(wěn)定的驗(yàn)算盡管規(guī)范中有計(jì)算公式，但由于經(jīng)驗(yàn)性較強(qiáng)，往往計(jì)算值與實(shí)際情況有較大偏差，使得可靠度計(jì)算方法在整治建筑物中的應(yīng)用缺乏必要的理論基礎(chǔ)．目前整治建筑物的設(shè)計(jì)仍采用傳統(tǒng)的單一安全系數(shù)法，未考慮設(shè)計(jì)變量中客觀存在的變異性，因而不能全面準(zhǔn)確地反映整治建筑物的可靠性．基于現(xiàn)階段國內(nèi)外對(duì)整治建筑物可靠性及其使用壽命的研究寥寥無幾，因此如果能夠從整治建筑物的破壞機(jī)理出發(fā)，在充分考慮設(shè)計(jì)變量隨機(jī)性的基礎(chǔ)上分析整治建筑物的可靠度顯得尤為重要．

本文以整治建筑物中最常見的丁壩為研究對(duì)象，通過原型觀測(cè)和水槽試驗(yàn)，發(fā)現(xiàn)丁壩水毀主要是由于洪水漲退時(shí)壩體遭受洪水沖擊破壞，以及水流挾沙力增強(qiáng)導(dǎo)致壩頭附近沖刷坑增大使得壩體由于自身重力作用而導(dǎo)致的失穩(wěn)坍塌破壞．洪水作用下丁壩的破壞，從幾何結(jié)構(gòu)角度來看，原壩體結(jié)構(gòu)發(fā)生變形破壞并形成新的穩(wěn)定結(jié)構(gòu);從力學(xué)性能上來看，新結(jié)構(gòu)的結(jié)構(gòu)抗力、阻水能力等均有所降低．借鑒損傷力學(xué)中的宏觀分析方法，引進(jìn)損傷變量(本文指水毀體積)作為內(nèi)部變量［4］，通過宏觀現(xiàn)象分析壩體水毀過程，并由此建立丁壩水毀的本構(gòu)關(guān)系．

1 洪水條件下丁壩可靠度分析

1．1 洪水循環(huán)破壞試驗(yàn)設(shè)計(jì)

每次洪水的洪峰流量、洪水總量、洪水過程線形狀等各不相同，每年發(fā)生的洪水次數(shù)、歷時(shí)又各不相同．考慮到實(shí)際洪水過程與丁壩水毀之間的關(guān)系十分復(fù)雜、可變因素較多，因此本次水槽試驗(yàn)采用概化高斯曲線來模擬洪水過程．

根據(jù)長(zhǎng)江重慶航運(yùn)工程勘察設(shè)計(jì)院提供的“長(zhǎng)江中游荊江河段已建航道整治建筑物現(xiàn)狀調(diào)查及損毀原因分析”專題研究報(bào)告中的三峽水庫蓄水后壩下游宜昌、枝城、沙市、監(jiān)利、螺山、漢口等站歷年洪水資料，分析各次洪水的最大洪峰流量及其對(duì)應(yīng)的實(shí)測(cè)最大流速，并將其按照發(fā)生的頻率及流量大小分成不同的洪水段，再通過比尺換算到水槽試驗(yàn)中來，可以算出試驗(yàn)中洪峰流量在38．7～100 L/s范圍內(nèi)的分布，找出幾次發(fā)生頻率較大且有代表性的洪水，其洪峰流量分別為38．7，50，68和100 L/s，并將洪水過程概化成高斯曲線(見圖1)．水槽試驗(yàn)是在重慶交通大學(xué)國家內(nèi)河航道整治工程技術(shù)研究中心航道整治試驗(yàn)大廳長(zhǎng)30 m、寬2 m、高1 m的矩形玻璃水槽中進(jìn)行，最大可產(chǎn)生近0．2 m3/s的流量．本次試驗(yàn)的壩體結(jié)構(gòu)是長(zhǎng)江中下游常見的正挑梯形斷面拋石丁壩，迎水坡坡度為1∶1．5，背水坡坡度為1∶2，向河坡為1∶5，試驗(yàn)采用1∶60的正態(tài)模型．

圖1 水槽試驗(yàn)概化洪水過程曲線Fig．1 Generalized flood curves given by flume experiment

圖2 洪峰流量S與水毀循環(huán)次數(shù)N的關(guān)系Fig．2 Relationship between flood peak discharge S and water-damage cycle times N

天然河道中每一次的洪水歷時(shí)都不相同，但本次試驗(yàn)為了找出丁壩在不同洪峰流量下的水毀過程，將不同洪峰流量下每次洪水循環(huán)歷時(shí)視為相同，即ti－ti－1=C(C為常數(shù))．設(shè)S為最大洪峰流量，N為使丁壩水毀以至于無法滿足整治要求時(shí)的洪水循環(huán)次數(shù)(丁壩洪水循環(huán)使用壽命)．從水槽試驗(yàn)中得到，當(dāng)最大洪峰流量為38．7，50，68，100和120 L/s(依次為工況1～5)時(shí)，經(jīng)洪水循環(huán)作用丁壩水毀體積達(dá)到壩頭總體積30%時(shí)洪水循環(huán)次數(shù)分別為100，12，5，2和1．試驗(yàn)中當(dāng)洪峰流量為38．7 L/s時(shí)，無論循環(huán)多少次丁壩水毀體積都達(dá)不到30%，因此以100次來表示;這里的120 L/s是為了應(yīng)對(duì)將來有可能發(fā)生的超大洪水而進(jìn)行的模擬性試驗(yàn)．由此得到壩體水毀S－N關(guān)系曲線，如圖2所示．

這里需要說明的是丁壩作為最常見的一種河道整治建筑物，在建成投入使用一段時(shí)間后，壩體會(huì)出現(xiàn)不同程度的水毀，但以往經(jīng)驗(yàn)表明并不是只要丁壩發(fā)生損壞就一定失去其整治功能，有時(shí)壩體雖部分水毀但仍可以滿足整治要求．究竟丁壩水毀程度為多少時(shí)認(rèn)為其失效，航道部門一直沒有給出明確的規(guī)定．本文通過在清水沖刷試驗(yàn)中對(duì)壩頭斷面平均流速與壩體水毀體積的跟蹤測(cè)量，提出以丁壩水毀體積V1與壩頭總體積V0的比值V1/V0為表征丁壩水毀程度的指標(biāo)，并認(rèn)為當(dāng)丁壩水毀體積達(dá)到壩頭總體積的30%時(shí)丁壩失效，即當(dāng)V1=0．3V0時(shí)的水毀體積為臨界失效水毀體積．

選擇30%作為丁壩失效的臨界值，其原因在于清水沖刷試驗(yàn)中，對(duì)某一工況下丁壩壩軸線斷面上遠(yuǎn)離壩頭10 cm處的斷面平均流速進(jìn)行跟蹤，從水流調(diào)穩(wěn)時(shí)開始，到壩體水毀體積大約達(dá)到壩頭總體積的30%左右時(shí)結(jié)束，此時(shí)壩頭最大平均流速約減小50%(圖3)．而在整治工程中丁壩的主要作用為束水攻沙﹑穩(wěn)定航槽，隨著壩頭處流速的減小，丁壩漸漸達(dá)不到整治要求，因此認(rèn)為當(dāng)丁壩水毀體積達(dá)到壩頭總體積的30%時(shí)，丁壩失效．試驗(yàn)中對(duì)水毀體積的測(cè)量是在每一次洪水周期過后，在ti(i=1，2，3，…)時(shí)刻將每一次洪水循環(huán)后的水毀體積Vi疊加，求出=0．3時(shí)，記下在此洪峰流量下的洪水循環(huán)次數(shù)N，其中為累計(jì)水毀體積．

1．2 洪水循環(huán)可靠度模型

圖3 某一工況時(shí)的跟蹤流速Fig．3 Tracking velocity diagram under a working condition

由于散拋石壩的結(jié)構(gòu)松散度較大且壩基處泥沙運(yùn)動(dòng)的隨機(jī)性導(dǎo)致丁壩的水毀S－N關(guān)系也有明顯的隨機(jī)性．假設(shè)f(N)為在確定洪峰流量S下，丁壩水毀體積達(dá)到壩頭總體積的30%時(shí)的洪水循環(huán)次數(shù)N的密度函數(shù)，由于天然條件下丁壩周圍水沙條件復(fù)雜，壩體多以拋石為主，密實(shí)度較差，加之其工作環(huán)境一般位于水下，難以進(jìn)行實(shí)際觀測(cè)，所以丁壩受洪水沖擊作用下的失效模型究竟服從何種分布，目前還沒有充足的論證，也難以從理論上推斷究竟服從何種分布．對(duì)工況3和4(最大洪峰流量為68和100 L/s)分別進(jìn)行12次洪水循環(huán)作用下壩體破壞試驗(yàn)，得到的試驗(yàn)結(jié)果如圖4所示．

圖4曲線類似于某種指數(shù)形式或極值形式的分布函數(shù)，但由于本文中的數(shù)據(jù)是從概化水槽模型試驗(yàn)中得到的，加上樣本數(shù)量較少，而建模成敗的關(guān)鍵是要選擇恰當(dāng)?shù)膲勖植碱愋?，盲目地依賴于假設(shè)檢驗(yàn)是不可取的［4］，因此通過圖4曲線并不能準(zhǔn)確地得到洪水循環(huán)作用下丁壩的水毀失效分布模型．但考慮到曲線的形狀及工程上常采用的壽命分布形式［6］，本文借鑒損傷力學(xué)中通過定義損傷變量研究損傷演化規(guī)律來預(yù)測(cè)疲勞壽命的方法以及混凝土凍融循環(huán)作用下的失效分布模型［7］，提出一種應(yīng)用比較普遍的三參數(shù)威布爾(Weibull)分布模型對(duì)洪水循環(huán)作用下丁壩可靠度進(jìn)行探究．之所以選擇威布爾分布是由于威布爾分布對(duì)于各種類型的試驗(yàn)數(shù)據(jù)擬合的能力很強(qiáng)，在各個(gè)領(lǐng)域中有許多現(xiàn)象近似地符合于威布爾分布，它的適用范圍很廣，是在可靠性工程中廣泛使用的連續(xù)型分布形式［8］．選擇Weibull分布的另一個(gè)原因是［4］:無論指數(shù)分布還是極值分布都與Weibull分布密切相關(guān);在大量仿真試驗(yàn)的基礎(chǔ)上發(fā)現(xiàn)正態(tài)分布和三參數(shù)威布爾分布之間存在解析關(guān)系;工程中常?？梢杂猛紶柗植加行У?cái)M合對(duì)數(shù)正態(tài)數(shù)據(jù)，用三參數(shù)威布爾分布擬合正態(tài)和對(duì)數(shù)正態(tài)數(shù)據(jù)，相對(duì)誤差ε＜5%．由此可見，用三參數(shù)威布爾分布擬合試驗(yàn)數(shù)據(jù)，可以提供工程上能夠接受的統(tǒng)一結(jié)果．

通過以上的分析，本文擬采用比較普遍的三參數(shù)威布爾分布對(duì)丁壩壩體可靠度進(jìn)行推導(dǎo)，首先建立洪水循環(huán)作用下丁壩的壽命分布函數(shù):

圖4 丁壩失效時(shí)洪水循環(huán)次數(shù)－發(fā)生頻率曲線Fig．4 Flood cycles and frequency of occurrencecurve during spur dike failure

式中:X為壽命隨機(jī)變量;m為最小保證壽命;a為特征壽命或尺度參數(shù);b為形狀參數(shù)．

相應(yīng)的丁壩受洪水循環(huán)作用壽命分布密度函數(shù)可表示為

已知洪峰流量和丁壩受洪水循環(huán)作用壽命分布密度函數(shù)，則可得丁壩在洪水循環(huán)作用下的可靠度為

式中:n為給定的洪水循環(huán)次數(shù);Pr(N)為在洪水循環(huán)次數(shù)n下的可靠度．

由式(2)和(3)可知洪水條件下丁壩壩體可靠度為

1．3 不同洪峰流量時(shí)洪水循環(huán)作用下各階段丁壩失效概率

如圖2所示，在洪水歷時(shí)相同情況下，最大洪峰流量對(duì)丁壩水毀程度影響很大，當(dāng)洪峰流量為38．7 L/s時(shí)，壩體失效(水毀體積達(dá)到總體積的30%)洪水循環(huán)次數(shù)為100次，而當(dāng)洪峰流量為120 L/s時(shí)，壩體失效洪水循環(huán)次數(shù)為1次．由于在真實(shí)河道斷面上每一次洪峰流量、洪水總量等可能都各不相同，為了簡(jiǎn)化計(jì)算，假設(shè)將洪峰流量、歷時(shí)相近的洪水循環(huán)歸并為一個(gè)階段，這樣就把不同洪峰流量的洪水隨機(jī)作用下丁壩的水毀問題簡(jiǎn)化為不同洪水段的洪水順序作用下丁壩水毀失效問題．設(shè)第i階段中洪峰流量為Si，洪水循環(huán)次數(shù)為ni，在這樣的洪水循環(huán)作用下丁壩的使用壽命為Ni．為了便于分析，本文作了如下假設(shè):(1)對(duì)于任一給定的洪峰流量，在每一次洪水周期作用后，壩體水毀體積相同;(2)不同洪峰流量的洪水施加順序?qū)Χ问褂脡勖鼪]有影響．

圖5 洪水循環(huán)等效水毀關(guān)系Fig．5 Equivalent relationship of water-damage under the cyclic action by flood

由假設(shè)(1)可得，在水平Si下，洪水循環(huán)ni次產(chǎn)生的損傷因子(水毀程度)可以定義為Di=ni/Ni，當(dāng)ni=Ni時(shí)，Di=1，壩體水毀體積達(dá)到壩頭總體積的30%時(shí)，認(rèn)為丁壩失效．在Si水平洪水循環(huán)作用下，洪水循環(huán)次數(shù)與丁壩壩體水毀程度成線性關(guān)系(如圖5)．在Si水平洪水循環(huán)作用ni次，可以等效為在Sj水平下作用ni，je次洪水循環(huán)(如圖5)．

由假設(shè)(2)可知，丁壩累計(jì)損傷程度只與洪水大小有關(guān)，而與施加順序無關(guān)，雖然這與實(shí)際情況有所出入，但通過試驗(yàn)結(jié)果來看，在不同的洪水施加順序下壩體累計(jì)水毀體積相差并不大(試驗(yàn)中小于7%)，因此我們認(rèn)為此假設(shè)能夠滿足實(shí)際需要．對(duì)于丁壩來說，隨著使用年限的增加，如果沒有進(jìn)行及時(shí)修補(bǔ)，其失效概率逐年遞增，根據(jù)三參數(shù)威布爾分布的性質(zhì)可得［6］，失效概率遞增時(shí)，形狀參數(shù)b＞1．

根據(jù)文獻(xiàn)［7］對(duì)混凝土結(jié)構(gòu)在凍融作用下各階段失效概率的推導(dǎo)方法，本文嘗試對(duì)洪水循環(huán)作用下，丁壩在各階段的失效概率進(jìn)行推導(dǎo)．由失效概率和可靠度的關(guān)系可得受洪水循環(huán)作用的丁壩失效概率為

在第一階段洪水循環(huán)作用下丁壩壽命的密度函數(shù)為f1(N1)，洪水循環(huán)次數(shù)為n1，可知這時(shí)丁壩失效概率為

假設(shè)在第二階段洪水循環(huán)作用下丁壩壽命的密度函數(shù)為f2(N2)，洪水循環(huán)次數(shù)為n2，那么可根據(jù)上述丁壩水毀等效原理將第一階段洪水循環(huán)次數(shù)n1化成第二階段的等效洪水循環(huán)次數(shù)n1，2e，使得在第一階段洪水循環(huán)作用n1次壩體結(jié)構(gòu)失效概率與在第二階段洪水循環(huán)作用n1，2e次的失效概率相等，即有

因此，在相鄰的兩階段不同洪峰流量洪水循環(huán)作用下壩體結(jié)構(gòu)的失效概率為

同樣可將n1，2e+n2等效成第三階段洪水循環(huán)作用下的等效循環(huán)次數(shù)n1+2，3e，類似于式(8)，可得

在相鄰的3個(gè)階段不同洪峰流量洪水循環(huán)作用下丁壩的失效概率為

以此類推可得到在第k階段作用下前k－1階段洪水循環(huán)作用的等效循環(huán)次數(shù)為

在連續(xù)的k個(gè)階段洪水循環(huán)作用下丁壩的失效概率為

1．4 洪水作用下丁壩可靠度分析實(shí)例

通過上述理論來計(jì)算一個(gè)工程實(shí)例．某河段上一丁壩竣工后共經(jīng)歷了3次洪水作用，經(jīng)分析該丁壩經(jīng)歷的3次洪水過程可以歸結(jié)為兩階段洪水循環(huán)作用，試分析該丁壩的失效概率．本例中最小保證壽命m取0，各階段參數(shù)見表1．

表1 丁壩受洪水循環(huán)作用各階段分布函數(shù)的參數(shù)Tab．1 Parameters of distribution function in each stage under the cyclic action by flood

對(duì)于丁壩可靠度模型中的參數(shù)，由于實(shí)測(cè)資料有限，加上本文是首次提出洪水條件下丁壩的三參數(shù)Weibull可靠度模型，因此無法從其他的Weibull分布參數(shù)表中查找．根據(jù)表1的數(shù)據(jù)，這里我們選擇一種針對(duì)小樣本情況下對(duì)威布爾分布參數(shù)進(jìn)行估計(jì)的一種方法——最優(yōu)線性無偏估計(jì)法(BLUE)［9］．

此時(shí)對(duì)于式(1)中的丁壩受洪水循環(huán)作用的壽命分布函數(shù)，令X=lnN，則X服從極值分布，其分布函數(shù)為

式中:μ，σ分別為極值分布的位置參數(shù)和尺度參數(shù)，且μ=lna，σ=1/b．當(dāng)m=0時(shí)，有

利用BLUE方法對(duì)圖4中的數(shù)據(jù)進(jìn)行處理，獲得母體分布參數(shù)和可靠性指標(biāo)的估計(jì)，由BLUE無偏系數(shù)表，經(jīng)計(jì)算后可得到壩體水毀分布參數(shù)的估計(jì)值(見表1)．

對(duì)于其他洪峰流量下的洪水循環(huán)分布參數(shù)，我們可以通過今后的工作，進(jìn)一步進(jìn)行統(tǒng)計(jì)，以建立起三參數(shù)威布爾丁壩可靠度模型的參數(shù)表，進(jìn)而通過式(11)和(13)計(jì)算三階段與k階段的洪水循環(huán)作用下的丁壩失效概率．

2 經(jīng)一定洪水循環(huán)作用后丁壩剩余壽命預(yù)測(cè)

實(shí)際河道中的丁壩，在服役期間一直遭受水流沖擊作用，但并不是壩體一旦遭受洪水作用就會(huì)失去其整治能力，而是壩體先部分破壞，剩下的部分又重新形成新的穩(wěn)定結(jié)構(gòu)，仍能繼續(xù)工作和承受下一次洪水的作用，直到壩體水毀體積達(dá)到一定程度．考慮到在實(shí)際工程中維持航道整治的效果比重新對(duì)其進(jìn)行治理更為重要，因此在丁壩承受一定次數(shù)的洪水作用后，估算其還能再承受多少次洪水作用從而來預(yù)測(cè)丁壩的剩余壽命以便及時(shí)對(duì)其進(jìn)行維修，對(duì)維護(hù)航道整治效果有重要的作用．

下面以兩階段洪水循環(huán)作用為例說明丁壩剩余壽命的預(yù)測(cè)過程．設(shè)在第1階段，洪峰流量為S1下的洪水循環(huán)次數(shù)為n1，第2階段洪峰流量為S2下的洪水循環(huán)次數(shù)為n2．在這兩階段不同洪峰流量的洪水循環(huán)作用下丁壩的剩余壽命計(jì)算步驟如下:(1)計(jì)算n1相當(dāng)于S2水平下的洪水循環(huán)次數(shù)n1，2e，可由式(8)得到;(2)計(jì)算丁壩經(jīng)歷第二階段洪水循環(huán)作用后的剩余壽命nr，可寫成

對(duì)于更多階段的洪水循環(huán)情況下的剩余壽命預(yù)測(cè)的計(jì)算過程可依此類推．仍以本文1．4節(jié)的例子進(jìn)行剩余壽命預(yù)測(cè)，在經(jīng)歷了第2階段的洪水循環(huán)作用后，丁壩的剩余壽命為nr=5－2．87－2=0．13，即此丁壩還可以承受S2水平的洪水作用0．13次，也就是1次．基于目前國家水文預(yù)報(bào)已趨于成熟，可以通過對(duì)河道中每年的洪水預(yù)報(bào)來計(jì)算丁壩的剩余壽命．如預(yù)測(cè)到明年S2水平的洪水將會(huì)發(fā)生，則此丁壩剩余使用壽命為1年;如果為Si水平的洪水，則可以通過等效循環(huán)次數(shù)的計(jì)算公式，計(jì)算出壩體還可以承受Si水平的洪水循環(huán)次數(shù)．經(jīng)實(shí)地勘測(cè)，在經(jīng)過兩階段洪水循環(huán)作用后，該丁壩的水毀體積達(dá)到壩頭總體積的25．48%，而在經(jīng)過下一年一次S2水平的洪水作用后水毀體積達(dá)到壩頭總體積的32．6%，通過前文對(duì)壩體失效的定義可以看出該丁壩已經(jīng)失效，從預(yù)測(cè)結(jié)果來看，用丁壩剩余壽命計(jì)算公式得到的結(jié)果與實(shí)際情況較為相符．

3 結(jié)語

(1)借鑒損傷力學(xué)中的宏觀分析方法，基于三參數(shù)Weibull分布，建立了洪水作用下丁壩可靠度分析模型，利用水毀等效原理，推導(dǎo)出丁壩受洪水循環(huán)作用下失效概率計(jì)算式，該計(jì)算式形式簡(jiǎn)單，便于應(yīng)用．

(2)結(jié)合洪水作用下丁壩可靠度分析模型，通過論證和推導(dǎo)，得出了洪水作用下丁壩剩余壽命計(jì)算式，結(jié)合現(xiàn)代水文學(xué)中的洪水預(yù)測(cè)機(jī)制，可以得出在遭受一定洪水作用后丁壩的剩余使用年限，對(duì)維護(hù)航道有重要意義．

(3)本文通過水槽試驗(yàn)對(duì)拋石丁壩壩頭斷面平均流速與壩體水毀體積進(jìn)行跟蹤測(cè)量，首次提出以壩體水毀體積達(dá)到壩頭總體積的30%為臨界失效水毀體積，為拋石丁壩安全性的判別提供依據(jù)．

(4)對(duì)于采用Weibull分布來描述洪水作用下丁壩的失效模型是否正確還有待于通過大量的原型資料來進(jìn)一步驗(yàn)證．由于本文的數(shù)據(jù)都是從水槽試驗(yàn)中得到的，雖然模型通過一定的比尺來控制，但畢竟概化水槽無法真實(shí)模擬天然河道中的水沙及邊界條件，而且試驗(yàn)中的概化洪水過程也不能完全反映天然的洪水過程，因此接下來可以進(jìn)一步分析在實(shí)際洪水過程下丁壩的水毀情況，嘗試找出其水毀規(guī)律．

［1］王平義，程昌華，榮學(xué)文，等．航道整治建筑物水毀理論及模擬技術(shù)［M］．北京:人民交通出版社，2005．(WANG Pingyi，CHENG Chang-hua，RONG Xue-wen，et al．The water damage theory of waterway regulation structure and simulation technique［M］．Beijing:China Communications Press，2005．(in Chinese))

［2］YU Tao，WANG Ping-yi，YANG Cheng-yu，et al．Experimental study of water force acting on spur dike［J］．Advances in Civil Engineering，2011:3558-3562．

［3］XU Ge-ning，YANG Rui-gang，ZHOU Ke，et al．Methodology to estimate remaining service life of steel structure by possibilistic reliability theory［J］．Chinese Journal of Mechanical Engineering，2010，23(6):780-787．

［4］鄔華芝，郭海丁，高德平．疲勞破壞壽命的概率統(tǒng)計(jì)方法研究綜述［J］．強(qiáng)度與環(huán)境，2002，29(4):38-43．(WU Huazhi，GUO Hai-ding，GAODe-ping．Generalization of life characteristic investigation for probabilistic fatigue failure［J］．Structure＆Environment Engineering，2002，29(4):38-43．(in Chinese))

［5］曹艷敏，張華慶，蔣昌波，等．丁壩沖刷坑及下游回流區(qū)流場(chǎng)和紊動(dòng)特性試驗(yàn)研究［J］．水動(dòng)力學(xué)研究與進(jìn)展:A輯，2008，23(5):560-569．(CAO Yan-min，ZHANG Hua-qing，JIANG Chang-bo，et al．The experimental research of the flow field and turbulence characteristics in the scour and backflow region around a groin［J］．Chinese Journal of Hydrodynamics(SerA)，2008，23(5):560-569．(in Chinese))

［6］PRABHAKAR MURTHY D N，BULNER M，ECCLESTON JA．Weibull model selection for reliability modeling［J］．Reliability Engineering＆System Safety，2004，86(3):257-267．

［7］宋玉普，冀曉東．混凝土凍融損傷可靠度分析及剩余壽命預(yù)測(cè)［J］．水利學(xué)報(bào)，2006，37(3):259-263．(SONG Yu-pu，JI Xiao-dong．Analysis of reliability of concrete under freezing thawing action and evaluation of residual life［J］．Journal of Hydraulic Engineering，2006，37(3):259-263．(in Chinese))

［8］幺亮，馮蘊(yùn)雯．三參數(shù)威布爾分布疲勞壽命分散系數(shù)確定方法［J］．科學(xué)技術(shù)與工程，2009，9(6):1488-1493．(YAO Liang，F(xiàn)ENGYun-wen．Fatigue life scatter factor conformation method for sample follows the three parameters Webull distribution［J］．Science Technology and Engineering，2009，9(6):1488-1493．(in Chinese))

［9］張志華．可靠性理論及工程應(yīng)用［M］．北京:科學(xué)出版社，2012．(ZHANG Zhi-hua．Reliability theory and engineering application［M］．Beijing:Science Press，2012．(in Chinese))

［11］LEE E T，WANG J．Statistical methods for survival data analysis［M］．America:Wiley，2003．

［12］茆詩松，周紀(jì)薌，陳穎．試驗(yàn)設(shè)計(jì)［M］．北京:中國統(tǒng)計(jì)出版社，2004．(MAO Shi-song，ZHOU Ji-xiang，CHEN Ying．Experiment design［M］．Beijing:China Statistics Press，2004．(in Chinese))

［13］長(zhǎng)江航道局．航道工程手冊(cè)［M］．北京:人民交通出版社，2004．(Changjiang Waterway Bureau．Waterway engineering handbook［M］．Beijing:China Communications Press，2004．(in Chinese))