结构动力学多自由度线性体系Wilson-θ法程序编写
多自由度线性体系Wilson-θ法程序编写
【摘要】本文主要介绍了通过使用Matlab软件,Wilson-θ法编写多自由度线性体系的程序的原理、流程图、具体算例以及使用注意事项。通过该程序可以得到剪切型结构在任意函数荷载作用下各质点的位移函数。
【关键词】Matlab;多自由度;Wilson-θ法
1.wilson-法原理
wilson-法中最主要的步骤就是推导由t时刻的状态求时刻的状态的递推公式,现推导如下:
对积分
解出
代入
整理,得
其中
本程序的核心就是对以上公式的循环使用。
2.程序流程图
3.具体应用算例
如图所示,两自由度框架结构,其中
初始静止,求各层位移。
将相应的数据输入到程序中,得出各层位移关于时间的图像。图1为第一层,图2为第二层。
将所得数值解与精确解相比较,图中实线为数值解,虚线为精确解。
由两张图,我们可以看出数值解大致是与精确相近的,但是仍然有些许的不同,这可能是算法中仍然有缺陷,说明程序仍然有待改善。
4.程序使用注意事项
(1)本程序针对于剪切型刚架结构,对于其他结构无法使用。
(2)本程序中各质点的荷载必须是函数的形式(包括常数),即对于只有某些点的荷载无法使用,且荷载函数输入时,必须采用inline语句。例如荷载为常数10,则输入inline(’10’);
如荷载函数为sin(at),则输入inline(’sin(a*t)’,’t’)。
(3)本程序主要针对无阻尼情况,若有阻尼,只需输入阻尼矩阵即可。
(4)θ的值应大于1.37,通常取1.4,优化值为1.420815。
(5)从第一层开始为m1,m2……
【参考文献】
[1] 王焕定. 结构力学(第3版)[M]. 北京出版社.
[2] Anil K.Chopra. 结构动力学理论及其在地震工程中的应用(第2版)[M]. 北京出版社.
多自由度线性体系Wilson-θ法程序编写
2
图1
图2