基于matlab算法的可靠度分析,参考基于matlab算法的可靠度分析

结构可靠度基于Matlab算法性能比较

研究生课程考核试卷

科 目： 工程结构可靠度 教 师： 范文亮

姓 名： 李亚勇 学 号： 20111602138

专 业： 岩土工程 类 别： 学术

上课时间： 2012 年 4 月至2012 年 6 月

考 生 成 绩：

卷面成绩

平时成绩

课程综合成绩

阅卷评语：

阅卷教师 (签名)

重庆大学研究生院目录

1 绪论1

1.1概述1

1.2 国内外研究现状1

2 结构可靠度基本理论2

2.1 可靠性概念2

2. 2 结构可靠度与失效概率4

3结构可靠度分析方法5

3.1 一次二阶矩法5

3.1.1 验算点法(JC法)5

3.2 Monte Carlo 抽样法6

3.2.1 直接蒙特卡罗法6

3.2.2 重要抽样蒙特卡洛法7

3.3 一次渐近积分法8

3.4 Breitung法8

4 可靠度计算在Matlab环境下的实现与比较9

5 结论15

参考文献16

2

结构可靠度基于Matlab算法性能比较

摘要：本文对只有一个失效模式的构件可靠度问题的各类算法进行比较，并对各种算法(JC法、Breitung法、一次渐进法、重要抽样蒙特卡洛法、蒙特卡洛法)发展历史和推导过程做出论述。基于结构可靠度分析方法(一次二阶矩法，二次二阶矩法，蒙特卡洛法)的基本原理, 采取不同的算例利用MATLAB软件编制了相应计算程序，并对计算结果进行了比较。本文算例对比主要分为三大类：1.线性功能函数与非线性功能函数的对比；2.相同的功能函数相同的分布不同的参数；3.相同的功能函数相同的参数不同的分布。并通过对算例的比较得出每种方法的适用范围和精确程度，为以后工程中解决可靠度问题提供了理论和实例依据。

关键词：结构可靠度；Matlab;JC法；Breitung法；一次渐进法；重要抽样蒙特卡罗法；蒙特卡罗法；

1 绪论

1.1 概述

