pre.m

薄壁结构有限元计算

% 固有配率及振型计算
[xx,dd]=eig(K,M,'qz');
dd=sqrt(dd);
ww=sort(diag(dd));

% 阻尼矩阵形成
b=2*0.05/(ww(1)+ww(5));
a=b*ww(1)*ww(5);
C=a*M+b*K;

% 地震波输入
xg1=[(Zinput(1:500))' zeros(1,500);(Xinput(1:500))' zeros(1,500);(Yinput(1:500))' zeros(1,500);zeros(4,1000)];
xg=xg1;
for i=1:(n-1)
     xg=[xg;xg1];
end

% 广义位移初始化
x=zeros(7*n,1000);
v=x;
a=x;