The generation of arbitrary order curved meshes for 3D finite element analysis

A procedure for generating curved meshes, suitable for high-order finite element analysis, is described. The strategy adopted is based upon curving a generated initial mesh with planar edges and faces by using a linear elasticity analogy. The analogy employs boundary loads that ensure that nodes representing curved boundaries lie on the true surface. Several examples, in both two and three dimensions, illustrate the performance of the proposed approach, with the quality of the generated meshes being analysed in terms of a distortion measure. The examples chosen involve geometries of particular interest to the computational fluid dynamics community, including anisotropic meshes for complex three dimensional configurations.     

Free surface flow through rock-fill dams analyzed by FEM with level set approach

 A stabilized-finite element formulation is coupled with a level set technique for computations of incompressible non-linear flow with interfaces between two immiscible fluids. An interface capturing formulation (ICF) for non-linear, free surface, seepage flow in rock-fill dams is proposed. The formulation is derived for two- and three-dimensional flow within a fixed mesh domain. The resulting formulation is general and applicable for various steady and transient two-phase flow problems. FE-refinement is processed for the entire fixed mesh domains. A general solver is also reviewed for large and non-symmetric non-positive definite linear system of equations with the GMRES-update technique based on a Newton-iterative method. The computational procedure has been implemented in MATLAB. A comparison is performed between the 2-D computed test problem for coarse and refined meshes together with some proposed analytical solutions for nonlinear seepage flow with free surface in rock-fill dams. An expansion of the 2-D program code to a 3-D one for a rectangular rock-fill dam is also developed and simulated in MATLAB. The performance of the computations in 3-D is very promising and its opening the future for possible industrial applications using the same simple technique. Computations for a simple 3-D seepage flow problem with free surface in rock-fill dam are included in present paper. A general mesh generator and solver for large scale and complex 3-D flow problems in a real embankment dam is also under construction in C++.     

B spline-based method for 2-D large deformation analysis

In this paper, quadratic cardinal B spline functions are used for solution of 2-D large deformation problems. Because the B spline functions are directly used in function approximation, no meshes are needed and the mesh distortion issues in nonlinear analyses are avoided in this method. Using the B spline functions, the solution formulations based on total Lagrangian (TL) approach for two dimensional large deformation problems are established. The numerical examples of 2-D large deformation problems indicate that the B spline method is effective and stable for solving complicated problems.     

Simultaneous untangling and smoothing of moving grids

In this work, a technique for simultaneous untangling and smoothing of meshes is presented. It is based on an extension of an earlier mesh smoothing strategy developed to solve the computational mesh dynamics stage in fluid–structure interaction problems. In moving grid problems, mesh untangling is necessary when element inversion happens as a result of a moving domain boundary. The smoothing strategy, formerly published by the authors, is defined in terms of the minimization of a functional associated with the mesh distortion by using a geometric indicator of the element quality. This functional becomes discontinuous when an element has null volume, making it impossible to obtain a valid mesh from an invalid one. To circumvent this drawback, the functional proposed is transformed in order to guarantee its continuity for the whole space of nodal coordinates, thus achieving the untangling technique. This regularization depends on one parameter, making the recovery of the original functional possible as this parameter tends to 0. This feature is very important: consequently, it is necessary to regularize the functional in order to make the mesh valid; then, it is advisable to use the original functional to make the smoothing optimal. Finally, the simultaneous untangling and smoothing technique is applied to several test cases, including 2D and 3D meshes with simplicial elements. As an additional example, the application of this technique to a mesh generation case is presented. Copyright © 2008 John Wiley & Sons, Ltd.     

AUTOMATIC ADAPTIVE 3-D FINITE ELEMENT REFINEMENT USING DIFFERENT-ORDER TETRAHEDRAL ELEMENTS

Automatic refinement finite element analyses were carried out employing three different-order tetrahedral solid elements for the solution of 3-D stress analysis problems. Numerical results indicated that the adaptive refinement procedure could eliminate effectively the effect of singularities and the optimal convergence rate was achieved in all the examples tested. The preconditioned conjugate gradient technique was used for the solution of the large system of simultaneous equations. By interpolating the initial guess of the iteration solver from the previous converged solution, the number of iterations needed for the solution is lower than expected. Furthermore, when the mesh density distribution pattern has converged, it became even more efficient and independent of the number of degrees of freedom in the finite element mesh. The relative efficiency of the three different-order tetrahedral elements has also been compared in terms of storage and computational cost needed for achieving a certain accuracy. It is found that although the cubic T20 element can achieve the highest convergence rate, the T10 element is the most competitive and effective element in terms of storage and computational cost needed. © 1997 John Wiley & Sons, Ltd.     

