三维有限元网格尺寸场光滑化的优化模型和算法

针对单元尺寸值过渡剧烈会导致有限元网格包含低质量单元的问题,提出基于优化原理的单元尺寸场光滑化理论及对应的几何自适应四面体网格生成算法.首先输入CAD模型,生成一套覆盖模型内部的非结构背景网格;然后结合用户参数计算背景网格点上的曲率和邻近特征,以获得自适应CAD模型几何特征的初始单元尺寸场;再以最小化初始单元尺寸场的改变为目标,以单元尺寸值过渡受控为约束,通过求解一类凸优化问题光滑初始尺寸场;最后以光滑后的尺寸场为输入,先后在CAD模型表面与内部生成曲面网格和实体网格.实验结果表明,文中算法仅需5个用户参数,即可在给定CAD模型内部全自动生成高质量的四面体网格.     

二维自适应前沿推进网格生成

针对二维平面问题,通过曲率计算和基于中轴理论的邻近特征计算控制区域边界曲线的离散;修改经典的前沿推进算法,利用边界驱动的单元尺寸控制方式在区域内部布置疏密过渡合理的三角网格;结合几何和拓扑策略提升网格质量。实验表明,上述算法可生成单元质量高、尺寸过渡合理的计算网格。     

一种网格和节点同步生成的二维Delaunay网格划分算法

应用Lawson算法对网格的Delaunay性质进行维护,利用单元尺度场控制生成网格的疏密分布;找到任一不满足尺度场要求的单元,在其可插度最大的边上按一定法则插入新节点,加密网格,实现内节点的生成与网格划分同步进行.该算法避免了搜寻包含三角形的过程,提高了效率.通过多次划分实验表明,该算法的时间复杂度约为O(N1.2).同时,由于在不满足单元尺寸要求的单元边上插入新节点,直接对单元的边长进行控制,使得网格的质量和自适性更加良好.     

三维实体有限元自适应网格规划生成

为实现三维实体有限元网格自适应生成,设计了中心点、沿指定曲线和基于实体表面等网格加密生成方式;并根据分析对象几何特征和物理特性经验估计,以规划的方式构造自适应网格单元尺寸信息场．在此基础上,提出基于Delaunay剖分的动态节点单元一体化算法,生成几何特征和物理特性整体自适应的有限元网格．     

可控的高度规整三角网格生成算法

本文提出了一种新的用户可控的高度规整三角网格生成算法.通过在网格表面上构造3个标量场,利用其等值线相交生成高度规整的三角网格.算法借助N-对称方向场来指导生成网格的边方向,在网格表面指定密度场来控制采样密度,同时还提供了特征对齐和对带边界模型的处理能力.所有的控制需求都被纳入标量场求解框架中统一优化.实验表明,本文的方法能够满足多种用户控制需求,生成高度规整的三角网格.     

二维几何特征自适应有限元网格生成(二)--算法描述

以Delaunay三角剖分为基础,构造几何特征自适应有限元网格单元尺寸信息场,给出动态节点一单元一体化生成算法,实现二维形体几何特征自适应有限元网格的自动生成,并使分析对象力学特性得到一定程度的自适应．     

基于特征识别的多尺度构件网格自动生成算法

随着结构力学领域待解决问题复杂程度不断提高,多尺度构件的高质量网格生成对于其数值模拟的计算精度起着至关重要的作用。本文提出一种基于特征识别的网格自动生成技术方法,该方法将多尺度构件的不同量级尺度几何特征识别出来,根据其不同尺度尺寸设置相关区域的网格尺寸值,利用Delaunay三角化算法和前沿推进法生成能够反映不同尺度几何特征的网格单元,再对小尺度区域周围进行加密处理,最后通过几何指数控制函数将不同尺度网格过渡连接起来,形成多尺度构件的整体网格划分模型。通过2个几何模型的测试表明该方法生成的整体网格质量好,不同尺度区域网格过渡合理,自动化程度较高。     

任意多介质区域变尺寸有限元网格自动剖分

以Bowyer／watson算法为基础，结合Lawson算法等对有限元网格进行Delaunay剖分，引入节点间距函数，在区域内通过布置约束点和边来控制特定位置网格的尺寸，并实现网格的尺寸均匀过渡．采用了多种优化技术对生成的网格进行优化。对Bowyer／Watson算法和Delaunay剖分的一些缺陷及问题作了探讨，并提出了相应的处理办法。编写的网格自动剖分软件具有良好。的可视化界面，其图形生成、显示及修改功能十分完善，极大地提高了有限元网格剖分的效率和准确性.     

基于改进波前法的曲面网格生成算法

为提高传统波前法(Advancing Front Method,AFM)的网格生成效率,利用多维搜索二叉树数据结构实现临近前沿和节点的快速查找,使整个网格生成的时间复杂度接近线性.针对周期曲面网格的生成,提出2种点修正算子,避免传统算法添加虚边界导致局部网格单元质量较差和虚边界计算复杂的问题.网格生成实例表明:多维搜索二叉树提高网格生成速度,引进点修正算子的波前法改善周期曲面网格质量.     

