Shreya

%Bearing Capacity check

D_p= zeros;
avg= zeros;
Cu= zeros;
no_layers=input('Enter the total number of layers ');
text=cell(1,no_layers);
for i=1:no_layers
text1=input('Enter the layer type ','s');
text(i)=cellstr(text1);
s=strfind(text1,'Clay');
Gamma(i)=input('Enter the unit weight for this layer type in t/m^3 ');
Thk(i)=input('Enter the thickness for this layer type in m ');
if s~=0
C(i)=input('Enter the C value for this layer type in t/m^2 ');
mv(i)=input('Enter the mv for this layer type ');
N_value(i)=0;
phi(i)=0;
else
C(i)=0;
mv(i)=0;
N_value(i)=input('Enter the N value for this layer type ');
phi(i)=input('Enter the phi value for this layer type ');
end
end
Nc=input('Enter the Nc value');
FOS=input('Enter the FOS');
Problem_Type=input('Enter Problem Type(Dimension/Surcharge)','s');
PT=strcmp(Problem_Type,'Dimension');
if PT==1
load=input('Enter the load in t/m^2 ');
height=input('Enter the height in m ');
dia=input('Enter the diameter in m ');
spth=input('Enter the sand pad thickness in m ');
spuw=input('Enter the sand pad unit weight in t/m^2 ');
dia_fond=dia+spth;
Qn=(load*height*dia^2)/((dia+1)^2) + (spuw*spth);
else
surchrge=input('Enter the surcharge in t/m^2');
dia=input('Enter the diameter in m');
spth=input('Enter the sand pad thickness in m');
spuw=input('Enter the sand pad unit weight in t/m^2');
dia_fond=dia+spth;
Qn=(surchrge*dia^2)/((dia+1)^2) + (spuw*spth);
end
for i=1:no_layers
if Gamma(i)<1.9
C_BC(i)=C(i);
Thk_BC(i)=Thk(i);
else
C_BC(i)=0;
Thk_BC(i)=0;
end
end
for i=1:no_layers
avge(i)=C_BC(i)*Thk_BC(i);
end
Cu_avge=sum(avge)/sum(Thk_BC);
Qall= Cu_avge*Nc/FOS;
disp('Allowable bearing capacity of soil,qult in t/m^2 is :');
disp(Qall);
disp('Net foundation pressure,qn in t/m^2 is :');
disp(Qn);
if Qall>Qn
disp('Design is safe in bearing capacity!!');
else
disp('Design is not safe in B.C, Ground improvemnet is required.....Preloading and Sand drains will be used.');
end

%---------------------------------------------------

