fid = fopen('el.txt');
[Accelerate,count] = fscanf(fid,'%g'); %count 读入的记录的量
time=0:0.02:(count-1)*0.02;
% ***********线性加速度法计算各反应***********
%初始化各储存向量
Displace=zeros(1,count); %相对位移
Velocity=zeros(1,count); %相对速度
大学班长职责
损伤容限AbsAcce=zeros(1,count);%绝对加速度
烤箱烤肉串
for Damp=0.01:0.02:0.1 %结构阻尼比取为0.05
Tc=0.0:0.02:15; %结构自振周期
Dt=0.02; %地震记录的步长
%记录计算得到的反应，MDis为最大相对位移，MVel为最大相对速度%MAcc某阻尼时最大绝对加速度，用于画图
MDis=zeros(1,length(Tc));
居住证明在哪里开MVel=zeros(1,length(Tc));
MAcc=zeros(1,length(Tc));
t=1; %在下一个循环中控制不同的结构自振周期
for T=0.0:0.02:15
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;
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*(((Damp/(Frcy*Dt)+1)*s/DamFrcy)+(1/(Frcy^2*Dt))*c)+1/(Frcy^2*Dt);
背酸
B(2,2)=e_t*((Damp/(Frcy*Dt)*s/DamFrcy)+(1/(Frcy^2*Dt))*c)-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(1,t)=max(abs(Displace));
MVel(1,t)=max(abs(Velocity));
if T==0.0
MAcc(1,t)=max(abs(Accelerate)); %当结构的自振周期为0时，其绝对加速度应等于地面加速度
el
MAcc(1,t)=max(abs(AbsAcce));
end
Displace=zeros(1,count);%初始化各储存向量避免下次计算时引用到前面结果
徘徊什么意思Velocity=zeros(1,count);
老姥姥AbsAcce=zeros(1,count);
t=t+1; %t=length(Tc)，即所求结构自振周期有多少个，对应就运行多少次。
end
% ***********绘图***********
figure %绘制位移反应谱
plot(Tc,MDis)
title('相对位移反应谱')
xlabel('周期(s)')
ylabel('位移(m)')
legend('ζ=0.05')
grid
figure %绘制速度反应谱
plot(Tc,MVel)
title('相对速度反应谱')
xlabel('周期(s)')
ylabel('速度')
legend('ζ=0.05')
grid
figure %绘制速度反应谱
plot(Tc,MAcc)
title('绝对加速度反应谱')
xlabel('T_{n}(s)')
ylabel('加速度 ')
legend('ζ=0.05')
蚕蛾的功效与作用grid
end
