四面体网格质量度量准则的研究

简要回顾了度量四面体单元质量的几个常用准则。长期以来，这些准则被认为在一定意义下是相互等价的。数值试验表明，对单元形状的变化用不同的准则度量有可能产生矛盾的评判结果，从而影响到网格质量优化的进程和结果。因此在进行网格剖分或网格优化时，应对所采用的度量准则做仔细的研究和测试。     

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

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

三维六面体有限元网格自动划分中的一种单元转换优化算法

本文采用十节点曲边四面体转换为六面体网格,并采用非线性约束优化算法取Ｌａｐｌａｃｉａｎ光滑处理算法有效地提高六面体单元的质量,实现了对任意实体的六面体网格自动划分。     

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

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

两级分块结构化网格在SK型静态混合器中的应用

给出一种适用于SK型静态混合器流场数值模拟的两级分块结构化网格划分方案和分块圆弧区点的极坐标方程。与常用的非结构四面体网格划分对比,两种方法划分网格数量相近时,两级分块结构化网格达到的迭代精度高于非结构化网格一个等级,而非结构化网格收敛所用时间约为结构化网格收敛所用时间的1.25倍;运用控制体积容积变化值对两种方法划分网格进行质量评估,结构化网格划分方法接近99%网格质量指标在1.0~2.0范围内,而非结构化网格仅接近60%网格质量指标在1.0~2.0范围内。通过实验测量对两种方法划分网格数值模拟结果进行对比验证,结构化网格计算结果与实验值更为接近,表明其可行性和准确性均优于非结构四面体网格划分。     

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

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

基于点球弹簧修匀法的混合网格变形并行算法A parallel algorithm for hybrid mesh deformation based on Vertex-Ball Spring Smoothing

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

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

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

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

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

基于分区拼接网格技术高升力装置流场数值模拟

针对高升力装置构型模型结构复杂、流场变化剧烈等特点,本文采用分区拼接网格的思想分别按照流场和结构拓扑特点对高升力装置进行了网格分区。在分区的基础上逐块生成网格,减小了增升装置网格的生成难度,提高了网格质量,减少了网格数目。首先,研究了高升力装置的分区策略及流场特点;接着以MD30P-30N多段翼型为研究算例研究了网格比例和插值方法对计算结果的影响。经过分析对比可知:外部区域与近壁区域之间的比例不宜小于1:5;内部域网格比例不宜超过1:1.8,最好保持在1:1左右;计算中应该采用高阶精度插值以保证计算精度。采用某四段翼型进行了验证;最后采用NASA标准高升力装置进行了三维高升力装置流场数值模拟并与相应风洞实验数据及对接网格计算结果进行了比较与分析,验证了拼接网格技术的高效性与可靠性。同时分析研究了绕三维增升装置的流动及其周围复杂的粘性流动现象。     

判别二维有限元网格图的简便算法

根据图论中欧拉公式推得二维有限元网格图中单元数、节点数、边数及边界数之间的关系式,获得判别网格图中是否存在节点号重复和单元遗漏的简便算法。     

三维有限元六面体网格几何自适应再生成方法

基于几何自适应的六面体网再生成方法,提出了一种以栅格法为基本方法,对有限元网格再划分过程中的网格再划分标准、几何形状的获取以及新旧网格之间物理量的传递等关键问题进行了研究.重点介绍了基于几何自适应的六面体网格再生成算法,首先对旧网格实体模型进行识别,根据实体模型几何特征建立加密源点信息场;然后采用栅格法生成核心网格并对...     

覆盖件冲压仿真中的网格形态优化技术

覆盖件冲压仿真计算模型中网格密度分布的合理性与网格单元形态的优劣,对仿真结果的准确性有很大的影响.提出了一种改进的全四边形网格细分方法,使网格的密度分布适于覆盖件冲压分析计算要求,可保证细分后网格的协调性,并将算法推广以处理非结构化四边形网格和三角形四边形混合网格的细分.提出的网格细分策略,有助于提高细分后网格的质量.提出了适用于细分后四边形网格和非结构四边形网格的拓扑形态优化操作,可有效的提高网格模型的形态质量.     

一类多线程并行四面体网格优化算法

提出了一类多线程并行、组合了光滑化与局部重连操作的四面体网格优化算法。采用传统的基于拓扑的数据分解策略实现了并行光滑化算法,利用图染色算法将待光滑化的点分解成多个独立点集。同时提出将一类基于几何的数据分解策略应用于局部重连操作的并行化,在每个局部重连操作涉及的几何区域中定义一个特征点,然后沿希尔伯特(Hilbert)曲线对特征点进行排序,曲线的均匀分解对应局部重连操作在各线程的分配。这一分配策略的优点是使并行执行局部重连操作时重连区域相互干涉的情形极少出现。因此,当干涉情形出现时,可选择放弃产生干涉的操作,并行优化效率和效果并无明显的负面影响。最后,数值实验验证了本文算法的效率和有效性。