Automatic finite element mesh generation for maxillary second premolar.

Developing three dimensional finite element mesh models for irregular geometric objects requires a large amount of manual efforts, hence limiting the three dimensional approach for dental structure analyses. An automatic procedure which can be used to generate a three dimensional finite element mesh for the maxillary second premolar was developed in this study. Firstly, a embedded second premolar was sliced and scanned parallel to the occlusal surface. A self-developed image processing system was employed to detect the boundaries of different materials within each section. An automatic mesh generation program was used on these boundaries to create tetrahedral elements based on moving nodes of uniform cube approach. Six mesh models of the second premolar with different element sizes using linear and quadratic elements were analyzed. Strain energy and von Mises stresses were reviewed for convergence in the crown regions.     

Feature-based modeling for automatic mesh generation

Automatic meshing algorithms for finite element analysis are based on a computer understanding of the geometry of the part to be discretized. Current mesh generators understand the part as either a boundary representation, an octree, or a point set. A higher-level understanding of the part can be achieved by associating engineering significance and engineering data, such as loading and boundary conditions, with generic shapes in the part. This technique, called feature-based modeling, is a popular approach to integrating computer-aided design (CAD) and computer-aided manufacturing through the use of machinable shapes in the CAD model. It would seem that feature-based design also could aid in the finite element mesh generation process by making engineering information explicit in the model.This paper describes an approach to feature-based mesh generation. The feature representation of a fully functioning feature-based system that does automatic process planning and inspection was extended to include finite element mesh generation. This approach is based on a single feature representation that can be used for design, finite element analysis, process planning, and inspection of prismatic parts. The paper describes several advantages that features provide to the meshing process, such as improved point sets and a convenient method of simplifying the geometry of the model. Also discussed are possible extensions to features to enhance the finite element meshing process.     

Constrained mesh optimization on boundary

This paper presents a new mesh optimization approach aiming to improve the mesh quality on the boundary. The existing mesh untangling and smoothing algorithms (Vachal et al. in J Comput Phys 196: 627–644, 2004; Knupp in J Numer Methods Eng 48: 1165–1185, 2002), which have been proved to work well to interior mesh optimization, are enhanced by adding constrains of surface and curve shape functions that approximate the boundary geometry from the finite element mesh. The enhanced constrained optimization guarantees that the boundary nodes to be optimized always move on the approximated boundary. A dual-grid hexahedral meshing method is used to generate sample meshes for testing the proposed mesh optimization approach. As complementary treatments to the mesh optimization, appropriate mesh topology modifications, including buffering element insertion and local mesh refinement, are performed in order to eliminate concave and distorted elements on the boundary. Finally, the optimization results of some examples are given to demonstrate the effectivity of the proposed approach.     

Generation of finite element mesh with variable size over an unbounded 2D domain

With the advance of the finite element, general fluid dynamic and traffic flow problems with arbitrary boundary definition over an unbounded domain are tackled. This paper describes an algorithm for the generation of finite element mesh of variable element size over an unbounded 2D domain by using the advancing front circle packing technique. Unlike the conventional frontal method, the procedure does not start from the object boundary but starts from a convenient point within an open domain. The sequence of construction of the packing circles is determined by the shortest distance from the fictitious centre in such a way that the generation front is more or less a circular loop with occasional minor concave parts due to element size variation. As soon as a circle is added to the generation front, finite elements are directly generated by properly connecting frontal segments with the centre of the new circle. In contrast to other mesh generation schemes, the domain boundary is not considered in the process of circle packing, this reduces a lot of geometrical checks for intersection with frontal segments, and a linear time complexity for mesh generation can be achieved. In case the boundary of the domain is needed, simply generate an unbounded mesh to cover the entire object. As the element adjacency relationship of the mesh has already been established in the circle packing process, insertion of boundary segments by neighbour tracing is fast and robust. Details of such a boundary recovery procedure are described, and practical meshing problems are given to demonstrate how physical objects are meshed by the unbounded meshing scheme followed by the insertion of domain boundaries.     

An algorithm for automatic 2D finite element mesh generation with line constraints

