An unsymmetric 8‐node hexahedral element with high distortion tolerance

Among all 3D 8‐node hexahedral solid elements in current finite element library, the ‘best’ one can produce good results for bending problems using coarse regular meshes. However, once the mesh is distorted, the accuracy will drop dramatically. And how to solve this problem is still a challenge that remains outstanding. This paper develops an 8‐node, 24‐DOF (three conventional DOFs per node) hexahedral element based on the virtual work principle, in which two different sets of displacement fields are employed simultaneously to formulate an unsymmetric element stiffness matrix. The first set simply utilizes the formulations of the traditional 8‐node trilinear isoparametric element, while the second set mainly employs the analytical trial functions in terms of 3D oblique coordinates (R, S, T). The resulting element, denoted by US‐ATFH8, contains no adjustable factor and can be used for both isotropic and anisotropic cases. Numerical examples show it can strictly pass both the first‐order (constant stress/strain) patch test and the second‐order patch test for pure bending, remove the volume locking, and provide the invariance for coordinate rotation. Especially, it is insensitive to various severe mesh distortions. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐performance geometric nonlinear analysis with the unsymmetric 4‐node, 8‐DOF plane element US‐ATFQ4

A recent unsymmetric 4‐node, 8‐DOF plane element US‐ATFQ4, which exhibits excellent precision and distortion‐resistance for linear elastic problems, is extended to geometric nonlinear analysis. Since the original linear element US‐ATFQ4 contains the analytical solutions for plane pure bending, how to modify such formulae into incremental forms for nonlinear applications and design an appropriate updated algorithm become the key of the whole job. First, the analytical trial functions should be updated at each iterative step in the framework of updated Lagrangian formulation that takes the configuration at the beginning of an incremental step as the reference configuration during that step. Second, an appropriate stress update algorithm in which the Cauchy stresses are updated by the Hughes‐Winget method is adopted to estimate current stress fields. Numerical examples show that the new nonlinear element US‐ATFQ4 also possesses amazing performance for geometric nonlinear analysis, no matter whether regular or distorted meshes are used. It again demonstrates the advantages of the unsymmetric finite element method with analytical trial functions.     

An unsymmetric 4‐node, 8‐DOF plane membrane element perfectly breaking through MacNeal's theorem

Among numerous finite element techniques, few models can perfectly (without any numerical problems) break through MacNeal's theorem: any 4‐node, 8‐DOF membrane element will either lock in in‐plane bending or fail to pass a C₀ patch test when the element's shape is an isosceles trapezoid. In this paper, a 4‐node plane quadrilateral membrane element is developed following the unsymmetric formulation concept, which means two different sets of interpolation functions for displacement fields are simultaneously used. The first set employs the shape functions of the traditional 4‐node bilinear isoparametric element, while the second set adopts a novel composite coordinate interpolation scheme with analytical trail function method, in which the Cartesian coordinates (x,y) and the second form of quadrilateral area coordinates (QACM‐II) (S,T) are applied together. The resulting element US‐ATFQ4 exhibits amazing performance in rigorous numerical tests. It is insensitive to various serious mesh distortions, free of trapezoidal locking, and can satisfy both the classical first‐order patch test and the second‐order patch test for pure bending. Furthermore, because of usage of the second form of quadrilateral area coordinates (QACM‐II), the new element provides the invariance for the coordinate rotation. It seems that the behaviors of the present model are beyond the well‐known contradiction defined by MacNeal's theorem. Copyright © 2015 John Wiley & Sons, Ltd.     

A novel unsymmetric 8‐node plane element immune to mesh distortion under a quadratic displacement field

