一种基于金相图片的材料微结构几何模型建立方法与流程

本发明涉及集成计算材料工程领域，具体地涉及一种基于金相图片的材料微结构几何模型的程式化建立方法。

背景技术：

集成计算材料工程(integratedcomputationalmaterialsengineering，icme)将复杂的材料冶金过程和材料力学性能联系起来，通过虚拟实验为材料设计和加工工艺制定提供可靠的辅助信息。材料微结构几何模型作为icme仿真的初始输入数据之一，对最终仿真结果的可靠性以及准确性有着至关重要的作用。

材料微观组织几何模型建立方法主要分为两类：一类是基于微观特征统计数据和随机算法建立的伪微观组织几何模型，此类模型并不能完全反映材料的微结构特征，对于不具有明显统计特性的复杂微观组织结构，过度简化的微结构模型的可靠性将大打折扣；另一类是基于真实微观组织图片建立材料微结构几何模型，此类模型能够真实反映材料微观组织形貌特征，在预测结果可靠性上具有明显的优势。

由于腐蚀或拍摄等原因，材料微观组织金相图片中的晶界轮廓可能并不连续。因此，采用常规的图片处理方法无法将组织晶粒和晶界从图像中分离出来。

技术实现要素：

本发明针对现有技术中存在的缺陷，提出一种基于金相图片的材料微结构几何模型建立方法。该方法能够适应材料金相图片中可能存在的不连续晶界，实现双相多晶材料微结构几何模型的程式化建立。

本发明是通过以下技术方案实现的：

一种基于金相图片的材料微结构几何模型建立方法，包括如下步骤：

步骤1，获取材料微观组织金相图片并将其灰度化，使得图片中任意像素点的像素值为0～1之间的小数；

步骤2，以图片像素作为正方形元胞，灰度值作为元胞状态变量，将灰度图片转换为相同尺寸的元胞自动机模型，模型中单个正方形元胞的边长由原始金相图片标尺以及图片中标尺长度方向包含的像素数目计算得出：lc＝ls/ns，其中lc为元胞边长，ls为金相图标标尺，ns为标尺沿长度方向所包含的像素数目；

步骤3，阈值化与噪点过滤处理，具体过程如下：

(1)基于金相图片中不同相以及晶界之间的灰度差异，采用局部自适应阈值法将灰度图像二值化；

(2)设定一个临界噪点尺寸nc1，当二值化图像中某一像素团内像素点的数目小于nc1时，判定该像素团为图像噪点，并将其从图像中剔除；

步骤4，晶界抽取与第二相分离，具体过程如下：

(1)根据晶界在厚度上所占据的最大元胞数确定元胞自动机算法的邻域范围nb；

(2)对于每一个黑色元胞i，计算其邻居内的黑色元胞数目ni，若ni小于预设的临界值nc2时，当前黑色元胞被判定为晶界，否则为第二相；

步骤5，基体晶粒识别，具体过程如下：

(1)在未识别的基体组织元胞中随机选取一个元胞作为新晶粒的生长核心，赋予新的晶粒序号，同时将其邻居内的未识别的基体组织元胞转变为新晶粒前沿元胞；

(2)历遍所有的前沿元胞，判断其是否满足晶粒生长停止规则。若不满足，该前沿元胞将转变为当前晶粒元胞，并将其邻居内未识别的基体组织元胞转变为晶粒生长前沿元胞；否则前沿元胞转变为未识别基体组织元胞。循环该步骤直至所有前沿元胞都转变为当前晶粒元胞；

(3)重复(1)和(2)，直至所有基体组织元胞都转变为相应的晶粒元胞；

步骤6，晶界消除，具体过程如下：

(1)历遍所有晶界元胞，计算其邻居内不同晶粒元胞的数目，并将其转变为具有最大数目的晶粒；

(2)重复(1)直至所有晶界元胞都转变为相应的晶粒元胞。

优选地，步骤3(1)中所述的局部自适应阈值法，其局部自适应阀值由目标像素点局部邻域块像素值的高斯加权和确定，局部邻域块的大小nl由晶界元胞厚度ngb确定，nl＝(2ngb+1)²。

优选地，步骤3(2)中所述的临界噪点尺寸nc1由用户根据过滤后的质量确定，初始默认值为10。

优选地，步骤4(2)中所述的预设临界值nc2的定义如下：

nc2＝nb(1+2nb)-1。

优选地，步骤5(2)所述的晶粒生长停止规则定义如下：

对于当前前沿元胞2层厚度邻居内已识别的当前晶粒元胞数目ni，当ni小于预设的临界值nc3时，当前前沿元胞将不会转变为当前晶粒元胞。临界值nc3根据平直晶界特征确定，nc3＝10。