一种快速的可逆累进网格算法

在分析已有累进网格生成算法的基础上,构造了一种新的网格简化信息记录表示法,并提出一种基于“边折叠”网格简化方法的累进网格生成算法。此算法不仅消除了累进网格技术中的二义性,而且能够较大地提高累进网格的运算速度。     

一种快速的三维有限元网格生成方法

针对三维有限元网格的生成的速度较慢并且网格质量不高的问题,提出了一种基于约束波前法的三维有限元网格生成算法。算法的主要思想是用背景网格提高网格单元的可控性,避免网格单元生成时验证有效性的计算量,从而快速生成高质量的三维有限元网格。算法首先借助八叉树方法生成背景网格,其次利用背景网格的密度对模型表面进行三角剖分得到初始波前,然后依据背景网格的特征生成实体网格单元,最后对得到的结果进行优化。实验证明结合了八叉树和推进波前法的三维网格生成算法降低了波前法的时间复杂度,将其效率提高了20%,而且能得到更高质量的网格。     

基于Delaunay四面体剖分的网格分割算法

为了构建有意义曲面分片,提出一种基于Delaunay四面体剖分的网格分割算法.首先根据Delaunay四面体剖分得到多边形网格内部的四面体,求出每个面上反映网格内部信息的Delaunay体距离;然后对Delaunay体距离进行平滑处理,再对网格上面的Delaunay体距离进行聚类,用高斯混合模型对Delaunay体距离作柱状图的拟合,利用期望最大化算法来快速求得拟合结果;最后结合图切分技术,同时考虑聚类的结果、分割区域的边界平滑和视觉认知中的最小规则,得到最终的网格分割结果.实验结果表明,采用文中算法可以有效地实现有意义的网格分割.     

Adaptive remeshing process with quadrangular finite elements

Since the quality of FEM analysis directly depends on the quality of meshes, various mesh adaptation schemes have been researched. There are two stages on adaptive finite element analysis; to derive error measure and to control meshes based on error measure. The former has been well researched among applied mathematicians. However, the importance of the latter aspect wasn't considered enough. Even if the error measures were well estimated, the total performance of mesh adaptation might be poor with a poor mesh control. This paper proposes an effective mesh control scheme for h-adaptation, or adaptive remeshing scheme with the explicit relation between interpolation theory based on error measure and desirable mesh size. Total mesh adaptation is controlled by introducing Quality Index, or the ratio between the total error norm and the total energy norm which represents the quality of the total meshes; specifying the desirable value of Quality Index, then the adaptive remeshing process can handle it and Quality Index is almost converged to the given value. Since the full automatic feature of the mesh generator is a prerequisite for adaptive remeshing, the author also discusses the algorithm of the quadrangular mesh generator for arbitrary domains. After evaluation on a linear problem, it's confirmed that the proposed mesh control scheme and the proposed error measure-mesh size relations are acceptable. The incompatible case for mesh adaptation is also discussed in this paper.     

从表面重构的体数据实现3维网格剖分

目的 针对有限元分析中网格最优化问题,提出一种改进的生成四面体网格的自组织算法。方法 该算法首先应用几何方法将三角形表面模型重新构造成规定大小的分类体数据,同时由该表面模型建立平衡八叉树,计算用以控制网格尺寸的3维数组;然后将体数据转换成邻域内不同等值面的形态一致的边界指示数组;结合改进的自组织算法和相关3维数据的插值函数,达到生成四面体网格的目的。结果 实验结果对比表明,该方法能够生成更高比例的优质四面体,增强了对扁平面体的抑制能力,同时很好地保证了边界的一致。结论 在对封闭的3维表面网格进行有限元建模时,本文算法为其提供了一种有效、可靠的途径。     

边界简化与多目标优化相结合的高质量四边形网格生成

目的 高质量四边形网格生成是计算机辅助设计、等几何分析与图形学领域中一个富有挑战性的重要问题。针对这一问题,提出一种基于边界简化与多目标优化的高质量四边形网格生成新框架。方法 首先针对亏格非零的平面区域,提出一种将多连通区域转化为单连通区域的方法,可生成高质量的插入边界;其次,提出"可简化角度"和"可简化面积比率"两个阈值概念,从顶点夹角和顶点三角形面积入手,将给定的多边形边界简化为粗糙多边形;然后对边界简化得到的粗糙多边形进行子域分解,并确定每个子域内的网格顶点连接信息;最后提出四边形网格的均匀性和正交性度量目标函数,并通过多目标非线性优化技术确定网格内部顶点的几何位置。结果 在同样的离散边界下,本文方法与现有方法所生成的四边网格相比,所生成的四边网格顶点和单元总数目较少,网格单元质量基本类似,计算时间成本大致相同,但奇异点数目可减少70% 80%,衡量网格单元质量的比例雅克比值等相关指标均有所提高。结论 本文所提出的四边形网格生成方法能够有效减少网格中的奇异点数目,并可生成具有良好光滑性、均匀性和正交性的高质量四边形网格,非常适用于工程分析和动画仿真。     