An 8‐node quadrilateral plane finite element is developed based on a novel unsymmetric formulation which is characterized by the use of two sets of shape functions, viz., the compatibility enforcing shape functions and completeness enforcing shape functions. The former are chosen to satisfy exactly the minimum inter‐ as well as intra‐element displacement continuity requirements, while the latter are chosen to satisfy all the (linear and higher order) completeness requirements so as to reproduce exactly a quadratic displacement field. Numerical results from test problems reveal that the new element is indeed capable of reproducing exactly a complete quadratic displacement field under all types of admissible mesh distortions. In this respect, the proposed 8‐node unsymmetric element emerges to be better than the existing symmetric QUAD8, QUAD8/9, QUAD9, QUAD12 and QUAD16 elements, and matches the performance of the quartic element, QUAD25. For test problems involving a cubic or higher order displacement field, the proposed element yields a solution accuracy that is comparable to or better than that of QUAD8, QUAD8/9 and QUAD9 elements. Furthermore, the element maintains a good accuracy even with the reduced 2× 2 numerical integration. Copyright © 2003 John Wiley & Sons, Ltd.     

4‐node unsymmetric quadrilateral membrane element with drilling DOFs insensitive to severe mesh‐distortion

The unsymmetric finite element method is a promising technique to produce distortion‐immune finite elements. In this work, a simple but robust 4‐node 12‐DOF unsymmetric quadrilateral membrane element is formulated. The test function of this new element is determined by a concise isoparametric‐based displacement field that is enriched by the Allman‐type drilling degrees of freedom. Meanwhile, a rational stress field, instead of the displacement one in the original unsymmetric formulation, is directly adopted to be the element's trial function. This stress field is obtained based on the analytical solutions of the plane stress/strain problem and the quasi‐conforming technique. Thus, it can a priori satisfy related governing equations. Numerical tests show that the presented new unsymmetric element, named as US‐Q4θ, exhibits excellent capabilities in predicting results of both displacement and stress, in most cases, superior to other existing 4‐node element models. In particular, it can still work very well in severely distorted meshes even when the element shape deteriorates into concave quadrangle or degenerated triangle.     

Quadrilateral membrane elements with analytical element stiffness matrices formulated by the new quadrilateral area coordinate method (QACM‐II)

A novel strategy for developing low‐order membrane elements with analytical element stiffness matrices is proposed. First, some complete low‐order basic analytical solutions for plane stress problems are given in terms of the new quadrilateral area coordinates method (QACM‐II). Then, these solutions are taken as the trial functions for developing new membrane elements. Thus, the interpolation formulae for displacement fields naturally possess second‐order completeness in physical space (Cartesian coordinates). Finally, by introducing nodal conforming conditions, new 4‐node and 5‐node membrane elements with analytical element stiffness matrices are successfully constructed. The resulting models, denoted as QAC‐ATF4 and QAC‐ATF5, have high computational efficiency since the element stiffness matrices are formulated explicitly and no internal parameter is added. These two elements exhibit excellent performance in various bending problems with mesh distortion. It is demonstrated that the proposed strategy possesses advantages of both the analytical and the discrete method, and the QACM‐II is a powerful tool for constructing high‐performance quadrilateral finite element models. Copyright © 2008 John Wiley & Sons, Ltd.     

An unsymmetric 8-node hexahedral solid-shell element with high distortion tolerance: Geometric nonlinear formulations

A recent distortion-tolerant unsymmetric 8-node hexahedral solid-shell element US-ATFHS8, which takes the analytical solutions of linear elasticity as the trial functions, is successfully extended to geometric nonlinear analysis. This extension is based on the corotational (CR) approach due to its simplicity and high efficiency, especially for geometric nonlinear analysis where the strain is still small. Based on the assumption that the analytical trial functions can properly work in each increment during the nonlinear analysis, the incremental corotational formulations of the nonlinear solid-shell element US-ATFHS8 are derived within the updated Lagrangian (UL) framework, in which an appropriate updated strategy for linear analytical trial functions is proposed. Numerical examples show that the present nonlinear element US-ATFHS8 possesses excellent performance for various rigorous tests no matter whether regular or distorted mesh is used. Especially, it even performs well in some situations that other conventional elements cannot work.     

Hybrid displacement function element method: a simple hybrid‐Trefftz stress element method for analysis of Mindlin–Reissner plate

