第七讲蒙特卡洛模拟应用.ppt

第七讲 蒙特卡洛模拟的应用 一、蒙特卡罗方法在积分计算中的应用 计算多重积分是蒙特卡罗方法的重要应用领域之一 如欲求s重积分 P＝P(x1，x2，…，xs)为s维向量空间点 Vs为积分区域，P在Vs内 若S很大，维数很高？？ 定义取Vs上任意P点的概率密度函数为f(P) 令 则被积函数化为： 求积转化为求函数g在概率分布f(p)下的数学期望值。 给定f(p)分布随机抽样 抽取随机向量 P 的 N 个样本：Pi，i＝1，2，…，N则 就是积分的近似估计，其结果的方差为： 一般情况下，概率密度函数取为积分区域的均匀分布密度函数。 为所围的n维体积 二、解非线性方程 MC方法也是求解非线性代数方程组的有效方法。 用MC如何求？？ 构造目标函数 搜索X使得目标函数接近0 解法步骤 a.选出发点。在自变量的取值范围内[a,b]内选一组随机抽样值X0为出发点。 b.选初始搜索参数。初始步长 ，补偿改变参数 ，初始增量 ，以及搜索失败控制次数M等。 c.进行第k+1步的搜索。此时已知前k-1步搜索中目标函数的最小值且令S=minQ和搜索失败次数m。 d.确定搜索步长 。 如果 ，说明第k步搜索成功，令 ，取搜索步长 。 否则，第K步搜索失败，m增加1，若m<M则其他搜索参数不变；若m>M,则减小搜索步长，取 。 e.确定搜索增量 取一组[0,1]上的随机数{r1, r2,… rn}得到正负随机数R={2r1-1, 2r2-1,…2rn-1} 则k+1步的随机搜索点可定为 如果Q(Xk+1)小于一个小的整数，得到方程组的解。 否则，继续下一步搜索。 三、统计物理中的应用 微正则系综NVE 正则系综NVT 等温等压系综NPT 巨正则系综u（化学势）VT 满足系综条件的粒子（原子）分布是很多的，每一种存在的概率密度为其配分数。因此，我们在某一系综内求材料的一个物理量（密度）的时候，可采用Metropolis抽样法，抽取满足原子n个分布抽样，MC模拟后，得到物理量的期望值。 四、材料的晶粒生长 有限元冷轧 冷轧中通过有限元计算得到 五、复合材料强度的MC模拟 纤维增强金属基复合材料 纤维 基体（金属 或树脂） 复合材料的破坏行为 界面脱开 纤维断裂 基体屈服 纤维断裂引起应力集中，进一步导致另外两种破坏。 纤维的统计强度分析 在室温下，纤维的强度往往具有较大的分散性，这主要是由于生产工艺、环境因素等方面的影响，使得纤维的表面和内部存在各种不同的缺陷。 Coleman在1958年证实了单根纤维的强度应该服从weibull分布.这一结论得到了实验的支持，目前已经在实际工程中得到了广泛的应用。 纤维强度的威布尔分布 纤维的长度 材料参数 长度为L的纤维其强度在0- 之间的概率 分布函数 MC时纤维强度的抽样 在[0,1]之间取均匀随机数 得到长度L的纤维的强度抽样 直接抽样法 复合材料应力场的解析解法 加权系数由裂纹断面力自由边界条件求出。 纤维断裂 复合材料应力场的解析解法 纤维断裂引起了基体屈服。 复合材料应力场的解析解法 纤维断裂引起了界面脱开。 MC模拟过程