newmarkseismicresponse.m

that procedure has good adaptability

K=100*10^3;
M=128.5;
C=1125;
 %位移 速度 加速度初值设定
d1=0;
v1=0;
a1=0;

%外荷载向量
dt=0.002;
tend=4;
omg=25;
P=500;
t0=0;
f=zeros(1,floor(tend/dt)+1);
for i=2:1:(floor(tend/dt)+1);
    t(1,i)=t0+dt*(i-1);
    f(1,i)=P*sin(omg*t(1,i));
end

%调用NewMark-beta法求解位移 速度 加速度
[d,v,a]=Newmark(K,M,C,f,d1,v1,a1,dt,tend);
%图形显示
figure(1)
plot(t(1,:),d(1,:))
figure(2)
plot(t,v(1,:))
figure(3)
plot(t,a(1,:)) %加速度显示



nodedd=d(1,:)';
nodevv=v(1,:)';
nodeaa=a(1,:)';

d