Dear Jorgen,

I have to bother you again, on the problem of Large volume changes with geometric nonlinearity when I writing a UMAT of Drucker Prager constitutive model. My umat can run but the convergence rate is slow, so that some times I have to modify the convergence criterion a little bit. I suppose that the problem is due to the poor Jacobian Matrix when the stresses of an integration point approaches the apex, where I set a tensile strength. Thus there is only adjustment of the hydrastatic stresses. However, I am afraid that the Jacobian Matrix possibly casuses converange problem, if the large defromation exists at the same time.

The manual said that If the material model allows large volume changes and geometric nonlinearity is considered, the exact definition of the consistent Jacobian should be used to ensure rapid convergence. For the rate-form constitutive laws, which is commonly encountered by using pressure-dependent plasticity, the exact consistent Jacobian is given by

C=(1/J) * [Par (delta(J*stress)/Par (delta(strain))]

Based on manual, DFGRD0(3,3) and DFGRD1(3,3) provide deformation gradient at the beginning and the end of the increment. But how to calculate C using the DFGRD0(3,3) and DFGRD1(3,3)? I never read the discussion on this topic.

Thanks and waiting for your inspiration.

Best regards,

Willness