Hi

I am trying to look at material instabilities in Lee Fenves non-associated plasticity model (aka Concrete damaged plasticity in abaqus) in form of $\det(E_{t})=0$

looking at element stiffness matrix, I am suspicious that the values provided by abaqus is really  a tangential stiffness matrix because:

  • I expect a lower rank stiffness matrix in peak regime
  • In softening, the stiffness matrix shouldn't be positive-definite
  • the stiffness matrix shouldn't be symmetric
  • what I get from abaqus is in contradiction with my expectations. I get a symmetric positive definite matrix with the highest rank possible (18 for case of C3D8 element) which is not the case that i am looking for. That is why i am trying to get DDSDDE matrix from abaqus, which is in constitutive level, but I don't have idea about extracting it as an output. I would appreciate it if anybody can help.

    PS: I attached my abaqus input file

    Similar questions and discussions