采用Cellularautomaton法模拟动态再结晶过程的研究-机械工程学报.PDF

第4 1 卷第2 期 机 械 工 程 学 报 Vol. 41 No. 2 2005 年 2 月 CHINESE JOURNAL OF MECHANICAL ENGINEERING Feb. 2 0 0 5 采用Cellular automaton 法模拟 动态再结晶过程的研究* 肖 宏 (燕山大学机械工程学院 秦皇岛 066004) 柳本润 (东京大学生产技术研究所 东京 153-8505) 摘要：结合金属塑性成形过程的冶金学原理，给出了一种新的模拟动态再结晶过程的元胞自动机(Cellular automaton ，CA)模型，用来模拟塑性变形过程中动态再结晶过程。再结晶晶粒的生长速度与再结晶驱动力成正比， 再结晶驱动力取决于晶界能和位错密度。在一个增量步内，根据动态再结晶晶粒生长速度确定与其相邻原胞转变 概率。采用本模型对不同温度、应变速度、应变条件下的动态再结晶过程进行了模拟，模拟结果与试验结果以及 经验公式得到的结果相一致。本模型也可用于再结晶动力学、微观组织和织构变化过程的模拟。 关键词：组织模拟 塑性成形 动态再结晶 CA 法 中图分类号：TB115 有联系，无法与有限元等宏观模拟方法相耦合。 * 0 前言 近年来一种新的计算方法元胞自动机(Cellular automaton ，简称CA)法得到了人们的重视，它是基 金属材料的微观组织对材料的力学性能有很 于生物体发育中的细胞自我复制原理，将一个长方 大的影响，控制塑性成形过程微观组织的变化对优 形平面(或立方体)划分成若干个网格，每个网格表 化产品的力学性能具有重要的意义。动态再结晶是 示一个细胞或系统的基元，用某一数值表示其所处 热塑性成形过程中的一种普遍冶金物理现象，通常 的状态(如对再结晶问题表示晶格位向) ，基元的演 认为动态再结晶具备以下性质：①应变大于某一临 化用状态的变动来表示。CA 法具有 5 个特性：① 界值是动态再结晶产生的必要条件。②动态再结晶 基元均匀排列。②基元的状态随时间演变。③基元 晶粒是等轴的，名义尺寸取决于变形状态。③动态 的数值是有限的。④基元状态的演变规则是确定或 再结晶通常在晶界或晶内变形带上形核。以往对动 随机的。⑤每个基元的演化规则是局部的，仅同周 态再结晶过程的研究主要是采用试验的方法，构造 围的基元有关系。该方法抓住了简单性与复杂性这 [1~3] 一对主要矛盾，从而触及并体现了其他有关矛盾， 经验或半经验公式 ，此领域日本学者做了大量的 [4~5] 如局部与整体、宏观与微观、线性与非线性、决定 工作，并已将此方法用于轧制过程的在线控制 ， 但此方法不能揭示热塑性成形过程复杂的微观物理 性与随机性及数学模型与物理本质之间的矛盾。因 现象，如非均质形核、各向异性等。20 世纪 80 年 此它具有利用简单的、局部的规则和离散的方法描 代一些学者开始采用 Monte Carlo 随机统计方法对 述复杂的、全局的、连续系统的能力，且容易与其 [6,7] 他宏观方法( 如有限元法)相耦合。目前，一些学者 动态再结晶过程进行分析 。该方法以自由能最小 原理为基础模拟动态再结晶过程晶粒形核和生长过