生成DDA块体相邻面上刚性点对的Voronoi图方法-北京大学学报医学.PDF

　 第２８卷第１１期 长　江　科　学　院　院　报 Ｖｏｌ．２８　Ｎｏ．１１　　 　　２０１１年１１月 ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＲｉｖｅｒＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ Ｎｏｖ．２０１１　 　 　　文章编号：１００１－５４８５（２０１１）１１－００９６－０４ 生成 ＤＤＡ块体相邻面上刚性点对的Ｖｏｒｏｎｏｉ图方法 １ ２ ２ 张　杨 ，邬爱清 ，林绍忠 （１．扬州大学 水利科学与工程学院，江苏 扬州　２２５００９；２．长江科学院，武汉　４３００１０） 摘要：在高阶ＤＤＡ用于连续结构的计算分析中，需要在块体之间的相邻面上设置刚性的“连接点对”将相互独立的块 体连接成一个连续的整体。对于多构件的复杂结构，快速均匀布点是一个难点。基于Ｖｏｒｏｎｏｉ图理论，改进了几何法 生成Ｖｏｒｏｎｏｉ图的算法，并基于质心法对点的位置进行优化。据此编制了相关程序，并采用算例进行了验证。计算结 果表明：该方法可以快速在任意凸多边形内获得分布均匀的一系列点，能方便地应用于三维高阶ＤＤＡ计算分析。 关　键　词：高阶ＤＤＡ；点对生成；Ｖｏｒｏｎｏｉ图；几何法；质心法 中图分类号：Ｏ２４２　　　文献标识码：Ａ 家Ｍ．Ｇ．Ｖｏｒｏｎｏｉ明确提出Ｖｏｒｏｎｏｉ图的概念，同时将 １　ＤＤＡ计算网格 其扩展到高维空间。１９３４年，Ｂ．Ｄｅｌａｕｎａｙ从Ｖｏｒｏｎｏｉ ［２］ 图中演化出其对偶图Ｄｅｌａｕｎａｙ三角形网 。Ｖｏｒｏｎｏｉ 与有限元方法相比，ＤＤＡ以体为研究对象，具 图和它的对偶图Ｄｅｌａｕｎａｙ三角形的意义在不同的科 有以下特点： 系统由一系列体构成，不存在节点的 ① 学领域被发现和重视，已有近１００年的历史。 概念； 体的形状任意，各体具有相对独立的位移场 ② 随着计算机技术的发展，２０世纪８０年代以来 ［１］ 函数 。前一个特点决定了ＤＤＡ的计算网格不需 Ｖｏｒｏｎｏｉ图及其对偶图Ｄｅｌａｕｎｅｙ三角剖分被广泛研 要像有限元那样满足网格相容性条件（相容网格见 究和应用于越来越多的领域。它们在几何形体重 图１，非相容网格见图２），若将 ＤＤＡ方法应用于复 构、计算机图形学、图形处理与模式识别、物理、化学 杂连续结构的计算，前处理将变得非常简便、快捷。 和分子生物学以及机器人、ＧＩＳ、城市规划、油气藏数 后一个特点是每个 ＤＤＡ块体有自己独立的位移函 值模拟、编辑排版、监控网络、经济决策过程、多边形 数，这一优势决定了可以方便地通过提高模拟函数 网格单元剖分和数值计算、有限单元网格剖分、无单 阶次来提高解题精度，无需增加节点，仅仅需要在体 元法布点等方面都有应用。 与体之间满足一定的变形协调条件。因此，位移函 下面介绍一下Ｖｏｒｏｎｏｉ图的概念和特性： 数升阶更加方