三维域四面体网格生成的推进网阵方法和网格光顺

本文讨论了应用推进网阵方法在任意形状三维块体内的四面体元网格的生成过程，建议了这个过程中一个典型四面体元生成的一般方案，提出了一个简便的在三维非平面表面域上三角形网格的广义Ｌａｐｌａｃｉａｎ光顺方案和一个块体域内四面体元网格的光顺方案，给出的一些图形结果表明了所提出的网格生成器及其程序的功能和质量。本文工作已作为前处理机在英国有限元分析软件公司ＦＥＡＬｔｄ的通用有限元程序ＬＵＳＡＳ中实现。     

薄元分解与Laplacian光顺相结合的四面体有限元网格优化方法

提出一种有效的三维实体四面体有限元网格质量优化方法以满足有限元分析对网格质量的要求。对薄元分解法进行改进,改进的薄元分解法更全面地考虑了各种劣质单元类型,能够对三维实体网格剖分中产生的各种类别的孤立劣质单元进行有效的分解;将改进的薄元分解法与Laplacian光顺优化方法相结合以解决某些网格剖分算法如推进波前法和Delaunay三角化方法产生的非孤立劣质单元问题。经过实例检验,本文提出的四面体单元网格优化算法健壮有效、易于实现,能够显著提高最差单元的质量。     

用逐点插入法生成Delaunay四面体自适应网格

介绍一种基于Delaunay算法的四面体自适应网格的自动划分方法。该方法用单元尺度场控制生成网格的疏密分布,在不满足尺度场要求的单元面形心处插入新节点,同时计算新节点单元尺寸参数,实现三维实体的Delaunay四面体自动划分。此方法具有几个特点:一是表面网格与体内网格同步划分,无需区分两者;二是结点与单元同时生成;三是生成网格自适应性好,疏密分布任意。另外,还介绍了三维网格划分中两个相关算法:一个是约束面恢复算法,该算法基于约束面不允许有单元边与之相交的性质而提出的;另一个是将二维射线法推广至三维空间,判断一个点是否在一多面体内,实现了凹多面体的划分。最后通过算例对单元质量进行了评价。本文所述方法是一种有效的四面体自适应单元生成算法。     

一种四面体网格局部自适应方法在有限元中的应用

通过提出一种针对三维四面体网格的区域整形技术，给出了一种简便实用的四面体网格局部自适应加密方法，并将其与流量修正有限元法结合，应用于三维高速流的计算     

任意多面体边界一致四面体网格生成方法

提出一种能够满足多面体边界几何与拓扑约束的边界一致恢复算法,解决了任意多面体的边界一致四面体网格生成问题。在恢复多面体的几何约束时,边界上可能会引入Steiner点,这样就不满足拓扑约束。对此,本文采用动态规划方法将Steiner点从边界上消除,修复与其相关四面体单元的拓扑关系,以保持原多面体边界的拓扑完整性,并采用扩展的Laplacian光顺算法优化劣质单元。在理论上,本文算法能够保证完整地恢复任意多面体的边界。算例表明,本文提出的边界一致恢复算法鲁棒性高,可应用于复杂多面体模型。     

面向复杂几何模型的并行四面体网格生成方法

面向大规模工程计算等数值模拟领域,提出了一种支持复杂几何模型的大规模四面体网格并行生成方法。该方法以复杂几何模型作为输入,首先采用串行网格生成方法生成初始四面体网格,然后通过两级区域分解方法将初始网格分解为多个子网格并分配到相应的进程中,进程间并行地提取出子网格的表面网格,并基于几何模型对面网格进行贴体加密,最后对加密后的面网格采用Delaunay方法重新生成四面体网格,该方法可以更好地适应高性能计算机体系结构,较好地克服了并行方法中并行性能和网格质量不能兼顾的问题。对三峡大坝模型进行测试和验证,证明该方法具有良好的并行效率和可扩展性,可以在数万处理器核上并行生成数十亿高质量四面体网格。     

网格生长动力学模型的建立和求解及其混合黏性网格生成中的应用

真实飞行器外形复杂,涉及的流动结构和现象对数值模拟中所需的网格提出了较高的要求。本文首先研究了各向异性黏性网格的生成方法,提出了网格生长动力学模型,并建立了松弛算法求解该模型,从而得到合理的层状网格;提出并建立了各向同性四面体网格的阵面推进Delaunay网格细化算法。最后通过数个算例对所建立的混合黏性网格生成方法的适用性和可靠性进行考察。结果表明,本文建立的混合黏性网格生成方法适用于复杂飞行器的网格生成,具有较高的可靠性。     

