第3章结构可靠度计算.pptVIP

东南大学 第3章 结构可靠度计算 杨春侠 长沙理工大学土木工程学院 3.1 中心点法 一、计算公式 二、公式的应用 例题3-1 一圆截面直杆，承受拉力P=100kN，基本变 量及其统计参数如下，求此杆的可靠指标。 材料强度设计值 fy 杆的直径d 功能函数为 解： 三、中心点法的评价 1、中心点法的优点 直接给出?与随机变量统计参数之间的关系，不必知道基本变量的的真实概率分布，只需知道基本变量的统计参数即可计算可靠指标值； 若值β较小，即Ｐf 值较大时，Ｐf 值对基本变量联合概率分布类型很不敏感，由各种合理分布计算出的Ｐf 值大致在同一个数量级内； 对正常使用极限状态尤为适用 ( ? ＝1～2)。 2、中心点法的缺点： 算得的?与结构功能函数的形式有关； 没有考虑基本变量的实际分布； 对于非线性功能函数，近似在基本变量的平均值处按泰勒级数展开不尽合理； 当失效概率pf 10-3时(相应的? 3.09)， Z的分布类型不同将会导致pf 值在几个数量级的范围内波动，中心点法失效。 3.2 验算点法 一、验算点法的两点改进 R. Rackwitz和 B.Fiessler 等人提出验算点法。 （１）当功能函数Z为非线性时，不以通过中心点的超切平面作为线性近似，而以通过Z＝0上的某一点P* (x1*, x2*,····· , xn*)超切平面作为线性近似，以避免中心点方法中的误差。 （２）当基本变量xi 具有分布类型的信息时，将Xi 的分布在 验算点P* 处以与正态分布等价的条件，变换为当量正态分布，从而在β中合理地反映了分布类型的影响。 二、计算公式 3、U空间的可靠指标 可靠指标的几何意义 4、X 空间的验算点 三、当量正态化 2、对数正态分布当量正态化 常见概率分布类型 2、正态分布 3、标准正态分布 4、对数正态分布 5、伽马分布 6、极值Ⅰ型分布 7、极值Ⅱ型分布 8、极值Ⅲ型分布 9、极值Ⅲ型分布（威布尔分布） 四、验算点法计算步骤 习题 3-2 3.3 蒙特卡洛方法 蒙特卡洛法又称为随机抽样法，是一种基于“随机数”的计算方法。 源于美国在第二次世界大战中研制原子弹的“曼哈顿计划”。该计划的主持人之一、数学家冯·诺伊曼用驰名世界的赌城Monte Carlo来命名。 一、蒙特卡洛法的基本思路 针对实际问题建立一个简单且便于实现的概率统计模型，使所求的量恰好是该模型某个指标的概率分布或者数字特征。 对模型中的随机变量建立抽样方法，在计算机上进行模拟测试，抽取足够多的随机数，对有关事件进行统计。 对模拟试验结果加以分析，给出所求解的估计及其精度(方差)的估计。 二、随机数的生成 蒙特卡罗模拟的关键是生成优良的随机数。 在计算机实现中, 是通过确定性的算法生成随机数,所以这样生成的序列在本质上不是随机的,只是很好的模仿了随机数的性质。我们通常称之为伪随机数。 在模拟中，我们需要产生各种概率分布的随机数，而大多数概率分布的随机数产生均基于均匀分布U(0,1)的随机数。 可靠度计算中采用一般抽样法和重要抽样法。 三、蒙特卡洛方法计算可靠度的步骤 1、列出功能函数及基本变量的概率分布函数 2、估算抽样的次数 3、产生的随机数 4、计算基本变量 5、计算功能函数值，记录功能函数小于0的次数 ； 6、计算失效概率估计值 对于结构极限状态函数中包含多个正态或非正态基本变量，只要知道其概率分布类型及统计参数，就可采用迭代法计算?和验算点P*的坐标值。 已知:Xi(i=1,…，n)的分布类型及统计参数?Xi,σXi，极限状态方程g(X1,…,Xn)=0 假定验算点P*的坐标值初值 对于非正态变量Xi，根据公式(4),(5)或(6),(7)进行当量正态化，计算 以代替 将设计验算点P*的坐标Xi*代入极限状态方程 以求出? 计算Xi*并取代上一次验算点坐标，返回进行下一轮迭代 本次求得的?和Xi*即为所求的可靠指标和设计验算点P*的坐标值 是 否 |上次求出的?－上次求出的?| ≤允许误差? 承受恒载和楼面活荷载的钢筋混凝土轴心受压短柱，已知恒载产生的轴向力NG为正态分布，活载产生的轴向力NL为极值Ⅰ型分布，截面承载能力（抗力）R为对数正态分布，统计参数分别为?NG＝1159.1kN， ?NG＝81.1kN， ?NL＝765.5kN， ?NL＝222kN，?R＝4560kN，?R＝729.6kN，极限状态方程为 Z＝R－NG－NL＝0，求可靠指标?和设计验算点。 ?＝3.96 R*=3009.8 NG*=1194.1 NL*=1815.6 g(·)= R*- NG*- NL*=0.01≈0 结构可靠度计算 Chan