initial commit / add all books

10 files changed, 170 insertions, 0 deletions

diff --git a/2744/CH8/EX8.1/Ex8_1.sce b/2744/CH8/EX8.1/Ex8_1.sce
new file mode 100755
index 000000000..a0c4b9e37
--- /dev/null
+++ b/2744/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,14 @@
+clear;
+clc;
+d = 6;//inches
+b = 3/4;//inch
+P = 18;//tons
+e = 1/8;//inch
+A = b*d;//sq.in
+M = P*e;//ton-in
+Z = (1/6)*b*d^2;//in^3
+p_0 = P/A;// tons/in^2
+p_b = M/Z;// ton/in^2
+p_max = p_0+p_b;// tons/in^2
+p_min = p_0-p_b;// tons/in^2
+printf('p_max = %.1f tons/in^2.,tensile\n p_min = %.1f tons/in^2.,tensile',p_max,p_min);
diff --git a/2744/CH8/EX8.10/Ex8_10.sce b/2744/CH8/EX8.10/Ex8_10.sce
new file mode 100755
index 000000000..e37f3c44b
--- /dev/null
+++ b/2744/CH8/EX8.10/Ex8_10.sce
@@ -0,0 +1,29 @@
+clear;

+clc;

+h = 20;// feet

+b = 4;// feet

+d = 12;// feet

+h1 = 18;// feet

+//density of masonary

+rho_m = 140;// lb-ft^3

+//density of water

+rho_w = 62.5;// lb-ft^3

+W = rho_m*0.5*(b+d)*h;// lb-wt

+//to locate its line of action divide the dam section into a rectangle and a triangle 

+x1 = b/2;//feet

+a1 = b*h;// sq. feet

+a2 = 0.5*(d-b)*h;// sq. feet

+x2 = b+((d-b)/3);// feet

+x_bar = (a1*x1+a2*x2)/(a1+a2);// feet

+P = rho_w*h1^2/2;// lb-wt

+z = x_bar + (h1/3)*(P/W);// feet

+e = z - d/2;// feet

+p_0 = W/d;// lb/ft^2

+M = W*e;// lb-feet

+Z = 1*d^2*1/6;// ft^3

+p_b = M/Z;// lb-ft^2

+p_max = p_0 + p_b;// lb-ft^2

+p_min = p_0 - p_b;// lb-ft^2

+printf('p_max = %.1f lb/ft^2., compressive at B.\n p_min = %.1f lb/ft^2., compressive at A.',p_max,p_min);

+

+//there is an error in the answer given in text book

diff --git a/2744/CH8/EX8.2/Ex8_2.sce b/2744/CH8/EX8.2/Ex8_2.sce
new file mode 100755
index 000000000..8013b8cf9
--- /dev/null
+++ b/2744/CH8/EX8.2/Ex8_2.sce
@@ -0,0 +1,15 @@
+clear;

+clc;

+d1 = 12;//inches

+t = 1;//inch

+d2 = d1-2*t;//inches

+P = 5;//tons

+e = 12;//inch

+A = 0.25*%pi*(d1^2-d2^2);//sq.in

+M = P*e;//ton-in

+Z = %pi*(d1^4-d2^4)/(32*d1);//in^3

+p_0 = P/A;// tons/in^2

+p_b = M/Z;// ton/in^2

+p_max = p_0+p_b;// tons/in^2

+p_min = p_0-p_b;// tons/in^2

+printf('p_max = %.4f ton/in^2.,cmopressive\n p_min = %.4f ton/in^2.,\n i.e., %.4f ton/in^2., tensile',p_max,p_min,-p_min);

diff --git a/2744/CH8/EX8.3/Ex8_3.sce b/2744/CH8/EX8.3/Ex8_3.sce
new file mode 100755
index 000000000..22808198c
--- /dev/null
+++ b/2744/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,24 @@
+clear;

+clc;

+l = 6;//inches

+b = 4;//inches

+d = 1/2;//inch

+P = 10;//tons

+r = 1.5;//inches

+A = 4.771;// in^2

+J = 0.968;// inches

+I_xx = 6.07;// in^4

+I_yy = 8.64;// in^4

+e = r-J;// inches

+M = P*e;// ton-inches

+y_t = J;//inches 

+y_c = b-y_t;//inches

+//compressive

+p_c = M*y_c/I_xx;// tons/in^2

+//tensile

+p_t = M*y_t/I_xx;// tons/in^2

+//compressive

+p_0 = P/A ;// tons/in^2

+p_max = p_0+p_c;// tons/in^2

+p_min = p_0-p_t;// tons/in^2

+printf('p_max = %.3f tons/in^2.,compressive\n p_min = %.3f tons/in^2.,compressive',p_max,p_min);

