机械可靠性设计分析方法.ppt

机械可靠性设计 南京航空航天大学机电学院 于明礼 (a)根据提出的问题确定各变量之间的确定性函数关系。 (b)根据提出的问题构造一个简单、适用的概率模型或随机模型，使问题的解对应于该模型中随机变量的某些特征(例如概率、均值和方差等) 。 (c)根据模型中各个随机变量的分布，在计算机上产生随机数，实现一次模拟过程所需的足够数量的随机数。通常先产生均匀 分布的随机数，然后生成服从某一分布的随机数。 二、蒙特卡洛模拟求解步骤: (d) 根据概率模型的特点和随机变量的分布特性，设计和选取合适的抽样方法，并对每个随机变量进行随机抽样。这里的抽样方法有直接抽样、分层抽样、相关抽样、重要抽样等。 (e) 按所建立的模型进行仿真计算，求出问题的一个随机解。 (f) 统计分折模拟试验结果，给出问题的概率解以及解的精度估计。 在可靠性分析和设计中，用蒙特卡洛方法可以确定复杂随机 变量的概率分布和数字特征；可以通过随机模拟估计系统和零件 的可靠度；也可以模拟随机过程、寻求系统最优参数等。 二、蒙特卡洛模拟求解步骤: 对数正态分布 是标准正态分布抽样 正态分布 标准正态分布 指数分布 均匀分布 1 [0，1]均匀分布 或 密度函数f(x)或f(t) 分布名称 常见分布函数随机变量的随机抽样公式 例4-5某铝合金板的形状如图所示。受弯矩作用，其尺寸H、h、ρ均服从正态分布，分布参数为： 试确定理论应力集中系数ασ 的分布类型及分布参数。 M M 受弯矩作用的铝合金板 三、蒙特卡洛模拟算例 蒙特卡洛模拟算例的程序框图 开始 输入H、h、 ρ分布的类型及参数;N j=1 分别从H、h和ρ的分布中产生随机数Hf , hf , ρf 计算ασ 的随机数ασj j=N 进行分布类型判断、估计分布参数 输出ασ 的分布类型和分布参数 结束 j=j+1 解：理论应立集中系数ασ的计算公式为 蒙特卡洛模拟算例 * * §4.3 一次二阶矩方法求可靠度_工程方法 第四章 机械可靠性设计分析方法 §4.1干涉面积法 §4.2 分布代数 §4.4 蒙特卡洛模拟方法 §4.5 变异系数传递规律 从可靠度计算的普遍方程可以看出，对于应力和强度比较复杂的分布，由于积分困难，往往难以得出问题的解析解。 因此如何采用一些较好的近似方法，能比较方便地求得满足工程精度要求的零件可靠度的近似解，一直是人们探讨的一个问题。 3-1应力和强度概率密度曲线的干涉面积 O r, s f(s) g(r) s0= r0 f(s) g(r) a1 a2 应力和强度两个概率密度函数的交叉区，即干涉区阴影面积的大小，反映了零件或结构可靠度的高低。 该面积越小，可靠度越高，反之，可靠度越低。 干涉面积的大小是由两个概率密度函数平均值的相对位置及方差决定的。 可靠度可否通过计算该面积的大小给出？ 3-1应力和强度概率密度曲线的干涉面积 O r, s f(s) g(r) s0= r0 f(s) g(r) a1 a2 §4.1干涉面积法 设应力、强度两概率密度函数曲线的交点横坐标为s0=r0，并令 在应力s、强度r相互独立的情况下，零件的失效概率(不可靠度)可表示为 另一方面，零件的可靠度可表示为 即有 可见，失效概率数值上不等于干涉区的阴影面积。 由于可靠度R(t)总是小于(1-a1a2)，所以(1-a1a2)可作为零件可靠度的上限，成为衡量可靠性的一种指标，称为零件的非失效保证度。 若已知应力和强度的概率密度函数f(s)、g(r)，便可求出干涉面积a1和a2，由此便可估计出零件的可靠度。 例4-1 设某一零件的强度r和作用在该零件上的应力s均为正态分布。强度的均值和标准差分别为μr =180 Mpa， σr ＝ 8 Mpa，应力的均值和标准差分别为 μs ＝150 Mpa，σs ＝6 Mpa，试计算该零件的可靠度和非失效保证度。 解：由于应力和强度均服从正态分布，所以有 则可靠度为 现用干涉面积法估算零件的可靠度， 因s0=r0处有f(s0)=g(r0) 所以有 解得s0=r0=163.5 MPa，因此求得a1和a2分别为 可靠度的上限RU=0.9976比理论值高万分之十一，同样，可靠度下限 所以有 经验公式 该式结果与理论值相比误差约为0.14%，可见，干涉面积法得出的零件可靠度近似值 ，精度还是比较高的。 该例应力与强度均为正态分布，是为了将近似值与理论值比较，该法对其他任何形式的应力和强度的分布均适用 强度、应力和它们的干涉变量及其他许多随机事件往往需要用两个、三个或更多随机变量的函数Z=f(x1, x2, …,