求助，单自由度弹塑性时程分析程序

aken

本人编了一个单自由度的弹塑性时程分析程序，但是在程序运行后出来的骨架曲线与我图中的骨架曲线差别非常大，程序运行的骨架曲线只有弹性段，没有塑性段，而且位移、速度、加速度结果感觉也不对，请大家帮忙指点一下我的错误在哪里。我的email:
xiaojingtan@yahoo.com.cn，请与我联系。