earthquakeinput.m

薄壁结构有限元计算

% 输入地震波的数据（纵向、横向、垂向）
load X.dat;
load Y.dat;
load Z.dat;

% 将所输入的地震波数据重新排列为一维矩阵
Xinput=reshape(X',510,1);
Yinput=reshape(Y',510,1);
Zinput=reshape(Z',510,1);

% 寻找地震横波与纵波中的最大幅值
XYinput=(Xinput.^2+Yinput.^2).^0.5;
factorXY=max(XYinput);
factorZ=max(abs(Zinput));

% 将地震横波按规定把幅值缩减为0.1g，纵波幅值为横波的2/3
Xinput=Xinput/factorXY;
Yinput=Yinput/factorXY;
Zinput=(Zinput*2)/(factorZ*3);

% 画出每一向地震波的时程输入图
t=0.02:0.02:10;

subplot(221),plot(t,Xinput(1:500));
title('Longitude Input');
grid on;

subplot(222),plot(t,Yinput(1:500));
title('Tranveral Input');
grid on;

subplot(223),plot(t,Zinput(1:500));
title('Vertical Input');
grid on;

% 画出底部震动平台的三维地震波加速度轨迹图
subplot(224),plot3(Xinput,Yinput,Zinput);
title('3-D Input');
grid on;