运动学的硬化 fortran 代码

代码说明：

此代码为各向同性材料在平面应力状态的金属一样的随动强化你可以使用更精细的随动强化模型，所以用 fortran 语言写一个子程序来实现他们的本构关系。