一次二阶矩法ppt课件

ch4一次二阶矩法 引言一次二阶矩法是求解非时变荷载作用下结构可靠度问题的行之有效的近似方法 它既有较高的精度 又有较高的计算效率 一次二阶矩法是在基本变量xi i 1 2 n 的概率分布尚不清楚中时 采用只有均值 又称为一阶原点矩 和标准差 又称为二阶中心矩 的数学模型去求解结构可靠度的方法 Ch4一次二阶矩法 4 1一次二阶矩中心点法4 2一次二阶矩验算点法4 3JC法 4 1一次二阶矩中心点法 4 1 1 简单问题可靠指标的求法4 1 2 中心点法基本原理4 1 3 中心点法举例4 1 4 对中心点法的评述 在介绍中心点法之前 先回顾功能函数为线性函数 随机变量为正态分布情况下的可靠指标问题 I随机变量为正态分布 且功能函数为线性函数 假定抗力R和荷载效应S均服从正态分布 对于功能函数Z R S 由于Z是R S的线性函数 根据正态随机变量的特性 Z也服从正态分布 其平均值及标准差分别为 则定义可靠指标 为什么可以定义来衡量结构的可靠程度 4 1 1简单问题可靠指标 或失效概率 的求法 正态分布功能函数Z 其失效概率与可靠指标之间的精确关系 这种情况下可靠指标是唯一的 且与失效概率之间有精确的对应关系 II随机变量为对数正态分布 功能函数为线性函数 对数正态分布的定义 若lnR lnS服从正态分布 则称R S服从对数正态分布 假定抗力R和荷载效应S均服从对数正态分布 且结构功能函数可表示为Z lnR lnS 由于lnR lnS均服从正态分布 Z也服从正态分布 是lnR lnS的线性函数 根据正态随机变量的特性 其平均值及标准差也可以精确得到 可靠指标为 根据x的统计参数 求lnx的统计参数 结论 随机变量为正态分布 或对数正态分布 且功能函数为线性函数的条件下 结构的可靠指标或失效概率可根据随机变量的统计参数方便的求出 随机变量为正态分布 由其形成的线性函数也服从正态分布 4 1 2中心点法 引言 假定随机变量x1 x2 xn服从任意分布 功能函数不是线性函数 这时 精确求解Z的平均值和标准差是非常困难的 即便能够求得 Z也不服从正态分布 也不能用上面方法来计算结构的可靠指标 若将非线性功能函数作为泰勒级数展开 并取其一次展开式 前两项 但一次展开式已不是原来的功能函数 所计算可靠指标与结构失效概率之间不再存在精确的对应关系 在这种情况下如何选择展开点 从而使近似计算结果与精确失效概率的误差最小 成为一次二阶矩法要研究的问题 4 1 2中心点法 基本原理 设结构的极限状态方程为将极限状态函数在中心点M 处展开为泰勒级数 并作线性化处理 得根据概率论中随机变量参数估计 Z 的统计参数为 结构的可靠指标 在中心点处展开为泰勒级数 4 1 3中心点法举例 钢梁承受确定性弯矩 截面塑性抵抗矩W和屈服强度f都是随机变量 已知统计参数为 W f 试用中心点法计算该钢梁抗弯的可靠指标 4 1 3中心点法举例 解法1取功能函数极限状态方程为 4 1 3中心点法举例 解法2取功能函数极限状态方程为 4 1 4对中心点法的评述 与随机变量为正态分布 功能函数为线性函数的情况下求可靠指标 形式一样 却有本质的区别 1 中心点法不论R S的概率分布如何 直接用其平均值 标准差求可靠指标 事实上 平均值 标准差相同的随机变量有无穷多个 每一种情况的可靠度都是不同的 而这里是同一个结果 2 前者的可靠指标是唯一的 并与失效概率有精确的对应关系 而中心点法可靠指标与建立的功能函数表达式的形式有关 4 1 3对中心点法的评述 中心点法的主要弱点没有考虑基本变量的概率分布均值 方差及可靠指标的计算式是误差传递公式同一个结构往往可以列出几种等价的极限状态方程 不同的极限状态函数在运用中心点法计算时 其结果可能不一致 将非线性功能函数在随机变量的平均值处展开不合理 由于平均值不在极限状态曲面上 展开后的线性极限状态面可能会较大程度地偏离原来的极限状态曲面 基本变量不服从正态分布和对数正态分布时 计算出的结构可靠指标与结构的实际情况出入较大 中心点法的线性示意图 该法选用的线性化点 即平均值点 不在失效边界上 4 2验算点法 改进的一次二阶矩法 4 2 1验算点法基本思想运用泰勒级数进行线性化处理时 略去高阶项的误差随线性化点到失效边界的距离增加而增大 中心点法运用泰勒级数在中心点处展开进行线性化处理 选择在失效边界Z 0的某点处将极限状态函数展开为泰勒级数 以避免这种形式的误差 验算点法 4 2 1验算点法基本思想 关于设计验算点设点P x1 x2 xn 为极限状态方程Z 0所对应的曲面上的点 d P M 为点P到中心点M 的距离 则能使mind P M 的点P 称为设计验算点 简称为验算点 记为P x 1 x 2 x n 显然验算点的坐标满足 验算点法示意 M 验算点法 4 2 2设计验算点法求可靠指标 理论推导 确定验算点位置进而求可靠指标当线性化点选在设计验算点xi i 1 2 n 上时Z的均值为 求Z的数学期望得 由于设计验算点在失效边界上 故有则有 假设各随机变量独立 则可求解Z的方差 引入分离函数式 将上面的根式线性化 得表示第 个随机变量对整个标准差的相对影响 因此称为灵敏系数 4 2 2设计验算点法求可靠指标 4 2 2设计验算点法求可靠指标 根据可靠指标的定义 有由于 必有 对于所有i 将乘到分母上整理后得 4 2 2设计验算点法求可靠指标 式 13 代表n个方程 再加上 14 共有n 1个方程 未知数有和 也是n 1个 可联立求解 解出设计验算点P 及相应的 常用迭代法求解 4 2 2设计验算点法求可靠指标 计算步骤 1 选取设计验算点坐标的初值 一般取 2 由式 11 计算的值 其中包括 3 由式 13 得到和的关系 4 由式 14 解出值 5 将该值代入式 13 求出新值以该新重复进行 2 5 计算 直到值与上次相等或误差不超过允许值 此时即为所求的可靠指标 即为所确定的设计验算点坐标 4 2 3验算点法举例 设极限状态功能函数为 均为随机变量 并具有如下统计信息 试用验算点法计算结构的可靠指标 1 取均值作为验算点初值 2 计算 0 2639 0 7037 0 6597 3 与的关系 从而有新的 4 由 14 求及新的 4 2 3验算点法举例 5 以新进行第二次迭代 重复4直到值与上次相等或误差不超过允许值 此时即为所求的可靠指标 即为所确定的设计验算点坐标 本例第三次迭代求得 第四次迭代求得 即为最终结果 设计验算点法求可靠指标评述 验算点法对极限状态方程中服从正态分布的随机变量计算结果尚可 而非正态分布误差较大 JC法求可靠指标 问题的提出 验算点法对极限状态方程中服从正态分布的随机变量计算结果尚可 而非正态分布误差较大 在工程结构可靠度分析中 永久荷载一般服从正态分布 风压 雪压 楼面活载服从其它类型 如极值 型等 截面抗力R服从对数正态分布 因此在极限状态方程中 常包含非正态变量分布的基本变量 对于这种极限状态方程的可靠度分析 一般要把非正态变量进行当量正态化 JC法的基本概念就是在应用前面所述方法 验算点法 时 将非正态的随机变量先行 当量正态化 JC法是由Rackwitz Fiessler Hasofer Lind等人先后提出来的 因为国际安全度联合委员会 JCSS 推荐采用这个方法而得名 JC法求可靠指标计算步骤 设结构的极限状态方程为Z g x1 x2 xn 0 基本变量xi的统计参数均值为标准差为 i 1 n 则JC法的具体实施步骤加下 l 假定基本随机变量的设计验算点P 的坐标值 一般取为 2 对基本变量xi进行当量正态化处理 计算其当量正态分布yi的统计参数 并用来代替和并记为 1 n 3 计算的值 其中包括 JC法求可靠指标计算步骤 4 确定和的关系 5 求 6 求出新值 7 以该新重复进行 2 5 计算 直到值与上次相等或误差不超过允许值 此时即为所求的可靠指标 即为所确定的设计验算点坐标 JC法计算实例 某轴心受压短柱 承受永久荷载产生的压力SG 汽车 人群可变荷载产生的压力SQ 截面抗力为R 各变量的统计信息如下表 试用JC法求可靠指标及对应的失效概率 JC法计算实例 1 建立极限状态方程Z g R SG SQ R SG SQ 02 假定初始的设计验算点 由已知条件得 x1 R 365 8 x2 SG 63 x3 SQ 84 03 将非正态