Automatic generation of transitional meshes

Advances in commercial computer‐aided design software have made finite element analysis with three‐dimensional solid finite elements routinely available. Since these analyses usually confine themselves to those geometrical objects for which particular CAD systems can produce finite element meshes, expanding the capability of analyses becomes an issue of expanding the capability of generating meshes. This paper presents a method for stitching together two three‐dimensional meshes with diverse elements that can include tetrahedral, pentahedral and hexahedral solid finite elements. The stitching produces a mesh that coincides with the edges which already exist on the portion of boundaries that will be joined. Moreover, the transitional mesh does not introduce new edges on these boundaries. Since the boundaries of the regions to be stitched together can have a mixture of triangles and quadrilaterals, tetrahedral and pyramidal elements provide the transitional elements required to honor these constraints. On these boundaries a pyramidal element shares its base face with the quadrilateral faces of hexahedra and pentahedra. Tetrahedral elements share a face with the triangles on the boundary. Tetrahedra populate the remaining interior of the transitional region. Copyright © 2001 John Wiley & Sons, Ltd.     

A face‐based smoothed finite element method (FS‐FEM) for 3D linear and geometrically non‐linear solid mechanics problems using 4‐node tetrahedral elements

This paper presents a novel face‐based smoothed finite element method (FS‐FEM) to improve the accuracy of the finite element method (FEM) for three‐dimensional (3D) problems. The FS‐FEM uses 4‐node tetrahedral elements that can be generated automatically for complicated domains. In the FS‐FEM, the system stiffness matrix is computed using strains smoothed over the smoothing domains associated with the faces of the tetrahedral elements. The results demonstrated that the FS‐FEM is significantly more accurate than the FEM using tetrahedral elements for both linear and geometrically non‐linear solid mechanics problems. In addition, a novel domain‐based selective scheme is proposed leading to a combined FS/NS‐FEM model that is immune from volumetric locking and hence works well for nearly incompressible materials. The implementation of the FS‐FEM is straightforward and no penalty parameters or additional degrees of freedom are used. The computational efficiency of the FS‐FEM is found better than that of the FEM. Copyright © 2008 John Wiley & Sons, Ltd.     

Linear finite elements with orthogonal discontinuous basis functions for explicit earthquake ground motion modeling

The dynamic explicit finite element method is commonly used in earthquake ground motion modeling. In this method, the element mass matrix is approximately lumped, which may lead to numerical dispersion. On the other hand, the orthogonal finite element method, based on orthogonal polynomial basis functions, naturally derives a lumped diagonal mass matrix and can be applied to dynamic explicit finite element analysis. In this paper, we propose finite elements based on orthogonal discontinuous basis functions, the element mass matrices of which are lumped without approximation. Orthogonal discontinuous basis functions are used to improve the accuracy and reduce the numerical dispersion in earthquake ground motion modeling. We present a detailed formulation of the 4‐node tetrahedral and 8‐node hexahedral elements. The relationship between the proposed finite elements and conventional finite elements is investigated, and the solutions obtained from the conventional explicit finite element method are compared with analytical solutions to verify the numerical dispersion caused by the lumping approximation. Comparison of solutions obtained with the proposed finite elements to analytical solutions demonstrates the usefulness of the technique. Examples are also presented to illustrate the effectiveness of the proposed method in earthquake ground motion modeling in the actual three‐dimensional crust structure. Copyright © 2010 John Wiley & Sons, Ltd.     

Achieving finite element mesh quality via optimization of the Jacobian matrix norm and associated quantities. Part II—A framework for volume mesh optimization and the condition number of the Jacobian matrix

Three‐dimensional unstructured tetrahedral and hexahedral finite element mesh optimization is studied from a theoretical perspective and by computer experiments to determine what objective functions are most effective in attaining valid, high‐quality meshes. The approach uses matrices and matrix norms to extend the work in Part I to build suitable 3D objective functions. Because certain matrix norm identities which hold for 2×2 matrices do not hold for 3×3 matrices, significant differences arise between surface and volume mesh optimization objective functions. It is shown, for example, that the equality in two dimensions of the smoothness and condition number of the Jacobian matrix objective functions does not extend to three dimensions and further, that the equality of the Oddy and condition number of the metric tensor objective functions in two dimensions also fails to extend to three dimensions. Matrix norm identities are used to systematically construct dimensionally homogeneous groups of objective functions. The concept of an ideal minimizing matrix is introduced for both hexahedral and tetrahedral elements. Non‐dimensional objective functions having barriers are emphasized as the most logical choice for mesh optimization. The performance of a number of objective functions in improving mesh quality was assessed on a suite of realistic test problems, focusing particularly on all‐hexahedral ‘whisker‐weaved’ meshes. Performance is investigated on both structured and unstructured meshes and on both hexahedral and tetrahedral meshes. Although several objective functions are competitive, the condition number objective function is particularly attractive. The objective functions are closely related to mesh quality measures. To illustrate, it is shown that the condition number metric can be viewed as a new tetrahedral element quality measure. Published in 2000 by John Wiley & Sons, Ltd.     

