Recovering Mullins damage hyperelastic behaviour with physics augmented neural networks

The aim of this work is to develop a neural network for modelling incompressible hyperelastic behaviour with isotropic damage, the so-called Mullins effect. This is obtained through the use of feed-forward neural networks with special attention to the architecture of the network in order to fulfil several physical restrictions such as objectivity, polyconvexity, non-negativity, material symmetry and thermodynamic consistency. The result is a compact neural network with few parameters that is able to reconstruct the hyperelastic behaviour with Mullins-type damage. The network is trained with artificially generated plane stress data and even correctly captures the full 3D behaviour with much more complex loading conditions. The energy and stress responses are correctly captured, as well as the evolution of the damage. The resulting neural network can be seamlessly implemented in widely used simulation software. Implementation details are provided and all numerical examples are performed in Abaqus.