一种三维表面网格自动生成方法

提出了一种改进的区域分解法,用于三维表面三角形和四边形网格的自动生成。基于区域分解法的三维表面网格生成的主要技术包括三维表面三角形面片拓扑关系的建立,表面三角形面片的分组及特征面的确定,空间曲线上网格节点的生成,区域内外边界的合并,最佳剖分面的确定,以及网格的光滑处理等。本文对这些技术进行了详细的介绍。本文提出的表面网格生成方法具有算法可靠、生成网格效率高以及生成网格质量高等优点。最后给出三维表面网格生成实例验证了本文提出方法的可靠性。     

有限元网格全自动生成的四分法

本文论述了有限元网格全自动生成的四分法。这种方法适用于任意形状的平面结构,且形体内部可以带有若干个任何形状的孔洞;自动生成的单元以四点元为主,边界辅以三角形元;该法还允许有限元网格局部加密,加密区交接处引入五点元,网格加密部位可以有若干处。有限元网格全自动划分算法在CAD工作站上实现,用人机交互会话方式控制输入/输出操作,方便、灵活、高效,可做为前处理器与任何一个有限元分析程序系统相联接。本文所采用的方法和数据结构,可以推广到带孔洞的在意形状空间形体,实现对更广泛的三维结构的有限元网格全自动生成,     

适应复杂外形的三维粘性混合网格生成算法

提出一类适应复杂外形的粘性混合网格生成算法。表面网格由前沿推进三角形曲面网格程序获得,边界层布置各向异性的三棱柱体网格,远物面区域采用Delaunay方法生成四面体网格。针对模型的复杂几何特征,综合采用了各种网格处理技术,以保证边界层网格的质量,并避免算法失效问题。网格实例及计算结果表明了本文算法的实用性及和效性。     

参数投影法表面网格生成

本文采用投影方法化三维网格生成为两维网格生成，大大精简了计算量，且通用性强，生成的网格平面性好，正交性强，编制的软件为用户提供了良好的界面，使用方便。     

黏性边界层网格自动生成

高雷诺数黏性流动在壁面附近存在边界层,在计算模拟中自动生成可靠且有效的计算网格仍然是计算流体力学存在的瓶颈.三棱柱/四面体混合网格技术在一定程度上缓解了这个困难.然而,对于复杂外形的情况,在边界层内自动高效生成高质量的三棱柱单元仍然十分困难.常用的层推进法在凹凸区域及角点处生成的边界层网格单元质量较差,边界层网格最外层尺寸不均匀.针对这些问题,发展了一种黏性边界层网格自动生成方法,通过顶点周围边的二面角识别物面网格特征确定多生长方向,预估并调整生长高度处理相交情况.同时提出一种三维前沿尺寸调节方式,提高了边界层网格单元的正交性,保证了边界层网格与远场网格尺寸的光滑过渡.通过ONERA M6翼型以及带发动机短舱的DLR-F6翼身组合体等外形的网格生成实例及绕流数值模拟,将计算值与标准实验值进行对比,结果表明:该方法能够自动高效地生成满足数值计算需求的混合网格.     

基于叉树法的自适应有限元局部网格加密研究

本文以叉树法网格生成技术为基础，提出了效率高，易实现的局部网格加密步骤，同时设计了相互匹配的网格数据结构和树数据结构以及局部网格修正算法用以支持本文所提出的局部网格加密步骤，最后，给出了局部网格加密例子，用以证明本文所提出的方法可以很容易地生成任意密度分布的网格，对于实施自适应ｈ－方案是十分有利的。     

非结构网格变形方法研究进展

动网格技术中的非结构网格变形是计算流体力学和计算固体力学的关键技术之一.本文在总结现有非结构网格变形方法的基础上, 提出了一种网格变形方法的详细分类, 将现有方法归类为虚拟结构法、偏微分方程法和代数法. 本文综述了各类方法的最新研究进展, 分析并比较了各类方法的特性, 评述了当前网格变形研究的几个主要方向:复杂结构外形在不规则变形下的动网格生成、三维动网格生成、并行动网格生成和动节点技术. 最后简要地探讨了该领域的发展趋势.      

