作业1
绘制1940 El Centro，N—S分量地震动的绝对加速度、相对速度和相对位移反应谱。地震动:在PEER Ground Motion Databa自行下载经典的1940 El Centro，N-S分量。 要求：在此模板内完成，A4纸打印.自编程序与软件(Bispec或Seismosigna等）计算反应谱进行对比。提交自编写程序。
Matlab程序：
clear
fid = fopen（'E：\Earthquake\El centro。txt’）；
［Accelerate，count] = fscanf（fid，’％g'）； ％count 读入的记录的量
Accelerate=9.8*Accelerate’; %单位统一为 m和s
time=0:0.02：(count-1）*0。02; ％单位 s
Displace=zeros(1，count)； ％相对位移
Velocity=zeros(1,count)； ％相对速度
AbsAcce=zeros（1，count); ％绝对加速度
DampA=[0.00,0。02，0.05]; %三个阻尼比
TA=0.0：0.02：4； ％TA=0.000001:0。02：4;结构周期
Dt=0.02； %地震记录的步长
%记录计算得到的反应，MDis为某阻尼时最大相对位移
％MVel为某阻尼时最大相对速度，MAcc某阻尼时最大绝对加速度
MDis=zeros(3,length(TA))；
MVel=zeros(3，length(TA))；
MAcc=zeros(3,length（TA)）；
j=1；
for Damp=[0.00，0.05，0。1]
t=1；
for T=0。0：0.02：4
Frcy=2＊pi/T ；
DamFrcy=Frcy＊sqrt(1-Damp＊Damp)；
e_t=exp(—Damp*Frcy*Dt)；
s=sin(DamFrcy＊Dt)；
c=cos（DamFrcy*Dt）；
毕业设计英文翻译 A=zeros（2,2）;
A（1,1)=e_t＊（s＊Damp/sqrt(1-Damp＊Damp)+c)；
A（1，2）=e_t*s/DamFrcy；
auspice
A(2,1)=-Frcy*e_t*s/sqrt（1-Damp*Damp)；
A（2,2)=e_t*(—s＊Damp/sqrt（1—Damp＊Damp)+c)；
d_f=（2*Damp^2-1）/（Frcy^2*Dt)；
d_3t=Damp/(Frcy^3*Dt);
B=zeros(2,2)；
B（1,1）=e_t＊((d_f+Damp/Frcy）*s/DamFrcy+（2*d_3t+1/Frcy^2)*c）—2*d_3t；
B（1，2）=—e_t＊（d_f＊s/DamFrcy+2*d_3t＊c)—1/Frcy^2+2*d_3t; B（2，1）=e_t*（(d_f+Damp/Frcy)*（c—Damp/sqrt(1-Damp^2）＊s)-(2＊d_3t+1/Frcy^2）＊
(DamFrcy*s+Damp*Frcy*c)）+1/（Frcy^2＊Dt）;
B(2，2)=e_t*(1/(Frcy^2＊Dt)*c+s*Damp/（Frcy*DamFrcy＊Dt))-1/(Frcy^2＊Dt)；
for i=1：(count—1）
Displace（i+1)=A（1，1）*Displace(i)+A(1,2)＊Velocity(i）+B(1，1)＊Accelerate(i）+B（1，2）＊
Accelerate（i+1）；
Velocity（i+1)=A(2,1）*Displace（i)+A（2，2）*Velocity(i)+B（2,1）*Accelerate(i）+B（2,2)*
Accelerate(i+1);
AbsAcce（i+1)=-2＊Damp＊Frcy＊Velocity(i+1）—Frcy^2＊Displace(i+1）；
end
MDis（j，t)=max（abs(Displace));
MVel（j,t）=max（abs（Velocity))；
if T==0.0
MAcc(j，t)=max（abs(Accelerate）)；
el
MAcc（j，t)=max(abs（AbsAcce））；
end
Displace=zeros(1，count）;
Velocity=zeros(1,count);
AbsAcce=zeros(1，count）；
t=t+1;
end
j=j+1;
end
clo all
figure %绘制位移反应谱
plot（TA，MDis（1,：),'—b’,TA,MDis（2，：),’-r'，TA,MDis（3，:),’：k')
title(’Displacement’)
xlabel（’Tn(s)'）
ylabel（’Displacement（m)’)
legend('ζ=0'，’ζ=0。02','ζ=0。05')
grid
figure %绘制速度反应谱
plot（TA,MVel（1,:)，'-b'，TA，MVel（2,:）,'-r’,TA，MVel(3，:）,’:k’）
深圳公共营养师title（'Velocity'）
xlabel（’Tn(s)')
ylabel('velocity（m/s)’）
legend（’ζ=0',’ζ=0.02'，'ζ=0.05’）
grid
figure ％绘制绝对加速度反应谱
plot（TA，MAcc（1，：)，’-b’,TA，MAcc（2,:）,’—r',TA,MAcc（3，：），':k’)
title('Absolute Acceleration’）
xlabel（'Tn(s)')
ylabel（'absolute acceleration（m/s^2）'）
legend（'ζ=0’,'ζ=0.02'，'ζ=0。05')
grid
由此得到下图反应谱与由Seismosigna得到的反应谱作出比较
basket ca相对位移反应谱（上）与Seismosigna计算相对位移反应谱（下）的比较
相对速度反应谱（上)与Seismosigna计算相对速度反应谱(下）的比较
绝对加速度反应谱（左图)与Seismosigna计算绝对加反应谱（右图）的比较
作业2、
mondial使用中心差分方法和Newmark方法分别计算下面单自由度结构的地震反应.题：考虑下面的具有2个自由度的一个简单系统，分析时间步长Δt=0.28c的情况,求每时刻的加速度、速度、位移
求：
要 求:在此模板内完成，A4纸打印。
日语字典
MATLAB程序：
clear
％求自振圆频率与振型
m=[2 0;0 1];
k=[6 —2；—2 4]；
［v，d]=eig(inv(m)＊k);
w=sqrt（d）；
fai1=v（:，1)。/v(1,1）;
fai2=v(:，2)./v（1，2)；
fai=[fai1 fai2］；
mm=fai’＊m＊fai；kk=fai'*k*fai； %等效质量矩阵及等效刚度矩阵
p0=［0;10]；
％中心差分法
s0=[0；0］;u0=[0;0］； ％初速度及初位移均为0
％以下过程均采用q、q1、q2的形式分别表示位移、速度、加速度
q0=［fai1'＊m＊s0./(fai1'＊m*fai1）;fai2'*m＊s0./(fai2'＊m＊fai2）］;
q01=［fai1’*m＊u0。/（fai1’＊m*fai1)；fai2'＊m＊u0./(fai2’*m*fai2)］；
P0=fai'＊p0;
q02=mm^（—1)*(P0-kk*q0);
t=0.28；
q（:，1）=q0-t*q01+t^2＊q02/2； %q（:，1）表示q-1，即初始时刻前一时刻的位移
q（：,2）=[0；0]； ％q(：，2）表示初始时刻的位移
Kh=mm/（t^2）；
felixstowea=mm/（t^2）;
b=kk—2*mm/(t^2）;
P=fai’＊p0;
for i=2:7;
web proxy rver
Ph（：，i)=P—a*q（:，i—1）—b＊q（：，i)；
q(:,i+1）=Kh^(—1）＊Ph（:，i）；
q1(:,i)=(q(：,i+1）—q(：，i-1))/(2＊t);
q2（:,i）=（q(：，i+1)—2*q（：,i）+q(：,i-1））/（t^2);
end
u=fai*q；
u1=fai＊q1；
u2=fai*q2;
for i=1:6
U(：,i）=u(:，i+1）;
U1（:，i)=u1(:,i+1）;
U2（:，i）=u2（：，i+1）;
end
U
保持距离英文兽欲U1
U2
%Newmark法（采用线加速度法）
gamma=1/2;
beta=1/6;
Kh1=kk+mm/（beta*t^2）；
a1=mm/（beta＊t);
b1=mm/(2＊beta);
detaP=[0 0 0 0 0;0 0 0 0 0］；
q1（：,1)=［0;0］；
q11（:,1）=[0;0]；
q12(:，1)=q02;
for i=1：5
detaPh(：,i）=detaP（:,i）+a1＊q11(：,i）+b1*q12(:，i)；
detaq(:，i）=Kh1^(-1)*detaPh(:,i）；
detaq1（：，i）=gamma*detaq（:，i）/(beta＊t)—gamma＊q11（：，i）/beta+t＊(1—gamma/(2＊beta））*q12(:,i）;
