initial commit / add all books

7 files changed, 275 insertions, 0 deletions

diff --git a/1898/CH7/EX7.1/Ex7_1.sce b/1898/CH7/EX7.1/Ex7_1.sce
new file mode 100755
index 000000000..921a97e4d
--- /dev/null
+++ b/1898/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,40 @@
+clear all; clc;

+

+disp("Scilab Code Ex 7.1 : ")

+

+//Given:

+V = 3; //kN

+h = 125; //mm

+b = 100; //mm

+y_top = 50; //mm

+x_right = 37.5; //mm

+

+//Part (a):

+

+//Section Properties:

+I = (b*h^3)/12;

+y_dash_1 = ((h-y_top)-(h/2));

+A = y_top*b;

+Q = (y_dash_1+(y_top/2))*A;

+

+//Shear Stress:

+tou_p = (V*Q)/(I*b); //tou = VQ/It

+tou_p = tou_p*10^3;

+

+//Part (b):

+

+//Section Properties:

+y_dash_2 = (y_dash_1+(y_top));

+a_dash= b*y_dash_2;

+Q_dash =(y_dash_2*a_dash)/2;

+

+//Shear Stress:

+tou_max = (V*Q_dash)/(I*b);

+tou_max = tou_max*10^3;

+

+//Display:

+

