[交通运输]5_机械可靠性的蒙特卡罗方法及应用_2.pptVIP

机械可靠度蒙特卡罗仿真及应用 张建国 工程系统工程系 内容简介 蒙特卡罗方法介绍 蒙特卡罗方法基本步骤 直接抽样方法及实例 抽样方差减缩技术 对偶抽样 条件期望抽样 重要抽样 重要抽样方法及工程实例 蒙特卡罗方法介绍 蒙特卡罗法(Monte Carlo)是一种以概率统计理论和方法为基础的数值计算方法，是最基本，适用范围最广的方法 随机模拟(Random simulation ) 随机抽样技术(Random sampling) 统计实验（Statistical testing) 计算机应用技术的发展与进步，以Monte Carlo法为基础的各类抽样方差减缩技术，成为解决非线性机械、结构系统可靠性分析的重要手段 对偶抽样 条件期望抽样 重要抽样 蒙特卡罗方法介绍 蒙特卡罗方法基本思想 当所求的问题是某个事件的概率时，可以通过抽样试验的方法得到这种事件出现的频率，把它作为问题的解 用基本随机变量的联合概率密度函数进行抽样，用落入失效域内样本点的个数与总样本点的个数之比作为失效概率的无偏估计值 由大数定理可知，当模拟次数足够大时， 依概率1收敛于 。 蒙特卡罗方法特点 Monte Carlo法优点 方法及其程序结构简单，容易实现 对功能函数的形式、随机变量的维数及其分布型式均无特殊要求 收敛的概率和收敛的速度与问题的维数无关 只要抽样次数足够多，可靠度精度满足要求，可检验其他方法结果 计算精度高，数值模拟的误差、次数和精度容易确定 通用性强 缺点 效率太低，对于可靠度要求普遍比较高的机械、结构可靠性问题，工程实用性较差 蒙特卡罗方法基本步骤 构造概型(基本变量的概率分布模型) 概型的概率分布、分布的分位点、数学期望或事件出现的概率即为所求的解 定义随机变量(定义随机变量以及功能函数Z) 随机变量的分布或数字特征恰好就是问题的解，根据概率统计模型的特点和计算实践的需要，尽量改进模型，以便减小方差和降低费用，提高计算效率 通过模拟获得子样(计算功能函数抽样(随机数生成)) 根据概型找出对随机变量的抽样方法，实现从已知概率分布抽样 统计计算 直接抽样法 在整个区间上进行均匀抽样，通过大量模拟抽样，统计失效次数 由大数定理得出，在大量随机模拟的情况下，失效概率近似等于其失效发生的频率 进行可靠性分析的一种最基本的方法 直接抽样方法几乎不需要做任何前期准备工作和特殊处理 直接抽样法的效率和精度也较低，适用于结构可靠度不是很高的情况 直接抽样法 直接抽样方法 随机数的产生和检验 乘同余法 加同余法 混合同余法 正态随机变量的抽样方法 哈斯汀（Hasting）有理逼近法 统计近似法 变换法 一般随机变量的抽样方法 反函数方法 舍选法 直接抽样法流程 定义随机变量概率分布类型以及定义功能函数 产生随机数 用抽样方法对随机变量抽样,代入随机数计算 带入极限状态方程，计算 循环叠代计算,统计 直接抽样方法算例 例1：结构的功能函数为 服从正态分布 服从极值I型分布 试用直接抽样估算结构的失效概率。 直接抽样方法算例 思路 例题已经将随机变量的分布以及功能函数的定义给出,所以可以直接利用直接抽样方法对 和 进行抽样. 随机数自动生成,带入随机变量进行抽样计算,再将抽样值带入功能函数 后，判断 的符号，对于模拟次数 ，统计 的次数 根据公式 估算结构的失效概率。 在下表中给出了 次时的模拟结果 直接抽样方法算例 失效概率的估算值： 算法比较 用JC方法计算的结果为 ，相应的可靠度指标为 ； 用实用分析法计算的结果为 ， 蒙特卡罗方法模拟5000次的结果 ， 实用分析法的结果比JC法更为准确 直接抽样方法算例 例2：某一串联结构体系有三个主要失效模式，相应的功能函数为 其中 为相互独立的标准正态随机变量。试用直接模特卡罗方法估算结构体系的失效概率。 直接抽样方法算例 思路：本例为计算结构体系可靠度问题，结构的功能函数可以表示为 下表给出模拟次数100～1000000次时的结构体系失效概率的估算结果。 抽样方差减缩技术 直接抽样方法中 这就说明 与 成反比;当 为一小量, 估算精度 即 才能满足对失效率的抽样估计。 而工程结构破坏的概