Converting a tetrahedral mesh to a prism–tetrahedral hybrid mesh for FEM accuracy and efficiency

This paper presents a computational method for converting a tetrahedral mesh to a prism–tetrahedral hybrid mesh for improved solution accuracy and computational efficiency of finite element analysis. The proposed method performs this conversion by inserting layers of prism elements and deleting tetrahedral elements in sweepable sub‐domains, in which cross‐sections remain topologically identical and geometrically similar along a certain sweeping path. The total number of finite elements is reduced because roughly three tetrahedral elements are converted to one prism element. The solution accuracy of the finite element analysis improves since a prism element yields a more accurate solution than a tetrahedral element due to the presence of higher‐order terms in the shape function. Only previously known method for creating such a prism–tetrahedral hybrid mesh was to manually decompose a target volume into sweepable and non‐sweepable sub‐volumes and mesh each of the sub‐volumes separately. Unlike the previous method, the proposed method starts from a cross‐section of a tetrahedral mesh and replaces the tetrahedral elements with layers of prism elements until prescribed quality criteria can no longer be satisfied. A series of computational fluid dynamics simulations and structural analyses have been conducted, and the results verified a better performance of prism–tetrahedral hybrid mesh. Copyright © 2009 John Wiley & Sons, Ltd.     

A displacement‐based finite element formulation for general polyhedra using harmonic shape functions

A displacement‐based continuous‐Galerkin finite element formulation for general polyhedra is presented for applications in nonlinear solid mechanics. The polyhedra can have an arbitrary number of vertices or faces. The faces of the polyhedra can have an arbitrary number of edges and can be nonplanar. The polyhedra can be nonconvex with only the mild restriction of star convexity with respect to the vertex‐averaged centroid. Conforming shape functions are constructed using harmonic functions defined on the undeformed configuration, thus requiring the use of a total‐Lagrangian finite element formulation in large deformation applications. For nonlinear applications with computationally intensive constitutive models, it is important to minimize the number of element quadrature points. For this reason, an integration scheme is adopted in which the number of quadrature points is equal to the number of vertices. As a first step toward verifying the element behavior in the general nonlinear setting, several linear verification examples are presented using both random Voronoi meshes and distorted hexahedral meshes. For the hexahedral meshes, results for the polyhedral formulation are compared to those of the standard trilinear hexahedral formulation. The element behavior in the nearly incompressible regime is also examined. Copyright © 2013 John Wiley & Sons, Ltd.     

A practical rule for the derivation of transition finite elements

The derivation of shape functions of two‐dimensional regular and transition finite element using Pascal triangle concept is presented. A practical rule for the trial function selection from Pascal triangle is developed. The derived elements are tested using some case studies. The results are validated by ANSYS. The developed rule is also generalized for the three‐dimensional finite elements. Copyright © 2000 John Wiley & Sons, Ltd.     

Polyhedral elements by means of node/edge‐based smoothed finite element method

The node‐based or edge‐based smoothed finite element method is extended to develop polyhedral elements that are allowed to have an arbitrary number of nodes or faces, and so retain a good geometric adaptability. The strain smoothing technique and implicit shape functions based on the linear point interpolation make the element formulation simple and straightforward. The resulting polyhedral elements are free from the excessive zero‐energy modes and yield a robust solution very much insensitive to mesh distortion. Several numerical examples within the framework of linear elasticity demonstrate the accuracy and convergence behavior. The smoothed finite element method‐based polyhedral elements in general yield solutions of better accuracy and faster convergence rate than those of the conventional finite element methods. Copyright © 2016 John Wiley & Sons, Ltd.     

A transition element for uniform strain hexahedral and tetrahedral finite elements

