Discontinuous Galerkin methods for non‐linear elasticity

This paper presents the formulation and a partial analysis of a class of discontinuous Galerkin methods for quasistatic non‐linear elasticity problems. These methods are endowed with several salient features. The equations that define the numerical scheme are the Euler–Lagrange equations of a one‐field variational principle, a trait that provides an elegant and simple derivation of the method. In consonance with general discontinuous Galerkin formulations, it is possible within this framework to choose different numerical fluxes. Numerical evidence suggests the absence of locking at near‐incompressible conditions in the finite deformations regime when piecewise linear elements are adopted. Finally, a conceivable surprising characteristic is that, as demonstrated with numerical examples, these methods provide a given accuracy level for a comparable, and often lower, computational cost than conforming formulations. Stabilization is occasionally needed for discontinuous Galerkin methods in linear elliptic problems. In this paper we propose a sufficient condition for the stability of each linearized non‐linear elastic problem that naturally includes material and geometric parameters; the latter needed to account for buckling. We then prove that when a similar condition is satisfied by the discrete problem, the method provides stable linearized deformed configurations upon the addition of a standard stabilization term. We conclude by discussing the complexity of the implementation, and propose a computationally efficient approach that avoids looping over both elements and element faces. Several numerical examples are then presented in two and three dimensions that illustrate the performance of a selected discontinuous Galerkin method within the class. Copyright © 2006 John Wiley & Sons, Ltd.     

A locking‐free selective smoothed finite element method using tetrahedral and triangular elements with adaptive mesh rezoning for large deformation problems

A novel finite element (FE) formulation with adaptive mesh rezoning for large deformation problems is proposed. The proposed method takes the advantage of the selective smoothed FE method (S‐FEM), which has been recently developed as a locking‐free FE formulation with strain smoothing technique. We adopt the selective face‐based smoothed/node‐based smoothed FEM (FS/NS‐FEM‐T4) and edge‐based smoothed/node‐based smoothed FEM (ES/NS‐FEM‐T3) basically but modify them partly so that our method can handle any kind of material constitutive models other than elastic models. We also present an adaptive mesh rezoning method specialized for our S‐FEM formulation with material constitutive models in total form. Because of the modification of the selective S‐FEMs and specialization of adaptive mesh rezoning, our method is locking‐free for severely large deformation problems even with the use of tetrahedral and triangular meshes. The formulation details for static implicit analysis and several examples of analysis of the proposed method are presented in this paper to demonstrate its efficiency. Copyright © 2014 John Wiley & Sons, Ltd.     

Locking in the incompressible limit for the element‐free Galerkin method

Volumetric locking (locking in the incompressible limit) for linear elastic isotropic materials is studied in the context of the element‐free Galerkin method. The modal analysis developed here shows that the number of non‐physical locking modes is independent of the dilation parameter (support of the interpolation functions). Thus increasing the dilation parameter does not suppress locking. Nevertheless, an increase in the dilation parameter does reduce the energy associated with the non‐physical locking modes; thus, in part, it alleviates the locking phenomena. This is shown for linear and quadratic orders of consistency. Moreover, the biquadratic order of consistency, as in finite elements, improves the locking behaviour. Although more locking modes are present in the element‐free Galerkin method with quadratic consistency than with standard biquadratic finite elements. Finally, numerical examples are shown to validate the modal analysis. In particular, the conclusions of the modal analysis are also confirmed in an elastoplastic example. Copyright © 2001 John Wiley & Sons, Ltd.     

Adaptive methods for thermomechanical coupled contact problems

In this paper an adaptive method for the analysis of thermomechanical coupled multi‐body contact problems is presented. The method is applied to non‐linear elastic solids undergoing finite (thermal) deformations. The contact model considers non‐linear pressure‐dependent heat flux as well as frictional heating in the interface. A time–space‐finite element discretization of the governing equations is formulated including unilateral constraints due to contact. A staggered solution algorithm has been constructed that allows an independent spatial discretization of the coupled subproblems. A posteriori projection‐based error estimators, which enforce implicitly the special boundary conditions due to thermal contact, are used to control the spatial discretization as well as the adaptive time stepping. Numerical examples are presented to corroborate the applicability of the adaptive algorithm to the considered problem type. Copyright © 2004 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.     

A finite deformation brick element with inhomogeneous mode enhancement

In this paper we describe a new enhanced assumed strain finite element for finite deformations. The element is based on the split of the deformation of an element into a homogeneous and inhomogeneous part. The enhancement is applied to the inhomogeneous part only. For the homogeneous part a compressible Neo‐Hooke material is used, while for the inhomogeneous part linear elasticity is assumed. In several examples it is shown that the element is locking and hourglassing free as well as insensitive to initial element distortion. Copyright © 2008 John Wiley & Sons, Ltd.     

