工程结构可靠度中非正态分布转为正态分布课件

1非正态分布基本变量的情况如果极限状态方程中的基本变量Xi是非正态随机变量，则需首先将非正态变量在一定的条件下等效为正态变量，即进行当量（或等效）正态化。11非正态分布基本变量的情况如果极限状态方程中的基本变量Xi是2

当量正态化条件：在设计验算点P*处非正态变量和当量正态变量的概率分布函数取值相等。（尾部面积相等）在设计验算点P*处非正态变量和当量正态变量的概率密度函数取值相等。（纵坐标相等）22当量正态化条件：23如果随机变量Xi为极值I型分布变量：等效转换后的当量正态随机变量Xi的平均值和标准差分别为Xi′和Xi′，其概率分布函数和概率密度函数分别为Fxi′(x)和fXi′(x)

。由条件①，33如果随机变量Xi为极值I型分布变量：等效转换后的当量正态4从而求得当量正态分布的平均值Xi′为由条件②，44从而求得当量正态分布的平均值Xi′为由条件②，455556对于非正态随机变量，以从公式(1),(2)求得的Xi′和Xi′分别代替Xi和Xi后，所有的随机变量现在都变成了正态分布随机变量，所以前述正态分布基本变量情况下求和设计验算点P*的公式和方法也就均可应用了。注意：

