2、018,正態(tài)分布）；2.上游水位Hs=（119m,0.06,正態(tài)分布）；3.下游水位Hx=（25m,0.055,正態(tài)分布）；4.揚(yáng)壓力系數(shù)a=（0.35,0.15,極值I型分布）；5.摩擦系數(shù)=（0.72,0.25,對(duì)數(shù)正態(tài)分布）；6.凝聚力c=（600kN/m2,0.6,極值I型分布）；7.混凝土抗壓強(qiáng)度Rc=（9980kN/m2,0.23,對(duì)數(shù)正態(tài)分布）。用JC法計(jì)算該混凝土重力壩段的可靠度。作業(yè)要求：要求給出計(jì)算程序、計(jì)算步驟。1.給出求解可靠指標(biāo)的程序框圖；2.編制計(jì)算程序并計(jì)算出滿足精度要求的結(jié)果；3.試用截尾JC法計(jì)算該混凝土重力壩段的可靠
3、度并進(jìn)行分析比較。一.混凝土重力壩可靠度計(jì)算1.程序流程2.可靠度方程的建立混凝土重力壩失穩(wěn)形式有:壩段抗滑失穩(wěn)和壩段抗剪短失穩(wěn)混凝土重力壩壩段的可靠度有：（1）壩段抗滑穩(wěn)定的可靠度；（2）壩段抗壓強(qiáng)度穩(wěn)定的可靠度?，F(xiàn)計(jì)算混凝土重力壩壩段的可靠度如下，取單位寬度的壩段進(jìn)行計(jì)算。（1）壩段抗滑穩(wěn)定的可靠度1.建立壩段抗滑穩(wěn)定極限狀態(tài)方程壩段抗滑穩(wěn)定的極限狀態(tài)方程為。運(yùn)用抗剪斷強(qiáng)度公式計(jì)算抗滑穩(wěn)定的極限狀態(tài)方程為：1。豎向力（不包含）壩體自重：上游水重：下游水重：豎向力：2。揚(yáng)壓力3。上游水壓力4。下游水壓力說(shuō)明：a.計(jì)算上游水重時(shí)，近似認(rèn)為上游水位低于上
4、游折坡處，故水重按三角形計(jì)算；b.揚(yáng)壓力折減處距壩踵L。L根據(jù)SL319-2005《混凝土重力壩設(shè)計(jì)規(guī)范》屮對(duì)廊道的規(guī)定，再結(jié)合排水孔，近似取L二10】m綜上，抗滑穩(wěn)定的極限狀態(tài)方程為：2.對(duì)于非正態(tài)分布的的變量先進(jìn)行當(dāng)量正態(tài)化：在設(shè)計(jì)驗(yàn)算點(diǎn)X*對(duì)數(shù)正態(tài)分布變量f與Rc的
5、1與6標(biāo)準(zhǔn)差：均值：其中：在設(shè)計(jì)驗(yàn)算點(diǎn)x客極值I型分布變量a與c的參數(shù)a與k：則其概率密度和分布函數(shù)分別為由當(dāng)量正態(tài)化條件得1.列出Matlab程序語(yǔ)言編制的JC法程序所需相關(guān)公式2.編寫(xiě)Matlab程序語(yǔ)言編制的JC法程序，程序代碼：clear;clc;muX=l24;l19;25
6、;0.35;0.72;600J;%均值sigmaX=[0.018;0.06;0.055;0.15;0.25;0.6];%方差%將對(duì)數(shù)正態(tài)分布轉(zhuǎn)化為正態(tài)分布，求當(dāng)量正態(tài)化后的方差和均值sLn=sqrt(log(1+(sigmaX(5)/muX(5))A2));mLn=muX(5)-sLnA2/2;%將極值1型分布當(dāng)量正態(tài)化，求當(dāng)量正態(tài)化后的方差和均值aEv4=sqrt(6)*sigmaX(4)/pi;uEv4=-psi(l)*aEv4-muX(4);aEv6=sqrt(6)*sigmaX(6)/pi;uEv6=-psi(1)*aEv6-muX(6);mu
7、Xl=muX;sigmaX1=sigmaX;x=muX;normX=eps;whileabs(norm(x)-normX)/normX>le-6normX=norm(x);g=(6807.81*x(l)+0.735*x(2)A2+3.528*x(3)A2-49*x(2)-528.612*x(4)*x(2)-100&224*x(3)+52&612*x(4)*x(3))*x(5)+107.88*x(6)-4.9*x(2)A2+4.9*x(3)A2;gX=[6807.81*x(5);(l.47恢⑵?49?52&612*x(4))*x(5)-9.8*x(2);
8、(7.056*x(3)?100&224+528.612*x(4))*x(5)+9.8*x(3);(-528.612*x(2)+52&612*x(3))*x(5);6807.81*x(l)+0.735*x(2)A2+3.528*x(3)A2-49*x(2)-1008.224^x(3)-528.612*x(2)*x(4)+528.612*x(3)*x(4);107.88];%logncdf求對(duì)數(shù)分布累積分布函數(shù),evcdf求極值分布累積分布函數(shù)cdfX=[1-evcdf(-x(4),uEv4,aEv4);logncdf(x(5),mLn,sLn);1-ev
9、cdf(-x(6),uEv6,aEv6)];%gnpdf求對(duì)數(shù)分布概率密度函數(shù)，evpdf求極
