结构动力学MATLAB编程例题,华工结构动力学编程 MATLAB

华工结构动力学编程

一．所编程序及编译：

clear%***********读入地震记录***********

fid=fopen('ei.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);%绝对加速度

Damp=0.05;%结构阻尼比取为0.05

Tc=0.0:0.05:10;%结构自振周期

Dt=0.02;%地震记录的步长

%记录计算得到的反应，MDis为最大相对位移，MVel为最大相对速度%MAcc某阻尼时最大绝对加速度，用于画图

MDis=zeros(1,length(Tc));

MVel=zeros(1,length(Tc));

MAcc=zeros(1,length(Tc));

t=1;%在下一个循环中控制不同的结构自振周期

forT=0.0:0.05:10

Frcy=2*pi/T;%结构自振频率

DamFrcy=Frcy*sqrt(1-Damp*Damp);%计算公式化简

e_t=exp(-Damp*Frcy*Dt);

s=sin(DamFrcy*Dt);