由区域生长算法实现四边形网格划分

合集下载

区域算法的分割方法

区域算法的分割方法
区域算法的分割方法是一种将图像分成若干个区域的图像分割技术。

这些区域通常具有相似的特征，例如颜色、纹理、亮度等，并且在图像内部具有连续性。

常见的基于区域的分割方法包括以下几种：
1. 区域生长法（Region Growing）：从某个起始点开始，不断地向周围生长，直到达到一定条件为止，例如颜色、纹理或灰度值等。

这种方法需要手动选择起始点，并且计算复杂度较高。

2. 分裂与合并法（Split and Merge）：将整张图像看作一个区域，将其分割成若干个子区域，然后对相邻的区域进行比较，并将相似的区域合并成更大的区域。

这种方法可以自动选择区域，并且可以产生较平滑的分割结果。

3. 基于能量函数的方法：通过定义一个能量函数来分割图像。

能量函数是一种用于衡量分割质量的函数，通常包括区域的特征和区域之间的相似度。

通过最小化能量函数来得到最佳的分割结果。

4. 基于图论的方法：将图像表示为一个图，其中每个像素表示一个节点，节点之间的边表示像素之间的相似度。

然后使用图分割算法来将图分割成若干个子图，每个子图对应一个区域。

这种方法可以自动选择区域，并且通常具有较高的分割质量。

这些基于区域的分割方法可以根据具体的应用场景和需求进行选择和应用。

复杂区域四边形网格生成的一种改进方法

年来，非结构化网格的生成方法取得了快速发展。
行波法（ｄｎｎｒｔｎｕ也称作推进ａｉｆｎｔｈｉｅ（ｖｃｇｅｑ）Ａｏｃ波前最初是由Ｌｌ法，ｏ］ｌ提出，并用于平面区域三角形网动生后来，ｅｉ［将这种格自成。Ｐｒｒ１ａｅ等２方法进
（找Ｎ足件满条礁Ｐ集Ｎ使４出中）的合２，
Ａ
得Ｐ＜，＜，Ｄ与生成波上除ＡＡ２Ｐｎ且ＰｎＬＢ２ＬＢ外的任一线段都不相交，本文取ｎ１；２．＝５（）若Ｎ不是空集，则选取其中相对于Ａ５２Ｂ的张角ＬＰＡＢ最大的点Ｐ作为新生成点，构造单元ＡＢＡＰ；若丛是空集，ＡＣ中不包含点集凡但ＡＢ
一步完善。Ｐｒｒ等发展的技术中，在ｅｉａｅ区域中的内部节点与单元是在剖分过程中同步生成的，并且通过引入背景网格来控制内部节点的生成，这一改进便于将网格与离散误差联系起来，有利于寻求最优
网随后，ｅｉ及Ｌｈｅ３格。Ｐｒｒｏｎｒ同时成功地把这种ａｅ１】
技术推广到任意形状三维区域的四面体网格生成。与其它各种网格生成方法相比，行波法对复杂几何形状与边界的二维区域（或三维区域）中的三角形（或
２三角形网格生成方法
一般来说，对一个给定的区域进行四边形网格划分可以分两步：首先，在整个区域内生成三角形网格；然后，在三角形网格的基础上划分四边形网格。利用波前法生成三角形网格的技术己经相对成
熟，这里只简要介绍一下其生成步骤（见图１）０
、、
‘．‘ ．沪．口．口产口
３区域描述．１通常直接对复杂区域生成网格比较困难，较为方便的作法是根据材料介质分区、网格密度分布、边界条件等进行区域分解，把复杂区域分解成若干个较为简单的子域。每个子域分别生成网最后，格，将各个子域的网格连接在一起，就构成了区域的网

基于带约束前沿推进的四边形网格生成方法

基于带约束前沿推进的四边形网格生成方法汪攀;张见明;韩磊;鞠传明;池宝涛【摘要】In view of the limitation of the current various Q-Morph algorithms that can generate isolated triangles in the process of generating quadrilateral mesh,a new advancing front method with constraint was proposed.The method is a kind of Q-Morph algorithm based on advancing front.First,the two side edges and one top edge are generated according to the current front.Then the triangles in the quadrangle are eliminated from the triangular domain.Finally a quadrilateral mesh is generated.Moreover,some constraining conditions are imposed on the advancing procedure to avoid generating isolated triangles,which can improve algorithm efficiency and meshquality.Numerical experiments are presented to demonstrate that the new method could generate high-quality quadrilateral meshes for complex surface automatically.%针对当前各种Q-Morph算法在生成四边形网格的过程中容易产生残余三角形这一缺陷,提出了一种带约束的前沿推进算法.该方法是一种基于前沿推进思想的Q-Morph算法,由当前前沿生成两条侧边和一条顶边,并删除其内部的三角形,从而将三角形网格合并生成一个四边形网格,并在前沿推进的过程中加入了约束条件,从而可以有效地避免残余三角形的产生,提高了算法的效率和最终生成网格的质量.数值实验表明,该算法能够全自动生成质量较好的四边形网格.【期刊名称】《湖南大学学报（自然科学版）》【年(卷),期】2017(044)008【总页数】6页(P29-34)【关键词】网格生成;四边形网格;Q-Morph;前沿推进【作者】汪攀;张见明;韩磊;鞠传明;池宝涛【作者单位】湖南大学汽车车身先进设计制造国家重点实验室,湖南长沙 410082;湖南大学汽车车身先进设计制造国家重点实验室,湖南长沙 410082;湖南大学汽车车身先进设计制造国家重点实验室,湖南长沙 410082;湖南大学汽车车身先进设计制造国家重点实验室,湖南长沙 410082;湖南大学汽车车身先进设计制造国家重点实验室,湖南长沙 410082【正文语种】中文【中图分类】TP391在求解各类工程和科学计算问题上，有限元法或边界元法[1-2]已成为最有效的数值分析方法之一，网格划分是数值分析前处理的重要部分，已成为当前有限元或边界研究的重要领域.四边形网格既可以直接用于壳分析，也可用于六面体网格的输入，而且四边形网格在计算精度、计算效率等方面均优于三角形网格，因此研究四边形网格生成算法很有必要.四边形网格生成算法分为间接法和直接法.其中直接法包括映射法、子域映射法[3-4]和铺砖法[5-6]等，这些算法往往比较复杂或自适应性较差.为了减少四边形网格生成的难度，通常采用间接法[7-8].Lee等提出的Q-Morph算法[9]是典型的间接法，间接法即先生成三角形网格，再将三角形网格转化为四边形网格，常见的转化技巧有2种：分解法，将1个三角形分解为3个四边形；合并法，将两个相邻的三角形合并生成1个四边形[10].合并法在合并操作时会考虑生成四边形网格的质量，因此其生成的初始网格要优于分解算法.合并算法最大的难题在于最后会留下很多残余三角形，阻碍了全四边形网格的生成，此时常见的解决办法是用分解算法，将残余三角形分解为3个四边形，再通过拉普拉斯优化进行局部光顺处理，但是处理过多的残余三角形，不仅耗时较多，算法效率较低，也会影响最终生成网格的质量.为了减少残余三角形的数量，Lo[11]尝试设置一些探索性规则来确定三角形网格的合并顺序；胡向红等[10]在三角形合并的过程中加入了约束条件，从而避免了残余三角形的产生，但是该算法生成的网格质量较差.基于前沿推进的思想，Owen等提出了改进的Q-Morph算法[12-13]，该算法改进Lee等提出算法中存在较多不规则点的缺陷，目前这一算法得到广泛的使用，并已集成到大型的商用有限元软件中.该方法可以生成质量较好的四边形网格，但是在合并的过程中也容易产生残余三角形.针对上述算法的不足，本文提出了一种基于带约束前沿推进的四边形网格生成方法，该方法以Owen提出的Q-Morph算法为基础，以保证在前沿推进的过程中，生成的四边形质量较好，并且在合并操作的过程中加入了约束条件，以保证前沿推进的过程中不会产生残余三角形，该算法能够自动高效地生成质量较好的四边形网格. 本文提出的算法是一种间接的四边形网格生成方法.先利用推动波前法生成三角形网格，然后基于Q-Morph算法以及前沿推进的思想，将三角形的外边界作为初始前沿，从外向里推进，将三角形网格合并生成四边形网格，其主要步骤如下，流程图如图1所示.1)用推动波前法生成初始三角形网格；2)取三角形网格的外边界，将其加入前沿链表，作为初始前沿；3)从链表中选取一条前沿作为当前前沿，将其作为当前四边形的底边；4)在底边的两个端点处生成两条侧边，侧边可以是已经存在的边，也可以是通过对角交换获取的边；5)通过交换对角线操作获取顶边；6)将底边、侧边、顶边生成一个四边形，并删除里面包含的三角形；7)通过局部优化改善四边形网格的质量和周围三角形网格的质量；8)更新前沿.基于前沿推进的Q-Morph算法主要包括：生成侧边、生成顶边、删除四边形内的三角形和局部光顺等几个主要的步骤(如图2所示).2.1 侧边的生成下面介绍以AB为前沿边，生成侧边AD和BC，步骤如下：1)确定当前前沿和相邻前沿的角平分线，如图3(a)所示的虚线；2)计算与前沿结点相连的边与角平分的夹角，如图3(a)所示的α1，α2，并取其中的最小值α=min(α1,α2)；3)根据α确定侧边的生成方式，如果α≤ε(本文中ε取π/6)则直接以α对应的边作为侧边，如图3(b)所示的BC.否则，采取交换对角线的方式获取侧边，如图3(d)所示的AD.从侧边的生成过程可以看出，该方式生成的网格质量相比传统的直接删除一条公共边生成的四边形网格质量要好.当α较小时，直接取该边做侧边，当α较大时，则必然存在分居角平分线两侧的边，此时通过交换对角线得到的新边，则必然与角平分线的夹角较小，即适合做侧边，因此，此方法总能找到一条合适的侧边.2.2 顶边的生成在以AB为前沿边，生成侧边AD和BC以后，连接CD，CD即为顶边，但是很可能CD并不是三角形中已经存在的边，此时需要通过交换对角线，恢复边CD.对于任意的三角网，理论可以证明，可以通过有限次对角线交换恢复给定的任意一条边，具体步骤如下：1)连接CD，得到所有与CD相交的三角形.如图4所示的T0,T1,T2,T3；2)由C点开始，依次交换T0,T1,T2,T3的对角线，交换对角线时，如果相邻两个三角形构成的四边形是凹四边形，则将对角线交换的过程延后至下一次循环中进行；3)当出现CD边时，循环终止，即得到了顶边CD，否则进入步骤2).直至恢复了顶边CD.2.3 删除内部三角形生成四边形ABCD以后，需删除四边形内部的三角形，本文删除内部三角形的方法也是基于前沿推进的思想提出来的，具体步骤如下(如图5所示)：1)将四边形内部的三角形加入前沿链表中；2)以当前前沿为基础，搜索与其相邻的三角形，删除该三角形，并判断该三角形另外两条边，若该边为前沿边，则将其从前沿链表中删除，否则将其加入前沿链表；3)当前沿链表为空时，该过程结束，否则进入步骤2)，直至前沿链表为空.2.4 约束条件Owen提出的Q-Morph算法[13]是以前沿推进的方式生成四边形，但是在前沿推进的过程中，无法保证待合并的三角形区域是一个单连通域，也就是说三角形区域可能被前沿分割成几个相互独立的区域，这就大大增加了产生残余三角形的概率.为了解决此问题，可以在前沿推进的过程中加入约束条件，以确保待合并的三角形始终是单连通域，从而避免了残余三角形的产生.首先引入几个定义.前沿边：位于四边形区域与三角形区域分界线上的边.前沿结点：位于前沿边上的结点.在以当前前沿(如图6所示的AB)生成四边形单元的过程中，可能出现的所有情况如图6所示.从图6中可以看出，在(a)，(d)，(f)，(i)和(j)所示的情形下，三角形区域在前沿AB生成四边形以后仍然是一个单连通区域，此时可以以AB为前沿生成四边形ABCD，除此之外，在生成四边形ABCD后，前沿会将三角形区域分割成多个块，从而可能会导致残余三角形的形成，此时不能以AB为前沿生成四边形ABCD.用nFrontNode表示ABCD中前沿结点的个数，nFrontEdge表示ABCD中前沿边的个数，则前沿推进的约束条件可以归纳如下：1)nFrontNode=2；2)nFrontNode=3，nFrontEdge=2；3)nFrontNode=4，nFrontEdge>2.因此，只有满足这三个条件之一，才可以生成四边形.此外，除了上述10种情况以外，还有一种特殊情况，如图7所示.在常规的前沿推进的过程中，若出现如图8(a)所示的情形时(粗线段为前沿)，以AB为前沿，在形成侧边时(如图8(b)所示)，若α1满足要求，则BD是一条侧边，生成的如图8(b)所示的四边形ABDE，更新前沿后会产生如图8(c)所示的残余三角形BCD.若在前沿推进的过程中，加入上述约束条件，则可以避免残余三角形的产生，因为生成图8(b)所示的四边形ABDE时，nFrontNode=4，nFrontEdge=2，不满足约束条件，此时不能生成四边形，然后循环前沿链表中的下一前沿BC，对于可能生成的四边形BCDE(如图9(a)所示)，nFrontNode=4，nFrontEdge=3，满足约束条件，则可生成四边形单元，更新前沿后如图9(c)所示，同理，再以EF为前沿时生成四边形EFAB，更新前沿后如图9(d)所示，显而易见，经过了带约束的前沿推进操作，生成了四边形ABCD和EFAB后，三角形区域仍是一个单连通域，没有残余三角形的产生.从图9可以看出，带约束条件的前沿推进不会产生残余三角形，从而避免了处理残余三角形这一复杂程序带来的时间消耗，以保证生产高质量的全四边形网格，提高了算法效率和网格质量.本文采用几何准则对生成的网格质量进行评价.如图10所示的四边形ABCD，连接其对角线AC，BD，得到ABC，ABD，CDB，CDA等4个三角形，先计算这4个三角形的质量因子，然后根据这4个三角形的质量因子，得到四边形ABCD的质量因子.其中三角形的质量因子α[9]定义如下：α的意义就是三角形的面积与边长的平方和之比，0<α≤1，α越大表示质量越好.令4个三角形△ABC，△ABD，△CDB，△CDA的α因子依次为α1，α2，α3，α4，假设α1≥α2≥α3≥α4，则四边形ABCD的质量因子β[14]定义如下：其中0<β≤1，β越大表示网格质量越好.为了说明本文提出算法的优越性，下面给出一个带小孔的平板的算例(因为对称性，在分析时取其1/4即可).并和Owen的算法[13]以及胡向红等人的算法[10]进行对比，如图11所示.从图11和表1中可以看出：Owen的算法[13]生成的四边形网格质量较好，但是有残余三角形的存在.胡向红等人的算法[10]生成的网格没有残余三角形，但是网格质量较差.本文提出的算法可以有效地解决上述两种算法的不足，生成高质量的全四边形网格.下面给出一个复杂多连通域实体的表面网格划分(如图12所示).从图12中可以看出，对于复杂的多连通域,本文的算法也可以生成高质量的四边形网格,从而进一步说明该算法的可行性.本文提出的算法是一种间接的四边形网格生成方法，该方法基于前沿推进的思想，在单个前沿推进的过程中，取合适侧边和顶边，控制当前前沿生成的网格质量，并在推进的过程中加入约束条件，避免了残余三角形的产生，数值算例表明该算法能够生成质量较好的四边形网格.相对于常规的前沿推进算法，本文方法可以有效避免残余三角形的产生，提高了算法的效率和网格生成的质量.相对于胡向红的区域生长算法[10]，本文方法可以有效控制当前四边形的质量，从而使最终生成的网格质量有显著的提高.†通讯联系人，E-mail:***************.cn【相关文献】[1] ZHANG Jianming，QIN Xianyun，HAN Xu.A boundary face method for potential problems in three dimensions[J].International Journal for Numerical Method inEngineering，2009，80(3)：320-337.[2] 张见明，李湘贺，陆陈俊，等.多域边界面法在稳态热传导问题中的应用[J].湖南大学学报：自然科学版，2014，41(5)：58-64.ZHANG Jianming，LI Xianghe，LU Chenjun，et al.Application of multi-domain boundary face method in steady state heat conduction[J].Journal of Hunan University：Natural Sciences，2014，41(5)：58-64.(In Chinese)[3] RUIZ-GIRONS E，SARRATE J.Generation of structured meshes in multiply connected surfaces using subsub-mapping[J].Advances in Engineering Software，2010，41(2)：379-387.[4] MUKHERJEE N.An art gallery approach to submap meshing[J].Procedia Engineering，2014，82(82)：313-324.[5] BLACKER T D，STEPHENSON M B.Paving：a new appro-ach to automated quadrilateral mesh generation[J].International Journal for Numerical Methods in Engineering，1991，32(4)：811-847.[6] 梅中义，范玉青.一种用于自动生成二维全四边形有限元网格的改进的铺设算法[J].计算机辅助设计与图形学学报，2000，12(6)：428-434.MEI Zhongyi，FAN Yuqing.A modified paving technique for quadrilateral mesh generation in planar regions[J].Journal of Computer-Aided Design & Computer Graphics，2000，12(6)：428-434.(In Chinese)[7] PELLENARD B，ORBAY G，CHEN J，et al.QMCF：QMorph cross field-driven quad-dominant meshing algorithm[J].Procedia Engineering，2014，82：338-350.[8] MERHOF D，GROSS R，TREMEL U，et al.Anisotropic quadrilateral meshing generation：anindirect approach[J].Engineering Computational Technology，2007，38(11/12)：860-867.[9] LEE C K，LO S H.A new scheme for the generation of a graded quadrilateralmesh[J].Computers & Structures，1994，52(5)：847-857.[10]胡向红，陈康宁.由区域生长算法实现四边形网格划分[J].计算机辅助设计与图形学学报，2004，16(1)：29-34.HU Xianghong，CHEN kangning.Generating quadrilateral meshes by using region growth algorithm[J].Journal of Computer Aided Design & Computer Graphics，2004，16(1)：29-34.(In Chinese)[11]LO S H.Generating quadrilateral elements on plane and over curvedsurfaces[J].Computers & Structures，1989，31(3)：421-426.[12]OWEN S J，STATEN M L，CANANN S A，et al.Q-Morph：an indirect approach to advancing front quad meshing[J].International Journal for Numerical Methods in Engineering，1999，44(9)：1317-1340.[13]OWEN S J，STATEN M L，CANANN S A.Advancing front quadrilateral meshing usingtriangle transformations[C]//Proceedings of the 9th International Meshing Roundtable.New Orleans，2000：409-428.[14]PARK C，NOH J S，JANG I S，et al.A new automated scheme of quadrilateral mesh generation for randomly distributed line constraints[J].Computer-Aided Design，2007，39(4)：258-267.。

