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.     

A shape‐free 8‐node plane element unsymmetric analytical trial function method

The unsymmetric FEM is one of the effective techniques for developing finite element models immune to various mesh distortions. However, because of the inherent limitation of the metric shape functions, the resulting element models exhibit rotational frame dependence and interpolation failure under certain conditions. In this paper, by introducing the analytical trial function method used in the hybrid stress‐function element method, an effort was made to naturally eliminate these defects and improve accuracy. The key point of the new strategy is that the monomial terms (the trial functions) in the assumed metric displacement fields are replaced by the fundamental analytical solutions of plane problems. Furthermore, some rational conditions are imposed on the trial functions so that the assumed displacement fields possess fourth‐order completeness in Cartesian coordinates. The resulting element model, denoted by US‐ATFQ8, can still work well when interpolation failure modes for original unsymmetric element occur, and provide the invariance for the coordinate rotation. Numerical results show that the exact solutions for constant strain/stress, pure bending and linear bending problems can be obtained by the new element US‐ATFQ8 using arbitrary severely distorted meshes, and produce more accurate results for other more complicated problems. Copyright © 2012 John Wiley & Sons, Ltd.     

Rotation‐free triangular shell element using node‐based smoothed finite element method

In this paper, a triangular thin flat shell element without rotation degrees of freedom is proposed. In the Kirchhoff hypothesis, the first derivative of the displacement must be continuous because there are second‐order differential terms of the displacement in the weak form of the governing equations. The displacement is expressed as a linear function and the nodal rotation is defined using node‐based smoothed finite element method. The rotation field is approximated using the nodal rotation and linear shape functions. This rotation field is linear in an element and continuous between elements. The curvature is defined by differentiating the rotation field, and the stiffness is calculated from the curvature. A hybrid stress triangular membrane element was used to construct the shell element. The penalty technique was used to apply the rotation boundary conditions. The proposed element was verified through several numerical examples.     

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.     

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.     

A 3‐node C0 triangular element for the modified couple stress theory based on the smoothed finite element method

In this paper, a 3‐node C⁰ triangular element for the modified couple stress theory is proposed. Unlike the classical continuum theory, the second‐order derivative of displacement is included in the weak form of the equilibrium equations. Thus, the first‐order derivative of displacement, such as the rotation, should be approximated by a continuous function. In the proposed element, the derivative of the displacement is defined at a node using the node‐based smoothed finite element method. The derivative fields, continuous between elements and linear in an element, are approximated with the shape functions in element. Both the displacement field and the derivative field of displacement are expressed in terms of the displacement degree of freedom only. The element stiffness matrix is calculated using the newly defined derivative field. The performance of the proposed element is evaluated through various numerical examples.     

On inter‐element forces in the FEM‐displacement formulation,and implications for stress recovery

This paper describes a new technique for the determination of the inter‐element forces and tractions, as well as stress state at nodes, as a post‐processing step after the solution of standard FE‐displacement calculation. The work is motivated in the context of a broader development of a procedure to simulate fracture processes using a discrete approach without the need of double‐noded interface elements. The technique, easily implementable, is based on the double minimization of an objective function, representing the error between the inter‐element stress tractions and the projection of the best‐fit stress tensor T along the planes of the interfaces converging at an element corner node. The formulation is illustrated with some basic examples in which the resulting stress tensors and inter‐element forces are compared to theoretical solutions and to the results obtained by using a traditional stress average smoothing method. Copyright © 2005 John Wiley & Sons, Ltd.     

A systematic construction of B‐bar functions for linear and non‐linear mixed‐enhanced finite elements for plane elasticity problems

In a previous paper a modified Hu–Washizu variational formulation has been used to derive an accurate four node plane strain/stress finite element denoted QE2. For the mixed element QE2 two enhanced strain terms are used and the assumed stresses satisfy the equilibrium equations a priori for the linear elastic case. In this paper an alternative approach is discussed. The new formulation leads to the same accuracy for linear elastic problems as the QE2 element; however it turns out to be more efficient in numerical simulations, especially for large deformation problems. Using orthogonal stress and strain functions we derive B̄ functions which avoid numerical inversion of matrices. The B̄ ‐strain matrix is sparse and has the same structure as the strain matrix B obtained from a compatible displacement field. The implementation of the derived mixed element is basically the same as the one for a compatible displacement element. The only difference is that we have to compute a B̄ ‐strain matrix instead of the standard B ‐matrix. Accordingly, existing subroutines for a compatible displacement element can be easily changed to obtain the mixed‐enhanced finite element which yields a higher accuracy than the Q4 and QM6 elements. Copyright © 1999 John Wiley & Sons, Ltd.     

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 large deformation solid‐shell concept based on reduced integration with hourglass stabilization

