Newmark b method for bi-linear SDOF structure with iteration：用Newmark b法求解带迭代的双线性SDOF结构。 在文本中重现示例表。-matlab开发

最初的程序是由 Christopher Wong 先生开发的。 感谢他的紧凑工作。 我从原始程序中学到了很多东西。 在我的案例中，参考书是 Chopra 博士的《结构动力学》，2001 年第二版，第 191 页。 第 2 版和第 4 版之间的模型参数似乎有所不同。 我试图重现表 E5.6，p192。 原始程序中的收敛过程不清楚。 我是Matlab的初学者。 于是，为了解决非线性DSOF和练习Matlab，我开始修改Christopher Wong先生的原始程序。 以下是修改项目。 1) 使用的参数来自 Chopra 博士的《结构动力学》第二版。 2) 使用全波长正弦波产生双向屈服。 3) 考虑兼容性条件来计算加速度。 4) 引入变量flg 来指定力-变形关系的状态，例如flg=1 弹性，flg=2 屈服，flg=3 回弹弹性。 5) 当flg=2的du*du0<0时产生回弹弹性，其中d