基于区域分解技术的并行四面体网格生成算法

Ｐａｒａｌｌｅｌｔｅｔｒａｈｅｄｒａｌｍｅｓｈｇｅｎｅｒａｔｉｏｎａｌｇｏｒｉｔｈｍｂａｓｅｄｏｎｄｏｍａｉｎｄｅｃｏｍｐｏｓｉｔｉｏｎ
ＸＵＱｕａｎ，ＣＵＩＴａｏ。，ＬＩＵＱｉｎｇ－ｋａｉ，ＣＡＯＸｉａｏ — ｌｉｎ
２０１４年１月
计算机工程与设计
ＣＯＭＰＵＴＥＲＥＮＧＩＮＥＥＲＩＮＧＡＮＤＤＥＳＩＧＮ
Ｊａｎ．２０１４
Ｖｏ１．３５Ｎｏ．１
第３５卷
第１期
基于区域分解技术的并行四面体网格生成算法
徐权，崔涛，刘青凯，曹小林
法相比，不仅大大降低了网格生成的时间，同时保证了子区域之间交界面上网格的一致性和协调性，生成了高质量的四面
体网格。
关键词：网格生成；约束Ｄｅｌａｕｎａｙ三角化；并行四面体网格生成；区域分解；并行算法中图法分类号：ＴＰ３１１．１文献标识号：Ａ文章编号：１０００ — ７０２４（２０１４）０１ — ０１５３ — ０５
Ａｂｓｔｒａｃｔ：Ａｐａｒａｌｌｅｌｔｅｔｒａｈｅｄｒａｌｍｅｓｈｇｅｎｅｒａｔｉｏｎｍｅｔｈｏｄｂａｓｅｄｏｎｄｏｍａｉｎｄｅｃｏｍｐｏｓｉｔｉｏｎｉｓｐｒｅｓｅｎｔｅｄ．Ｉｔｃａｎｂｅａｐｐｌｉｅｄｆｏｒｔｈｅ３一ｄｉｍｅｎｓｉｏｎｃｏｍｐｌｅｘｇｅｏｍｅｔｒｉｅｓ．Ｂａｓｅｄｏｎｔｈｅｄｉｖｉｄｅａｎｄｃｏｎｑｕｅｒｐｒｉｎｃｉｐｌｅ，ｔｈｅ３Ｄｃｏｍｐｕｔａｔｉｏｎａｌｄｏｍａｉｎｉｓｄｉｖｉｄｅｄｉｎｔｏｍａｎｙｓｕｂ－ｄｏｍａｉｎｓ，ｗｈｉｃｈａｒｅｄｉｓｔｒｉｂｕｔｅｄｉｎｔｏｅａｃｈｐｒｏｃｅｓｓｏｒ．Ｆｉｎａｌｌｙ，ｆｏｒｓｕｂ－ｄｏｍａｉｎｓ，ｃｏｎｓｔｒａｉｎｅｄＤｅｌａｕｎａｙｔｒｉａｎｇｕｌａｔｉｏｎａｎｄｉｔｅｒａｔｉｖｅｔｅｃｈｎｏｌｏｇｙａｒｅｕｓｅｄｔｏｃｏｎｓｔｒｕｃｔｔｅｔｒａｈｅｄｒａｌｍｅｓｈｅｓｓｉｍｕｌｔａｎｅｏｕｓｌｙ．Ｔｈｅｎｕｍｅｒｉｃａｌｅｘｐｅｒｉｍｅｎｔｓｓｈｏｗｔｈｉｓａｌｇｏｒｉｔｈｍｉｓｓｃａｌａｂｌｅａｎｄｓｔａｂｌｅ，ａｎｄｉｔｃａｎｇｕａｒａｎｔｅｅｔｈｅｃｏｎｓｉｓｔｅｎｃｙｏｆｔｈｅｍｅｓｈｅｓｏｎｔｈｅｉｎｔｅｒｆａｃｅａｎｄｃａｎｇｅｎｅｒａｔｅｈｉｇｈｑｕａｌｉｔｙｔｅｔｒａｈｅｄｒａｌｍｅｓｈｅｓａｔ

四边形网格——精选推荐

四边形⽹格四边形⽹格定义四边形⽹格，顾名思义，它的每个⽹格⾯⽚是⼀个四边形。

有时候，四边形⽹格⾥会掺杂⼀些三⾓形⾯⽚，我们把这类⽹格也都叫做四边形⽹格。

三⾓形⽹格常见于逆向建模领域，⽐如通过三维扫描仪扫描得到的⽹格。

四边形⽹格常见于正向建模系统，如3dsMax，ZBrush等。

这主要是因为点云或者三⾓形⽹格转成四边形⽹格有⼀定的难度，特别是⾼质量的四边形⽹格。

正则点：内点-度数为4；边界点(⾮拐点(Corner))-度数为3；边界点(凸拐点)-度数为2；边界点(凹拐点)-度数为4分类：主要是根据顶点的正则度来进⾏分类。

如下图所⽰，第⼀类⽹格为正则⽹格，所有顶点度数为4，只有特殊拓扑结构的⽹格能达到正则；第⼆类是半正则⽹格，它是分⽚正则的；第三类是度数半正则⽹格，它的顶点度数绝⼤部分是4；最后⼀类是⽆序的四边形⽹格，它有很多⾮正则点。

四边形⽹格的优缺点与三⾓形⽹格相⽐，四边形⽹格有⼀些优点：特征边对齐：四边形⽹格的边可以很⾃然的与特征边进⾏对齐，边⾛向也可以很⾃然的与模型的⼏何特征⾛向对齐。

样条曲⾯和细分曲⾯：样条曲⾯和Catmull-Clark细分曲⾯常见的定义域就是四边形纹理贴图：半正则的四边形⽹格，每个正则⽚可以很好的与图⽚对齐，有利于图⽚的采样精确性同时四边形⽹格也有⼀些缺点：四边形可能不共⾯⾼质量的四边形⽹格奇异点个数尽量少，布局合理⾯⽚的边⾛向要与⼏何特征⾛向对齐边长尽量均匀化，或者⾃适应⼏何特征⽣成四边形⽹格的⽅法正向建模软件可以直接创建四边形⽹格。

逆向建模的⽹格⼀般是三⾓形⽹格，需要⽅法把三⾓形⽹格转成四边形⽹格Catmull-Clark细分三⾓形可以得到⼀个四边形⽹格，它的质量⽐较低，奇异点多，边⾛向不好。

后续可以应⽤⼀些四边形⽹格优化的⽅法参数化⽅法Morse-Smale complex ⽅法⽹格分割后，分⽚参数化的⽅法四边形⽹格的处理⽹格简化⼏何优化⾼阶曲⾯拟合，细分曲⾯拟合⽹格拓扑优化有兴趣的读者，欢迎参考。

平面任意区域四边形网格自动生成的一种方法

平面任意区域四边形网格自动生成的一种方法
顾元宪;马正阳
【期刊名称】《计算机辅助设计与图形学学报》
【年(卷),期】1998(010)005
【摘要】在改进节点连接法的基础上，提出了一种平面任意区域的有限元网格全自动剖分方法，既能快速生成四边形单元网格，也能生成三角形单元网格；【总页数】8页(P432-439)
【作者】顾元宪;马正阳
【作者单位】大连理工大学工程力学系;大连理工大学工程力学系
【正文语种】中文
【中图分类】TP391.72
【相关文献】
1.一种平面区域四边形网格生成方法 [J], 但文蛟;张平;胡树根;王耘
2.一种平面区域四边形网格生成方法 [J], 但文蛟;张平;胡树根;王耘
3.任意平面区域三角形网格的全自动生成算法 [J], 张修忠;金生
4.任意平面区域的变密度四边形网格生成方法 [J], 杨晓东;申长雨;陈静波;刘春太
5.任意四边形网格上解辐射输运方程的一种数值方法 [J], 冯庭桂
因版权原因，仅展示原文概要，查看原文内容请购买。

平面四边形网格自动生成方法研究

第4类，封闭单元操作，在铺砌过程的最后阶段只剩下6个节点时，根据铺砌边界上大鱼150º 的节点内角数目和相对位臵，产生封闭单元。
网格中单元的大小由铺砌边界上的节点的空间大小所决定。在铺砌过程中，维持这一空间大小不变。这样，我们可以通过修改固定节点的空间大小来控制单元的大小。各个节点的空间大小可以不相同，因为大小不一的单元组成的网格，有利于分析单元，便于网格的疏密分布。
②平面四边形网格的生成方法有两类主要的方法。一类是间接法，即在区域内部先生成三角形网格，然后分别将两个相邻的三角形合并成为一个四边形。生成的四边形的内角很难保证接近直角。所以再采用一些相应的修正方法加以修正。间接法优点是首先就得到了区域内的整体的网格尺寸的信息，对四边形网格尺寸梯度的控制一直是四边形网格生成技术的难点。缺点是生成的网格质量相对比较差，需要多次的修正，同时需要首先生成三角形网格，生成的速度也比较慢，程序的工作量大。另外一类是直接法，二维的情况称为铺砖法。采用从区域的边界到区域的内部逐层剖分的方法。这种方法到现在已经逐渐替代间接法而称为四边形网格的主要生成方法。它的优点是生成的四边形的网格质量好，对区域边界的拟合比较好，最适合流体力学的计算。缺点是生成的速度慢，程序设计复杂。
3、边界节点数的调整
第一种方法：将节点1、2合并成节点4，节点4的位臵取节点1、2、3的平均位臵。
第二种方法；删除单元a以后，将节点1、2合并成节点4，节点4的位臵取节点1、2、3的平均位臵。
以上两种方法每次使栅格边界减少2个节点，如果要减少一个可删除单元1,、2、3、4以后.，合并相应的节点。
基区划分任意一个几何轮廓若干个子区域直线段子区域曲线段子区域四边形基区类四边形基区类四边形基区是指在一定条件下可以把曲线段近似看作一条直线来处理把它看成是四边形的一条基区网格生成21四边形基区的网格生成由于基区形状与所要生成的单元具有相似性这种相似性的程度越高就越容易处理

