[机械可靠性设计0704分析方法.pptVIP

机械可靠性设计 南京航空航天大学机电学院 于明礼 (a)根据提出的问题确定各变量之间的确定性函数关系。 (b)根据提出的问题构造一个简单、适用的概率模型或随机模型，使问题的解对应于该模型中随机变量的某些特征(例如概率、均值和方差等) 。 (c)根据模型中各个随机变量的分布，在计算机上产生随机数，实现一次模拟过程所需的足够数量的随机数。通常先产生均匀 分布的随机数，然后生成服从某一分布的随机数。 二、蒙特卡洛模拟求解步骤: Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. (d) 根据概率模型的特点和随机变量的分布特性，设计和选取合适的抽样方法，并对每个随机变量进行随机抽样。这里的抽样方法有直接抽样、分层抽样、相关抽样、重要抽样等。 (e) 按所建立的模型进行仿真计算，求出问题的一个随机解。 (f) 统计分折模拟试验结果，给出问题的概率解以及解的精度估计。 在可靠性分析和设计中，用蒙特卡洛方法可以确定复杂随机 变量的概率分布和数字特征；可以通过随机模拟估计系统和零件 的可靠度；也可以模拟随机过程、寻求系统最优参数等。 二、蒙特卡洛模拟求解步骤: Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 分布名称 [0，1]均匀分布 均匀分布 指数分布 标准正态分布 正态分布 对数正态分布 密度函数f(x)或f(t) 1 或 是标准正态分布抽样 常见分布函数随机变量的随机抽样公式 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例4-5某铝合金板的形状如图所示。受弯矩作用，其尺寸H、h、ρ均服从正态分布，分布参数为： 试确定理论应力集中系数ασ 的分布类型及分布参数。 M M 受弯矩作用的铝合金板 三、蒙特卡洛模拟算例 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 蒙特卡洛模拟算例的程序框图 开始 输入H、h、 ρ分布的类型及参数;N j=1 分别从H、h和ρ的分布中产生随机数Hf , hf , ρf 计算ασ 的随机数ασj j=N 进行分布类型判断、估计分布参数 输出ασ 的分布类型和分布参数 结束 j=j+1 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 解：理论应立集中系数ασ的计算公式为 蒙特卡洛模拟算例 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 由于H、h、ρ均服从正态分布，所以根据正态分布 的抽样公式以及ασ 的计算公式编制计算机程序，上 机运行。输入参数为： 模拟次数N=1000。 输出结果为： ασ 服从正态分布， 均值为： 标准差： 即： 蒙特卡洛模拟算例 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 用蒙特卡洛仿真计算应力和强度为任意分布时的可靠度 任意分布的应力—强度模型都可以用蒙特卡洛模拟法求可靠度的近似值，结果的精度随模拟的次数的增多而增高。模拟程序的流程图如右图所示。 开始 输入应力和强度分布类型和 参数，模拟次数N,置j=1 对应力和强度各产生 一个随机数xsj和xSj 比较xsj和xSj并记下xsjxSj的次数N1 j=N ? 输出R=N1/N 结束 j=j+1 Evaluation only. Created w