某四层框架结构在顶部受一个简谐荷载
�
�
�
�
�
�
�
�
�
1
0
4
sin
t
t
FF
�
的作用，力的作用时间 t
1
=5.0s，
F
0
=100N。试计算响应的时间为 20s（计算分 2000 步完成）。阻尼矩阵由 Rayleigh 阻尼构
造。
结构基本计算简图
计算 Newmark-beta 法的 matlab 程序
m=[1,2,3,4];
m=diag(m);
k= [800 -800 0 0;
-800 2400 -1600 0;
0 -1600 4800 -3200;
0 0 -3200 8000];
c=0.05*m+0.02*k;
f0=100;
t1=5;
nt=2000;
dt=0.01;
gama=0.5;
beta=0.25;
评论2