initial commit / add all books

8 files changed, 228 insertions, 0 deletions

diff --git a/2858/CH9/EX9.1/Ex9_1.sce b/2858/CH9/EX9.1/Ex9_1.sce
new file mode 100755
index 000000000..dfb3ea991
--- /dev/null
+++ b/2858/CH9/EX9.1/Ex9_1.sce
@@ -0,0 +1,39 @@
+//example 9.1

+clc; funcprot(0);

+sall=30;//sigma allowed

+pi=%pi;

+zbar=12.1;

+L1=10;

+L2=20;

+Gamma=0.12;

+Gamma1=0.1294-0.0624;

+phi=40*pi/180;

+Ka=(tan(pi/4-phi/2))^2;

+Kp=(tan(pi/4+phi/2))^2;

+s1=Gamma*L1*Ka;//sigma1

+s2=Gamma*L1*Ka+Gamma1*L2*Ka;//sigma2

+L3=s2/(Gamma1*(Kp-Ka));

+disp(L3,"length in ft");

+P=1/2*s1*L1+s1*L2+1/2*(s2-s1)*L2+1/2*s2*L3;

+disp(P,"force in kip/ft");

+s5=Gamma*L1*Kp+Gamma1*L2*Kp+Gamma*L3*(Kp-Ka);//sigma5

+disp(s5,"pressure in kip/ft");

+A1=s5/(Gamma1*(Kp-Ka));

+A2=8*P/(Gamma1*(Kp-Ka))

+A3=6*P*(2*zbar*(Gamma1*(Kp-Ka))+s5)/(Gamma1*(Kp-Ka))^2

+A4=P*(6*zbar*s5+4*P)/(Gamma1*(Kp-Ka))^2

+printf("A1,A2,A3,A4 respectively in %.2f,%.2f,%.2f,%.2f \n",A1,A2,A3,A4);

+disp("slight error due to rounding off error")

+//partb

+deff('y=f(x)','y=x^4+41.7*x^3-270.5*x^2-13363*x-106863');

+[x]=fsolve(20,f);

+D=1.88+x;

+disp(D,"value of D, ft")

+TL=L1+L2+1.3*D;

+disp(TL,"total length in ft");

+//partc

+z=sqrt(2*P/(Gamma1*(Kp-Ka)));//zdash

+Mmax=P*(z+zbar)-1/2*(Gamma1*(Kp-Ka))*z^2*z/3;

+S=Mmax*12/sall;

+disp(S,"section modulus in in^3/ft")

+

diff --git a/2858/CH9/EX9.2/Ex9_2.sce b/2858/CH9/EX9.2/Ex9_2.sce
new file mode 100755
index 000000000..271c01569
--- /dev/null
+++ b/2858/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,32 @@
+//example 9.2

+clc; funcprot(0);

+sall=172.5*1000;//sigma allowed

+pi=%pi;

+c=47;

+zbar=1.78;

+L1=2;

+L2=3;

+Gamma=15.9;

+Gamma1=19.33-9.81;

+phi=32*pi/180;

+Ka=(tan(pi/4-phi/2))^2;

+Kp=(tan(pi/4+phi/2))^2;

+s1=Gamma*L1*Ka;//sigma1

+s2=Gamma*L1*Ka+Gamma1*L2*Ka;//sigma2

+P=1/2*s1*L1+s1*L2+1/2*(s2-s1)*L2;

+disp(P,"force in kN/ft");

+deff('y=f(x)','y=127.4*x^2-104.4*x-357.15');

+[x]=fsolve(2,f);

+D=x;

+disp(D,"value of D, m")

+disp(D*1.5,"actual D in m")

+L4=D*(4*c-(Gamma*L1+Gamma1*L2)-P/D)/4/c;

+disp(L4,"length in m");

+s6=4*c-(Gamma*L1+Gamma1*L2);//sigma6

+s7=4*c+(Gamma*L1+Gamma1*L2);//sigma7

+z=P/s6;//zdash

+Mmax=P*(z+zbar)-1/2*s6*z^2;

+S=Mmax*12/sall;

+disp(S,"section modulus in m^3/m")

