Radial Return Mapping Algorithm for Plasticity

Hi guys, I am currently working on building my own model of radial return mapping algorithm for plasticity of a beam, but encountered some issue (which I don’t know why, the values just doesn’t match with the one I got from commercial software like Abaqus or Comsol)…… If any of you are willing to enlighten me / share about how your ‘successful code’ looks like, please do so……..

For helping you, we need a description (in terms of equations) of the problem you want to solve and/or a code try.