Eduardo Kausel-Fundamental solutions in elastodynamics_ a compendium-Cambridge University Press (2006)

222 MATLAB programs
if k >=1; D2(1:k)=0; end
S1 = S1.*D1;
S2 = S2.*D2;
D1 = 2*D1.*s1;
D2 = 2*D2.*s2;
% Displacements due to impulsive line load
sgn = fac*sign(x)*sign(z);
Uxx = fac*real(p2.*S2-Q1.*D1);
Uzz = fac*real(p1.*S1-Q2.*D2);
Uzx = sgn*imag(q1.*p1.*D1-q2.*S2);
Uxz = sgn*imag(q1.*S1-q2.*p2.*D2);
end
T = [T0,T];
Uxx = [Uxx0,Uxx];
Uzz = [Uzz0,Uzz];
Uxz = [Uxz0,Uxz];
Uzx = [Uzx0,Uzx];
% Plot results
plot (T,Uxx);
tmax = max(T);
tit = sprintf(...
‘Horizontal displacement due to horizontal line load, \\nu=%5.2f’,...
pois);
titx = ‘Time’;
title (tit);
grid on;
%axis([0 tmax -2 2]);
pause;
plot (T,Uzz);
tit = sprintf((...
‘Vertical displacement due to vertical line load, \\nu=%5.2f’,pois);
title (tit);
grid on;
%axis([0 tmax -2 2]);
pause;
plot (T,Uzx);
tit = sprintf((...
‘Vertical displacement due to horizontal line load, \\nu=%5.2f’,pois);
title (tit);
grid on;
%axis([0 tmax -0.5 0.5]);
if theta > mang
hold on
plot ([tr tr], [0 -wr]);
end
hold off;
pause;