+disp("slight error due to rounding off error")

+

diff --git a/2858/CH9/EX9.3/Ex9_3.sce b/2858/CH9/EX9.3/Ex9_3.sce
new file mode 100755
index 000000000..9304a6e68
--- /dev/null
+++ b/2858/CH9/EX9.3/Ex9_3.sce
@@ -0,0 +1,34 @@
+//example 9.3

+clc; funcprot(0);

+pi=%pi;

+zbar=2.23;

+L1=2;

+L2=3;

+Gamma=15.9;

+Gamma1=19.33-9.81;

+phi=32*pi/180;

+Ka=(tan(pi/4-phi/2))^2;

+Kp=(tan(pi/4+phi/2))^2;

+s1=Gamma*L1*Ka;//sigma1

+s2=Gamma*L1*Ka+Gamma1*L2*Ka;//sigma2

+L3=s2/(Gamma1*(Kp-Ka));

+disp(L3,"length in m");

+P=1/2*s1*L1+s1*L2+1/2*(s2-s1)*L2+1/2*s2*L3;

+disp(P,"force in kN/m");

+deff('y=f(x)','y=x^3+6.99*x^2-14.55');

+[x]=fsolve(1.4,f);

+D=L3+x;

+disp(D,"value of D in m")

+AL=1.4*D;

+disp(AL,"actual length in m");

+//partb

+L4=1.4;

+F=P-1/2*(Gamma1*(Kp-Ka)*L4^2);

+disp(F,"anchor force in kN/m");

+//partc

+deff('y=f(x)','y=x^2+6.682*x-14.44');

+[x]=fsolve(1.7,f);

+z=x+2;

+Mmax=-1/2*s1*L1*(x+2/3)+F*(x+1)-s1*x*x/2-1/2*Ka*Gamma1*x^3/3;

+disp(Mmax,"maximum moment in kN-m/m")

+

diff --git a/2858/CH9/EX9.4/Ex9_4.sce b/2858/CH9/EX9.4/Ex9_4.sce
new file mode 100755
index 000000000..7ed1658ac
--- /dev/null
+++ b/2858/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,25 @@
+//example 9.4

+clc; funcprot(0);

+Gamma=15.9;

+Gamma1=19.33-9.81;

+GD=0.23;

+CDL1=1.172;

+L1=2;

+L2=3;

+Dth=(L1+L2)*GD*CDL1;

+disp(Dth,"theoritical depth in m");

+Dac=1.4*Dth;

+disp(Dac,"actual depth in m");

+printf("approximate it as 2 m");

+//partb

+CFL1=1.073;

+GF=0.07;

+Gammaa=(Gamma*L1^2+Gamma1*L2^2+2*Gamma*L1*L2)/(L1+L2)^2;

+F=Gammaa*(L1+L2)^2*GF*CFL1;

+disp(F,"force in kN/m");

+//partc

+GM=0.021;

+CML1=1.036;

+Mmax=Gammaa*(L1+L2)^3*GM*CML1;

+disp(Mmax,"maximum moment in kN-m/m");

+

diff --git a/2858/CH9/EX9.5/Ex9_5.sce b/2858/CH9/EX9.5/Ex9_5.sce
new file mode 100755
index 000000000..f68f2e92f
--- /dev/null
+++ b/2858/CH9/EX9.5/Ex9_5.sce
@@ -0,0 +1,19 @@
+//example 9.5

+clc; funcprot(0);

+Mmax=43.72;

+sp=["PSA-31","PSA-23"];//sheet pile

+H=[7.9,7.9]

+I=[4.41e-6,5.63e-6];

+p=[0.00466,0.00365];

+S=[10.8e-5,12.8e-5];

+Md=[18.63,22.08];

+printf("SheetPile    I(m^4/m)  H(m)  p\t     Logp  S(m^3/m)  Md(kN-m/m) Md/Mmax \n")

+for i=1:2

+    Logp(i)=log10(p(i));

+    k(i)=Md(i)/Mmax;

+    printf("%s\t   %.2e    %.2f %f %.2f %.2e   %.2f\t %.3f\n",sp(i),I(i),H(i),p(i),Logp(i),S(i),Md(i),k(i));