An Eulerian finite‐volume scheme for large elastoplastic deformations in solids

Conservative formulations of the governing laws of elastoplastic solid media have distinct advantages when solved using high‐order shock capturing methods for simulating processes involving large deformations and shock waves. In this paper one such model is considered where inelastic deformations are accounted for via conservation laws for elastic strain with relaxation source terms. Plastic deformations are governed by the relaxation time of tangential stresses. Compared with alternative Eulerian conservative models, the governing system consists of fewer equations overall. A numerical scheme for the inhomogeneous system is proposed based upon the temporal splitting. In this way the reduced system of non‐linear elasticity is solved explicitly, with convective fluxes evaluated using high‐order approximations of Riemann problems locally throughout the computational mesh. Numerical stiffness of the relaxation terms at high strain rates is avoided by utilizing certain properties of the governing model and performing an implicit update. The methods are demonstrated using test cases involving large deformations and high strain rates in one‐, two‐, and three‐dimensions. Copyright © 2009 John Wiley & Sons, Ltd.     

Extension of the ‘solid‐shell’ concept for application to large elastic and large elastoplastic deformations

In the present contribution we extend a previously proposed so‐called solid–shell concept which incorporates only displacement degrees of freedom to the simulation of large elastic and large elastoplastic deformations of shells. Therefore, the modifications necessary for hyper‐elastic or elastoplastic material laws are discussed. These modifications concern the right Cauchy–Green tensor for large elastic deformations, respectively, the deformation gradient for elastoplasticity which then are consistent to the modified Green–Lagrange strains that are necessary for transverse shear and membrane locking free solid–shell element formulations. However, in addition to the locking mentioned above especially in the range of plasticity incompressibility locking becomes important. Thus, the second major aspect of this contribution is the discussion of several ways to avoid incompressibility locking also including the investigation of eigenmodes. Finally, a selective reduced integration scheme with reduced integration for the volumetric term is employed and described in detail, although it is limited to material laws which allow the decomposition into a volumetric and a deviatoric part. Some numerical examples show the range of application for the proposed elements. Copyright © 2000 John Wiley & Sons, Ltd.     

Non‐linear analysis of structures using high performance hybrid elements

In this work, we describe the formulation and implementation for stress‐based hybrid elements for conducting non‐linear analysis of elastic structures. The motivation behind developing these elements is that they should be as simple to use as standard displacement‐based isoparametric brick elements, but at the same time, be relatively immune to the shortcoming that these elements suffer from, namely, ‘locking’ problems which occur when they are used to model plate/shell geometries, almost incompressible materials or when the elements are distorted, and so on. The formulation is based on a two‐field mixed variational principle. Numerical examples are presented to demonstrate the excellent performance of the proposed elements on a variety of challenging problems involving very large deformations, buckling, mesh distortions, almost incompressible materials, etc. Copyright © 2006 John Wiley & Sons, Ltd.     

An enriched meshless method for non‐linear fracture mechanics

This paper presents an enriched meshless method for fracture analysis of cracks in homogeneous, isotropic, non‐linear‐elastic, two‐dimensional solids, subject to mode‐I loading conditions. The method involves an element‐free Galerkin formulation and two new enriched basis functions (Types I and II) to capture the Hutchinson–Rice–Rosengren singularity field in non‐linear fracture mechanics. The Type I enriched basis function can be viewed as a generalized enriched basis function, which degenerates to the linear‐elastic basis function when the material hardening exponent is unity. The Type II enriched basis function entails further improvements of the Type I basis function by adding trigonometric functions. Four numerical examples are presented to illustrate the proposed method. The boundary layer analysis indicates that the crack‐tip field predicted by using the proposed basis functions matches with the theoretical solution very well in the whole region considered, whether for the near‐tip asymptotic field or for the far‐tip elastic field. Numerical analyses of standard fracture specimens by the proposed meshless method also yield accurate estimates of the J‐integral for the applied load intensities and material properties considered. Also, the crack‐mouth opening displacement evaluated by the proposed meshless method is in good agreement with finite element results. Furthermore, the meshless results show excellent agreement with the experimental measurements, indicating that the new basis functions are also capable of capturing elastic–plastic deformations at a stress concentration effectively. Copyright © 2003 John Wiley & Sons, Ltd.     

The enhanced assumed strain method for the isogeometric analysis of nearly incompressible deformation of solids