+printf("\n\nThe shear stress in the beam at point P    = %1.3f MPa',tou_p);

+printf('\nThe maximum shear stress in the beam       = %1.3f MPa',tou_max);

+

+//----------------------------------------------------------------------END--------------------------------------------------------------------------------

diff --git a/1898/CH7/EX7.2/Ex7_2.sce b/1898/CH7/EX7.2/Ex7_2.sce
new file mode 100755
index 000000000..8b624ffb1
--- /dev/null
+++ b/1898/CH7/EX7.2/Ex7_2.sce
@@ -0,0 +1,50 @@
+clear all; clc;

+

+disp("Scilab Code Ex 7.2 : ")

+

+//Given:

+V = 80; //kN

+thick_1 = 20/1000; //m

+thick_2 = 15/1000; //m

+l = 300/1000; //m

+y = 100/1000; //m

+h = 2*y;

+y_dash = y +thick_1/2;

+

+//Part(a):

+

+I1 = (thick_2*(h^3))/12;

+I2 = (l*(thick_1^3))/12;

+I3 = (l*thick_1*(y_dash)^2);

+I = I1+2*(I2+I3); //Moment of inertia

+

+Q_b = y_dash*l*thick_1;

+//At B'

+tou_b_dash = (V*Q_b)/(I*l*1000);

+//At B

+tou_b = (V*Q_b)/(I*thick_2*1000);

+

+//At C:

+Q_c = (y_dash*l*thick_1)+(y*thick_2*y/2);

+tou_c = (V*Q_c)/(I*thick_2*1000);

+

+//Part(b)

+

+

+y0 = -0.1;

+y1 = 0.1;

+

+function Q =f(y),Q = ((0.735 - (7.5*y*y))*10^-3),

+endfunction

+Int =intg(y0,y1,f)

+

+V_w = (V*Int*thick_2)/(I*thick_2);

+

+//Display:

+

+printf("\n\nThe shear stress at B dash             = %1.2f MPa',tou_b_dash);

+printf("\nThe shear stress at B                  = %1.1f MPa',tou_b);

+printf("\nThe shear stress at C                  = %1.1f MPa',tou_c);

+printf("\nThe shear force resisted by the web    = %1.1f kN',V_w);

+

+//------------------------------------------------------------------------END----------------------------------------------------------------------------------

diff --git a/1898/CH7/EX7.3/Ex7_3.sce b/1898/CH7/EX7.3/Ex7_3.sce
new file mode 100755
index 000000000..f236accae
--- /dev/null
+++ b/1898/CH7/EX7.3/Ex7_3.sce
@@ -0,0 +1,39 @@
+clear all; clc;

+

+disp("Scilab Code Ex 7.3 : ")

+

+//Given:

+udl = 6.5; //kN

+l_bc = 8; //m

+l = 150/1000;//m

+t = 30/1000;//m

+

+//Internal Shear:

+w = udl*l_bc/2;

+l_wc = l_bc/4;

+l_bw = l_bc - l_wc;

+V = (w*l_bw)/l_bc;

+R_b = w - V;

+

+//Section Properties:

+y1= l/2;

+A = (l*t);

+y2= l+(t/2);

+y_dash = (y1*A + y2*A)/(2*A);

+I1 = (t*l^3)/12;

+I2 = (A*(y_dash-y1)^2);

+I3 = (l*t^3)/12;

+I4 = (A*(y2 - y_dash)^2);

+I = I1+I2+I3+I4;

+

+Q = ((l+t)-(t/2)-y_dash)*A;

+

+//Shear Stress:

+tou_max = (V*Q)/(I*t*1000);

+

+//Display:

+

+printf("\n\nThe maximum shear stress in the glue necessary to hold the boards together    = %1.2f MPa',tou_max);

+

+

+//----------------------------------------------------------------------END--------------------------------------------------------------------------------

diff --git a/1898/CH7/EX7.4/Ex7_4.sce b/1898/CH7/EX7.4/Ex7_4.sce
new file mode 100755
index 000000000..73c377ca9
--- /dev/null
+++ b/1898/CH7/EX7.4/Ex7_4.sce
@@ -0,0 +1,48 @@
+clear all; clc;

+

+disp("Scilab Code Ex 7.4 : ")

+

+//Given:

+

+V = 850; //kN

+l1 =250/1000; //m

+l2 = 300/1000; //m

+l3 = 125/1000;//m

+t = 10/1000; //m

+h = 200/1000; //m

+

+A1 = l1*t;

+A2 = l2*t;

+A3 = l3*t;

+

+y1 = l2+(t/2);

+y2 = l2/2;

+y3 = h+(t/2);

+

+y_dash = (2*y2*A2 + A1*y1 + A3*y3)/(2*A2 + A1 + A3);

+

+I1 = ((l1*t^3)/12) +(A1 * (l2+(t/2)-y_dash)^2);

+I2 = ((t*l2^3)/12) +(A2 * (y_dash - (l2/2))^2);

+I3 = ((l3*t^3)/12) +(A1 * (h+(t/2)-y_dash)^2);

+I = 2*I2 + I1 + I3;

+

+Q_b = (l2+(t/2) - y_dash)*A1; //Q = y'A'

+Q_c = (h+(t/2) - y_dash)*A3; //Q = y'A'

+

+//Shear Flow:

+

+q_b = (V*Q_b)/I;

+q_c = (V*Q_c)/I;

+

+q_b = q_b/(2*1000);

+q_c = q_c/(2*1000);

+

+//Display:

+

+printf("\n\nThe shear flow at B, resisted by the glue is    = %1.2f MN/m',q_b);

+printf("\nThe shear flow at C, resisted by the glue is    = %1.4f MN/m',q_c);

+

+

+

+

+//----------------------------------------------------------------------END--------------------------------------------------------------------------------

diff --git a/1898/CH7/EX7.5/Ex7_5.sce b/1898/CH7/EX7.5/Ex7_5.sce
new file mode 100755
index 000000000..de724633b
--- /dev/null
+++ b/1898/CH7/EX7.5/Ex7_5.sce
@@ -0,0 +1,32 @@
+clear all; clc;

+

+disp("Scilab Code Ex 7.5 : ")

+

+//Given:

+V = 80; //N

+t = 1.5; //cm

+a = 7.5; //cm

+b = a-2*t; //cm

+F_nail= 30; //N

+

+//Section Properties:

+I = (a*a^3 - b*b^3 )/12;

+Q_b = (((a-2*t)/2)+(t/2))*a*t; //Q = y'A'

+Q_c = (((a-2*t)/2)+(t/2))*(a-2*t)*t; //Q = y'A'

+

+//Shear Flow:

+q_b = (V*Q_b)/I;

+q_c = (V*Q_c)/I;

+

+s_b = F_nail/(q_b/2);

+s_c = F_nail/(q_c/2);

+

+//Display:

+

+

+printf("\n\nThe maximum spacing of nails required at B is    = %1.0f cm',s_b);

+printf("\nThe maximum spacing of nails required at C is    = %1.1f cm',s_c);

+

+

+

+//----------------------------------------------------------------------END--------------------------------------------------------------------------------

diff --git a/1898/CH7/EX7.6/Ex7_6.sce b/1898/CH7/EX7.6/Ex7_6.sce
new file mode 100755
index 000000000..db891753f
--- /dev/null
+++ b/1898/CH7/EX7.6/Ex7_6.sce
@@ -0,0 +1,33 @@
+clear all; clc;

+

+disp("Scilab Code Ex 7.6 : ")

+

+//Given:

+F = 40; //N

+s = 9; //cm

+h = 5; //cm

+t = 0.5; //cm

+w = 3; //cm

+w_3 = w/3; //cm

+

+//Calculations:

+

+I = (w*h^3)/12 - (2*w_3*(h - 2*t)^3)/12;

+

+//Case 1:

+

+Q1 = ((h-t)/2)*(w*t);

+V1 =((F/s)*I)/Q1 ;  //q = VQ/I

+

+//Case2:

+

+Q2 = ((h-t)/2)*(w_3*t);

+V2 =((F/s)*I)/Q2 ;  //q = VQ/I

+

+//Display:

+

+

+printf("\n\nThe largest vertical shear that can be supported in Case 1    = %1.1f N',V1);

+printf("\nThe largest vertical shear that can be supported in Case 2    = %1.1f N',V2);

+

+//-------------------------------------------------------------------------END---------------------------------------------------------------------------------------

diff --git a/1898/CH7/EX7.7/Ex7_7.sce b/1898/CH7/EX7.7/Ex7_7.sce
new file mode 100755
index 000000000..6e752ac58
--- /dev/null
+++ b/1898/CH7/EX7.7/Ex7_7.sce
@@ -0,0 +1,33 @@
+clear all; clc;

+

+disp("Scilab Code Ex 7.7 : ")

+

+//Given:

+V = 10; //kN

+b1 = 6; //cm

+h1 = 8; //cm

+t = 1; //cm

+b2 = b1-2*t;

+h2 = h1-2*t; //cm

+b3 = 4; //cm

+

+//Calculations:

+I = ((b1*h1^3)/12) - ((b2*h2^3)/12);

+

+q_b = 0;

+

+Q_c = ((b1/2)+(t/2))*(b3+(t))*t; 

+q_c = (V*Q_c*100)/(I); //Q = VQ/I

+

+Q_d = (2*h1/4*t*b3) + ((b1/2)+(t/2))*b3*t;

+q_d = (V*Q_d*100)/(I); //Q = VQ/I

+

+//Display:

+

+

+printf("\n\nVariation of shear flow at B  = %1.1f N/mm',q_b);

+printf("\nVariation of shear flow at C  = %1.1f N/mm',q_c);

+printf('\nVariation of shear flow at D  = %1.1f N/mm',q_d);

+

+//-------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+