三维晶粒生长计算机模拟教案分析.pptVIP

Monte Carlo法模拟三维晶粒生长 小组介绍 华中科技大学电子科学与技术系006班二组 成员名单 颜庙青 李进 李怡然 赵伟民 朱俊锋 刘栋 陈峙 游州 余洋 苏丽琼 周丽 引言 在研究材料晶粒组织,亚结构等微观组织尺度的材料力学问题, 首先要建立细观结构 。材料的显微组织模型.材料模型要反应细观结构的几何形状,分布和界面条件.通常采用代表性体积元法,即设计初一定的体积元,来表示具有规则的,周期性分布的 典型体积元法的缺点是无法研究颗粒空间分布对材料性能的影响;对于多晶体的材料,也不能表示具有不同尺寸,形状等参数的晶粒组织. 本文的目的 一 基本物理模型与处理方法 1．1 基本理论 通常，对多晶材料，在晶核形成之后，其晶粒生长现象分为两种不同类型：正常晶粒生长和异常晶粒生长（也称重结晶）。正常晶粒生长时，晶粒尺寸一致增长，归一化晶粒尺寸F(R/)和拓扑分布函数P()不随时间改变。其中，R是晶粒半径，是平均晶粒半径，是晶粒边数。这种情况下，生长有如下规律：R=K获 其中，t为时间，K,B为常数，n=0.5.当时间很长，即t较大时，若R=Ro (Ro为初始的晶粒半径)，式（1），（2）是等效的且m=1/n. 二 多晶体组织三维建模及可视化 2.1 三维单相多晶体晶粒组织建模 2.2 晶粒组织模型的三维可视化 2.1 三维单相多晶体晶粒组织建模 不论单相材料还是多相复合材料,其基体组织一般是多晶聚集体,而晶粒的形状,尺寸及晶体取向等参数对其性能有重要的影响,因而,在晶粒尺度对多晶体组织进行表征和模型化是非常重要的.材料的显微组织模型化不仅要符合材料显微组织的几何形态要求,还要尽可能符合材料显微组织形成和演变的物理过程,已经有工作利用Monte Carlo方法(Q态Potts模型)及其改经的方法来模拟三维多晶体组织的晶粒长大过程,研究表明Mont Carlo方法可以成功魔力曲率驱动的晶粒长大过程,仿真结果符合世界材料中多晶体组织形成和演化过程的基本规律. 采用具有周期性边界条件,边长为32的简单立方点阵,系统的点阵位置总数N=32*32*32=32768,点阵中的每个阵点被赋予一个初始取向随机数 S( S E[1,32] ) ,代表相应体积元的晶体取向,具有相同取向数的26个相邻点阵为之属于相同的晶粒,而不同取向数的相邻点阵位置之间形成晶界,每个阵点的能量为: 其中,J为单位面积上的晶界能;S,S为第i个点位置和其邻居j的取向数;N是计算中考虑的近邻数目(N=26,包括第1,第2和第3近邻); 为Kronecker 函数. 在模拟过程中,随机选取一个点阵位置S,随机赋予它任意一个近邻的取向数S”,从S向S”的转变按照如下概率进行: 其中能量差为：E= 通常,每N个单元完成再取向尝试后计为一个Monte Carlo Step(简写为MCS),定义为模拟的时间单位,N是系统的点阵位置总数,通过改变总的MCS步数,可以获得不同平均晶粒尺寸的多晶体组织.通过上述方法获得的多晶体组织模型包含的细观力学结构参量有:晶粒的尺寸和尺寸分布,晶粒的形状,晶粒的取向,晶粒的各向异性和晶界的材料特性等,可以应用于材料特性等,可以应用于材料的塑性变形,损伤以及断裂过程等的细观力学研究. 2.2 晶粒组织模型的三维可视化 单晶多晶体组织三维建模的计算机实现,考虑了这种显微组织模型三维可视化的需要,用1个三维数组来存储多晶体组织模型数据结构,数组的元素可以根据需要包含多晶体组织体积元的空间坐标,晶体取向,形变组织的存储能等信息.从科学三维可视化的角度看,整个三维数组就是1个体数据结构.我们可以根据不同的需要,来实现多晶体组织体数据的三维可视化 三 程序仿真过程 3．1 模拟晶粒生长过程如下： a.随机选取一个试验格点； b.该格点的一个新的试验取向在其余Q-1种可能的取向中任选其一； c.计算可能改变取向有关的能量； d.通过计算改变几率W,判断取向改变是否可以实现 式中，为Boltzman常数，T为转变温度。 即若 =0,则认为改变取向式允许的，若=0,则取向概率为，则改变取向允许，否则不允许，即只有体系能量减少的转变才是允许的。 经试验，发现这种基于Monte Carlo 的模拟方法计算量相当大，需要在大，巨型计算机上才能得到较好的结果。基于上述原因，此仿真方法中仿真的是32的简单立方体。 3.2 程序说明 语言：FORTRAN 描述： 程序从一个名叫Grain