工程结构可靠度中非正态分布转为正态分布课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《工程结构可靠度中非正态分布转为正态分布课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 工程 结构 可靠 中非 正态分布 转为 课件
- 资源描述:
-
1、1非正态分布基本变量的情况如果极限状态方程中的基本变量Xi是非正态随机变量,则需首先将非正态变量在一定的条件下等效为正态变量,即进行当量(或等效)正态化。12 当量正态化条件:在设计验算点P*处非正态变量和当量正态变量的概率分布函数取值相等。(尾部面积相等)在设计验算点P*处非正态变量和当量正态变量的概率密度函数取值相等。(纵坐标相等)23如果随机变量Xi为极值I型分布变量:等效转换后的当量正态随机变量Xi的平均值和标准差分别为Xi和Xi,其概率分布函数和概率密度函数分别为Fxi(x)和fXi(x)。由条件,*1*iiiiiXFXXXFXFXiXiXiXiXiXiXi得34从而求得当量正态分布
2、的平均值Xi为*1*1iiiXiXiXXFX由条件,2exp2112exp212*22*iXiiXiXiiXiiXiXiiXiXXXXfXf45 *1*1*1*1111 2iXXiiiXiiiXXXiiiXiiiXiXXiiiXXiXiiXiiXiXXFXXfXFXFXfX又由公式的推导过程知所以有即得56 对于非正态随机变量,以从公式(1),(2)求得的Xi和Xi分别代替Xi和Xi后,所有的随机变量现在都变成了正态分布随机变量,所以前述正态分布基本变量情况下求和设计验算点P*的公式和方法也就均可应用了。注意:当X*中仅有部分基本变量为非正态分布时,只需将这部分基本变量当量正态化。如果随机变量
3、Xi为对数正态分布基本变量:将对数正态分布的将对数正态分布的Xi 直接根据当量化处理直接根据当量化处理的两个条件转化为当量正态分布。的两个条件转化为当量正态分布。67 iiiXiXiiiXiXiiXXyxXXxXXiysdyeyxXPxFxexxfXlnln2ln0ln2lnlnln 21 0 21 2ln2ln2ln2ln令:有服从对数正态分布,则因为78 分布。态分布来计算对数正态从正因此可以利用式正态分布的分布函数,分布的分布函数化为了此式实际上将对数正态代入上式得则axdsexFdsydyiiiXiXiiXXxsXXa ln 21 lnlnln2lnlnln289 22*lnln22l
4、nln1*ln1*ln*ln*lnlnln22*ln*2 ln ln 11 22 XiiiiiiiiiXiXiiXXiiiXiiXiiXXiXiXXiXXXXiiaFXfXXfXXfXeeXX由公式和 ln 3iX910 *1*ln*1ln*ln*ln*ln*lnln*ln1 ln ln ln 1 ln iiiiiiiiiiiiiXiXiXXiXXXiXXXiiXXiiXaXFXXXXXXXXXX由公式、和上式 41011 现在以从公式(3),(4)求得的Xi和Xi分别代替Xi和Xi后,即可将对数正态随机变量变成了正态分布随机变量,接着可按前述公式和方法求和设计验算点P*。根据以上的讨论,对于
5、结构极限状态函数中包含多个正态或非正态基本变量的一般情况,只要知道了各基本变量的概率分布类型及统计参数,就可采用迭代法计算和设计验算点P*的坐标值。其计算框图如下:1112已知:Xi(i=1,n)的分布类型及统计参数xi,xi,极限状态方程g(x1xn)=0假定设计验算点P*的坐标值初值:Xi*(可取Xi*xi)对于非正态变量Xi,根据Xi*和公式(1),(2)求出xi、xi以代替xi、xi2121*|cosniXPiXPiXiiiXgXg求A1213的坐标值即求出验算点代入坐标转换公式:将求得的*cos,.cosPXi xxixiixixiix将设计验算点P*的坐标值Xi*代入极限状态方程
6、以求出0,*2*1nXXXg|上次求出的上次求出的|允许误差以本次求得的Xi*作为下次的取用值 本次求得的和Xi*即为所求的可靠指标和设计验算点P*的坐标值是否A1314由于是以Z的一阶原点矩和二阶中心矩表达的,且在计算时考虑了基本变量的分布类型,并采用了线性化的近似手段,因此这种结构可靠度的计算方法通常称为“考虑变量分布类型的一次二阶矩方法”。上面介绍的验算点方法是国际安全度联合委员会(JCSS)推荐采用的拉克维茨菲斯勒法(RackwitzFiessler),所以简称JC法或R-F法。实际上不同的研究者提出了很多种验算点法,它们各有优缺点,其中我国大连理工大学的赵国藩院士也提出了一种验算点法
7、,计算较JC法简单,计算精度也很高。1415极限状态方程为:EI-78.12P=0 E=2*107;E=0.5*107;I=10-4;I=0.2*10-4;P=4;P=1;=1.2825 k=3.5499;15161、EI正态分布、PI型分布m1=E=2*107;n1=E=0.5*107;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;
8、m31=m3-n31*norminv(Fm,0,1);1617G=m1*m2-78.12*m3;K0=G/sqrt(m2*n1)2+(m1*n2)2+(n31*78.12)2);1718R=-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);1819m12=m1+K0*n1*R;m22=m2+K0*S*n2;m32=m31+K0*T*n31;VE=(m12-m1)/n1;VI=(
9、m22-m2)/n2;VP=(m32-m31)/n31;第一次迭代结束1920第二次迭代首先将P的I型变量转为正态变量m1=m12;n1=0.5*107;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+7
10、8.12*n31*VP;K1=Q/sqrt(m2*n1)2+(m1*n2)2+(n31*78.12)2);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*107+K1*0.5*107*R;m22=10-4+K1*S*0.2*10-4;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m2)/n2;
11、VP=(m32-m31)/n31;第1、2次迭代结果为2122最后根据可靠度验算条件编程while abs(K1-K0)0.003*K0m1=m12;n1=0.5*107;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
12、*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*107+K1*0.5*107*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)
13、/n31;end2223经过循环迭代结果为23242、E正态、PI I型分布这种组合先把PII型分布转为正太变量,用计算后的平均值、标准差求可靠度,因为内容太多 就只写程序码和运算结果,不详细写了。m1=2*107;n1=0.5*107;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(-e
14、xp(-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*
15、n21)2+(n31*78.12)2);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*107;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-n2
16、1*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
17、(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*107+K1*0.5*107*R;m22=m21+K1*S*n21;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m21)/n21;VP=(m32-m31)/n31;2526循环迭加while abs(K1-K0)0.003*K0m1=m12;n1=0.5*107;m2=m22;n2=0.2*10-4;m3=m32;n3=1;c=1.2825;d=3.5499;K0=K1;Fm=exp
18、(-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
19、/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*107+K1*0.5*107*R;m22=m21+K1*S*n21;m32=m31+K1*T*n31;VE=(m12-m1)/n1;VI=(m22-m21)/n21;VP=(m32-m31)/n31;end2627最终运行
20、结果为:27283、EIP都为I型分布这种组合应该先把EIP都转化为正态分布,在利用均值、标准差求可靠度 第一次迭加:m1=2*107;n1=0.5*107;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(-ex
21、p(-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+(n3
22、1*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;2829第二次迭加:m1=m12;n1=0.5*107;m2=m22;n2=0.2*10-4;m3=m32;n3=1;c=1.2825;d=3.5499;Fm=exp(-e
23、xp(-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);
24、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+
25、(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;2930可靠度精度要求循环迭加:while abs(K1-K0)0.003*K0m1=m12;n1=0.5*107;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=normpd
展开阅读全文