In order to develop robust finite element models for analysis of thin and moderately thick plates, a simple hybrid displacement function element method is presented. First, the variational functional of complementary energy for Mindlin–Reissner plates is modified to be expressed by a displacement function F, which can be used to derive displacement components satisfying all governing equations. Second, the assumed element resultant force fields, which can satisfy all related governing equations, are derived from the fundamental analytical solutions of F. Third, the displacements and shear strains along each element boundary are determined by the locking‐free formulae based on the Timoshenko's beam theory. Finally, by applying the principle of minimum complementary energy, the element stiffness matrix related to the conventional nodal displacement DOFs is obtained. Because the trial functions of the domain stress approximations a priori satisfy governing equations, this method is consistent with the hybrid‐Trefftz stress element method. As an example, a 4‐node, 12‐DOF quadrilateral plate bending element, HDF‐P4‐11 β, is formulated. Numerical benchmark examples have proved that the new model possesses excellent precision. It is also a shape‐free element that performs very well even when a severely distorted mesh containing concave quadrilateral and degenerated triangular elements is employed. Copyright © 2014 John Wiley & Sons, Ltd.     

Hyperelastic finite deformation analysis with the unsymmetric finite element method containing homogeneous solutions of linear elasticity

A recent unsymmetric 4-node, 8-DOF plane finite element US-ATFQ4 is generalized to hyperelastic finite deformation analysis. Since the trial functions of US-ATFQ4 contain the homogenous closed analytical solutions of governing equations for linear elasticity, the key of the proposed strategy is how to deal with these linear analytical trial functions (ATFs) during the hyperelastic finite deformation analysis. Assuming that the ATFs can properly work in each increment, an algorithm for updating the deformation gradient interpolated by ATFs is designed. Furthermore, the update of the corresponding ATFs referred to current configuration is discussed with regard to the hyperelastic material model, and a specified model, neo-Hookean model, is employed to verify the present formulation of US-ATFQ4 for hyperelastic finite deformation analysis. Various examples show that the present formulation not only remain the high accuracy and mesh distortion tolerance in the geometrically nonlinear problems, but also possess excellent performance in the compressible or quasi-incompressible hyperelastic finite deformation problems where the strain is large.     

An incompatible and unsymmetric four-node quadrilateral plane element with high numerical performance

Recent studies show that the unsymmetric finite element method exhibits excellent performance when the discretized meshes are severely distorted. In this article, a new unsymmetric 4-noded quadrilateral plane element is presented using both incompatible test functions and trial functions. Five internal nodes, one at the elemental central and four at the middle sides, are added to ensure the quadratic completeness of the elemental displacement field. Thereafter, the total nine nodes are applied to form the shape functions of trial function, and the Lagrange interpolation functions are adopted as the incompatible test shape functions of the internal nodes. The incompatible test displacements are then revised to satisfy the patch test. Numerical tests show that the present element can provide very good numerical accuracy with badly distorted meshes. Unlike the existing unsymmetric four-node plane elements in which the analytical stress fields are employed, the present element can be extended to boundary value problems of any differential equations with no difficulties.     

Analysis of a rotation‐free 4‐node shell element

In this paper, a shell element for small and large deformations is presented based on the extension of the methodology to derive triangular shell element without rotational degrees of freedom (so‐called rotation‐free). As in our original triangular S3 element, the curvatures are computed resorting to the surrounding elements. However, the extension to a quadrilateral element requires internal curvatures in order to avoid singular bending stiffness. The quadrilateral area co‐ordinates interpolation is used to establish the required expressions between the rigid‐body modes of normal nodal translations and the normal through thickness bending strains at mid‐side. In order to propose an attractive low‐cost shell element, the one‐point quadrature is achieved at the centre for the membrane strains, which are superposed to the bending strains in the centred co‐rotational local frame. The membrane hourglass control is obtained by the perturbation stabilization procedure. Free, simply supported and clamped edges are considered without introducing virtual nodes or elements. Several numerical examples with regular and irregular meshes are performed to show the convergence, accuracy and the reasonable little sensitivity to geometric distortion. Based on an updated Lagrangian formulation and Newton iterations, the large displacements of the pinched hemispherical shell show the effectiveness of the proposed simplified element (S4). Finally, the deep drawing of a square box including large plastic strains with contact and friction completes the ability of the rotation‐free quadrilateral element for sheet‐metal‐forming simulations. Copyright © 2005 John Wiley & Sons, Ltd.     

