结构可靠度编程作业.docVIP

下载本文档

466
0
约6.75千字
约 8页
2017-07-22 发布于重庆
举报
版权申诉

结构可靠度编程作业.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

结构可靠度编程作业

Matlab编程作业姓名：余振学号专业：结构工程指导老师：杨春侠习题一：承受恒载和楼面活荷载的钢筋混凝土轴心受压短柱，已知恒载产生的轴向力NG为正态分布，活载产生的轴向力NL为极值Ⅰ型分布，截面承载能力（抗力）R为对数正态分布，统计参数分别为(NG＝1159.1kN， (NG＝81.1kN， (NL＝765.5kN， (NL＝222kN，(R＝4560kN，(R＝729.6kN，极限状态方程为 Z＝R－NG－NL＝0，求可靠指标(和设计验算点。方法一 muX=[4560,765.5,1159.1]; %均值 cvX=[729.6/4560,222/765.5,81.1/1159.1]; %变异系数 bbeta=JC_3(muX,cvX) M-file ： function bbeta=JC_3(muX,cvX) sigmaX=cvX.*muX; %标准差 sLn=sqrt(log(1+(sigmaX(1)/muX(1))^2)); mLn=log(muX(1)/sqrt(1+cvX(1)^2)); aEv=pi/(sqrt(6)*sigmaX(2)); %求极值I型分布的参数 psi=0.5772 uEv=muX(2)-psi/aEv; sigmaX1=sigmaX; x=muX; m=0; %记录循环次数 b0=0; %可靠度的初始值 b1=3.2; %可靠度的第二个值 while abs((b1-b0)/b1)1e-3 sigmaX1(1)=sLn*muX(1); %抗力R当量正态化 muX1(1)=muX(1)*(1+mLn-log(muX(1))); t=exp((-aEv)*(muX(2)-uEv)); %活载L当量正态化 f1=aEv*exp((-aEv)*(muX(2)-uEv)-t); f2=exp(-t); a=norminv(f2); %标准正态概率密度反函数求值 y=normpdf(a); %标准正态概率密度求值 sigmaX1(2)=y/f1; muX1(2)=muX(2)-a*sigmaX1(2); sigmaX1=[sigmaX1(1);sigmaX1(2);sigmaX1(3)]; w=norm(sigmaX1); bbeta=(muX1(1)-muX1(2)-muX1(3))/w; %求bbeta值 Alphar=-sigmaX1(1)/w; %求方向余弦 Alphal=sigmaX1(2)/w; Alphag=sigmaX1(3)/w; x(1)=muX1(1)+Alphar*bbeta*sigmaX1(1); %求循环后验算点的坐标值 x(2)=muX1(2)+Alphal*bbeta*sigmaX1(2); x(3)=muX1(3)+Alphag*bbeta*sigmaX1(3); x=[x(1);x(2);x(3)]; muX=x; m=m+1; b0=b1; b1=bbeta; end disp(结果：) fprintf(循环次数m：m=%d\n,m); fprintf(可靠度指标贝塔：bbeta=%1.2f\n,bbeta); fprintf(最后验算点坐标：muX=[%1.2f;%1.2f;%1.2f]\n,muX); 结果：循环次数m：m=6 可靠度指标贝塔：bbeta=3.959 最后验算点坐标：muX=[3011.84;1817.98;1193.86] 方法二 muX=[4560,765.5,1159.1]; cvX=[729.6/4560,222/765.5,81.1/1159.1]; bbeta=JC_3(muX,cvX) M-file如下： function bbeta=JC_3(muX,cvX) sigmaX=cvX.*muX; sLn=sqrt(log(1+(sigmaX(1