结构可靠度设计原理与应用.docx
《结构可靠度设计原理与应用.docx》由会员分享，可在线阅读，更多相关《结构可靠度设计原理与应用.docx（6页珍藏版）》请在冰豆网上搜索。
结构可靠度设计原理与应用
结构可靠度设计原理与应用
结构可靠度设计原理与应用华中科技大学研究生课程考试答题本考生姓名考生学号系、年级结构工程硕1401班类别学术型考试科目结构可靠度设计原理与应用考试日期2015年1月20日评分题号得分题号得分总分：
评卷人：
注：
1、无评卷人签名试卷无效。
2、必须用钢笔或圆珠笔阅卷，使用红色。
用铅笔阅卷无效。
题号回答内容得分中心点法1.如图所示圆截面直杆，承受拉力，已知材料的强度设计值的均值，标准差，杆直径的均值，标准差，在功能函数为：
1）；2），在这两种情况下，试用中心点法求其可靠度指标和可靠度。
（5分）解：
（1）clearall;clc;mufy=310;sigmafy=25;mud=30;sigmad=3;P=120000;symsfyd;%定义符号变量fy和dZ=（pi*d/4）*fy-P;%定义目标函数pdfy=diff（Z,fy,1）;%Z对fy求一阶偏导pdd=diff（Z,d,1）;%Z对gmafy=25;mud=30;sigmad=3;P=120000;symsfyd;%定义符号变量fy和dZ=fy-4*P/（pi*d）;%定义功能函数pdfy=diff（Z,fy,1）;%Z对fy求一阶偏导pdd=diff（Z,d,1）;%Z对d求一阶偏导fy=mufy;d=mud;%将均值赋给fy和da=subs（pdfy）;%求在均值点处的偏导数b=subs（pdd）;c=subs（Z）;%求功能函数在均值点的值muZ=c;sigmall;clc;muw=100;sigmaw=20;mufai=0.6109;sigmafai=0.0873;tao=52;%角度输入要用弧度制，注意单位要统一symswfai;%定义变量类型为符号变量Z=w*tan（fai）-tao;pdw=diff（Z,w,1）;pdfai=diff（Z,fai,1）;w=muw;fai=mufai;a=subs（pdw）;b=subs（pdfai）;c=subs（Z）;muZ=c;sigmaZ=（a*sigma分布；恒载,服从正太分布；活载Q服从极致I型分布，。
试用JC法求目标可靠指标时，构件截面的抗力平均值（20分）（提示：
;）解：
JC法clearall;clc;muX=[319.52;53;70];cuX=[0.17;0.07;0.29];sigmaX=cuX.*muX;sLn=sqrt（log（1+（sigmaX
（1）/muX
（1））））;mLn=log（muX
（1））-sLn/2;aEv=sqrt（6）*sigmaX（3）/pi;uEv=0.5772*aEv-muX（3）;muX1=muX;sigmaX1=sigmaX;x=muX;normX=eps;count=0;whileabs（norm（x）-normX）/normX>1e-6normX=norm（x）;g=x
（1）-x
（2）-x（3）;gX=[1;-1;-1];cdfX=[logncdf（x
（1）,mLn,sLn）;1-evcdf（-x（3）,uEv,aEv）];pdfX=[lognpdf（x
（1）,mLn,sLn）;evpdf（-x（3）,uEv,aEv）];nc=norminv（cdfX）;sigmaX1（1:
2:
3）=normpdf（nc）./pdfX;muX1（1:
2:
3）=[x（1:
2:
3）-nc.*sigmaX1（1:
2:
3）];gs=gX.*sigmaX1;ccosX=-gs/norm（gs）;count=count+1;x=muX1+3.7*sigmaX1.*ccosXendxR=x
（1）cosR=ccosX
（1）;sigmaX1=x.*sqrt（log（1+cuX））;sigmaRs=sigmaX1
（1）;muX1=x-3.7*sigmaX1.*ccosX;muRs=muX1
（1）;muR=sqrt（1+cuX
（1））.*exp（log（x
（1））-1+muX1
（1）./x
（1））displayString1=[抗力平均值,num2str（muR）]displayString2=[迭代次数为,num2str（count）]%结论：
JC法求解，当目标可靠指标[β]=3.7时，经过8次迭代，抗力平均值为519.7522。
蒙特卡罗法5.设某构件正截面强度计算的极限状态方程为Z=R-S=0。
其中R和S分别为正态和极值I型分布的随机变量，其统计量为R（100,20）和S（80,24），20和24为标准差。
试用JC法和蒙特卡罗模拟分别求解构件失效概率。
（20分）解：
（1）JC法clearall;clc;symsRS;muR=100;sigmaR=20;muS=80;sigmaS=24;Z=R-S;alpha=0.78*sigmaS;k=muS-0.5772*alpha;%alpha和k为极值I型分布的两个参数symsx;Scdf=exp（-exp（-（x-k）/alpha））;%S概率分布函数Spdf=exp（-（x-k）/alpha-exp（-（x-k）/alpha））/alpha;%S概率密度函数r=muR;s=muS;%初始验算点x=s;a=subs（Scdf）;b=subs（Spdf）;%求1pf=normcdf（-beta）%pf=0.2221
（2）蒙特卡罗法clearall;clc;symsRS;muR=100;sigmaR=20;muS=80;sigmaS=24;Z=R-S;r=normrnd（muR,sigmaR,1000,1000）;s=rand（1000,1000）;alpha=0.78*sigmaS;k=muS-0.5772*alpha;%极值I型分布的两个参数symsx;y=exp（-exp（-（x-k）/alpha））;%y为S的概率分布函数f=finverse（y）;x=s;s1=subs（f）;R=r;S=s1;z=subs（Z）;m=0;n=0;fori=1:
1000forj=1:
1000ifz（i,j）。
式中，，服从对数正态分布；，为正态分布；，为正态分布；，为正态分布。
试用蒙特卡洛法计算该结构构件的可靠度。
（10分）解：
clearall;clc;symsx1x2x3x4x111;%x111=log（x1）mux1）;x2=normrnd（mux2,sigmax2,1000,1000）;x3=normrnd（mux3,sigmax3,1000,1000）;x4=normrnd（mux4,sigmax4,1000,1000）;z=subs（Z）;m=0;n=0;fori=1:
1000forj=1:
1000ifz（i,j）试用蒙特卡洛法计算该结构构件的可靠度。
（10分）解：
clearall;clc;symsx1x2x3x4x111x222;3000,3000）;x2=lognrnd（mux222,sigmax222,3000,3000）;x3=normrnd（mux3,sigmax3,3000,3000）;x=rand（3000,3000）;x4=subs（f）;z=subs（Z）;m=0;n=0;fori=1:
3000forj=1:
3000ifz（i,j）基于随机域的卡亨南-洛维变换展开理论，本文利用埃尔米特多项式混沌展开理论来建立随机响应面函数，其未知量可以通过概率搭配方法计算得到，结构的可靠度可以通过几何方法求得。
此方法和蒙特卡罗法比较，可以得出：
在分析二阶和三阶随机响应面的性质时，前者可以得到较高的计算精度，同时计算效率也比较高。
关键词：
随机响应面法；结构可靠度；随机域；蒙特卡罗模拟1引言在结构工程中，可靠度分析已经被越来越广泛的使用。
对于大型而又复杂的工程结构，可靠度代表结构对荷载和环境的承载能力。
因此，如何高效的计算结构的可靠度有很重要的实践意义。
各种各样的可靠度分析方法，如蒙特卡罗模拟法和扰动随机有限元法等等，已经在文化领域得到认可。
一般情况下，第一步就是建立结构的极限状态方程。
响应面法是计算结构可靠度的有效方法中的一种。
法拉韦利【1】根据大量的实验研究建立了确定性响应面法；巴赫尔【2】建立了插值迭代技术，主要用来模拟极限状态面，尽管计算量较小，检查点也是近似的，但仍具有很重要的意义；卡伊马兹【3】提出了一种加权回归法，此方法的关键就是要求所选的试验点在最可能的地方。
随机响应面法（SRSM）是传统响应面法的扩展。
两种方法最显著的区别在于建立极限状态函数的方式不同。
本文所用的随机响应面法被用来模拟输入的不确定性和系统反应之间的复杂关系，并将输出结果通过埃尔米特多项式混沌展开理论转化为标准随机变量。
2谱随机有限元法2.1随机域离散化一般地，随机域不仅是空间坐标轴的函数，而且也是随机变量的函数。
零平均值随机域表示为：
（1）其中和是随机域的均值和协方差函数。
通常情况下，是有界限的、对称的和正定的，同时也不局限于同一领域，因此它可以分解为：
（2）其中特征值和相关函数可以从下面的积分等式的解中得到。
（3）因此随机域可以被近似地离散为随机变量，如下式所示：
（4）其中是一组不相关的高斯随机变量，其均值为0，当其展开项达到第三项时，其计算结果就会比较精确。
此外，随机变量也满足下面这个等式。
（5）2.2随机刚度方程的扩展对于随机结构而言,它的刚度方程可以表示为（6）其中刚度矩阵可以表示为（7）考虑材料属性的随机性，弹性刚度矩阵可以用下式表示：
（8）其中是一个确定的矩阵，表示材料属性的随机域，其均值和正确的结构模型都是已知的。
因此，联立方程（4）和（8）可以得到方程（9），其表达式如下：
（9）将式（4）和式（8）代入式（7），则可得到扩展刚度矩阵，其表达式如下：
（10）其中是单元刚度矩阵的平均值，是随着的变化而有微小变化的单元刚度矩阵。
将上式与随机域的离散化方程联合，可以得到随机结构有限元方程，其表达式如下：
（11）3基于随机响应面法（SRSM）的结构可靠度计算在运用随机响应面法之前，控制点P的随机域可以通过一组不相关的随机变量根据式（4）近似得到，其表达式可表示为：
（12）其中，是未知系数，多维数组P服从埃尔米特多项式，且能用下式表示出来：
（13）因此，式（12）可以被改写为：
（14）其中是待定系数。
和与和不同。
在的样本空间里，每一组都相当于一个配置点。
将每一个配置点代入式（11），可以得到结点位移矩阵。
所以，随机响应面法的关键步骤是如何合理地选择配置点。
为了得到更高的计算精度，配置点应该从下一个高阶埃尔米特多项式的根中进行选择。
与此同时，配置点应该接近于原点，且与原点对称【4】。
一般地，基于回归的方法可以通过选择系数数量两倍多的点来得到更高的计算精度【5】。
将和代入式（14），可以得到一个线性的代数方程，如下式所示：
（15）其中是待定系数，可以通过线性最小平方法得到。
将代入式（12），可以得到指定点P的极限状态方程，如下所示：
（16）其中代表
升级会员 