initial commit / add all books

14 files changed, 466 insertions, 0 deletions

diff --git a/1898/CH13/EX13.1/Ex13_1.sce b/1898/CH13/EX13.1/Ex13_1.sce
new file mode 100755
index 000000000..1c67b528d
--- /dev/null
+++ b/1898/CH13/EX13.1/Ex13_1.sce
@@ -0,0 +1,25 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.1 : ")

+

+//Given:

+l = 7.2*1000; //mm

+E = 200; //GPa

+ro = 75; //mm

+ri = 70; //mm

+sigma_y = 250; //MPa

+

+//Calculations:

+I = (%pi/4)*(ro^4 - ri^4)

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

+

+Pcr = (%pi^2*(E*10^6)*I*(1000)^-2)/(l^2); //Pcr = (%pi^2*EI)/(l^2)

+

+sigma_cr = (Pcr*1000)/A;

+

+if(sigma_cr<sigma_y)

+   

+    printf("\n\nThe maximum allowable axial load that the column can support    = %1.1f kN',Pcr);

+end

+

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

diff --git a/1898/CH13/EX13.10/Ex13_10.sce b/1898/CH13/EX13.10/Ex13_10.sce
new file mode 100755
index 000000000..0bc807691
--- /dev/null
+++ b/1898/CH13/EX13.10/Ex13_10.sce
@@ -0,0 +1,37 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.10 : ")

+

+//Given:

+L = 750; //mm

+P = 60; //kN

+sigma = 195; //N/mm^2

+K = 1;

+

+//Calculations:

+b2 = (P*1000)/(2*sigma);

+b = sqrt(b2);

+

+A = 2*b*b;

+Iy = (1/12)*(2*b*b^3);

+ry = sqrt(Iy/A);

+

+sl_ratio = (K*L)/(ry);

+

+

+

+if(sl_ratio>12)

+    b4 = (P*1000*2598.1^2)/(2*378125); //Eqn 13.26

+    b = b4^(1/4);

+   

+    sl_ratio = (2598.1)/(b);

+    w = 2*b;

+    

+    if(sl_ratio>55)

+        printf('\n\nThe thickness of the bar  = %1.0fmm',b);

+       printf('\nThe width of the bar      = %1.0fmm',w);

+   end

+end

+

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

+   

diff --git a/1898/CH13/EX13.11/Ex13_11.sce b/1898/CH13/EX13.11/Ex13_11.sce
new file mode 100755
index 000000000..1139abcf7
--- /dev/null
+++ b/1898/CH13/EX13.11/Ex13_11.sce
@@ -0,0 +1,25 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.11 : ")

+

+//Given:

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

+y1 = 150; //mm

+x1 = 40; //mm

+A = (x1*y1);

+d = 40;

+K = 1;

+

+//Eqn 13.29

+

+L2 = (3718*A*d^2)/(P);

+L = sqrt(L2);

+KL_d = (K*L)/(d);

+

+if(KL_d>26 & KL_d<=50)

+     printf('\n\nThe greatest allowable length L as specified by the NFPA  = %1.0f mm',L);

+   

+end 

+

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

+

diff --git a/1898/CH13/EX13.12/Ex13_12.sce b/1898/CH13/EX13.12/Ex13_12.sce
new file mode 100755
index 000000000..4eb3846c8
--- /dev/null
+++ b/1898/CH13/EX13.12/Ex13_12.sce
@@ -0,0 +1,32 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.12 : ")

+

+//Given:

+L = 1600; //mm

+K = 2;

+l = 80; //mm

+b = 40; //mm

+e = 20; //mm

+c = 40; //mm

+

+//Calculations:

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

+A = l*b;

+r = sqrt(I1/A);

+sl_ratio = (K*L)/(r);

+

+//Eqn 13.26:

+sigma_allow = (378125)/(sl_ratio^2);

+

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

+coefficient = (1/A) + (e*c)/I2;

+sigma_max = sigma_allow;

+P = sigma_max/coefficient;

+P = P/1000;

+

+//Display:

+

+printf('\n\nThe load that can be supported if the column is fixed at its base  = %1.2f kN',P);

+

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

diff --git a/1898/CH13/EX13.13/Ex13_13.sce b/1898/CH13/EX13.13/Ex13_13.sce
new file mode 100755
index 000000000..39e1aa598
--- /dev/null
+++ b/1898/CH13/EX13.13/Ex13_13.sce
@@ -0,0 +1,42 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.13 : ")

+

+//Given:

+sigmaB_allow = 160; //MPa

+E = 200; //GPa

+sigma_y = 250; //MPa

+K= 1;

+A = 3790; //mm^2

+Ix = 17.1*10^6; //mm^4

+ry = 38.2; //mm

+d = 157; //mm

+c= d/2; 

+e = 750; //mm

+L = 4000; //mm

+

+sl_ratio = (K*L)/(ry);

+sl_ratio_c = sqrt((2*%pi^2*E*1000)/(sigma_y));

+

+

+

+if(sl_ratio<sl_ratio_c)

+    num = (1 - (sl_ratio^2/(2*sl_ratio_c^2)))*sigma_y;

+    denom1 = (5/3) + ((3/8)*sl_ratio/sl_ratio_c);

+    denom2 = (sl_ratio^3)/(8*sl_ratio_c^3);

+    sigmaA_allow = num/(denom1 - denom2);

+    

+    coeffP = 1/(sigmaA_allow*A) + (e*c)/(Ix*sigmaB_allow);

+    P = 1/coeffP;

+    

+    sigA = (P/A)/(sigmaA_allow);

+    P = P/1000; //in kN

+   

+  

+    if(sigA < 0.15)

+        printf('\n\nThe maximum allowable value of eccentric load  = %1.2f kN',P);

+    end

+end

+

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

+    

diff --git a/1898/CH13/EX13.14/Ex13_14.sce b/1898/CH13/EX13.14/Ex13_14.sce
new file mode 100755
index 000000000..4467cfb1e
--- /dev/null
+++ b/1898/CH13/EX13.14/Ex13_14.sce
@@ -0,0 +1,37 @@
+

+clear all; clc;

+

+disp("Scilab Code Ex 13.14 : ")

+

+//Given:

+K = 2;

+d= 60; //mm

+L = 1200; //mm

+e = 80; //mm

+c = d;

+A = 60*120; //mm^2

+l = 60; //mm

+b = 120;//mm

+

+

+//Calculations:

+sl_ratio = (K*L)/(d);

+

+if(sl_ratio>26 & sl_ratio<50)

+    sigma_allow = (3718)/(sl_ratio^2);

+    sigma_max = sigma_allow;

+    

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

+    coeffP = (1/A) + (e*c)/(I);

+    P = sigma_max/coeffP;

+    P = P/1000; //kN

+    

+     printf('\n\nThe eccentric load that can be supported = %1.2f kN',P);

+end

+

+//Answer given in the textbook varies.

+

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

+    

+    

+    

diff --git a/1898/CH13/EX13.2/Ex13_2.sce b/1898/CH13/EX13.2/Ex13_2.sce
new file mode 100755
index 000000000..12909c52c
--- /dev/null
+++ b/1898/CH13/EX13.2/Ex13_2.sce
@@ -0,0 +1,27 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.2 : ")

+

+//Given:

+E = 200; //GPa

+I = 15.3*10^6; //mm^4

+l= 4*1000; //mm

+A = 5890; //mm^2

+sigma_y = 250; //MPa

+

+//Calculations:

+

+Pcr = ((%pi^2)*E*10^6*I*1000^-2)/(l^2); //Pcr = (%pi^2*EI)/(l^2)

+

+sigma_cr = (Pcr*1000)/A;

+

+if(sigma_cr>sigma_y)

+    Pcr = (sigma_y*A);

+    Pcr = Pcr/1000; //in kN

+end

+

+//Display:

+

+printf("\n\nThe maximum allowable axial load that the column can support    = %1.1f kN',Pcr);

+

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

diff --git a/1898/CH13/EX13.3/Ex13_3.sce b/1898/CH13/EX13.3/Ex13_3.sce
new file mode 100755
index 000000000..f06bc7d5c
--- /dev/null
+++ b/1898/CH13/EX13.3/Ex13_3.sce
@@ -0,0 +1,33 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.3 : ")

+

+//Given:

+E = 200;//GPa

+Ix = 13.4*10^-6;

+Iy = 1.83*10^-6;

+l = 8;

+KLx = 0.5*l; //m

+KLy = 0.7*(l/2); //m

+rx = 66.2; //mm

+ry = 24.5; //mm

+

+Pcrx = (%pi^2*E*10^6*Ix)/(KLx^2); //Pcr = (%pi^2*EI)/(l^2)

+Pcry = (%pi^2*E*10^6*Iy)/(KLy^2); //Pcr = (%pi^2*EI)/(l^2)

+

+Pcr = min(Pcrx,Pcry);

+A = 3060; //mm^2

+sigma_cr = Pcr/A;

+

+sl_ratio_x = (KLx*1000)/(rx);

+sl_ratio_y = (KLy*1000)/(ry);

+s_ratio = max(sl_ratio_x, sl_ratio_y);

+

+//Display:

+

+printf("\n\nThe maximum load that the column can support without buckling    = %1.0f kN',Pcr);

+printf("\nThe largest slenderness ratio                                    = %1.1f N/mm^2',s_ratio);

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

+

+

+

diff --git a/1898/CH13/EX13.4/Ex13_4.sce b/1898/CH13/EX13.4/Ex13_4.sce
new file mode 100755
index 000000000..bf866097e
--- /dev/null
+++ b/1898/CH13/EX13.4/Ex13_4.sce
@@ -0,0 +1,31 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.4 : ")

+

+//Given:

+E = 70;//GPa

+Ix = 61.3*10^-6;

+Iy = 23.2*10^-6;

+l = 5;

+KLx = 2*l; //m

+KLy = 0.7*(l); //m

+FS = 3; //Factor of safety

+sigma_y = 215; //MPa

+

+

+Pcrx = (%pi^2*E*10^6*Ix)/(KLx^2); //Pcr = (%pi^2*EI)/(l^2)

+Pcry = (%pi^2*E*10^6*Iy)/(KLy^2); //Pcr = (%pi^2*EI)/(l^2)

+

+Pcr = min(Pcrx,Pcry);

+A = 7.5*10^-3; //mm^2

+P_allow = Pcr/FS;

+sigma_cr = (Pcr*10^-3)/A;

+

+

+if(sigma_cr<sigma_y)

+

+  printf("\n\nThe largest allowable load that the column can support     = %1.0f kN',P_allow);

+end

+

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

+

diff --git a/1898/CH13/EX13.5/Ex13_5.sce b/1898/CH13/EX13.5/Ex13_5.sce
new file mode 100755
index 000000000..f0b7eadd5
--- /dev/null
+++ b/1898/CH13/EX13.5/Ex13_5.sce
@@ -0,0 +1,38 @@
+clear all; clc;
+
+disp("Scilab Code Ex 13.5 : ")
+
+//Given:
+E = 200*10^3; //MPa
+sigma_y = 250; //MPa
+x1 = 50; //mm
+y1 = 75; //mm
+z1 = 4.5; //m
+e = 25; //mm
+
+Ix = (1/12)*x1*(y1*2)^3;
+A = x1*2*y1;
+rx = sqrt(Ix/A);
+L = z1*1000;
+KL = 1*L;
+
+sl_ratio = KL/rx;
+c = y1;
+ec_r = e*c/(rx^2);
+P_a = 83; //MPa
+A = 7500; //mm^2
+P = P_a*A;
+P = P/1000; //in kN
+
+k = (L/(2*rx))*(sqrt(P/(E*A)));
+sigma_max = (P*1000/A)*(1+ec_r*sec(k)); //Secant Formula
+
+l = sqrt((P*1000)/(E*Ix));
+nu_max = e*(sec(l*L/2)-1);
+
+//Display:
+
+printf('\n\nThe allowable eccentric load that can be applied on the column  = %1.1fkN',P);
+printf('\nThe maximum deflection of the column due to the loading         = %1.0f mm',nu_max);
+
+//--------------------------------------------------------------------------END------------------------------------------------------------------------------------
diff --git a/1898/CH13/EX13.6/Ex13_6.sce b/1898/CH13/EX13.6/Ex13_6.sce
new file mode 100755
index 000000000..001e38a67
--- /dev/null
+++ b/1898/CH13/EX13.6/Ex13_6.sce
@@ -0,0 +1,36 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.6 : ")

+

+//Given:

+z1 = 4*1000; //mm

+e = 200; //mm

+KLy = 0.7*z1;

+Iy = 20.4*10^6;

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

+sigma_y =250; //MPa

+

+//y-y Axis Buckling:

+Pcry = (%pi^2*E*10^6*Iy)/(KLy^2); //Pcr = (%pi^2*EI)/(l^2)

+Pcry = Pcry/1000;

+

+//x-x Axis Yielding:

+Kx= 2;

+KLx = Kx*z1;

+c = (z1-KLy)/2;

+rx = 89.9;

+

+//Solved by applying the Secant Formula and then finding Px by trial and error:

+

+trial_Px = 419.4; //kN

+

+A = 7850;//mm^2

+sigma = (trial_Px*1000)/(A);

+

+if(Pcry>trial_Px & sigma<sigma_y)

+printf('\n\nThe maximum eccentric load that the column can support = %1.1fkN',trial_Px);

+printf('\nFailure will occur about the x-x axis.');

+

+end

+

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

diff --git a/1898/CH13/EX13.7/Ex13_7.sce b/1898/CH13/EX13.7/Ex13_7.sce
new file mode 100755
index 000000000..2ee337dac
--- /dev/null
+++ b/1898/CH13/EX13.7/Ex13_7.sce
@@ -0,0 +1,39 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.7 : ")

+

+//Given:

+d = 30; //mm

+r = d/2;

+L = 600; //mm

+sigma_pl = 150;//MPa

+

+//Calculations:

+I = (%pi/4)*(r^4);

+A = %pi*r^2;

+r_gyr = sqrt(I/A);

+K = 1;

+sl_ratio = (K*L)/(r_gyr);

+flag1 = 0;

+

+//Assuming the critical stress is elastic:

+E = 150/0.001;

+sigma_cr1 = (%pi^2*E)/(sl_ratio^2); //Pcr = (%pi^2*EI)/(l^2)

+

+

+if(sigma_cr1 > sigma_pl)

+    Et = (270 - 150)/(0.002 - 0.001);

+    sigma_cr2 = (%pi^2*Et)/(sl_ratio^2); //Pcr = (%pi^2*EI)/(l^2)

+    

+   if(sigma_cr2>150 & sigma_cr2<270) 

+       Pcr = sigma_cr2*A;

+       Pcr = Pcr/1000; //in kN

+        printf('\n\nThe critical load when used as a pin supported column = %1.0fkN',Pcr);

+   

+   end

+    

+       

+end

+

+

+

diff --git a/1898/CH13/EX13.8/Ex13_8.sce b/1898/CH13/EX13.8/Ex13_8.sce
new file mode 100755
index 000000000..73ccfc378
--- /dev/null
+++ b/1898/CH13/EX13.8/Ex13_8.sce
@@ -0,0 +1,30 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.8 : ")

+

+//Given:

+E = 200*10^3; //MPa

+sigma_y = 250; //MPa

+L = 5*1000; //mm

+K = 1;

+A = 19000; //mm^2

+rx = 117; //mm

+ry = 67.4; //mm

+

+//Calculations:

+sl_ratio = (K*L)/(ry);

+sl_ratio_c = sqrt((2*%pi^2*E)/(sigma_y));

+

+if(sl_ratio>0 & sl_ratio<sl_ratio_c)

+    num = (1 - (sl_ratio^2/(2*sl_ratio_c^2)))*sigma_y;

+    denom1 = (5/3) + ((3/8)*sl_ratio/sl_ratio_c);

+    denom2 = (sl_ratio^3)/(8*sl_ratio_c^3);

+    sigma_allow = num/(denom1 - denom2);

+    

+    P = sigma_allow*A;

+    P = P/1000;

+       printf('\n\nThe largest load the pin supported column can safely bear = %1.0f kN',P);

+       

+   end

+   

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

diff --git a/1898/CH13/EX13.9/Ex13_9.sce b/1898/CH13/EX13.9/Ex13_9.sce
new file mode 100755
index 000000000..17411e206
--- /dev/null
+++ b/1898/CH13/EX13.9/Ex13_9.sce
@@ -0,0 +1,34 @@
+clear all; clc;

+

+disp("Scilab Code Ex 13.9 : ")

+

+//Given:

+P = 80; //kN

+E = 210*10^3; //MPa

+sigma_y = 360; //MPa

+L = 5000; //mm

+K = 0.5;

+

+//Calculations:

+I_by_d = (1/4)*(%pi)*(0.5^4);

+A_by_d = (1/4)*(%pi);

+r_by_d = sqrt(I_by_d/A_by_d);

+

+sl_ratio_c = sqrt((2*%pi^2*E)/(sigma_y));

+sigma_allow = (P*1000)/A_by_d;

+

+d4 = (sigma_allow*23*(K*L)^2*16)/(12*%pi^2*E);

+d = d4^(1/4);

+

+//Check:

+d = ceil(d);

+r = d/4;

+KL_r = (K*L)/r;

+

+

+if(KL_r>sl_ratio_c & KL_r<200)

+        printf('\n\nThe smallest diameter of the rod as allowed by AISC specification  = %1.0fmm',d);

+   

+end 

+

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