【正文】 (norminv(Fm,0,1),0,1)/fm。n2=*10^4。m22=10^4+K1*S**10^4。 m31=m3n31*norminv(Fm,0,1)。d=。VP=(m32m31)/n31。 R=m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。m2=10^4。 m12=m11+K1*n11*R。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。n2=*10^4。m22= m21+K1*S*n21。 m31=m3n31*norminv(Fm,0,1)。 m11=m1n11*norminv(Fm,0,1)。n2=*10^4。m22=m21+K0*S*n21。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。m2=10^4。 m12=2*10^7+K1**10^7*R。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。d=。VP=(m32m31)/n31。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 m21=m2n21*norminv(Fm,0,1)。n2=*10^4。m22=m21+K0*S*n21。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 c=。VP=(m32m31)/n31。 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 K0=K1。 第 2次迭代結(jié)果為 21 最后根據(jù)可靠度驗(yàn)算條件編程 while abs(K1K0)*K0 m1=m12。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 Fm=exp(exp(c*(m3d)))。 第一次迭代結(jié)束 19 第二次迭代 首先將 P的 I型變量轉(zhuǎn)為正態(tài)變量 m1=m12。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。d=k=。 α= k=。 上面介紹的驗(yàn)算點(diǎn)方法是國際安全度聯(lián)合委員會(huì)(JCSS)推薦采用的拉克維茨－菲斯勒法 (Rackwitz－ Fiessler)，所以簡(jiǎn)稱 JC法或 RF法。ln*ln*39。 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 對(duì)于非正態(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)用了 。39。39。 由條件 ① ， ? ? ? ?? ?? ?*139。 （ 縱坐標(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) 。 1i i iX i X i XX F X?? ? ??? ? ? ??由條件 ② ， ? ? ? ?? ????????????????????????? ?????????????????2exp2112exp21239。39。39。*ln* 1 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)算點(diǎn) P*的坐標(biāo)值 是 否 A 13 由于 ?是以 Z的一階原點(diǎn)矩和二階中心矩表達(dá)的，且在計(jì)算 ?時(shí)考慮了基本變量的分布類型，并采用了線性化的近似手段，因此這種結(jié)構(gòu)可靠度的計(jì)算方法通常稱為“考慮變量分布類型的一次二階矩方法”。бP=1。 c=α=。 17 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。VP=(m32m31)/n31。d=。 20 最后計(jì)算新設(shè)計(jì)驗(yàn)算點(diǎn)的均值 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。VP=(m32m31)/n31。d=。 K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。VI=(m22m2)/n2。n3=1。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 24 第二次迭加 m12=m1+K0*n1*R。m2=m22。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。VI=(m22m21)/n21。 c=。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 T=n3*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。n1=*10^7。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m12=m11+K0*n11*R。m2=m22。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m12=m11+K1*n11*R。m2=m22。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 T=n31*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。n1=*10^7。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。 K0=G/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。VI=(m22m2)/n2。 c=。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 m12=m11+K1*n11*R。m2=m22。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。 Q=Gm2*n11*VEm11*n2*VI+*n31*VP。 VE=(m12m11)/n11。n3=1。m32=m3+K0*T*n3。m3=m32。 m12=2*10^7+K1**10^7*R。m2=m22。 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n3*)^2)。VP=(m32m3)/n3。d=k=。 B=log(m1/sqrt(1+(n1/m1)^2))。 m31=m3n31*norminv(Fm,0,1)。 41 第二次迭代： m1=m12。 m11=m1*(1log(m1)+B)。 K1=(m11m2m31)/sqrt(n11^2+n2^2+n31^2)。n1=*m1。 m11=m1*(1log(m1)+B)。 K1=(m11m2m31)/sqrt(n11^2+n2^2+n31^2)。n1=*m1。n11=m1*A。 R=n11/sqrt(n11^2+n21^2+n3^2)。m2=m22。 c=。G=n21/sqrt(n11^2+n21^2+n3^2)。n2=*m2。 c=。G=n21/sqrt(n11^2+n21^2+n3^2)。n2=*m2。d=。Q=n3/sqrt(n11^2+n21^2+n3^2)。m3=m32。 Fm=exp(exp(c*(m2d)))。 m12=m11+n11*K1*R。n3=33*m3。 Fm=exp(exp(c*(m2d)))。 m12=m11+n11*K1*R。n3=*m3。 B=log(m3/sqrt(1+(n3/m3)^2))。m22=m2+G*K0*n2。 c=。 m31=m3*(1log(m3)+B)。m32=m31+n31*K1*Q。 c=。 m31=m3*(1log(m3)+B)。m32=m31+n31*K1*Q。d=。 n31=m3*A。 56 第二次迭代： m1=m12。 Fm=exp(exp(c*(m2d)))。 K1=(m1m21m31)/sqrt(n1^2+n21^2+n31^2)。n1=*m1。 Fm=exp(exp(c*(m2d)))。 K1=(m1m21m31)/sqrt(n1^2+n21^2+n31