|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
一个弹性时程分析程序设计例题:某三层钢筋混凝土结构地震波采用200gal el centre波,采样周期为0.02s.用弹性时程分析发求解结构地震反应.程序如下:
clear;
h=[4000,3300,3300]
m=[2.762 2.760 2.300]*1e+3;
k0=[2.485 1.921 1.522]*1e+5;
cn=length(m) %总层数
%earthquake parameter
ct=1.4 %Wilson-法的值
dt=0.02
load ('dzhbo1')
xs=200/max(abs(dzhbo1)) %调整地震输入加速度幅值200gal
ag=dzhbo1*0.01*xs; %地震波,单位m/s2
ndzh=400
ag1=ag(1:ndzh)
ag2=ag(2:ndzh+1)
agtao=ct*(ag2-ag1) %initial value
chsh=zero(cn,1)
wyi1=chsh
sdu1=chsh
jsdu1=chsh
wyimt=chsh
sdumt=chsh
jsdumt=chsh
unit=one(cn,1)
m=diag(m)
[ik]=matrixju(k0,cn)
[x,d]=eig(ik,m)
d=sqrt(d)
w=sort(diag(d))
a=2*w(1)*w(2)*(0.05*w(2)-0.07*w(1))/(w(2)^2-w(1)^2);
b=2*(0.07*w(2)-0.05*w(1))/(w(2)^2-w(1)^2);
c0=a*m+b*ik
运行后说红字部分错误?那个地震波不是个.m文件吗?是不是就是地震的数据呢? |
评分
-
1
查看全部评分
-
|