元 胞 自 动 机

一、元胞自动机的定义元胞自动机，即CellularAutomaton(CA)，也称为细胞自动机、点格自动机、分子自动机或单元自动机。它是一种利用简单编码与仿细胞繁殖机制的非数值算法的空间分析模式。散布在规则格网(LatticeGrid)中的每一元胞(Cell)取有限的离散状态，遵循同样的作用规则，依据确定的局部规则作同步更新。大量元胞通过简单的相互作用而构成动态系统的演化。不同于一般的动力学模型，元胞自动机不是由严格定义的物理方程或函数确定，而是用一系列模型构造的规则构成。凡是满足这些规则的模型都可以算作是元胞自动机模型。因此，元胞自动机是一类模型的总称，或者说是一个方法框架。其特点是时间、空间、状态都离散，每个变量只取有限多个状态，且其状态改变的规则在时间和空间上都是局部的。二、元胞自动机的分类1.按动力学行为平稳型元胞自动机周期型元胞自动机混沌型元胞自动机复杂型元胞自动机2.按元胞空间分类一维元胞自动机二维元胞自动机三维元胞自动机多维元胞自动机三、元胞自动机的应用•1.对凝固组织的模拟（枝晶的生长、定向凝固组织，焊接接头微观组织）•2.对热压缩行为的研究(位错组态等)•3.对再结晶组织的模拟（静态再结晶、动态再结晶）元胞自动机自产生以来，被广泛地应用到社会、经济、军事和科学研究的各个领域。应用领域涉及社会学、生物学、生态学、信息科学、计算机科学、数学、物理学、化学、地理、歹境、军事学等。四、元胞自动机模型•1.元胞和元胞的划分•2.元胞的状态空间•3.邻居类型和边界条件•4.形核规则及核心生长规则•5.模型算法划分元胞优点缺点三角形邻居数目较少不易表达和显示，需要转换为四边形正方形易于表达和显示不能较好的模拟各向同性的现象六边形可较好的模拟各向同性，使模拟结果显示更接近于真实情况表达与显示困难表1-二维元胞划分类型及其优缺点二维邻居类型（1）Neumann型（4个近邻元胞）（2）Moore型（8个近邻元胞）（3）交替Moore型边界条件(1)周期边界条件：指相对边界连接起来的元胞空间.对于一维空间，表现为一个首尾相连的圈，对于二维空间，上下相连，左右相接，形成一个拓扑圆环面。(2)反射边界条件：指在边界外邻居的元胞状态是以边界为轴的镜面反射。(3)固定边界条件：指所有边界外元胞均取某一固定值，如0，1等。五、模型的建立过程•1.选择形核方式：位置过饱和(SS)•2.元胞的划分：正方形，200×200•3.边界条件：周期性•4.邻居类型：交替Moore型•5.设定元胞的状态：0—未再结晶，1—再结晶•6.初始晶粒组织：（1）颜色变量：101~150的正整数（2）晶粒取向：1~180的正整数，随机表示（3）晶界变量：存放晶界元胞•7.动态再结晶过程：（1）初始条件：热变形条件和相关材料参数，初始位错密度等obDTm(K)μ(Gpa)b(nm)Qb(KJ/mol)()Qa(KJ/mol)θm(°)3/ms2(/)mJm（2）位错密度其中，：硬化系数，：软化系数（母相晶粒中每个元胞的初始位错密度相同，新再结晶晶粒中每个元胞的初始位错密度为零。）临界位错密度：式中，：界面能；m：晶界迁移率；：单位长度位错线的能量，：位错平均自由程12()ddkkdtdt1k2k1/3220[]3icblmi22cbl（3）形核：形核率：其中，C：常数；T：绝对温度；R：气体常数；：激活能（4）晶粒长大：长大速率：晶界迁移率：（—材料的晶界厚度；—绝对零度时的晶界扩散系数；—晶界扩散激活能）(,)exp[]aQNTCRTaQiivmfexp()obbbDQmkTRTobDbQ驱动力：—第i个再结晶晶粒的位错密度；—与之相邻晶粒的位错密度；—半径；—界面能；（—再结晶晶粒与相邻晶粒间的取相差；—大角度晶界的取向差；—大角度晶界的晶界能）()2/imiiifrimiri[1ln]iiimmmimm（5）时间步长：（—晶界迁移速率修正系数，取值1~10）再结晶时间：再结晶晶粒尺寸：20202max1GBdkdtvkmkGBk433RtNv4032222RtRvdRvdtvtN选择形核方式从界面输入初始条件：网格规模、边界条件、邻居类型、形核率等时间步t=1以选定的形核方式形核核心或晶粒生长确定再结晶体积分数Φ，元胞状态，晶粒取向等输出数据和图像t=t+1直到再结晶体积分数Φ=100%输出晶界图输出再结晶曲线和Avrami曲线计算晶粒尺寸分布表2-再结晶CA模拟的N-S流程图输入初始状态定义第二相粒子分布判断晶界，位错密度演变DRX形核晶粒长大更新元胞状态应变叠加，变形完成结束是否是否否*包含粒子的动态再结晶CA程序流程图精品课件！精品课件！目前的研究情况1.在定性模拟不同应变速率下动态再结晶晶粒的形核即长大过程的基础上，考虑温度对DRX过程的影响，改进位错密度演化机制、形核机制和再结晶晶粒长大机制，建立模型。2.以动态再结晶的金属学原理为基础，综合考虑加工硬化、动态回复、再结晶形核率、位错密度、以及变形温度、应变速率的影响。3.采用晶体塑形有限元（CPFEM）和元胞自动机（CA）耦合的方法模拟热变形的静态再结晶；通过增加形变规则，得到不同变形程度时的形变微观组织，在此基础上建立模型，观察微观组织演变规律及其动力学行为，并分析了储存能的变化规律。4.将宏观热变形参数与位错密度变化相耦合，只能模拟低应变速率下的动态再结晶过程。若采用真实时间替代模拟时间步，可模拟个应变速率下动态再结晶过程的微观组织结构变化。