Mean‐strain 10‐node tetrahedron with energy‐sampling stabilization

In this study, a new mean‐strain 10‐node tetrahedral element is developed using energy‐sampling stabilization. The proposed 10‐node tetrahedron is composed of several four‐node linear tetrahedral elements, four tetrahedra in the corners and four tetrahedra that tile the central octahedron in three possible sets of four‐node linear tetrahedra, corresponding to three different choices for the internal diagonal. The assumed strains are calculated from mean ‘basis function gradients.’ The energy‐sampling technique introduced previously for removing zero‐energy modes in the mean‐strain hexahedron is adapted for the present element: the stabilization energy is evaluated on the four‐corner tetrahedra. The proposed element naturally leads to a lumped‐mass matrix and does not have unphysical low‐energy vibration modes. For simplicity, we limit our developments to linear elasticity with compressible and nearly incompressible material. The numerical tests demonstrate that the present element performs well compared with the classical 10‐node tetrahedral elements for shell and plate structures, and nearly incompressible materials. Copyright © 2016 John Wiley & Sons, Ltd.     

Stabilization of mixed tetrahedral elements at large deformations

This paper presents stabilized mixed finite element formulations for tetrahedral elements at large deformations using volume and area bubble functions. To this end, the corresponding weak formulations are derived for the standard two‐field method, the method of incompatible modes and the enhanced strain method. Then, the weak formulations will be linearized. Furthermore, the matrix formulations for the weak formulations and its linearizations are summarized. The numerical results for incompressible rubber‐like materials using a Neo‐Hookean material law show the locking‐free performance and the drastic damping of the stresses for the new stabilized tetrahedral elements in finite deformation problems. This paper is an extension of the works published by the authors regarding small deformation problems for linear elasticity and plasticity. Copyright © 2011 John Wiley & Sons, Ltd.     

Enhanced transverse shear strain shell formulation applied to large elasto‐plastic deformation problems

In this work, a previously proposed Enhanced Assumed Strain (EAS) finite element formulation for thin shells is revised and extended to account for isotropic and anisotropic material non‐linearities. Transverse shear and membrane‐locking patterns are successfully removed from the displacement‐based formulation. The resultant EAS shell finite element does not rely on any other mixed formulation, since the enhanced strain field is designed to fulfil the null transverse shear strain subspace coming from the classical degenerated formulation. At the same time, a minimum number of enhanced variables is achieved, when compared with previous works in the field. Non‐linear effects are treated within a local reference frame affected by the rigid‐body part of the total deformation. Additive and multiplicative update procedures for the finite rotation degrees‐of‐freedom are implemented to correctly reproduce mid‐point configurations along the incremental deformation path, improving the overall convergence rate. The stress and strain tensors update in the local frame, together with an additive treatment of the EAS terms, lead to a straightforward implementation of non‐linear geometric and material relations. Accuracy of the implemented algorithms is shown in isotropic and anisotropic elasto‐plastic problems. Copyright © 2004 John Wiley & Sons, Ltd.     

A meshfree‐enriched finite element method for compressible and near‐incompressible elasticity

In this paper, a two‐dimensional displacement‐based meshfree‐enriched FEM (ME‐FEM) is presented for the linear analysis of compressible and near‐incompressible planar elasticity. The ME‐FEM element is established by injecting a first‐order convex meshfree approximation into a low‐order finite element with an additional node. The convex meshfree approximation is constructed using the generalized meshfree approximation method and it possesses the Kronecker‐delta property on the element boundaries. The gradient matrix of ME‐FEM element satisfies the integration constraint for nodal integration and the resultant ME‐FEM formulation is shown to pass the constant stress test for the compressible media. The ME‐FEM interpolation is an element‐wise meshfree interpolation and is proven to be discrete divergence‐free in the incompressible limit. To prevent possible pressure oscillation in the near‐incompressible problems, an area‐weighted strain smoothing scheme incorporated with the divergence‐free ME‐FEM interpolation is introduced to provide the smoothing on strains and pressure. With this smoothed strain field, the discrete equations are derived based on a modified Hu–Washizu variational principle. Several numerical examples are presented to demonstrate the effectiveness of the proposed method for the compressible and near‐incompressible problems. Copyright © 2012 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.     

P1‐nonconforming quadrilateral finite element for topology optimization