一种城市路网多层次复合网格模式识别方法

㊀㊀第５２卷㊀第１１期测㊀绘㊀学㊀报V o l．５２,N o．１１㊀２０２３年１１月A c t aG e o d a e t i c ae tC a r t o g r a p h i c aS i n i c a N o v e m b e r,２０２３引文格式:王安东,武芳,巩现勇,等．一种城市路网多层次复合网格模式识别方法[J]．测绘学报,２０２３,５２(１１):１９９４Ｇ２００６．D O I:１０．１１９４７/j．A G C S．２０２３．２０２２０５２８．WA N G A n d o n g,WUF a n g,G O N GX i a n y o n g,e t a l．Ar e c o g n i t i o n a p p r o a c h f o r c o m p o u n d g r i d p a t t e r no f u r b a n r o a d n e t w o r k s [J]．A c t aG e o d a e t i c a e tC a r t o g r a p h i c aS i n i c a,２０２３,５２(１１):１９９４Ｇ２００６．D O I:１０．１１９４７/j．A G C S．２０２３．２０２２０５２８．一种城市路网多层次复合网格模式识别方法王安东,武㊀芳,巩现勇,翟仁健,刘呈熠,邱㊀越,张寒雪信息工程大学地理空间信息学院,河南郑州４５０００１A r e c o g n i t i o na p p r o a c h f o r c o m p o u n d g r i d p a t t e r no f u r b a n r o a dn e t w o r k s W A N G A n d o n g,W U F a n g,G O N G X i a n y o n g,Z H A IR e n j i a n,L I U C h e n g y i,Q I U Y u e,Z H A N GH a n x u eI n s t i t u t eo fG e o s p a t i a l I n f o r m a t i o n,I n f o r m a t i o nE n g i n e e r i n g U n i v e r s i t y,Z h e n g z h o u４５０００１,C h i n aA b s t r a c t:A s t h es k e l e t o no fu r b a nc i t i e s,t h es p a t i a l p a t t e r nr e c o g n i t i o no f r o a dn e t w o r k s i so f g r e a t s i g n i f i c a n c ef o r m a p g e n e r a l i z a t i o n,s p a t i a ld a t a m i n i n g,a n d m u l t iＧs c a l er e p r e s e n t a t i o n．T h i s p a p e r p r e s e n t sa na p p r o a c h t o r e c o g n i z i n g t h ec o m p o u n d g r i d p a t t e r no f r o a dn e t w o r k sw i t h l o c a l h e t e r o g e n e i t y b a s e do n r o a dm e s h e s．F i r s t l y,t h em u l t i l e v e l c o g n i t i v e c h a r a c t e r i s t i c s o f t h e l i n e a r a n d g r i d p a t t e r n o f r o a d m e s h e sa r ea n a l y z e d,a n dt h e m u l t i l e v e lc o g n i t i v eo r d e r,w h i c hf r o m b a s i c m e s h,c o m p o u n d m e s ht o r e g u l a r p a t t e r n,i s p r o p o s e d．S e c o n d l y,t h e r e c o g n i t i o n m e t h o d s o fi n c l u s i o n r e l a t i o n s h i p,p a r a l l e l r e l a t i o n s h i p,a n dl i n e a r p a t t e r n b e t w e e nr o a d m e s h e s a r e d e s i g n e d c o n s i d e r i n g t h e c o m p o s a b i l i t y, l i n e a r i t y,a n de x t e n s i b i l i t y o fc o m p o u n dl i n e a r p a t t e r n．F i n a l l y,t h el i n e a r p a t t e r n sa r ec o m b i n e d a n d d e c o m p o s e d t oe x t r a c t t h ec o m p o u n d g r i d p a t t e r no f r o a d m e s h e s．E x p e r i m e n t ss h o wt h a t t h e p r o p o s e d m e t h o d i s e f f e c t i v e f o r c o m p o u n d g r i d p a t t e r n r e c o g n i t i o nw i t h t h ea g r e e m e n t so f h u m a n s p a t i a l c o g n i t i v e c h a r a c t e r i s t i c s．K e y w o r d s:c a r t o g r a p h i c g e n e r a l i z a t i o n;r o a d n e t w o r k;r o a d m e s h e s;p a t t e r nr e c o g n i t i o n;m u l t i l e v e l c o g n i t i o n;g r i d p a t t e r n摘㊀要:道路网作为城市骨架,其模式识别对于地图综合㊁空间数据挖掘与多尺度表达具有重要意义.针对大比例尺数据中局部异质性明显的道路网格模式识别问题,提出基于网眼的城市道路多层次复合网格模式识别方法.首先分析了道路网眼直线和网格模式的多层次认知特点,提出了基础网眼➝复合网眼➝规则模式的多层次认知顺序;然后考虑复合直线模式的组合性㊁延伸性和直线性约束,设计了道路网眼直线模式㊁包含关系和并列关系的识别方法;最后通过对直线模式的组合分解,提取道路网眼的网格模式.试验表明本文方法能有效识别路网数据中的复合网格模式,识别结果符合人类认知特点.关键词:制图综合;道路网;网眼;模式识别;多层次认知;网格模式中图分类号:P２０８㊀㊀㊀㊀文献标识码:A㊀㊀㊀㊀文章编号:１００１Ｇ１５９５(２０２３)１１Ｇ１９９４Ｇ１３㊀㊀城市道路网是城市范围内不同功能㊁等级㊁区位的道路,以一定密度和适当形式组成的网络结构[１].作为城市的基础骨架,其结构模式体现了城市的主要结构和空间格局,反映出城市的地形地貌特点㊁功能结构和规划治理情况,蕴含着大量城市形成和发展的内在机制[２].对其结构模式的挖掘和识别是地图综合㊁城市形成㊁更新和扩张㊁交通规划设计等领域的研究热点和难点[３Ｇ５].相关研究从不同研究重点出发,将道路网分为不同结构模式,如网格模式㊁环型模式㊁放射型模式㊁复杂道路交叉口等显式模式[６Ｇ１１],以及城市中心㊁热点区域㊁城市建成区等隐式模式[１２].网格模式作为城市道路网的典型结构模式,在城市布局中十分常见,在长达几千年的城市发展史上,都得到广泛的采用[１３].道路网格模式的识别对城市空间特征挖掘㊁交通规划及地图自动综合具第１１期王安东,等:一种城市路网多层次复合网格模式识别方法有重要意义.根据识别的基本模式单元,现有道路网格模式的识别方法可大体分为两类.(１)基于路段的识别方法.此类方法多将道路网抽象为图结构,以图中顶点作为基本处理单元,从顶点的几何㊁上下文关系特征中抽象出与路网结构相关的特征项,借助图论㊁统计学或机器学习方法进行处理或学习,实现网格模式识别.例如,基于道路结点和改进的霍夫变换策略来实现规则格网的识别[１４];通过构建道路网对偶图,采用交㊁并㊁联合等图运算来提取基础格网模式[１５];基于道路结点的几何㊁拓扑特征,利用多项式评定模型识别道路网中的典型结构模式[１６];基于道路网的线性单元剖分,提出５种特征参量,采用支持向量机分类来提取网格模式[１７Ｇ１８];在构建道路网原始图的基础上,利用图卷积神经网络模型,通过学习人工标注样本,实现网格模式的识别[１９]等. (２)基于网眼的识别方法.此类方法将道路网中路段围成的闭合区域转化为面,即道路网眼,通过计算网眼与邻近网眼的形状㊁方向㊁尺寸相似性及排列特征,采用邻近搜索㊁任务分类㊁隶属度计算或自组织映射聚类等方法完成网格模式的识别.例如,通过计算网眼的几何特征相似性,采用区域生长算法识别网格模式[２０Ｇ２１];基于网眼的形状和关系描述参量,采用机器学习算法识别网格模式,以减少参数阈值设置的人工干预[２２Ｇ２５],等.然而,当前研究中至少存在如下问题有待解决:①根据定义,道路网格模式的基本特征是由两组几乎平行的道路垂直相交构成,网眼形状大多为近似矩形或平行四边形.然而,基于路段的识别方法大多以道路网结构中的正交性原则为依据,设计网格模式的特征因子,忽略了网格模式中相邻网格间的尺寸㊁形状相似性和分布的延伸性.从结构模式的定义来看,这些研究的部分识别结果更接近于道路网的正交模式或方格模式 ,而非网格模式 .②基于网眼的已有识别方法均以单个网眼面要素及相邻网眼间一对一的邻近关系作为研究对象,对于整体规则㊁局部破碎的网眼群组,无法将其作为一个整体参与模式构建.为了解决具有局部异质性的道路网格模式识别问题,本文基于道路网眼,顾及视觉认知的层次性,提出一种城市道路网多层次复合网格模式的识别方法,主要解决两个问题:①道路网眼分布模式的多层次认知特征和定义;②网眼复合直线模式和复合网格模式的识别.１㊀道路网眼分布模式的多层次认知特点道路网眼是指道路网中路段围成的闭合区域.与其他面状地图要素相似,道路网眼群组具有丰富的分布模式,直线模式和网格模式是两种典型的分布模式.直线模式中道路网眼有规律地呈直线分布,网眼间具有相似的几何特征,如方向㊁尺寸等;网格模式由若干组近似平行的直线模式与另外若干组近似平行的直线模式,以近似正交的方式相交构成.作为道路对空间划分的结果,其分布模式与道路结构模式有着密切联系.道路网眼的直线和网格模式是道路网格模式的两种表现形式,图１中模式１㊁２分别为网眼的直线和网格模式,两组模式均表现为道路网的网格模式.因此,本文以道路网眼为基本模式单元,通过提取其直线和网格模式,实现道路网格模式的识别.图１㊀道路网结构模式与网眼分布模式F i g．１㊀T h es t r u c t u r a l p a t t e r no f r o a dn e t w o r k sa n dt h ed i s t r i b u t i o n p a t te r no fm e s h依据格式塔认知准则[２６]和大范围优先的视知觉认知理论[２７Ｇ２８],人类更倾向于以主体ң细节的顺序来认知事物.道路网眼群组具有丰富复杂的几何㊁拓扑特征,人类的认知过程也必然遵循一定的顺序,从而形成了其空间关系的层次性.观察者观看地图时,首先关注道路网眼的整体特征,如图２(a)红色方框中网眼整体具有明显的网格模式特征;然后,才会注意到局部网眼间的细节特征,如图２(b)所示,蓝色网眼间的几何形态差异和复杂拓扑关系会被进一步感知.然而,当前相关研究仅考虑相邻单个网眼间的一对一关系,难以识别由多个不规则的网眼多边形拼接而成的网格网眼,对于大比例尺地图中局部异质性明显的道路网格模式,识别结果并不符合人类认知[２９].５９９１N o v e m b e r ２０２３V o l ．５２N o ．１１A G C Sh t t p :ʊx b ．c h i n a s m p ．c om 图２㊀整体到局部的认知过程F i g ．２㊀T h e c o gn i t i o n p r o c e s s f r o m w h o l e t o p a r t ㊀㊀为解决网格模式中局部网眼破碎的问题,本文引入复合道路网眼的概念.复合道路网眼的直线和网格模式具有多层次认知的特点:宏观尺度下,道路网眼整体呈直线或网格模式分布;中观尺度下,模式由几何特征相似㊁排列规律相近的简单或复合矩形网眼构成;微观尺度下,根据简单网眼的组合方式,复合矩形网眼进一步划分为包含关系和并列关系复合网眼.复合道路网眼直线和网格模式的多层次认知关系如图３所示.以图４(a)中道路网为例,各层次的具体含义如下.图３㊀道路网眼多层次认知关系F i g ．３㊀T h e r e l a t i o n s h i p o fm u l t i l e v e l c o gn i t i on 图４㊀多层次认知过程F i g ．４㊀T h em u l t i l e v e l c o gn i t i v e p r o c e s s ㊀㊀(１)整体层,包括直线模式和网格模式,其中直线模式的识别是网格模式识别的前提和保障.网眼直线模式具有以下的表现形式:①模式内部的各个网眼具有相似的形状㊁大小和方向特征;②模式内相邻网眼方向一致,且模式的全局方向与各网眼组件方向近似相同或正交;③模式内相邻网眼的公共边近似为两最小面积外界矩形的长边或短边,图４(b )中红色虚线分别表示呈直线模式分布的网眼.网格模式由多组直线模式近似垂直相交构成,处于更高的认知层次,如图４(b )中９组直线模式以近似正交的方式相交构成网格模式.(２)组件层,构成单元为近似矩形的道路网眼.根据矩形网眼所中的基础网眼数量,分为简单矩形网眼和复合矩形网眼,分别如图４(c )中网眼３㊁９和网眼１㊁２㊁４㊁５㊁６㊁７㊁８.其中,简单矩形网眼形状为近似矩形;复合矩形网眼包含多个任意形状的道路网眼,组合后形状为近似矩形.从矩形网眼与基础网眼的空间对应关系的角度来看,组件层中包含１ʒ１(简单网眼)和１ʒn (复合网眼)的空间对应关系.(３)原子层,构成单元为由道路网结点㊁路段直接围成的封闭区域多边形,即简单网眼.根据邻接关系,将组件层复合网眼中简单网眼间关系划分为包含关系和并列关系.包含关系复合网眼由一个主体网眼和若干次要网眼组成,如图４(d)中灰色网眼５㊁６,网眼８㊁９,网眼１０㊁１１和网眼１７㊁１８.主次网眼间空间邻近,整体轮廓互补.其中,主体网眼的面积相对较大,在视觉认知中占主导地位,反映该复合网眼的主要形状特征,如图４(d )中网眼６㊁８㊁１０和１８;次要网眼的面积相对较小,在视觉认知中占从属地位,如图４(d )中网眼５㊁９㊁１１和１７.并列关系复合网眼由若干简单网眼组成,网眼为任意形状多边形,组合后形状为近似矩形,并与相邻网眼构成直线模式.网眼间并列关系难以通过自底向上的组合方法进行探测,６９９１第１１期王安东,等:一种城市路网多层次复合网格模式识别方法其关系的识别依赖于复合网眼邻域的模式特征.图４(d)中蓝色网眼１㊁２㊁３㊁４,网眼１２㊁１３和网眼１４㊁１５㊁１６分别为具有并列关系的简单网眼.２㊀道路网多层次网格模式识别结合道路网眼的多层次认知特点,本文采用自底向上与自顶向下相结合的策略,提出一种多层次道路网眼直线和网格模式的识别方法,整体框架如图５所示,基本思想和关键步骤如下.(１)根据视知觉感知理论中完整性㊁规则性等心理倾向,采用自底向上的策略,合并具有包含关系的网眼,将整体规则㊁局部不规则的相邻网眼组合为视觉感知上更高级的复合网眼.(２)考虑相邻网眼间的尺寸㊁形状相似性和分布的直线性,构建直线模式结构化参数,以此为约束提取直线模式.(３)根据直线模式分布的延伸性,自顶向下构建直线模式连续匹配模板,搜索合并直线模式两端具有并列关系的道路网眼,实现复合直线模式的提取.(４)采用降维的思想,将直线模式网眼组以二维的线段表示,对其进行分解和组合,实现道路网眼网格模式的提取.图５㊀本文方法整体框架F i g．５㊀T h e f r a m e w o r ko f t h e p r o p o s e dm e t h o d道路网眼形状为近似矩形是其作为直线和网格模式的组成单元的必要条件[２０Ｇ２４].本文采用矩形度(R e c)和凹凸度(C o n v)[２９]作为网眼矩形相似度的度量参数,具体含义与计算方法见表１.表１㊀网眼矩形相似度的度量参数T a b．１㊀T h em e t r i c s o f r o a dm e s h r e c t a n g l e s i m i l a r i t y参数参数含义计算方法矩形度R e c描述多边形呈矩形的程度网眼自身面积与其最小面积外接矩形面积的比值凹凸度C o n v描述道路网眼多边形的凹凸程度网眼面积与其凸包面积的比值２．１㊀简单直线模式识别简单直线模式由当前数据中的简单网眼构成,其正确识别是复合直线模式识别的基础.考虑直线模式中网眼的相似性㊁直线性和延伸性等结构特征,结合格式塔认知准则,从道路网眼的大小相似性㊁直线性和对齐程度３个方面引入识别道路网眼直线模式的结构化参数(表２).表２㊀网眼直线模式结构化参数T a b．２㊀T h e s t r u c t u r a l p a r a m e t e r s o f r o a dm e s h l i n e a r p a t t e r n参数参数含义计算方法面积比R a r e a描述两邻接道路网眼之间面积大小差异相邻网眼间较小网眼与较大网眼面积的比值方向差异D o r i e n t描述３个相邻道路网眼之间的直线性相邻３个网眼间质心连线的夹角公共边长度比R c e描述相邻道路网眼之间的对齐程度相邻网眼公共边与网眼最小外接矩形边长度比值的较小值㊀㊀根据直线模式的组织规律,同一模式内部要素间具有相似的几何形态结构.对于道路网眼,其几何形态结构主要表现为形状和大小.其中形状依靠上文中矩形相似度参数进行约束;大小相似度利用网眼面积比进行度量,面积比越大,网眼间大小相似程度越高.为保证模式中的网眼沿直线分布,以表２中方向差异参数作为约束,方向差异D o r i e n t越接近于１８０ʎ,模式的直线性越强.考虑到直线模式相邻网眼间具有相互对齐的特点,引入公共边长度比R c e对模式内相邻网眼间对齐程度进行约束.如图６所示,对于相邻网眼M１㊁M２,定义其公共边长度比R(c e)１,２为公共边(P１P２)长度与公共边对应网眼最小外接矩形边(E M１㊁E M２)长度的比值的较小值,公共边长度比７９９１N o v e m b e r ２０２３V o l ．５２N o ．１１A G C Sh t t p :ʊx b ．c h i n a s m p ．c o m 越接近于１,网眼对齐程度越高.图６㊀公共边长度比F i g ．６㊀T h e l e n g t h r a t i oo f c o mm o ne d ge 综上,本文识别网眼简单直线模式的步骤如下.步骤１:根据网眼间是否具有公共边,构建相邻网眼间邻近关系,同时提取网眼群中矩形度和凹凸度分别大于阈值δR e c 和δC o n v 的简单矩形网眼,加入列表L i s t S G M .步骤２:选取L i s t S G M 中任一网眼M i 及其邻接网眼M j ,计算M i 与M j 的面积比(R a r e a )i ,j 和公共边长度比(R c e )i ,j ,若(R a r e a )i ,j ＞δa r e a 且(R c e )i ,j ＞δc e (式中δa r e a 和δc e 分别为人工设定的面积比和公共边长度比参数阈值),则将网眼M i 与M j 的邻近边e i ,j 加入直线模式临时列表t L i s t L P ,否则,返回步骤２.步骤３:设当前直线模式一侧搜索方向为{i ,j },即S e a r c h L e f t ＝{i ,j },另一侧搜索方向为{j ,i },即S e a r c h R i gh t ＝{j ,i },以S e a r c h L e f t ＝{i ,j }为起始搜索方向,选取M j 的邻接网眼M k ,若M j 除M i 不存在其他邻近网眼,则此搜索方向终止,以S e a r c h R i gh t ＝{j ,i }方向继续搜索.步骤４:计算M j 与M k 的面积比(R a r e a )j ,k 和公共边长度比(R c e )j ,k ,若(R a r e a )j ,k ＞δa r e a 且(R c e )j ,k ＞δc e ,则执行步骤５,否则,返回步骤３.步骤５:计算M i ㊁M j ㊁M k 的方向差异D o r i e n t ,若D o r i e n t ＞δo r i e n t ,则将e j ,k 添加至当前直线模式列表t L i s t L P 中,并令j ＝k ,否则,返回步骤３.步骤６:若当前直线模式向两侧搜索均终止,则该组直线模式识别结束,将当前直线模式列表t L i s t L P 加入直线模式识别结果列表L i s t L P 中,返回步骤２.循环步骤２６,直至L i s t S G M 中全部网眼均被遍历.以图７(a)中道路网数据为例,经上述步骤识别的网眼简单直线模式如图７(b )中红色线段所示.受道路网中较低等级道路影响,一些在大尺度上认知为整体的网眼被分割成若干小网眼,呈现出局部破碎的现象,如图７(b )中网眼１㊁２㊁３.由图７(c )中网眼邻近关系可以看出,简单直线模式的提取方法仅利用简单网眼间一对一的邻近关系(蓝色线段),难以反映道路网的整体结构模式,需要进一步利用复合直线模式识别方法,提取局部破碎㊁整体规则的直线模式.图７㊀道路网眼简单直线模式识别结果F i g ．７㊀T h e r e c o g n i t i o n r e s u l t s o f s i m pl e l i n e a r p a t t e r n ２．２㊀复合网眼构建与复合直线模式识别复合直线模式的识别是解决网眼局部异质性,实现由低级基础网眼到高级认知模式过渡的关键,其难点在于组件层中复合矩形网眼的识别和构建.对于复合矩形网眼中的包含关系和并列关系,本文分别采用自底向上和自顶向下的策略对其进行识别.２．２．１㊀包含关系识别根据包含关系网眼间轮廓互补的特点,参考文献[３０]中对相离面要素主次关系识别的方法,引入公共边长周长比(R l e n g t h )和约束面积比(R c a )两个参数,分别从一维和二维两个维度描述相邻网眼之间的包含程度.以图８中道路网眼为例,红色线段P ２P ３P ４表示相邻网眼M １㊁M ２间公共边,虚线矩形为M ２最小面积外接矩形S M B R M ２,灰色多边形P １P ２P ３P ４为S M B R M ２与M １的公共区域多边形,参数含义及计算方法见表３.８９９１第１１期王安东,等:一种城市路网多层次复合网格模式识别方法图８㊀包含关系参数F i g ．８㊀T h e i n c l u s i o n r e l a t i o n s h i ppa r a m e t e r 表３㊀包含关系识别参数T a b ．３㊀T h e r e c o g n i t i o n p a r a m e t e r s o f i n c l u s i o n r e l a t i o n s h i p参数参数含义计算方法公共边周长比R L 描述两相邻网眼边界的包含程度相邻网眼公共边长度与网眼周长的比值约束面积比R c a描述两相邻网眼区域的包含程度网眼与相邻网眼最小面积外接矩形交集面积与网眼面积的比值根据包含关系参数和网眼矩形相似度参数,识别包含关系复合矩形网眼的步骤如下.步骤１:计算网眼矩形度R e c 和凹凸度C o n v ,将R e c ＜δR e c 或C o n v ＜δC o n v 的网眼加入列表m L i s t 中.从中选取网眼M i ,计算其与邻近网眼M j 的公共边周长比(R l e n g t h )i ,j 和约束面积比(R c a )i ,j ,若(R l e n g t h )i ,j ＞δL 且(R c a )i ,j ＞δc a ,则将网眼M i 和M j 记为包含关系组,其中δl e n g t h 和δc a 分别为人工设定的公共边周长比和约束面积比参数阈值.步骤２:合并网眼M i 与M j ,记新网眼为M n ,若R e c n ＞δR e c 或C o n v n ＞C o n v i ,说明次要网眼对主要网眼的规则程度具有补充作用,将M n 加入列表m L i s t 中,并从中删除M i 和M j ,否则删除M n .步骤３:循环步骤１㊁２,直至列表m L i s t 中的元素数量不再减少为止,此时全部具有包含关系的复合矩形网眼均被识别.以图７中道路网眼为例,利用上述步骤识别㊁合并包含关系网眼的过程如图９所示.图９(a )为步骤１包含关系的识别结果,粉色线段表示网眼间的包含关系.图９(b )为步骤２包含关系第一次合并结果(深色网眼),粉色线段表示所产生的新的包含关系.经数轮迭代,具有包含关系网眼的最终合并结果如图９(c )所示,其中深色网眼为合并后的包含关系复合网眼.图９㊀包含关系复合网眼识别与合并过程F i g ．９㊀T h e r e c o g n i t i o na n d c o m b i n a t i o n p r o c e s s o f c o m p o u n dm e s h e sw i t h i n c l u s i o n r e l a t i o n s h i p２．２．２㊀并列关系识别及复合直线模式提取从认知角度来看,并列关系复合网眼产生于邻域内直线模式的延伸,例如对于图９(c )中邻近关系相似的网眼对１㊁２和２㊁３,网眼１㊁２更容易被组合为复合网眼.因此,本文利用网眼邻域的模式特征,采用自顶向下的匹配策略,识别具有并列关系的复合网眼.基本思想为:首先根据直线模式两端网眼确定初始匹配模板的几何特征;然后结合直线模式延伸性构建连续匹配模板,向两端搜索㊁匹配待识别网眼;最后结合直线模式约束条件,判定待识别网眼组合后能否构成直线模式,实现复合直线模式的提取.结合图９(a)中道路网眼,说明并列关系及复合直线模式识别方法的具体步骤.步骤１:合并具有包含关系的矩形网眼(图９(c )),提取简单直线模式(图１０(a)中红色线段),存入S L P L i s t 中.步骤２:选取任意一组直线模式S L P i (图１０(a)中红色加粗线段),提取其首㊁末端网眼的最小面积外接矩形,以首㊁末端网眼与其邻接网眼几何中心的相对距离d ㊁方向o 为约束,沿直线模式两端延伸方向计算匹配模板的位置,匹配模板分别记为T l ㊁T r (图１０(a)中蓝色矩形).步骤３:搜索与T l 和T r 存在面状交集的网９９９１N o v e m b e r ２０２３V o l ．５２N o ．１１A G C Sh t t p :ʊx b ．c h i n a s m p ．c o m 眼,记为M s ,如果其模板重叠度(R t o )s ,l ＞δt o 或(R t o )s ,r ＞δt o ,则将M s 存入并列关系候选列表M L i s t 中,执行步骤４,否则,则执行步骤２,式中δT C 为人为设定的模板重叠度参数阈值.步骤４:若M l i s t 中网眼数量大于１,合并M L i s t 中全部网眼,记为M n ,如图１０(b )中蓝色网眼,若R e c (M n )＞δR e c ,则根据２．１节中方法,判断其能否满足直线模式结构化参数约束,若满足,执行步骤５,否则终止该侧搜索.步骤５:当直线模式S L P i 向两侧搜索均终止时,将其存入复合直线模式列表M L P L i s t,并从S L P L i s t 中移除,执行步骤２.循环步骤１５,直至S L P L i s t 为空时,结束循环.识别结果如图１０(c)所示,其中红色线段表示网眼直线模式.图１０㊀并列关系网眼与直线模式识别过程F i g ．１０㊀T h e r e c o gn i t i o n p r o c e s s o f p a r a l l e l r e l a t i o nm e s h e s a n d l i n e a r p a t t e r n ２．３㊀网格模式提取道路网眼的网格模式由近似正交的直线模式相交构成,处于更高的认知层次.由网格模式概念可知,组成网格模式的直线模式之间需满足以下３项条件:①各组直线模式近似平行或正交;②正交的直线模式间具有相交关系;③各直线模式构成闭合回路.对于条件①,由于在直线模式中,网眼为方向一致的近似矩形,网眼构成的直线模式方向基本确定,若任意两组直线模式包含同一网眼,则其关系为近似正交;若任意两组不相交直线模式间,存在其他直线模式同时包含以上两组直线模式中的网眼,则两组直线模式近似平行.故条件①可由条件②代替.另外,若多组直线模式构成闭合回路,则相互正交的直线模式间必然相交,故条件③为条件②的充分条件.综上,本文通过对直线模式网眼构成闭合回路进行识别,提取其中的网格模式.当前研究大多采用图论中算法识别多组直线模式中的闭合回路[３１Ｇ３２],算法实现较为复杂.本文从几何角度出发,通过对直线模式邻近图中由结点和线段形成的封闭多边形进行聚类,实现网格模式的提取.以图１１(a )道路数据为例,说明算法具体步骤.步骤１:以各直线模式中网眼几何中心为结点,依次连接相邻网眼,构建直线模式关系图G l(图１１(b)中红色线段).步骤２:提取G l 中封闭区域多边形P (图１１(b )中晕线多边形P １㊁P ２㊁P ３),将其存入列表T P L i s t.步骤３:遍历T P L i s t ,若两多边形P i ㊁P j 间存在公共边(图１１(b )中多边形P ２㊁P ３),则将其聚类为一组,存入列表C l u s t e r L i s t 中.图１１㊀网格模式提取过程F i g．１１㊀T h e e x t r a c t i o n p r o c e s s o f g r i d p a t t e r n ０００２第１１期王安东,等:一种城市路网多层次复合网格模式识别方法㊀㊀步骤４:根据聚类结果,提取每组聚类内多边形顶点所对应的道路网眼,记为一组网格模式网眼G r i d m ,存入G r i d _l i s t 中,如图１１(c )中蓝色和黄色线段对应网眼分别为两组网格模式.３㊀试验与讨论３．１㊀试㊀验本文基于P y t h o n 和Q G I S 编程实现以上算法.试验分为两组,数据分别采用不同空间结构模式的国内外道路网数据.试验中网眼矩形度和凹凸度阈值的设置影响简单矩形网眼以及包含关系㊁并列关系的识别,结合前人的研究成果[２９,３３],取δR e c ＝０．９㊁δC o n v ＝０．９５.算法中所涉及其他阈值设置的指导思想为:面积比和公共边长度比的阈值越大,方向差异的阈值越小,模式的直线性越强.在识别包含关系矩形网眼时,公共边周长比和约束面积比的阈值越大,主要网眼对次要网眼的包含程度越强;在识别并列关系矩形网眼时,模板重合度的阈值越大,模式的同质性越强.试验区域１为国外某地区道路网,数据来自O p e n S t r e e t M a p.如图１２所示,所选试验区域内道路分布密度较高,密度不均匀,没有全局平稳的特征,存在支离破碎的短小路段,但大部分呈垂直交错结构分布,视知觉上呈现出明显的网格模式特征.试验区域包含２４６９个道路简单网眼.经反复测试,试验参数设置见表４.图１２(b )中灰色和蓝色网眼分别为识别出的具有包含和并列关系的复合网眼,提取的直线模式和网格模式分别如图１２(c )㊁(d)所示.图１２㊀试验１直线模式和网格模式识别结果F i g ．１２㊀T h e r e c o g n i t i o n r e s u l t s o f l i n e a r p a t t e r na n d g r i d p a t t e r no f e x pe r i m e n t １表４㊀参数阈值设置T a b ．４㊀A d v i s a b l e p a r a m e t e r s e t t i n g试验简单直线模式识别包含关系识别并列关系识别面积比阈值方向差阈值公共边长度比阈值公共边周长比阈值约束面积比阈值模板重合度阈值试验１０．５１０ʎ０．９０．４０．５０．８试验２０．４１０ʎ０．８０．４０．５０．７１００２。