In this study, an algorithm is designed specifically for automatic finite element (FE) mesh generation on the transverse structure of hulls reinforced by stiffeners. Stiffeners attached to the transverse structure are considered as line constraints in the geometry boundary. For the FE mesh generation used in this study, the line constraints are treated as boundaries and by that means the geometry domain attached to the line constraints is decomposed into sub-domains, constrained only by the closed boundaries. Then, the mesh can be generated directly on those sub-domains by the traditional approach. The performance of the proposed algorithm is evaluated and the quality of the generated mesh meets expectations.     

Asymptotic boundary conditions with immersed finite elements for interface magnetostatic/electrostatic field problems with open boundary

Many of the magnetostatic/electrostatic field problems encountered in aerospace engineering, such as plasma sheath simulation and ion neutralization process in space, are not confined to finite domain and non-interface problems, but characterized as open boundary and interface problems. Asymptotic boundary conditions (ABC) and immersed finite elements (IFE) are relatively new tools to handle open boundaries and interface problems respectively. Compared with the traditional truncation approach, asymptotic boundary conditions need a much smaller domain to achieve the same accuracy. When regular finite element methods are applied to an interface problem, it is necessary to use a body-fitting mesh in order to obtain the optimal convergence rate. However, immersed finite elements possess the same optimal convergence rate on a Cartesian mesh, which is critical to many applications. This paper applies immersed finite element methods and asymptotic boundary conditions to solve an interface problem arising from electric field simulation in composite materials with open boundary. Numerical examples are provided to demonstrate the high global accuracy of the IFE method with ABC based on Cartesian meshes, especially around both interface and boundary. This algorithm uses a much smaller domain than the truncation approach in order to achieve the same accuracy.     

Numerical comparison between two cost functions in contact shape optimization

A finite element approach to shape optimization in a 2D frictionless contact problem for two different cost functions is presented in this work. The goal is to find an appropriate shape for the contact boundary, performing an almost constant contact-stress distribution. The whole formulation, including the mathematical model for the unilateral problem, sensitivity analysis and geometry definition is treated in a continuous form, independently of the discretization in finite elements. Shape optimization is performed by a direct modification of the geometry throughB-spline curves and an automatic mesh generator is used at each new configuration to provide the finite element input data. Augmented-Lagrangian techniques (to solve the contact problem) and an interior-point mathematical-programming algorithm (for shape optimization) are used to obtain numerical results.     

A mesh generation method based on set theory

A method is proposed to study automatic generation of input data to finite element programs by the computer with minimal human involvement. Based on a set of objective criteria, a consistent procedure is developed to scan through the geometric idealization of a boundary to form regions for mesh generation. This procedure is supported by a rigorously built mathematical model. A cyclic set is defined according to the ordering of segments on a boundary. Corresponding to this is a quadrilateral set whose elements represent the partitioning of the geometry into regions for mesh generation. Theorems are proven establishing the uniqueness and completeness of this procedure. By means of this mathematical model, a computer program is written for automatic mesh generation together with evaluation of results by a given set of objective criteria.Examples are presented with simply connected and with multiply connected geometry.     

Eulerian shape design sensitivity analysis and optimization with a fixed grid

Conventional shape optimization based on the finite element method uses Lagrangian representation in which the finite element mesh moves according to shape change, while modern topology optimization uses Eulerian representation. In this paper, an approach to shape optimization using Eulerian representation such that the mesh distortion problem in the conventional approach can be resolved is proposed. A continuum geometric model is defined on the fixed grid of finite elements. An active set of finite elements that defines the discrete domain is determined using a procedure similar to topology optimization, in which each element has a unique shape density. The shape design parameter that is defined on the geometric model is transformed into the corresponding shape density variation of the boundary elements. Using this transformation, it has been shown that the shape design problem can be treated as a parameter design problem, which is a much easier method than the former. A detailed derivation of how the shape design velocity field can be converted into the shape density variation is presented along with sensitivity calculation. Very efficient sensitivity coefficients are calculated by integrating only those elements that belong to the structural boundary. The accuracy of the sensitivity information is compared with that derived by the finite difference method with excellent agreement. Two design optimization problems are presented to show the feasibility of the proposed design approach.     

A new front updating solution applied to some engineering problems

