Observations regarding assumed-stress hybrid plate elements

Edge displacement of assumed-stress hybrid plate elements are commonly chosen as either linear, quadratic or cubic functions of the edge co-ordinate. The present paper examines the effect of this choice on the acceptability, accuracy, and possible zero energy modes of elements based on a linear stress field. It is concluded that triangular elements may be preferable to quadrilaterals, and that linear edge displacements are not desirable.     

Stability and accuracy of finite element methods for flow acoustics. II: Two‐dimensional effects

This is the second of two articles that focus on the dispersion properties of finite element models for acoustic propagation on mean flows. We consider finite element methods based on linear potential theory in which the acoustic disturbance is modelled by the convected Helmholtz equation, and also those based on a mixed Galbrun formulation in which acoustic pressure and Lagrangian displacement are used as discrete variables. The current paper focuses on the effects of numerical anisotropy which are associated with the orientation of the propagating wave to the mean flow and to the grid axes. Conditions which produce aliasing error in the Helmholtz formulation are of particular interest. The 9‐noded Lagrangian element is shown to be superior to the more commonly used 8‐noded serendipity element. In the case of the Galbrun elements, the current analysis indicates that isotropic meshes generally reduce numerical error of triangular elements and that higher order mixed quadrilaterals are generally less effective than an equivalent mesh of lower order triangles. Copyright © 2005 John Wiley & Sons, Ltd.     

The effects of element shape on the critical time step in explicit time integrators for elasto‐dynamics

In this paper, the effects of element shape on the critical time step are investigated. The common rule‐of‐thumb, used in practice, is that the critical time step is set by the shortest distance within an element divided by the dilatational (compressive) wave speed, with a modest safety factor. For regularly shaped elements, many analytical solutions for the critical time step are available, but this paper focusses on distorted element shapes. The main purpose is to verify whether element distortion adversely affects the critical time step or not. Two types of element distortion will be considered, namely aspect ratio distortion and angular distortion, and two particular elements will be studied: four‐noded bilinear quadrilaterals and three‐noded linear triangles. The maximum eigenfrequencies of the distorted elements are determined and compared to those of the corresponding undistorted elements. The critical time steps obtained from single element calculations are also compared to those from calculations based on finite element patches with multiple elements. Copyright © 2014 John Wiley & Sons, Ltd.     

Determination of coefficients of the crack tip asymptotic field by fractal hybrid finite elements

This paper applies the fractal finite element method (FFEM) together with 9-node Lagrangian hybrid elements to the calculation of linear elastic crack tip fields. An explicit stabilization scheme is employed to suppress the spurious kinematic modes of the sub-integrated Lagrangian element. An extensive convergence study has been conducted to examine the effects of the similarity ratio, reduced integration and the type of elements on the accuracy and stability of the numerical solutions. It is concluded that (i) a similarity ratio close to unity should be used to construct the fractal mesh, (ii) sub-integrated Lagrangian elements with occasionally unstable behaviour should not be used, and (iii) the good accuracy (with differences less than 0.5% with existing available solutions) and stability over a wide range of numerical tests support the use of fractal hybrid finite elements for determining crack tip asymptotic fields.     

Optimization of stress modes by energy compatibility for 4‐node hybrid quadrilaterals

Optimal hybrid stress quadrilaterals can be obtained by adopting appropriate stresses and displacements, and satisfying the energy compatibility condition is shown to be an ultimate key to obtaining optimal stress modes. By using compatible isoparametric bilinear (Q₄) displacements and 5‐parameter energy compatible stresses of the combined hybrid finite element CH(0‐1), a robust 4‐node plane stress element ECQ₄ is derived. Equivalence to another hybrid stress element LQ₆ with 9‐parameter complete linear stresses based on a modified Hellinger–Reissner principle is established. A convergence analysis is given and numerical experiments show that elements ECQ₄/LQ₆ have high performance, i.e. are accurate at coarse meshes, insensitive to mesh distortions and free from locking. Copyright © 2003 John Wiley & Sons, Ltd.     

A theoretically improved and easily implemented version of the Ahmad Thick Shell element