基于射线穿透法的GPU并行阶梯型有限差分网格生成算法

三维大规模有限差分网格生成技术是三维有限差分计算的基础,网格生成效率是三维有限差分网格生成的研究热点。传统的阶梯型有限差分网格生成方法主要有射线穿透法和切片法。本文在传统串行射线穿透法的基础上,提出了基于GPU （graphic processing unit）并行计算技术的并行阶梯型有限差分网格生成算法。并行算法应用基于分批次的数据传输策略,使得算法能够处理的数据规模不依赖于GPU内存大小,平衡了数据传输效率和网格生成规模之间的关系。为了减少数据传输量,本文提出的并行算法可以在GPU线程内部相互独立的生成射线起点坐标,进一步提高了并行算法的执行效率和并行化程度。通过数值试验的对比可以看出,并行算法的执行效率远远高于传统射线穿透法。最后,通过有限差分计算实例可以证实并行算法能够满足复杂模型大规模数值模拟的需求。     

基于点球弹簧修匀法的混合网格变形并行算法

针对二维/三维混合网格,提出基于点球弹簧修匀法的并行网格变形算法。按特定模板将混合网格中的非三角形/四面体单元分解成三角形/四面体单元。针对每个内部节点及其相邻节点建立相应的子弹簧系统,并通过增加Ball-Vertex弹簧避免弹簧系统的塌陷问题。由于点球弹簧法在计算中逐点对网格内部节点进行计算,在计算过程中具有良好的弱耦合性质,因此有利于算法并行化。在并行化时仅需对网格进行虚拟分区操作,不必进行复杂的几何分区操作,同时避免了混合网格不同单元之间的兼容性问题。该方法适用于具有复杂外形的大规模混合网格的变形问题,能够显著提高网格变形的效率,同时具有良好的适应性。     

基于转换模板的三维实体全六面体网格生成方法

以四面体一六面体基本转换模板为基础,提出一系列具有伸缩性的扩展转换模板．可根据需要选择不同模板及其组合,将四面体分解为不同数量、不同密度过渡形式的六面体单元。这样,初始四面体网格不需要划分得很细,生成的六面体单元数量也可以通过采用不同规格的扩展转换模板而得到控制。提出了基于CAD几何造型的边界节点坐标修正方法．从而使边界网格更好地拟合几何模型边界。     

大量椭球颗粒随机分布三维区域的模拟及其四面体网格快速生成算法

针对颗粒为椭球形状并且随机分布的多相复合材料给出了三维情形下的计算机模拟方法,并对此复杂区域提供了一种四面体网格快速生成算法。此方法的基本思想:首先是根据材料的要求,给出随机颗粒材料的计算机模拟,然后对模拟后的颗粒边界初始化,得到一些初始剖分点,并给出剖分区域的初始网格(正四面体),接着把颗粒的初始边界点投入到剖分区域中,替换区域中靠此点近的点,并且采用加点寻找边界的方法和修正初始点方法来解决四面体穿刺边界问题,从而得到整个区域的四面体网格。最后用实例说明了算法的有效性。     

用混合网格求解三维可压雷诺平均Navier—Stokes方程

用混合网格求解了三维紊流 N-S方程。在物面附近采用三棱柱网格 ,其它区域则采用四面体网格。方程的求解采用 Jamson的有限体积法 ,紊流模型采用两层 Baldwin-Lomax代数紊流模型。数值算例表明 ,用混合网格求解三维紊流 Navier-Stokes是非常有效的。     

三维随机凹凸型混凝土骨料细观建模方法研究

为了提高混凝土细观力学中数值模拟预处理的要求,提出了一种网格预先生成的三维随机凹凸型混凝土骨料建模方法。该方法利用已生成的网格状态进行骨料投放,有效解决了在细观三相介质中网格质量问题。采用随机凹凸型骨料使骨料模型更加接近真实形态,避免了该方法对骨料形貌描述不准确。同时,提出一种四面体分割法来判断空间三维骨料与单元节点的位置关系,进行细观三相介质单元判断,该方法可保证投放的随机性以及级配分布的合理性,三维投放骨料含量可达到65%。采用本文提出的预处理方法,结合混凝土塑性损伤模型,分析了钢筋锈蚀引起的保护层开裂破坏过程,结果表明:受随机骨料的影响,顺筋裂纹的产生具有不确定性。