为了便于进行一些概念性研究，特编写了一个糖葫芦串层模型程序。

程序下载地址

http://www.luxinzheng.net/download/MDOF_Public.f90

程序非常简单，主要部件只有两个：

1 时程积分算法。

这里用了最简单的中心差分法。公式为

虽然中心差分法精度稍差，但是对于层模型来说已经足够了

2 楼层的滞回模型

一共包括四个滞回模型：

(1)    平行四边形弹塑性模型

(2)    Clough模型

(3)    万能塑性铰模型（陆新征－曲哲10参数塑性铰模型）

注意，前两个模型只有硬化、没有软化，所以是Happy Model，不适用于研究倒塌问题。万能塑性铰模型有软化模块，可以用于进行倒塌模拟。以一6层结构为例，计算得到的总位移和层间位移如下

供大家学习，参考。程序刚写完，可能还有一些错误，请大家及时指正。

程序下载地址

http://www.luxinzheng.net/download/MDOF_Public.f90