A transition element is presented for meshes containing uniform strain hexahedral and tetrahedral finite elements. It is shown that the volume of the standard uniform strain hexahedron is identical to that of a polyhedron with 14 vertices and 24 triangular faces. Based on this equivalence, a transition element is developed as a simple modification of the uniform strain hexahedron. The transition element makes use of a general method for hourglass control and satisfies first‐order patch tests. Example problems in linear elasticity are included to demonstrate the application of the element. Copyright © 1999 John Wiley & Sons, Ltd. This paper was produced under the auspices of the U.S. Government and it is therefore not subject to copyright in the U.S.     

A spline wavelet finite‐element method in structural mechanics

The wavelet‐based methods are powerful to analyse the field problems with changes in gradients and singularities due to the excellent multi‐resolution properties of wavelet functions. Wavelet‐based finite elements are often constructed in the wavelet space where field displacements are expressed as a product of wavelet functions and wavelet coefficients. When a complex structural problem is analysed, the interface between different elements and boundary conditions cannot be easily treated as in the case of conventional finite‐element methods (FEMs). A new wavelet‐based FEM in structural mechanics is proposed in the paper by using the spline wavelets, in which the formulation is developed in a similar way of conventional displacement‐based FEM. The spline wavelet functions are used as the element displacement interpolation functions and the shape functions are expressed by wavelets. The detailed formulations of typical spline wavelet elements such as plane beam element, in‐plane triangular element, in‐plane rectangular element, tetrahedral solid element, and hexahedral solid element are derived. The numerical examples have illustrated that the proposed spline wavelet finite‐element formulation achieves a high numerical accuracy and fast convergence rate. Copyright © 2005 John Wiley & Sons, Ltd.     

A general 3D BEM/FEM coupling applied to elastodynamic continua/frame structures interaction analysis

This paper presents a procedure for coupling general finite element models with three‐dimensional bodies modelled by the Boundary Element Method (BEM). Shells, plates and frames are modelled by the Finite Element Method (FEM) and coupled to the BEM domain directly or by means of rigid blocks. The coupling is used for the analysis of buildings connected to half‐space by means of rigid footings, piles or plates in bending and other problems where combinations of different types of sub‐domains are required, composite domains for instance. Several numerical examples are analysed to demonstrate the robustness and accuracy of the proposed scheme. Copyright © 1999 John Wiley & Sons, Ltd.     

A finite element method based on C0‐continuous assumed gradients

This article presents an alternative approach to assumed gradient methods in FEM applied to three‐dimensional elasticity. Starting from nodal integration (NI), a general C0‐continuous assumed interpolation of the deformation gradient is formulated. The assumed gradient is incorporated using the principle of Hu‐Washizu. By dual Lagrange multiplier spaces, the functional is reduced to the displacements as the only unknowns. An integration scheme is proposed where the integration points coincide with the support points of the interpolation. Requirements for regular finite element meshes are explained. Using this interpretation of NI, instabilities (appearance of spurious modes) can be explained. The article discusses and classifies available strategies to stabilize NI such as penalty methods, SCNI, α‐FEM. Related approaches, such as the smoothed finite element method, are presented and discussed. New stabilization techniques for NI are presented being entirely based on the choice of the assumed gradient interpolation, i.e. nodal‐bubble support, edge‐based support and support using tensor‐product interpolations. A strategy is presented on how the interpolation functions can be derived for various element types. Interpolation functions for the first‐order hexahedral element, the first‐order and the second‐order tetrahedral elements are given. Numerous examples illustrate the strengths and limitations of the new schemes. Copyright © 2010 John Wiley & Sons, Ltd.     

A partition‐of‐unity‐based finite element method for level sets

Level set methods have recently gained much popularity to capture discontinuities, including their possible propagation. Typically, the partial differential equations that arise in level set methods, in particular the Hamilton–Jacobi equation, are solved by finite difference methods. However, finite difference methods are less suited for irregular domains. Moreover, it seems slightly awkward to use finite differences for the capturing of a discontinuity, while in a subsequent stress analysis finite elements are normally used. For this reason, we here present a finite element approach to solving the governing equations of level set methods. After a review of the governing equations, the initialization of the level sets, the discretization on a finite domain, and the stabilization of the resulting finite element method will be discussed. Special attention will be given to the proper treatment of the internal boundary condition, which is achieved by exploiting the partition‐of‐unity property of finite element shape functions. Finally, a quantitative analysis including accuracy analysis is given for a one‐dimensional example and a qualitative example is given for a two‐dimensional case with a curved discontinuity. Copyright © 2008 John Wiley & Sons, Ltd.     