基于区域递归分解算法的四边形网格生成技术的研究

第28卷　第4期2004年8月武汉理工大学学报(交通科学与工程版)Journal of Wuhan University of Technology(Tr ansportat ion Science&Engineer ing)Vol.28　No.4Aug.2004基于区域递归分解算法的四边形网格生成技术的研究*¹刘　阳　张　武　唐锦春(浙江大学建筑工程学院　杭州　310027)摘要:基于区域递归分解算法是一种有效的四边形网格生成算法,能较好地运用于复杂的单连通和多连通图形区域,自动化程度高.文中修正Sar ra te提出的算法,描述程序实现的具体过程.原算法对六节点子区域的剖分效果不理想,文中就此提出了区分和剖分各类六节点子区域的方法,提高了算法效率,并给出网格生成实例.关键词:四边形;网格划分;区域递归分解中图法分类号:T P399;T V131.4 网格生成是有限元前处理阶段的关键一步,网格生成的质量直接影响有限元分析结果的精度.三角形网格的计算性能不如四边形网格,因此四边形网格生成技术就成为当前研究的热点.基于区域递归分解的网格生成算法最早由Sluiter[1]提出,T albert[1]和Sarr ate[2]分别作出了改进.文中即采用Sarrate提出的算法,对其作出修正,并提出一种区分六节点子区域类型的有效方法,给出程序实现的具体过程.该算法具有自动化程度高,网格疏密控制较好,能较好地适用于复杂的图形区域和多连通区域的特点.1　算法描述运用区域递归分解法,首先在初始边界线段的端点处指定单元尺寸信息,根据单元尺寸信息,在线段内部渐变的布点,从而形成初始的边界节点,将初始边界转化为由边界节点形成的多边形区域.连接两个边界节点,形成剖分线,在剖分线上布置新节点,把初始区域剖分成两个子区域.重复这个过程,对子区域不断进行剖分,直至最后形成四边形网格.区域递归分解算法的关键在于剖分线的选取,新节点生成等方面.1.1　剖分线的选取为得到质量较好的网格,应选取该区域内最优的剖分线.定义如下形式的目标函数,最优的剖分线使该目标函数的值达到最小.SL=c1<+c2R+c3l+c4A(1)式中:c1,c2,c3,c4为4个系数;<,R,l,A为4个几何参数,其定义如下.<=1A1+A2<P/2andA3+A4<P/2L(F1,F2)P/2≤A1+A2≤2P/3orP/2≤A3+A4≤2P/3 7(A1,A2,A3,A4)A1+A2>2P/3andA3+A4>2P/3(2)L(F1,F2)=1+ûF1F27(A1,A2,A3,A4)û(3)7(A1,A2,A3,A4)=ûA1-A2û+ûA3-A4û(A1+A2+A3+A4)(4) F1=(A1+A2)-2P/3P/6P/2≤A1+A2≤2P/31otherwise(5)¹收稿日期:20040423 刘　阳:男,25岁,硕士生,主要研究领域为有限元计算及有限元前处理研究 *浙江省教委回国人员基金资助(G59990)F 2=(A 3+A 4)-2P /3P /6P /2≤A 3+A 4≤2P /31otherwise(6)R =R i +R j200　R i =800≤Hi ≤2P /302P /3<H i ≤2PR j =800≤H j ≤2P /302P /3<H j ≤2P (7)l =l splitl char(8)A =ûa 1-a 2û(a 1+a 2)(9) <是角度控制参数,使得剖分线尽可能靠近角平分线,避免内角小于90°,使得生成的单元内角过于尖锐.式(2)中的A 1,A 2,A 3,A 4分别是剖分线与边界线段组成的角度值,如图1所示.文中对式(3)做了修正,原文的公式如下[2].L (F 1,F 2)=(1-F 1F 2)+F 1F 27(A 1,A 2,A 3,A 4) 很显然原文的公式与对称性要求是矛盾的,当F 1<0,F 2=1的时候,剖分线越接近角平分线,L (F 1,F 2)值越大,该剖分线被选中的几率减小.运用式(3),在相同情况下,L (F 1,F 2)值越小,该剖分线被选中的几率增大.R 是结构化控制参数,H i ,H j 分别是点i ,j 的内角值,即H i =A 1+A 2,H j =A 3+A 4.l 为长度控制参数,l split 是剖分线的长度,l char 定义为待剖分区域的特征线长度,l char =(x max -x min )2+(y max -y min )2.A 是对称性参数,A1,A 2分别是两个子区域的面积值.式(1)中4个系数的取值为:c 1=0.52,c 2=0.17,c 3=0.17,c 4=0.14.图1　剖分线选取示意1.2　剖分线上新节点的布置剖分线上节点的布置需要满足以下两个要求:(1)点在剖分线上的间隔值必须线性地分布,以保证单元尺寸渐变过渡;(2)剖分线上新插入的节点个数,应保证两个子区域的边界节点个数均为偶数.假设剖分线端点的单元尺寸为d A ,d B 不失一般性,可以假设d A <d B ,剖分线长度为l .内部节点数目可按下式确定[3].N =int(l (d B -d A )ln(d Bd A))-1d A <d B int(ld A)-1d A =d B(10) 式(10)中的int ()为取整函数,根据准则b 的奇偶性要求,适当增减N 的值,得到剖分线内部节点数目.然后可根据文献[1]中所用的方法布置新节点.1.3　六节点子区域的剖分按上述算法进行剖分时,最后会形成六节点和四节点子区域,四节点子区域即是所求的四边形网格,对于六节点子区域如果仍旧按照该算法进行剖分,则可能出现网格质量较差的网格.文中对每一类六节点子区域分别处理,以获得质量较好的网格.实际编程时,可根据六节点子区域中内角为180°的顶点个数及它们之间的相对位置关系进行区分.1)没有内角为180°的顶点　对于这种情况,可直接在区域内部寻找最短的剖分线,把区域剖分成两个四边形网格,如图2a )所示.a )没有内角为　180°的顶点 b)有一个内角为　180°的顶点图2　六节点子区域处理2)只有一个内角为180°的顶点　与情况1)类似,可直接将此顶点与区域内另一顶点相联,得到两个四边形网格,如图2b)所示.3)有两个内角为180°的顶点　根据两个顶点的位置关系,可进一步细分为三类.(1)两个顶点相邻,如图3a)所示.顶点5,6的内角为180°,可在区域内部插入7,8两点,形成右边所示的4个四边形网格.点7可取点2,3,6所组成的三角形的重心,点8可取点2,3,5所组成的三角形的重心.(2)两个顶点之间隔着一个顶点,如图3b )所示.顶点2,6的内角为180°,可在区域内部插入点7,形成3个四边形网格.点7的可取点2,4,6所组成的三角形的重心.(3)两个顶点之间隔着两个顶点,如图3c )·569·　第4期刘　阳等:基于区域递归分解算法的四边形网格生成技术的研究a )两个顶点相邻b )两个顶点之间隔着一个顶点c )两个顶点之间隔着两个顶点图3　有2个内角为180°顶点的六节点子区域处理所示.顶点2,5的内角为180°,连接顶点2,5形成2个四边形网格. 4)有3个内角为180°的顶点　根据3个顶点的位置关系,可进一步细分为3种情况.(1)3个顶点中,有2个相邻,如图4a )所示.顶点3,5,6的内角为180°,可连接点3,6形成两个四边形网格,由点3,4,5,6组成的畸形网格可在网格修匀阶段得到改善.(2)3个顶点都相邻,如图4b)所示,一般这种情况较少出现.顶点3,4,5内角为180°,连线点1,4,并均匀插入7,8两点,形成两个六节点区域,分别按照情况(1)进行剖分,最后得到4个四边形网格.(3)3个顶点两两都不相邻,如图4c )所示.顶点2,4,6内角为180°,可在区域内部插入点7,形成3个四边形网格,点7可取点2,4,6所组成的三角形的重心.a)3个顶点中,有2个相邻b)3个顶点都相邻c)3个顶点两两都不相邻图4　有3个内角为180°顶点的六节点子区域处理2　网格生成的程序实现基于区域递归分解算法的网格生成方法类似于一个二叉树结构.文中建立单向链表(domain list ),链表中每个元素是待剖分的区域,不断循环遍历链表进行区域剖分,直至链表中每个节点都是四边形区域.可建立以下两个类. 节点类(class point ),存储节点坐标值、节点的单元尺寸、是否为边界节点、节点编号、相邻节点的节点号; 区域类(class domain),存储组成待剖分区域的节点.2.1　初始区域的建立按逆时针方向输入边界形成初始待剖分区域.多连通区域,可在内外边界之间添加一条切割线,将多连通区域转化为单连通区域.类似的,对于n 次多连通区域,进行n -1次切割就可转化为单连通区域.将边界按上文方法进行离散,形成初始的边界节点.把初始边界节点存入区域类中,作为链表的初始节点.2.2　区域剖分反复遍历链表中每个元素所表示的区域,在该区域内,寻找最优的剖分线,生成两个子区域,存入链表内,同时删除该元素.直至每个元素都表示四边形网格.2.3　剖分线的合理性检查在进行剖分时,剖分线必须满足以下两个条件,才是合理的剖分:(1)剖分线不能与边界线段相交;(2)剖分线不能出现在区域以外.为此在选取剖分线时,首先要检查剖分线的合理性.具体实现时,可将剖分线与边界线段一一做求交运算,如果剖分线与边界线段相交,且交点位于线段内部,则该剖分线不合理,不能选取.文献[4]对此有具体的描述.3　为提高算法效率而做的改进区域递归分解算法的关键一步是在待剖分区域中寻找最优的剖分线.严格地讲,必须对所有剖分线计算公式(1)的值,选取其中最优一条.对于包含有N 个边界节点的待剖分区域,则需要进行的计算次数为:(N -1)+(N -2)+…+2+1=0.5×(N 2-N ).如果边界节点较多,计算量无疑是巨大的.为了改善这个状况,对于具有较长直线段边界(相对于该边界上的单元尺寸而言)的待剖分区域,考虑到剖分线应具有的特性(式(1)),可在直线段边界上,只选取位于中部的节点作为剖分线的起点,而不考虑其他节点.经过上述处理,可大·570·武汉理工大学学报(交通科学与工程版)2004年　第28卷大减少待选剖分线的数量,使得计算次数明显减少.当待剖分区域具有大量的直线段边界时,算法效率提高尤为明显.4　算例和结论初始区域被划分成四边形网格后,还需要进行网格修匀才能生成质量较高的网格.由于区域递归分解算法充分考虑了几何图形的对称性和网格尺寸的渐变分布,每次剖分都达到局部最优,因此文中运用Laplacian 修匀方法即可获得高质量的网格.以下是几个四边形网格的实例,如图5所示. a)圆周40等分 b)圆环40等分 c)箱梁截面局部 d)网格疏密过渡图5　网格生成算例从以上算例可见,区域递归分解法是一种有效的网格生成方法,在较规则区域可生成结构化网格,在复杂二维区域自动生成渐变非结构四边形网格,生成的网格质量好.文中改进了区域递归分解算法:修正文献[2]中所含的错误;对六节点子区域分别进行处理,提出区分各类六节点区域的方法;最后在提高算法效率方面做了努力.所开发的网格自动划分程序,具有图形交互输入、网格自动划分、材料定义和边界条件添加等功能,在固体力学和流体力学计算中,有较强的实用价值.参考文献1　Talbert J A,Pa rkinson A R.Development of an au-tomat ic t wo-dimensional finite elem ent mesh gener a-tor using quadr ilatera l elements and Bezier cur ve boundar y definition.International Jour nal for Numer -ical M ethods in Engineering,1990,29:1551～15672　Sar rate J,Huerta A.Eff icient unstr uct ur ed quardri-lat eral mesh generation .Inter national Jour nal for Numerical Methods in Engineer ing,2000,49:1327～13503　田春松,胡健伟.平面区域渐变三角网格自动生成.数值计算与计算机应用,1993(4):303～3114　周培德.计算几何——算法分析与设计.北京:清华大学出版社,1999.15～26Quadrilateral Mesh Generation Based on Recursive Decomposition AlgorithmLiu Yang Zhang Wu Tang Jinchun(College of Civil Engineer ing and Architectur e ,Zhe jiang Univer sity ,H a ngzhou ,310027)AbstractThis paper pr esents an efficient algorithm for creating automatic quadrilateral mesh based on a re-cursive decomposition domain of quadr ilateral elements .This technique can generate quadrilateral meshes over complex simple or multiple connected domain.T he paper describes the programe imple-mentation of the algorithm by Sarrate,and redevelopes it by classifying and treating six node sub-do-mains one by one .T here are several example meshes at the end of this paper .Key wor ds :quadrilateral ;mesh ;domain decomposition·571·　第4期刘　阳等:基于区域递归分解算法的四边形网格生成技术的研究。

