A dual element based geometric element transformation method for all-hexahedral mesh smoothing

Due to their increased complexity hexahedral elements are more challenging with respect to mesh generation and mesh improvement techniques than tetrahedral elements. In particular, there is a lack of geometry-based all-hexahedral smoothing methods for mesh quality improvement being easy to implement, practicable, and efficient. The recently introduced geometric element transformation method represents a new promising element oriented smoothing concept to resolve this deficiency. By giving a dual octahedron based regularizing transformation this new approach is adapted in order to smooth all-hexahedral meshes. First numerical tests indicate that the resulting smoothing method yields high quality results at least comparable to those of a state of the art global optimization-based approach while being significantly faster.     

Hexahedral mesh smoothing via local element regularization and global mesh optimization

The quality of finite element meshes is one of the key factors that affect the accuracy and reliability of finite element analysis results. In order to improve the quality of hexahedral meshes, we present a novel hexahedral mesh smoothing algorithm which combines a local regularization for each hexahedral mesh, using dual element based geometric transformation, with a global optimization operator for all hexahedral meshes. The global optimization operator is composed of three main terms, including the volumetric Laplacian operator of hexahedral meshes and the geometric constraints of surface meshes which keep the volumetric details and the surface details, and another is the transformed node displacements condition which maintains the regularity of all elements. The global optimization operator is formulated as a quadratic optimization problem, which is easily solved by solving a sparse linear system. Several experimental results are presented to demonstrate that our method obtains higher quality results than other state-of-the-art approaches.     

多边形单元网格自动生成技术

近年来兴起的多边形有限元方法,在有限元计算中采用多边形单元划分网格,不仅可以更好地适应求解区域的几何形状,而且增加了网格划分的灵活性。为了更方便有效地生成多边形单元网格,在Delaunay三角形的基础上,通过将共圆Delaunay三角形合并为一个圆内接多边形,首先提出了Delaunay多边形的概念,进而提出了一种多边形网格自动生成的Delaunay多边形化算法。利用该Delaunay多边形化技术,对工程中常见的几何形状进行网格划分的具体算例表明,Delaunay多边形化方法可以生成性质优良的多边形单元网格。     

An Efficient Optimization Procedure for Tetrahedral Meshes by Chaos Search Algorithm

A simple and efficient local optimization-based procedure for node reposition-ing/smoothing of three-dimensional tetrahedral meshes is presented.The initial tetrahedral mesh is optimized with respect to a specified element shape measure by chaos search algorithm,which is very effective for the optimization problems with only a few design variables.Examples show that the presented smoothing procedure can provide favorable conditions for local transformation approach and the quality of mesh can be significantly improved by the combination of these two procedures with respect to a specified element shape measure.Meanwhile,several commonly used shape measures for tetrahedral element,which are considered to be equivalent in some weak sense over a long period of time,are briefly re-examined in this paper.Preliminary study indicates that using different measures to evaluate the change of element shape will probably lead to inconsistent result for both well shaped and poorly shaped elements.The proposed smoothing approach can be utilized as an appropriate and effective tool for evaluating element shape measures and their influence on mesh optimization process and optimal solution.     

面向法向域网格建模的非线性引导滤波

网格建模是数字几何处理领域的基础性研究问题.为了提高网格建模的简便性和鲁棒性,首先提出了一种非线性的引导滤波算法.滤波过程在法向域进行,滤波后的法向是引导网格法向的局部二次变换;然后,应用上述算法研究了建模方面的2个重要问题:网格去噪和网格平滑,其中的难点在于如何构造合适的引导网格.针对去噪问题,每次迭代时利用双边法向滤波得到引导网格;针对平滑问题,引导网格以高斯滤波结果作为初始值,进而结合原始网格不断进行更新;最后,在形状复杂或特征丰富的网格模型上进行了去噪、平滑等实验,结果表明,该算法简单实用、鲁棒,去噪时能够有效地去除强噪声,保持模型的几何特征;平滑时能够提取出中小尺度的特征,保留大尺度的特征.     

基于LEPP递变的四面体网格自适应剖分算法

为了更合理地进行四面体网格剖分,提出了一种根据待剖分对象形态不同进行网格密度自适应调整的四面体网格剖分方法。该方法首先采用BCC(body-centered cubic)网格初始化网格空间,并根据表面曲率的大小以及距离物体表面的远近,采用LEPP(longest edge propagation path)算法由外至内对初始化后的网格空间进行不同尺度的细分;然后对横跨表面的网格进行调整,以形成对象的表面形态;最后采用以质量函数引导的拉普拉斯平滑与棱边收缩(edge collapse)的方法对网格的质量进行优化来最终得到待剖分对象的四面体网格。结果表明,该方法所生成的网格不仅具有自适应的网格密度,而且网格质量比常用的Advancing Front算法也有所提高。对于基于3维断层图像或表面模型进行有限元建模,该方法不失为一种行之有效的好方法。     

