【正文】 m12=m1+n1*K1*R。 A=。n3=3*m3。 m12=m1+n1*K1*R。 A=。m3=m32。Q=n31/sqrt(n1^2+n21^2+n31^2)。 m21=m2n21*norminv(Fm,0,1)。n2=*m2。G=n2/sqrt(n11^2+n2^2+n31^2)。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。n2=*m2。G=n2/sqrt(n11^2+n2^2+n31^2)。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。m2=m22。 R=n11/sqrt(n11^2+n2^2+n31^2)。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。n1=*m1。 K1=(m11m21m3)/sqrt(n11^2+n21^2+n3^2)。 m11=m1*(1log(m1)+B)。n1=*m1。 K1=(m11m21m3)/sqrt(n11^2+n21^2+n3^2)。 m11=m1*(1log(m1)+B)。 49 第二次迭代： m1=m12。 m21=m2n21*norminv(Fm,0,1)。 B=log(m1/sqrt(1+(n1/m1)^2))。m3=m3+n3*K1*Q。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 A=。m3=m3+n3*K1*Q。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 A=。m22=m21+G*K0*n21。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。n3=*m3。 m12=m11+n11*K1*R。 Fm=exp(exp(c*(m3d)))。n3=*m3。 m12=m11+n11*K1*R。 Fm=exp(exp(c*(m3d)))。m3=m32。Q=n31/sqrt(n11^2+n2^2+n31^2)。d=。n2=*m2。 n2=бI=*m2。m22=10^4+K1*S**10^4。 K0=K1。VI=(m22m2)/n2。 K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 m1=m12。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n3*)^2)。n1=*10^7。 T=n3*((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 c=。VI=(m22m2)/n2。 K1=Q/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。m2=m22。 m12=m1+K0*n1*R。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。n3=1。 VE=(m12m11)/n11。 Q=Gm21*n11*VEm11*n21*VI+*n31*VP。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 c=。VI=(m22m21)/n21。 K1=Q/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 c=。VI=(m22m21)/n21。 K0=G/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。n3=1。 VE=(m12m1)/n1。 Q=Gm21*n1*VEm1*n21*VI+*n31*VP。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。m2=m22。 m12=2*10^7+K1**10^7*R。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 c=。VI=(m22m21)/n21。 K0=G/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。n1=*10^7。 m12=2*10^7+K1**10^7*R。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。n2=*10^4。m22=10^4+K1*S**10^4。 m31=m3n31*norminv(Fm,0,1)。n2=*10^4。m22=m2+K0*S*n2。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。m2=цI=10^4。б E=*10^7。其計算框圖如下： 11 已知 :Xi(i=1,??n) 的分布類型及統(tǒng)計參數?xi,σ xi，極限狀態(tài)方程 g(x1??x n)=0 假定設計驗算點 P*的坐標值初值： Xi*(可取 Xi*＝ ?xi) 對于非正態(tài)變量 Xi， 根據 Xi*和公式 (1),(2)求出 ?xi′ 、 σ xi′ 以代替 ?xi、 σ xi 2121?**||cos??????????????????????????niXPiXPiXiiiXgXg???求A 12 的坐標值即求出驗算點代入坐標轉換公式：將求得的*?39。 6 ? ?? ?? ?? ? ? ?? ?iiiXiXiiiXiXiiXXyxXXxXXiysdyeyxXPxFxexxfXlnln2ln0ln2lnlnln 21 0 21 2ln2ln2ln2ln?????????????????????????令：有服從對數正態(tài)分布，則因為7 ? ?? ?? ?分布。*39。239。39。非正態(tài)分布基本變量的情況 如果極限狀態(tài)方程中的基本變量 Xi是非正態(tài)隨機變量 ， 則需首先將非正態(tài)變量在一定的條件下等效為正態(tài)變量 ， 即進行當量 （ 或等效 ） 正態(tài)化 。**39。*39。* 1 *39。態(tài)分布來計算對數正態(tài)從正因此可以利用式正態(tài)分布的分布函數，分布的分布函數化為了此式實際上將對數正態(tài)代入上式得則axdsexFdsydyiiiXiXiiXXxsXXa ln 21 lnlnln2lnlnln2??????????? ????????????????8 ? ? ? ?? ?? ?? ? ? ?? ?? ?? ? ? ?22**l n l n22l n l n39。39。цI=10^4。 n2=бI=*10^4。 m31=m3n31*norminv(Fm,0,1)。m32=m31+K0*T*n31。m3=m32。 然后求 KI（可靠度） G=m1**m3。m32=m31+K1*T*n31。m3=m32。 m31=m3n31*norminv(Fm,0,1)。m22=10^4+K1*S**10^4。m2=10^4。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。VP=(m32m31)/n31。d=。 m31=m3n31*norminv(Fm,0,1)。m22= m21+K1*S*n21。n2=*10^4。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。VI=(m22m21)/n21。 c=。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 R=m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。VP=(m32m31)/n31。d=。 m21=m2n21*norminv(Fm,0,1)。 R=m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。VP=(m32m31)/n31。d=。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 K1=Q/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。VI=(m22m21)/n21。 c=。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。m22=m2+K0*S*n2。n2=*10^4。 m11=m1n11*norminv(Fm,0,1)。 R=m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。VP=(m32m31)/n31。d=。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m12=m11+K1*n11*