In this paper a new eight‐node (brick) solid‐shell finite element formulation based on the concept of reduced integration with hourglass stabilization is presented. The work focuses on static problems. The starting point of the derivation is the three‐field variational functional upon which meanwhile established 3D enhanced strain concepts are based. Important additional assumptions are made to transfer the approach into a powerful solid‐shell. First of all, a Taylor expansion of the first Piola–Kirchhoff stress tensor with respect to the normal through the centre of the element is carried out. In this way the stress becomes a linear function of the shell surface co‐ordinates whereas the dependence on the thickness co‐ordinate remains non‐linear. Secondly, the Jacobian matrix is replaced by its value in the centre of the element. These two assumptions lead to a computationally efficient shell element which requires only two Gauss points in the thickness direction (and one Gauss point in the plane of the shell element). Additionally three internal element degrees‐of‐freedom have to be determined to avoid thickness locking. One important advantage of the element is the fact that a fully three‐dimensional stress state can be modelled without any modification of the constitutive law. The formulation has only displacement degrees‐of‐freedom and the geometry in the thickness direction is correctly displayed. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

A simple unsymmetric 4-node 12-DOF membrane element for the modified couple stress theory

In this work, the recently proposed unsymmetric 4-node 12-DOF (degree-of-freedom) membrane element (Shang and Ouyang, Int J Numer Methods Eng 113(10): 1589-1606, 2018), which has demonstrated excellent performance for the classical elastic problems, is further extended for the modified couple stress theory, to account for the size effect of materials. This is achieved via two formulation developments. Firstly, by using the penalty function method, the kinematic relations between the element's nodal drilling DOFs and the true physical rotations are enforced. Consequently, the continuity requirement for the modified couple stress theory is satisfied in weak sense, and the symmetric curvature test function can be easily derived from the gradients of the drilling DOFs. Secondly, the couple stress field that satisfies a priori the related equilibrium equations is adopted as the energy conjugate trial function to formulate the element for the modified couple stress theory. As demonstrated by a series of benchmark tests, the new element can efficiently capture the size-dependent responses of materials and is robust to mesh distortions. Moreover, as the new element uses only three conventional DOFs per node, it can be readily incorporated into the standard finite element program framework and commonly available finite element programs.     

A novel singular node‐based smoothed finite element method (NS‐FEM) for upper bound solutions of fracture problems

It is well known that the lower bound to exact solutions in linear fracture problems can be easily obtained by the displacement compatible finite element method (FEM) together with the singular crack tip elements. It is, however, much more difficult to obtain the upper bound solutions for these problems. This paper aims to formulate a novel singular node‐based smoothed finite element method (NS‐FEM) to obtain the upper bound solutions for fracture problems. In the present singular NS‐FEM, the calculation of the system stiffness matrix is performed using the strain smoothing technique over the smoothing domains (SDs) associated with nodes, which leads to the line integrations using only the shape function values along the boundaries of the SDs. A five‐node singular crack tip element is used within the framework of NS‐FEM to construct singular shape functions via direct point interpolation with proper order of fractional basis. The mix‐mode stress intensity factors are evaluated using the domain forms of the interaction integrals. The upper bound solutions of the present singular NS‐FEM are demonstrated via benchmark examples for a wide range of material combinations and boundary conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

An enhanced cell‐based smoothed finite element method for the analysis of Reissner–Mindlin plate bending problems involving distorted mesh

In this work, an enhanced cell‐based smoothed finite element method (FEM) is presented for the Reissner–Mindlin plate bending analysis. The smoothed curvature computed by a boundary integral along the boundaries of smoothing cells in original smoothed FEM is reformulated, and the relationship between the original approach and the present method in curvature smoothing is established. To improve the accuracy of shear strain in a distorted mesh, we span the shear strain space over the adjacent element. This is performed by employing an edge‐based smoothing technique through a simple area‐weighted smoothing procedure on MITC4 assumed shear strain field. A three‐field variational principle is utilized to develop the mixed formulation. The resultant element formulation is further reduced to a displacement‐based formulation via an assumed strain method defined by the edge‐smoothing technique. As the result, a new formulation consisting of smoothed curvature and smoothed shear strain interpolated by the standard transverse displacement/rotation fields and smoothing operators can be shown to improve the solution accuracy in cell‐based smoothed FEM for Reissner–Mindlin plate bending analysis. Several numerical examples are presented to demonstrate the accuracy of the proposed formulation.Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

