I am trying to predict cutting forces in orthogonal cutting using Oxley's model using power law in flow stresses. However, I can't get out of the loop where temperature is just changing iteratively. In fact, at some moment, temperature increase in shear zone is negative due to negative value of beta. I can't figure out where the problem is.