一种改进的全四边形网格生成方法及优化

一种改进的全四边形网格生成方法及优化李晓辉;李昌华【摘要】在任意二维区域的全四边形网格生成方法中,铺砌法是目前较好的一种方法.但是当边界出现不规则区域时,网格生成完毕后内部会产生一些质量较差的单元,且判断网格交叉现象是否发生和解决网格交叉问题比较困难.为了提高该方法的适应性和可靠性,提出一种新的改进算法,它对边界不规则区域网格生成和解决网格交叉问题等关键技术进行改进,并加入四边形网格优化方法,将改进后的铺砌法应用于船舶有限元网格划分中,取得了较好的应用效果,最后给出算例进行了验证.【期刊名称】《现代电子技术》【年(卷),期】2010(033)008【总页数】5页(P119-122,126)【关键词】网格生成;四边形单元;网格优化;铺砌法【作者】李晓辉;李昌华【作者单位】西安建筑科技大学,信息与控制工程学院,陕西,西安,710055;西安建筑科技大学,信息与控制工程学院,陕西,西安,710055【正文语种】中文【中图分类】TP391.70 引言随着计算机科学技术的快速发展，以有限元技术为代表的数值方法得到了广泛应用。

然而，制约其进一步发展的主要因素之一是网格划分技术和优化技术,特别是在船舶板架的有限元网格划分中，由于精度和质量上的考虑，要求网格单元尽量是四边形。

