可靠性和可靠性灵敏度分析的MonteCarlo数值模拟法.ppt

《可靠性和可靠性灵敏度分析的MonteCarlo数值模拟法.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、机械系统可靠性和可靠性灵敏度分析的MonteCarlo数值模拟法授课教师：黄贤振图1机械系统的分类机械系统是指由多个机械基本要素（如机械零、部件，机构，机器等）组成的，可以完成所需动作(或动作过程)，实现动作的传递和力的变换以及机械能的转化和利用的装置。1机械系统2机械可靠性应力分析应力分布参数强度分析强度分布参数载荷分布参数机械系统构件分布参数机械系统性能分布参数机械系统构件分布参数f(x)fR(xR)fS(xS)x机械可靠性模型应力是指对产品功能有影响的各种因素，比如机械应力、变形、磨损等。强度是指产品承受应力的能力，比如机械强度、刚度、抗裂度等。假定机械系统的强度随机变量xR，应力随

2、机变量为xS，机械系统极限状态函数（功能函数）为极限状态的定义为：整个机械系统的一部分超过某一特定状态就不能满足设计规定的某一功能要求，此特定状态为该功能的极限状态。2机械可靠性f(x)fR(xR)fS(xS)xxSdxS应力xS落在附近dxS小区间内的概率为在应力xS落在小区间dxS内的条件下，强度xR小于应力xS的概率为应力xS落在小区间dxS内，同时强度xR小于应力xS的概率为根据可靠度的定义，对于应力xS所有的可能值强度xR均小于应力xS的概率，及事件(xR

3、变量的机械系统的功能函数可以表示为则极限状态方程g(x1,x2,…,xn)=0将结构的基本随机变量空间分为失效域和可靠区域两部分。2机械可靠性式中,f(xi)(i=1,2,…,n)为随机变量xi的概率密度函数。式中，f(x1,x2,…,xn)是基本随机变量x=[x1,x2,…,xn]T的联合概率密度函数。若基本随机变量是相互独立的，则有失效概率Pf可表示为MonteCarlo可靠性分析方法又称随机抽样法、概率模拟法或统计试验法。该方法是通过随机模拟或者说统计试验来进行结构可靠性分析的。由于它是以概率和数理统计理论为基础的，故被无理学家以赌城MonteCarlo来命名。3MonteCarlo

4、可靠性分析3MonteCarlo可靠性分析根据概率论的大数定律，若有来自同一母体且有相同分布的n个相互独立的随机样本x1,x2,…,xn,他们具有相同的有限均值μ和方差σ2，则对于任意的ε>0,有上式表明样本均值是依概率收敛于母体的均值μ的。3.1MonteCarlo模拟的理论基础3MonteCarlo可靠性分析另外，设随机事件A发生的概率P(A),在n次独立试验中，事件A发生的频率为m，则随机事件A发生的频率W(A)=m/n,对于任意ε>0，有上式表明事件发生的频率是依概率收敛于事件发生的概率的。3MonteCarlo可靠性分析3.2MonteCarlo可靠性分析的原理与计算公式1.Mo

5、nteCarlo方法求解失效概率估计值的计算公式失效概率可以改写成下式所示的失效域指示函数IF(x)的数学期望形式：式中，为失效域指示函数；Rn为n维变量空间；E[.]为数学期望算子。3MonteCarlo可靠性分析失效概率为失效域指示函数的数学期望，依据大数定律，失效域指示函数的数学期望可以有失效域指示函数的样本均值来近似。即随机变量的联合概率密度函数fX(x)抽取N个样本xj(j=1,2,…,N),落入失效域F内样本点的个数Nf与总样本点的个数N之比即为失效概率的估计值。(1)3MonteCarlo可靠性分析2.MonteCarlo失效概率估计值的方差分析由式(1)可以看出，失效概率估

6、计值的随机样本xj(j=1,2,…,N)的函数，因此也是一个随机变量。为了对计算出的的收敛性有一个清楚的认识，有必要对的方差进行分析。对式(1)两边求数学期望，可得失效概率估计值的期望如下所示：样本xj与母体x独立同分布3MonteCarlo可靠性分析由上式可知，，即为Pf的无偏估计在数值模拟的过程中，以指示函数IF(x)的样本均值近似代替E[IF(x)],则失效概率估计值的期望可以近似表达为3MonteCarlo可靠性分析失效概率估计值的方差可以通过对式(1)两边求方差得如下：由于样本方差依概率收敛于母体的方差，所以可以用IF(.)的样本方差(2)xj与母体x独立同分布3MonteCar

7、lo可靠性分析将上式代人式(2)，可得失效概率估计值的方差估计为进而得到估计值的变异系数为指示函数IF(x)方差的无偏估计可以进一步表达为？3MonteCarlo可靠性分析IF(xj)方差的无偏估计xj相互独立xj与母体x独立同分布依据随机变量的分布形式和参数，由随机样本的产生方法产生N组随机样本向量的样本xj=(xj1,xj2,…,xjn)(j=1,2,…,N)。将随机向量样本xj代人极限状态方程，并根据状态指示函数I