On a fully three-dimensional finite-strain viscoelastic damage model: Formulation and computational aspects

A fully three-dimensional finite-strain viscoelastic model is developed, characterized by: (i) general anisotropic response, (ii) uncoupled bulk and deviatoric response over any range of deformations, (iii) general relaxation functions, and (iv) recovery of finite elasticity for very fast or very slow processes; in particular, classical models of rubber elasticity (e.g. Mooney-Rivlin). Continuum damage mechanics is employed to develop a simple isotropic damage mechanism, which incorporates softening behavior under deformation, and leads to progressive degradation of the storage modulus in a cyclic test with increasing amplitude (Mullins' effect). A numerical integration procedure is proposed which trivially satisfies objectivity and bypasses the use of midpoint configurations. The resulting algorithm can be exactly linearized in closed form, and leads to symmetric tangent moduli. Quasi-incompressible response is accounted for within the context of a three-field variational formulation of the Hu-Washizu type.