结构可靠度分析与设计的编程实践1.docx

UndergraduateCourse"Loads&StructuralDesignMethods"
Project#2

结构可靠度分析与设计的编程实践
姓名:杨冬冬学号：1103300607班级：1033119

摘要：本文使用MATLAB软件，运用设计点法、RF法、蒙特卡洛方法等几种求可靠度指标的方法对相关习题进行了编程计算，得到了与预期结果十分接近的结果。
引言
科学实验和力学的发展是土木工程由完全凭直觉和经验走向了科学与经验相结合的道路。土木工程的设计理论和方法的发展经历了几个阶段：容许应力设计法到极限强度设计法到极限状态设计法。极限状态设计法引入可靠性的概念，应用概率论和数理统计的知识进行设计，经历了四个水准：半经验半概率、近似概率、全概率、最优概率，使结构的设计更加科学和可控。今天，国际上结构设计理论的发展趋势是概率极限状态设计法。在结构建造和使用过程中，结构可靠与不可靠是不可预知的，这是因为建造和使用中存在了诸多不确定性。而结构的可靠度关乎人们的生命财产安全，至关重要，因此学习计算结构的可靠度很有必要。

结构可靠度分析的基本原理
结构可靠度指的是结构在规定的时间，在规定的条件下，完成预定功能的概率，而通过极限状态的概率计算可以得到安全概率和失效概率，而且。而结构可靠度以结构可靠指标计量，定义为。通过结构可靠指标对结构的可靠度进行判别有着较为科学的依据，可靠度指标与失效概率相关，表明其具有概率意义。在取值中，可靠指标取值越大，则结构的失效概率越小；结构的安全概率越大，结构的保证率越大，更加安全。
当用概率描述结构的可靠性时，结构可靠度是一个概率计算问题，即计算结构可靠度就是计算结构在规定时间内、规定条件下结构能够完成预定功能的概率。从简单到复杂或精确程度的不同，先后提出的可靠度计算方法有一次二矩法，二次二矩法，蒙特卡洛方法及其他方法。一次二矩法又分为中心点法和设计点法，其中验算点法是目前可靠度分析最常用的方法。
对于一般形式的结构功能函数。式中为n个相互独立随机变量，其平均值和标准差是已知的。在计算=/时，需要计算功能函数Z的平均值和标准差，如果进一步已知各随机变量概率密度函数，精确计算这两个值要进行高维数值积分，这又回到了与直接根据随机变量概率分布计算可靠度相同的问题；如果各随机变量的概率密度函数是未知的，则不能给出和的近似结果，然而展开点的选择却是一个值得考虑的问题，一个简单的方法是将展开点选为平均值处，即中心点。功能函数Z的平均值和方差近似为

从而，按中心点法计算的可靠指标为

由上式可以看出，中心点法使用了结构功能函数的一次泰勒级数展开式和随机变量的前两阶矩（平均值和方差），故称为一次二阶矩方法。但是该方法存在很多缺点：1.功能函数在平均值处展开不尽合理；2.对于力学意义相同但数学表达式不同的结构功能函数，由此计算的β可能不同；3.没有考虑随机变量的概率分布。
于是产生了改进的一次二阶矩方法也就是设计点法：在极限状态曲面上寻找验算点，并在此基础上进行泰勒级数展开，应用随机变量的前二阶矩，采用非正态随机变量的当量正态化，迭代求解结构的失效概率的一种方法，该方法简称验算点法。
RF法：RF法用于处理随机变量非正态分布的情况，为JCSS（JointCommitteeofStructuralSafety）推荐的方法，所计算的可靠指标又称Rackwitz—Fiessler可靠指标。

JC法的基本原则：
–按照等效正态化原则将非正态随机变量转化为当量正态化随机变量
–解决由于非正态随机变量导致的可靠指标与失效概率不一一对应的不足
–转化之后可以利用验算点法求解
蒙特卡洛方法：对功能函数中所有随机变量进行随机抽样，根据随机变量的抽样值进行结构功能函数值Z=g(x)的计算，根据结构可靠度基本理论可知当功能函数Z>0时，结构可靠，当功能函数Z=0时，结构处于极限状态，当功能函数Z<0时，结构失效；若进行N次随机模拟抽样计算，功能函数小于零的次数为，则结构的失效概率的估计值为：
结构可靠性设计的基本原理
如果建立了结构或者结构构件的功能函数，并且已知荷载效应及抗力的概率分布及平均值和变异系数，可计算结构或者结构构件的可靠指标β。同样，如果已知荷载效应及抗力的概率分布、平均值和变异系数及目标可靠度指标βT，则可计算抗力平均值，进而确定结构或者结构构件的抗力标准值，对结构或者结构构件进行设计。
结构可靠度设计实际上是结构可靠度分析的逆过程，可靠度分析是已知荷载和抗力的概率分布和统计参数求可靠指标，而可靠度设计是已知和载荷目标可靠度求抗力。对于可靠度设计，除抗力服从对数正态分布的情况外，一般不能直接迭代计算抗力平均值，这时可以先假定抗力平均值，按结构可靠度分析的方法计算可靠指标β，当计算的可靠度指标β大于目标可靠度指标βT