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

非正态分布基本变量的情况如果极限状态方程中的基本变量Xi是非正态随机变量，则需首先将非正态变量在一定的条件下等效为正态变量，即进行当量（或等效）正态化。2021/7/5 当量正态化条件：在设计验算点P*处非正态变量和当量正态变量的概率分布函数取值相等。（尾部面积相等）在设计验算点P*处非正态变量和当量正态变量的概率密度函数取值相等。（纵坐标相等）2021/7/5如果随机变量Xi为极值I型分布变量 ：等效转换后的当量正态随机变量Xi的平均值和标准差分别为?Xi′和?Xi′，其概率分布函数和概率密度函数分别为Fxi′(x)和fXi′(x) 。由条件①，2021/7/5从而求得当量正态分布的平均值?Xi′为由条件②，2021/7/52021/7/5 对于非正态随机变量，以从公式(1),(2) 求得的?Xi′和?Xi′分别代替?Xi和?Xi后，所有的随机变量现在都变成了正态分布随机变量，所以前述正态分布基本变量情况下求?和设计验算点P*的公式和方法也就均可应用了。注意： 当X*中仅有部分基本变量为非正态分布时，只需将这部分基本变量当量正态化。 如果随机变量Xi为对数正态分布基本变量：将对数正态分布的Xi 直接根据当量化处理的两个条件转化为当量正态分布。2021/7/52021/7/52021/7/52021/7/52021/7/5 现在以从公式(3),(4) 求得的?Xi′和?Xi′分别代替?Xi和?Xi后，即可将对数正态随机变量变成了正态分布随机变量，接着可按前述公式和方法求?和设计验算点P*。 根据以上的讨论，对于结构极限状态函数中包含多个正态或非正态基本变量的一般情况，只要知道了各基本变量的概率分布类型及统计参数，就可采用迭代法计算?和设计验算点P*的坐标值。其计算框图如下：2021/7/5已知:Xi(i=1,……n)的分布类型及统计参数?xi,σxi，极限状态方程g(x1……xn)=0假定设计验算点P*的坐标值初值：Xi*(可取Xi*＝ ?xi) 对于非正态变量Xi，根据Xi*和公式(1),(2)求出?xi′、σxi′以代替?xi、σxiA2021/7/5A将设计验算点P*的坐标值Xi*代入极限状态方程 以求出?以本次求得的Xi*作为下次的取用值 否|上次求出的?－上次求出的?| ≤允许误差?是本次求得的?和Xi*即为所求的可靠指标和设计验算点P*的坐标值2021/7/5由于?是以Z的一阶原点矩和二阶中心矩表达的，且在计算?时考虑了基本变量的分布类型，并采用了线性化的近似手段，因此这种结构可靠度的计算方法通常称为“考虑变量分布类型的一次二阶矩方法”。上面介绍的验算点方法是国际安全度联合委员会(JCSS)推荐采用的拉克维茨－菲斯勒法(Rackwitz－Fiessler)，所以简称JC法或R-F法。实际上不同的研究者提出了很多种验算点法，它们各有优缺点，其中我国大连理工大学的赵国藩院士也提出了一种验算点法，计算较JC法简单，计算精度也很高。2021/7/5RF法算例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.2825 k=3.5499;2021/7/51、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);2021/7/5G=m1*m2-78.12*m3;K0=G/sqrt((m2*n1)^2+(m1*n2)^2+(n31*78.12)^2);2021/7/5R=-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);2021/7/5m12=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;第一次迭代结束2021/7/5第二次迭代首先将P的I型变量转为正态变量m1=m12;n1=0.5*10^