当X*中仅有部分基本变量为非正态分布时，只需将这部分基本变量当量正态化。如果随机变量Xi为对数正态分布基本变量：将对数正态分布的Xi直接根据当量化处理的两个条件转化为当量正态分布。66对于非正态随机变量，以从公式(1),(2)求得的7777888899991010101011现在以从公式(3),(4)求得的Xi′和Xi′分别代替Xi和Xi后，即可将对数正态随机变量变成了正态分布随机变量，接着可按前述公式和方法求和设计验算点P*。根据以上的讨论，对于结构极限状态函数中包含多个正态或非正态基本变量的一般情况，只要知道了各基本变量的概率分布类型及统计参数，就可采用迭代法计算和设计验算点P*的坐标值。其计算框图如下：1111现在以从公式(3),(4)求得的Xi′和Xi12已知:Xi(i=1,……n)的分布类型及统计参数xi,σxi，极限状态方程g(x1……xn)=0假定设计验算点P*的坐标值初值：Xi*(可取Xi*＝xi)对于非正态变量Xi，根据Xi*和公式(1),(2)求出xi′、σxi′以代替xi、σxiA1212已知:Xi(i=1,……n)的分布类型及统计参数xi,13将设计验算点P*的坐标值Xi*代入极限状态方程以求出|上次求出的－上次求出的|≤允许误差以本次求得的Xi*作为下次的取用值本次求得的和Xi*即为所求的可靠指标和设计验算点P*的坐标值是否A1313将设计验算点P*的坐标值Xi*代入极限状态方程14由于是以Z的一阶原点矩和二阶中心矩表达的，且在计算时考虑了基本变量的分布类型，并采用了线性化的近似手段，因此这种结构可靠度的计算方法通常称为“考虑变量分布类型的一次二阶矩方法”。上面介绍的验算点方法是国际安全度联合委员会(JCSS)推荐采用的拉克维茨－菲斯勒法(Rackwitz－Fiessler)，所以简称JC法或R-F法。实际上不同的研究者提出了很多种验算点法，它们各有优缺点，其中我国大连理工大学的赵国藩院士也提出了一种验算点法，计算较JC法简单，计算精度也很高。1414由于是以Z的一阶原点矩和二阶中心矩表达的，且在计算时15RF法算例1：极限状态方程为：EI-78.12P=0цE=2*10^7;бE=0.5*10^7;цI=10^-4;бI=0.2*10^-4;цP=4;бP=1;α=1.2825k=3.5499;1515RF法算例1：极限状态方程为：EI-78.12P=0161、EI正态分布、PI型分布m1=цE=2*10^7;n1=бE=0.5*10^7;m2=цI=10^-4;n2=бI=0.2*10^-4;m3=цP=4;n3=бP=1;c=α=1.2825;d=k=3.5499;Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);16161、EI正态分布、PI型分布m1=цE=2*10^7;n17G=m1*m2-78.12*m3;K0=G/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);1717G=m1*m2-78.12*m3;K0=G/sqrt((18R=-m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);T=n31*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);1818R=-m2*n1/sqrt((m2*n1)^2+(m1*19m12=m1+K0*n1*R;m22=m2+K0*S*n2;m32=m31+K0*T*n31;VE=(m12-m1)/n1;VI=(m22-m2)/n2;VP=(m32-m31)/n31;第一次迭代结束1919m12=m1+K0*n1*R;m22=m2+K0*S*n20第二次迭代首先将P的I型变量转为正态变量m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);然后求KI（可靠度）G=m1*m2-78.12*m3;Q=G-m2*n1*VE-m1*n2*VI+78.12*n31*VP;K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);2020第二次迭代首先将P的I型变量转为正态变量2021最后计算新设计验算点的均值R=-m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);S=-m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);m12=2*10^7+K1*0.5*10^7*R;m22=10^-4+K1*S*0.2*10^-4;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m2)/n2;VP=(m32-m31)/n31;第1、2次迭代结果为2121最后计算新设计验算点的均值R=-m2*n1/sqrt((22最后根据可靠度验算条件编程whileabs(K1-K0)>0.003*K0m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;K0=K1;Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m2*n1*VE-m1*n2*VI+78.12*n31*VP;K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);R=-m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);S=-m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);m12=2*10^7+K1*0.5*10^7*R;m22=10^-4+K1*S*0.2*10^-4;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m2)/n2;VP=(m32-m31)/n31;end2222最后根据可靠度验算条件编程2223经过循环迭代结果为2323经过循环迭代结果为23242、E正态、PII型分布这种组合先把PII型分布转为正太变量，用计算后的平均值、标准差求可靠度，因为内容太多就只写程序码和运算结果，不详细写了。m1=2*10^7;n1=0.5*10^7;m2=10^-4;n2=0.2*10^-4;m3=4;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;K0=G/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);R=-m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);S=-m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);T=n31*78.12/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);24242、E正态、PII型分布2425第二次迭加m12=m1+K0*n1*R;m22=m21+K0*S*n21;m32=m31+K0*T*n31;VE=(m12-m1)/n1;VI=(m22-m21)/n21;VP=(m32-m31)/n31;m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m21*n1*VE-m1*n21*VI+78.12*n31*VP;K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);R=-m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);S=-m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);m12=2*10^7+K1*0.5*10^7*R;m22=m21+K1*S*n21;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m21)/n21;VP=(m32-m31)/n31;2525第二次迭加2526循环迭加whileabs(K1-K0)>0.003*K0m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;K0=K1;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m21*n1*VE-m1*n21*VI+78.12*n31*VP;K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);R=-m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);S=-m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);m12=2*10^7+K1*0.5*10^7*R;m22=m21+K1*S*n21;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m21)/n21;VP=(m32-m31)/n31;end2626循环迭加2627最终运行结果为:2727最终运行结果为:27283、EIP都为I型分布这种组合应该先把EIP都转化为正态分布，在利用均值、标准差求可靠度第一次迭加：m1=2*10^7;n1=0.5*10^7;m2=10^-4;n2=0.2*10^-4;m3=4;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m1-d)));fm=c*(exp(-c*(m1-d))*exp(-exp(-c*(m1-d))));n11=normpdf(norminv(Fm,0,1),0,1)/fm;m11=m1-n11*norminv(Fm,0,1);Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;K0=G/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);R=-m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);S=-m1*n21/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);T=n31*78.12/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);m12=m11+K0*n11*R;m22=m21+K0*S*n21;m32=m31+K0*T*n31;VE=(m12-m1)/n11;VI=(m22-m21)/n21;VP=(m32-m31)/n31;28283、EIP都为I型分布2829第二次迭加：m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m1-d)));fm=c*(exp(-c*(m1-d))*exp(-exp(-c*(m1-d))));n11=normpdf(norminv(Fm,0,1),0,1)/fm;m11=m1-n11*norminv(Fm,0,1);Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m21*n11*VE-m11*n21*VI+78.12*n31*VP;K1=Q/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);R=-m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);S=-m1*n21/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);T=n31*78.12/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);m12=m11+K1*n11*R;m22=m21+K1*S*n21;m32=m31+K1*T*n31;VE=(m12-m11)/n11;VI=(m22-m21)/n21;VP=(m32-m31)/n31;2929第二次迭加：2930可靠度精度要求循环迭加：whileabs(K1-K0)>0.003*K0m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;K0=K1;Fm=exp(-exp(-c*(m1-d)));fm=c*(exp(-c*(m1-d))*exp(-exp(-c*(m1-d))));n11=normpdf(norminv(Fm,0,1),0,1)/fm;m11=m1-n11*norminv(Fm,0,1);Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m21*n11*VE-m11*n21*VI+78.12*n31*VP;K1=Q/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);R=-m2*n11/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);S=-m1*n21/sqrt((m21*n11)^2+(m11*n21)^2+(n31*78.12)^2);T=n31*78.12/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);m12=m11+K1*n11*R;m22=m21+K1*S*n21;m32=m31+K1*T*n31;VE=(m12-m11)/n11;VI=(m22-m21)/n21;VP=(m32-m31)/n31;end3030可靠度精度要求循环迭加：3031迭加结果为：3131迭加结果为：31324、E、P为I型分布.I为正态分布这种组合应该先把E、P都转化为正态分布，在利用均值、标准差求可靠度第一次迭加：m1=2*10^7;n1=0.5*10^7;m2=10^-4;n2=0.2*10^-4;m3=4;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m1-d)));fm=c*(exp(-c*(m1-d))*exp(-exp(-c*(m1-d))));n11=normpdf(norminv(Fm,0,1),0,1)/fm;m11=m1-n11*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;K0=G/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);R=-m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);S=-m1*n2/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);T=n31*78.12/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);m12=m1+K0*n1*R;m22=m2+K0*S*n2;m32=m31+K0*T*n31;VE=(m12-m11)/n11;VI=(m22-m2)/n2;VP=(m32-m31)/n31;32324、E、P为I型分布.I为正态分布3233

