Hi Mabin,
I assume you want to know why your code is not performing as you had expected.
Your code has a problem since you calculate the "sigprimeh" value inside the innermost loop, overwriting it after each iteration and only charting the final value. After computing the total for each "zii", you must relocate the assignment statement outside the innermost loop in order to plot the cumulative sum for each "zii" value. Also you should use an array to store the values of "sigmaprimeh" as that will preserve all the calculated values. For example:
clear;
clc;
gamma = 95;
zwater = 0;
H = 10;
TYP = 1;
x = 8;
v = 0.49;
array = 2;
arraysize = 6;
Q = 250;
Q = (Q/array^2);
nax = array;
nay = array;
xmin = x - arraysize/2;
xmax = x + arraysize/2;
yleft = -(arraysize/2);
yright = (arraysize/2);
z = 0:0.1:H;
sigma_h = zeros(1,length(z)); % Preallocate memory for sigma_h
xa = linspace(xmin+(1/(2*array)*arraysize),xmax-(arraysize/(2*array)),nax);
ya = linspace(yleft+(1/(2*array)*arraysize),yright-(arraysize/(2*array)),nay);
for ii = 1:length(z)
sum =0;
for xii = xa
for yii = ya
R = sqrt(xii^2+yii^2+z(ii)^2);
sigprimeh = ((Q/(R*pi))*((3*z(ii)*xii^2)/R^3)*((R*(1-2*v))/(R+z(ii))));
sum = sum +sigprimeh;
end
end
sigma_h(ii) = sum;
end
plot(sigma_h,z)
title('Stress Due to Surcharge Load')
xlabel('\sigma`h (ksf)');
ylabel('Depth "z"(ft)');
set(gca,'YDir','reverse')
This will compute the accumulated sum for each "zii" value and store it in the "sigma_h" array, which can then be plotted against "z".
I hope this answers your question.
Regards,
Tushar
