initial commit / add all books

13 files changed, 427 insertions, 0 deletions

diff --git a/1898/CH14/EX14.1/Ex14_1.sce b/1898/CH14/EX14.1/Ex14_1.sce
new file mode 100755
index 000000000..9534b5795
--- /dev/null
+++ b/1898/CH14/EX14.1/Ex14_1.sce
@@ -0,0 +1,35 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.1 : ")

+

+//Given:

+sigma_y = 310; //N/mm^2

+db =18; //mm

+rb = db/2;

+Ab = %pi*(rb^2);

+E = 210*10^3; //N/mm^2

+da1 = 20; //mm

+ra1 = da1/2;

+Aa1 = %pi*(ra1^2);

+La1 = 50;//mm

+La2= 6; //mm

+da2 =18; //mm

+ra2 = da2/2;

+Aa2 = %pi*(ra2^2);

+Lb = 56; //mm

+

+

+//Bolt A:

+P_max = sigma_y*Ab;

+Uia = (P_max^2/(2*E))*(La1/Aa1 + La2/Aa2); //Ui = (N^2L)/(2AE)

+Uia = Uia/1000;

+

+//Bolt B:

+Uib = (P_max^2/(2*E))*(Lb/Ab);

+Uib = Uib/1000;

+

+//Display:

+    printf('\n\nThe greatest amount of strain energy absorbed by bolt A    = %1.3f J',Uia);

+    printf('\nThe greatest amount of strain energy absorbed by bolt B    = %1.3f J',Uib);

+    

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

diff --git a/1898/CH14/EX14.10/Ex14_10.sce b/1898/CH14/EX14.10/Ex14_10.sce
new file mode 100755
index 000000000..bd4ca6544
--- /dev/null
+++ b/1898/CH14/EX14.10/Ex14_10.sce
@@ -0,0 +1,30 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.10 : ")

+

+//Given:

+m = 80*1000; //kg

+v = 0.2; //m/s

+l_ac = 1.5; //m

+E = 200*10^9; //N/m^2

+w = 0.2; //m

+I = (1/12)*(w^4);

+l_ab = 1000; //mm

+

+//Calculations:

+del_Amax = sqrt((m*v^2*l_ac^3)/(3*E*I));

+

+P_max = (3*E*I*del_Amax)/(l_ac^3);

+theta_A = (P_max*l_ac^2)/(2*E*I);

+del_Amax = del_Amax*1000;

+del_Bmax = del_Amax + (theta_A*l_ab);

+

+

+//Display:

+

+    printf('\n\nThe maximum horizontal displacement of the post at B due to impact   = %1.1f mm',del_Bmax);

+    

+

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

+

+

diff --git a/1898/CH14/EX14.11/Ex14_11.sce b/1898/CH14/EX14.11/Ex14_11.sce
new file mode 100755
index 000000000..1ddc41e5f
--- /dev/null
+++ b/1898/CH14/EX14.11/Ex14_11.sce
@@ -0,0 +1,30 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.11 : ")

+

+//Given:

+A = 400*10^-6; //m^2

+E = 200*10^6; //kN/m^2

+P = 100; //kN

+

+//Virtual Work Equation:

+

+n = [0 0 -1.414 1];

+N = [-100 141.4 -141.4 200];

+L = [4 2.828 2.828 2];

+del_cv = 0;

+

+for i=1:4

+    del_cv = del_cv + (n(i)*N(i)*L(i))/(A*E);

+end

+

+del_cv = del_cv*1000;

+

+//Display:

+

+    printf('\n\nThe vertical displacement of joint C of the steel truss   = %1.1f mm',del_cv);

+    

+

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

+

+

diff --git a/1898/CH14/EX14.12/Ex14_12.sce b/1898/CH14/EX14.12/Ex14_12.sce
new file mode 100755
index 000000000..3ac25a7cb
--- /dev/null
+++ b/1898/CH14/EX14.12/Ex14_12.sce
@@ -0,0 +1,43 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.12 : ")

+

+//Given:

+A = 300*10^-6; //m^2

+E = 210*10^6; //kN/m^2

+P = 60; //kN

+F_ac = 1.25; //kN

+

+//Part a:

+

+//Virtual Work Equation:

+

+n = [0 1.25 0 -0.75];

+N = [0 75 -60 -45];

+L = [1.5 2.5 2 1.5];

+del_ch = 0;

+

+for i=1:4

+    del_ch = del_ch + (n(i)*N(i)*L(i))/(A*E);

+end

+

+del_chA = del_ch*1000;

+

+//Part b:

+

+del_L = -6; //mm

+del_chB = F_ac*del_L;

+

+if(del_chB<0)

+    

+

+//Display:

+   

+    printf('\n\nThe horizontal displacement of joint C if a force is applied to B   = %1.3f mm',del_chA);

+    printf('\nThe horizontal displacement of joint C if AC is fabricated short    = %1.1f mm',del_chB);

+end

+

+    

+

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

+