第二次迭加:m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m1-d)));fm=c*(exp(-c*(m1-d))*exp(-exp(-c*(m1-d))));n11=normpdf(norminv(Fm,0,1),0,1)/fm;m11=m1-n11*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m2*n11*VE-m11*n2*VI+78.12*n31*VP;K1=Q/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);R=-m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);S=-m1*n2/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);m12=m11+K1*n11*R;m22=10^-4+K1*S*0.2*10^-4;m32=m31+K1*T*n31;VE=(m12-m11)/n11;VI=(m22-m2)/n2;VP=(m32-m31)/n31;3333第二次迭加:3334

可靠度精度要求验算迭加:whileabs(K1-K0)>0.003*K0m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;K0=K1;Fm=exp(-exp(-c*(m1-d)));fm=c*(exp(-c*(m1-d))*exp(-exp(-c*(m1-d))));n11=normpdf(norminv(Fm,0,1),0,1)/fm;m11=m1-n11*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m2*n11*VE-m11*n2*VI+78.12*n31*VP;K1=Q/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);R=-m2*n11/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);S=-m1*n2/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m2*n11)^2+(m11*n2)^2+(n31*78.12)^2);m12=m11+K1*n11*R;m22=10^-4+K1*S*0.2*10^-4;m32=m31+K1*T*n31;VE=(m12-m11)/n11;VI=(m22-m2)/n2;VP=(m32-m31)/n31;end3434可靠度精度要求验算迭加:3435迭加结果为：3535迭加结果为：35365、EIP都为正态分布这种组合最简单，直接利用EIP均值、标准差求可靠度第一、2次迭加：m1=2*10^7;n1=0.5*10^7;m2=10^-4;n2=0.2*10^-4;m3=4;n3=1;G=m1*m2-78.12*m3;K0=G/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);R=-m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);S=-m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);T=n3*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);m12=m1+K0*n1*R;m22=m2+K0*S*n2;m32=m3+K0*T*n3;VE=(m12-m1)/n1;VI=(m22-m2)/n2;VP=(m32-m3)/n3;m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;G=m1*m2-78.12*m3;Q=G-m2*n1*VE-m1*n2*VI+78.12*n3*VP;K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);R=-m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);S=-m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);T=n3*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);m12=2*10^7+K1*0.5*10^7*R;m22=10^-4+K1*S*0.2*10^-4;m32=m3+K1*T*n3;VE=(m12-m1)/n1;VI=(m22-m2)/n2;VP=(m32-m3)/n3;36365、EIP都为正态分布3637可靠度验算循环迭加为：whileabs(K1-K0)>0.003*K0m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;K0=K1;G=m1*m2-78.12*m3;Q=G-m2*n1*VE-m1*n2*VI+78.12*n3*VP;K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);R=-m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);S=-m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);T=n3*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n3*78.12)^2);m12=2*10^7+K1*0.5*10^7*R;m22=10^-4+K1*S*0.2*10^-4;m32=m3+K1*T*n3;VE=(m12-m1)/n1;VI=(m22-m2)/n2;VP=(m32-m3)/n3;end3737可靠度验算循环迭加为：3738迭代结果为：3838迭代结果为：3839