An implementation of the Ahmad Thick Shell element using vector manipulation, reduced integration and the incorporation of the missing term in the approximating polynomial is presented. The first two aspects are merely applications. The theory is presented elsewhere^6,7,11,18 and will not be repeated for the sake of conciseness. The third aspect was achieved by adding to the original nodal configuration a central node having only one degree-of-freedom. That node is then eliminated to preserve the original number of degrees-of-freedom. The Ahmad thick shell element is adequately presented in finite element literature.^1-3 The co-ordinate definition and the displacement field are retained and will not be dealt with henceforth. However, a different concept of the strain definition due to Irons is used and the bending terms are included in the modulus matrix. The element like all second generation isoparametric elements is easily implemented through a shape function subroutine suitable for quadrilaterals. It reproduces exactly rigid body motion for any combination of elements of any geometry even for elements withcurved sides and variable thickness. Although Ahmad's thick shell element passed the patch test for unequal parallelograms but not for quadrilaterals the version presented passes that test for quadrilaterals. In spite of its credentials, the new version denoted by A has a major setback. the presence of one too many spurious mechanism is reported.     

A simple cubic displacement element for plate bending

Early attempts to construct a triangular finite element for plate bending problems from a compatible cubic displacement field are not entirely satisfactory. The present paper shows how an accurate plate element can be achieved using independent cubic polynomial assumptions for the internal and boundary displacements in conjunction with a modified potential energy principle. This approach yields a simple algebraic formulation with favourable connection quantities at the element vertices which will appeal to practical users of the conventional finite element displacement method. Moreover, in Appendix I it is shown that the cubic element is identical to a previous hybrid stress element with linear internal bending and twisting moments and cubic boundary displacements. The stresses obtained from the former hybrid finite element solution therefore satisfy the strain compatibility conditions exactly. This remarkable result has an important significance in the theory of hybrid finite elements.     

The hybrid-Trefftz finite element model and its application to plate bending

This paper presents a new hybrid element approach and applies it to plate bending. In contrast to more conventional models, the formulation is based on displacement fields which fulfil a priori the non-homogeneous Lagrange equation (Trefftz method). The interelement continuity is enforced by using a stationary principle together with an independent interelement displacement. The final unknowns are the nodal displacements and the elements may be implemented without any difficulty in finite element libraries of standard finite element programs. The formulation only calls for integration along the element boundaries which enables arbitrary polygonal or even curve-sided elements to be generated. Where relevant, known local solutions in the vicinity of a singularity or stress concentration may be used as an optional expansion basis to obtain, for example, particular singular corner elements, elements presenting circular holes, etc. Thus a high degree of accuracy may be achieved without a troublesome mesh refinement. Another important advantage of the formulation is the possibility of generating by a single element subroutine a large number of various elements (triangles, quadrilaterals, etc.), presenting an increasing degree of accuracy. The paper summarizes the results of numerical studies and shows the excellent accuracy and efficiency of the new elements. The conclusions present some ideas concerning the adaptive version of the new elements, extension to nonlinear problems and some other developments.     

Theoretical aspects of the smoothed finite element method (SFEM)

This paper examines the theoretical bases for the smoothed finite element method (SFEM), which was formulated by incorporating cell‐wise strain smoothing operation into standard compatible finite element method (FEM). The weak form of SFEM can be derived from the Hu–Washizu three‐field variational principle. For elastic problems, it is proved that 1D linear element and 2D linear triangle element in SFEM are identical to their counterparts in FEM, while 2D bilinear quadrilateral elements in SFEM are different from that of FEM: when the number of smoothing cells (SCs) of the elements equals 1, the SFEM solution is proved to be ‘variationally consistent’ and has the same properties with those of FEM using reduced integration; when SC approaches infinity, the SFEM solution will approach the solution of the standard displacement compatible FEM model; when SC is a finite number larger than 1, the SFEM solutions are not ‘variationally consistent’ but ‘energy consistent’, and will change monotonously from the solution of SFEM (SC = 1) to that of SFEM (SC → ∞). It is suggested that there exists an optimal number of SC such that the SFEM solution is closest to the exact solution. The properties of SFEM are confirmed by numerical examples. Copyright © 2006 John Wiley & Sons, Ltd.     

Improved stress evaluation under thermal load for simple finite elements

The method of initial strain is replaced by a method of initial displacements and applied to simple four-noded membrane and shell elements. The thermal stresses in nodes are improved for quadrilaterals and hexahedrons, which have a discrepancy between interpolation functions for temperature (linear) and strains (constant). To get initial displacements, the singular element stiffness matrix has to be inverted. Minor changes in subroutine MINV³ (enclosed version STDMNV) were made to do this.     

