A generic anisotropic continuum damage model integration scheme adaptable to both ductile damage and biological damage-like situations

• We present an integration for anisotropic continuum damage in finite strains.
• It is a staggered scheme with an analytical consistent tangent operator.
• Its main feature is its versatility, it can be coupled to any damage evolution.
• The model is particularised to represent biologically induced stiffness variation.
• It is thus a fully coupled remodelling model applied in orthodontic tooth movement.

This paper aims at presenting a general versatile time integration scheme applicable to anisotropic damage coupled to elastoplasticity, considering any damage rate and isotropic hardening formulations. For this purpose a staggered time integration scheme in a finite strain framework is presented, together with an analytical consistent tangent operator. The only restrictive hypothesis is to work with an undamaged isotropic material, assumed here to follow a J2 plasticity model. The only anisotropy considered is thus a damage-induced anisotropy. The possibility to couple any damage rate law with the present algorithm is illustrated with a classical ductile damage model for aluminium, and a biological damage-like application. The later proposes an original bone remodelling law coupled to trabecular bone plasticity for the simulation of orthodontic tooth movements. All the developments have been considered in the framework of the implicit non-linear finite element code Metafor (developed at the LTAS/MN2L, University of Liège, Belgium – www.metafor.ltas.ulg.ac.be).