RF法算例2：极限状态方程为：R-SG-SQ=0假定:m1=цE=373.92;n1=бE=0.154*m1;m2=цI=127.2;n2=бI=0.07*m2;m3=цP=51.54;n3=бP=0.233*m3;c=α=1.2825;d=k=3.5499;方法一:R对数G正态分布QI值分布m1=373.92;n1=0.154*m1;m2=127.2;n2=0.07*m2;m3=51.54;n3=0.233*m3;A=sqrt(log(1+(n1/m1)^2));B=log(m1/sqrt(1+(n1/m1)^2));m11=m1*(1-log(m1)+B);n11=m1*A;3939RF法算例2：极限状态方程为：R-SG-SQ=040

c=1.2825/n3;d=m3-0.5772/c;Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);4040c=1.2825/n3;d=m3-0.5772/c;441

K0=(m11-m2-m31)/sqrt(n11^2+n2^2+n31^2);R=n11/sqrt(n11^2+n2^2+n31^2);G=n2/sqrt(n11^2+n2^2+n31^2);Q=n31/sqrt(n11^2+n2^2+n31^2);m12=m11+n11*K0*R;m22=m2+G*K0*n2;m32=m31+n31*K0*Q;4141K0=(m11-m2-m31)/sqrt(n11^2+42

第二次迭代：m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;A=0.1531;B=5.9123;m11=m1*(1-log(m1)+B);n11=m1*A;c=0.1068;d=46.1355;Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);K1=(m11-m2-m31)/sqrt(n11^2+n2^2+n31^2);R=n11/sqrt(n11^2+n2^2+n31^2);G=n2/sqrt(n11^2+n2^2+n31^2);Q=n31/sqrt(n11^2+n2^2+n31^2);m12=m11+n11*K1*R;m22=m2+G*K1*n2;m32=m31+n31*K1*Q;4242第二次迭代：4243

可靠度精度验算循环程序：whileabs(K1-K0)>(1e-1)*K0m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;K0=K1;A=0.1531;B=5.9123;m11=m1*(1-log(m1)+B);n11=m1*A;c=0.1068;d=46.1355;Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);K1=(m11-m2-m31)/sqrt(n11^2+n2^2+n31^2);R=n11/sqrt(n11^2+n2^2+n31^2);G=n2/sqrt(n11^2+n2^2+n31^2);Q=n31/sqrt(n11^2+n2^2+n31^2);m12=m11+n11*K1*R;m22=m2+G*K1*n2;m32=m31+n31*K1*Q;end4343可靠度精度验算循环程序：4344

迭代结果：4444迭代结果：4445

方法2：R对数GI值分布Q正态布这种方法先将R、G转为正太分布，在利用平均值、标准差求可靠度第一次迭代:m1=373.92;n1=0.154*m1;m2=127.2;n2=0.07*m2;m3=51.54;n3=0.233*m3;A=sqrt(log(1+(n1/m1)^2));B=log(m1/sqrt(1+(n1/m1)^2));m11=m1*(1-log(m1)+B);n11=m1*A;c=1.2825/n2;d=m2-0.5772/c;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);K0=(m11-m21-m3)/sqrt(n11^2+n21^2+n3^2);R=n11/sqrt(n11^2+n21^2+n3^2);G=n21/sqrt(n11^2+n21^2+n3^2);Q=n3/sqrt(n11^2+n21^2+n3^2);m12=m11+n11*K0*R;m22=m21+G*K0*n21;m32=m3+n3*K0*Q;4545方法2：R对数GI值分布Q正态布4546