Comparison of implicit and explicit hybridizable discontinuous Galerkin methods for the acoustic wave equation

We describe implicit and explicit formulations of the hybridizable discontinuous Galerkin method for the acoustic wave equation based on state‐of‐the‐art numerical software and quantify their efficiency for realistic application settings. In the explicit scheme, the trace of the acoustic pressure is computed from the solution on the two elements adjacent to the face at the old time step. Tensor product shape functions for quadrilaterals and hexahedra evaluated with sum factorization are used to ensure low operation counts. For applying the inverse mass matrix of Lagrangian shape functions with full Gaussian quadrature, a new tensorial technique is proposed. As time propagators, diagonally implicit and explicit Runge–Kutta methods are used, respectively. We find that the computing time per time step is 25 to 200 times lower for the explicit scheme, with an increasing gap in three spatial dimensions and for higher element degrees. Our experiments on realistic 3D wave propagation with variable material parameters in a photoacoustic imaging setting show an improvement of two orders of magnitude in terms of time to solution, despite stability restrictions on the time step of the explicit scheme. Operation counts and a performance model to predict performance on other computer systems accompany our results. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Hybrid and mixed-penalty finite elements for 3-D analysis of soft hydrated tissue

Solution of biomechanics problems involving three-dimensional (3-D) behaviour of soft tissue on geometries representative of such tissue in vivo will require the use of numerical methods. Toward this end, a pair of tetrahedral finite elements has been developed. The equations which are used to model the tissue behaviour for both elements are those commonly known as the linear biphasic equations. This model assumes that hydrated soft tissue is a mixture of two incompressible, immiscible phases, and employs mixture theory to derive governing equations for its mechanical behaviour. The finite element techniques applied to these equations for the two elements are the mixed-penalty method and the hybrid method. Both elements are described here, and the special requirements for 3-D analysis are discussed. Results obtained by solving canonical problems in two and three dimensions using both elements are presented and compared. Both elements are found to produce excellent results. The hybrid element is also noted to have advantages for non-linear analyses involving finite deformation which will require solution in the future.     

On using degenerated solid shell elements in adaptive refinement analysis

Three different degenerated shell elements are studied in an adaptive refinement procedure for the solution of shell problems. The stress recovery procedure expressed in a convective patch co‐ordinate system is used for the construction of continuous smoothed stress fields for the a posteriori error estimation. The performance of the stress recovery procedure, the error estimator and the adaptive refinement strategy are tested by solving three benchmark shell problems. It is found that when adaptive refinement is used, the adverse effects of boundary layers and stress singularities are eliminated and all the elements tested are able to achieve their optimal convergence rates. It is also found that the accuracy of the shell elements increases with the number of polynomial terms included in the stress and strain approximations. In addition, if complete Lagrangian polynomial terms are used, the element will be less sensitive to shape distortion than the one in which only complete polynomial terms are employed. Copyright © 1999 John Wiley & Sons, Ltd.     

A coupled two‐scale shell model with applications to layered structures

In this paper a coupled two‐scale shell model is presented. A variational formulation and associated linearization for the coupled global–local boundary value problem is derived. For small strain problems, various numerical solutions are computed within the so‐called FE ² method. The discretization of the shell is performed with quadrilaterals, whereas the local boundary value problems at the integration points of the shell are discretized using 8‐noded or 27‐noded brick elements or so‐called solid shell elements. At the bottom and top surface of the representative volume element stress boundary conditions are applied, whereas at the lateral surfaces the in‐plane displacements are prescribed. For the out‐of‐plane displacements link conditions are applied. The coupled nonlinear boundary value problems are simultaneously solved within a Newton iteration scheme. With an important test, the correct material matrix for the stress resultants assuming linear elasticity and a homogeneous continuum is verified.Copyright © 2013 John Wiley & Sons, Ltd.     

A partition of unity‐based ‘FE–Meshfree’ QUAD4 element for geometric non‐linear analysis