A unified 3D‐based technique for plate bending analysis using scaled boundary finite element method

This paper presents a unified technique for solving the plate bending problems by extending the scaled boundary finite element method. The formulation is based on the three‐dimensional governing equation without enforcing the kinematics of plate theory. Only the in‐plane dimensions are discretised into finite elements. Any two‐dimensional displacement‐based elements can be employed. The solution along the thickness is expressed analytically by using a matrix function. The proposed technique is consistent with the three‐dimensional theory and applicable to both thick and thin plates without exhibiting the numerical locking phenomenon. Moreover, the use of higher order spectral elements allows the proposed technique to better represent curved boundaries and to achieve high accuracy and fast convergence. Numerical examples of various plate structures with different thickness‐to‐length ratios demonstrate the applicability and accuracy of the proposed technique. Copyright © 2012 John Wiley & Sons, Ltd.     

Unstructured,curved elements for the two‐dimensional high order discontinuous control‐volume/finite‐element method

Quadrilateral and triangular elements with curved edges are developed in the framework of spectral, discontinuous, hybrid control‐volume/finite‐element method for elliptic problems. In order to accommodate hybrid meshes, encompassing both triangular and quadrilateral elements, one single mapping is used. The scheme is applied to two‐dimensional problems with discontinuous, anisotropic diffusion coefficients, and the exponential convergence of the method is verified in the presence of curved geometries. Copyright © 2014 John Wiley & Sons, Ltd.     

A suitable low‐order,tetrahedral finite element for solids

To use the all‐tetrahedral mesh generation capabilities existing today, we have explored the creation of a computationally efficient eight‐node tetrahedral finite element (a four‐node tetrahedral finite element enriched with four mid‐face nodal points). The derivation of the element's gradient operator, studies in obtaining a suitable mass lumping and the element's performance in applications are presented. In particular, we examine the eight‐node tetrahedral finite element's behavior in longitudinal plane wave propagation, in transverse cylindrical wave propagation, and in simulating Taylor bar impacts. The element samples only constant strain states and, therefore, has 12 hourglass modes. In this regard, it bears similarities to the eight‐node, mean‐quadrature hexahedral finite element. Comparisons with the results obtained from the mean‐quadrature eight‐node hexahedral finite element and the four‐node tetrahedral finite element are included. Given automatic all‐tetrahedral meshing, the eight‐node, mean‐quadrature tetrahedral finite element is a suitable replacement for the eight‐node, mean‐quadrature hexahedral finite element and meshes requiring an inordinate amount of user intervention and direction to generate. Copyright © 1999 John Wiley & Sons, Ltd. This paper was produced under the auspices of the U.S. Government and it is therefore not subject to copyright in the U.S.     

Clustered generalized finite element methods for mesh unrefinement,non‐matching and invalid meshes

In spite of significant advancements in automatic mesh generation during the past decade, the construction of quality finite element discretizations on complex three‐dimensional domains is still a difficult and time demanding task. In this paper, the partition of unity framework used in the generalized finite element method (GFEM) is exploited to create a very robust and flexible method capable of using meshes that are unacceptable for the finite element method, while retaining its accuracy and computational efficiency. This is accomplished not by changing the mesh but instead by clustering groups of nodes and elements. The clusters define a modified finite element partition of unity that is constant over part of the clusters. This so‐called clustered partition of unity is then enriched to the desired order using the framework of the GFEM. The proposed generalized finite element method can correctly and efficiently deal with: (i) elements with negative Jacobian; (ii) excessively fine meshes created by automatic mesh generators; (iii) meshes consisting of several sub‐domains with non‐matching interfaces. Under such relaxed requirements for an acceptable mesh, and for correctly defined geometries, today's automated tetrahedral mesh generators can practically guarantee successful volume meshing that can be entirely hidden from the user. A detailed technical discussion of the proposed generalized finite element method with clustering along with numerical experiments and some implementation details are presented. Copyright © 2006 John Wiley & Sons, Ltd.     

Three‐dimensional finite elements with embedded strong discontinuities to model material failure in the infinitesimal range