A 20‐node hexahedron element with enhanced distortion tolerance

This paper presents a distortion resistant 20‐node hexahedron element that employs two different sets of shape functions for the trial and test functions. The formulation seeks to satisfy the continuity and completeness requirements by exploiting the intrinsic properties of these two sets of shape functions. Several test problems are used to assess the performance of the element under various mesh distortions. The ability of the proposed as well as the classical 20‐node element to maintain solution accuracy under severe mesh distortions has been studied. The proposed element exhibits a very high tolerance to mesh distortions. In particular, for problems involving linear and quadratic displacement fields, the element is capable of reproducing exact solution under all admissible geometrical distortions of the mesh. For test problems involving higher‐order displacement fields, the performance of the present element is in general better than that of the classical element. Copyright © 2004 John Wiley & Sons, Ltd.     

A hybrid‐mixed four‐node quadrilateral plate element based on sampling surfaces method for 3D stress analysis

The hybrid‐mixed assumed natural strain four‐node quadrilateral element using the sampling surfaces (SaS) technique is developed. The SaS formulation is based on choosing inside the plate body N not equally spaced SaS parallel to the middle surface in order to introduce the displacements of these surfaces as basic plate variables. Such choice of unknowns with the consequent use of Lagrange polynomials of degree N–1 in the thickness direction permits the presentation of the plate formulation in a very compact form. The SaS are located at Chebyshev polynomial nodes that allow one to minimize uniformly the error due to the Lagrange interpolation. To avoid shear locking and have no spurious zero energy modes, the assumed natural strain concept is employed. The developed hybrid‐mixed four‐node quadrilateral plate element passes patch tests and exhibits a superior performance in the case of coarse distorted mesh configurations. It can be useful for the 3D stress analysis of thin and thick plates because the SaS formulation gives the possibility to obtain solutions with a prescribed accuracy, which asymptotically approach the 3D exact solutions of elasticity as the number of SaS tends to infinity. Copyright © 2016 John Wiley & Sons, Ltd.     

A four‐node,shear‐deformable shell element developed via explicit Kirchhoff constraints

An efficient, four‐node quadrilateral shell element is formulated using a linear, first‐order shear deformation theory. The bending part of the formulation is constructed from a cross‐diagonal assembly of four three‐node anisoparametric triangular plate elements, referred to as MIN3. Closed‐form constraint equations, which arise from the Kirchhoff constraints in the thin‐plate limit, are derived and used to eliminate the degrees‐of‐freedom associated with the ‘internal’ node of the cross‐diagonal assembly. The membrane displacement field employs an Allman‐type, drilling degrees‐of‐freedom formulation. The result is a displacement‐based, fully integrated, four‐node quadrilateral element, MIN4T, possessing six degrees‐of‐freedom at each node. Results for a set of validation plate problems demonstrate that the four‐node MIN4T has similar robustness and accuracy characteristics as the original cross‐diagonal assembly of MIN3 elements involving five nodes. The element performs well in both moderately thick and thin regimes, and it is free of shear locking. Shell validation results demonstrate superior performance of MIN4T over MIN3, possibly as a result of its higher‐order interpolation of the membrane displacements. It is also noted that the bending formulation of MIN4T is kinematically compatible with the existing anisoparametric elements of the same order of approximation, which include a two‐node Timoshenko beam element and a three‐node plate element, MIN3. 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.     

Parallel computing of high‐speed compressible flows using a node‐based finite‐element method