%gamma_sat calculation
text_lyr=char('DenseSand','FineSand','MediumSand','CoarseSand','SiltySand','FirmClay','SoftClay','MarineClay','
StiffClay','GreyClay');
Gamma_sat=[20.95172414 19.6 19.31176471 19.04 18.31052632 18.78333333 17.885 19.90625 19.60625
18.54054054];
C1=[39 52 39 43 24 0 0 0 0 0];
C2=[4.5 3.3 4.5 11.8 5.3 0 0 0 0 0];

for i=1:no_layers
E_sand(i)=C1(i)*C2(i)*N_value(i);
end

E=[1065 1016 1065 2200 1035 600*Cu_avge 600*Cu_avge 600*Cu_avge 600*Cu_avge 600*Cu_avge];

for i=1:no_layers
for j=1:10
Comp(i)=strcmp(text(i),text_lyr(j));
if Comp(i)==1
gammasat(i)=Gamma_sat(j);
E_lyr(i)=E(j);
else
end
 end
end
%--------------------------------------------------------------------

for i=1 : no_layers

if Gamma(i)< 1.9
chck(i)= 1;
else
chck(i)= 0;
end
end
Layer_consi=sum(chck) ;
% % --------------------------------------------------------------------
% % Preloading Calculations
SB_gamma= input ('Enter the unit weight of sand bag in t/m^3 ');
%mid layer depth
cum_sum=cumsum(Thk);
mid_layer_depth(1)=Thk(1)/2;
for i=2:no_layers
mid_layer_depth(i) = cum_sum(i-1)+Thk(i)/2;
end
Qult_imp=0;
c=C;
l=1;
%D_p;
while Qult_imp<Qn
B1=dia;
B2=dia;
pi=3.14;
Layer_consi;
for i=1:Layer_consi %Finding Factor for embankment loading
a1= atan((B1+B2)/(mid_layer_depth(i)))-atan(B1/mid_layer_depth(i));
a2= atan(B1/mid_layer_depth(i));
sf=2;
Iz(i)=(1/pi)*((B1+B2/B2)*(a1+a2)-(B1/B2*(a2))); %#ok<*SAGROW>
D_p(i)=2*SB_gamma*Iz(i);
end

H1_Embk=Qall/D_p(1);
for i=1:Layer_consi
delta_P(i)=H1_Embk*D_p(i);
delta_Cu(i)=0.3*delta_P(i);
Cu(i)= c(i)+delta_Cu(i);
end
Thk_layrs=0;
for i=1:Layer_consi
avg(i)=Cu(i)*Thk(i);
 Thk_layrs=Thk_layrs+Thk(i);
end
Cu_avg=sum(avg)/sum(Thk_layrs);
Qult_imp=Cu_avg*Nc/FOS;
c=Cu;
H_all_layers(1)=H1_Embk;
Qall=Qult_imp -(SB_gamma*sum(H_all_layers));
l= l+1;
end
disp('Number of preloading stages is : ')
disp(l-1)
disp('Height of preloading for each step, in meters is : ')
disp(H_all_layers)
disp('Final allowable bearing capacity after preloading in t/m^2 is : ')
disp(Qult_imp)
% % ----------------------------------------------------------------
% % Settlement Calculations without preloading
disp('Calculating Settlements without preloading applied.........')
% Immediate settlement of clay layers
Eclay= 600*Cu_avge;
S_imm_clay_cen=0.75*Qn*dia_fond/Eclay;
S_imm_clay_cor=S_imm_clay_cen*0.64;
disp('Immediate settlement of clay layers in mm is : ')
disp(S_imm_clay_cen*1000)
% Immediate settlement of Sand layers
% Finding p
for i = 1:no_layers
Pre_p(i) = Thk(i)*Gamma_sat(i);
end
Cum_sum_p=cumsum(Pre_p);
P(1)= Pre_p(1)/2;
for i=2 :no_layers
P(i) = Cum_sum_p(i-1)*Gamma_sat(i)*Thk(i)/2;
end
% finding delta P
mid_layer_thk= dia_fond+mid_layer_depth;
for i = 1: no_layers;
del_P (i) =Qn*((dia_fond^2)/mid_layer_thk(i)^2);
end
S_Imm_sand(i)=2.3*P(i)*Thk(i)*log10((P(i)+del_P(i))/P(i))/E(i);
for i=1:no_layers
if C(i)==0
index =i;
else
end
end
Imm_sett_sand= S_Imm_sand(i);
disp('Immediate settlement of sand layers in mm is : ')
 disp(1000*Imm_sett_sand)
Total_Imm_sett= Imm_sett_sand+S_imm_clay_cen;
disp('Total immendiate settlement of all layers in mm is : ')
disp(1000*Total_Imm_sett)
%Consolidation settlement of clay layers

for i=1:no_layers
S_Consi_set(i)= mv(i)*Thk(i)*del_P(i);
end

Total_consi_set=sum(S_Consi_set); %how many layers???

disp('Total consolidation of all layers in mm is : ')
disp(1000*Total_consi_set)

Total_Sett= Total_consi_set+Total_Imm_sett;
disp('Total settlement before preloading is applied, in mm, is: ')
disp(Total_Sett)

% %-----------------------------------------------------
% % Degree of consolidation & sand drains requirement
% fixing length of sand drain

max_conset=max(S_Consi_set);
for i=1:no_layers
H_factor(i)=S_Consi_set(i)/max_conset;
if H_factor(i)>0.25
H_drain_prt(i)=Thk(i);
else
H_drain_prt(i)=0;
end
end

H_drain=sum(H_drain_prt);
if H_drain>30
H_drain=30;
end
disp('Calculating time that will be required to achieve 90% degree of consolidation...')
H=0;
Cv=input('Enter the Cv value in m^2/day '); %0.02%m2/day
Cr=input('Enter the Cr value in m^2/day '); %0.04%m2/day

for i=1:Layer_consi
H=H+Thk(i);
end
time_cons= (0.848*H^2)/Cv;
if time_cons>365
 time_cons_year= time_cons/365;
disp('time as calculated to achieve 90% degree of consolidation, in year, is : ')
disp(time_cons_year)
else
time_cons_year= time_cons;

disp('Time as calculated to acheive 90% degree of consolidation, in days, is : ')

disp( time_cons_year)
end
if time_cons<200 % setting limit below which sand drain is not required
disp('No sand drains is required because time needed for 90% consolidation is very less! ');
else
disp('Since time required to aceive 90% degree of consolidation is very high, sand drains is required..');
an = 0;
numberOfElements = length(H_drain_prt);
for i=1:numberOfElements
ie(i) = isequal(H_drain_prt(i),an);
end
lyr = numberOfElements - sum(ie);
disp('Applying all conditions for fixing the height of sand drain and considering that since most of the
consoldation settlement is taking place upto layer ');
disp(lyr)
disp('Height of sand drain is choosen as ')
disp(H_drain)
end
%-----------------------------------------------------------------------------------------------------------------%sand
%Sand drains calculation
Drain_type = input('Enter the drain type(Sand Drains/Sand wicks) ', 's');
st = strfind(Drain_type,'Drains');
if st~=0
Ini_dia_SD = input('Enter the diameter of the drain in mm. ' );
else
Section = input('Enter the section type(Circular/Rectangular) ', 's');
st2 = strfind(Section, 'Circular');
if st~=0
Ini_dia_SW_C = input('Enter the diameter of the drain in mm. ');
else
Sub_type = input('Enter the sub type(PVC, Geo Drain, Membrane drain, Acid drain, Col band)','s');
Sub_type_cat = char('PVC', 'Geo drain', 'Membra Drain', 'Acid Drain', 'Col band');
Sub_type_Catalog = cellstr(Sub_type_cat);
Deff_wick = [64.93521678 63.02535746 62.38873769 67.48169587 67.48169587];

for i = 1:5
Comp(i) = strcmp(Sub_type, Sub_type_Catalog(i)); %#ok<*SAGROW>
if Comp(i)==1
Drain_dia_WR = Deff_wick(i);
else
 end
end
end
end

if st~=0
Drain_dia = Ini_dia_SD;
elseif st2~=0
Drain_dia = Ini_dia_SW_C;
else
Drain_dia = Drain_dia_WR;
end

Pattern = input('Enter the pattern (Triangular/Rectengular/Hexagonal) ','s');

Spacing = input('Enter the spacing for sand drains in mm ');
Min_spc = input('Enter the minimum spacing requirement in mm ');
time_fix = input('Enter the time period for applying GIT in days ');
U = 0;
m = 0;

while U <0.9
st3 = strfind(Pattern,'Triangular');
if st3~=0
D_inf = 1.05*Spacing;
else
st4 = strfind(Pattern,'Rectengular');
if st4~=0
D_inf = 1.13*Spacing;
else
D_inf = 1.29*Spacing;
end
end

n = D_inf/Drain_dia;
Fn = (n^2/((n^2)-1))*log(n)-((3*(n^2))-1)/(4*(n^2));
Tr = Cr*time_fix/(D_inf/1000)^2;
Ur = 1-exp(-(8*Tr/Fn));
Tv = Cv*time_fix/H^2;
Uv = sqrt(4*Tv/pi);
U = 1-(1-Ur)*(1-Uv);
Spacing = Spacing-100;
m = m+1;
if U>0.9
break
end
if Spacing < Min_spc
 disp('minimum spacing is reached with initial diameter , diameter will be increased')
Drain_dia= Drain_dia+25;
Spacing = (m)*100+Spacing;
end
end
Drain_dia;
Final_U=U;
disp('The final degree of consolidation achieved is :')
disp(100*Final_U)
Final_Spacing= Spacing +100;
disp('Revised diameter of the drain ,in mm, is')
disp(Drain_dia)
disp('Final spacing provided is :')
disp(Final_Spacing)
%_______________________________________________________
%Post treatment settlement calculations
disp('Calculating post treatment settlement calculations....')
Thk_layrs=0;
for i=1:Layer_consi
avgw(i)=C(i)*Thk(i);
end
C_org_avg=sum(avgw)/sum(Thk_layrs);
PT_Imm_set=0.75*SB_gamma*sum(H_all_layers)*dia_fond/(600*C_org_avg);
disp('Total immediate settlement in mm is:')
disp(1000*(PT_Imm_set))
%Post Treatment Consolidation Settlement
PT_Consi_sett = Final_U*Total_consi_set;
disp('Total Consolidation settlement in mm. is: ')
disp(1000*PT_Consi_sett)
PT_total_sett=PT_Imm_set+PT_Consi_sett;
disp('total settlement after ground improvement , in m,is :')
disp(PT_total_sett)
Sett_left=Total_Sett-PT_total_sett;
disp('settlement after removal of preload and constructiojn in mm is :')
disp(1000*Sett_left)
if Sett_left*1000<125
disp('Since settlement left is less than 125 mm so it is oK for residential buildings')
else
disp('Since settlement left is more than 125 mm so it is not safe for residential building ')
end
%____________________________________

Error using input

Cannot call INPUT from EVALC.

Error in githw2 (line 8)

no_layers=input('Enter the total number of layers ');

Published with MATLAB R2016a