【正文】 2)。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。d=。n2=*10^4。VP=(m32m31)/n31。m22=10^4+K1*S**10^4。 R=m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 m31=m3n31*norminv(Fm,0,1)。 m11=m1n11*norminv(Fm,0,1)。d=。n2=*10^4。VP=(m32m31)/n31。m22=m2+K0*S*n2。 R=m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 c=。m2=10^4。VI=(m22m21)/n21。 m12=m11+K1*n11*R。 K1=Q/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。d=。n2=*10^4。VP=(m32m31)/n31。m22= m21+K1*S*n21。 R=m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 m31=m3n31*norminv(Fm,0,1)。 m21=m2n21*norminv(Fm,0,1)。 m11=m1n11*norminv(Fm,0,1)。d=。n2=*10^4。VP=(m32m31)/n31。m22=m21+K0*S*n21。 R=m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 c=。m2=10^4。VI=(m22m21)/n21。 m12=2*10^7+K1**10^7*R。 K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。d=。n2=*10^4。VP=(m32m31)/n31。m22= m21+K1*S*n21。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 m31=m3n31*norminv(Fm,0,1)。 m21=m2n21*norminv(Fm,0,1)。d=。n2=*10^4。VP=(m32m31)/n31。m22=m21+K0*S*n21。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 c=。m2=10^4。VP=(m32m31)/n31。m22=10^4+K1*S**10^4。 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 m31=m3n31*norminv(Fm,0,1)。 K0=K1。m3=m32。 第 2次迭代結(jié)果為 21 最后根據(jù)可靠度驗(yàn)算條件編程 while abs(K1K0)*K0 m1=m12。m32=m31+K1*T*n31。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 然后求 KI（可靠度） G=m1**m3。 Fm=exp(exp(c*(m3d)))。m3=m32。 第一次迭代結(jié)束 19 第二次迭代 首先將 P的 I型變量轉(zhuǎn)為正態(tài)變量 m1=m12。m32=m31+K0*T*n31。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 m31=m3n31*norminv(Fm,0,1)。d=k=。 n2=бI=*10^4。 α= k=。цI=10^4。 上面介紹的驗(yàn)算點(diǎn)方法是國際安全度聯(lián)合委員會(JCSS)推薦采用的拉克維茨－菲斯勒法 (Rackwitz－ Fiessler)，所以簡稱 JC法或 RF法。39。ln*ln*39。態(tài)分布來計(jì)算對數(shù)正態(tài)從正因此可以利用式正態(tài)分布的分布函數(shù)，分布的分布函數(shù)化為了此式實(shí)際上將對數(shù)正態(tài)代入上式得則axdsexFdsydyiiiXiXiiXXxsXXa ln 21 lnlnln2lnlnln2??????????? ????????????????8 ? ? ? ?? ?? ?? ? ? ?? ?? ?? ? ? ?22**l n l n22l n l n39。 1 * *1111 2iXXiiiXiiiX X Xi i iXiiiXiXX i iiXX i X i iX i i X iXXFXXf X F XF X f X?????????? ? ????????????????????? ? ???????－－－又 由 公 式 的 推 導(dǎo) 過 程 知＝所 以 有即 得5 對于非正態(tài)隨機(jī)變量 ， 以從公式 (1),(2) 求得的 ?Xi′ 和 ?Xi′ 分別代替 ?Xi和 ?Xi后 ， 所有的隨機(jī)變量現(xiàn)在都變成了正態(tài)分布隨機(jī)變量 ， 所以前述正態(tài)分布基本變量情況下求 ?和設(shè)計(jì)驗(yàn)算點(diǎn) P*的公式和方法也就均可應(yīng)用了 。* 1 *39。39。*39。39。**39。 由條件 ① ， ? ? ? ?? ?? ?*139。非正態(tài)分布基本變量的情況 如果極限狀態(tài)方程中的基本變量 Xi是非正態(tài)隨機(jī)變量 ， 則需首先將非正態(tài)變量在一定的條件下等效為正態(tài)變量 ， 即進(jìn)行當(dāng)量 （ 或等效 ） 正態(tài)化 。 （ 縱坐標(biāo)相等 ） 2 如果隨機(jī)變量 Xi為 極值 I型 分布 變量 ： 等效轉(zhuǎn)換后的當(dāng)量正態(tài)隨機(jī)變量 Xi的平均值和標(biāo)準(zhǔn)差分別為 ?Xi′ 和 ?Xi′ ， 其概率分布函數(shù)和概率密度函數(shù)分別為 Fxi′ (x)和 fXi′ (x) 。39。 1i i iX i X i XX F X?? ? ??? ? ? ??由條件 ② ， ? ? ? ?? ????????????????????????? ?????????????????2exp2112exp21239。239。39。*39。39。 6 ? ?? ?? ?? ? ? ?? ?iiiXiXiiiXiXiiXXyxXXxXXiysdyeyxXPxFxexxfXlnln2ln0ln2lnlnln 21 0 21 2ln2ln2ln2ln?????????????????????????令：有服從對數(shù)正態(tài)分布，則因?yàn)? ? ?? ?? ?分布。*ln* 1 39。其計(jì)算框圖如下： 11 已知 :Xi(i=1,??n) 的分布類型及統(tǒng)計(jì)參數(shù)?xi,σ xi，極限狀態(tài)方程 g(x1??x n)=0 假定設(shè)計(jì)驗(yàn)算點(diǎn) P*的坐標(biāo)值初值： Xi*(可取 Xi*＝ ?xi) 對于非正態(tài)變量 Xi， 根據(jù) Xi*和公式 (1),(2)求出 ?xi′ 、 σ xi′ 以代替 ?xi、 σ xi 2121?**||cos??????????????????????????niXPiXPiXiiiXgXg???求A 12 的坐標(biāo)值即求出驗(yàn)算點(diǎn)代入坐標(biāo)轉(zhuǎn)換公式：將求得的*?39。?cos,.cosPX ixxixiixixiix?????????將設(shè)計(jì)驗(yàn)算點(diǎn) P*的坐標(biāo)值 Xi*代入極限狀態(tài)方程 以求出 ? ? ? 0, **2*1 ?nXXXg ?|上次求出的 ?－ 上次求出的 ?| ≤ 允許誤差 ? 以本次求得的 Xi*作為下次的取用值 本次求得的 ?和 Xi*即為所求的可靠指標(biāo)和設(shè)計(jì)驗(yàn)算