单自由度matlab时程分析程序,matlab单自由度的时程分析程序

clear;clc;

% 结构模型初始参数---------------------------------------------------------- m=3e3; %质量(单位：kg)

k=1e6; %刚度((单位：N/m))

kesai=0.05; %阻尼比取0.05

c=2*kesai*sqrt(k*m); %阻尼系数

% 读取地震波数据------------------------------------------------------------ acc=textread('D:\处理后的smc文件\51WCW_90_chnua370295.smc_090501.a','%f','headerlines',56); PGA_Max=max(abs(acc)) %最大地面加速度绝对值

% Newmark-beta法的基本参数--------------------------------------------------

beta=1/6; gama=0.5; %按线性加速度法计算更接近真实结果，故取此组参数

dt=0.02; %地震加速度时程波记录时间间隔

b1=1/(beta*dt^2); b2=1/(beta*dt); b3=1-1/(2*beta); %计算参数

b4=gama/(beta*dt); b5=gama/beta-1; b6=(1-gama/(2*beta)) *dt;

ke=k+m*b1+c*b4; %等效刚度

% 设定结构初始状态为零,生成向量空间存储计算值---------------------------------