montecarlo计算结构可靠度的matlab程序-matlab work.rar

%问题2、按原混凝土设计规范设计的轴心受压构件极限状态方程为：R－SG－SQ = 0，ρ= 2.0； %求β。 %其中：KR = 1.33，δR = 0.17； %KSG = 1.06，δSG = 0.07； %KSQ = 0.70，δSQ = 0.29。 function [bb]=kkdmc mm=1000; nn=mm^2; pp=0.2;%pp=sq/sg荷载效应比值 sg=10; sq=sg*pp; r=1.55*(sg+sq);%按一个分项系数设计 u=[1.33 1.06 0.7];%矩阵U第一列代表抗力R,第二列代表SG,第三列代表SQ p=[0.17 0.07 0.29];%矩阵p的三列分别为R,SG,SQ的标准差 a1=log(lognrnd(u(1)*r,p(1)*r,mm,mm));%对数正态分布随机抽样 a3=normrnd(u(2)*sg,p(2)*sg,mm,mm);%正态分布随机抽样 aaa=pi/(p(3)*sq*sqrt(6));%m为极值I型当量正态化过程中的一个常数 bbb=u(3)*sq-0.5772/aaa;%n为极值I型当量正态化过程中的另一个常数 a=rand(mm); a4=(aaa*bbb-log(log(1./a)))/aaa; l=a1-a3-a4;%功能函数 o=find(l<0);%失效的次数的位置 s=length(o);%失效的次数 bb=norminv(1-s/mm^2)%可靠度指标 end %结果3.8