Summary An automatic mesh generation dealing with domains of an arbitrary shape could be realized by an advancing front method. The mesh generator based on this method creates triangle elements inside a domain starting with the polygonal (polyhedral in 3D) discretisation of its border. In this paper an original algorithm for the front updating procedures as a part of the mesh generator is presented. The proposed algorithm provides an efficient mesh generation procedure. It has been verified on the various domains with complex geometry and with nonuniform distribution of edge nodes such as the discretisation of the switched reluctance motor and power cable configuration, respectively. The related finite element calculations are carried out.     

Adaptive finite element mesh triangulation using self-organizing neural networks

The finite element method is a computationally intensive method. Effective use of the method requires setting up the computational framework in an appropriate manner, which typically requires expertise. The computational cost of generating the mesh may be much lower, comparable, or in some cases higher than the cost associated with the numeric solver of the partial differential equations, depending on the application and the specific numeric scheme at hand.The aim of this paper is to present a mesh generation approach using the application of self-organizing artificial neural networks through adaptive finite element computations. The problem domain is initially constructed using the self-organizing neural networks. This domain is used as the background mesh which forms the input for finite element analysis and from which adaptive parameters are calculated through adaptivity analysis. Subsequently, self-organizing neural network is used again to adjust the location of randomly selected mesh nodes as is the coordinates of all nodes within a certain neighborhood of the chosen node. The adjustment is a movement of the selected nodes toward a specific input point on the mesh. Thus, based on the results obtained from the adaptivity analysis, the movement of nodal points adjusts the element sizes in a way that the concentration of elements will occur in the regions of high stresses. The methods and experiments developed here are for two-dimensional triangular elements but seem naturally extendible to quadrilateral elements.     

A new scheme for mesh generation and mesh refinement using graph theory

There are an extensive number of algorithms available from graph theory, some of which, for problems with geometric content, make graphs an attractive framework in which to model an object from its geometry to its discretization into a finite element mesh. This paper presents a new scheme for finite element mesh generation and mesh refinement using concepts from graph theory. This new technique, which is suitable for an interactive graphical environment, can also be used efficiently for fully automatic remeshing in association with self-adaptive schemes. Problems of mesh refinement around holes and local mesh refinement are treated. The suitability of the algorithms presented in this paper is demonstrated by some examples.     

Generation of a finite element MESH from stereolithography (STL) files

The aim of the method proposed here is to show the possibility of generating adaptive surface meshes suitable for the finite element method, directly from an approximated boundary representation of an object created with CAD software. First, we describe the boundary representation, which is composed of a simple triangulation of the surface of the object. Then we will show how to obtain a conforming size-adapted mesh. The size adaptation is made considering geometrical approximation and with respect to an isotropic size map provided by an error estimator. The mesh can be used “as is” for a finite element computation (with shell elements), or can be used as a surface mesh to initiate a volume meshing algorithm (Delaunay or advancing front). The principle used to generate the mesh is based on the Delaunay method, which is associated with refinement algorithms, and smoothing. Finally, we will show that not using the parametric representation of the geometrical model allows us to override some of the limitations of conventional meshing software that is based on an exact representation of the geometry.     

A hybrid and automated approach to adapt geometry model for CAD/CAE integration

Traditional adaption of CAD geometry, which plays an important role in generating effective and fit-for-purpose finite element models, is usually carried out manually and optionally with excessive dependence on engineer’s experience. Automatic and efficient geometry modification before simulation evidently improves design efficiency and quality, and cuts down product life cycle. This paper represents an automatic approach to generate simplified and idealized geometry models for CAE simulation, which consists of hybrid model simplification criteria, feature-based model simplification, and simulation intent driven geometry modification. Hybrid adaption criteria takes detailed features geometric dimension, geometry topology, design intent into consideration synthetically. Simulation intent-driven modification with the help of virtual topology operators helps to deal with geometry at a higher level to get an ameliorative boundary for mesh without perturbing the original design model with constructing history for down-stream manufacture-oriented application, such as machining feature recognition and process planning. Development of plug-in toolkit guarantees automation of the simplification process and helps generate simulation-fitted geometry for subsequent analysis and simulation process. Prototype system and cases are implemented to demonstrate the result and efficiency of the proposed approach.     

