工程结构可靠度中非正态分布转为正态分布.pptVIP

非正态分布基本变量的情况; 当量正态化条件： 在设计验算点P*处非正态变量和当量正态变量的概率分布函数取值相等。（尾部面积相等） 在设计验算点P*处非正态变量和当量正态变量的概率密度函数取值相等。（纵坐标相等）;如果随机变量Xi为极值I型分布变量 ： ;从而求得当量正态分布的平均值?Xi′为;; 对于非正态随机变量，以从公式(1),(2) 求得的?Xi′和?Xi′分别代替?Xi和?Xi后，所有的随机变量现在都变成了正态分布随机变量，所以前述正态分布基本变量情况下求?和设计验算点P*的公式和方法也就均可应用了。;;;;; 现在以从公式(3),(4) 求得的?Xi′和?Xi′分别代替?Xi和?Xi后，即可将对数正态随机变量变成了正态分布随机变量，接着可按前述公式和方法求?和设计验算点P*。;已知:Xi(i=1,……n)的分布类型及统计参数?xi,σxi，极限状态方程g(x1……xn)=0;将设计验算点P*??坐标值Xi*代入极限状态方程 以求出?;由于?是以Z的一阶原点矩和二阶中心矩表达的，且在计算?时考虑了基本变量的分布类型，并采用了线性化的近似手段，因此这种结构可靠度的计算方法通常称为“考虑变量分布类型的一次二阶矩方法”。;RF法算例1：;1、EI正态分布、PI型分布;;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=n31*78.12/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);;m12=m1+K0*n1*R;m22=m2+K0*S*n2;m32=m31+K0*T*n31;;第二次迭代;最后计算新设计验算点的均值;;经过循环迭代结果为 ;2、E正态、PI I型分布 这种组合先把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);;第二次迭加 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*n