This investigation focuses on an alternative approach to topology optimization problems involving incompressible materials using the P1‐nonconforming finite element. Instead of using the mixed displacement‐pressure formulation, a pure displacement‐based approach can be employed for finite element formulation owing to the Poisson locking‐free property of the P1‐nonconforming element. Moreover, because the P1‐nonconforming element has linear shape functions that are defined at element vertices, it has considerably fewer degrees of freedom than other quadrilateral nonconforming elements and its implementation is as simple as that of the conforming bilinear element. Various problems dealing with incompressible materials and pressure‐loaded structures found in published works are solved to verify the applicability of the proposed method. The application of the method is extended to the optimal design of fluid channels in the Stokes flow. This is done by expressing pressure in terms of volumetric strain rates and developing a velocity‐field‐only finite element formulation. The optimization results obtained from all the problems considered in this study are in close agreement with those found in the literature. Copyright © 2010 John Wiley & Sons, Ltd.     

Development of shear locking‐free shell elements using an enhanced assumed strain formulation

The degenerated approach for shell elements of Ahmad and co‐workers is revisited in this paper. To avoid transverse shear locking effects in four‐node bilinear elements, an alternative formulation based on the enhanced assumed strain (EAS) method of Simo and Rifai is proposed directed towards the transverse shear terms of the strain field. In the first part of the work the analysis of the null transverse shear strain subspace for the degenerated element and also for the selective reduced integration (SRI) and assumed natural strain (ANS) formulations is carried out. Locking effects are then justified by the inability of the null transverse shear strain subspace, implicitly defined by a given finite element, to properly reproduce the required displacement patterns. Illustrating the proposed approach, a remarkably simple single‐element test is described where ANS formulation fails to converge to the correct results, being characterized by the same performance as the degenerated shell element. The adequate enhancement of the null transverse shear strain subspace is provided by the EAS method, enforcing Kirchhoff hypothesis for low thickness values and leading to a framework for the development of shear‐locking‐free shell elements. Numerical linear elastic tests show improved results obtained with the proposed formulation. Copyright © 2001 John Wiley & Sons, Ltd.     

Enrichment of enhanced assumed strain approximations for representing strong discontinuities: addressing volumetric incompressibility and the discontinuous patch test

We present a geometrically non‐linear assumed strain method that allows for the presence of arbitrary, intra‐finite element discontinuities in the deformation map. Special attention is placed on the coarse‐mesh accuracy of these methods and their ability to avoid mesh locking in the incompressible limit. Given an underlying mesh and an arbitrary failure surface, we first construct an enriched approximation for the deformation map with the non‐linear analogue of the extended finite element method (X‐FEM). With regard to the richer space of functions spanned by the gradient of the enriched approximation, we then adopt a broader interpretation of variational consistency for the construction of the enhanced strain. In particular, in those elements intersected by the failure surface, we construct enhanced strain approximations which are orthogonal to piecewise‐constant stress fields. Contrast is drawn with existing strong discontinuity approaches where the enhanced strain variations in localized elements were constructed to be orthogonal to constant nominal stress fields. Importantly, the present formulation gives rise to a symmetric tangent stiffness matrix, even in localized elements. The present modification also allows for the satisfaction of a discontinuous patch test, wherein two different constant stress fields (on each side of the failure surface) lie in the solution space. We demonstrate how the proposed modifications eliminate spurious stress oscillations along the failure surface, particularly for nearly incompressible material response. Additional numerical examples are provided to illustrate the efficacy of the modified method for problems in hyperelastic fracture mechanics. Copyright © 2003 John Wiley & Sons, Ltd.     

Assumed strain formulation for the four-node quadrilateral with improved in-plane bending behaviour

A new assumed strain quadrilateral element with highly accurate in-plane bending behaviour is presented for plane stress and plane strain analysis. The basic idea of the formulation consists in identification of various modes of deformation and then in proper modification of the strain field in some of these modes. In particular, the strain operator corresponding to the in-plane bending modes is modified to simulate the strain field resulting from the assumptions usually made in structural mechanics. The modification of the strain field leads to the assumed strain operator on the element level. As a result, the so-called shear and membrane locking phenomena are alleviated. The element exhibits remarkable success in bending- dominated problems even when severely distorted and high aspect ratio meshes are used. Another advantage of the present assumed strain element is that locking for nearly incompressible materials is also mitigated. While this assumed strain element passes the patch test only for the parallelogram shapes, the element provides convergent solutions as long as the initially general form of the element approaches a parallelogram shape with the refinement of the mesh.     

Topology optimization with incompressible materials under small and finite deformations using mixed u/p elements