diff --git a/1898/CH14/EX14.13/Ex14_13.sce b/1898/CH14/EX14.13/Ex14_13.sce
new file mode 100755
index 000000000..dbacc100f
--- /dev/null
+++ b/1898/CH14/EX14.13/Ex14_13.sce
@@ -0,0 +1,23 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.13 : ")

+

+//Given:

+del_T = 60; //degree celcius

+alpha = 12*10^-6; //per degree celcius

+E = 200*10^6; //kN/m^2

+A = 250*10^-6; //m^2

+L = 4; //m

+

+//Virtual Work Equation:

+n = 1.155; //kN

+N = -12; //kN

+

+del_bh = (n*N*L)/(A*E) + (n*alpha*del_T*L);

+del_bh = del_bh*1000;

+

+//Display:

+

+printf('\n\nThe horizontal displacement of joint B of the truss   = %1.2f mm',del_bh);

+

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

diff --git a/1898/CH14/EX14.16/Ex14_16.sce b/1898/CH14/EX14.16/Ex14_16.sce
new file mode 100755
index 000000000..a5ae25568
--- /dev/null
+++ b/1898/CH14/EX14.16/Ex14_16.sce
@@ -0,0 +1,41 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.16 : ")

+

+//Given:

+I = 175.8*10^-6; //m^4

+E = 200*10^6; //kN/m^2

+Ra = 1;//kN

+l_ab = 3; //m

+l_bc = 6; //m

+

+

+//Virtual Work Equation:

+m1 = -1; //*x1

+M1 = -2.5; //*x1^3

+m2 = -0.5; //*x2

+

+x10 = 0;

+x11 = l_ab;

+I1 = integrate('m1*M1*(x1^4)','x1',x10,x11);

+

+x20 = 0;

+x21 = l_bc;

+I2 = integrate('m2*123.75*(x2^2)','x2',x20,x21);

+

+x20 = 0;

+x21 = l_bc;

+I3 = integrate(' -m2*22.5*(x2^3)','x2',x20,x21);

+

+In = I1 + I2 + I3;

+del_A = (In)/(E*I);

+del_A = del_A*1000;

+

+

+//Display:

+

+    printf('\n\nThe displacement of point A of the steel beam   = %1.1f mm',del_A);

+    

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

+

+

diff --git a/1898/CH14/EX14.17/Ex14_17.sce b/1898/CH14/EX14.17/Ex14_17.sce
new file mode 100755
index 000000000..3b9396fb1
--- /dev/null
+++ b/1898/CH14/EX14.17/Ex14_17.sce
@@ -0,0 +1,37 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.17 : ")

+

+//Given:

+E = 210*10^3; //N/mm^2

+P = 40*10^3;//N

+A_ab = 1250; //mm^2

+A_ac = 625; //mm^2

+A_cd = 1250; //mm^2

+A_bc = 625; //mm^2

+

+N_by_P = [0 0 1.67 -1.33];

+L = [4000 3000 5000 4000];

+A = [A_ab A_bc A_ac A_cd];

+N = zeros(4);

+sum = 0;

+

+

+for i =1:4

+    N(i) = N_by_P(i)*P;

+    num(i) = N(i)*N_by_P(i)*L(i);

+   

+end

+

+for i = 1:4

+    sum = sum + (num(i)/(A(i)*E)); //By Castigliano's Second theorem.

+end

+

+del_ch = sum;

+

+//Display:

+    printf('\n\nThe horizontal displacement of joint C of the steel truss   = %1.2f mm',sum);

+    

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

+

+

diff --git a/1898/CH14/EX14.18/Ex14_18.sce b/1898/CH14/EX14.18/Ex14_18.sce
new file mode 100755
index 000000000..065bafbce
--- /dev/null
+++ b/1898/CH14/EX14.18/Ex14_18.sce
@@ -0,0 +1,31 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.18 : ")

+

+//Given:

+E = 200*10^6; //kN/m^2

+P = 0;//N

+A = 400*10^-6; //m^2

+

+N_by_P = [0 0 -1.414 1];

+L = [4 2.828 2.828 2];

+N = [-100 141.4 -141.4 200];

+sum = 0;

+

+

+for i =1:4

+    num(i) = N(i)*N_by_P(i)*L(i); 

+end

+

+for i = 1:4

+    sum = sum + (num(i)/(A*E)); //By Castigliano's Second theorem.

+end

+

+del_ch = sum*1000;

+

+//Display:

+  printf('\n\nThe vertical displacement of joint C of the steel truss   = %1.1f mm',del_ch);

+    

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

+

+

diff --git a/1898/CH14/EX14.21/Ex14_21.sce b/1898/CH14/EX14.21/Ex14_21.sce
new file mode 100755
index 000000000..f7ab470e6
--- /dev/null
+++ b/1898/CH14/EX14.21/Ex14_21.sce
@@ -0,0 +1,39 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.21 : ")

+

+//Given:

+I = 125*10^-6; //m^4

+E = 200*10^6; //kN/m^2

+Rc = 5;//kN

+l_ac = 6; //m

+l_cb = 4; //m

+

+

+//Castigliano's Second Theorem:

+m = 0.4/9;

+

+x10 = 0;

+x11 = l_ac;

