用MATLAB进行平面杆系模型的弹性时程分析时遇到的问题

glwh

我自己用MATLAB编的框架结构的平面杆系模型，用的是一致质量矩阵。然后按照徐赵东的《MATLAB在建筑抗震工程中的应用》一书中的弹性时程分析程序(wilson-θ法)来对结构进行计算，结果随时间推移，节点响应越来越大，最后出现无穷大响应。我又把模型换成层剪切模型用同样的弹性时程分析程序进行计算，结果却正常。