三维实体仿真建模的网格自动生成方法

有限元网格模型的生成与几何拓扑特征和力学特性有直接关系。建立网格模型时，为了更真实地反映原几何形体的特征，在小特征尺寸或曲率较大等局部区域网格应加密剖分；为提高有限元分析精度和效率，在待分析的开口、裂纹、几何突变、外载、约束等具有应力集中力学特性的局部区域，网格应加密剖分。为此，该文提出了基于几何特征和物理特性相结合的网格自动生成方法。该方法既能有效地描述几何形体，又能实现应力集中区域的网格局部加密及粗细网格的均匀过渡。实例表明本方法实用性强、效果良好。     

一个高效可靠的三维AFT四面体网格生成算法

针对三维推进波前算法(AFT-Advancing Front Technique)存在的效率与收敛性问题,文中提出了一整套改进方案,给出了基于拓扑连接的网格数据结构和基于Hash表的网格元素的插入、查找、删除算法,提高了整个算法的效率.通过在网格生成过程中动态维护前沿的尺寸信息,提高四面体单元的整体质量.在内核回退求解时通过引入前沿优先因子,改变前沿推进的路径,大大增加了回退求解的成功概率;对于极少数不能回退求解的内核采用基于线性规划的插点方法加以解决,这样就基本保证了整个算法的收敛.在网格生成以后,通过删除不必要的内部节点、合并相关四面体单元以及对所有内部节点进行基于角度的优化,从而进一步有效提高了网格质量.数值算例表明,文中提出的改进算法具有接近线性的时间复杂度,生成网格质量好.该算法已经得到工程应用.     

基于拉普拉斯算子的三角网格模型的平滑与压缩算法

三角网格模型被广泛应用于各个领域并迅速发展,为了既保留网格模型的局部几何特征,又在平滑去噪的同时能够较好地保持边缘、纹理等细节信息,提出一种基于拉普拉斯算子的偏微分方程平滑方法,可以得到较好的去噪效果;为了更加方便三维数据的传输与操作,采用了一种在拉普拉斯算子的基础上,对三角网格模型进行特征分解进而进行光谱压缩的方法,可以实现对模型的压缩.     

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

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

Block Cartesian abstraction of a geometric model and its application in hexahedral mesh generation

In this work, a fuzzy logic approach is proposed to transform a geometric model of arbitrary shape to its block Cartesian abstraction. This abstraction is topologically similar to the original model and it contains geometric sub-entities which are all aligned in the Cartesian directions. This is achieved by calculating the modifications made to the face normal vectors as a result of the influences of the adjacent faces. A fuzzy logic inference engine is developed by combining heuristics to emulate the local changes in face normal vectors with respect to the changes in the global space. A three-dimensional field morphing algorithm is used to position the features of this block Cartesian abstraction so that a congruent geometric model can be reconstructed. Such a model is useful for the generation of structured quadrilateral boundary element meshes or structured hexahedral meshes based on grid-based meshing method, mesh mapping or sweeping. This approach is also able to overcome the traditional problem of having poorly shaped elements at the boundary using the grid-based method of mesh generation. As the topology of the block Cartesian abstraction is congruent to the original model, the mesh can be mapped back to the original model by employing an inverse operation of the transformation.     

A cascadic geometric filtering approach to subdivision

A new approach to subdivision based on the evolution of surfaces under curvature motion is presented. Such an evolution can be understood as a natural geometric filter process where time corresponds to the filter width. Thus, subdivision can be interpreted as the application of a geometric filter on an initial surface. The concrete scheme is a model of such a filtering based on a successively improved spatial approximation starting with some initial coarse mesh and leading to a smooth limit surface.

In every subdivision step the underlying grid is refined by some regular refinement rule and a linear finite element problem is either solved exactly or, especially on fine grid levels, one confines to a small number of smoothing steps within the corresponding iterative linear solver. The approach closely connects subdivision to surface fairing concerning the geometric smoothing and to cascadic multigrid methods with respect to the actual numerical procedure. The derived method does not distinguish between different valences of nodes nor between different mesh refinement types. Furthermore, the method comes along with a new approach for the theoretical treatment of subdivision.     

基于混合频谱信号编码的网格纹理平滑

