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

function [ ] = Torsion Half(r, z, z0) 231
W2 = i*2*pi*R2/cs*f;
%H = S1*(1+W1).*exp(-W1)+S2*(1+W2).*exp(-W2);
H = S1*(1+W1).*exp(-W1)+S2*(1+W2).*exp(-W2);
plot (f, real(H));
hold on;
plot(f,imag(H),‘r’);
tit = sprintf(...
‘Transfer function at r=%5.3f, z=%5.3f due to torsional load at z=%5.3f’,...
r, z, z0);
title (tit);
xlabel(titx);
grid on;
pause;
hold off;
EasyPlot(‘TorHalf FD.ezp’, tit, titx, f, H, ‘c’);
A2 = abs(H);
plot(f, A2);
grid on;
tit = ‘Abs. Val. transfer function, direct’;
title (tit);
xlabel(titx);
pause;
ang2 = unwrap(angle(H));
plot(f, ang2);
grid on;
title (‘Angle transfer function, direct’);
xlabel(titx);
pause;
plot(f, A1./A2);
grid on;
title (‘A1/A2’);
xlabel(titx);
pause;
ang = (ang2-ang1)/pi;
plot(f, ang);
xlabel(titx);
grid on;
title (‘Diff. Angle / pi’);
xlabel(titx);
pause;
close all;
return