Finite element formulation for modelling nonlinear viscoelastic elastomers

Nonlinear viscoelastic response of reinforced elastomers is modeled using a three-dimensional mixed finite element method with a nonlocal pressure field. A general second-order unconditionally stable exponential integrator based on a diagonal Padé approximation is developed and the Bergström–Boyce nonlinear viscoelastic law is employed as a prototype model. An implicit finite element scheme with consistent linearization is used and the novel integrator is successfully implemented. Finally, several viscoelastic examples, including a study of the unit cell for a solid propellant, are solved to demonstrate the computational algorithm and relevant underlying physics

Arbitrary order 2D virtual elements for polygonal meshes: Part II, inelastic problem

The present paper is the second part of a twofold work, whose first part is reported in [3], concerning a newly developed Virtual Element Method (VEM) for 2D continuum problems. The first part of the work proposed a study for linear elastic problem. The aim of this part is to explore the features of the VEM formulation when material nonlinearity is considered, showing that the accuracy and easiness of implementation discovered in the analysis inherent to the first part of the work are still retained. Three different nonlinear constitutive laws are considered in the VEM formulation. In particular, the generalized viscoplastic model, the classical Mises plasticity with isotropic/kinematic hardening and a shape memory alloy (SMA) constitutive law are implemented. The versatility with respect to all the considered nonlinear material constitutive laws is demonstrated through several numerical examples, also remarking that the proposed 2D VEM formulation can be straightforwardly implemented as in a standard nonlinear structural finite element method (FEM) framework

Mixed finite element methods for linear elasticity with weakly imposed symmetry

In this paper, we construct new finite element methods for the approximation of the equations of linear elasticity in three space dimensions that produce direct approximations to both stresses and displacements. The methods are based on a modified form of the Hellinger--Reissner variational principle that only weakly imposes the symmetry condition on the stresses. Although this approach has been previously used by a number of authors, a key new ingredient here is a constructive derivation of the elasticity complex starting from the de Rham complex. By mimicking this construction in the discrete case, we derive new mixed finite elements for elasticity in a systematic manner from known discretizations of the de Rham complex. These elements appear to be simpler than the ones previously derived. For example, we construct stable discretizations which use only piecewise linear elements to approximate the stress field and piecewise constant functions to approximate the displacement field.Comment: to appear in Mathematics of Computatio

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