I have been trying to simulate a triaxial test of rockfill from a research paper for a while now. I have made the parameters of the duncan chang model i.e. the young's modulus, bulk modulus and the ultimate deviatoric stress functions of stress so that i can use the values of k, k_ur, n, m, phi, Rf from the paper but i have not been able to get it to converge. If anyone has any experience or advice regarding using duncan chang model with comsol I would greatly appreciate your advice. I am uploading my model with this question.