结构可靠度matlab,基于MATLAB的结构可靠度分析

·181· 居 舍 研究探讨 2018年6月(上) 结构可靠度即是结构在规定时间内、规定条件下完成预定功能的概率 [1]。规定时间是指结构的设计基准期，规定条件是指预先确定的施工条件和适用条件。本文利用 MATLAB 软件，以某薄壁型钢梁为例，编制了用于结构可靠度计算的改进一次二阶矩法即设计验算点法计算程序，运行结果表明计 算速度较快，结果合理，可供工程技术人员参考使用。 1 结构可靠度的理论 现有结构可靠度计算方法中，关于一次二阶矩方法主要包括中心点法和设计验算点法，而对于非正态变量的处理方面则有 JC 法、映射变换法、实用分析法等，其中设计验算点法由于其随机变量为正态分布，对于功能函数线性或非线性程度不高的单一结构，计算简便、效果理想 [2，3]。 1.1 结构可靠度与失效概率 为了便于理解结构可靠度与失效概率，选取简单的功能函数来进行说明。假定一个结构的功能函数仅与荷载效应 S 和结构抗力 R 有关，且假定随机变量 S、R 都为正态分布。那么该结构的功能函数为 Z=g(R，S)=R-S，对应的极限状态方程为 ：Z=R-S=0。 若结构抗力 R 大于荷载效应 S，即 Z=R-S ＞0。则结构处于可靠状态，反之结构处于失效状态。由于 R、S 的均值和标准差分别为μR、μS 和σ R、σS，则 Z 的均值为μZ=μR-μS，标准差 σZ= 。令结构的失效率为 Pf，结构的可靠度 Pr。由概率论可知 P(Z<0)+P(Z>0)=1，即失效概率和可靠度的关系为 Pf+Pr=1。 对于非正态分布的随机变量，可以根据随机变量的概率分布函数通过积分得到结构的可靠度和失效概率。 1.2 结构可靠度与可靠指标 同样为了方便说明，我们以正态分布随机变量为例进行说明。首先把 Z 的正态分布 N(μ，σZ)转换为标准正态分布 N(0，1)，之后定义可靠指标β(无量纲)，并令β=μZ/ σZ。得到失效概率 Pf=φ(-β)。利用可靠度与失效概率的关系，得到可靠度与可靠指标之间的关系为：Pr=1-Pf=1-φ(-β)=φ(β)。由此可知，可靠指标β与可靠度 Pr 成正相关，与失效率为 Pf 成反相关。 若随机变量为非正态分布，则可以通过计算 Z 的均值 μZ 和标准差 σZ，再由 β=μZ/σZ 计算出近似的可靠指标。 2 在MATLAB中实现设计验算点法 在运用设计验算点法进行结构可靠度分析之前，为了减少大量繁琐的有限元计算工作，更为了说明设计验算点法的应用，因此我们选用功能函数已知的算例来分析。 受永久荷载作用的薄壁型钢梁，Z=Wf-M=0。已知截面 抵抗矩W 的均值μW=54.27，变异系数Vw=0.05 ；钢材强度 f 的均值μf=3800，变异系数Vf=0.08 ；弯矩 M 的均值μM=1.3×105，变异系数 VM=0.07。所有变量均服从正态分布。 求钢梁的失效概率 pf。 解：功能函数的梯度为▽ g(W，f，M)=(f，W，-1)T。 MATLAB 源程序如下 ： clear ；clc ；(清除工作空间里的变量数据，清除命令窗口的内容)muX=[54.72 ；3800 ；1.3e5] ；cvX=[5 ；8 ；7]/100 ；(输入均值、变异系数) sigmaX=cvX.*muX ；(计算标准差)x=muX ；normX=eps ；whileabs(norm(x)-normX)/normX>1e-6(循环语句)normX=norm(x)；g=x(1)*x(2)-x(3)；(功能函数)gX=[x(2)；x(1)；-1] ；(梯度函数) gs=gX.*sigm