较之现有技术而言，本发明针对材料真实组织金相图片提出了一种双相多晶材料微结构几何模型程式化建立方法，并且针对组织金相图片中存在的不连续晶界特征提出了一种新的晶粒分离策略。通过本发明方法建立的材料微结构几何模型可以准确有效地反映双相多晶材料微观组织组成相以及晶粒之间的几何拓扑关系。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：

图1为本发明所述方法的流程图。

图2为本发明所述方法的执行过程示例图。

图3是晶界抽取与第二相分离算法示意图。

图4是晶粒识别与晶界消除算法示意图。

具体实施方式

下面结合说明书附图和实施例对本发明内容进行详细说明：

一种基于金相图片的材料微结构几何模型建立方法，实现流程如图1所示，算法执行过程示例如图2所示，初始金相图片输入如图2中(a)所示，最终输出结果如图2中(f)所示，具体包括如下步骤：

步骤1，获取材料微观组织金相图片并将其灰度化，使得图片中任意像素点的像素值为0～1之间的小数，如图2中(b)所示；

步骤2，以图片像素作为正方形元胞，灰度值作为元胞状态变量，将灰度图片转换为相同尺寸的元胞自动机模型，模型中单个正方形元胞的边长由原始金相图片标尺以及图片中标尺长度方向包含的像素数目计算得出：lc＝ls/ns，其中lc为元胞边长，ls为金相图标标尺，ns为标尺沿长度方向所包含的像素数目；

步骤3，阈值化与噪点过滤处理，具体过程如下：

(1)基于金相图片中不同相以及晶界之间的灰度差异，采用局部自适应阈值法将灰度图像二值化，如图2中(c)所示；

(2)设定一个临界噪点尺寸nc1，当二值化图像中某一像素团内像素点的数目小于nc1时，判定该像素团为图像噪点，并将其从图像中剔除，如图2中(d)所示。

步骤4，晶界抽取与第二相分离，算法示意图如图3所示，算法执行结果如图2中(e)所示，具体过程如下：

(1)根据晶界在厚度上所占据的最大元胞数确定元胞自动机算法的邻域范围nb；

(2)对于每一个黑色元胞i，计算其邻居内的黑色元胞数目ni，若ni小于预设的临界值nc2时，当前黑色元胞被判定为晶界，否则为第二相；

步骤5，基体晶粒识别，具体过程如下：

(1)在未识别的基体组织元胞中随机选取一个元胞作为新晶粒的生长核心，赋予新的晶粒序号，同时将其邻居内的未识别的基体组织元胞转变为新晶粒前沿元胞；

(2)历遍所有的前沿元胞，判断其是否满足晶粒生长停止规则。若不满足，该前沿元胞将转变为当前晶粒元胞，并将其邻居内未识别的基体组织元胞转变为晶粒生长前沿元胞；否则前沿元胞转变为未识别基体组织元胞。循环该步骤直至所有前沿元胞都转变为当前晶粒元胞；

(3)重复(1)和(2)，直至所有基体组织元胞都转变为相应的晶粒元胞，如图4中(a)、(b)、(c)和(d)所示；

步骤6，晶界消除，具体过程如下：

(1)历遍所有晶界元胞，计算其邻居内不同晶粒元胞的数目，并将其转变为具有最大数目的晶粒；

(2)重复(1)直至所有晶界元胞都转变为相应的晶粒元胞，如图4中(e)、(f)所示。

优选地，步骤3(1)中所述的局部自适应阈值法，其局部自适应阀值由任意像素点局部邻域块的高斯加权和确定，局部邻域块的大小nl由晶界元胞厚度ngb确定，nl＝(2ngb+1)²。

优选地，步骤3(2)中所述的临界噪点尺寸nc1由用户根据过滤后的质量确定，初始默认值为10。

优选地，步骤4(2)中所述的预设临界值nc2的定义如下：

nc2＝nb(1+2nb)-1。

优选地，步骤5(2)所述的晶粒生长停止规则定义如下：

对于当前前沿元胞2层厚度邻居内已识别的当前晶粒元胞数目ni，当ni小于预设的临界值nc3时，当前前沿元胞将不会转变为当前晶粒元胞。临界值nc3根据平直晶界特征确定，nc3＝10。

以上所述仅为本发明的优选实施例而已，其并非因此限制本发明的保护范围，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，通过常规的替代或者能够实现相同的功能在不脱离本发明的原理和精神的情况下对这些实施例进行变化、修改、替换、整合和参数变更均落入本发明的保护范围内。