【正文】 2+(n31*)^2)。 K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 然后求 KI（可靠度） G=m1**m3。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 c=。m3=m32。m2=m22。 第一次迭代結(jié)束 19 第二次迭代 首先將 P的 I型變量轉(zhuǎn)為正態(tài)變量 m1=m12。VI=(m22m2)/n2。m32=m31+K0*T*n31。 18 m12=m1+K0*n1*R。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 K0=G/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 m31=m3n31*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。d=k=。n3=бP=1。 n2=бI=*10^4。n1=б E=*10^7。 α= k=。 цP=4。цI=10^4。 14 RF法算例 1： 極限狀態(tài)方程為： =0 цE=2*10^7。 上面介紹的驗算點方法是國際安全度聯(lián)合委員會(JCSS)推薦采用的拉克維茨－菲斯勒法 (Rackwitz－ Fiessler)，所以簡稱 JC法或 RF法。39。39。 根據(jù)以上的討論 ， 對于結(jié)構極限狀態(tài)函數(shù)中包含多個正態(tài)或非正態(tài)基本變量的一般情況 ， 只要知道了各基本變量的概率分布類型及統(tǒng)計參數(shù) ，就可采用迭代法計算 ?和設計驗算點 P*的坐標值 。ln*ln*39。 * 1 * 39。態(tài)分布來計算對數(shù)正態(tài)從正因此可以利用式正態(tài)分布的分布函數(shù)，分布的分布函數(shù)化為了此式實際上將對數(shù)正態(tài)代入上式得則axdsexFdsydyiiiXiXiiXXxsXXa ln 21 lnlnln2lnlnln2??????????? ????????????????8 ? ? ? ?? ?? ?? ? ? ?? ?? ?? ? ? ?22**l n l n22l n l n39。 如果隨機變量 Xi為對數(shù)正態(tài)分布基本變量： 將對數(shù)正態(tài)分布的 Xi 直接根據(jù)當量化處理的兩個條件轉(zhuǎn)化為當量正態(tài)分布 。 1 * *1111 2iXXiiiXiiiX X Xi i iXiiiXiXX i iiXX i X i iX i i X iXXFXXf X F XF X f X?????????? ? ????????????????????? ? ???????－－－又 由 公 式 的 推 導 過 程 知＝所 以 有即 得5 對于非正態(tài)隨機變量 ， 以從公式 (1),(2) 求得的 ?Xi′ 和 ?Xi′ 分別代替 ?Xi和 ?Xi后 ， 所有的隨機變量現(xiàn)在都變成了正態(tài)分布隨機變量 ， 所以前述正態(tài)分布基本變量情況下求 ?和設計驗算點 P*的公式和方法也就均可應用了 。 39。* 1 *39。1*39。39。*iXiiXiXiiXiiXiXiiXiXXXXfXf????????4 ? ?? ?? ?? ? ? ?? ?? ?? ?? ?? ? ? ? ? ?*39。*39。239。39。 * 1 * 39。**39。*39。 由條件 ① ， ? ? ? ?? ?? ?*139。 （ 尾部面積相等 ） ?在設計驗算點 P*處非正態(tài)變量和當量正態(tài)變量的概率密度函數(shù)取值相等 。非正態(tài)分布基本變量的情況 如果極限狀態(tài)方程中的基本變量 Xi是非正態(tài)隨機變量 ， 則需首先將非正態(tài)變量在一定的條件下等效為正態(tài)變量 ， 即進行當量 （ 或等效 ） 正態(tài)化 。 1 當量正態(tài)化條件： ?在設計驗算點 P*處非正態(tài)變量和當量正態(tài)變量的概率分布函數(shù)取值相等 。 （ 縱坐標相等 ） 2 如果隨機變量 Xi為 極值 I型 分布 變量 ： 等效轉(zhuǎn)換后的當量正態(tài)隨機變量 Xi的平均值和標準差分別為 ?Xi′ 和 ?Xi′ ， 其概率分布函數(shù)和概率密度函數(shù)分別為 Fxi′ (x)和 fXi′ (x) 。39。39。* iiiiiXFXXXFXFXiXiXiXiXiXiXi???????????? ????????得3 從而求得當量正態(tài)分布的平均值 ?Xi′ 為 ? ? ? ?39。 1i i iX i X i XX F X?? ? ??? ? ? ??由條件 ② ， ? ? ? ?? ????????????????????????? ?????????????????2exp2112exp21239。*39。239。*39。39。*39。*39。 39。39。 注意： 當 X*中僅有部分基本變量為非正態(tài)分布時 ，只需將這部分基本變量當量正態(tài)化 。 6 ? ?? ?? ?? ? ? ?? ?iiiXiXiiiXiXiiXXyxXXxXXiysdyeyxXPxFxexxfXlnln2ln0ln2lnlnln 21 0 21 2ln2ln2ln2ln?????????????????????????令：有服從對數(shù)正態(tài)分布，則因為7 ? ?? ?? ?分布。 1 * **ln1*ln*ln *lnl n l n22*ln*2 ln ln 11 22 Xiiiiiiiii X i XiiXXiiiX i i X iiXXiXiXXiXXXXiiaF X f XXfXXfXeeXX????????????? ? ?????????????? ???? ? ??????????????????? ?? ?????????－－由 公 式 和? ?ln 3iX9 ? ? ? ?? ?? ?39。*ln* 1 39。ln*ln**lnln**ln1 ln ln ln 1 l n i i iiiiiiiiiiiX i X i XXiXXXiXXXi i XXi i XaX F XXXXXXXXXX????????????????? ? ????????? ? ? ??????????????? ??? ???????? ?? ? ? ???????? ? ?由 公 式 、 和 上 式? ? 410 現(xiàn)在以從公式 (3),(4) 求得的 ?Xi′ 和 ?Xi′ 分別代替 ?Xi和 ?Xi后 ， 即可將對數(shù)正態(tài)隨機變量變成了正態(tài)分布隨機變量 ， 接著可按前述公式和方法求 ?和設計驗算點 P*。其計算框圖如下： 11 已知 :Xi(i=1,??n) 的分布類型及統(tǒng)計參數(shù)?xi,σ xi，極限狀態(tài)方程 g(x1??x n)=0 假定設計驗算點 P*的坐標值初值： Xi*(可取 Xi*＝ ?xi) 對于非正態(tài)變量 Xi， 根據(jù) Xi*和公式 (1),(2)求出 ?xi′ 、 σ xi′ 以代替 ?xi、 σ xi 2121?**||cos??????????????????????????niXPiXPiXiiiXgXg???求A 12 的坐標值即求出驗算點代入坐標轉(zhuǎn)換公式：將求得的*?39。*39。?cos,.cosPX ixxixiixixiix?????????將設計驗算點 P*的坐標值 Xi*代入極限狀態(tài)方程