第二次迭代m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;A=0.1531;B=5.9123;m11=m1*(1-log(m1)+B);n11=m1*A;c=0.1440;d=123.1927;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);K1=(m11-m21-m3)/sqrt(n11^2+n21^2+n3^2);R=n11/sqrt(n11^2+n21^2+n3^2);G=n21/sqrt(n11^2+n21^2+n3^2);Q=n3/sqrt(n11^2+n21^2+n3^2);m12=m11+n11*K1*R;m22=m21+G*K1*n21;m3=m3+n3*K1*Q;4646第二次迭代4647

可靠度验算迭代whileabs(K1-K0)>(1e-1)*K0m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;K0=K1;A=0.1531;B=5.9123;m11=m1*(1-log(m1)+B);n11=m1*A;c=0.1440;d=123.1927;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);K1=(m11-m21-m3)/sqrt(n11^2+n21^2+n3^2);R=n11/sqrt(n11^2+n21^2+n3^2);G=n21/sqrt(n11^2+n21^2+n3^2);Q=n3/sqrt(n11^2+n21^2+n3^2);m12=m11+n11*K1*R;m22=m21+G*K1*n21;m3=m3+n3*K1*Q;end4747可靠度验算迭代4748

迭代结果：4848迭代结果：4849

方法3：RI值G对数分布Q正态分布这种方法先将R、Q转为正态分布，在利用平均值、标准差求可靠度第一次迭代:m1=373.92;n1=0.154*m1;m2=127.2;n2=0.07*m2;m3=51.54;n3=0.233*m3;A=sqrt(log(1+(n1/m1)^2));B=log(m1/sqrt(1+(n1/m1)^2));m11=m1*(1-log(m1)+B);n11=m1*A;c=1.2825/n2;d=m2-0.5772/c;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);K0=(m11-m21-m3)/sqrt(n11^2+n21^2+n3^2);R=n11/sqrt(n11^2+n21^2+n3^2);G=n21/sqrt(n11^2+n21^2+n3^2);Q=n3/sqrt(n11^2+n21^2+n3^2);m12=m11+n11*K0*R;m22=m21+G*K0*n21;m32=m3+n3*K0*Q;4949方法3：RI值G对数分布Q正态分布4950

第二次迭代：m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;A=0.1531;B=5.9123;m11=m1*(1-log(m1)+B);n11=m1*A;c=0.1440;d=123.1927;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);K1=(m11-m21-m3)/sqrt(n11^2+n21^2+n3^2);R=n11/sqrt(n11^2+n21^2+n3^2);G=n21/sqrt(n11^2+n21^2+n3^2);Q=n3/sqrt(n11^2+n21^2+n3^2);m12=m11+n11*K1*R;m22=m21+G*K1*n21;m3=m3+n3*K1*Q;5050第二次迭代：5051

验算迭代：whileabs(K1-K0)>(1e-1)*K0m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;K0=K1;A=0.1531;B=5.9123;m11=m1*(1-log(m1)+B);n11=m1*A;c=0.1440;d=123.1927;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);K1=(m11-m21-m3)/sqrt(n11^2+n21^2+n3^2);R=n11/sqrt(n11^2+n21^2+n3^2);G=n21/sqrt(n11^2+n21^2+n3^2);Q=n3/sqrt(n11^2+n21^2+n3^2);m12=m11+n11*K1*R;m22=m21+G*K1*n21;m3=m3+n3*K1*Q;end5151验算迭代：5152

方法4：RI值G正态分布Q对数分布这种方法先将R、Q转为正态分布，在利用平均值、标准差求可靠度第一次迭代:m1=373.92;n1=0.154*m1;m2=127.2;n2=0.07*m2;m3=51.54;n3=0.233*m3;c=1.2825/n1;d=m1-0.5772/c;Fm=exp(-exp(-c*(m1-d)));fm=c*(exp(-c*(m1-d))*exp(-exp(-c*(m1-d))));n11=normpdf(norminv(Fm,0,1),0,1)/fm;m11=m1-n11*norminv(Fm,0,1);A=sqrt(log(1+(n3/m3)^2));B=log(m3/sqrt(1+(n3/m3)^2));m31=m3*(1-log(m3)+B);n31=m3*A;K0=(m11-m2-m31)/sqrt(n11^2+n2^2+n31^2);R=n11/sqrt(n11^2+n2^2+n31^2);G=n2/sqrt(n11^2+n2^2+n31^2);Q=n31/sqrt(n11^2+n2^2+n31^2);m12=m11+n11*K0*R;m22=m2+G*K0*n2;m32=m31+n31*K0*Q;5252方法4：RI值G正态分布Q对数分布5253

