工程结构可靠度.doc

结构可靠度分析与设计的编程实践
摘要:本文基于结构可靠度的分析方法(改进一次二阶矩法,JC法,蒙特卡洛法)的基本原理,采取不同的算例利用MATLAB软件编制了相应计算程序分析结构可靠度。同时,考虑随机变量的分布信息,还对结构构件可靠度校准及分项系数校准等进行了简单分析。通过此次编程实践,加深了对结构可靠度重要性的理解。
引言
结构的可靠度是结构在规定时间内,规定条件下结构能够完成预定功能的概率。由于材料生产的不均匀性和结构作用荷载的随机性,导致抗力和荷载都属于符合某种情形的分布,通过概率论知识并依托Matlab数学分析软件我们可以计算出结构的可靠指标和失效概率,确定结构的安全程度。从简单到复杂或精确程度的不同,先后提出的可靠度计算方法有一次二阶矩方法,高次二阶矩方法,蒙特卡洛方法及其他方法。一次二阶矩又分为中心点法和验算点法,其中验算点法是目前可靠度分析最常用的方法。
2. 结构可靠度分析的基本原理
结构的可靠度通常受各种荷载、材料性能、几何参数、计算公式精确性等因素的影响。这些因素一般具有随机性,特别是当结构可能受到多种因素的影响,而且每一种影响不一定完全服从正态分布时,需要预先对它们进行当量正态化处理。用失效概率来描述结构的可靠度,在理论上是合理的,但在具体计算上又非常困难,因此,必须寻求近似的或简化的计算方法,以达到实用的目的。在实际工程中,结构的功能函数往往是由多个随机变量组成的非线性函数,而且这些随机变量并不都服从正态分布或对数正态分布,因此不能直接采用前述公式计算可靠指标,而需要作出某些近似简化后进行计算。下面本文将介绍分析结构可靠度的几种常用方法。
中心点法
若结构功能函数为线性功能函数,即
式中:为已知常数(i=1,2,3⋅⋅⋅⋅⋅⋅n)。
由此可知功能函数的统计参数为
根据概率论中心极限定理,当随机变量的数量n较大时,可以认为Z近似服从于正态
分布,则可靠指标直接按下式计算
然后再求得得结构的失效概率:
若结构的功能函数为非线性情况,即
将Z在随机变量的平均值(即中心点)处按泰勒级数展开,并仅取线性项,即
接下来的算法就参照线性功能函数的情况。
验算点法
JC法是国际“结构安全度联合委员会(JCSS)”推荐的一种计算结构可靠指标更为一般的方法。起主要特点是,对于非线性的功能函数,线性化近似不是选在中心点处,而是选在
失效边界上,即以通过极限状态方程上得某一点的切平面做线性近似,以提高可靠指标的计算精度。这里特定的点称为设计验算点,它与结构最大可能的失效概率相对应,因此也叫验算点法。
当功能函数为线性函数时:
式中:为已知常数(i=1,2,3⋅⋅⋅⋅⋅⋅n)。
为进一步在标准正太坐标系中研究可靠指标的几何意义,按下式将随机变量变换为标准正态随机变量
则结构功能函数可由表示为:
从而功能函数的平均值和标准差为:
这样可靠度指标为:
由上式计算的可靠指标与结构失效概率的对应关系:
当功能函数为非线性函数时:
为计算可靠指标,将非线性功能函数在验算点处泰勒级数展开,则结构功能函数的一次展开式为:
则Z的平均值和标准差可分别近似表示为
由此可求得可靠指标:
JC法
当随机变量不服从正态分布时,则需要将其进行当量正态化将随机变量等效为正态随机变量,当量正态化的条件是,在验算点处使非正态随机变量的概率分布函数值与当量正态随机变量的概率分布函数值相等,的概率密度函数值与当量正态随机变量的概率密度函数值相等,即
则:
当量正态化后:
上述三式构成非正态随机变量情况下可靠指标的迭代计算公式。
3. 结构可靠性设计的基本原理
结构设计首先要保证结构的安全性,其次是保证结构的适用性、耐久性和偶然作用下的整体稳定性。为实现这些目的,就要对结构进行合理的设计。结构的设计由最初的完全经验设计法,逐步向以科学试验和理论计算为基础的方向发展。从力学计算理论的发展过程看,经历了容许应力设计法,破损阶段设计法和极限状态设计法;从概率理论的应用过程看,经历了定制设计法、半概率设计法和近似概率设计法。现在。国际上结构设计理论的发展趋势是概率极限状态设计法。
结构可靠度校准
在现实工程中,按照结构设计规范规定的设计表达式进行可靠度计算,按规范设计的结构或结构构件的可靠度指标。在这种情况下反映的是一种结构或一类结构构件的可靠指标,代表了设计规范的可靠水平,这一过程称为结构可靠度校准。
结构可靠度校准法计算步骤如下:
(1) 假设一个荷载效应比;
(2) 确定构件抗力的特征值:
(3) 确定基本变量的均值和标准差:
均值:
标准差:
(4) 确定极限状态方程:
(5) 利用PF法求解可靠度指标;