I'm performing Cell model analyses which need to maintain the triaxiality and lode parameter constant same as the analyses by C. Tekoglu, with the title "Representative volume element calculations under constant stress triaxiality, Lode parameter, and shear ratio".

Actually, I have written the code but the problem is that the triaxiality is not getting controlled. If anyone has done this work or has knowledge of this field please respond.

I can share the INP file of the model and the code that I have written. The paper attached herewith of which subroutine I'm trying to write.