基于图像的自适应有限元网格划分方法技术

一种基于图像的自适应有限元网格划分方法，属于材料结构有限元分析的网格划分方法。方法为：一、利用图像分析法区分图像中不同的相或组分；二、利用计算机语言程序，将所有像素的位置信息读入计算机内存中；三、根据空间结构的精细程度，对结构粗大，成分单一的区域降低图像分辨率；四、将像素直接构建成有限元单元，不同分辨率造成不同的网格密度；五、以映射网格过渡不同密度的网格。本发明专利技术利用图像直接构建节点和单元，跳过实体建模，因而能够节省系统资源。并且能够在保证计算精度得前提下优化网格结构，显著降低节点数量，节省计算时间和内存消耗。全过程自动完成，对于节点数量巨大的模型非常实用。

【技术实现步骤摘要】

本专利技术涉及一种材料结构有限元分析的网格划分方法，特别是一种基于图像的自适应有限元网格划分方法。
技术介绍
材料的性能对于其使用效能有决定性的作用，因此测试材料的性能是在使用它前必做的一项功课。然而对于结构材料设计来说，常常要频繁的改变材料的结构以期能得到想要的性能，这使得需要大量的测试工作。通过材料结构预测材料的性能，能够极大的缩短材料设计的周期，因此大量应用在复合材料设计等领域。传统的预测方法一般都使用两种方法：1)基于理想假设的理论计算和2)基于实验的经验公式，或者将两者结合起来。对于第一种方法，实际情况往往比理想状态复杂的多，因此理论计算的结果和实际状态可能会有较大偏差。而对于第二种方法，由于实验条件的不同，由经验公式计算的值也可能会与真实结果有较大的差别。这使得无论是以上哪种方法都很难适用于普遍场合。近些年随着计算机技术的发展，人们可以利用计算机程序对于简单的材料结构进行虚拟构建，然后进行有限元等数值建模，再设置合理的边界条件和初始条件进行计算，能够对材料的性能进行分析，这种方法即为计算机辅助设计CAD(ComputerAidedDesign)。由于该方法是利用材料的真实结构进行建模的，因此计算结果有较好的准确率。目前，这种方法广泛应用于复合材料结构设计领域。然而，对于一些没有周期结构的、非均质的复合材料，例如热喷涂WC-Co涂层，由于其结构太复杂，很难进行构建，因此该方法也难有用武之地。数字摄影方法是将连续的图像信息存储为离散的像素，堆叠起来逼近原图像。这和有限元计算中将连续的空间离散成单元的思路很相近，因此便诞生了基于图像的数值建模方法。基于图像的数值建模是指利用能够代表材料结构的二维或三维数码照片，建立数值模型的一种方法。该方法直接以材料的图像作为其结构，常以像素为单元，构建材料的数值计算模型，因此能够适用于各种复杂的材料结构。在近十年来，有学者利用该方法对矿物材料、建筑材料和涂层材料进行建模，成功预测了这些材料的力学、热学等性能。当材料的结构精细且复杂无规律时，需要有较高图像的分辨率来显示材料的结构。然而当计算区域大小一定时，分辨率越高，像素数就越多，以像素为单元的模型就越大，计算耗时和对系统资源的消耗就会指数上升。特别是对三维图像建模时，由于节点和单元体数量是同边长二维模型的1.5次方倍左右，使得分析时间和所需计算机资源极度增加，这可能会失去数值模拟简便快捷的优势。因此，在精细复杂的组织区域保持高的单元密度，在简单粗大的组织区域降低单元密度，建立对于结构精细程度的自适应网格，能够在保证计算精度的同时，提高计算效率。目前的有限元网格划分方法，都需要先建立实体模型，因而并不适用于图像有限元建模方法。因此本专利技术提出了一种基于图像像素的自适应网格划分方法。
技术实现思路
为了在兼顾计算精度和提高计算效率，保持数值模拟在材料性能预测周期上的优势，本专利技术提供一种根据材料结构，自动控制网格的粗细程度，基于图像像素的自适应网格划分方法，该方法对于节点量大的三维问题尤其适用。本专利技术解决其技术问题所采用的技术方案是：该划分方法：首先以适宜的分辨率摄取材料的数值图像，在组织结构精细的区域保持高密度单元，在组织结构粗大的区域降低节点和单元的密度，通过合理的过渡单元连接上述两种尺度的单元。具体包括以下步骤：(a)获取能代表材料结构的二维或三维图像；(b)以图像分析法区分图像中不同的相或组分；(c)将所有像素的位置信息读入计算机内存，标注出他们属于哪个相或组分；(d)降低成分均一、结构简单的区域的图像分辨率；(e)除了不同分辨率之间的过渡区，把其余区域每个像素直接建成一个单元；(f)建立不同分辨率之间过渡区域的单元；(g)构建单元模型。所述的步骤(d)中，对应实体结构中成分均一、结构简单的区域，只要保证像素对应的材料单一性，则不断循环地减半该区域的图像分辨率。所述的步骤(e)和(f)中，直接根据像素的材料归类和位置信息构建有限元单元，节点的位置对应像素的顶点。所述的步骤(f)中，不同分辨率的区域之间，二维网格采用映射型的多边形过渡连接，三维网格采用映射型的多面体进行过渡；根据材料各部分区域不同的结构精细程度，选择不同节点网格密度，越精细的区域节点密度越高，反之节点密度越低。有益效果，由于采用了上述方案，本专利技术为一种利用材料图像的有限元建模方法，其作用为改善了以图像像素直接构建成有限元单元的图像有限元建模方法，建模思想是基于可调节的图像分辨率原理，对于成分单一或组织粗大的区域，降低该区域的图像分辨率，也就能够降低该区域的网格密度。采用多级图像分辨率，即可有效降低有限元模型的节点数量，节省了系统资源和计算时间。优点是：基于图像像素直接构建有限元单元，避免构建实体模型，降低系统资源消耗，提高运算速度。应用于组织结构较为复杂的材料时，相比以每个像素为单元的建模方式，能够明显减少节点数量，降低系统资源消耗，在保证计算精度的前提下提高运算速度。网格划分采用计算机程序自动完成，操作简单方便。附图说明图1是本专利技术基于图像的自适应有限元网格划分方法的流程示意图。图2是本专利技术冷喷涂SiC增强的镁合金复合材料的二维断面图像(1024×768像素)。图3是图2经过图像处理后的照片，其中黑色是孔隙，白色为SiC，其余为镁合金。图4是以图2为建模对象的二维自适应有限元网格图。图5是本专利技术等离子喷涂氧化钇稳定氧化锆涂层的三维重建图像(80×80×50像素)。图6是以图4为建模对象的三维自适应有限元网格图。具体实施方式该划分方法：首先以适宜的分辨率摄取材料的数值图像，在组织结构精细的区域保持高密度单元，在组织结构粗大的区域降低节点和单元的密度，通过合理的过渡单元连接上述两种尺度的单元。具体操作步骤如下：第一步，获取能代表材料结构的二维或三维图像；第二步，以图像分析法区分图像中不同的相或组分；第三步，利用计算机语言程序，将所有像素的位置信息读入计算机内存中，标注出他们属于哪个相或组分；第四步，降低成分均一结构简单的区域的图像分辨率；第五步，除了不同分辨率之间的过渡区，把其余区域每个像素直接建成一个单元；第六步，建立不同分辨率之间过渡区域的单元；第七步，构建有限元模型。所述第一步，摄取图像要选择合适的分辨率和面域或体域大小；图像的分辨率选择要以能清楚展示材料的结构为准，不宜过高或过低；摄取图像的面域或体域大小要能代表材料的真实结构，对于周期性排本文档来自技高网...

