微观至介观尺度的模拟方法概述.pptVIP

  1. 1、本文档共65页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
6.4 元胞自动机的分类 一维细胞自动机在随机的初始条件下所产生的花纹可以归结为4类: 1、固定值型:细胞自动机演化到一定时刻就变成了一种颜色的方格; 2、周期型:细胞自动机固定在具有一定循环结构中不再改变; 3、混沌型(或叫随机类型):结构在不停的变化,但是它们没有确定的变化规律; 4、复杂型:这类结构介于完全秩序与完全混沌之间,会产生一些局部的复杂结构,但整体似乎又不是完全混沌随机。 6.5 概率性元胞自动机方法 将确定性元胞自动机变为非确定性的基本方法有两种: 第一种方法就是随机地选择所研究的晶格格点,而不是系统化地按顺序选择,但是要使用确定性变换规则; 第二种方法就是用概率性变换代替确定性变换,但要系统地研究所有格点。 主要讨论第二种方法,并将之归为概率性或随机性元胞自动机。 6.5 概率性元胞自动机方法 概率性元胞自动机,就其基本过程和要素方面而言,非常相似于普通的元胞自动机,只不过转变规则由确定性的换成了随机性的。 设有N个格点组成一个一维链,其中每个格点有k个可能的状态Sv=0, 1, 2, ..., k-1。从而整个链共有kN个不同的排列方式。由(S1, S2, …, SN)描述的某给定晶格状态用下式整数标记: (6.3) 6.5 概率性元胞自动机方法 假设每个状态i的存在概率为Pi。作为时间的函数的Pi(t),按照其转变概率以离散时间步t=0, 1, 2, …的方式变换发展。如果只考虑邻近的时间步(t-1),这一规则可用下式给出: 转移概率Tij就表示由前一时刻的j状态转变到i状态的概率。 对于离散型元胞自动机方法,转移矩阵Tij 是由局部规则决定: (6.4) (6.5) 6.5 概率性元胞自动机方法 和 分别表示状态j和i的格点变量;因而,变量的转换只有其最近邻及其自己的状态有关。 虽然概率性元胞自动机与Metropolis蒙特卡洛算法之间具有一定的相似性,但二者之间还是有差别的。这种差别主要表现在两个方面: 第一,蒙特卡洛方法每个时间步只更新一个格点,而概率元胞自动机像大多数自动机一样,每次要全部一起更新; 第二,元胞自动机没有本征的长度或时间标度。 元胞自动机的标定参数主要是由构成物理模型的基础来决定,而不是由所采用的元胞自动机算法来决定。 6.5 概率性元胞自动机方法 元胞自动机与波茨蒙特卡罗自旋模型的区别: (1)元胞自动机对微观体系不存在内禀标度,如果选择合适的基础单元,并且建立与场变量相匹配的代数、微分或积分方程,那么元胞自动机可以用于任意空间和时间尺度上的问题处理。蒙特卡罗方法对于微观体系是有内禀标度的。 (2)在蒙特卡罗方法中,广义自旋格座是随机抽样顺序考察的,而元胞自动机则是同步一起更新。 (3)元胞自动机比多态波茨模型使用了更多的确定性或概率性变换规则。 6.6 非平衡现象的模拟 6.6.1热力学模拟 在金属的热变形过程中,会发生诸如再结晶、连续与非连续型晶粒生长和不连续沉淀等非平衡转变现象和微结构瞬态问题。按照微结构的观点,这些转变现象都是由于高角晶界的运动引起的。由于吉布斯自由焓存在梯度,原子或原子团将从一个晶粒越迁转移到其邻近晶粒。对同相界面,其净驱动压强 (6.6) 6.6.1热力学模拟 在实际材料中,各种贡献都将影响到局域自由焓的值。 (1)在冷加工金属中,位错密度ρ的增加对所储存的弹性能的贡献,在对驱动压强中占最大的份额。 Δρ是界面两边的位错密度差;μ为各向同性极限下的体剪切模量;b表示伯格矢量的大小。 如果将存在于元胞壁的位错(ρw)和元胞内的位错(ρi)的贡献分别表述。 (6.7) (6.8) ρw只能用亚晶粒尺寸D和亚晶粒壁的界面能γsub表述 6.6.1热力学模拟 (2)作用于各晶粒上的拉普拉斯压强或毛细压强的贡献。 对于常见的晶粒粒度分布和球形晶粒 α为2~3的常数;γ为界面能;1/R为曲率。 对于薄膜,还有来自表面能梯度的贡献: 式中,B表示薄膜宽度;h为膜厚;Δγ代表表面能变化量。 (6.9) (6.10) 6.6.1热力学模拟 (3)在过饱和态,对驱动压强还有一项化学贡献。其对应的转变称为非连续沉淀。 kB为波耳兹曼常数; Ω为原子体积; T1为(数值)实验中的实际温度; T0为相应于T1时过饱和浓度的平衡温度;c0为浓度。 (6.11) 6.6.1热力学模拟 (4)其他因素对总压强的贡献 在总的驱动压强中,还要考虑冷加工或硬化金属间化合物中由于损失长程有序而产生的贡献。更进一步还应当考虑来自于磁性、弹性及温度场等梯度的贡献,但是这类贡献在实际应用中意义不大。 (5)可能的反驱动力 杂质阻力以及在有序化合金中高角晶界运动在远处产生畴的结构。 6.6.2 动力学模拟 为使原级再结晶能够启动,要在热力学、力学和

文档评论(0)

181****4529 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档