STL文件上有限元网格的生成.doc

STL文件上有限元网格的生成 摘要： 这些被提议的方法的是为了展示一个适合有限元方法的且关于物体近似边界表示的直接由CAD软件生成的曲面网格产生的可能性。首先，我们将描述由一个简单的物体表面三角剖分组成的边界表示法。接着我们将展示如何获得一个相容的形状适应网格。形状适应是考虑到几何逼近和由一个误差估计量的各向同性的尺寸映射做出的。网格可以被用来用于有限元计算（通过壳体元素），或者可以用来作为一个启动体积网格算法（Delaunay 或advancing front）曲面的网格。这个用来产生网格的原则是基于与求精算法相关的Delaunay方法和光滑化原则。最后，我们将展示不用用于克服常规的基于一个几何特定表示法的网格化软件的限制的几何模型的参数表示法。 关键字：网格产生;STL文件格式Bisection算法;Delaunay 三角剖分。 CAD接口 1.1基于CAD的数据集 市面上的大多数CAD软件可以产生STL文件，且这些大体上用作实现样机研究和绘制图形的目的。这些文件表现了立体的三角剖分的边界。STL三角剖分的生产算法都是高效的，而且如果可以接受大规模数据集的话平面可以被精确地逼近。但是，这不是一个真正几何模型是因为STL文件格式仅由一个扩展的三角平面列组成。这些小平面除了一般立体外部定向的座标外，还由三个三角顶点座标组成。这种三角剖分法是建立用来最小化一个关于立体(图1)真实边界的几何逼近准则的。（见图1） 1.2 STL三角剖分的特性 STL三角剖分不能直接用在有限元方法(FEM)中，主要是因为它要求在计算域的几何描述上的特殊性。在FEM中，几何和函数支持由元（三角形或其他）提供，且它们必须有一特殊形状，也就是就计算误差估算量的固有大小和固有品质因子应该越低越好。在这项工作中，我们利用如下三角形的品质因子，在此di,i=0…2代表三角形每边的长度。我们假定: 品质因子就是: 这个品质因子在0（对于所有退化三角形而言）和1（对等边三角形而言）之间。 当然，这些要求是和应用相关的。比如，在流体力学中很普遍地应用在边界层的各向异性元和冲激波当中。很显然一个由最小化几何准则得到的网格不能适合FEM的要求，因为它依赖于表面的的曲率和拓扑，且生成的三角形在某方向上会大大伸长（见图2－4）。但是，得到的网格一般都是一致的。 1.3 STL网格中的几何恢复 STL文件内容包含每个三角形顶角座标和相应法线。为了达到曲面的完全网格化，我们需要获得关于被网格化曲面的拓扑和曲率的有关数据。三角形间的拓扑和连通度可以通过避免在STL文件间的顶角冗余得到。这可以通过利用一个利用字典法排序和储存了的顶角的二叉树做出的。在这个过程结束时，我们得到了一物体表面的网格及所有类型的可被重新网格化所需的连通度，特别是用来作邻接查找的连通度（见小节三）。 由于许多工业制品有棱与边，我们需要决定它们定位的位置。这可以通过一个边测定程序考查在三角剖分中的实体间的连通度得到。在保持它们原貌以用于重新网格化的过程中模型边的查找显得十分重要。也就是说它们不能被Delaunauy重新网格化法切割。这种查找是基于邻接三角形的角的计算。两三角形的公共部分被认为是一个边，如果在此两个三角形的法向量间的角比一个特定的值（一般为20。）大的话。另外还应该注意的是我们可以处理诸如在一个集合（比如为了应用边界条件）中的两个不同部分间的接触部分的面积。这些部分将被认为是边，且将被保留在最后的网格中（图5）。 接着曲率通过寻找其它与之连接的三角剖分的每个元素来求值且检查它们的法向量（图6）。 接下来的程序是用来：对于每个三角形Ti,我们在它的三个顶点（Vj(j=0…2)进行迭代。在每个顶点，有两个边Sjk(k=0…1),与之相关的三角与Ti，Tjk邻接。如果Sjk是个陡沿，我们定义一个变量为0;不然就是1。曲率通过下式计算出来： j是 在Sj0和Sj1间的角，’ j是法向量Tj0和Tj1间的角。 在这之后，曲率在每个项点Vi外推(Cvi),应用连接Vi的n个三角形的最大曲率： 我们应该注意到在这里考虑的曲率是关于曲面的曲率张量范数（也就是说我们不用在意曲率产生的效应方向），要注意到得到的曲率不是十分准确的。这是因为几何曲面表示（平面集）是离散。 实际上，这些曲率不仅驱动一任选网格的细化。在较少弯曲平面上的结果是粗糙的，在其它地方可得到改善。图六描绘了利用上面提到的公式得到的一个穿过每个三角形的线性插值的结果。 在这方面，我们需要的用于网格生成的程序的信息是关于目标网格的各向同性尺寸映射的信息。既然这不是当前研究工作的目的，我们假定或者常尺寸映射是由用户提供的，或者一个尺寸映射是基于一个已被提供的误差估计量（从先前的FEM算得）。我们也考虑当密度映射是基于一个几何准则（控制在STL网格和三角