【正文】 注意： 當(dāng) X*中僅有部分基本變量為非正態(tài)分布時 ，只需將這部分基本變量當(dāng)量正態(tài)化 。 實際上不同的研究者提出了很多種驗算點法，它們各有優(yōu)缺點，其中我國大連理工大學(xué)的趙國藩院士也提出了一種驗算點法，計算較 JC法簡單，計算精度也很高。 Fm=exp(exp(c*(m3d)))。n1=*10^7。 T=n3*((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 Fm=exp(exp(c*(m3d)))。 end 22 經(jīng)過循環(huán)迭代結(jié)果為 23 E正態(tài)、 PI I型分布 這種組合先把 PII型分布轉(zhuǎn)為正太變量，用計算后的平均值、標(biāo)準差求可靠度，因為內(nèi)容太多 就只寫程序碼和運算結(jié)果，不詳細寫了。 m31=m3n31*norminv(Fm,0,1)。m3=m32。 S=m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 K0=K1。m22= m21+K1*S*n21。 m11=m1n11*norminv(Fm,0,1)。m32=m31+K0*T*n31。 Fm=exp(exp(c*(m2d)))。m32=m31+K1*T*n31。 m11=m1n11*norminv(Fm,0,1)。m22= m21+K1*S*n21。 m11=m1n11*norminv(Fm,0,1)。 32 第二次迭加 : m1=m12。 G=m1**m3。m3=m32。 R=m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 K0=G/sqrt((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 G=m1**m3。m3=m32。n1=б E=*m1。n11=m1*A。m2=m22。G=n2/sqrt(n11^2+n2^2+n31^2)。 c=。n2=*m2。Q=n3/sqrt(n11^2+n21^2+n3^2)。 Fm=exp(exp(c*(m2d)))。n3=33*m3。 m12=m11+n11*K1*R。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 A=。m3=m3+n3*K1*Q。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。d=。 52 第二次迭加： m1=m12。 K1=(m11m2m31)/sqrt(n11^2+n2^2+n31^2)。 Fm=exp(exp(c*(m1d)))。n1=*m1。 R=n1/sqrt(n1^2+n21^2+n31^2)。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。n2=*m2。G=n21/sqrt(n1^2+n2^21+n31^2)。 m31=m3*(1log(m3)+B)。m32=m31+n31*K1*Q。 c=。 B=log(m3/sqrt(1+(n3/m3)^2))。 m12=m11+n11*K1*R。n3=*m3。 A=。Q=n31/sqrt(n11^2+n2^2+n31^2)。n2=*m2。 c=。G=n21/sqrt(n11^2+n21^2+n3^2)。m2=m22。n11=m1*A。 K1=(m11m21m3)/sqrt(n11^2+n21^2+n3^2)。n1=*m1。 m11=m1*(1log(m1)+B)。 m21=m2n21*norminv(Fm,0,1)。m32=m31+n31*K1*Q。 A=。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。m22=m2+G*K0*n2。n3=*m3。 VE=(m12m1)/n1。 36 可靠度驗算循環(huán)迭加為： while abs(K1K0)*K0 m1=m12。m2=m22。n2=*10^4。 m31=m3n31*norminv(Fm,0,1)。 33 可靠度精度要求驗算迭加 : while abs(K1K0)*K0 m1=m12。 Fm=exp(exp(c*(m3d)))。m32=m31+K0*T*n31。d=。 R=m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 K0=K1。 S=m1*n21/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 Fm=exp(exp(c*(m1d)))。 S=m1*n21/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。d=。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。m3=m32。 G=m1**m3。 m1=m12。 m21=m2n21*norminv(Fm,0,1)。m32=m31+K1*T*n31。n3=1。 Q=Gm2*n1*VEm1*n2*VI+*n31*VP。 VE=(m12m1)/n1。m3=цP=4。*39。 39。* iiiiiXFXXXFXFXiXiXiXiXiXiXi???????????? ????????得3 從而求得當(dāng)量正態(tài)分布的平均值 ?Xi′ 為 ? ? ? ?39。*39。1*39。 根據(jù)以上的討論 ， 對于結(jié)構(gòu)極限狀態(tài)函數(shù)中包含多個正態(tài)或非正態(tài)基本變量的一般情況 ， 只要知道了各基本變量的概率分布類型及統(tǒng)計參數(shù) ，就可采用迭代法計算 ?和設(shè)計驗算點 P*的坐標(biāo)值 。n1=б E=*10^7。 18 m12=m1+K0*n1*R。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。m2=m22。 T=n3*((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 Fm=exp(exp(c*(m2d)))。 VE=(m12m1)/n1。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。n1=*10^7。 G=m1**m3。m3=4。 G=m1**m3。n3=1。 Q=Gm21*n11*VEm11*n21*VI+*n31*VP。n3=1。 G=m1**m3。m3=4。 T=n31*((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。 VE=(m12m11)/n11。 Fm=exp(exp(c*(m3d)))。 end 34 迭加結(jié)果為： 35 EIP都為正態(tài)分布 這種組合最簡單，直接利用 EIP均值、標(biāo)準差求可靠度 第一、 2次迭加： m1=2*10^7。VP=(m32m3)/n3。 VE=(m12m1)/n1。 m12=2*10^7+K1**10^7*R。m2=。G=n2/sqrt(n11^2+n2^2+n31^2)。d=。m3=m32。Q=n31/sqrt(n11^2+n2^2+n31^2)。 Fm=exp(exp(c*(m2d)))。n3=3*m3。m22=m21+G*K1*n21。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 A=sqrt(log(1+(n1/m1)^2))。m32=m3+n3*K0*Q。 m21=m2n21*norminv(Fm,0,1)。 B=。 end 51 方法 4： RI值 G正態(tài)分布 Q對數(shù)分布 這種方法先將 R、 Q轉(zhuǎn)為正態(tài)分布，在利用平均值、標(biāo)準差求可靠度 第一次迭代 : m1=。 K0=(m11m2m31)/sqrt(n11^2+n2^2+n31^2)。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。m2=m22。 R=m11/sqrt(n11^2+n2^2+n31^2)。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。n2=*m2。Q=n31/sqrt(n1^2+n21^2+n31^2)。 m21=m2n21*norminv(Fm,0,1)。m22=m21+G*K1*n21。 K0=K1。 B=。 m12=m1+n1*K0*R。m3=。 m11=m1n11*norminv(Fm,0,1)。Q=n31/sqrt(n11^2+n2^2+n31^2)。n2=*m2。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 R=n11/sqrt(n11^2+n21^2+n3^2)。m2=m22。n11=m1*A。 K0=(m11m21m3)/sqrt(n11^2+n21^2+n3^2)。 end 47 迭代結(jié)果： 48 方法 3： RI值 G對數(shù)分布 Q正態(tài)分布 這種方法先將 R、 Q轉(zhuǎn)為正態(tài)分布，在利用平均值、標(biāo)準差求可靠度 第一次迭代 : m1=。 B=。 m21=m2n21*norminv(Fm,0,1)。m32=m3+n3*K0*Q。 A=sqrt(log(1+(n1/m1)^2))。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。m22=m2+G*K1*n2。n3=*m3。 Fm=exp(exp(c*(m3d)))。m3=цP=。 G=m1**m3。 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 T=n3*((m2*n1)^2+(m1*n2)^2+(n3*)^2)。 m12=m11+K1*n11*R。d=。 R=