结构系统可靠性课程设计.docVIP

核科学与技术学院 结构系统可靠性分析课程设计 学 号： 专 业：核科学与技术 学生姓名： 任课教师： 教授 2015年4月 问题的提出 理想情况下，若元件的功能函数为，各自变量的联合概率密度函数为，则元件失效概率可以通过解析的方法求得，即： 然而，由于数学，上式不仅无法通过解析法求解，利用数值积分求解时面临计算效率的问题。 对功能函数进行标准化变换后，改进一次二阶矩法消除了功能函数的形式对可靠性指标的计算的影响。一般认为，标准化后，从原点到的最短距离即为可靠性指标β，因此改进一次二阶矩法的主要问题是寻找这一距离，或者说寻找上距离原点最近的点。 我对课堂上所学到的方法的理解是，若在极限状态方程M=0 基本假设 数学模型 when 构造拉格朗日函数，有： 将上述函数对zi与求λ偏导后令其等于0，得到非线性方程组： 上述方程组可以写为： 即F(x)=0。 向量函数F(x)的雅可比矩阵为： 则根据牛顿迭代法，有迭代式： 算例与程序实现 σ，直径为D，所受拉力为T，则拉杆的功能函数可以写为g=σ-(4T/πD2)；假设各变量的均值与标准差如表4.1所示。 表4.1 算例输入数据表 参数 均值 标准差 σ 170MPa 10.8MPa D 0.0294m 0.0005m T 40kN 4.5kN 所编制的MATLAB程序如下： syms g z1 z2 z3 lamd; g(z1,z2,z3)=170+10.8*z1-4*(0.04+0.0045*z2)/(pi*(0.0294+0.0005*z3)^2); L=z1^2+z2^2+z3^2+lamd*g; f1=diff(L,z1); f2=diff(L,z2); f3=diff(L,z3); f4=diff(L,lamd); F(z1,z2,z3,lamd)=[f1;f2;f3;f4]; Jac(z1,z2,z3,lamd)=[diff(f1,z1) diff(f1,z2) diff(f1,z3) diff(f1,lamd);diff(f2,z1) diff(f2,z2) diff(f2,z3) diff(f2,lamd);diff(f3,z1) diff(f3,z2) diff(f3,z3) diff(f3,lamd);diff(f4,z1) diff(f4,z2) diff(f4,z3) diff(f4,lamd);]; x0=[0;0;0;0];delta=zeros(1,3); for i=1:3 x1=x0-inv(Jac(x0(1),x0(2),x0(3),x0(4)))*F(x0(1),x0(2),x0(3),x0(4)); x0=x1; delta(i)=abs(vpa(g(x0(1),x0(2),x0(3)))); end delta 所输出的结果为： delta = 1.509157317239845 0.069487094565626 0.000009862142070 即，迭代三次后，所得到的坐标能够使g(z1,z2,z3)=0.000009862142070，可以近似认为这一点位于极限状态曲面上。 该点坐标为： z1= -6.8622094451667668002649151995478（对应σ） z2= 4.5377285391719876172778985750505（对应T） z3= -2.1510122276458808636524699657479（对应D） 则该杆件的可靠性指标β=8.503396519511131265425973854947。 讨论与心得