I tried to write VUMAT for implementation of Moooney-Rivlin hyperelastic. My VUMAT calculate stresses correctly (I check it with MATLAB). I model a cube that pull with negative pressure boundary condition in one direction. It works correctly for one element but for more than 8 element my model distort excessively.
it's the same for displacement boundary condition.
Can you help me please??