Finite element analysis on implicitly defined domains: An accurate representation based on arbitrary parametric surfaces

In this paper, we present some novel results and ideas for robust and accurate implicit representation of geometric surfaces in finite element analysis. The novel contributions of this paper are threefold: (1) describe and validate a method to represent arbitrary parametric surfaces implicitly; (2) represent arbitrary solids implicitly, including sharp features using level sets and boolean operations; (3) impose arbitrary Dirichlet and Neumann boundary conditions on the resulting implicitly defined boundaries. The methods proposed do not require local refinement of the finite element mesh in regions of high curvature, ensure the independence of the domain’s volume on the mesh, do not rely on boundary regularization, and are well suited to methods based on fixed grids such as the extended finite element method (XFEM). Numerical examples are presented to demonstrate the robustness and effectiveness of the proposed approach and show that it is possible to achieve optimal convergence rates using a fully implicit representation of object boundaries. This approach is one step in the desired direction of tying numerical simulations to computer aided design (CAD), similarly to the isogeometric analysis paradigm.     

A strategy of automatic hexahedral mesh generation by using an improved whisker-weaving method with a surface mesh modification procedure

One of the demands for three dimensional (3D) finite element analyses is the development of an automatic hexahedral mesh generator. For this problem, several methods have been proposed by many researchers. However, reliable automatic hexahedral mesh generation has not been developed at present. In this paper, a new strategy of fully automatic hexahedral mesh generation is proposed. In this strategy, the prerequisite for generating a hexahedral mesh is a quadrilateral surface mesh. From the given surface mesh, combinatorial dual cycles (sheet loops for the whisker-weaving algorithm) are generated to produce a hexahedral mesh. Since generating a good quality hexahedral mesh does not depend only on the quality of quadrilaterals of the surface mesh but also on the quality of the sheet loops generated from it, a surface mesh modification method to remove self-intersections from sheet loops is developed. Next, an automatic hexahedral mesh generator by the improved whisker-weaving algorithm is developed in this paper. By creating elements and nodes on 3D real space during the weaving process, it becomes possible to generate a hexahedral mesh with fewer bad-quality elements. Several examples will be presented to show the validity of the proposed mesh generation strategy.     

Automatic Delaunay mesh generation method and physically-based mesh optimization method on two-dimensional regions

Delaunay mesh generation method is a common method for unstructured mesh (or unstructured grid) generation. Delaunay mesh generation method can conveniently add new points to the existing mesh without remeshing the whole domain. However, the quality of the generated mesh is not high enough if compared with some mesh generation methods. To obtain high-quality mesh, this paper developed an automatic Delaunay mesh generation method and a physically-based mesh optimization method on two-dimensional regions. For the Delaunay mesh generation method, boundary-conforming problem was ensured by create nodes at centroid of mesh elements. The definition of node bubbles and element bubbles was provided to control local mesh coarseness and fineness automatically. For the physically-based mesh optimization method, the positions of boundary node bubbles are predefined, the positions of interior node bubbles are adjusted according to interbubble forces. Size of interior node bubbles is further adjusted according to the size of adjacent node bubbles. Several examples show that high-quality meshes are obtained after mesh optimization.

     

Generation of triangular mesh with specified size by circle packing

This paper describes an algorithm for the generation of a finite element mesh with a specified element size over an unbound 2D domain using the advancing front circle packing technique. Unlike the conventional frontal method, the procedure does not start from the object boundary but starts from a convenient point within the open domain. As soon as a circle is added to the generation front, triangular elements are directly generated by properly connecting frontal segments with the centre of the new circle. Circles are packed closely and in contact with the existing circles by an iterative procedure according to the specified size control function. In contrast to other mesh generation schemes, the domain boundary is not considered in the process of circle packing, this reduces a lot of geometrical checks for intersection between frontal segments. If the mesh generation of a physical object is required, the object boundary can be introduced. The boundary recovery procedure is fast and robust by tracing neighbours of triangular elements. The finite element mesh generated by circle packing can also be used through a mapping process to produce parametric surface meshes of the required characteristics. The sizes of circles in the pack are controlled by the principal surface curvatures. Five examples are given to show the effectiveness and robustness of mesh generation and the application of circle packing to mesh generation over curved surfaces.