第二次迭加：m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;c=0.0223;d=348.004;Fm=exp(-exp(-c*(m1-d)));fm=c*(exp(-c*(m1-d))*exp(-exp(-c*(m1-d))));n11=normpdf(norminv(Fm,0,1),0,1)/fm;m11=m1-n11*norminv(Fm,0,1);A=0.2299;B=3.9159;m31=m3*(1-log(m3)+B);n31=m3*A;K1=(m11-m2-m31)/sqrt(n11^2+n2^2+n31^2);R=n11/sqrt(n11^2+n2^2+n31^2);G=n2/sqrt(n11^2+n2^2+n31^2);Q=n31/sqrt(n11^2+n2^2+n31^2);m12=m11+n11*K1*R;m22=m2+G*K1*n2;m32=m31+n31*K1*Q;5353第二次迭加：5354

可靠度验算迭加：whileabs(K1-K0)>(1e-1)*K0m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;K0=K1;c=0.0223;d=348.004;Fm=exp(-exp(-c*(m1-d)));fm=c*(exp(-c*(m1-d))*exp(-exp(-c*(m1-d))));n11=normpdf(norminv(Fm,0,1),0,1)/fm;m11=m1-n11*norminv(Fm,0,1);A=0.2299;B=3.9159;m31=m3*(1-log(m3)+B);n31=m3*A;K1=(m11-m2-m31)/sqrt(n11^2+n2^2+n31^2);R=m11/sqrt(n11^2+n2^2+n31^2);G=n2/sqrt(n11^2+n2^2+n31^2);Q=n31/sqrt(n11^2+n2^2+n31^2);m12=m11+n11*K1*R;m22=m2+G*K1*n2;m32=m31+n31*K1*Q;end5454可靠度验算迭加：5455

迭加结果：5555迭加结果：5556

方法5：R正太分布GI值分布Q对数分布这种方法先将G、Q转为正态分布，在利用平均值、标准差求可靠度第一次迭代:m1=373.92;n1=0.154*m1;m2=127.2;n2=0.07*m2;m3=51.54;n3=0.233*m3;c=1.2825/n2;d=m2-0.5772/c;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);A=sqrt(log(1+(n3/m3)^2));B=log(m3/sqrt(1+(n3/m3)^2));m31=m3*(1-log(m3)+B);n31=m3*A;K0=(m1-m21-m31)/sqrt(n1^2+n21^2+n31^2);R=n1/sqrt(n1^2+n21^2+n31^2);G=n21/sqrt(n1^2+n21^2+n31^2);Q=n31/sqrt(n1^2+n21^2+n31^2);m12=m1+n1*K0*R;m22=m21+G*K0*n21;m32=m31+n31*K0*Q;5656方法5：R正太分布GI值分布Q对数分布5657

第二次迭代：m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;c=0.01440;d=123.1927;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);A=0.2299;B=3.9159;m31=m3*(1-log(m3)+B);n31=m3*A;K1=(m1-m21-m31)/sqrt(n1^2+n21^2+n31^2);R=n1/sqrt(n1^2+n21^2+n31^2);G=n21/sqrt(n1^2+n2^21+n31^2);Q=n31/sqrt(n1^2+n21^2+n31^2);m12=m1+n1*K1*R;m22=m21+G*K1*n21;m32=m31+n31*K1*Q;5757第二次迭代：5758