网格纹理平滑技术要求既能保持模型大尺度结构特征又能去除模型小尺度纹理.然而当模型小尺度纹理与噪声相差较大时,大多数网格光顺算法会将网格纹理识别为特征加以保持,而无法有效将其去除;现有的基于谱分析的网格光顺方法尽管能有效去除网格纹理,但又无法同时保持模型大尺度结构特征.为解决该问题,本文提出一种基于混合频谱信号编码的低通过滤网格纹理平滑算法.首先采用基于视觉感知的特征识别方法,准确区分模型大尺度与小尺度特征.然后,基于顶点特征尺度,采用差异性频谱信号编码的方式进行几何信息重建,最终实现在保持网格模型大尺度结构特征的同时有效去除小尺度纹理.算法解决了现有网格光顺方法在模型小尺度纹理特征与噪声有明显区别时,无法有效去除纹理的问题;并且也解决了现有基于谱分析的网格光顺方法无法在去除模型小尺度纹理的同时,又能保持模型大尺度特征的矛盾.实验结果验证了算法的有效性.     

基于邻域相似性的三角网格光顺算法

为了提高三角网格模型的质量,满足模型后续处理的要求,提出了一种基于邻域相似性的网格光顺算法。类比图像像素的灰度值,首先构造顶点的双边滤波微分算子,作为其几何灰度值;然后计算顶点邻域之间的相似性,作为顶点几何灰度值的权值,并对顶点的邻域顶点的几何灰度值进行加权平均,得到该顶点的最终几何灰度值;最后将顶点沿着其法矢量方向移动几何灰度值大小的距离,得到光顺后的三角网格模型。实验证明,该算法在光顺模型的同时有效地保持了网格的几何特征。     

Adaptive generation of hexahedral element mesh using an improved grid-based method

An improved grid-based algorithm for the adaptive generation of hexahedral finite element mesh is presented in this paper. It is named as the inside-out grid-based method and involves the following four steps. The first step is the generation of an initial grid structure which envelopes the analyzed solid model completely. And the elements size and density maps are constructed based on the surface curvature and local thickness of the solid model. Secondly, the core mesh is generated through removing all the undesired elements using even and odd parity rules. The third step is to magnify the core mesh in an inside-out manner through a surface node projection process using the closest position approach. To match the mesh to the characteristic boundary of the solid model, a minimal Scaled Jacobian criterion is employed. Finally, in order to handle the degenerated elements and improve the quality of the resulting mesh, two comprehensive techniques are employed: the insertion technique and collapsing technique. The present method was applied in the mesh construction of different engineering problems. Scaled Jacobian and Skew metrics are used to evaluate the hexahedral element mesh quality. The application results show that all-hexahedral element meshes which are well-shaped and capture all the geometric features of the original solid models can be generated using the inside-out grid-based method presented in this paper.     

A New Approach to Fully Automatic Mesh Generation

Automatic mesh generation is one of the most important parts in CIMS (Computer Integrated Manufacturing System).A method based on mesh grading propagation which automatically produces a triangular mesh in a multiply connected planar region is presented in this paper.The method decomposes the planar region into convex subregions,using algorithms which run in linear time.For every subregion,an algorithm is used to generate shrinking polygons according to boundary gradings and form delaunay triangulation between two adjacent shrinking polygons,both in linear time.It automatically propagates boundary gradings into the interior of the region and produces satisfactory quasi-uniform mesh.     

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

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

隐式曲面的多分辨率法向网格逼近

法向网格是一种新型的曲面多分辨率描述方式，其中每个层次都可以表示为其前一个粗糙层次的法向偏移．文中提出一种基于法向网格表示的隐式曲面多分辨率网格逼近算法．首先通过基于空间剖分技术的多边形化算法获得隐式曲面的粗糙逼近网格，并利用网格均衡化方法对粗糙网格进行优化，消除其中的狭长三角形；然后利用法向细分规则迭代地对网格中的三角面片进行细分，并利用区间算术技术沿法向方向对隐式曲面进行逼近．最终生成的隐式曲面分片线性逼近网格为法向网格．该逼近网格为隐式曲面提供了一种多分辨率表示，网格具有细分连通性，其数据量较传统的多边形化算法所生成的网格有大幅度的压缩．该算法可用于隐式曲面的多级绘制、累进传输及相关数字几何处理．     

多孔平面的快速边界元划分

在 3D VL SI互连寄生电容的边界元素法计算中 ,多孔平面的边界元划分是十分困难的问题 .文中提出一种快速划分多孔平面边界元的方法 ,它可高效处理非正交几何边界形状 ,形成规则的梯形元 .与全局扫描线法相比 ,有较高的划分速度、计算速度与精度