From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 2858/CH9/EX9.1/Ex9_1.sce | 39 +++++++++++++++++++++++++++++++++++++++
 2858/CH9/EX9.2/Ex9_2.sce | 32 ++++++++++++++++++++++++++++++++
 2858/CH9/EX9.3/Ex9_3.sce | 34 ++++++++++++++++++++++++++++++++++
 2858/CH9/EX9.4/Ex9_4.sce | 25 +++++++++++++++++++++++++
 2858/CH9/EX9.5/Ex9_5.sce | 19 +++++++++++++++++++
 2858/CH9/EX9.6/Ex9_6.sce | 31 +++++++++++++++++++++++++++++++
 2858/CH9/EX9.7/Ex9_7.sce | 26 ++++++++++++++++++++++++++
 2858/CH9/EX9.8/Ex9_8.sce | 22 ++++++++++++++++++++++
 8 files changed, 228 insertions(+)
 create mode 100755 2858/CH9/EX9.1/Ex9_1.sce
 create mode 100755 2858/CH9/EX9.2/Ex9_2.sce
 create mode 100755 2858/CH9/EX9.3/Ex9_3.sce
 create mode 100755 2858/CH9/EX9.4/Ex9_4.sce
 create mode 100755 2858/CH9/EX9.5/Ex9_5.sce
 create mode 100755 2858/CH9/EX9.6/Ex9_6.sce
 create mode 100755 2858/CH9/EX9.7/Ex9_7.sce
 create mode 100755 2858/CH9/EX9.8/Ex9_8.sce

(limited to '2858/CH9')

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")
+
-- 
cgit