【正文】 RF法算例 2： 極限狀態(tài)方程為： RSGSQ =0 假定 :m1=цE=。VI=(m22m2)/n2。m32=m3+K1*T*n3。 m12=2*10^7+K1**10^7*R。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 G=m1**m3。n3=1。n2=*10^4。n1=*10^7。VP=(m32m3)/n3。 VE=(m12m1)/n1。m22=10^4+K1*S**10^4。 T=n3*((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 Q=Gm2*n1*VEm1*n2*VI+*n3*VP。n3=1。n2=*10^4。n1=*10^7。VP=(m32m3)/n3。 VE=(m12m1)/n1。m22=m2+K0*S*n2。 T=n3*((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 G=m1**m3。m3=4。m2=10^4。 end 34 迭加結(jié)果為： 35 EIP都為正態(tài)分布 這種組合最簡(jiǎn)單，直接利用 EIP均值、標(biāo)準(zhǔn)差求可靠度 第一、 2次迭加： m1=2*10^7。VI=(m22m2)/n2。m32=m31+K1*T*n31。 m12=m11+K1*n11*R。 S=m1*n2/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 K1=Q/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 G=m1**m3。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m1d)))。d=。n3=1。n2=*10^4。n1=*10^7。VP=(m32m31)/n31。 VE=(m12m11)/n11。m22=10^4+K1*S**10^4。 T=n3*((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 R=m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 Q=Gm2*n11*VEm11*n2*VI+*n31*VP。 m31=m3n31*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m11=m1n11*norminv(Fm,0,1)。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。d=。n3=1。n2=*10^4。n1=*10^7。VP=(m32m31)/n31。 VE=(m12m11)/n11。m22=m2+K0*S*n2。 T=n31*((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 R=m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 G=m1**m3。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m1d)))。 c=。m3=4。m2=10^4。 end 30 迭加結(jié)果為： 31 E、 P為 I型分布 .I 為正態(tài)分布 這種組合應(yīng)該先把 E、 P都轉(zhuǎn)化為正態(tài)分布，在利用均值、標(biāo)準(zhǔn)差求可靠度 第一次迭加： m1=2*10^7。VI=(m22m21)/n21。m32=m31+K1*T*n31。 m12=m11+K1*n11*R。 S=m1*n21/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 K1=Q/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 G=m1**m3。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m2d)))。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m1d)))。d=。n3=1。n2=*10^4。n1=*10^7。VP=(m32m31)/n31。 VE=(m12m11)/n11。m22= m21+K1*S*n21。 T=n31*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 R=m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 Q=Gm21*n11*VEm11*n21*VI+*n31*VP。 m31=m3n31*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m21=m2n21*norminv(Fm,0,1)。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 m11=m1n11*norminv(Fm,0,1)。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。d=。n3=1。n2=*10^4。n1=*10^7。VP=(m32m31)/n31。 VE=(m12m1)/n11。m22=m21+K0*S*n21。 T=n31*((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 R=m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 G=m1**m3。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m2d)))。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m1d)))。 c=。m3=4。m2=10^4。 end 26 最終運(yùn)行結(jié)果為 : 27 EIP都為 I型分布 這種組合應(yīng)該先把 EIP都轉(zhuǎn)化為正態(tài)分布，在利用均值、標(biāo)準(zhǔn)差求可靠度 第一次迭加： m1=2*10^7。VI=(m22m21)/n21。m32=m31+K1*T*n31。 m12=2*10^7+K1**10^7*R。 S=m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 G=m1**m3。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m2d)))。d=。n3=1。n2=*10^4。n1=*10^7。VP=(m32m31)/n31。 VE=(m12m1)/n1。m22= m21+K1*S*n21。 T=n3*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 Q=Gm21*n1*VEm1*n21*VI+*n31*VP。 m31=m3n31*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m21=m2n21*norminv(Fm,0,1)。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。d=。n3=1。n2=*10^4。n1=*10^7。VP=(m32m31)/n31。 VE=(m12m1)/n1。m22=m21+K0*S*n21。 T=n31*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 G=m1**m3。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m2d)))。 c=。m3=4。m2=10^4。 m1=2*10^7。VP=(m32m31)/n31。 VE=(m12m1)/n1。m22=10^4+K1*S**10^4。 T=n3*((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 Q=Gm2*n1*VEm1*n2*VI+*n31*VP。 m31=m3n31*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 K0=K1。 c=。m3=m32。m2=m22。 第 2次迭代結(jié)果為 21 最后根據(jù)可靠度驗(yàn)算條件編程 while abs(K1K0)*K0 m1=m12。VI=(m22m2)/n2。m32=m31+K1*T*n31。 m12=2*10^7+K1**10^7*R。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^