循环加载问题

小弟最近想模拟拉－压循环加载来研究包辛格效应，不过遇到一些问题还请大虾能够帮忙解答，小弟不胜感激！
我选用的材料具体参数为弹性模量210E9 ，泊松比为0.3 ，密度为7800，同时采用混合硬化模型：C＝25.5E9，gama＝81，初始屈服应力为200E6，Q=2000E6，b＝0.26。
采用explicit进行计算每步长度为0.01秒共分4步（拉－压－拉－压）
拉压两次循环得到如图结果，感觉应力单位的数量级有问题，同时循环加载曲线好像和资料上的也有一些出入，不知是为什么。