水利工程风险分析(全).pptVIP

作业 用蒙特卡洛方法编程，计算本章例题的可靠度。 作业形式：程序源代码电子文档，程序说明和计算结果。 交作业时间：10月31日 (7)计算 值 ： （二）重新选取β=2.5，有： 在此循环内有： 第一次 第二次 第三次 第四次 29.12 29.184 / / 51.60 51.489 / / 0.928 0.928 / / 0.372 0.374 / / β 2.5 2.5 / / 计算 值 ： （三）令 有： 在此循环内有： 第一次 第二次 第三次 第四次 23.19 22.85 22.81 / 49.76 50.44 50.53 / 0.949 0.952 0.951 / 0.314 0.306 0.305 / β 4.2 4.2 4.2 / 计算 值 ： （四）令 有： 在此循环内有： 第一次 第二次 第三次 第四次 22.60 22.57 22.57 / 50.49 50.52 50.54 / 0.953 0.953 / / 0.303 0.301 / / β 4.26 4.26 4.26 / 计算 值 ： 所以可靠度指标为β=4.26，验算点为： 总结一下： 1、计算出正态分布下的数学特征。 2、给出 和 的表达式。 3、选取β值和定义 4、进入一个小循环 （1）由 公式计算 （2）由 公式计算 （3）循环（1）~ (2)直至 的误差满足要求，并计算此时的 （4）计算 ，若 则进入下一循环。 5、再次假定一个β（若 则β应增大，若 则β应减小）。并取上一循环的 作为此次循环开始的 值，并进入循环。 （1）由 计算 。 （2）由 计算 。 （3）循环（1）~(2)，直至 满足误差要求，并计算出此时的 值。 （4）计算 ，若 则进入下一循环。 6、令 ，并取上一循环的最后的 作为本次初始 （1）由 计算 。 （2）由 计算 。 （3）循环（1）~(2)，直至 满足误差要求，并计算出此时的 值。 （4）计算 ，若 则进入下一循环。（仍为循环6）。 （5）若 ，则β= βn+1，验算点为最后得到的 JC法 JC法是拉克维茨和非斯莱等人提出来的，它适用于随机变量为任意分布下结果可靠度的求解，该法通俗易懂，计算精度又能满足工程需要，该法已经为国际安全度联合委员会（JCSS）所采用，故又称为JC法。 基本原理是：首先把随机变量Xi原来的非正态分布用正态分布替代，但对于代替的正态分布函数要求在设计验算点X*处的累积概率分布函数（CDF）值和概率密度函数（PDF）值都和原来的CDF和PDF值相同，然后根据这两个条件求得等效正态分布的均值和标准差，最后采用验算点法求解结构的可靠性指标。 JC法中正态分布替代原任意分布 下面讨论如何利用上述当量正态化来求解等效的正态分布的均值 和标准差 。 设原来验算点的累积概率为： 用来代替的正态分布验算点累积概率为： （1） 根据替代条件，则以上两个概率相等，即： 同时设原来分布验算点的概率密度函数为： 替代的正态分布验算点处概率密度函数为， 要求PDF相等，即上述两式相等，有： （2） 由（1）有： （3） 代入（2）有： 即： （4） 由（3）可得： （5） 其中 和 分别代表变量Xi原来的累积概率分布函数和概率密度函数，Φ和φ分别代表标准正态分布下的累积概率分布函数和概率密度函数，其值分别可由查表和下式求得： （6） （4） （5） （6） 该三式即为对一般分布的计算公式。 当分布为正态分布时，不必用上述公式进行转换，而可以直接将该变量的均值和标准差作为变换后的均值和标准差。 当分布为对数正态分布时，则可对上式进行简化。 由前可知： （7） （8） 又： （9） 式中： （11） 依概率论有： （10） （9）和（11）分别为对应正态分布下的Xi的累积概率分布和概率密度函数，代入（4）和（5）得到对数分布下分布变量Xi的替代正态分布变量的均值和标准差。 化简后得： （12） （13） 总结一下： （12） （13） （10） 求得替代的正态分布的均值 和标准差 后，即可用与一次二阶矩法大致相同的解法进行求解。 下面给出用JC法求解可靠