diff --git a/2744/CH8/EX8.4/Ex8_4.sce b/2744/CH8/EX8.4/Ex8_4.sce
new file mode 100755
index 000000000..12568c949
--- /dev/null
+++ b/2744/CH8/EX8.4/Ex8_4.sce
@@ -0,0 +1,24 @@
+clear;
+clc;
+b = 5;//inches
+t = 1/2;//inch
+P = 12;//tons
+d = 1/2;//inch
+r = 3/2;//inch
+A = (b-d)*t;// in^2
+p_0 = P/A;// tons/in^2
+a1 = b*d;// in^2
+x1 = 0;
+a2 = d*t;// in^2
+x2 = -3/2;
+e = (a1*x1-a2*x2)/(a1-a2);//inches
+M = P*e;// ton-inches
+y_c = a1-e;// inches
+y_t = a1+e;// inches
+I_yy = (t*b^3)/12 - ((t*d^3)/12 + 0.5*d*r^2);// in^4
+I_GG = I_yy - (a1-a2)*(e^2);// in^4
+p_c = M*y_c/I_GG;// tons/in^2
+p_t = M*y_t/I_GG;// tons/in^2
+p_max = p_0 + p_t;// tons/in^2
+p_min = p_0 - p_c;// tons/in^2
+printf('p_max = %.1f tons/in^2., tensile\n p_min = %.2f tons/in^2., tensile',p_max,p_min);
diff --git a/2744/CH8/EX8.5/Ex8_5.sce b/2744/CH8/EX8.5/Ex8_5.sce
new file mode 100755
index 000000000..1f8576741
--- /dev/null
+++ b/2744/CH8/EX8.5/Ex8_5.sce
@@ -0,0 +1,15 @@
+clear;

+clc;

+h = 20;// feet

+b = 12;// feet

+d = 4;// feet

+p = 30;// lb. per sq.foot

+rho = 140;// lb. per cubic foot

+p_0 = rho*h;// lb-ft^2

+P = p*b*h;// lb-wt

+M = P*h/2;//lb-ft

+Z = b*d^2/6;// ft^3

+p_b = M/Z;// lb/ft^2

+p_max = p_0 + p_b;// lb/ft^2

+p_min = p_0 - p_b;// lb/ft^2

+printf('p_max = %d lb/ft^2., compressive\n p_min = %d lb/ft^2., compressive',p_max,p_min);

diff --git a/2744/CH8/EX8.6/Ex8_6.sce b/2744/CH8/EX8.6/Ex8_6.sce
new file mode 100755
index 000000000..458650568
--- /dev/null
+++ b/2744/CH8/EX8.6/Ex8_6.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+h = 80;// feet

+p = 28;// lb. per sq.foot

+rho = 126;// lb. per cubic foot

+p_0 = rho*h/2240;// tons-ft^2

+p_max = 7;// tons/ft^2

+d = 4;// feet

+p_b = p_max - p_0;// tons/ft^2

+D = sqrt(3*p*h^2 /(2*p_b*2240) +sqrt(d^4 + (3*p*h^2 /(2*p_b*2240))^2));// feet

+t = 0.5*(D-d);// feet

+printf('The necessary thickness is, t = %d feet',round(t));

diff --git a/2744/CH8/EX8.7/Ex8_7.sce b/2744/CH8/EX8.7/Ex8_7.sce
new file mode 100755
index 000000000..26a476185
--- /dev/null
+++ b/2744/CH8/EX8.7/Ex8_7.sce
@@ -0,0 +1,17 @@
+clear;

+clc;

+h = 60;//feet

+rho = 130;// lb. per cubic foot

+D = 12;//feet

+d = 5;//feet

+P_h = 24;// lb. per sq. foot

+p_0 = rho*h;// lb-ft^2

+P = P_h*D*h;// lb-wt

+M = P*h/2;// lb-feet

+Z = %pi*(D^4 - d^4)/(32*D);//lb-ft^3

+p_b = M/Z;// lb/ft^2

+p_max = p_0 + p_b;// lb/ft^2

+p_min = p_0 - p_b;// lb/ft^2

+printf('p_max = %d lb/ft^2., compressive\n p_min = %d lb/ft^2., compressive',p_max,p_min);

+

+//there is an error in the answer given in text book

diff --git a/2744/CH8/EX8.8/Ex8_8.sce b/2744/CH8/EX8.8/Ex8_8.sce
new file mode 100755
index 000000000..1b30c86f5
--- /dev/null
+++ b/2744/CH8/EX8.8/Ex8_8.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+h = 120;// feet

+d = 5;// feet

+h1 = 49;// feet

+p = 42;// lb. per square foot

+c = 0.6;

+//p = k*sqrt(x)

+k = p/sqrt(h1);

+M = integrate('18*x^(3/2)','x',0,120);

+printf('Bending moment at the foot of the chimney is, M = %d lb-ft',M);

+//there is an error in the answer given in text book

diff --git a/2744/CH8/EX8.9/Ex8_9.sce b/2744/CH8/EX8.9/Ex8_9.sce
new file mode 100755
index 000000000..3a9a5642a
--- /dev/null
+++ b/2744/CH8/EX8.9/Ex8_9.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+h = 100;// feet

+d = 4;// feet

+p = 50;// lb. per square foot

+c = 2/3;

+M = integrate('(100*x/3)*(10-(6*x/100))','x',0,100);

+printf('Bending moment at the foot of the chimney is, M = %d lb-feet',M);