可靠度验算迭代:whileabs(K1-K0)>(1e-1)*K0m1=m12;n1=0.154*m1;m2=m22;n2=0.07*m2;m3=m32;n3=0.233*m3;K0=K1;c=0.01440;d=123.1927;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);A=0.2299;B=3.9159;m31=m3*(1-log(m3)+B);n31=m3*A;K1=(m1-m21-m31)/sqrt(n1^2+n21^2+n31^2);R=n1/sqrt(n1^2+n21^2+n31^2);G=n21/sqrt(n1^2+n2^21+n31^2);Q=n31/sqrt(n1^2+n21^2+n31^2);m12=m1+n1*K1*R;m22=m21+G*K1*n21;m32=m31+n31*K1*Q;end5858可靠度验算迭代:5859

迭代结果：5959迭代结果：59放映结束！无悔无愧于昨天，丰硕殷实的今天，充满希望的明天。60放映结束！无悔无愧于昨天，丰硕殷实的今天，充满希望的明天。661非正态分布基本变量的情况如果极限状态方程中的基本变量Xi是非正态随机变量，则需首先将非正态变量在一定的条件下等效为正态变量，即进行当量（或等效）正态化。611非正态分布基本变量的情况如果极限状态方程中的基本变量Xi是62

当量正态化条件：在设计验算点P*处非正态变量和当量正态变量的概率分布函数取值相等。（尾部面积相等）在设计验算点P*处非正态变量和当量正态变量的概率密度函数取值相等。（纵坐标相等）622当量正态化条件：263如果随机变量Xi为极值I型分布变量：等效转换后的当量正态随机变量Xi的平均值和标准差分别为Xi′和Xi′，其概率分布函数和概率密度函数分别为Fxi′(x)和fXi′(x)

。由条件①，633如果随机变量Xi为极值I型分布变量：等效转换后的当量正态64从而求得当量正态分布的平均值Xi′为由条件②，644从而求得当量正态分布的平均值Xi′为由条件②，465655566对于非正态随机变量，以从公式(1),(2)求得的Xi′和Xi′分别代替Xi和Xi后，所有的随机变量现在都变成了正态分布随机变量，所以前述正态分布基本变量情况下求和设计验算点P*的公式和方法也就均可应用了。注意：