+I11 = integrate('4.4*(x1^2)','x1',x10,x11);

+I12 = integrate('-m*(x1^4)','x1',x10,x11);

+I1 = I11 + I12;

+

+x20 = 0;

+x21 = l_cb;

+I21 = integrate('6*0.6*(x2^2)','x2',x20,x21);

+I22 = integrate('18*0.6*(x2)','x2',x20,x21);

+I2 = I21+I22;

+

+In = I1 + I2 ;

+del_cv = (In)/(E*I);

+del_cv = del_cv*1000;

+

+

+//Display:

+

+    printf('\n\nThe vertical displacement of point C of the steel beam   = %1.1f mm',del_cv);

+    

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

+

+

diff --git a/1898/CH14/EX14.5/Ex14_5.sce b/1898/CH14/EX14.5/Ex14_5.sce
new file mode 100755
index 000000000..6dd06de2e
--- /dev/null
+++ b/1898/CH14/EX14.5/Ex14_5.sce
@@ -0,0 +1,30 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.5 : ")

+

+//Given:

+G = 75*10^9; //N/m^2

+ro = 80/1000; //m

+t = 15/1000; //m

+ri = ro - t;

+l1 = 750/1000; //m

+l2 = 300/1000; //m

+T1 = 40; //Nm

+T2 =15; //Nm

+

+//Calculations:

+

+J = (%pi/2)*(ro^4 - ri^4);

+

+//Eqn 14-22

+U1 = (T1^2*l1)/(2*G*J); 

+U2 = (T2^2*l2)/(2*G*J);

+Ui = U1 + U2;

+Ui = Ui*10^6; //in micro Joule

+

+//Display:

+

+   printf('\n\nThe strain energy stored in the shaft    = %1.0fX10^-6 J',Ui);

+    

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

+    

diff --git a/1898/CH14/EX14.6/Ex14_6.sce b/1898/CH14/EX14.6/Ex14_6.sce
new file mode 100755
index 000000000..26ad594d5
--- /dev/null
+++ b/1898/CH14/EX14.6/Ex14_6.sce
@@ -0,0 +1,36 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.6 : ")

+

+//Given:

+l_ab = 1; //m

+l_bc = 2; //m

+N_ab = 11.547*1000; //N

+Nb = 20*1000; //N

+Nc = -23.094*1000; //N

+N_ac = -20*1000; //N

+A = 100/(1000^2); //mm^2

+E = 200*10^9; //N/m^2

+P = 20*10^3;//N

+

+//Eqn 14-26

+P_by_2 = P/2;

+l_ac = sqrt(l_bc^2 - l_ab^2);

+del = 0;

+

+N2= [N_ab^2 Nc^2 N_ac^2];

+L = [l_ab l_bc l_ac];

+

+for i = 1:3

+    del = del + (N2(i)*L(i))/(2*A*E);

+end

+

+del_bh = del/P_by_2;

+del_bh = del_bh*1000;

+

+//Display:

+

+printf('\n\nThe horizontal displacement at point B   = %1.2fmm',del_bh);

+    

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

+    

diff --git a/1898/CH14/EX14.8/Ex14_8.sce b/1898/CH14/EX14.8/Ex14_8.sce
new file mode 100755
index 000000000..fff23bee6
--- /dev/null
+++ b/1898/CH14/EX14.8/Ex14_8.sce
@@ -0,0 +1,27 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.8 : ")

+

+//Given:

+ro = 60; //mm

+ri = 50; //mm

+E = 70; //kN/mm^2

+W = 600; //kN

+L = 240; //mm

+h = 0;

+

+//Part a:

+

+A = (%pi)*(ro^2 - ri^2);

+del_st= (W*L)/(A*E);

+

+//Part b:

+

+del_max = del_st*(1 + sqrt(1 + 2*(h/del_st)));

+

+//Display:

+

+  printf('\n\nThe maximum displacement at the top of the pipe for gradually applied load   = %1.4f mm',del_st);

+    printf('\nThe maximum displacement at the top of the pipe for suddenly applied load    = %1.4f mm',del_max);

+

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

diff --git a/1898/CH14/EX14.9/Ex14_9.sce b/1898/CH14/EX14.9/Ex14_9.sce
new file mode 100755
index 000000000..ec0583916
--- /dev/null
+++ b/1898/CH14/EX14.9/Ex14_9.sce
@@ -0,0 +1,25 @@
+clear all; clc;

+

+disp("Scilab Code Ex 14.9 : ")

+

+//Given:

+W = 6000; //N

+h = 50; //mm

+E = 210*1000; //N/mm^2

+L = 5000; //mm

+I = 87.3*10^6; //mm^2

+

+//Calculations:

+

+del_st = (W*L^3)/(48*E*I);

+del_max = del_st*(1 + sqrt(1 + 2*(h/del_st)));

+

+c = 252/2;

+sigma_max = (12*E*del_max*c)/(L^2);

+

+//Display:

+

+     printf('\n\nThe maximum bending stress in the steel beam   = %1.2f N/mm^2',sigma_max);

+    printf('\nThe maximum deflection in the beam             = %1.3f mm',del_max);

+

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