原计划是由克里斯托弗·黄先生。谢谢他的紧凑的工作。我学到了很多从原来的计划。我的参考书是结构动力学Chopra博士第二ed.2001 p191。模型参数似乎不同的第二和第四版之间。我试图复制表E5.6 p192。原始程序的收敛过程还不清楚。我是一个初学者的Matlab。所以,解决非线性DSOF和练习Matlab的缘故,我开始修改原程序的Christopher Wong先生。以下是修改的项目。

1)第二版的参数使用的动态结构Chopra博士。
2)全波长度的正弦波是用于生成产生在两个方向上。
3)相容性条件被认为是计算加速度。
4)flg被介绍给一个变量指定force-deformation关系状态,如flg = 1弹性,flg = 2产量、flg = 3弹性反弹。
5)反弹弹性时生成du * du0 < 0 flg = 2,杜是增量变形。