目前，网格生成方法[1]主要有Delaunay三角化方法、铺砌法(Paving)、四(八)叉树法、行波法等[2]。

其中,铺砌法是由T.D.Blacker 和M.B.Stephenson 提出的[3],该算法重复地在区域边界内部放一层或者铺一层单元，从区域边界开始向内部形成四边形单元。

White[4]重新设计了铺砌算法,单元不是生成一排,而是一个接一个的生成，它生成的网格质量和灵活性要优于其他算法,尤其是边界单元(接近正方形)[5]。

方兴、张武等[6]在White的基础上对铺砌法做了改进，主要是提出一种节点计算方法，但是在复杂区域的相交和缝合处理等方面仍有不足。

金属成形有限元模拟四边形网格自动划分方法

金属成形有限元模拟四边形网格自动划分方法
张清萍;赵国群;王鹏程
【期刊名称】《锻压技术》
【年(卷),期】2002(27)3
【摘要】对四边形网格自动生成方法进行了研究。

将目前网格划分效果较好的铺路法进行了改进 ,重新设计了网格生成的总体算法 ,单元由边界向区域内部不再逐层生成 ,而是逐段或逐个生成 ,直至网格布满整个区域。

本算法还解决了有关的关键技术 ,其中包括网格交叉处理 ,边界的光滑处理 ,网格的整合和光滑等。

最后给出了网格划分实例 ,验证了算法的可靠性和稳定性。

【总页数】3页(P47-49)
【关键词】金属成形;有限元模拟;四边形网格;自动划分;网络划分
【作者】张清萍;赵国群;王鹏程
【作者单位】山东大学模具工程技术研究中心;内蒙古工业大学材料工程系
【正文语种】中文
【中图分类】TG301;TB115
【相关文献】
1.塑性成形有限元模拟六面体网格划分和再划分技术 [J], 张向;陈军;阮雪榆
2.金属成形有限元模拟网格再划分技术 [J], 赵新中;张向
3.金属塑性成形有限元六面体网格自动划分 [J], 虞松;赵国群;王广春;张存生
4.网格再划分技术在金属成形过程有限元模拟中的应用 [J], 吴伯杰;陈森灿;王少林
5.特别适用于金属成形模拟的四边形网格划分方法 [J], 郭晓霞
因版权原因，仅展示原文概要，查看原文内容请购买。

现代建筑物点云平面特征识别方法

河南科技Henan Science and Technology电气与信息工程总第875期第4期2024年2月收稿日期：2023-09-06作者简介：王新静（1984—），男，博士，副教授，研究方向：土地复垦与生态重建；段晨鑫（1998—），女，硕士生，研究方向：摄影测量与遥感。

现代建筑物点云平面特征识别方法王新静段晨鑫姚怡烨（华北水利水电大学，河南郑州　450046）摘　要：【目的】基于现代建筑物点云数据面片特征，提出一种基于随机抽样一致算法的平面分割识别方法。

【方法】该方法先利用三维格网划分来建立空间格网单元，再根据随机采样点来确定局部格网单元，通过随机机制来拟合平面模型，经过局部打分来确定候选模型集，利用法向约束和共面分割来解决过分割和欠分割的问题。

【结果】采用该方法可获取当前最优模型和一致集，并完成点云分割。

【结论】试验结果表明，该方法能对富有平面特征的建筑物进行有效分割。

关键词：点云；分割；局部采样；一致集中图分类号：TP391 文献标志码：A 文章编号：1003-5168（2024）04-0004-05DOI ：10.19968/ki.hnkj.1003-5168.2024.04.001The Plane Feature Recognition Method for Modern Building PointCloudWANG Xinjing DUAN Chenxin YAO Yiye(North China University of Water Resources and Electric Power, Zhengzhou 450046, China)Abstract: [Purposes ] According to the patch characteristics of modern building point cloud data, a planesegmentation recognition method based on random sampling consistent algorithm is proposed in this pa⁃per. [Methods ] In this method, the spatial grid element is established by three-dimensional grid divi⁃sion, and then the local grid element is determined according to the random sampling points, the plane model is fitted by random mechanism, and the candidate model set is determined by local scoring. Nor⁃mal constraints and coplanar segmentation are used to solve the problems of over-segmentation and in⁃sufficient segmentation. [Findings ] Finally, the current optimal model and consistent set are obtained, and the point cloud segmentation can be completed. [Conclusions ] The experimental results show thatthis method can segment the buildings with plane features effectively.Keywords: point cloud; segmentation; local sampling; consistent set0 引言激光雷达技术是一种新型的高精度、多领域三维数据采集技术，使用三维激光扫描技术能获取点云数据，并基于此来建立各种三维模型，从而被广泛应用于文物保护［1］、地形测绘［2］、城市建筑测量［3］、大型构筑物［4］等领域中。

基于Euler操作的四边形网格细分算法设计与实现

ｔｕｄｉａｅａｕｄｖｓｏｅｐｏｏｅ．Ｉｅｅｏｅａｉｎ，ｅｃｎｔｍｉＯｔａａｎｔｅｓｌ，ａｄｉｃａｈｉｒｔｎｉｏｑａｒｌｔｒｌｂｉｉｉｎａｒｐｓｄｎｔｓｐｒｔｓｓｒｈｏａｈｏｅｉａｏｃＳｔｔｎｏｉｎｆｔｃｅａｉｓｈｉｃｂｐｔｎａｅｔｏｓａｓｑｅｔｌｄｅｅｔｅｃｍｂｎｔｎｏｍｅｏｅｅｏｅａｉｎｎｓｂｉｉｉｇｐｏｅｓＢａｅｎｔｅＥｌｒｐｒｔｎ，ｂｔｆｅｅｕｎｉｆｃｉｏｉａｉｆｓａａｎｖｏｏｆｈｓｐｒｔｏｓｉｕｄｖｄｒｃｓ．ｔｎｓｄｏｕｅｅａｉｓｈｏｏｏｏｈｈｔｐｍａｎｕｌｕｄｖｓｏｒｕｄｉｔｒｌｓｅｅｓｃｅｓｕｌｉｒｌｄｄａｂｉｉｉｎｆａｒｌｅａａｓｏｑａｍｅｈｓａｃｓｆｌｉｌｍｅｔｄｂｅｅｔｇｈｌｅｇａａｓｒｃｕｅｒｕｙｍｐｅｎｅｙｓｌｃｉａｆｄｅｄｔｔｕｔｒ．ｎＫｅｒｓｓｂｉｉｉｎｓｒａｅ；ＥｕｅｌｓＥｕｅｐｒｔｎ；ｑａｒｌｔｒｌｕｄｖｓｏ；ｈｌｅｇａａｓｒｃｕｅｙｗｏｄ：ｕｒｌｒｅａｉｓｕｄｉｅａｂｉｉｎｏｏａｓｉａｆｄｅｄｔｕｔｒｔ
摘要：通过分析细分曲面在Ｅｌｒ则下的性质，出了适用于四边形网格细分的基本Ｅｌ操作方法。在这些方法中，ｕｅ规提ｕｅｒ任

四边形有限元网格生成方法关键技术

有限元是工程界应用最为广泛的一种数值模拟方法。在有限元法中，首先需要把分析域离散为一系列单元。另外，在大变形工艺的模拟过程中，网格也常常需经多次重新划分，因此，通用性强、自动程度高且能产生高质量单元的网格生成法对有限元模拟来说非常重要。它不仅关系到有限元模拟是否能成功地进行到底，而且对模拟的精度也有很大的影响。目前，典型的四边形网格生成法有：映射法、波前法、基于栅格法等。但映射法只能生成结构化的网格，且需要人工事先将被划分区域分割为一系列可映射的子区域，自动化程度相对低；基于栅格法所生成的网格与选择的初始网格及其取向有关，且网格边界单元质量较差；波前法中具有代表性的为 Paving[2]法，该方法涉及的技术较多，程序实现难度较大。本文提出的四边形网格生成法的基本原理已在前期的论文[3]中阐述，现对其中的一些关键技术作进一步改进，使该方法更加完善，程序的运行速度更快。
θ2 θ4 θi θj
(d A = d B ) （3）
Si = m * i + b (d − d B ) 式中， m = A N
N +1 i =1
（4）；b = dB − m 。
理论上讲，节点距离之和
∑S
i
应等于线段
的长度 l ，可实际计算出来的距离之和却不一定等于线段的实际长度。因此，对每个节点间距 Si 都乘以一个调整系数，即节点距离之和与实际线段长度之比。如图１所示，在长度为 10 的一条线段 AB 上插入 C，D，E 3 个节点，A 和 B 2 端点的单元尺寸分别为 2 和 5，由（4）式求得节点间距离 Si（i=1, 2, 3, 4）分别 2，3，4，5，但其距离之和不等于 10。为此，每个节点间距 Si 都乘以一个调整系数 14/10 ，这样节点间的间距变为 1.43，2.14，2.86，3.57。

【doc】平面任意区域四边形网格自动生成的一种方法