This paper focuses on topology optimization utilizing incompressible materials under both small‐ and finite‐deformation kinematics. To avoid the volumetric locking that accompanies incompressibility, linear and nonlinear mixed displacement/pressure (u/p) elements are utilized. A number of material interpolation schemes are compared, and a new scheme interpolating both Young's modulus and Poisson's ratio (E ‐ν interpolation) is proposed. The efficacy of this proposed scheme is demonstrated on a number of examples under both small‐ and finite‐deformation kinematics. Excessive mesh distortions that may occur under finite deformations are dealt with by extending a linear energy interpolation approach to the nonlinear u/p formulation and utilizing an adaptive update strategy. The proposed optimization framework is demonstrated to be effective through a number of representative examples.     

A new conforming quadrilateral plate bending element

A new four-node conforming quadrilateral element (NCQ) for plate bending is Described. The element is based on an earlier free formulation due to Bergan,^1,2 but extra twist variables are included to make it more accurate. The element geometry is in bilinear polynomials representation, while the displacement functions are described by slight modifications of bicubic polynomials which satisfy energy orthogonality and need not be force-orthogonal. It is shown that the coupling stiffness between the fundamental deformation modes and its higher-order deformation modes has been altered. Numerical examples indicate that the element gives a surprisingly accurate solution.     

Finite element formulation for modelling large deformations in elasto‐viscoplastic polycrystals

Anisotropic, elasto‐viscoplastic behaviour in polycrystalline materials is modelled using a new, updated Lagrangian formulation based on a three‐field form of the Hu‐Washizu variational principle to create a stable finite element method in the context of nearly incompressible behaviour. The meso‐scale is characterized by a representative volume element, which contains grains governed by single crystal behaviour. A new, fully implicit, two‐level, backward Euler integration scheme together with an efficient finite element formulation, including consistent linearization, is presented. The proposed finite element model is capable of predicting non‐homogeneous meso‐fields, which, for example, may impact subsequent recrystallization. Finally, simple deformations involving an aluminium alloy are considered in order to demonstrate the algorithm. Copyright © 2004 John Wiley & Sons, Ltd.     

Patch‐averaged assumed strain finite elements for stress analysis

A finite element model for linear‐elastic small deformation problems is presented. The formulation is based on a weighted residual that requires a priori the satisfaction of the kinematic equation. In this approach, an averaged strain‐displacement matrix is constructed for each node of the mesh by defining an appropriate patch of elements, yielding a smooth representation of strain and stress fields. Connections with traditional and similar procedure are explored. Linear quadrilateral four‐node and linear hexahedral eight‐node elements are derived. Various numerical tests show the accuracy and convergence properties of the proposed elements in comparison with extant finite elements and analytic solutions. Specific examples are also included to illustrate the ability to resist numerical locking in the incompressible limit and insensitive response in the presence of shape distortion. Furthermore, the numerical inf‐sup test is applied to a selection of problems to show the stability of the present formulation. Copyright © 2012 John Wiley & Sons, Ltd.     

Hybrid stress tetrahedral elements with Allman's rotational D.O.F.s

This paper presents two hybrid stress four‐node tetrahedron solid elements which are equipped with the rotational d.o.f.s proposed by Allman. Inasmuch Allman's rotation is employed, the elements are plagued by zero‐energy rotation modes which induce no strain. A modified Hellinger–Reissner functional that treats the rotation and the skew symmetric stress as independent fields is employed to formulate a stabilization scheme. Particular effort has been made to reduce the number of stress modes to minimum without sacrificing the frame invariance and proper rank of the element. The computational cost of the element is reduced by adopting orthogonal constant and non‐constant symmetric stress modes. Numerical benchmark tests indicate that accuracy of the element with the minimum number of stress modes is close to another multi‐field element which, however, is not frame invariant and exhibits unsuppressed zero‐energy deformation modes. Copyright © 2000 John Wiley & Sons, Ltd.     

Analytic preconditioners for the electric field integral equation

Since the advent of the fast multipole method, large‐scale electromagnetic scattering problems based on the electric field integral equation (EFIE) formulation are generally solved by a Krylov iterative solver. A well‐known fact is that the dense complex non‐hermitian linear system associated to the EFIE becomes ill‐conditioned especially in the high‐frequency regime. As a consequence, this slows down the convergence rate of Krylov subspace iterative solvers. In this work, a new analytic preconditioner based on the combination of a finite element method with a local absorbing boundary condition is proposed to improve the convergence of the iterative solver for an open boundary. Some numerical tests precise the behaviour of the new preconditioner. Moreover, comparisons are performed with the analytic preconditioner based on the Calderòn's relations for integral equations for several kinds of scatterers. Copyright © 2004 John Wiley & Sons, Ltd.     

