結(jié)構(gòu)可靠度編程作業(yè)

承受恒載和樓面活荷載的鋼筋混凝土軸心受壓短柱，已知恒載產(chǎn)方法一muX=[4560,765.5,1159.1];%均值cvX=[729.6/4560,222/765.5,81.1/1159.1];%變異系數(shù)mLn=log(muX(1)/sqrt(1+cvaEv=pi/(sqrt(6)*sigmaX(2));%求極值I型分布的參數(shù)psi=0.5772while%記錄循環(huán)次數(shù)%可靠度的初始值%可靠度的第二個(gè)值%抗力R當(dāng)量正態(tài)化muX1(1)=muX(1)*(1+mLn-log(mt=exp((-aEv)*(muX(2%活載L當(dāng)量正態(tài)化%標(biāo)準(zhǔn)正態(tài)概率密度反函數(shù)求值%標(biāo)準(zhǔn)正態(tài)概率密度求值muX1(2)=muX(2)-a*sigAlphar=-sigmaX1(1)/w;%求方向余弦x(1)=muX1(1)+Alphar*bbeta*sigmaX1(1):%求循環(huán)后驗(yàn)算點(diǎn)的坐標(biāo)值x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbeta*disp(結(jié)果：)fprintf(循環(huán)次數(shù)m:m=%d\n',m);fprintf(可靠度指標(biāo)貝塔：bbeta=%1.2f\n',bbeta);fprintf(最后驗(yàn)算點(diǎn)坐標(biāo)：muX=[%1.2f;%1.2f;%1.2f]\n’,muX);循環(huán)次數(shù)m:m=6cvX=[729.6/4560,222/7M-file如下：whileabs(norm(x)-normX)/normX>1e-6%記錄循環(huán)次數(shù)%向量x二范數(shù)的前后兩個(gè)值的相對(duì)誤差%計(jì)算x的二范數(shù)%抗力R的當(dāng)量正態(tài)化muX1(1)=muX(1)*(1+mLn-log(mut=exp((-1/aEv)*(muX(2)-uEv));%活載L的當(dāng)量正態(tài)化muX1(2)=muX(2)-a*sigx(1)=muX1(1)+Alphar*bbeta*x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbeta*fprintf(可靠度指標(biāo)貝塔：bbeta=%1.2f\n',bbeta);fprintf(最后驗(yàn)算點(diǎn)坐標(biāo)：muX=[%1.2f;%1.2f;%1.2f]\n',muX);循環(huán)次數(shù)m:m=7cvX=[729.6/4560,222/7-fmLn=log(muX(1)/sqrt(1+cv%計(jì)算標(biāo)準(zhǔn)差以便下面使用%求極值I型分布的參數(shù)muX1(1)=muX(1)*(1+mLn-log(mu%對(duì)抗力R進(jìn)行當(dāng)量正態(tài)化%對(duì)活載L進(jìn)行當(dāng)量正態(tài)化x(1)=muX1(1)+Alphar*bbeta*sigmaX1(1);%求循環(huán)后驗(yàn)算點(diǎn)的坐標(biāo)值x(2)=muX1(2)+Alphal*bbetamuX=x;%將x賦予給muX以便下次循環(huán)使用fprintf(可靠度指標(biāo)貝塔：bbeta=%1.2f\n',bbeta);fprintf(最后驗(yàn)算點(diǎn)坐標(biāo)：muX=[%1.2f;%1.2f;%1.2f]\n',muX);方法一.cvX=[0.17,0.29,0.07];%變異系數(shù)muX1(2)=muX(2)-a*sig%計(jì)算x的二范數(shù)%活載L的當(dāng)量正態(tài)化x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbetafprintf(循環(huán)次數(shù)m:m=%d\n',m);fprintf(最后驗(yàn)算點(diǎn)坐標(biāo)：muX=[%1.2f;%1.2f;%1.2f]\n’,muX);循環(huán)次數(shù)m:m=9方法二cvX=[0.17,0.29,0.07];%變異系數(shù)%求極值I型分布的參數(shù)%對(duì)活載L進(jìn)行當(dāng)量正態(tài)化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);%求循環(huán)后驗(yàn)算點(diǎn)的坐標(biāo)值x=[x(1);x(2);x(3)];%將計(jì)算的驗(yàn)算點(diǎn)坐標(biāo)放入數(shù)組xmuX=x;