Isogeometric analysis has recently become very popular for the numerical modeling of structures and fluids. Among other potential features, advantages of using a non‐uniform rational B‐splines (NURBS)‐based isogeometric analysis over the traditional finite element method include the possibility of using higher‐order polynomials for the basis functions of the approximation space, which may be easily built on a recursive (hierarchical) fashion as well as higher convergence ratio. Nevertheless, NURBS‐based isogeometric analysis suffers from the same problems depicted by other methods when it comes to reproduce isochoric deformations, that is, it shows volumetric locking, especially for low‐order basis functions. Similar remedies as those that have been proposed for the finite element method may be appropriate for integration in the NURBS‐based isogeometric analysis and some have already been tried with success. In this work, the analysis of the underlying space of incompressible deformations of a NURBS‐based isogeometric approximation is performed with the main objective of understanding the likelihood of volumetric locking. As a remedy, the enhanced assumed strain methodology is blended with the NURBS‐based isogeometric analysis to alleviate the volumetric locking associated with incompressible deformations. The solution includes a stabilization term derived directly from a penalized form of the classical Veubeke–Hu–Washizu three‐field variational principle. Copyright © 2012 John Wiley & Sons, Ltd.     

On the design of energy–momentum integration schemes for arbitrary continuum formulations. Applications to classical and chaotic motion of shells

The construction of energy–momentum methods depends heavily on three kinds of non‐linearities: (1) the geometric (non‐linearity of the strain–displacement relation), (2) the material (non‐linearity of the elastic constitutive law), and (3) the one exhibited in displacement‐dependent loading. In previous works, the authors have developed a general method which is valid for any kind of geometric non‐linearity. In this paper, we extend the method and combine it with a treatment of material non‐linearity as well as that exhibited in force terms. In addition, the dynamical formulation is presented in a general finite element framework where enhanced strains are incorporated as well. The non‐linearity of the constitutive law necessitates a new treatment of the enhanced strains in order to retain the energy conservation property. Use is made of the logarithmic strain tensor which allows for a highly non‐linear material law, while preserving the advantage of considering non‐linear vibrations of classical metallic structures. Various examples and applications to classical and non‐classical vibrations and non‐linear motion of shells are presented, including (1) chaotic motion of arches, cylinders and caps using a linear constitutive law and (2) large overall motion and non‐linear vibration of shells using non‐linear constitutive law. Copyright © 2004 John Wiley & Sons, Ltd.     

A parallel computational method in steady power‐law creep

A new approach to parallelization of materially non‐linear problems in solid mechanics is developed. It is based on approximating generalized models of subdomains. The procedure does not retain the same substructuring technique used in a linear version. The convergence proof of the single‐ and multilevel‐domain decomposition algorithms uses the principle of minimum potential energy dissipation and investigated properties of the substructural models. The high efficiency of the approach introduced is shown through the study of several examples. The method developed in this paper for steady creep can be used without modification to solve non‐linear elasticity problems and, at active loading, plasticity problems for bodies of the power‐law strain–stress diagrams. Copyright © 2001 John Wiley & Sons, Ltd.     

Efficient simulation of multi‐body contact problems on complex geometries: A flexible decomposition approach using constrained minimization

We consider the numerical simulation of non‐linear multi‐body contact problems in elasticity on complex three‐dimensional geometries. In the case of warped contact boundaries and non‐matching finite element meshes, particular emphasis has to be put on the discretization of the transmission of forces and the non‐penetration conditions at the contact interface. We enforce the discrete contact constraints by means of a non‐conforming domain decomposition method, which allows for optimal error estimates. Here, we develop an efficient method to assemble the discrete coupling operator by computing the triangulated intersection of opposite element faces in a locally adjusted projection plane but carrying out the required quadrature on the faces directly. Our new element‐based algorithm does not use any boundary parameterizations and is also suitable for isoparametric elements. The emerging non‐linear system is solved by a monotone multigrid method of optimal complexity. Several numerical examples in 3D illustrate the effectiveness of our approach. Copyright © 2008 John Wiley & Sons, Ltd.     

Polygonal finite elements for finite elasticity