The enriched space–time finite element method (EST) for simultaneous solution of fluid–structure interaction

The paper introduces a weighted residual‐based approach for the numerical investigation of the interaction of fluid flow and thin flexible structures. The presented method enables one to treat strongly coupled systems involving large structural motion and deformation of multiple‐flow‐immersed solid objects. The fluid flow is described by the incompressible Navier–Stokes equations. The current configuration of the thin structure of linear elastic material with non‐linear kinematics is mapped to the flow using the zero iso‐contour of an updated level set function. The formulation of fluid, structure and coupling conditions uniformly uses velocities as unknowns. The integration of the weak form is performed on a space–time finite element discretization of the domain. Interfacial constraints of the multi‐field problem are ensured by distributed Lagrange multipliers. The proposed formulation and discretization techniques lead to a monolithic algebraic system, well suited for strongly coupled fluid–structure systems. Embedding a thin structure into a flow results in non‐smooth fields for the fluid. Based on the concept of the extended finite element method, the space–time approximations of fluid pressure and velocity are properly enriched to capture weakly and strongly discontinuous solutions. This leads to the present enriched space–time (EST) method. Numerical examples of fluid–structure interaction show the eligibility of the developed numerical approach in order to describe the behavior of such coupled systems. The test cases demonstrate the application of the proposed technique to problems where mesh moving strategies often fail. Copyright © 2007 John Wiley & Sons, Ltd.     

A practical large‐strain solid finite element for sheet forming

An alternative approach for developing practical large‐strain finite elements has been introduced and used to create a three‐dimensional solid element that exhibits no locking or hourglassing, but which is more easily and reliably derived and implemented than typical reduced‐integration schemes with hourglassing control. Typical large‐strain elements for forming applications rely on reduced integration to remove locking modes that occur with the coarse meshes that are necessary for practical use. This procedure introduces spurious zero‐energy deformation modes that lead to hourglassing, which in turn is controlled by complex implementations that involve lengthy derivations, knowledge of the material model, and/or undetermined parameters. Thus, for a new material or new computer program, implementation of such elements is a daunting task. Wang–Wagoner‐3‐dimensions (WW3D), a mixed, hexahedral, three‐dimensional solid element, was derived from the standard linear brick element by ignoring the strain components corresponding to locking modes while maintaining full integration (8 Gauss points). Thus, WW3D is easily implemented for any material law, with little chance of programming error, starting from programming for a readily available linear brick element. Surprisingly, this approach and resulting element perform similarly or better than standard solid elements in a series of numerical tests appearing in the literature. The element was also tested successfully for an applied sheet‐forming analysis problem. Many variations on the scheme are also possible for deriving special‐purpose elements. Copyright © 2005 John Wiley & Sons, Ltd.     

Deviatoric hybrid model and multivariable elimination at element level for incompressible medium

A deviatoric hybrid element approach, in which the deviatoric stress σ′, the pressure p and the displacement u are independently dealt with as the element variables, is suggested. The present approach is naturally universal for compressible and fully incompressible mediums. Moreover, it can be extended to the simulation of Stokes flow directly. The resulting hybrid model is able to meet the zero volumetric strain constraint in terms of the incompatible displacement mode only. Therefore an incompressible elimination can be carried out within an individual element, and the complex system elimination for nodal displacements is then avoided. The present 3‐field hybrid model maintains the important features of current hybrid stress elements—finally resulting in a set of displacement‐type discrete equations which can be easily solved, while not a set of u ‐p mixed‐type equations resulted. Regarding the numerical stability of the element, an effective strategy is offered to suppress all the zero energy modes hidden in the model. Copyright © 1999 John Wiley & Sons, Ltd.     

A new enhanced assumed strain quadrilateral membrane element with drilling degree of freedom and modified shape functions

A new four‐node quadrilateral membrane finite element with drilling rotational degree of freedom based on the enhanced assumed strain formulation is presented. A simple formulation is achieved by five incompatible modes that are added to the Allman‐type interpolation. Furthermore, modified shape functions are used to improve the behaviour of distorted elements. Numerical results show that the proposed new element exhibits good numerical accuracy and improved performance, and in many cases, superior to existing elements. In particular, Poisson's locking in nearly incompressible elasticity fades and the element performs well when it becomes considerably distorted even when it takes almost triangular shape. Copyright © 2016 John Wiley & Sons, Ltd.