% fig_08_front
% Normalized displacement components, ux/b and uy/b,
% on a grid near the edge dislocation
% Weertman & Weertman (1964)
% equations (8.36) and (8.37)

clear all, clf reset; % clear memory and figures
b = 1; % Burgers vector
mu = 30000; pr = 0.25; lm = (2*mu*pr)/(1-2*pr); % Elastic moduli
c1 = (0.5*mu)/(lm+2*mu); c2 = (lm+mu)/(lm+2*mu);
x = linspace(-50*b,50*b,100)+eps;
y = linspace(-50*b,50*b,100);
[X,Y]=meshgrid(x,y); % define Cartesian grid
DEN = X.^2 + Y.^2;
UXDB = (-1/(2*pi))*(atan2(Y,X) + c2*(X.*Y)./DEN);
UYDB = (-1/(2*pi))*(-c1*log(DEN) + c2*(Y.^2)./DEN);
[T,R] = cart2pol(X,Y);
UXDB(find(R<(5*b))) = nan; UYDB(find(R<(5*b))) = nan;
surf(X,Y,UXDB)
xlabel('x-axis'), ylabel('y-axis'), zlabel('displacement ux')
figure, contourf(X,Y,UXDB,10), colorbar, title('displacement ux')
xlabel('x-axis'), ylabel('y-axis')
figure, surf(X,Y,UYDB)
xlabel('x-axis'), ylabel('y-axis'), zlabel('displacement uy')
figure, contourf(X,Y,UYDB,10), colorbar, title('displacement uy')
xlabel('x-axis'), ylabel('y-axis')