Wavelet-based multiresolution analysis of irregular surface meshes

We extend Lounsbery's multiresolution analysis wavelet-based theory for triangular 3D meshes, which can only be applied to regularly subdivided meshes and thus involves a remeshing of the existing 3D data. Based on a new irregular subdivision scheme, the proposed algorithm can be applied directly to irregular meshes, which can be very interesting when one wants to keep the connectivity and geometry of the processed mesh completely unchanged. This is very convenient in CAD (computer-assisted design), when the mesh has attributes such as texture and color information, or when the 3D mesh is used for simulations, and where a different connectivity could lead to simulation errors. The algorithm faces an inverse problem for which a solution is proposed. For each level of resolution, the simplification is processed in order to keep the mesh as regular as possible. In addition, a geometric criterion is used to keep the geometry of the approximations as close as possible to the original mesh. Several examples on various reference meshes are shown to prove the efficiency of our proposal.     

从表面重构的体数据实现三维网格剖分

针对有限元分析中网格最优化问题,本文提出一种改进的生成四面体网格的自组织算法。该算法首先应用几何方法将三角形表面模型重新构造成规定大小的分类体数据,同时由该表面模型建立平衡八叉树,计算用以控制网格尺寸的三维数组;然后将体数据转换成邻域内不同等值面的形态一致的边界指示数组;结合改进的自组织算法和相关三维数据的插值函数,达到生成四面体网格的目的。实验对比表明,该方法能够生成更高比例的优质四面体,同时很好地保证了边界的一致。在对封闭的三维表面网格进行有限元建模时,本文算法为其提供了一种有效、可靠的途径。     

All-hexahedral mesh generation via inside-out advancing front based on harmonic fields

The generation of hexahedral meshes is an open problem that has undergone significant research. This paper deals with a novel inside-out advancing front method to generate unstructured all-hexahedral meshes for given volumes. Two orthogonal harmonic fields, principal and radial harmonic fields, are generated to guide the inside-out advancing front process based on a few user interactions. Starting from an initial hexahedral mesh inside the given volume, we advance the boundary quadrilateral mesh along the streamlines of radial field and construct layers of hexahedral elements. To ensure high quality and uniform size of the hexahedral mesh, quadrilateral elements are decomposed in such a way that no non-hexahedral element is produced. For complex volume with branch structures, we segment the complex volume into simple sub-volumes that are suitable for our method. Experimental results show that our method generates high quality all-hexahedral meshes for the given volumes.     

A Parallel Implementation of ALE Moving Mesh Technique for FSI Problems using OpenMP

This paper investigates a high performance implementation of an Arbitrary Lagrangian Eulerian moving mesh technique on shared memory systems using OpenMP environment. Moving mesh techniques are considered an integral part of a wider class of fluid mechanics problems that involve moving and deforming spatial domains, namely, free-surface flows and Fluid Structure Interaction (FSI). The moving mesh technique adopted in this work is based on the notion of nodes relocation, subjected to a certain evolution as well as constraint conditions. A conjugate gradient method augmented with preconditioning is employed for solution of the resulting system of equations. The proposed algorithm, initially, reorders the mesh using an efficient divide and conquer approach and then parallelizes the ALE moving mesh scheme. Numerical simulations are conducted on the multicore AMD Opteron and Intel Xeon processors, and unstructured triangular and tetrahedral meshes are used for the 2D and 3D problems. The quality of generated meshes is checked by comparing the element Jacobians in the reference and current meshes, and by keeping track of the change in the interior angles in triangles and tetrahedrons. Overall, 51 and 72% efficiencies in terms of speedup are achieved for both the parallel mesh reordering and ALE moving mesh algorithms, respectively.     

Anisotropic quadrilateral mesh generation: An indirect approach

In this paper a new indirect approach is presented for anisotropic quadrilateral mesh generation based on discrete surfaces. The ability to generate grids automatically had a pervasive influence on many application areas in particularly in the field of Computational Fluid Dynamics. In spite of considerable advances in automatic grid generation there is still potential for better performance and higher element quality. The aim is to generate meshes with less elements which fit some anisotropy criterion to satisfy numerical accuracy while reducing processing times remarkably. The generation of high quality volume meshes using an advancing front algorithm relies heavily on a well designed surface mesh. For this reason this paper presents a new technique for the generation of high quality surface meshes containing a significantly reduced number of elements. This is achieved by creating quadrilateral meshes that include anisotropic elements along a source of anisotropy.