可靠性与可靠性灵敏度分析的Monte Carlo数值模拟法.ppt

3 Monte Carlo 可靠性分析 5 Monte Carlo 可靠性灵敏度分析 本文观看结束！！！ 5 Monte Carlo 可靠性分析 采用向量运算代替循环提高Matlab计算效率 4 可靠性灵敏度 可靠性灵敏度定义为基本随机变量分布参数的变化引起失效概率变化的比率，在数学上可靠性灵敏度是由失效概率Pf对基本随机变量分布参数θx的偏导数予以表达，即。可靠性灵敏度反应了基本变量分布参数对失效概率的影响程度，无量纲正则化的可靠性灵敏度可以给出基本变量分布参数对可靠度的重要排序。 5 Monte Carlo 可靠性灵敏度分析 可靠性灵敏度为失效概率 Pf 对基本随机变量 xi 的分布参数 θ (i=1, 2, …, n; k=1, 2, …, mi)，其中mi为第i个变量xi的分布参数的总数) 的偏导数，将失效概率的积分公式对分布参数 求导数，便可以得到可靠性灵敏度 如下所示： 5.1 Monte Carlo 可靠性灵敏度估计的计算公式 将可靠性灵敏度定义式做如下变换，可使可靠性灵敏度变成数学期望的形式，之后就可以采用 Monte Carlo 数值模拟来估计可靠性灵敏度。 5 Monte Carlo 可靠性灵敏度分析 采用样本均值代替总体均值 xj是按联合概率密度函数fX(x)抽取N个样本中的第j个样本。 5 Monte Carlo 可靠性灵敏度分析 当 x为相互独立的n维正态随机变量时 可靠性灵敏度为 5 Monte Carlo 可靠性灵敏度分析 2. Monte Carlo 可靠性灵敏度估计的方差分析 对可靠性灵敏度估计值计算表达式两边求数学期望 由于样本点xj与母体x独立同分布，所以可以用样本的均值近似代替总体的期望，于是可以得到Monte Carlo 可靠性灵敏度分析结果的数学期望为 故Monte Carlo 可靠性灵敏度估计值为无偏估计。 5 Monte Carlo 可靠性灵敏度分析 用样本方差近似代替母体方差，可得Monte Carlo 可靠性灵敏度估计值得方差如下 而估计值得变异系数为 5 Monte Carlo 可靠性灵敏度分析 5 Monte Carlo 可靠性灵敏度分析 5 Monte Carlo 可靠性灵敏度分析 5 Monte Carlo 可靠性灵敏度分析 5 Monte Carlo 可靠性灵敏度分析 5 Monte Carlo 可靠性灵敏度分析 谢 谢 欣 赏！ 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 东北大学机械设计及理论研究所 机械系统可靠性和可靠性灵敏度分析的Monte Carlo 数值模拟法 授课教师：黄贤振 图1 机械系统的分类 机械系统是指由多个机械基本要素（如机械零、部件，机构，机器等）组成的，可以完成所需动作(或动作过程)，实现动作的传递和力的变换以及机械能的转化和利用的装置。 1 机械系统 2 机械可靠性 应力分析 应力分布参数 强度分析 强度分布参数 载荷 分布参数 机械系统 构件分布参数 机械系统 性能分布参数 机械系统 构件分布参数 f(x) fR(xR) fS(xS) x 机械可靠性模型 应力是指对产品功能有影响的各种因素，比如机械应力、变形、磨损等。 强度是指产品承受应力的能力，比如机械强度、刚度、抗裂度等。 假定机械系统的强度随机变量 xR，应力随机变量为 xS，机械系统极限状态函数（功能函数）为 极限状态的定义为：整个机械系统的一部分超过某一特定状态就不能满足设计规定的某一功能要求，此特定状态为该功能的极限状态。 2 机械可靠性 f(x) fR(xR) fS(xS) x xS dxS 应力xS落在附近dxS小区间内的概率为 在应力xS落在小区间dxS内的条件下，强度xR小于应力xS的概率为 应力xS落在小区间dxS内，同时强度xR小于应力xS的概率为 根据可靠度的定义，对于应力xS所有的可能值强度xR均小于应力xS的概率，及事件(xRxS)的边缘概率就是机械系统的可靠度(对上式的求和) 2 机械可靠性 同理，也可以求得失效概率的另一种表达式 2 机械可靠性 具有n设计变量的机械系统的功能函数可以表示为 则极限状态方程g(x1, x2, …, xn)=0将结构的基本随机变量空间分为失效域和可靠区域两部分。 2 机械可靠性 式中, f(xi) (i=1, 2, …, n)为随机