31.08.2018 Views

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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

236 MATLAB programs<br />

xlabel(titx);<br />

grid on;<br />

title(tit);<br />

pause;<br />

EasyPlot(‘SH stratum TD.ezp’, tit, titx, T, U, ‘r’);<br />

% Frequency doma<strong>in</strong> solution via method of images<br />

ni = 100; % Number of images<br />

fac = 0.25/mu/i;<br />

fmax = N*ffund;<br />

nf = 600;<br />

df = fmax/nf;<br />

F = [df:df:fmax];<br />

W = 2*pi*F;<br />

H = zeros(size(F));<br />

for j=-ni:ni<br />

t1 = sqrt(x2+(dz-j*h2)ˆ2)/cs;<br />

t2 = sqrt(x2+(sz-j*h2)ˆ2)/cs;<br />

H = H + (besselh(0,2,t1*W) - besselh(0,2,t2*W))*(-1)ˆj;<br />

end<br />

H = fac*H;<br />

plot (F, real(H));<br />

hold on;<br />

plot (F, imag(H),‘r’);<br />

tit = spr<strong>in</strong>tf(...<br />

‘Green function at x=%5.3f, z=%5.3f due to SH load <strong>in</strong> stratum at z=%5.3f’,...<br />

x, z, z0);<br />

grid on;<br />

title(tit);<br />

titx = ‘Frequency (rad/s)’;<br />

xlabel(titx);<br />

pause;<br />

hold off;<br />

EasyPlot(‘SH stratum FD1.ezp’, tit, titx, F, H, ‘c’);<br />

% Frequency doma<strong>in</strong> solution via normal modes<br />

nm = 20; % number of normal modes<br />

fac = 1/mu/i;<br />

tit = spr<strong>in</strong>tf(...<br />

‘Green function at x=%5.3f, z=%5.3f due to SH load <strong>in</strong> stratum at z=%5.3f’,...<br />

x, z, z0);<br />

x = abs(x/h);<br />

z = z/h;<br />

z0 = z0/h;<br />

W = W*h*(1-i*xi)/cs;<br />

H = zeros(size(W));<br />

W2 = W.ˆ2;<br />

for j=1:nm<br />

j1 = pi*(j-0.5);

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!