平面任意区域四边形网格自动生成的一种方法第10卷第5期1998年9月计算机辅助设计与图形学J.CAD&CGV o1.10.No.5Sep..1998Z～平面任意区域四边形网格自动生成的一种方法/V顾元宪马正阳关振群0矗豇垦工夭学f茬豸大连l1∞24)弋,7Z摘要在改进节点连接法的基础上,提出了一种平面任意区域的有限元同格仝自动剖分方法.既能快速生成四边开j单元同格,也能生成三角形单元网格}而且可很好地控制单元的质量,特别是边界单元的质量;同时可以是活地进行局部同格加密关键词鬯,乏堡中图分类号TP391边形单元网格C AMETHODOFAUTOMATICGENERATIONOFQUADRILATERAL—ELEMENTMESHESONARBITRARYPLANEDOMAIN GUYuan—XianMAZheng——Y angGUANZhen——Qun(DepartmentofEngineeringMechanws?Dah'anunawrsityTechnfogy,Dalian116024~ AbstractAnewmethodforautomaticgenerationoffiniteelementmeshesonarbi—trary2Ddomainisproposedherebymeansofimprovingthenode—connectionmethod. Themethodcangeneratequadrilateral—elementmeshesaswellastriangular—element mesheseffectively+Thequalityofmeshes,particularlyofelementsONboundaties.can bewellcontrolled,andthelocalmeshescanberefinedeasily. Keywordsmeshgeneration.finiteelement.quadrilateralelement1引言有限元方法在应用中的一个主要困难是建立有限元模型,即剖分有限元网格,建立计算模型需耗费很大的工作量.NN~F4快速地剖分高质量的有限元网格是一个重要的研199791?07收稿r1997—04—20收到修改穑.本文得副杰出青年科学基金和国家自然科学基盎资助.嗣元宪.1954年7月生t教授,博士生导师t主要研究方向为结构优他,有限元方法及前后处理,科学计算可视他,计算工程软件马正阳t碗士研究生,主要研究方向为有限元阿格剖分,CAD软件.关振群,讲师,主要研究方向为CAD,计算工程软伴5期顾元宪等:平面任意区域四边形网格自动生成的一种方法433究课题,学者们已提出了许多网格自动生成的算法.在众多的网格生成算法中,三角形网格的自动生成算法最为成熟,如三角化方法" 等.但在有限元分析中三角形单元的精度低,为了满足精度分析的需要,又出现了一些直接生成四边形单元网格的算法.如四分法,铺路法:,以及中轴线法等.这些方法可直接生成质量较高的四边形单元网格,但算法比较复杂.还有一种较为简单易行的方法,是将已有的三角形单元网络通过每两个三角形的合并来得到四边形单元网格.如波前法.这种方法一方面存在着算法的复杂性问题,如在网格的生成过程中需不断地生成新的波前序列;另一方面网格的质量缺乏保证,因为两个质量较好的三角形单元合并而成的四边形单元质量并不一定好.本文提出一种自动生成三角形网格并转化为四边形网格的算法.由于在三角形网格生成时就考虑了合并后四边形网格的质量,围此不但可快速生成高质量的三角形网格,而且最终能生成高质量的四边形网格.整个算法分为布点,三角形网格生成,四边形网格生成三个部分.首先根据网格密度(亦即单元大sb)及分布的要求在边界与平面域内生成有限元网格的节点,然后将这些节点相互联接形成三角形单元.由于布点算法十分简单,同时节点分布规则,围而在生成三角形单元时可避免过多的相交判断,速度很快.在三角形单元生成后,一方面可通过光顺来得到高质量的三角形网格.另一方面,可将三角形单元两两合并来得到四边形单元.同样由于布点的规则性,大部分的相邻三角形单元合并后可得到一个近似于矩形的四边形单元,这就使生成的四边形网格具有较高的质量.另外本文还提供了一种新的数据结构来实现节点一单元关系的快速查寻,从而大大提高了网格的生成速度.由于上述算法三个部分的相互独立性,其中任一部分的改动都不会影响其它部分.因此,可以通过布点算法来实现加密,也可直接处理用任何其他方法生成的三角形网格来生成四边形网格.2任意平面域内的布点算法文献[6]提出一种简单的布点算法,它采用一组平行线与平面域的内外边界相交,在域内截取的部分平行线上均匀插入节点,布点的间隔与单元的平均大小相同,这样就在平面域内得到均匀分布的节点.这种方法对于与坐标轴垂直或平行的边界无疑效果很好,最终生成的边界单元形状接近于矩形,可以很好地拟合边界.但对于曲线边界,就会造成三角形边界单元大小不一,合并而成的四边形单元质量也不高.而边界单元的质量直接影响有限元分析的精度,因而保证边界单元的质量是至关重要的.本文改进了文献E6]的布点算法,使边界单元的形状接近于矩形,从而改善了边界外网格的质量.具体做法如下:(1)先沿平面域的内外边界进行布点,这些节点构成一组新的内外边界,然后把这组新的内外边界作为文献[6]方法中的内外边界进行布点.沿内外边界进行布点时,要根据边界折线的夹角来确定插入节点的方式和位置.从图1所示生成的四边形网格可以看出,用本文的布点算法形成的四边形网格在边界处的质量大大地提高了.计算机辅助设计与图形学)文献[6]的布点算法(6)本文改进的布点算法图1两种不同布点方式生成四边形网格的比较(2)用非均匀布点的方法实现对边界及其邻近区域的网格加密.根据网格加密的要求,通过边界点的非均匀分布来实现对某个区域的加密.具体做法是在沿边界布点时,布点数量及间距由边界点的间距决定,即边界点越密,布点数越多,附近区域的网格也越密.其网格生成实例可见图l】(c),(d).3三角形单元网格的生成本文算法采用逐点生成单元的方法,让每个点与周围其它点相互联接生成三角形单元.该方法充分利用了节点分布的规则性,生成单元时只要不与其周围局部范围内的单元相交就是有效单元,这样就大大减小了单元生成时的判断次数,提高了计算速度.3.1单元生成算法(1)取一节点为基点A(2)取与此节点最近的一点为第2点B;(3)找出在AB左侧与AB最近的3点7',7',7';;(4)在三点中取一个能与AB点组成质量较好且有效的三角形单元的点,记为C,如图2中的T:点;(5)如果ABC三点已经组成单元,贝9转(6);否则,在单元表中插入由ABC三点组成的一个新单元后转(6);(6)将C点做为新的B点,如果新的B点与初始的B点重合,则退出;否贝9转(3); 至此,一个点周围的单元全部形成完毕如图3,可转至下一个节点重复上述算法.直到遍历所有的节点,这样平面域就被全部剖分为三角形单元.B图2从,z,s中选取最好的单元ABT.图3在基点周围形成单元5期顾元宪等平面任意区域四边形网格自动生成的一种方法3.2三角形单元质量和有效性的判断三角形的质量可以用一个参数a来表示,一,a的取值是一个经验值,本文推荐的a值是≥0.1.几种典型三角形的值如图4.△图4几种典型三角形的值三角形单元的有效性是指生成的单元不与已有的单元相交.由于节点的规贝9性,只要判断ABC组成的单元不与A,B,T,,T相邻的单元相交,即可确定ABC是有效的单元3.3单元查找算法在三角形有效性判断以及后面要提到的光~N(SMOOTH)处理中,都要大量地使用以下几个查找算法:(1)查询与一点相连的所有单元号;(2)查询与一点相连的所有节点号;(3)查询与一个单元相连的所有单元号.其中算法(1)是最基本的,其它两种算法可由它来导出.这几种算法都是比较费时间的,因而也是决定网格生成速度的主要因素.本文提供了一种节点一单元邻接链表的数据结构来实现快速查询算法.以三角形单元为例.首先建立三个数组:Element[1..Nelem][1..3]存储已生成的单元的节点号;Nelem为已生成单元数;Einfo[1.Nelem][1…3]是一个初值为0的空数组,其每个元素位置与Element数组一一对应;Ninfo[1..Nnode]是一个初值为0的空数组,Nnode为总的节点数.Ninfo与Einfo组成一个节点一单元邻接关系链表.当一个新的单元形成后,设其单元号为E,三个节点为N(1,3).将三个节点号依次存入Element数组中,如果此时Ninfo[N~]值为0,则将E填入Ninfo[N,]中;如果此时Ninfo[N,]值不为0,则先将Ninfo[M]值填入Einfo[E][i]中.再将E填入Ninfo[N]中. 这样当所有的单元都形成后就可以得到完整的节点一单元邻接关系链表.欲查询与一个节点N相连的单元号时,先查找Ninfo[Ⅳ]中的值,此值即为与节点Ⅳ相连的第一个单元号,记为E1.然后在查找Ⅳ在Element[E1]E1..3]中的位置,设其为i(1一<<一3),此时查找Einfo[E1]l]中的值.此值即为与节点Ⅳ相连的第二个单元号E2.再重复E1的过程,直到目的值为0时则查找结束.下面给出节点一单元邻接关系链表形成及查询与一节点相连的所有单元号的C语言算法:(1)节点一单元邻接关系链表形成算法(2)查询与一节点相连的单元算法设新形成的单元号为E设节点号为Ⅳfor(/j1"<一3;++)E=NionfoEN]fwhile(El—O)Ⅳ.;Ele㈣t[E][]{if(Ninfo[N.3-0)printf(,E);计算机辅助设计与图形学(for(一1,f<一3;++)Ninfor[N~]--E}{}if(ElementEElJrf]一一)elsebreak;{)Ei~o[E][i]--Ninfo[M];E--Einfo[E]Ei]}Ninfo[N.]一E1))}由上述算法看出在生成单元的过程中可以很方便地建立这种节点一单元邻接链表数据结构,从而在查询节点与单元相互关系时用info与Einfo的链接关系代替了与所有单元节点号的比较,因此可显着地提高查询速度.由图5所示的单元节点关系得到的节点一单元邻接链表数据结构如下:单元号4E]ementl23l3451461621Einf0Ninfo节点号0001l00220334l452345圆图5节点1周围单元图6三角形单元的光顺处理图7四边形单元的光J匝处理3.4光顺处理如果需要生成三角形网格,在所有的三角形单元都生成后,可对全部单元进行光顺处理以提高网格质量.光顺处理就是对每一节点的位置根据其周围节点的位置进行调整,以期得到形状较好的三角形网格.如图6中'是节点^光顺后的新位置.这里我们采用拉普拉斯算法,见式(1),其中Ⅳ.为与点相连的所有节点.由于采用了 2.3节中的查找算法,使光顺处理的速度大大提高了..X一1∑X,,=∑(1)…H—l…Hl如果要生成四边形单元网格,这一步可以省略,因为对于一些形状不好的三角形单元5期顾元宪等:平面任意区域四边形网格自动生成的一种方法在合并后可能生成一个质量较好的四边形单元,尤其对于上述方法生成的三角形网格更是如此.4四边形单元网格的生成本文采用的四边形网格生成算法是基于一个三角形单元网格,将三角形单元逐个与其相邻的单元进行最优合并.这样一来合并的过程就与三角形网格的形成方式无关,因而可以适用于任何方法生成的三角形阿格.4.1单元生成算法由于一个三角形单元最多可以与3个三角形相邻,所在存在三种合并成四边形单元的可能,因此要选择一个最优的合并方式,使合并后的四边形单元形状最接近于矩形.由于前述方法布点的规则性,使三种合并方式中必然只有一个最接近于矩形,因此可以用最简单的内角判断来选择最佳的合并方式.即让四边形的4个内角大小控制在45.一135.之间.此外,还要处理下面的问题:(1)优先合并单元:在三角形阿格边界处存在着一些只有一个相邻单元的单元,我们把这些单元列入优先合并的单元,首先加以处理,以尽量减少后文所述的"孤岛三角形"的数量.注意在每合并一个四边形单元后,都有可能产生新的需优先合并的单元,因此要时刻判断并找出有无新的优先合并的单元,如有,应优先处理,如图8.优先台并单元腰由于^研Ⅱ且D合并后新形成的优先台并单元c凸EECF从与删日邻的三个单元中选择最好的进行台并图8四边形单元的生成(2)孤岛三角形单元:在合并的过程中,一些三角形单元周围的相邻单元都已被合并,或是周围的相邻单元都不能组成质量较好的四边形单元,这样的三角形单元称为"孤岛三角形".在所有能合并的单元都合并完后,在"孤岛三角形"内部插入一个新的节点.使其形成3个四边形单元;而其它四边形单元也要相应地分为4个四边形单元,以保证全部为四边形单元,如图9.(3)生成单元的数量:由于对孤岛三角形单元的处理,使生成单元的数量比预计的要.多一倍.为了避免这种问题,在生成三角形阿格时采用的边界点是间隔的实际边界点.如图9中的A点,当合并完三角形单元,进行单元细分处理时,边界上的单元细分要使用先吁雷图9孤岛三角形的处理计算机辅助设计与图形学前未使用的边界点(如图9中的a点)来作为插入的新节点.这样就保证了边界单元很好的拟合于原来的边界,同时降低了生成网格的密度.可以对整个理与三角形Einf0的维图10一…定警晶采用公式其中1兰1三x一(x+x一+2X一),Y(y+y4-2Y)(2)4.3四边形单元质量的判定文献[7]用参数D表示四边形单元的质量,见式(3),4个顶点n,b,c,d如图lO所示. O≤D≤1.0,当四边形abcd为正方形时D=1.0,当四边形abcd有3个以上的点在一条直线上时D一0.D值越大,单元质量越好,D的平均值可表示所有单元的平均质量.在图11所示阿格剖分实例中,给出了单元质量参数D的平均值和最小值D….其中最小的D值为0.705,表明单元质量是好的.露c)方孔边局部加密,D一0.825,D～=0762(d)凹点附近局部加密.=0846,Dm.=0.789 图11几个阿格生成实倒5期顾元宪等:平面任意区域四边形阿格自动生成的一种方法4395结论和展望本文给出的算法,经实例验证可以自动生成较高质量的四边形或三角形单元网格.其优点是方法简单,运算速度快,边界区域的单元质量好,便于局部丽格加密,适台于任意复杂形状的平面区域.三角形丽格生成的节点连接法可以推广到三维实体区域,用以自动生成四面体单元网格.参考文献lCavendishJC.Automatictriangu]atJonoarbitraryplanardomainsforthefiniteelementmet hod.Int盯【Ld0阻lJournalOrNumer~a1MethodsinEngineeringt1974,8|679—6982SbepardMS,Y erryMA.ApproachingtheautomaticgenerationoffinitedementmeshComp uterMechanicalEngineering,1983—1:49—563BlackerT.StephensonMB.Pavir~iAnewapproachtoautomaticalquadrilateralmeshgener ation.In々er【Ⅲti0【Ⅲl JournalforNumerlcMMethodsinEngineering,1991,32:811--8474TamTKH,ArmstrongOG.2Dfilfiteelementmeshgenerationbymedlaraidsubdivision.Ad vanceEngineeringSoftware,1991,13;313—3245ZhuJZ,ela1.Anewapproachtothedevelopmentofautomaticquadrilateralme*hgeneratio n.Inter~tiona]Jour-halforNumerica]MethodsinEngineering一1991,32:849—8666LoSH,AnewmehgenerationschedulearbitraryptanardomainstIntJNumerlcolMethodsi nEnginerring,l985,21:1403—14267李华.平面曲面区域的四边形阿格自动剖分[博士学位论文].大连;大连理工大学工程力学系,1995。

有限元四边形网格划分的两种算法

有限元四边形网格划分的两种算法
潘子杰;杨文通
【期刊名称】《机械设计与制造》
【年(卷),期】2002(000)002
【摘要】有限元四边形网格划分要比三角形网格划分难很多,这里通过两种方法解决四边形网格划分问题.一种是把已有三角形网格转化为三角形/四边形混合网格,另一种是通过映射法生成全四边形网格.两种方法都具有计算效率高,四边形网格质量好的特点.
【总页数】2页(P50-51)
【作者】潘子杰;杨文通
【作者单位】北京工业大学,北京,100022;北京工业大学,北京,100022
【正文语种】中文
【中图分类】O242.21
【相关文献】
1.基于特征约束的四边形网格划分算法 [J], 刘虎;翟建军;陈文亮
2.四边形网格划分过程中的边界提取与优化算法的实现 [J], 齐维毅;丁言镁;吴丽娟
3.基于特征约束的四边形网格划分算法 [J], 刘虎;翟建军;陈文亮
4.二维图形的四边形有限元网格划分 [J], 刘晓飞;金泉林
5.四边形有限元网格划分方法——二分法的改进 [J], 郭晓霞
因版权原因，仅展示原文概要，查看原文内容请购买。

基于Euler操作的四边形网格细分算法设计与实现

基于Euler操作的四边形网格细分算法设计与实现
郑立垠;周笑天;于萍;张晓伟
【期刊名称】《计算机工程与设计》
【年(卷),期】2007(28)13
【摘要】通过分析细分曲面在Euler规则下的性质,提出了适用于四边形网格细分的基本Euler操作方法.在这些方法中,任意一种操作都是不可分割原子操作,而每一次细分迭代实际上就是这些基本操作顺序有效的组合.在Euler操作的基础上,选择halfedge数据结构,四边形网格的初始细分和对偶细分得以成功实现.
【总页数】4页(P3151-3153,3156)
【作者】郑立垠;周笑天;于萍;张晓伟
【作者单位】中国石油大学(华东)计算机与通信工程学院,山东,东营,257061;中国石油大学(华东)计算机与通信工程学院,山东,东营,257061;中国石油大学(华东)计算机与通信工程学院,山东,东营,257061;中国石油大学(华东)教务处,山东,东
营,257061
【正文语种】中文
【中图分类】TP391
【相关文献】
1.基于四边形网格细分曲面的生成 [J], 杨俊清
2.基于四边形网格参数细分的平面与自由曲面求交算法 [J], 李慧莹;陈良骥
3.基于C-C细分的四边形网格上插值光滑Bézier曲面的生成 [J], 蒋跃华;陈志杨;
李根;张三元;叶修梓
4.基于四边形网格的新细分算法的研究 [J], 王其华;孙立镌
5.基于不分裂算法的Euler网格细分方法 [J], 宁建国;刘春;马天宝
因版权原因，仅展示原文概要，查看原文内容请购买。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

定义３前沿边．．位于兰角形区域和四边形区域分冲线上的边．定义４前沿结点．．位于前沿边上的顶点．定义５次前沿边．于三角形区域中的非边．位界边，其两个端点中至少有一个端点是前沿结点．定义６边界环．．由一系列边界边首尾依次相
法往往比较复杂或受到一定的限制．为了减少四边
等〔」 ‘提出了前沿算法的线性改进算法— 分级算。
法，该算法虽然简单，但是生成的四边形网格存在大
量的残余三角形；赵熠等［ＪＩ提出的泛前沿算法虽然１
克服了前沿算法的缺陷，但是在最初生成的四边形网格中仍存在一些残余三角形，需要通过对这些三角形进行移动来消除，增加了算法的复杂性针对上述算法的不足，本文提出了一种将三角形网格转
” 二ｔ
法达到理论最大值，只能获得一个最大近似值ｗ．
在实际算法中，。Ｅ用与其相对应的内。边Ｌ来代替．
为了获取命，本文提出了区域生长算法来生成四边形网基本思想如下：格，通过三角形两两合并生成一个连通的四边形区域和一个连通的三角形区域，每一次三角形合并操作都从生成的次前沿边中选择具
成的四边形网格存在少量的残余三角形；阂卫东
四边形网格［；。〕进行参数曲面（１在Ｕ中，３如ＮＲＳＢｚｒ－ｌｅ重构时也需要事先生成ＵＢ，，ｉ等）ｅｅＢＳｎｉｐ相应的四边形网格〔１与全自４．－５动生成三角形网格相比，前只有少数的方法能直接生成四边形网格，目如四叉树法［铺路法１、［６１７映射法［等，１ “ 而且这些算」
门户－一
内一ａｌ
＜ａ４
ｘ２ａ
图１四边形的４个外角
凸四Ｒ边形的值介于０１和之间，当四边形退化为三角形时ＰｏＲ＝；值越大，表明四边形质量越好，越接近于矩形．由于凹四边形的Ｒ值也大于。，因此无法通过Ｒ值来表明四边形的凹凸性．如图２
连接Ｆ，３Ｄ
Ｂ Leabharlann 图２对角线与４条边形成的４个夹角
刀＝１
刀－＿ｐ
刀＝．０２５
卢＝．０３；ｉ
刀＝１．０７４
图３一些典型平面四边形的Ｒ值
３四边形网格的生成
３１相关定义．
定义１度．网格中与某个顶点相邻的元素（．三
１引
言
间接算法的一个重要环节．研究人员已经提出了很多由三角形网格生成四
在有限元分析中广泛使用三角形网格，但有时为了提高分析精度或满足一定的特殊要求需要建立
边形网格的算法，ｈ等［提出了前沿算法，Ｚｕ［９，该算法需要进行烦琐的前沿拆分、合成和求交操作，且生
ａｄａｔｑａｒａｒｌｓｅａｄｍｃｎｔｉｎｃｎｉｏｓｉｐｓｏｔｅｒｉｇｏｑｌｙｄｉｔａｍｅｈｓｓｅｓａｉｏｄｉａｅｏｅｎｍｇｇｏｕｉｕｌｅｎｏｏｒｎｇｔｎｒｍｄｈｅｎｐｏｅｕｅａｏｇｎｒｔｇｌｅｔａｇｓＴｅｓｍｓｏｔｉｏｅｔｎａｔｅｔｒｉｔｅｒｃｄｒｔｖｉｅｅｉｉａｄｎｌ．ｎｅｏｈｇｒｉｓａｎｅｎｈｏｄａｎｓｔｒｅｈｏｍｏｉｎｐａｏｒｋｏｅｅｆｉｔｌｄｉｔｒｌｓｅ．ｎｉｑａｒａｅａｍｅｈｓｉａｕｌＫｙｒｓｅｈｅａｏ；ｄｉｔａｍｅ；ｎｕｒｓｅｗｄｕｒａｒｌｈｔａｇｌｍｅｈｏｍｓｇｎｒｉｑａｌｅｅｔｎｓｒｉａ
研究方向为逆向工程．ＡＣＤ＆Ｃ陈康宁，１３Ｇ．男．９５生．教授博士生导师主要研究方向为计算机集成质量系统
万方数据
计算机辅助设计与图形学学报
２００４年
Ｚ四边形质量因子ｐ
定义四边形质量因子１３的方法比较多［１，１４本２１文根据四边形的４条边所夹的４个外角来定义Ｒ值．
第１６卷第１期
２００４年１月
计算机辅助设计与图形学学报
ＪＯＵＲＬＦＮＡＯＣＯＭＰＵＴＥＡＩＥＤＳＧＮ＆ＣＲＤＤＥＩＯＭＰＥＧＲＨＩＳＵＴＲＡＰＣ
Ｖｏ．Ｎ．ｌ１，１６ｏＪ，２０－０４
结束．
３４三角形合并的约束条件．
在三角形合并过程可能发生的各种情形如图６
所示．
图５由三角形网格生成的带权无向图
要使合并后生成的四边形网格质量最优，理论上相当于从无向图中选择Ｋ／２条边进行三角形合
并作，得ｗ一艺ｃｉ大但际用无操使ａｎ．实应中；最，
角形或四边形）的个数．定义２边界边．．网格中只有一个元素（三角形
或四边形）与其相邻的边
行三角形合并操作之前需要对三角形网格做预处理．图４所示为网格中度为３４的顶点，和由于与这类顶点相邻的三角形在合并后很容易生成凹四边形，因此可以先将这类顶点从网格中删除，分别做如图４所示的操作．
ＡｓａｒｈｔｉｐｖｔｑａｔｏｇｎｒｉｑａｒｔａｍｓｓｂｔｃｒｔＡｉｇｗｈｉｍｍｒｅｕｌｆｅｔｇｄｌｅｌｈｆｍｒｏｒｔａｏｔｏｏｈｅｎｇｏｌｇｅｉｙｅａｎｕｉｒｅｅｒａｏｔａｇｌｏｅｉｐｐｓ．ｍｔｄｉｉｓｗｙｅｉｔｏａｅｔｎｌｔｇｎｒｅｒｎｕｒｓｒｏｅＴｅｈｏｔｚｔｅｏｍｒｎｗａｊｎｔａｇｓｅｅｔｉａｎｓｏｄｈｅｏｐｍｅｈａｆｇｇｄｃｒｅｏａｉ
为删除边，合并与其相邻的两个三角形生成一个新四边形．
Ｓｐ．Ｆｎｄｅｉ，Ｆｏｔｄｅｉ和Ｆｏｔｔ４更新ｒｔｇｓＳｂｒＥｇＬｓｒ－ｅｏＥＬｔｕｎｔｎ从ｕＬｓ如果Ｓｂｒｔｉ非空，ｌｉ，，ｔｕＦｏＬｓｎｔ则返回Ｓｐ；ｔ３否则算法ｅ
一争
几
岁
凡
弓
所示，对于四边形ＡＣ连接顶点Ａ，ＢＤ，Ｃ形成对角线，得到四边形的４条边与此对角线所形成的４个
夹角分别为Ｙｌ７和八；Ｙ十Ｙ李１ ’ ｌＹ，２３若，２８或０Ｙ十Ｙ１１＇３４８，＞０则表明该四边形是凹四边形；同理，
由区域生长算法实现四边形网格划分
胡向红
摘法
陈康宁
〔西安交通大学机械工程学院机械电子工程系西安７０４）１０９
要针对当前各种三角形网格转化为四边形网格算法的缺陷，出一种实现这种转化的算法— 区域生长算提该算法通过有选择地合并三角形来获得高质量的四边形网格，并且对两个三角形的合并增加约束条件以避免在
ＨＸａｇｏｇｈｎｎｎｇｕｎｈｎｉＣｅＫｇｉａｎ
（ｅａｔｅｔｈｒｎ，ｎｉｅｎ，－ＪＭｃｎａＥｇｎｒｇＸ＇ｎｏｍＵｉｒｔ，＇ｎ１０９ＤｐｒｎｏＭｎａｍｉＥｇｎｒｇＳ！１ｅａｉｌｉｅｎ，ｎＪｒｇｖｓＸｒ７０４）ｍｆ，ｅｉｃ月ｈｅｎｅｉｉｉ，ｎｅｉａｙｉ．
化为四边形网格的新算法— 区域生长算法．
形网格生成的难度，通常采用间接算法，即先生成三角形网格，然后再将其转化为四边形网格．因此，如
何将已生成的三角形网格转化为四边形网格是实施
原稿收到日期：０－１９修改稿收到日期：０－－本课题得到国家自然科学基金（９７０５资助胡向红男１７生、２２１３０２３４２０００５６５７）９８硕士，主要
生成四边形网格的过程中产生残余三角形；同时对生成的四边形网格进行一系列的质量改善操作
关扭词网格划分；四边形网格；三角形网格中图法分类号Ｔ３１Ｐ９
ＧｎｒｔｇａｒａｅａＭｅｈｓＵｉｇｇｎｏｔＡｇｒｈｅｅａｉＱｄｉｔｌｓｅｂｓＲｉＧｗｈｏｉｍｎｕｌｒｙｎｅｏｒｌｔ
如图ｔ所示，四边形ＡＣ〕４个外角分别为设ＢＬ的
可以作同样的判断，对于凹四边形的ａ可以赋值为一．１一些典型平面四边形的月值如图３所示，
ａ，ａ和ａ，ｌａ，４２３对其进行降序排列，ａｉａ％即，２ａｉａ，３４则四边形的质量因子Ｒ可定义为
并后生成的新前沿与某一边界环相交，并且该边界环不属于前沿的一部分，则将该边界环加人到前沿
非前沿结点
卿四形城边区口三形域角区
图６三角形合并过程中可能出现的１种基本情况７
Ｄ（－１２一，来表示，ｉ，，Ｋ）任意两个相邻的三角形
结点Ｄ和几用一条。且，边Ｅ连接，将这两个相邻三角形合并后生成的四边形质量因ｉ子Ｒ作为权赋给边Ｅ，；则原三，角形网格将形成一个带权无向图，如
图５所示．
四形量因边质子凡，将作权赋这两三并其为值给个角形的公共；从！｝选择边Ｌ，Ｌ中具有最大值的为边生成，权边作删除
连生成的一个封闭环．３２三角形网格的预处理．
ａ度为３的顶点的处理