Iím facing a problem in the implementation of a VUMAT for an hyperelastic anisotropic material, similar to the one implemented as *ANISOTROPIC HYPERELASTIC, HOLZAPFEL. In particular, Iíve noticed that varying the local material orientations using the *ORIENTATION keyword Iím obtaining incoherent results, and sometimes the error "Bad Material definition in element number 1: zero or negative initial dilatational modulus caused by bad material data".
It looks like that when ABAQUS is running the initial data check phase imposing a set of fictitious strains, these strains are applied in the local coordinate system, while the stretchNew variables passed into the subroutine are in the global coordinate system. So, it could happen to obtain a dilatation of the material (J=det(F)>1) even if the strains are compressive.
Iím working on two input files:
- the first one has a rectangular orientation in which the local x,y,z axes have the same direction of the global ones; this one is working fine.
- the second one has a rectangular orientation in which the local the y-axis and the z-axis have the opposite direction of the global ones; this one is giving the Bad Material definition error.
Is this an issue any of you has ever faced?
Thanks for your help,