为了保证结构、构件等物体的适用性和耐久性要求，工程人员要对构件物体的受力进行可靠度分析。可靠度分析采用容许应力法即对材料强度进行一系列折减(即当材料所能承受的最大力为F,在工程应用中我们将F除以一个大于1的系数n，采用除数作为构件所能承受的最大的力进行受力分析。但是这种方法有自身的缺陷性，例如对于一个受弯构建来讲，当力F增大2倍时，弯矩增大量却大于2倍。因此，用容许应力法做结构可靠度设计有自身的缺陷，现在工程中已经较少采用。随着概率论与数理统计的发展，分项系数法逐渐在可靠度分析中采用，分项系数法不像容许应力法那样只对应力进行一项折减，而是分别对不同的力不同的材料采取不同的分项系数进行可靠度计算，分项系数法现在在工程可靠度分析中广泛采用。由于材料生产的不均匀性和结构作用荷载的随机性，导致抗力和荷载都属于符合某种情形的分布，通过概率论知识并依托Matlab数学分析软件我们可以计算出结构的可靠指标和失效概率，确定结构的安全程度。

1.2 国内外研究现状

可靠度理论在国外开展的比较早，我国知道20世纪70年代才开始可靠度方面的研究，虽说我国开展的较晚但也做了大量的研究并提出了一些自己的理论，目前我国从事结构可靠度方面的研究仍然比较多。可靠度理论在结构工程方面的应用研究开始得较早，目前已广泛应用于水利、煤炭、通信等多种行业。早在1947年，前苏联的尔然尼钦【1】通过大量研究就提出了用一次二阶矩方法来估计结构的失效概率。但那以后的一段时间里，这种计算失效概率的方法并没有得到发展仍然停留在古典可靠性理论，而且只适用于随机变量符合标准正太分布的情况。其后，美国的C.A.Comell、A.H.s【2】和wH.Tang【3】在总结了大量的工程实践经验发展了工程技术中应用概率的概念和方法。其中CA.Comell于1969年提出了作为衡量结构安全度的统一标准—与结构失效概率相联系的可靠度指标，并建立了结构安全度的二阶矩模式。Lind【4】在总结前人经验的基础上提出了材料和构件分项系数的概念，将可靠度指标表达为设计人员习惯使用的分项系数形式。随后美国伊利诺斯大学A.H.S在结构可靠度方面做出了较大的进步，提出了广义可靠度概率法，并对各种结构的分布情况做了分析。1976年国际“结构安全度联合委员会”(Jcss)采用RackwitZ【5】和Fiessler【6】等提出的“当量正态”法以考虑随机变量实际分布的二阶矩模式，这种方法极大的提高了二阶距方法的精度。在此后一段时间可靠度理论越来越完善，在结构设计中的应用也越来越广泛。可靠度理论迎来的发展的黄金期。

20 世纪 70 年代以来，各个国家都开始了基于概率论与数理统计理论的可靠度分析体制和相应的设计规范的建立。加拿大【7】、联邦德国【8】、北欧五国、美国等国家先后建立了自己的可靠度设计规范，我国也在90年代建立了自己的可靠度设计规范。目前我国土木工作者依托大量的工作实践提出了自己对于可靠度的理解，为世界范围内可靠度理论的发展做出了自己的贡献。 从可靠度理论开始应用于结构工程方面，世界范围内的科学家依托自己对于可靠度理论的理解并结合相应的工程实践提出了许多根据可靠度理论计算可靠指标和失效概率的方法，例如：一次二阶距法(中心点法、JC法、实用分析法等)、二次二阶距法【9】、二次四阶矩法【10】、一次渐进法、蒙特卡洛法【11】(直接蒙特卡洛法、重要抽样蒙特卡洛法、方向抽样蒙特卡洛法等)等等方法。每种方法都有自身的实用范围和优缺点，要总结相应方法的适用范围，在不同的情况下采用不同的计算方法，获取更加精确的结果。

2 结构可靠度基本理论

2.1 可靠性概念

从事工程结构设计的基本目的，是在一定的经济条件下，赋予结构以适当的可靠度，使结构在预定的使用期限内，能满足设计所预期的各种功能要求。无论是房屋、桥梁、隧洞等结构，都必须满足以下四项基本功能要求【12】：

1.能承受在正常施工和正常使用时可能出现的各种作用；

2.在正常使用时具有良好的工作性能；

3.在正常维护下具有足够的耐久性能；

4.在偶然事件发生时(如地震、火灾等)及发生后，仍能保持必需的整体稳定性。

上述第1，4项为结构的安全性要求，第2项为结构的适用性要求，第3项是结构的耐久性要求。结构若同时满足安全性，适用性和耐久性要求，则称该结构可靠，即结构的可靠性是结构安全性、适用性和耐久性的统称。

所谓结构的耐久性能，是指在一定的期限内结构材料和主体结构没有发生过大的老化或者变形影响适用者的舒适度和房屋失效概率增加。

所谓结构的安全性，是指结构在遇到偶然荷载(例如50年一遇的地震、大风等)，结构不至于倒塌对人的生命财产造成一些损失。

可靠度的定义【13】：结构在规定的时间内，在规定的条件下，完成预定功能的能力，称为结构可靠性。对于特定情况而言，可具体解释成在规定的参照时期内，结构将不会达到某一特定极限状态的概率，即所谓的结构可靠度。结构可靠度是结构可靠性的概率量度。

上述“规定的时间”，一般指结构设计基准期，目前世界上大多数国家普通结构的设计基准期均为50年。由于荷载效应一般随设计基准期增长而增大(假如一个地区遇到7级地震的概率是20年一遇，那么设计基准期选为50年结构遇到7级地震的概率和次数就要比基准期为30年的结构要大，对结构造成的破坏也就更大)，抵抗能力随着基准期增长降低，而影响结构抗力的材料性能指标则随设计基准期的增大而减小，因此结构可靠度与“规定的时间”有关，“规定的时间”越长，结构的可靠度越低。

一般情况下，总可以将影响结构可靠度的因素归纳为两个综合量，即结构或结构构件的荷载效应S和抗力R。令

Z=g(R，S)=R—S

由于R和S都是随机变量，因此Z是一个随机变量，Z可能出现下列三种情况：

Z>0,结构可靠

Z<0,结构实效

Z=0,结构处于极限状态

结构可靠性设计既是求出Z大于0的概率即可靠度，Z小于0对应于结构的失效概率。对于结构抗力R由于构件材料、尺寸、制造工艺的差别，造成R并不是一个固定不变的数而是一个关于材料、尺寸、制造差异的一个随机变量。S也是一个随机变量(例如一个房屋所承受的雪荷载、风荷载等值并不是固定的一个数，我们可以通过对不同月份和年份的统计出雪荷载的大小，并根据数理统计原理统计出S的分布情况来)。结构的极限状态方程如下：

Z=R—S (1)

假设R和S为两个相互独立的正态随机变量，他们的均值和方差分别为μR，μS；σR，σS。根据概率论知识可以得到：

μZ=μR-μS (2)

σZ2=σR2+σS2 (3)

令β=μZ/σZ，称β为结构的可靠指标。β越小结构的失效概率越大；β越大，结构的失效概率越小，结构越安全。

2. 2 结构可靠度与失效概率

对于公式Z=R-S计算P(Z>0),在概率论与数理统计中先求出极限状态方程的联合分布密度再求出极限状态方程的Z>0的区域，对联合分布密度函数在该区域内积分求出极限状态方程的概率。对于Z<0即对联合分布密度在Z<0的区域内积分。对于结构可靠度也采用同上面相应算法表示为 Pr，fX(X1,X2,…Xn)为结构的联合密度函数，Z>0为联合分布函数的积分区域

(4)

相反，如果结构不能完成预定的功能，则称相应的概率为结构失效的概率，表示为Pf，即：

<