This paper presents new three‐dimensional finite elements with embedded strong discontinuities in the small strain infinitesimal range. The goal is to model localized surfaces of failure in solids, such as cracks at fracture, through enhancements of the finite elements that capture the propagating discontinuities of the displacement field in the element interiors. In this way, such surfaces of discontinuity can be sharply resolved in general meshes not necessarily related to the detailed geometry of the surface, unknown a priori. An important issue is also the consideration of general finite element formulations in the developments (e.g., basic displacement‐based, mixed or enhanced assumed strain finite element formulations), as needed to optimally resolve the continuum problem in the bulk. The actual modeling of the discontinuity effects, including the incorporation of the cohesive law defining the discontinuity constitutive response, is carried out at the element level with the proper enhancement of the discrete strain field of the element. The added elemental degrees of freedom approximate the displacement jumps associated with the discontinuity and are defined independently from element to element, thus allowing their static condensation at the element level without affecting the global mechanical problem in terms of the number and topology of the global degrees of freedom. In fact, this global‐local structure of the finite element methods developed in this work arises naturally from a multi‐scale characterization of these localized solutions, with the discontinuities understood to appear in the small scales, thus leading directly to these computationally efficient numerical methods for their numerical resolution, easily incorporated to an existing finite element code. The focus in this paper is on the development of finite elements incorporating a linear interpolation of the displacement jumps in the general three‐dimensional setting. These interpolations are shown to be necessary for hexahedral elements to avoid the so‐called stress locking that occurs with simpler constant approximations of the jumps (namely, a spurious transfer of stresses across the discontinuity not allowing its full release and, hence, resulting in an overstiff or locked numerical solution). The design of the new finite elements is accomplished in this work by a direct identification of the separation modes to be incorporated in the discrete strain field of the element, rather than from an assumed discontinuous interpolation of the displacements, assuring with this approach their locking‐free response by design. An additional issue addressed in the paper is the geometric characterization and propagation of the discontinuity surfaces in the general three‐dimensional setting of interest here. The paper includes a series of numerical simulations illustrating and evaluating the properties of the new finite elements. Copyright © 2012 John Wiley & Sons, Ltd.     

A stabilized nodally integrated tetrahedral

A stabilized, nodally integrated linear tetrahedral is formulated and analysed. It is well known that linear tetrahedral elements perform poorly in problems with plasticity, nearly incompressible materials, and acute bending. For a variety of reasons, low‐order tetrahedral elements are preferable to quadratic tetrahedral elements; particularly for nonlinear problems. But the severe locking problems of tetrahedrals have forced analysts to employ hexahedral formulations for most nonlinear problems. On the other hand, automatic mesh generation is often not feasible for building many 3D hexahedral meshes. A stabilized, nodally integrated linear tetrahedral is developed and shown to perform very well in problems with plasticity, nearly incompressible materials and acute bending. The formulation is analytically and numerically shown to be stable and optimally convergent for the compressible case provided sufficient smoothness of the exact solution u ∈ C² ∩ (H¹)³. Future work may extend the formulation to the incompressible regime and relax the regularity requirements; nonetheless, the results demonstrate that the method is not susceptible to locking and performs quite well in several standard linear and nonlinear benchmarks. Published in 2006 by John Wiley & Sons, Ltd.     

Boundary element formulation for 3D transversely isotropic cracked bodies

The boundary traction integral representation is obtained in elasticity when the classical displacement representation is differentiated and combined according to Hooke's law. The use of both traction and displacement integral representations leads to a mixed (or dual) formulation of the BEM where the discretization effort for crack problems is much smaller than in the classical formulation. A boundary element analysis of three‐dimensional fracture mechanics problems of transversely isotropic solids based on the mixed formulation is presented in this paper. The hypersingular and strongly singular kernels appearing in the formulation are regularized by using two terms of the displacement series expansion and one term of the traction expansion, at the collocation point. All the remaining integrals are analytically evaluated or transformed by means of Stokes' theorem into regular or weakly singular integrals, which are numerically computed. The method is general and can be used for elements of any shape including quarter‐point crack front elements. No change of co‐ordinates is required for the integration. The formulation as presented in this paper is something as clear, general and easy to handle as the classical BE formulation. It is used in combination with three‐dimensional quadratic and quarter‐point elements to obtain accurate results for several different crack problems. Cracks in boundless and finite transversely isotropic domains are studied. The meshes are simple and include only discretization of the crack and the external boundary. The obtained results are in good agreement with those existing in the literature. Copyright © 2004 John Wiley & Sons, Ltd.