transientresponse3.m

薄壁结构有限元计算

% Newmark 动力时程积分法

% 输入计算参数gamma、beta、叠代步长deltat
gamma=input('Please input the value of gamma coefficient(recommend 0.5): ');
beta=input('Please input the value of beta coefficient(recommend 0.25): ');
deltat=input('Please input the value of time interval: ');

% 计算八个系数0.
a1=1/beta/deltat/deltat;
a2=gamma/beta/deltat;
a3=1/beta/deltat;
a4=0.5/beta-1;
a5=gamma/beta-1;
a6=(gamma/beta/2-1)*deltat;
a7=deltat*(1-gamma);
a8=gamma*deltat;

% 叠代计算每一步长中的位移、速度、加速度
for i=1:999
    x(:,i+1)=(a1*M+a2*C+K)\(-M*xg(:,i+1)+M*(a1*x(:,i)+a3*v(:,i)+a4*a(:,i))+C*(a2*x(:,i)+a5*v(:,i)+a6*a(:,i)));
    a(:,i+1)=a1*(x(:,i+1)-x(:,i))-a3*v(:,i)-a4*a(:,i);
    v(:,i+1)=v(:,i)+a7*a(:,i)+a8*a(:,i+1);
end