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

218 MATLAB programs
T1 = conj(sqrt(T2-a2)); % make the imaginary part negative
q1 = c*T1+i*s*T;
% Complex horizontal slowness, P waves
p1 = c*T+i*s*T1; % ” vertical ” ” ”
Q1 = q1.ˆ2;
s1 = sqrt(Q1+1);
S1 = 2*Q1+1;
R1 = S1.ˆ2 - 4*Q1.*p1.*s1; % Rayleigh function, P waves
D1 = p1./T1./R1;
% Derivative of q1 divided by R1
Ux1 = 2*fac*sign(x)*imag(q1.*s1.*D1);
Uz1 = -fac*sign(z)*real(S1.*D1);
T = [T0,T];
Ux = [Ux0,Ux1];
Uz = [Uz0,Uz1];
% Plot results
plot (T,r*Ux);
tmax = max(T);
tit = sprintf(...
‘Horizontal displacement due to line blast load, \\nu=%5.2f’,pois);
title (tit);
grid on;
%axis([0 tmax -2 2]);
pause;
% Create EasyPlot file
fout = fopen(‘garvin-x.ezp’, ‘w’);
fprintf (fout, ‘/et g “%s”\n’, tit);
fprintf (fout, ‘/og on\n’);
fprintf (fout, ‘/sd off\n’);
fprintf (fout, ‘/sm off\n’);
fprintf (fout, ‘%15.5f %15.5f\n’, [T;r*Ux]);
fclose (fout);
plot (T,r*Uz);
tit = sprintf(...
‘Vertical displacement due to line blast load, \\nu=%5.2f’,pois);
title (tit);
grid on;
%axis([0 tmax -2 2]);
pause;
% Create EasyPlot file
fout = fopen(‘garvin-z.ezp’, ‘w’);
fprintf (fout, ‘/et g “%s”\n’, tit);
fprintf (fout, ‘/og on\n’);
fprintf (fout, ‘/sd off\n’);
fprintf (fout, ‘/sm off\n’);
fprintf (fout, ‘%15.5f %15.5f\n’, [T;r*Uz]);
fclose (fout);
close;