用matlab建立晶体模型,一种建立多晶体几何模型的方法与流程

本发明涉及材料结构建模技术领域，具体地说是一种基于voronoi图空间分割原理，利用ug/openapi编程构建多晶体微观结构几何建模的方法。

背景技术：

在对多晶材料进行塑性变形研究时，合适的多晶体几何模型的建立是尤为重要的。早期研究的模型是根据经典塑性理论建立的，经典塑性理论建立在均匀连续介质假设和宏观唯象实验的基础上，对于解决如大变形、织构形成和演化、塑性各向异性、变形局部化、高温和高应变等现代工程实际和高技术领域所面临的问题，有着很大局限性。因此要想解决这些问题就需要从微观层次构建更加符合实际的模型。而且随着计算机技术的迅速发展、各种实验测量技术的发展完善，为多晶材料的研究提供了更加符合实际、更加复杂理论的物质条件。

通常研究多晶变形的有限元模型多以正六边形或以每一个单元代表一个晶粒，这种类似的建模方法无法对细观结构下晶粒内部变形的不均匀性进行研究，而且不能够描述晶界对多晶材料的影响，导致无法研究晶粒演变之间的关系。

由于实际工程材料晶粒大小不一，形状各异，建立与实际组织相近的多晶几何模型是目前有限元分析的关键，其中包括建立符合实际形状的晶粒，以及适合的晶界有限元模型。

现有一些利用扫描电镜等电子设备，通过电子设备扫描样品的表面信息，经过后期处理信息合成相应的晶体模型。这种方法工作量大，最终构建出的模型与真实模型差距比较明显。而且通常由x射线测定的晶体学织构是从统计学的角度出发观察多晶体取向分布状况,因而是从宏观的角度分析问题.宏观织构往往不能确定微观结构和晶粒的单个取向。

目前多晶体建模应用较多的有voronoi图、蒙特卡罗和元胞自动机等，其中蒙特卡罗由于对基础风险因素的一些假设，使得存在模型风险，计算量大，准确性的提高速度慢，如果计算机模拟产生的是伪随机数，那么可能导致错误结果。

再结晶模拟的元胞自动机方法是hesselbarth等提出的，他们建立的模型由于结构简单，易于将复杂的物理现象用简单的转化规则表达出来，但其中也有一些问题没有解决，如再结晶分数-时间关系的定量描述、晶粒尺寸的分布规律及如何显示微观组织演化过程等。

voronoi图的方法可以与有限元方法相耦