【技术保护点】
一种基于图像的有限元自适应网格划分方法，其特征在于：该划分方法：首先以适宜的分辨率摄取材料的数值图像，在组织结构精细的区域保持高密度单元，在组织结构粗大的区域降低节点和单元的密度，通过合理的过渡单元连接上述两种尺度的单元；具体步骤如下：（a）获取能代表材料结构的二维或三维图像；（b）以图像分析法区分图像中不同的相或组分；（c）将所有像素的位置信息读入计算机内存，标注出他们属于哪个相或组分；（d）降低成分均一、结构简单的区域的图像分辨率；（e）除了不同分辨率之间的过渡区，把其余区域每个像素直接建成一个单元；（f）建立不同分辨率之间过渡区域的单元；（g）构建单元模型。

【技术特征摘要】
1.一种基于图像的有限元自适应网格划分方法，其特征在于：该划分方法：首先以适宜的分辨率摄取材料的数值图像，在组织结构精细的区域保持高密度单元，在组织结构粗大的区域降低节点和单元的密度，通过合理的过渡单元连接上述两种尺度的单元；具体步骤如下：
（a）获取能代表材料结构的二维或三维图像；
（b）以图像分析法区分图像中不同的相或组分；
（c）将所有像素的位置信息读入计算机内存，标注出他们属于哪个相或组分；
（d）降低成分均一、结构简单的区域的图像分辨率；
（e）除了不同分辨率之间的过渡区，把其余区域每个像素直接建成一个单元；
（f）建立不同分辨率之间过渡区域的单元；
（g）构建单元模型。
2.根据权利要求1所述的基于图像的有限元自适应网格划分方法，...

【专利技术属性】
技术研发人员：乔江浩，谭娜，刘洪涛，张德坤，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：江苏;32