【正文】 以求出 ? ? ? 0, **2*1 ?nXXXg ?|上次求出的 ?－ 上次求出的 ?| ≤ 允許誤差 ? 以本次求得的 Xi*作為下次的取用值 本次求得的 ?和 Xi*即為所求的可靠指標(biāo)和設(shè)計驗算點 P*的坐標(biāo)值 是 否 A 13 由于 ?是以 Z的一階原點矩和二階中心矩表達(dá)的，且在計算 ?時考慮了基本變量的分布類型，并采用了線性化的近似手段，因此這種結(jié)構(gòu)可靠度的計算方法通常稱為“考慮變量分布類型的一次二階矩方法”。 實際上不同的研究者提出了很多種驗算點法，它們各有優(yōu)缺點，其中我國大連理工大學(xué)的趙國藩院士也提出了一種驗算點法，計算較 JC法簡單，計算精度也很高。б E=*10^7。бI=*10^4。бP=1。 15 EI正態(tài)分布、 PI型分布 m1=цE=2*10^7。m2=цI=10^4。m3=цP=4。 c=α=。 Fm=exp(exp(c*(m3d)))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 16 G=m1**m3。 17 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 T=n31*((m2*n1)^2+(m1*n2)^2+(n31*)^2)。m22=m2+K0*S*n2。 VE=(m12m1)/n1。VP=(m32m31)/n31。n1=*10^7。n2=*10^4。n3=1。d=。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m31=m3n31*norminv(Fm,0,1)。 Q=Gm2*n1*VEm1*n2*VI+*n31*VP。 20 最后計算新設(shè)計驗算點的均值 R=m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 T=n3*((m2*n1)^2+(m1*n2)^2+(n31*)^2)。m22=10^4+K1*S**10^4。 VE=(m12m1)/n1。VP=(m32m31)/n31。n1=*10^7。n2=*10^4。n3=1。d=。 Fm=exp(exp(c*(m3d)))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 G=m1**m3。 K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*)^2)。 m12=2*10^7+K1**10^7*R。m32=m31+K1*T*n31。VI=(m22m2)/n2。 end 22 經(jīng)過循環(huán)迭代結(jié)果為 23 E正態(tài)、 PI I型分布 這種組合先把 PII型分布轉(zhuǎn)為正太變量，用計算后的平均值、標(biāo)準(zhǔn)差求可靠度，因為內(nèi)容太多 就只寫程序碼和運算結(jié)果，不詳細(xì)寫了。n1=*10^7。n2=*10^4。n3=1。d=。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 m21=m2n21*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m31=m3n31*norminv(Fm,0,1)。 K0=G/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 S=m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 24 第二次迭加 m12=m1+K0*n1*R。m32=m31+K0*T*n31。VI=(m22m21)/n21。 m1=m12。m2=m22。m3=m32。 c=。 Fm=exp(exp(c*(m2d)))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 G=m1**m3。 K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 S=m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 m12=2*10^7+K1**10^7*R。m32=m31+K1*T*n31。VI=(m22m21)/n21。 25 循環(huán)迭加 while abs(K1K0)*K0 m1=m12。m2=m22。m3=m32。 c=。 K0=K1。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 m21=m2n21*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m31=m3n31*norminv(Fm,0,1)。 Q=Gm21*n1*VEm1*n21*VI+*n31*VP。 R=m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*)^2)。 T=n3*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。m22= m21+K1*S*n21。 VE=(m12m1)/n1。VP=(m32m31)/n31。n1=*10^7。n2=*10^4。n3=1。d=。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。 m11=m1n11*norminv(Fm,0,1)。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 m21=m2n21*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m31=m3n31*norminv(Fm,0,1)。 K0=G/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 S=m1*n21/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 m12=m11+K0*n11*R。m32=m31+K0*T*n31。VI=(m22m21)/n21。 28 第二次迭加： m1=m12。m2=m22。m3=m32。 c=。 Fm=exp(exp(c*(m1d)))。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m2d)))。 n21= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 G=m1**m3。 K1=Q/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 S=m1*n21/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 m12=m11+K1*n11*R。m32=m31+K1*T*n31。VI=(m22m21)/n21。 29 可靠度精度要求循環(huán)迭加： while abs(K1K0)*K0 m1=m12。m2=m22。m3=m32。 c=。 K0=K1。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。 m11=m1n11*norminv(Fm,0,1)。 fm=c*(exp(c*(m2d))*exp(exp(c*(m2d))))。 m21=m2n21*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m31=m3n31*norminv(Fm,0,1)。 Q=Gm21*n11*VEm11*n21*VI+*n31*VP。 R=m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*)^2)。 T=n31*((m21*n1)^2+(m1*n21)^2+(n31*)^2)。m22= m21+K1*S*n21。 VE=(m12m11)/n11。VP=(m32m31)/n31。n1=*10^7。n2=*10^4。n3=1。d=。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))))。 m11=m1n11*norminv(Fm,0,1)。 fm=c*(exp(c*(m3d))*exp(exp(c*(m3d))))。 m31=m3n31*norminv(Fm,0,1)。 K0=G/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 S=m1*n2/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 m12=m1+K0*n1*R。m32=m31+K0*T*n31。VI=(m22m2)/n2。 32 第二次迭加 : m1=m12。m2=m22。m3=m32。 c=。 Fm=exp(exp(c*(m1d)))。 n11= normpdf(norminv(Fm,0,1),0,1)/fm。 Fm=exp(exp(c*(m3d)))。 n31= normpdf(norminv(Fm,0,1),0,1)/fm。 G=m1**m3。 K1=Q/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 S=m1*n2/sqrt((m2*n11)^2+(m11*n2)^2+(n31*)^2)。 m12=m11+K1*n11*R。m32=m31+K1*T*n31。VI=(m22m2)/n2。 33 可靠度精度要求驗算迭加 : while abs(K1K0)*K0 m1=m12。m2=m22。m3=m32。 c=。 K0=K1。 fm=c*(exp(c*(m1d))*exp(exp(c*(m1d))