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

(a) 根据提出的问题确定各变量之间的确定性函数关系。 (b) 根据提出的问题构造一个简单、适用的概率模型或 随机模型，使问题的解对应于该模型中随机变量的某 些特征 ( 例如概率、均值和方差等 ) 。 (c) 根据模型中各个随机变量的分布，在计算机上产生 随机数，实现一次模拟过程所需的足够数量的随机数。 通常先产生均匀 分布的随机数，然后生成服从某一 分布的随机数。 二、蒙特卡洛模拟求解步骤 : (d) 根据概率模型的特点和随机变量的分布特性，设计 和选取合适的抽样方法，并对每个随机变量进行随 机抽样。这里的抽样方法有直接抽样、分层抽样、 相关抽样、重要抽样等。 (e) 按所建立的模型进行仿真计算，求出问题的一个随 机解。 (f) 统计分折模拟试验结果，给出问题的概率解以及解 的精度估计。 在可靠性分析和设计中，用蒙特卡洛方法可以确定 复杂随机 变量的概率分布和数字特征；可以通过随 机模拟估计系统和零件 的可靠度；也可以模拟随机 过程、寻求系统最优参数等。 二、蒙特卡洛模拟求解步骤 : 分布名称 密度函数 f(x) 或 f(t) 或 [0 ， 1] 均匀分布 1 均匀分布 指数分布 标准正态分布 正态分布 是标准正态分布抽样 对数正态分布 ( ) F X ? ( ) F t ? ? 1 ( ) b a ? ( ) b a a ? ? ? t e ? ? ? 1 1 ln(1 ) ln ? ? ? ? ? ? ? 或 2 2 1 2 t e ? ? 1 2 1 2 2ln cos 2 2ln sin 2 ? ?? ? ?? ? ? 2 2 ( ) 2 1 2 t e ? ? ?? ? ? 01 N t ? ? ? ? 01 N t 2 2 1 exp 2 ( ) [ln( ) ] , ( ) 2 t a t a t a ?? ? ? ? ? ? ? ? ? ? ? ? ? ? 01 exp( ) N a t ? ? ? 常见分布函数随机变量的随机抽样公式 例 4-5 某铝合金板的形状如图所示。受弯矩作用，其 尺寸 H 、 h 、 ρ 均服从正态分布，分布参数为： 2 2 2 ~ (40,0.2 ), ~ (20,0.1 ), ~ (4,0.167 ) H N h N N ? 试确定理论应力集中系数 α σ 的分布类型及分 布参数。 H ρ h b M M 受弯矩作用的铝合金板 三、蒙特卡洛模拟算例 蒙 特 卡 洛 模 拟 算 例 的 程 序 框 图 开始 输入 H 、 h 、 ρ 分布的类型及参数 ; N j =1 分别从 H 、 h 和 ρ 的分布中产生随机数 H f , h f , ρ f 计算 α σ 的随机数 α σj j = N 进行分布类型判断、估计分布参数 输出 α σ 的分布类型和分布参数 结束 j = j +1 解：理论应立集中系数 α σ 的计算公式为 0.85 1 1 9.6 1.12 1 H h h H h ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 蒙特卡洛模拟算例 由于 H 、 h 、 ρ 均服从正态分布，所以根据正态分布 的抽样公式以及 α σ 的计算公式编制计算机程序，上 机运行。输入参数为： 模拟次数 N =1000 。 输出结果为： α σ 服从正态分布， 1.4795 ? ? 0.0183 S ? ? 2 (1.4795, 0.0183 ) N ? ? ? 均值为： 标准差： 即： 2 2 2 ~ (40,0.2 ), ~ (20,0.1 ), ~ (4,0.167 ) H N h N N ? 蒙特卡洛模拟算例 机械 可靠性设计 § 4.3 一次二阶矩方法求可靠度 _ 工程方法 第四章 机械可靠性设计分析方法 § 4.1 干涉面积法 § 4.2 分布代数 § 4.4 蒙特卡洛模拟方法 § 4.5 变异系数传递规律 从可靠度计算的普遍方程可以看出，对于应力和强度比较复 杂的分布，由于积分困难，往往难以得出问题的解析解。 因此如何采用一些较好的近似方法，能 比较方便地求得满足工程 精度要求的零件可靠度的近似解 ，一直是人们探讨的一个问题。 3-1 应力和强度概率密度曲线的干涉面积 O r , s f ( s ) g ( r ) s 0 = r 0 f ( s ) g ( r ) a 1 a 2 应力和强度两个概率密度函 数的交叉区，即干涉区阴影 面积的大小，反映了零件或 结构可靠度的高低。 该面积越小，可靠度越高， 反之，可靠度越低。 干涉面积的大小是由两个概率密度函数平均值的相对位置及方 差决定的。 可靠度可否通过计算该面积的大小给出？ 3-1 应力