8 files changed, 155 insertions, 0 deletions
diff --git a/2858/CH3/EX3.1/Ex3_1.sce b/2858/CH3/EX3.1/Ex3_1.sce
new file mode 100755
index 000000000..85351535e
--- /dev/null
+++ b/2858/CH3/EX3.1/Ex3_1.sce
@@ -0,0 +1,15 @@
+//example 3.1
+clc; funcprot(0);
+Nc=17.69;
+Nq=7.44;
+Ny=3.64;
+q=3*115;
+Gamma=115; //lb/ft^3
+c=320;
+B=5;//ft
+FS=4;//factor of safety
+qu=1.3*c*Nc+q*Nq+0.4*Gamma*B*Ny
+qall=qu/FS; //q allowed
+Q=qall*B^2;
+disp(Q,"allowable gross load in lb");
+
diff --git a/2858/CH3/EX3.2/Ex3_2.sce b/2858/CH3/EX3.2/Ex3_2.sce
new file mode 100755
index 000000000..b40e27e83
--- /dev/null
+++ b/2858/CH3/EX3.2/Ex3_2.sce
@@ -0,0 +1,28 @@
+//example 3.2
+clc; funcprot(0);
+Gamma=105;//lb/ft^3
+Gammasat=118;//lb/ft^3
+FS=3;
+pa=2014.125;//lb/ft^2
+Depth=[5,10,15,20,25];
+N60=[4,6,6,10,5];
+sigmao=[0,0,0,0,0];
+phi=[0,0,0,0,0]
+Gammaw=62.4;
+s=0;
+printf("depth (ft)\tN60\t  \tstress(lb/ft^2)\t phi(degrees)\n")
+for i=1:5
+    sigmao(i)=2*Gamma+(Depth(i)-2)*(Gammasat-Gammaw);
+    phi(i)=sqrt(20*N60(i)*sqrt(pa/sigmao(i)))+20;
+    printf(" %.2f\t       %.2f\t\t   %.2f  \t%.2f \n",Depth(i),N60(i),sigmao(i),phi(i));
+avgphi=phi(i)/5+s;
+s=avgphi;
+end
+disp(round(avgphi),"average friction angle in degrees");
+//using graph get the values of other terms in terms of B and solve for B
+deff('y=f(x)','y=-150000/x^2+5263.9+5527.1/x+228.3*x');
+[x, v, info ]=fsolve(4,f);
+disp(x,"the width in ft");
+
+
+
diff --git a/2858/CH3/EX3.3/Ex3_3.sce b/2858/CH3/EX3.3/Ex3_3.sce
new file mode 100755
index 000000000..591573558
--- /dev/null
+++ b/2858/CH3/EX3.3/Ex3_3.sce
@@ -0,0 +1,29 @@
+//example 3.3
+clc; funcprot(0);
+phi=25; //degrees
+Es=620; //kN/m^2
+Gamma=18;//kN/m^2
+Df=0.6;
+B=0.6;
+L=1.2;
+Fqc=0.347;
+Nq=10.66;
+Nc=20.72;
+Ngamma=10.88;
+mu=0.3;
+Fyd=1;
+c=48;//kN/m^2
+q=Gamma*(Df+B/2);
+Ir=Es/(2*(1+mu)*(c+q*tan(phi*%pi/180)));
+disp(Ir,"value of Ir");
+Fcc=Fqc-(1-Fqc)/(Nq*tan(phi*%pi/180));
+Fcs=1+Nq/Nc*B/L;
+Fqs=1+B/L*tan(phi*%pi/180);
+Fys=1-0.4*B/L;
+Fcd=1+0.4*Df/B;
+Fqd=1+2*tan(phi*%pi/180)*(1-sin(phi*%pi/180))^2*Df/B;
+q1=0.6*18;
+Fyc=Fqc;
+qu=c*Nc*Fcs*Fcd*Fcc+q1*Nq*Fqs*Fqd*Fqc+1/2*Gamma*Ngamma*Fys*Fyd*Fyc;
+disp(qu,"ultimate bearing capacity in kN/m^2");
+
diff --git a/2858/CH3/EX3.4/Ex3_4.sce b/2858/CH3/EX3.4/Ex3_4.sce
new file mode 100755
index 000000000..dd0aa235b
--- /dev/null
+++ b/2858/CH3/EX3.4/Ex3_4.sce
@@ -0,0 +1,21 @@
+//example 3.4
+clc; funcprot(0);
+q=110*4; //lb/ft^2
+Nq=33.3;
+phi=35;
+Df=4;
+B=6;
+Gamma=110;
+Ngamma=48.03;
+B1=6-2*0.5;
+Fqi=1;
+Fyi=1;
+Fyd=1;
+Fqs=1;
+Fys=1;
+Fqd=1+2*tan(phi*%pi/180)*(1-sin(phi*%pi/180))^2*Df/B;
+qu=q*Nq*Fqs*Fqd*Fqi+1/2*B1*Gamma*Ngamma*Fys*Fyd*Fyi;
+Qult=B1*1*qu;
+disp(Qult,"ultimate bearing capacity in lb/ft");
+disp(Qult/2000,"ultimate bearing capacity in ton/ft");
+
diff --git a/2858/CH3/EX3.5/Ex3_5.sce b/2858/CH3/EX3.5/Ex3_5.sce
new file mode 100755
index 000000000..60188d9ee
--- /dev/null
+++ b/2858/CH3/EX3.5/Ex3_5.sce
@@ -0,0 +1,14 @@
+//example 3.5
+clc; funcprot(0);
+e=0.5;
+B=6;
+k=e/B;
+Gamma=110;
+q=440;
+disp("get the values of Nqe and Nye from the figure from the value of e/B");
+Nye=26.8;
+Nqe=33.4;
+Qult=B*1*(q*Nqe+1/2*Gamma*B*Nye);
+disp(Qult,"ultimate bearing capacity in lb/ft");
+disp(Qult/2000,"ultimate bearing capacity in ton/ft");
+
diff --git a/2858/CH3/EX3.6/Ex3_6.sce b/2858/CH3/EX3.6/Ex3_6.sce
new file mode 100755
index 000000000..0bc1a627f
--- /dev/null
+++ b/2858/CH3/EX3.6/Ex3_6.sce
@@ -0,0 +1,22 @@
+//example 3.6
+clc; funcprot(0);
+Df=0.7;
+//from table
+Nq=18.4;
+Ny=22.4;
+q=12.6;
+phi=30; //angle
+L=1.5;
+Fyd=1;
+Gamma=18;
+L1=0.85*1.5;
+L2=0.21*1.5;
+B=1.5;
+A=1/2*(L1+L2)*B;
+B1=A/L1; //B'
+Fqs=1+B1/L1*tan(phi*%pi/180);
+Fys=1-0.4*B1/L1;
+Fqd=1+2*tan(phi*%pi/180)*(1-sin(phi*%pi/180))^2*Df/B;
+Qult=A*(q*Nq*Fqs*Fqd+1/2*Gamma*B1*Ny*Fys*Fyd);
+disp(Qult,"ultimate load in kN");
+
diff --git a/2858/CH3/EX3.7/Ex3_7.sce b/2858/CH3/EX3.7/Ex3_7.sce
new file mode 100755
index 000000000..443b4ad48
--- /dev/null
+++ b/2858/CH3/EX3.7/Ex3_7.sce
@@ -0,0 +1,15 @@
+//example 3.7
+clc; funcprot(0);
+e=0.15;
+B=1.5;
+Fqs=1;
+L=1.5;
+Gamma=18;
+q=0.7*18;
+//from table
+Nqe=18.4;
+Nye=11.58;
+Fys=1+(2*e/B-0.68)*(B/L)+(0.43-3/2*e/B)*(B/L)^2;
+Qult=B*L*(q*Nqe*Fqs+1/2*Gamma*Nye*Fys);
+disp(Qult,"ultimate load in kN");
+
diff --git a/2858/CH3/EX3.8/Ex3_8.sce b/2858/CH3/EX3.8/Ex3_8.sce
new file mode 100755
index 000000000..5da812803
--- /dev/null
+++ b/2858/CH3/EX3.8/Ex3_8.sce
@@ -0,0 +1,11 @@
+//example 3.8
+clc; funcprot(0);
+q=16;
+Nqei=14.2;
+Gamma=16
+B=1.5;
+Nyet=20;
+Qult=B*(Nqei*q+1/2*Gamma*B*Nyet);
+disp(Qult,"ultimate load in kN/m");
+
+