An efficient parallel computing method for high‐speed compressible flows is presented. The numerical analysis of flows with shocks requires very fine computational grids and grid generation requires a great deal of time. In the proposed method, all computational procedures, from the mesh generation to the solution of a system of equations, can be performed seamlessly in parallel in terms of nodes. Local finite‐element mesh is generated robustly around each node, even for severe boundary shapes such as cracks. The algorithm and the data structure of finite‐element calculation are based on nodes, and parallel computing is realized by dividing a system of equations by the row of the global coefficient matrix. The inter‐processor communication is minimized by renumbering the nodal identification number using ParMETIS. The numerical scheme for high‐speed compressible flows is based on the two‐step Taylor–Galerkin method. The proposed method is implemented on distributed memory systems, such as an Alpha PC cluster, and a parallel supercomputer, Hitachi SR8000. The performance of the method is illustrated by the computation of supersonic flows over a forward facing step. The numerical examples show that crisp shocks are effectively computed on multiprocessors at high efficiency. Copyright © 2003 John Wiley & Sons, Ltd.     

Distortion‐immune nine‐node displacement‐based quadrilateral thick plate finite elements that satisfy constant‐bending patch test

We employ the linked interpolation concept to develop two higher‐order nine‐node quadrilateral plate finite elements with curved sides that pass the constant bending patch test for arbitrary node positions. The linked interpolation for the plate displacements is expanded with three bubble parameters to get polynomial completeness necessary to satisfy the patch test. In contrast to some other techniques, the elements developed in this way retain a symmetric stiffness matrix at a marginal computational expense at the element level. The new elements generated using this concept are tested on several examples with curved sides or some other kind of geometric distortion. Copyright © 2014 John Wiley & Sons, Ltd.     

Application of piece‐wise linear weight functions for 2D 8‐node quadrilateral element in contact problems

The present study is a continuation of our previous work with the aim to reduce problems caused by standard higher order elements in contact problems. The difficulties can be attributed to the inherent property of the Galerkin method which gives uneven distributions of nodal forces resulting in oscillating contact pressures. The proposed remedy is use of piece‐wise linear weight functions. The methods to establish stiffness and/or mass matrix for 8‐node quadrilateral element in 2D are presented, i.e. the condensing and direct procedures. The energy and nodal displacement error norms are also checked to establish the convergence ratio. Interpretation of calculated contact pressures is discussed. Two new 2D 8‐node quadrilateral elements, QUAD8C and QUAD8D, are derived and tested in many examples, which show their good performance in contact problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Non‐linear exact geometry 12‐node solid‐shell element with three translational degrees of freedom per node

This paper presents the finite rotation exact geometry (EG) 12‐node solid‐shell element with 36 displacement degrees of freedom. The term ‘EG’ reflects the fact that coefficients of the first and second fundamental forms of the reference surface and Christoffel symbols are taken exactly at each element node. The finite element formulation developed is based on the 9‐parameter shell model by employing a new concept of sampling surfaces (S‐surfaces) inside the shell body. We introduce three S‐surfaces, namely, bottom, middle and top, and choose nine displacements of these surfaces as fundamental shell unknowns. Such choice allows one to represent the finite rotation higher order EG solid‐shell element formulation in a very compact form and to derive the strain–displacement relationships, which are objective, that is, invariant under arbitrarily large rigid‐body shell motions in convected curvilinear coordinates. The tangent stiffness matrix is evaluated by using 3D analytical integration and the explicit presentation of this matrix is given. The latter is unusual for the non‐linear EG shell element formulation. Copyright © 2011 John Wiley & Sons, Ltd.     

Development of a 4‐node hybrid stress tetrahedral element using a node‐based smoothed finite element method

In this paper, a new 4‐node hybrid stress element is proposed using a node‐based smoothing technique of tetrahedral mesh. The conditions for hybrid stress field required are summarized, and the field should be continuous for better performance of a constant‐strain tetrahedral element. Nodal stress is approximated by the node‐based smoothing technique, and the stress field is interpolated with standard shape functions. This stress field is linear within each element and continuous across elements. The stress field is expressed by nodal displacements and no additional variables. The element stiffness matrix is calculated using the Hellinger‐Reissner functional, which guarantees the strain field from displacement field to be equal to that from the stress field in a weak sense. The performance of the proposed element is verified by through several numerical examples.