[发明专利]基于三维随机重构的无机非金属材料性能分析方法

权利要求书

1.一种基于三维随机重构的无机非金属材料性能分析方法，其特征在于：该方法是在无机非金属材料分析常用手段即电子显微镜的基础上，通过连续多次不同深度的扫描，获得所有截面的二维图像；再将二维图像处理成晶体成分特征矩阵，多个二维晶体成分特征矩阵进行重构，生成三维晶体特征矩阵；同时提取晶体成分特征，形成特征库，为进一步的三维随机重构做准备；基于特征库，对三维晶体特征矩阵随机重构，形成新的矩阵，并转化为有限元软件的输入文件，最后进行计算分析，得到无机非金属材料的性能参数；

主要包括如下步骤：

1)二维图像获取：通过电子显微镜获得无机非金属材料的连续不同深度的二维图像，每一层的深度相同，采样区域保持固定不变，图像尺寸大小为X×Y像素，共采集n次；

2)二维图像参数化：对获取的二维图像进行处理，将基材的颜色去除，同时不同成分的晶体设定不同的颜色，通过MATLAB读取图像文件，获取RGB值，建立晶体成分与RGB值的映射关系，通过RGB值生成二维矩阵，大小为X×Y；对晶体成分赋值；基材的晶体赋值0，晶体成分1赋值1，晶体成分2赋值2，依次类推，晶体成分i赋值i；

3)三维重构：将每一个X×Y大小的二维晶体成分矩阵转化为新三维矩阵中的一层，形成X×Y×n的三维矩阵；此矩阵中的每一个数值i代表一个无机非金属材料单元，其大小为单个像素高度×单个像素宽度×取样深度；

4)三维结构样本随机重构：通过记录三维晶体成分矩阵中各个数值的特征参数，所述特征参数包括分布位置、大小、形态和体积占比，形成晶体特征库；以晶体特征库为基础，采用真实随机重构方法进行三维结构参数矩阵重构；

所述真实随机重构方法即以晶体特征库为基础，模拟晶体的大小、形状、分布状态的特征，建立符合真实构造的新三维晶体；

5)有限元计算分析：通过有限元分析软件，建立长度为X像素，宽度为Y像素，厚度为(n-1)×取样深度的长方体，按照三维结构参数矩阵的元素个数划分网格，得到X×Y×n个网格单元；通过修改对应网格单元的材料属性，实现无机非金属材料的微观建模，再进行计算分析，从而得到该无机非金属材料的性能参数。

2.根据权利要求1所述的基于三维随机重构的无机非金属材料性能分析方法，其特征在于：

具体步骤如下：

所述步骤1)中，选取需要分析的无机非金属材料样本，获得该无机非金属材料的成分特性，所述成分特性包括密度、弹性模量和主要成分；确定扫描区域，确定单次扫描递进深度，描完上一层后，将试样上一层去除再进行下一次扫描，以获得连续的无机非金属材料电子显微镜图像，共扫描n层；

所述步骤2)中，将二维图像参数化；保持扫描区域固定不变，将二维图像分割成X像素×Y像素的区域，通过MATLAB进行RGB读取，得到像素矩阵，也就是R、G、B信息分别保存在像素矩阵中；合成二维晶体成分矩阵；处理RGB矩阵，分别提取3层像素矩阵中位于(x，y)位置的数值：A1＝(x,y,1),A2＝(x,y,2),A3＝(x,y,3),生产数组Sx,y,1＝(A1,A2,A3)，通过模糊识别处理对应RGB-晶体成分映射库，将Sx,y,1转化为自然数1；无机非金属材料的基材映射数值设置为0，假设有i种晶体成分，就有i个数组Sx,y,i，Sx,y,i转化为自然数i；再将自然数i按其对应的x，y位置一一对应，生成一个新二维矩阵B1其元素为0～i的自然数，此矩阵表示该层电子显微镜图像的晶体成分；

所述步骤3)中，通过多个二维晶体成分矩阵，合成一个三维晶体成分矩阵C，其大小为1000×1000×100；提取晶体特征，包括晶体的占比、体积大小、分布形态；提取矩阵中相同元素数值i的占比：即统计数值i在整个X×Y×n矩阵中出现的次数，计算其在矩阵中的元素占比；具体意义为：晶体成分i在整个扫描无机非金属材料试样中的体积占比；

提取矩阵中相同元素数值i聚集体元素个数和、单个聚集体在元素i中的占比；单个聚集体之和：即扫描整个X×Y×n矩阵，当出现i元素时开始计数，假设第一个i元素的位置为(x_j,y_k,n_l)，扫描其相邻位置(x_j±1,y_k±1,n_l±1)的元素值，如出现数值i则扫描边界位置进行扩充，直至所有边界之外都没有元素i，E_ci为单个元素i聚集体的元素个数和；其物理意义为晶体成分i在无机非金属材料中形成的单个聚集体的大小；

提取单个元素数值i聚集体的体积占比：即单个聚集体E_ci在所有元素i中的占比；物理意义为，晶体成分i在无机非金属材料中形成的单个聚集体在晶体成分i中的占比，也即是分布情况；

提取单个元素数值i聚集体的形态特征值：最大元素距离L_ci，在最大元素距离方向的法向上提取聚集体的高度H_ci和宽度W_ci；物理意义为，单个元素数值i聚集体的长宽高；

所述步骤4)中，三维晶体成分矩阵随机重构，基于晶体成分特征库，使用真实随机重构方法，对三维晶体成分矩阵C进行重构，得到新的三维晶体成分矩阵C_f；物理意义为按照晶体成分特征库的规则，对无机非金属材料进行成分随机重构，从而得到f个新材料；

所述步骤5)中，生成有限元计算文件；以ABAQUS有限元软件为例，首先建立一个500μm×500μm×(5×99)μm的立方体，网格划分为1000×1000×99个，其中长度和宽度方向上均为1000个网格，厚度方向上为99个，设置边界条件，保存为inp格式文件；这样，ABAQUS有限元软件中建立的三维实体和网格单元与我们实际中的无机非金属材料是完全对应的，通过修改inp文件，将三维晶体成分矩阵C_f中不同元素的材料参数写入inp文件中相应单元中，就完成了有限元计算文件的生成；材料性能参数输出，通过ABAQUS有限元软件的并行处理功能，将f个计算文件提交处理，最终获得f个计算结果，即f个性能参数。