构可靠指标计算的蒙特卡罗法

《构可靠指标计算的蒙特卡罗法》由会员分享，可在线阅读，更多相关《构可靠指标计算的蒙特卡罗法（45页珍藏版）》请在文档大全上搜索。

1、ch5结构可靠度计算的蒙特卡罗法n5.1蒙特卡罗法概述n5.2蒙特卡罗法的优缺点n5.3抽样模拟总数与蒙特卡罗法的精度n5.4随机变量的抽样n5.5蒙特卡罗法计算可靠指标举例n5.6蒙特卡罗的重要抽样法5.1蒙特卡罗法概述n蒙特卡罗法是随机模拟法的代名词，用这种方法计算结构可靠指标的基本思想是：对于设计阶段的结构，其功能函数及所包含变量的统计特征都是已知的，通过某种方法,根据已知的概率特性(统计特征），产生大量设计变量的样本值，将其代入功能函数，“计算”结构的状态，并对计算结果进行分析统计,直接计算其失效概率。5.1.1对蒙特卡罗法简明的理解),(21nxxxgZ ;,;,2121nnxxxx

2、xxnxxx21,nZZZ21,1.已知2.用某种方法产生样本3.计算结构的状态当然，样本的统计特征应与已知值一致。失效极限状态可靠000),(21nxxxgZ统计时以 为失效。5.1.1对蒙特卡罗法简明的理解n一般地说，随机地对每个随机变量 进行抽样，以得到一个样本值 ，然后对是否出现功能函数 进行检查，如果超过了极限状态，则认为结构或构件已经“失效”了。n 针对具体问题，首先选定所需试验次数，然后进行随机抽样并记录试验结果，最后统计超越极限状态的试验次数。iXiX0)(iXZ5.1.2蒙特卡罗法的由来蒙特卡罗法的由来n蒙特卡罗原本是摩纳哥国的一个城市，该市有赌城之蒙特卡罗原本是摩纳哥国的一

3、个城市，该市有赌城之称。赌是一种不确定性行为，具有概率的特性，本质称。赌是一种不确定性行为，具有概率的特性，本质上是一种随机实验。上是一种随机实验。n1946年冯年冯.诺伊曼等人用诺伊曼等人用电子计算机模拟了裂变物质电子计算机模拟了裂变物质的中子连锁反应的中子连锁反应，由于研究涉及秘密工作（原子弹的，由于研究涉及秘密工作（原子弹的研制），将所使用的研制），将所使用的随机模拟随机模拟方法称为方法称为蒙特卡罗法蒙特卡罗法。5.1.3利用利用随机模拟随机模拟研究结构安全问题的数学基研究结构安全问题的数学基础础n利用利用随机模拟方法随机模拟方法研究结构安全问题是一种很自然的方研究结构安全问题是一种很自

4、然的方法，因为结构建造和使用本身就是一个随机实验。法，因为结构建造和使用本身就是一个随机实验。n在结构设计阶段，由于设计变量存在着不确定性，其在结构设计阶段，由于设计变量存在着不确定性，其具体的量值是未知的，只能通过对以往实验、实测和具体的量值是未知的，只能通过对以往实验、实测和调查资料的统计分析，从概率角度来推断结构未来的调查资料的统计分析，从概率角度来推断结构未来的性状；性状；n在结构建成并使用到设计规定期后，设计中所用的变在结构建成并使用到设计规定期后，设计中所用的变量都成了规定值，结构的最终状态也完全得以确定量都成了规定值，结构的最终状态也完全得以确定（完好或失效）。（完好或失效）。n

5、所以结构从建造到使用期内的表现，就是对所设计结构所以结构从建造到使用期内的表现，就是对所设计结构的一次随机实验结果。的一次随机实验结果。5.1.4蒙特卡罗法分析蒙特卡罗法分析可靠度的步骤可靠度的步骤n从数学的角度描述为：n1.利用随机抽样以获得每一个变量的样本值： ， ， ，n2.根据上述抽样值，计算功能函数的值Z：n3.进行了N次这样的试验（抽样），则失效概率可由下式近似给出： 1X2XnX ),(21nXXXgZ NZnPf)0(显而易见，在蒙特卡罗法中，失效概率就是结构失效次数占总试验次数的比例，这就是该方法的基本出发点。 5.1.5蒙特卡罗法表示结构可靠度的表示方法蒙特卡罗法表示结构可

6、靠度的表示方法n用蒙特卡罗法表示的失效概率也可用一个示性函数表示：NiifxZINp1)(10)(00)(1)(XZXZXZIi关于(0-1)分布:如果随机变量x只能取两个值0和1，它的分布律（离散型随机变量）是：) 10( 1 , 0,)1 (1pkppkXPkk则称x服从（0-1）分布。5.2蒙特卡罗法研究结构可靠度的优缺点蒙特卡罗法研究结构可靠度的优缺点n优点：回避了结构可靠度分析中的数学困难，不需要考虑极限状态曲面的复杂性。n缺点：计算工作量大（借助于计算机）n现状：不作为一种常规的结构可靠度分析方法来使用，只是用于一些复杂情况的可靠度分析（国防、航天领域）。如何随机抽样？如何保证样本

7、与实际情况大体相符合？抽样模拟多少次？5.3抽样模拟总数与蒙特卡罗法蒙特卡罗法精度)1 (12ffppNNppzppffff)1 (96. 12/当选取95的置信度来保证蒙特卡罗法的抽样误差时，有以相对误差 来表示，有fffffpNpppp12则上式可以近似表示为： fpN224fpN2 . 0fpN/100工程结构的失效概率是一个小量，可见， N足够大时才能给出正确的估计。 估计失效概率的方差：22)(ffppEpEf对于正态分布而言，具有95置信度的96. 12/Z估计失效概率的方差22)(ffppEpEf依据（1）方差的定义： （2）示性函数的性质 )1 (12ffppNNiifxZIN

8、p1)(1经过一系列推演得到失效概率估计值的方差：5.3抽样模拟总数与蒙特卡罗法蒙特卡罗法精度n结论：精度与抽样模拟总数有关。换句话说，要想提高精度，必须将抽样模拟总数提高24fpN5.4随机变量的抽样n要求：n（1）样本足够多n（2）样本统计特性应符合随机变量已知的统计特性。n思路：n首先产生在开区间（0，1）上的均匀样本值（随机数） ，在此基础上通过一定的计算再变换成给定分布变量的随机数。5.4.1随机数的产生随机数的产生n产生随机数的方法n随机数表：将利用某种方法（高速转盘、电子装置）产生的随机数记录于磁盘中，使用时输入计算机即可（一些数学手册中还附有随机数表）。n物理方法：由物理随机数

9、发生器（安装在计算机上）将具有随机性质的物理过程变换为随机数。是真正的随机数，不会出现循环现象，但不便于对结果复查，也不便于对不同方法进行对比，且发生器的稳定性检查和维护是一项繁琐的工作，该法不常用。n数学方法：根据数论方法通过数学递推公式运算来实现。速度快，即产即用，可重复生产。会出现循环现象且随机数之间存在一定的相关性，被称谓伪随机数。数学方法产生随机数数学方法产生随机数n用数学方法产生的“随机数”，由于是按确定的算法计算出来的，所以并不是真正的随机数，但如果计算方法选择得当，它们就近似地是相互独立和均匀分布的，经得起数理统计中的独立性检验和均匀分布检验。鉴于此，人们把这种数叫作伪随机数伪

10、随机数。n用数学方法产生的“随机数”，常用的方法是同余法，包括加同余法、乘同余法乘同余法和混合同余法。1.乘同余法乘同余法 n“同余同余”的概念 n（1）举例说明余数的意义：已知2011年4月16日是星期六，问2012年4月16日是星期几？n2012年为润年，从2011年4月16日至2012年4月16日共366天，365/752.2，2012年4月16日是星期六加2天等于星期一。1.乘同余法乘同余法n“同余同余”的概念n（2）关于同余的概念n两个正整数被一个正整数除有相同的余数，谓之同余。n有同余的两个数之差，一定能够被那个固定的正整数整除。n如果a和b都是整数，而m是一个固定的正整数，则当

11、（即m能够整除(a-b)）时，则称a,b对模m同余，记作 。)(bam)(modmba 举例：65和33对模16同余（0.0625）。65-333232/162所以)16(mod3365 Ki是比值 的整数部分，即 1.乘同余法 n乘同余法产生均匀分布随机序列的递推计算公式如下：axnimxxii 01, ), 2 , 1()(modmxi1)(1mxIntKiiiiimKxx11mxrii设 、m、a是选定的正整数，a为奇数。 则相应模m 的余数为 以模除上式得： 此即第i个均匀分布的随机数，反复迭代可以产生一个随机数序列r1、r2、r3例例5.1 取 有 乘同余法：产生区间(0,1)内均匀