The discussers ask why a hyperelastic model was employed. There are several reasons why using an hyper-elastic model is preferable. Probably the most important one is that closed cycles do not produce nor dissipate energy [17, 8]; instead, using a hypoelastic model, closed cycles may produce or dissipate energy, which is inconsistent with the definition of an elastic deformation path. As stressed by Borja et al., the use of an hypoelastic law with constant Poisson’s ratio in the context of the Modified Cam Clay may produce that closed-loop stress paths may result in a non-conservative elastic behavior; therefore, at the end of the loading path the elastic deformation might not be fully recovered. Also, hyper-elastic laws underpin the formulation of large-strain elastoplasticity . By using this theory, the constitutive model is inherently frame-indiferent  and algorithmic objective .