The recently published ‘FE–Meshfree’ QUAD4 element is extended to geometrical non‐linear analysis. The shape functions for this element are obtained by combining meshfree and finite element shape functions. The concept of partition of unity (PU) is employed for the purpose. The new shape functions inherit their higher order completeness properties from the meshfree shape functions and the mesh‐distortion tolerant compatibility properties from the finite element (FE) shape functions. Updated Lagrangian formulation is adopted for the non‐linear solution. Several numerical example problems are solved and the performance of the element is compared with that of the well‐known Q4, QM6 and Q8 elements. The results show that, for regular meshes, the performance of the element is comparable to that of QM6 and Q8 elements, and superior to that of Q4 element. For distorted meshes, the present element has better mesh‐distortion tolerance than Q4, QM6 and Q8 elements. Copyright © 2009 John Wiley & Sons, Ltd.     

Analytical integration formulae for linear isoparametric finite elements

Exact and approximate analytical expressions can be derived for integrals arising in finite element methods, employing isoparametric linear quadrilaterals in two space dimensions with bilinear basis functions. The formulae associated with rectangular elements, arbitrarily oriented in space, can be shown to be a special case. The proposed method provides considerable savings in computational effort, in comparison with a numerical method that employs Gaussian quadrature procedures. In addition, the method, when applied to a quadrilateral inscribable in a circle, can be shown to produce better accuracy than the associated (2 × 2) Gaussian quadrature formulae.     

A hybrid equilibrium element for folded plate and shell structures

This paper extends hybrid equilibrium formulation concepts, previously used with success for planar problems, to the analysis of folded plates and curved shells. A 2D hybrid equilibrium flat shell quadrilateral element is formulated for linear analysis, where detailed consideration is given to the implication of slope discontinuity when the element is used for non‐planar domains. Benchmark plate bending, folded plate and curved shell problems are modelled using equilibrium and conforming elements for comparison. In models of the latter two problems, torsional moments may be released along lines of slope discontinuity, and the effects of this assumption for the folded plate are studied by analysing a third type of model composed of 3D solid brick elements. The comparisons demonstrate an excellent performance from the new hybrid equilibrium analysis method for folded plates and curved shells. Copyright © 2013 John Wiley & Sons, Ltd.     

Structural‐acoustic coupling on non‐conforming meshes with quadratic shape functions

Fully coupled finite element/boundary element models are a popular choice when modelling structures that are submerged in heavy fluids. To achieve coupling of subdomains with non‐conforming discretizations at their common interface, the coupling conditions are usually formulated in a weak sense. The coupling matrices are evaluated by integrating products of piecewise polynomials on independent meshes. The case of interfacing elements with linear shape functions on unrelated meshes has been well covered in the literature. This paper presents a solution to the problem of evaluating the coupling matrix for interfacing elements with quadratic shape functions on unrelated meshes. The isoparametric finite elements have eight nodes (Serendipity) and the discontinuous boundary elements have nine nodes (Lagrange). Results using linear and quadratic shape functions on conforming and non‐conforming meshes are compared for an example of a fluid‐loaded point‐excited sphere. It is shown that the coupling error decreases when quadratic shape functions are used. Copyright © 2012 John Wiley & Sons, Ltd.     

A hybrid finite element formulation of the linear biphasic equations for hydrated soft tissue

A hybrid finite element method has been developed for application to the linear biphasic model of soft tissues. The biphasic model assumes that hydrated soft tissue is a mixture of two incompressible, immiscible phases, one solid and one fluid, and employs mixture theory to derive governing equations for its mechanical behaviour. These equations are time dependent, involving both fluid and solid velocities and solid displacement, and will be solved by spatial finite element and temporal finite difference approximation. The first step in the derivation of this hybrid method is application of a finite difference rule to the solid phase, thus obtaining equations with only velocities at discrete times as primary variables. A weighted residual statement of the temporally discretized governing equations, employing C° continuous interpolations of the solid and fluid phase velocities and discontinuous interpolations of the pore pressure and elastic stress, is then derived. The stress and pressure functions are chosen so that the total momentum equation of the mixture is satisfied; they are jointly referred to as an equilibrated stress and pressure field. The corresponding weighting functions are chosen to satisfy a relationship analogous to this equilibrium relation. The resulting matrix equations are symmetric. As an illustration of the hybrid biphasic formulation, six-noded triangular elements with complete linear, several incomplete quadratic, and complete quadratic stress and pressure fields in element local co-ordinates are developed for two dimensional analysis and tested against analytical solutions and a mixed-penalty finite element formulation of the same equations. The hybrid method is found to be robust and produce excellent results; preferred elements are identified on the basis of these results.