Feature based quadrilateral mesh generation for sculptured surface products

A successful approach to the generation of quadrilateral surface meshes for finite element analysis of sculptured surface products is presented. The free-form sculpted surface is divided into a feature based product anatomy which is used as the basis for defining regions of simpler shape that are related by the requirement for mesh continuity across their boundaries. Each surface region is meshed using a paving algorithm. Large surface areas are further subdivided to introduce new boundaries which improve contiol over mesh transitions and element distortion. A procedure, based on this approach, which generates well formed meshes on sculptured surface products is described using a hollow golf club head as an example. The approach is product specific and relies on initial user interaction to populate a feature based product data model. Thus, the data model contains, the attributes of a finite element mesh tailored to the product, which can be used to generate meshes on subsequent design iterations or other members of the product family with a high degree of automation, leading to reduced analysis cost.     

A distortion measure to validate and generate curved high‐order meshes on CAD surfaces with independence of parameterization

A framework to validate and generate curved nodal high‐order meshes on Computer‐Aided Design (CAD) surfaces is presented. The proposed framework is of major interest to generate meshes suitable for thin‐shell and 3D finite element analysis with unstructured high‐order methods. First, we define a distortion (quality) measure for high‐order meshes on parameterized surfaces that we prove to be independent of the surface parameterization. Second, we derive a smoothing and untangling procedure based on the minimization of a regularization of the proposed distortion measure. The minimization is performed in terms of the parametric coordinates of the nodes to enforce that the nodes slide on the surfaces. Moreover, the proposed algorithm repairs invalid curved meshes (untangling), deals with arbitrary polynomial degrees (high‐order), and handles with low‐quality CAD parameterizations (independence of parameterization). Third, we use the optimization procedure to generate curved nodal high‐order surface meshes by means of an a posteriori approach. Given a linear mesh, we increase the polynomial degree of the elements, curve them to match the geometry, and optimize the location of the nodes to ensure mesh validity. Finally, we present several examples to demonstrate the features of the optimization procedure, and to illustrate the surface mesh generation process. Copyright © 2015 John Wiley & Sons, Ltd.     

Parallel adaptive finite element computations with hierarchical preconditioning

A parallel implementation of an adaptive finite element program is treated which is characterized by an underlying parallel dynamic data structure based on linked lists and tree structures. In conjunction with a conjugate gradient solver an efficient methodology for treating adaptive finite element systems is shown. This is achieved by preconditioning using hierarchical bases with and without a coarse grid solver and by new methods of quasi-optimal load balancing. The different levels of nested meshes needed for preconditioning are governed either by global or by adaptive refinements. A termination algorithm based on the vector method is implemented for the non deterministic adaptive mesh refinement procedure. The problems concerning load balancing due to adaptive refinement are solved by a dynamic load balancing for the nodes.This work has been supported by Deutsche Forschungsgemeinschaft under grant no. Ste238/26.     

Variational generation of prismatic boundary‐layer meshes for biomedical computing

Boundary‐layer meshes are important for numerical simulations in computational fluid dynamics, including computational biofluid dynamics of air flow in lungs and blood flow in hearts. Generating boundary‐layer meshes is challenging for complex biological geometries. In this paper, we propose a novel technique for generating prismatic boundary‐layer meshes for such complex geometries. Our method computes a feature size of the geometry, adapts the surface mesh based on the feature size, and then generates the prismatic layers by propagating the triangulated surface using the face‐offsetting method. We derive a new variational method to optimize the prismatic layers to improve the triangle shapes and edge orthogonality of the prismatic elements and also introduce simple and effective measures to guarantee the validity of the mesh. Coupled with a high‐quality tetrahedral mesh generator for the interior of the domain, our method generates high‐quality hybrid meshes for accurate and efficient numerical simulations. We present comparative study to demonstrate the robustness and quality of our method for complex biomedical geometries. Copyright © 2009 John Wiley & Sons, Ltd.     

Mechanical analysis of 3-D braided composites by the finite multiphase element method

