任意多介质区域变尺寸.docVIP

任意多介质区域变尺寸 有限元网格自动剖分 吴振君 施斌 祁长青 (南京大学地球环境计算工程研究所 南京 210093) 【摘 要】以Bowyer/Watson算法为基础，结合Lawson算法等对有限元网格进行Delaunay剖分，引入节点间距函数，在区域内通过布置约束点和边来控制特定位置网格的尺寸，并实现网格的尺寸均匀过渡。采用了多种优化技术对生成的网格进行优化。对Bowyer/Watson算法和Delaunay剖分的一些缺陷及问题作了探讨，并提出了相应的处理办法。编写的网格自动剖分软件具有良好的可视化界面，其图形生成、显示及修改功能十分完善，极大地提高了有限元网格剖分的效率和准确性。 【关键词】有限元，Delaunay三角剖分，逐点插入法，多介质 ，网格尺寸，网格质量，可视化 ?1 引言 有限元法(FEM)广泛应用于土木建筑、岩土工程、地铁隧道、水利工程等各个领域中。网格剖分是有限元分析前处理的主要工作。对于大型的计算问题，有限元网格剖分是一项十分繁杂而艰苦的工作，这成为有限元技术应用的瓶颈。近年来，随着电子计算机及计算技术的飞速发展，网格自动剖分技术取得了卓有成效的进展。 由于三角形有描述方便、处理简单等特性 ,适用于对复杂区域简化处理,因此三角形单元在有限元分析中广泛被采用 。三角形网格剖分目前以Delaunay剖分方法和推进波前法(Advancing Front Technique)最为流行[1,2,3]。作者采用Delaunay剖分方法，用VC++编制了可视化的有限元网格自动剖分软件，实现了对多介质、限制数据、变尺寸、多连通域的任意平面区域三角形网格自动剖分。 2 Delaunay网格剖分的基本思想 定义一：给定平面上的点集,令线段集合 ，在上定义图形，其中，是线段集合的最大子集且满足除端点外的任意两线段互不相交，则称为上的三角剖分。 定义二：给定平面上的点集，对每一个，定义一个区域，对任意点，有不等式,成立，则称为点集的氏图，其中表示与之间的欧氏距离。 定义三：给定平面上的点集及其氏图，令线段集合，并且在氏图上的各自区域内具有共同的边界，则为集合的Delaunay三角剖分(其关系如图1示)。??????????? ? ??????? ?研究证明,同氏图互为对偶图的Delaunay三角剖分图具有如下性质。 (1)空外接圆性质。任何一个三角形的外接圆均不包含其它数据点(见图2)； (2)最小内角最大性质。在所有可能形成的三角剖分中, Delaunay三角剖分中三角形的最小内角之和是最大的。 这两个特性保证了Delaunay三角剖分能够尽可能地避免生成小内角的长薄单元，使三角形能够最接近等角或等边,这也是Delaunay三角剖分的算法依据。 3 Delaunay三角剖分的逐点插入法 Delaunay三角剖分的实现算法很多。Lawson提出了用逐点插入法实现Delaunay三角剖分的算法思想。Lee和Schachter, Bowyer, Watson, Sloan, Macedonio和 Pareschi,Floriani和Puppo,Tsai先后进行了发展和完善，本文是基于Bowyer/Watson算法实现网格剖分的，其基本思想是: 在已有网格中，每次插入一个新节点,在网格中找出包含的三角形(图3a),然后从开始递归搜寻外接圆包含点的三角形，并删除这些三角形，形成一个多边形(图3b阴影部分)，连接它的各个顶点和形成新的三角形(图3c)，这可以保证网格为Delaunay剖分。然后进行迭代，直至所有节点处理完。 　 (1)剖分区域处理。搜寻各个内、外边界及区域内的约束点、边(约束点是为控制网格尺寸而在区域内布置的点，两个约束点就构成了约束边)，生成网格剖分区域。定义介质信息及节点间距控制常数。剖分区域通过屏幕上鼠标绘图完成，参数以可视化方式输入和编辑。对于多种介质组成的区域,不同区域的边界线节点分别按逆时针方向编号。对于多连通域问题,外边界节点的编号为逆时针方向,内边界节点的编号为顺时针方向。这样使得计算区域处于边界线的左侧。区域内约束边若闭合，则同外边界一样，按逆时针方向编号。 (2)确定网格的最大包含区域，得到初始剖分三角形，这保证区域内所有的点都落在这个三角形内部。 (3)节点间距函数。 确定剖分区域后，为了控制区域内网格的尺寸及蔬密过渡，文中引入了节点间距函数。区域边界点和内部约束各点的节点间距控制信息，通过线性插值可以得到边界和约束边上各点的间距常数。区域内部加入新点后该点的间距函数S如下。 式中是该点的初始间距控制常数； 分别是三角形顶点的间距函数值；分别是节点和包含的三角形的顶点坐标。 (4)初始剖分。插入边界及区域内约束边上的点，根据节点间距控制信息线性插值确定插入节点的位置。每加入一点，以B