+    

+end

+plot(Logp,k);

+xtitle("","LogP","Md/Mmax");

+

diff --git a/2858/CH9/EX9.6/Ex9_6.sce b/2858/CH9/EX9.6/Ex9_6.sce
new file mode 100755
index 000000000..85972c5a6
--- /dev/null
+++ b/2858/CH9/EX9.6/Ex9_6.sce
@@ -0,0 +1,31 @@
+//example 9.6

+clc; funcprot(0);

+pi=%pi;

+R=0.6;

+L1=10;

+L2=20;

+Gammasat=122.4;

+l1=5;

+Gamma=110;

+C=0.68;

+L=L1+L2;;

+Gammaw=62.4;

+Gamma1=Gammasat-Gammaw;//gammadash

+Gammaav=(Gamma*L1+Gamma1*L2)/(L1+L2);

+phi=35*pi/180;

+Ka=(tan(pi/4-phi/2))^2;

+sa=C*Ka*Gammaav*L;//sigmaa

+sp=R*sa;//sigmap

+deff('y=f(x)', 'y=x^2+50*x-1000');

+[x]=fsolve(15,f);

+D=x;

+disp(D,"depth in ft");

+R=L/D*(L-2*l1)/(2*L+D-2*l1);

+printf("value of R=%.2f is OK\n",R);

+//partb

+F=sa*(L-R*D);

+disp(F,"Force in lb/ft")

+//partc

+Mmax=0.5*sa*L^2*((1-R*D/L)^2-2*l1/L*(1-R*D/L));

+disp(Mmax,"maximum moment lb-ft/ft");

+

diff --git a/2858/CH9/EX9.7/Ex9_7.sce b/2858/CH9/EX9.7/Ex9_7.sce
new file mode 100755
index 000000000..038803fb0
--- /dev/null
+++ b/2858/CH9/EX9.7/Ex9_7.sce
@@ -0,0 +1,26 @@
+//example 9.7

+clc; funcprot(0);

+pi=%pi;

+zbar=3.2;

+c=41;

+L1=3;

+L2=6;

+Gamma=17;

+Gamma1=20-9.81;

+phi=35*pi/180;

+Ka=(tan(pi/4-phi/2))^2;

+Kp=(tan(pi/4+phi/2))^2;

+s1=Gamma*L1*Ka;//sigma1

+s2=Gamma*L1*Ka+Gamma1*L2*Ka;//sigma2

+P=1/2*s1*L1+s1*L2+1/2*(s2-s1)*L2;

+disp(P,"Force in kN/m");

+s6=4*c-(Gamma*L1+Gamma1*L2);//sigma6

+disp(s6,"pressure in kN/m^2");

+deff('y=f(x)','y=x^2+15*x-25.43');

+[x]=fsolve(1.6,f);

+D=x;

+disp(D,"depth in m");

+F=P-s6*D;

+disp(F,"force in kN/m");

+printf("slight error due to rounding off")

+

diff --git a/2858/CH9/EX9.8/Ex9_8.sce b/2858/CH9/EX9.8/Ex9_8.sce
new file mode 100755
index 000000000..9ee170211
--- /dev/null
+++ b/2858/CH9/EX9.8/Ex9_8.sce
@@ -0,0 +1,22 @@
+//example 9.8

+clc; funcprot(0);

+pi=%pi;

+Gamma=105;

+Cov=14;

+B=15/12;

+Ka=0.26;

+phi=35*pi/180;

+H=37.5/12;

+h=15/12;

+t=6/12;

+Gc=150;//gamma concrete

+W=H*t*Gc;

+k=4.5;//kp*cos(delta)

+Pu=1/2*Gamma*H^2*(k-Ka*cos(phi));

+disp(Pu,"force in lb/ft");

+Pus=[(Cov+1)/(Cov+H/h)]*Pu;

+disp(Pus,"force in lb/ft");

+Be=0.227*(H+h)+B;

+Pu=Pus*Be;

+disp(Pu,"resistance of anchor plate i lb/ft")

+