结构可靠性作业.docVIP

华 中 科 技 大 学 研 究 生 课 程 考 试 答 题 本 考生姓名 考生学号 系、年级 类 别 考试科目 考试日期 年 月 日  评 分 题 号 得 分 题 号 得 分 总 分： 评 卷 人： 注：1、无评卷人签名试卷无效。 2、必须用钢笔或圆珠笔阅卷，使用红色。用铅笔阅卷无效。 题号 答 题 内 容 评分 1. 设极限状态方程为，式中，，为对数正态分布；，为对数正态分布；，为正态分布；，为极值I型分布。分别用中心点法、验算点法、蒙特卡洛法计算可靠指标并比较计算结果。 答：中心点法： 又称一次二阶矩法，β= UG/σG; 下面是相关程序及结果： function []=cenp() clear all;clc; mux1=2234.32;sigmax1=0.1;mux2=949.59;sigmax2=0.1; mux3=1521.9;sigmax3=0.109;mux4=496.1;sigmax4=0.292; %% beita=mu(G(x1,x2,x3,x4))/sigma(G(x1,x2,x3,x4)); %% G(x1,x2,x3,x4)=x1+x2-x3-x4; muG=mux1+mux2-mux3-mux4; sigmaG=sqrt(sigmax1^2+sigmax2^2+sigmax3^2+sigmax4^2); beita=muG/sigmaG %% Pf(G)=Φ(-beita); Pr(G)=1-Φ(-beita)=Φ(beita); Pr=normcdf(beita) 算得：β= 3406.5；有效概率Pr=1; 验算点法： 对于上述中心点法，实际上将不是正态分布的随机变量当成正态变量处理，结果并不准确。而且若极限状态方程为非线性，则要进行设计验算点处的泰勒展开，略去二阶以上变量，保留线性项，再做设计点处的正态化。一般工程中选均值处作为设计点。但是，在均值处进行泰勒展开及正态化是我们人为选择的，这种选择算出的β值是否准确并不知道，而且随着设计点的不同，β值并不相同。于是由β值是正则化空间中从原点到极限状态曲面的最短距离这一定义，人们得到一种迭代求解 β值的方法，从而使上述混乱得以解决，这种解法被称为验算点法，下面是相关程序及结果： function yansdf() clc;clear all; syms x1 x2 x3 x4; mux1=2234.32;sigmax1=0.1;mux2=949.59;sigmax2=0.1; mux3=1521.9;sigmax3=0.109;mux4=496.1;sigmax4=0.292; %% double tyx1=2;tyx2=2;tyx3=1;tyx4=3; %% norm:1;logn=2;ev=3; tyx=[tyx1 tyx2 tyx3 tyx4]; %% double g=x1+x2-x3-x4; %% sym x=[x1 x2 x3 x4]; %% sym for i = 1 : 4 dg(i) = diff(g,x(i)); %% sym end p=[mux1,mux2,mux3,mux4]; %% mu=p; mus=mu; %% double;for calculating mud=mu; %% double;for calculating sigma=[sigmax1,sigmax2,sigmax3,sigmax4]; %% double sigmas=sigma; %% double;for calculating sigmad=sigma; for i=1:4 if tyx(i)==1 elseif tyx(i)==2 s=log(1+sigma(i)/mu(i)); sigma(i)=sqrt(s); mu(i)=log(mu(i))-0.5*s; elseif tyx(i)==3 sigma(i)=sqrt(6)*sigma(i)/pi; mu(i)=mu(i)+psi(1)*sigma(i); % -psi(1) is eule