A solid‐like shell element allowing for arbitrary delaminations

In this contribution a new finite element is presented for the simulation of delamination growth in thin‐layered composite structures. The element is based on a solid‐like shell element: a volume element that can be used for very thin applications due to a higher‐order displacement field in the thickness direction. The delamination crack can occur at arbitrary locations and is incorporated in the element as a jump in the displacement field by using the partition of unity property of finite element shape functions. The kinematics of the element as well as the finite element formulation are described. The performance of the element is demonstrated by means of two examples. Copyright © 2003 John Wiley & Sons, Ltd.     

Bubble‐enhanced smoothed finite element formulation: a variational multi‐scale approach for volume‐constrained problems in two‐dimensional linear elasticity

This paper presents a bubble‐enhanced smoothed finite element formulation for the analysis of volume‐constrained problems in two‐dimensional linear elasticity. The new formulation is derived based on the variational multi‐scale approach in which unequal order displacement‐pressure pairs are used for the mixed finite element approximation and hierarchical bubble function is selected for the fine‐scale displacement approximation. An area‐weighted averaging scheme is employed for the two‐scale smoothed strain calculation under the framework of edge‐based smoothed FEM. The smoothed fine‐scale solution is shown to naturally contain the stress field jump of the smoothed coarse‐scale solution across the boundary of edge‐based smoothing domain and thus provides the possibility to stabilize the global solution for volume‐constrained problems. A global monolithic solution strategy is employed, and the fine‐scale solution is solved without the consideration of approximating the strong form of the fine‐scale equation. Several numerical examples are analyzed to demonstrate the accuracy of the present formulation. Copyright © 2014 John Wiley & Sons, Ltd.     

A mixed‐enhanced finite‐element for the analysis of laminated composite plates

This paper presents a new 4‐node finite‐element for the analysis of laminated composite plates. The element is based on a first‐order shear deformation theory and is obtained through a mixed‐enhanced approach. In fact, the adopted variational formulation includes as variables the transverse shear as well as enhanced incompatible modes introduced to improve the in‐plane deformation. The problem is then discretized using bubble functions for the rotational degrees of freedom and functions linking the transverse displacement to the rotations. The proposed element is locking free, it does not have zero energy modes and provides accurate in‐plane/out‐of‐plane deformations. Furthermore, a procedure for the computation of the through‐the‐thickness shear stresses is discussed, together with an iterative algorithm for the evaluation of the shear correction factors. Several applications are investigated to assess the features and the performances of the proposed element. Results are compared with analytical solutions and with other finite‐element solutions. Copyright © 1999 John Wiley & Sons, Ltd.     

Time‐domain hybrid formulations for wave simulations in three‐dimensional PML‐truncated heterogeneous media

We are concerned with the numerical simulation of wave motion in arbitrarily heterogeneous, elastic, perfectly‐matched‐layer‐(PML)‐truncated media. We extend in three dimensions a recently developed two‐dimensional formulation, by treating the PML via an unsplit‐field, but mixed‐field, displacement‐stress formulation, which is then coupled to a standard displacement‐only formulation for the interior domain, thus leading to a computationally cost‐efficient hybrid scheme. The hybrid treatment leads to, at most, third‐order in time semi‐discrete forms. The formulation is flexible enough to accommodate the standard PML, as well as the multi‐axial PML. We discuss several time‐marching schemes, which can be used à la carte, depending on the application: (a) an extended Newmark scheme for third‐order in time, either unsymmetric or fully symmetric semi‐discrete forms; (b) a standard implicit Newmark for the second‐order, unsymmetric semi‐discrete forms; and (c) an explicit Runge–Kutta scheme for a first‐order in time unsymmetric system. The latter is well‐suited for large‐scale problems on parallel architectures, while the second‐order treatment is particularly attractive for ready incorporation in existing codes written originally for finite domains. We compare the schemes and report numerical results demonstrating stability and efficacy of the proposed formulations. Copyright © 2014 John Wiley & Sons, Ltd.