当X*中仅有部分基本变量为非正态分布时，只需将这部分基本变量当量正态化。如果随机变量Xi为对数正态分布基本变量：将对数正态分布的Xi直接根据当量化处理的两个条件转化为当量正态分布。666对于非正态随机变量，以从公式(1),(2)求得的6767776868886969997070101071现在以从公式(3),(4)求得的Xi′和Xi′分别代替Xi和Xi后，即可将对数正态随机变量变成了正态分布随机变量，接着可按前述公式和方法求和设计验算点P*。根据以上的讨论，对于结构极限状态函数中包含多个正态或非正态基本变量的一般情况，只要知道了各基本变量的概率分布类型及统计参数，就可采用迭代法计算和设计验算点P*的坐标值。其计算框图如下：7111现在以从公式(3),(4)求得的Xi′和Xi72已知:Xi(i=1,……n)的分布类型及统计参数xi,σxi，极限状态方程g(x1……xn)=0假定设计验算点P*的坐标值初值：Xi*(可取Xi*＝xi)对于非正态变量Xi，根据Xi*和公式(1),(2)求出xi′、σxi′以代替xi、σxiA7212已知:Xi(i=1,……n)的分布类型及统计参数xi,73将设计验算点P*的坐标值Xi*代入极限状态方程以求出|上次求出的－上次求出的|≤允许误差以本次求得的Xi*作为下次的取用值本次求得的和Xi*即为所求的可靠指标和设计验算点P*的坐标值是否A7313将设计验算点P*的坐标值Xi*代入极限状态方程74由于是以Z的一阶原点矩和二阶中心矩表达的，且在计算时考虑了基本变量的分布类型，并采用了线性化的近似手段，因此这种结构可靠度的计算方法通常称为“考虑变量分布类型的一次二阶矩方法”。上面介绍的验算点方法是国际安全度联合委员会(JCSS)推荐采用的拉克维茨－菲斯勒法(Rackwitz－Fiessler)，所以简称JC法或R-F法。实际上不同的研究者提出了很多种验算点法，它们各有优缺点，其中我国大连理工大学的赵国藩院士也提出了一种验算点法，计算较JC法简单，计算精度也很高。7414由于是以Z的一阶原点矩和二阶中心矩表达的，且在计算时75RF法算例1：极限状态方程为：EI-78.12P=0цE=2*10^7;бE=0.5*10^7;цI=10^-4;бI=0.2*10^-4;цP=4;бP=1;α=1.2825k=3.5499;7515RF法算例1：极限状态方程为：EI-78.12P=0761、EI正态分布、PI型分布m1=цE=2*10^7;n1=бE=0.5*10^7;m2=цI=10^-4;n2=бI=0.2*10^-4;m3=цP=4;n3=бP=1;c=α=1.2825;d=k=3.5499;Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);76161、EI正态分布、PI型分布m1=цE=2*10^7;n77G=m1*m2-78.12*m3;K0=G/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);7717G=m1*m2-78.12*m3;K0=G/sqrt((78R=-m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);S=m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);T=n31*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);7818R=-m2*n1/sqrt((m2*n1)^2+(m1*79m12=m1+K0*n1*R;m22=m2+K0*S*n2;m32=m31+K0*T*n31;VE=(m12-m1)/n1;VI=(m22-m2)/n2;VP=(m32-m31)/n31;第一次迭代结束7919m12=m1+K0*n1*R;m22=m2+K0*S*n80第二次迭代首先将P的I型变量转为正态变量m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);然后求KI（可靠度）G=m1*m2-78.12*m3;Q=G-m2*n1*VE-m1*n2*VI+78.12*n31*VP;K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);8020第二次迭代首先将P的I型变量转为正态变量2081最后计算新设计验算点的均值R=-m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);S=-m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);m12=2*10^7+K1*0.5*10^7*R;m22=10^-4+K1*S*0.2*10^-4;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m2)/n2;VP=(m32-m31)/n31;第1、2次迭代结果为8121最后计算新设计验算点的均值R=-m2*n1/sqrt((82最后根据可靠度验算条件编程whileabs(K1-K0)>0.003*K0m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;K0=K1;Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m2*n1*VE-m1*n2*VI+78.12*n31*VP;K1=Q/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);R=-m2*n1/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);S=-m1*n2/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);m12=2*10^7+K1*0.5*10^7*R;m22=10^-4+K1*S*0.2*10^-4;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m2)/n2;VP=(m32-m31)/n31;end8222最后根据可靠度验算条件编程2283经过循环迭代结果为8323经过循环迭代结果为23842、E正态、PII型分布这种组合先把PII型分布转为正太变量，用计算后的平均值、标准差求可靠度，因为内容太多就只写程序码和运算结果，不详细写了。m1=2*10^7;n1=0.5*10^7;m2=10^-4;n2=0.2*10^-4;m3=4;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;K0=G/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);R=-m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);S=-m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);T=n31*78.12/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);84242、E正态、PII型分布2485第二次迭加m12=m1+K0*n1*R;m22=m21+K0*S*n21;m32=m31+K0*T*n31;VE=(m12-m1)/n1;VI=(m22-m21)/n21;VP=(m32-m31)/n31;m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m21*n1*VE-m1*n21*VI+78.12*n31*VP;K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);R=-m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);S=-m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);m12=2*10^7+K1*0.5*10^7*R;m22=m21+K1*S*n21;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m21)/n21;VP=(m32-m31)/n31;8525第二次迭加2586循环迭加whileabs(K1-K0)>0.003*K0m1=m12;n1=0.5*10^7;m2=m22;n2=0.2*10^-4;m3=m32;n3=1;c=1.2825;d=3.5499;K0=K1;Fm=exp(-exp(-c*(m2-d)));fm=c*(exp(-c*(m2-d))*exp(-exp(-c*(m2-d))));n21=normpdf(norminv(Fm,0,1),0,1)/fm;m21=m2-n21*norminv(Fm,0,1);Fm=exp(-exp(-c*(m3-d)));fm=c*(exp(-c*(m3-d))*exp(-exp(-c*(m3-d))));n31=normpdf(norminv(Fm,0,1),0,1)/fm;m31=m3-n31*norminv(Fm,0,1);G=m1*m2-78.12*m3;Q=G-m21*n1*VE-m1*n21*VI+78.12*n31*VP;K1=Q/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);R=-m2*n1/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);S=-m1*n21/sqrt((m21*n1)^2+(m1*n21)^2+(n31*78.12)^2);T=n3*78.12/sqrt((m21*n1)^2+(m1