第5行：对输出界面清空，对计算保留的数据清空。

第7~10行：设置计算所需要的动力参数。

第11~12行：读取txt文本的地震波。（注：也可以其他方式读取）

第14~15行：黏滞阻尼器所需要的计算参数。

第17~19行：设置双线性滞回的关键参数。

第20行：计算系统阻尼比。

第22~23行：计算前处理所需要整理的变量。

第25行：基于地震波步进迭代开始。

第26~34行：双线性滞回的未屈服阶段判断进行确定刚度。

第35~41行：双线性滞回的正向屈服阶段判断进行确定刚度。

第42~48行：双线性滞回的反向屈服阶段判断进行确定刚度。

第53~66行：进行Newmark-β的计算分析。（详情请阅读任意一本结构动力学）

第67~69行：结果进行整理储存。

第70~72行：结果进行整理，得到了相对加速度Acc，绝对加速度Aabs，位移Dish，恢复力Fsh。

到站了，希望各位看官有所收获！

源代码获得可关注“建源学堂”公众号，后台回复：Matlab双线性教学。