结构可靠度编程作业

承受恒载和楼面活荷载的钢筋混凝土轴心受压短柱，已知恒载产方法一muX=[4560,765.5,1159.1];%均值cvX=[729.6/4560,222/765.5,81.1/1159.1];%变异系数mLn=log(muX(1)/sqrt(1+cvaEv=pi/(sqrt(6)*sigmaX(2));%求极值I型分布的参数psi=0.5772while%记录循环次数%可靠度的初始值%可靠度的第二个值%抗力R当量正态化muX1(1)=muX(1)*(1+mLn-log(mt=exp((-aEv)*(muX(2%活载L当量正态化%标准正态概率密度反函数求值%标准正态概率密度求值muX1(2)=muX(2)-a*sigAlphar=-sigmaX1(1)/w;%求方向余弦x(1)=muX1(1)+Alphar*bbeta*sigmaX1(1):%求循环后验算点的坐标值x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbeta*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=6cvX=[729.6/4560,222/7M-file如下：whileabs(norm(x)-normX)/normX>1e-6%记录循环次数%向量x二范数的前后两个值的相对误差%计算x的二范数%抗力R的当量正态化muX1(1)=muX(1)*(1+mLn-log(mut=exp((-1/aEv)*(muX(2)-uEv));%活载L的当量正态化muX1(2)=muX(2)-a*sigx(1)=muX1(1)+Alphar*bbeta*x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbeta*fprintf(可靠度指标贝塔：bbeta=%1.2f\n',bbeta);fprintf(最后验算点坐标：muX=[%1.2f;%1.2f;%1.2f]\n',muX);循环次数m:m=7cvX=[729.6/4560,222/7EQ\*jc3\*hps46\o\al(\s\up4(b),M)-fEQ\*jc3\*hps46\o\al(\s\up4(et),ile)mLn=log(muX(1)/sqrt(1+cv%计算标准差以便下面使用%求极值I型分布的参数muX1(1)=muX(1)*(1+mLn-log(mu%对抗力R进行当量正态化%对活载L进行当量正态化EQ\*jc3\*hps36\o\al(\s\up15(w=nor),bbeta)EQ\*jc3\*hps36\o\al(\s\up15(sigma),muX)x(1)=muX1(1)+Alphar*bbeta*sigmaX1(1);%求循环后验算点的坐标值x(2)=muX1(2)+Alphal*bbetaEQ\*jc3\*hps36\o\al(\s\up13(x),x)EQ\*jc3\*hps36\o\al(\s\up13(+),2)EQ\*jc3\*hps36\o\al(\s\up13(l),;)EQ\*jc3\*hps36\o\al(\s\up13(ph),x)EQ\*jc3\*hps36\o\al(\s\up13(b),])muX=x;%将x赋予给muX以便下次循环使用fprintf(可靠度指标贝塔：bbeta=%1.2f\n',bbeta);fprintf(最后验算点坐标：muX=[%1.2f;%1.2f;%1.2f]\n',muX);方法一.cvX=[0.17,0.29,0.07];%变异系数muX1(2)=muX(2)-a*sig%计算x的二范数%活载L的当量正态化x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbetafprintf(循环次数m:m=%d\n',m);fprintf(最后验算点坐标：muX=[%1.2f;%1.2f;%1.2f]\n’,muX);循环次数m:m=9方法二cvX=[0.17,0.29,0.07];%变异系数%求极值I型分布的参数%对活载L进行当量正态化muX1(2)=muX(2)-a*sigsigmaX1=[sigmaX1(1)Alphar=-sigmaX1(1)/w;%求方向余弦x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbeta*x(1)=x(2)+x(3);%求循环后验算点的坐标值x=[x(1);x(2);x(3)];%将计算的验算点坐标放入数组xmuX=x;