材料微观原子排列的三维显示算法研究.docVIP

材料微观原子排列的三维显示算法研究 　　摘要：材料模拟系统常用到微观结构的三维模拟显示技术，以反映材料内部微观结构，并可提高计算材料研究结果的合理性。但在应用系统开发中，对于微观结构三维模拟显示较为困难，利用各种3D工具软件进行显示一般难于实现交互，细微调整及自编程序3D模型的建立也非常困难。因此，该文提出了球冠逐步覆盖算法，建立材料微观原子排列的3D描述模型，以此找到一种材料微观结构3D描述的方法，并以DELPHI7.0的OpenGL完成模拟图像测试输出，效果较好。该方法适合在实际应用项目中推广使用。 　　关键词：材料计算；微观模拟；3D显示；算法 　　中图分类号：TP274文献标识码：A文章编号：1009-3044(2011)12-2932-04 　　材料制造业是国民经济的物质基础、国家安全的重要保障和国家竞争力的主要体现。近年来，随着材料学的深入研究开展，新材料微观结构三维显示研究逐渐成为发展新材料需要解决的问题。 　　传统新材料的开发采用的是代表性体积元法，通过假设出的一定体积元，来代表具有规则的、周期性分布的微观结构[1]。这种方法一般用于非均质材料，对于多晶体材料则不能表示具有不同尺寸、形状等参数的晶粒组织。以上方法还要配合复杂的材料计算、反复的试验测试，这无疑将提高创造新材料的成本。随着计算机技术在材料学中的应用日益广泛，研究新材料已不再盲目地依靠实验室中的重复试验，而是运用科学的算法在日益强大的计算机中模拟材料的微观结构仿真模型、材料的抗拉伸及抗腐蚀等各种性能。这种通过电脑对材料进行显微分析的方法，使研究人员能够更好、更直观、快速的理解材料的深层原理，从而有效缩短研究周期，降低开发成本。 　　材料科学主要研究材料的晶粒组织、亚组织等微观的结构和分布,然后才能对结构和分子的分布进行分析和对比，进而开展材料科学方面的其他研究。建立材料的显微模型，能够有效地反映材料微观组织的晶体分布、结晶数量、聚合程度等，这些数据与材料各外在属性是密切相关的。所以，材料学中结构的可视化显示在材料学中占有很重要的作用。而3D可视化显示在材料计算中广泛需要，是三维空间图像表现即原子排列问题是其中一个主要问题[2]。本文即主要探讨球冠逐步覆盖算法及利用OPEN GL实现的方法[3-9]，研究球体表面密布算法问题的求解，实现在一个球体表面上布满小球且小球均匀分布不能叠加地进行显示，以及材料模拟原子排列三维微观组织结构的显示，拓展材料学研究的领域。 　　1 三维微观组织结构显示的覆盖算法分析 　　通常人们可以考虑利用空间几何思想求解原子间排列的密布问题，假设一个球体集，半径为R，在其表面密布着一层小球，小球半径为r（小球不得重叠且可以相切）。因此，问题就在于，要找出若干个半径为r的球体，能够覆盖半径为R的球体的表面。起初我们是将大球体按南北极经纬线的方法划分，并在其极点及各经纬线交点处放置小球（大球半径直接关系到小球密布的数量），但实际上是不可行的。原因在于处在极点附近的小球无论半径设定多小都会出现重叠问题，此法不能使用。 　　球内接正多边形的方法，试想球体内接一个正多边形，在正多边形与球面相交的每个交点上放置一个小球，方法可行。但带来的问题是交点灵活，小球的半径与正多边形的边长密切相关，而且还要考虑正几面体可以直观的显现小球密布大球最复杂的是内接多面体边数变小球半径随着变。该方案也不适合。 　　最后，我们收集各种资料整理思路决定采用宏观看问题的方法，既不锁定球的极点也不注重球面交点，从类似足球的表面图形着手，即球的内接正12面体，以每个面的中心为小球球心所在的轴线，逐步思考推算坐标，来解决球体的密布问题。 　　因为球内接正12面多边形酷似菠萝，故称这种算法传统被称为菠萝形算法。 　　菠萝型算法：在密布问题中考虑大球的下列内接多面体。它有两个面是正n边形A1A2...An和B1B2...Bn，其外接圆的半径都是R，并且这两个面彼此互相平行，另外的2n个面均是全等的五边形，外接圆的半径也是R（图1）。当n=5时，这个多面体就是球的内接正12面体。因为这个多面体形如菠萝，故称这种算法为菠萝形算法，这种算法共需2n+2个小球。 　　采用球面极坐标，并设R=1，设O1O2是垂直于面A1，A2，...，An的直径，它与该面的交点O3, 即为该正n边形的中心。取经过O1O2与五边形A1A2C2D1C1的顶点D1的半平面为半平面θ=0，而O4是这个五边形的外接圆圆心。那么,五边形A1A2C2D1C1关于平面θ=0是对称的，且A1、C1在同一条径线上，因此可设这些点的坐标为（设R=sinφ0）：O3（cosφ0,0,0）,A1(1, φ0, π/n)，O4（cosφ0, φ1,0）,C1(1, φ2, π/n ),D1(1, π－φ