A finite multiphase element method (FMEM), in which the element comprises more than one kind of material, has been proposed to predict the effective elastic properties of 3-D braided composites. This method is based on the variational principle and our previous geometric model that assumes the existence of different types of unit cells in the three regions in a 3-D braided composite, i.e. the interior, surface and corner. The numerical procedure involved two steps. First, a fine local mesh at the unit cell level is used to analyze the stress/strain of each unit cell. Then, a relatively coarse global mesh is used to obtain the overall responses of the composite at macroscopic level. By using the stress volume averaging method, the effective elastic properties of the composite can be calculated under the prescribed uniform strain boundary conditions. Finally, the predicted stress/strain curves are compared with experimental results, demonstrating the applicability of the FME method.     

Three dimensional automatic refinement method for transient small strain elastoplastic finite element computations

In this paper, the refinement strategy based on the “Non-Linear Localized Full MultiGrid” solver originally published in Int. J. Numer. Meth. Engng 84(8):947–971 (2010) for 2-D structural problems is extended to 3-D simulations. In this context, some extra information concerning the refinement strategy and the behavior of the error indicators are given. The adaptive strategy is dedicated to the accurate modeling of elastoplastic materials with isotropic hardening in transient dynamics. A multigrid solver with local mesh refinement is used to reduce the amount of computational work needed to achieve an accurate calculation at each time step. The locally refined grids are automatically constructed, depending on the user prescribed accuracy. The discretization error is estimated by a dedicated error indicator within the multigrid method. In contrast to other adaptive procedures, where grids are erased when new ones are generated, the previous solutions are used recursively to reduce the computing time on the new mesh. Moreover, the adaptive strategy needs no costly coarsening method as the mesh is reassessed at each time step. The multigrid strategy improves the convergence rate of the non-linear solver while ensuring the information transfer between the different meshes. It accounts for the influence of localized non-linearities on the whole structure. All the steps needed to achieve the adaptive strategy are automatically performed within the solver such that the calculation does not depend on user experience. This paper presents three-dimensional results using the adaptive multigrid strategy on elastoplastic structures in transient dynamics and in a linear geometrical framework. Isoparametric cubic elements with energy and plastic work error indicators are used during the calculation.     

Unified fluid–structure interpolation and mesh motion using radial basis functions

A multivariate interpolation scheme, using radial basis functions, is presented, which results in a completely unified formulation for the fluid–structure interpolation and mesh motion problems. The method has several significant advantages. Primarily, all volume mesh, structural mesh, and flow‐solver type dependence is removed, and all operations are performed on totally arbitrary point clouds of any form. Hence, all connectivity and user‐input requirements are removed from the computational fluid dynamics–computational structural dynamics (CFD–CSD) coupling problem, as only point clouds are required to determine the coupling. Also, it may equally well be applied to structured and unstructured grids, or structural and aerodynamic grids that intersect, again because no connectivity information is required. Furthermore, no expensive computations are required during an unsteady simulation, just matrix–vector multiplications, since the required dependence relations are computed only once prior to any simulation and then remain constant. This property means that the method is both perfectly parallel, since only the data relevant to each structured block or unstructured partition are required to move those points, and totally independent from the flow solver. Hence, a completely generic ‘black box’ tool can be developed, which is ideal for use in an optimization approach. Aeroelastic behaviour of the Brite–Euram MDO wing is analysed in terms of both static deflection and dynamic responses, and it is demonstrated that responses are strongly dependent on the exact CFD–CSD interpolation used. Mesh quality is also examined during the motion resulting from a large surface deformation. Global grid quality is shown to be preserved well, with local grid orthogonality also being maintained well, particularly at and near the moving surface, where the original orthogonality is retained. Copyright © 2007 John Wiley & Sons, Ltd.     

PCG solver and its computational complexity for implicit control-volume finite-element method of RTM mold filling simulation

For isoparametric element meshes, the control-volume finite-element method for resin transfer molding (RTM) mold filling generates an asymmetric matrix, and the performance of the pre-conditioner conjugate gradient (PCG) solver decreases by almost one order of magnitude, even for meshes with very few trivial asymmetric data points. In this paper, the asymmetric parts of the linear equations were transferred to the right-hand sides, and then the linear equations were transformed into an equivalent set of symmetric equations. The right-hand sides of the system of equations were updated only when the set of filled nodes changed. The time steps were controlled by the rule of “one time step, one element-size distance.” Based on the PCG solver and the time-step strategy, the computational complexity of the implicit control-volume method was analyzed and presented. Both analytical and case studies showed that the computational complexity of the PCG solver was of order N squared (where N is the number of nodes) for both 2.5D and 3D meshes. The proposed approach was very suitable for a 3D mesh and had the capability of simulating a mesh with 50,000 nodes in under one hour using a 2.0 GHz CPU, 512M RAM computer.     