Nonlinear elastic materials are of great engineering interest, but challenging to model with standard finite elements. The challenges arise because nonlinear elastic materials are characterized by non‐convex stored‐energy functions as a result of their ability to undergo large reversible deformations, are incompressible or nearly incompressible, and often times possess complex microstructures. In this work, we propose and explore an alternative approach to model finite elasticity problems in two dimensions by using polygonal discretizations. We present both lower order displacement‐based and mixed polygonal finite element approximations, the latter of which consist of a piecewise constant pressure field and a linearly‐complete displacement field at the element level. Through numerical studies, the mixed polygonal finite elements are shown to be stable and convergent. For demonstration purposes, we deploy the proposed polygonal discretization to study the nonlinear elastic response of rubber filled with random and periodic distributions of rigid particles, as well as the development of cavitation instabilities in elastomers containing vacuous defects. These physically‐based examples illustrate the potential of polygonal finite elements in studying and modeling nonlinear elastic materials with complex microstructures under finite deformations. Copyright © 2014 John Wiley & Sons, Ltd.     

Smooth finite strain plasticity with non‐local pressure support

The aim of this work is to introduce an alternative framework to solve problems of finite strain elastoplasticity including anisotropy and kinematic hardening coupled with any isotropic hyperelastic law. After deriving the constitutive equations and inequalities without any of the customary simplifications, we arrive at a new general elasto‐plastic system. We integrate the elasto‐plastic algebraico‐differential system and replace the loading–unloading condition by a Chen–Mangasarian smooth function to obtain a non‐linear system solved by a trust region method. Despite being non‐standard, this approach is advantageous, since quadratic convergence is always obtained by the non‐linear solver and very large steps can be used with negligible effect in the results. Discretized equilibrium is, in contrast with traditional approaches, smooth and well behaved. In addition, since no return mapping algorithm is used, there is no need to use a predictor. The work follows our previous studies of element technology and highly non‐linear visco‐elasticity. From a general framework, with exact linearization, systematic particularization is made to prototype constitutive models shown as examples. Our element with non‐local pressure support is used. Examples illustrating the generality of the method are presented with excellent results. Copyright © 2009 John Wiley & Sons, Ltd.     

An adjoint variable method for sensitivity analysis of non-linear elastic systems

An adjoint variable method for design sensitivity analysis of non-linear elastic systems is presented. The method uses domain parameterization and a mutual form of the Hu-Washizu energy principle, and extends results reported in a recent work for linear elastic systems to non-linear elasticity. Non-linearities due to finite deformations and non-linear, hyperelastic constitutive models are considered. In contrast to other methods for non-linear sensitivity analysis, the present formulation can be applied with force, displacement or mixed approximate solution methods. The mutual energy expression used in the adjoint sensitivity derivation is developed from a non-linear extension of the Hu-Washizu energy functional and yields a linear governing equation for the adjoint system. This has important ramifications for the computational cost of a sensitivity analyses of non-linear systems: excluding the cost of determining the response of the system, the cost of a sensitivity analysis for a non-linear system is essentially the same as that for a linear system. Finite element implementation of the resulting sensitivity expressions is discussed, and two numerical examples are presented. The first example involves large deformations of a Mooney-Rivlin body, while the second involves design sensitivity analysis for mixed solution methods.     

Efficient and accurate multilayer solid‐shell element: non‐linear materials at finite strain

We present in this paper an efficient and accurate low‐order solid‐shell element formulation for analyses of large deformable multilayer shell structures with non‐linear materials. The element has only displacement degrees of freedom (dofs), and an optimal number of enhancing assumed strain (EAS) parameters to pass the patch tests (both membrane and out‐of‐plane bending) and to remedy volumetric locking. Based on the mixed Fraeijs de Veubeke‐Hu‐Washizu (FHW) variational principle, the in‐plane and out‐of‐plane bending behaviours are improved and the locking associated with (nearly) incompressible materials is avoided via a new efficient enhancement of strain tensor. Shear locking and curvature thickness locking are resolved effectively by using the assumed natural strain (ANS) method. Two non‐linear 3‐D constitutive models (Mooney–Rivlin material and hyperelastoplastic material at finite strain) are applied directly without requiring the enforcement of the plane‐stress assumption. In particular, we give a simple derivation for the hyperelastoplastic model using spectral representations. In addition, the present element has a well‐defined lumped mass matrix, and provides double‐side contact surfaces for shell contact problems. With the dynamics referred to a fixed inertial frame, the present element can be used to analyse multilayer shell structures undergoing large overall motion. Numerical examples involving static analyses and implicit/explicit dynamic analyses of multilayer shell structures with both material and geometric non‐linearities are presented, and compared with existing results obtained from other shell elements and from a meshless method. It is shown that elements that did not pass the out‐of‐plane bending patch test could not provide accurate results, as compared to the present element formulation, which passed the out‐of‐plane bending patch test. The present element proves to be versatile and efficient in the modelling and analyses of general non‐linear composite multilayer shell structures. Copyright © 2005 John Wiley & Sons, Ltd.