Novel quadratic Bézier triangular and tetrahedral elements using existing mesh generators: Applications to linear nearly incompressible elastostatics and implicit and explicit elastodynamics

In this paper, we present novel techniques of using quadratic Bézier triangular and tetrahedral elements for elastostatic and implicit/explicit elastodynamic simulations involving nearly incompressible linear elastic materials. A simple linear mapping is proposed for developing finite element meshes with quadratic Bézier triangular/tetrahedral elements from the corresponding quadratic Lagrange elements that can be easily generated using the existing mesh generators. Numerical issues arising in the case of nearly incompressible materials are addressed using the consistent B -bar formulation, thus reducing the finite element formulation to one consisting only of displacements. The higher-order spatial discretization and the nonnegative nature of Bernstein polynomials are shown to yield significant computational benefits. The optimal spatial convergence of the B -bar formulation for the quadratic triangular and tetrahedral elements is demonstrated by computing error norms in displacement and stresses. The applicability and computational efficiency of the proposed elements for elastodynamic simulations are demonstrated by studying several numerical examples involving real-world geometries with complex features. Numerical results obtained with the standard linear triangular and tetrahedral elements are also presented for comparison.     

A coupled meshfree/finite element method for automotive crashworthiness simulations

Mesh distortion induced numerical instability is a major roadblock in automotive crashworthiness finite element simulations. Remedies such as wrapping elements with null shells and deletion of distorted meshes have been adopted but none of them seems robust enough to survive various scenarios. Meshfree methods have been developed over the past almost twenty years in view of their capabilities in dealing with large material deformation and separation, but have remained in academic research due to their unaffordable high computational cost in solving large-scale industrial applications. This paper presents a coupled meshfree/finite-element method which allows engineers to model the severe deformation area with the meshfree method while keeping the remaining area modeled by the finite element methods. The method is implemented into LS-DYNA version 971 and its later versions so that it is available for automotive crashworthiness simulations. In the paper, one linear patch test and three crash examples are presented to demonstrate the accuracy of the meshfree formulation, its effectiveness in resolving mesh distortion difficulty, and the efficiency of the coupled meshfree/finite element solver in handling large-scale models.     

OPTIMIZATION STRATEGIES IN UNSTRUCTURED MESH GENERATION

We propose a new optimization strategy for unstructured meshes that, when coupled with existing automatic generators, produces meshes of high quality for arbitrary domains in 3-D. Our optimizer is based upon a non-differentiable definition of the quality of the mesh which is natural for finite element or finite volume users: the quality of the worst element in the mesh. The dimension of the optimization space is made tractable by restricting, at each iteration, to a suitable neighbourhood of the worst element. Both geometrical (node repositioning) and topological (reconnection) operations are performed. It turns out that the repositioning method is advantageous with respect to both the usual node-by-node techniques and the more recent differentiable optimization methods. Several examples are included that illustrate the efficiency of the optimizer.     

An approach to refining three-dimensional tetrahedral meshes based on Delaunay transformations

A technique for refining three-dimensional tetrahedral meshes is proposed in this paper. The proposed technique is capable of treating arbitrary unstructured tetrahedral meshes, convex or non-convex with multiple regions resulting in high quality constrained Delaunay triangulations. The tetrahedra generated are of high quality (nearly equilateral). Sliver tetrahedra, which present a real problem to many algorithms are not produced with the new method. The key to the generation of high quality tetrahedra is the iterative application of a set of topological transformations based on the Voronoi–Delaunay theory and a reposition of nodes technique. The computational requirements of the proposed technique are in linear relationship with the number of nodes and tetrahedra, making it ideal for direct employment in a fully automatic finite element analysis system for 3-D adaptive mesh refinement. Application to some test problems is presented to show the effectiveness and applicability of the new method.     

An automatic adaptive refinement finite element procedure for 2D elastostatic analysis

An automatic adaptive refinement procedure for finite element analysis is presented. The procedure is applied to two-dimensional elastostatic problems to obtain solutions of prescribed accuracy. Through the combined use of new mesh generator using contour developed by Lo¹ and the concept of strain energy concentration, high-quality graded finite element meshes are generated. The whole process is fully automatic and no user intervention is required during the successive cycles of the mesh refinements. The Zienkiewicz and Zhu² error estimator is found to be effective and has been adopted for the present implementation. In the numerical examples tested, the error estimator gives an accurate error norm estimation and the effectivity index of the estimator converges to a value close to unity.