initial commit / add all books

14 files changed, 494 insertions, 0 deletions

diff --git a/1898/CH10/EX10.1/Ex10_1.sce b/1898/CH10/EX10.1/Ex10_1.sce
new file mode 100755
index 000000000..2b5ba7209
--- /dev/null
+++ b/1898/CH10/EX10.1/Ex10_1.sce
@@ -0,0 +1,29 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.1 : ")

+

+//Given:

+ep_x = 500; //Normal Strain

+ep_y = -300; //Normal Strain

+gamma_xy = 200; //Shear Strain

+theta = 30*(%pi/180);

+theta = theta*-1;

+

+ep_x_new = ((ep_x+ep_y)/2) + ((ep_x - ep_y)/2)*cos(2*theta) + (gamma_xy/2)*sin(2*theta);

+

+gamma_xy_new = -((ep_x - ep_y)/2)*sin(2*theta) + (gamma_xy/2)*cos(2*theta);

+gamma_xy_new = 2*gamma_xy_new;

+

+phi = 60*(%pi/180);

+ep_y_new = (ep_x+ep_y)/2 + ((ep_x - ep_y)/2)*cos(2*phi) + (gamma_xy/2)*sin(2*phi);

+

+//Display:

+

+

+printf('\n\nThe equivalent strain acting on the element in the x plain oriented at 30 degrees clockwise    = %1.1f *10^-6',ep_x_new);

+printf('\nThe equivalent strain acting on the element in the y plain oriented at 30 degrees clockwise    = %1.1f *10^-6',ep_y_new);

+printf('\nThe equivalent shear strain acting on the element                                              = %1.0f *10^-6',gamma_xy_new);

+

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

+

+

diff --git a/1898/CH10/EX10.10/Ex10_10.sce b/1898/CH10/EX10.10/Ex10_10.sce
new file mode 100755
index 000000000..05a8d3166
--- /dev/null
+++ b/1898/CH10/EX10.10/Ex10_10.sce
@@ -0,0 +1,38 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.10 : ")

+

+//Given:

+a = 300; //mm

+b = 50; //mm

+t = 20; //mm

+E_cu = 120*10^3; //MPa

+nu_cu = 0.34;// Poisson's ratio

+

+//By inspection:

+sigma_x = 800; //MPa

+sigma_y = -500; //MPa

+tou_xy = 0;

+sigma_z = 0;

+

+//By Hooke's Law:

+ep_x = (sigma_x/E_cu) - (nu_cu/E_cu)*(sigma_y + sigma_z);

+ep_y = (sigma_y/E_cu) - (nu_cu/E_cu)*(sigma_x + sigma_z);

+ep_z = (sigma_z/E_cu) - (nu_cu/E_cu)*(sigma_y + sigma_x);

+

+//New lengths:

+

+a_dash = a + ep_x*a;

+b_dash = b + ep_y*b;

+t_dash = t + ep_z*t;

+

+//Display:

+

+printf('\n\nThe new length     = %1.2fmm ',a_dash);

+printf('\nThe new width      = %1.2f mm ',b_dash);

+printf('\nThe new thickness  = %1.2f mm ',t_dash);

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

+

+

+

+

diff --git a/1898/CH10/EX10.11/Ex10_11.sce b/1898/CH10/EX10.11/Ex10_11.sce
new file mode 100755
index 000000000..19efbf42c
--- /dev/null
+++ b/1898/CH10/EX10.11/Ex10_11.sce
@@ -0,0 +1,34 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.11 : ")

+

+//Given:

+p = 20; //kPa

+E = 600; //kPa

+nu = 0.45

+a = 4; //cm

+b = 2; //cm

+c = 3; //cm

+

+//Dilatation:

+sigma_x = -p;

+sigma_y = -p;

+sigma_z = -p;

+

+e = ((1-2*nu)/E)*(sigma_x + sigma_y + sigma_z);

+

+//Change in Length:

+ep = (sigma_x - nu*(sigma_y + sigma_z))/E;

+

+del_a = ep*a;

+del_b = ep*b;

+del_c = ep*c;

+

+//Display:

+

+printf('\n\nThe change in length a      = %1.4fcm ',del_a);

+printf('\nThe change in length b      = %1.5fcm ',del_b);

+printf('\nThe change in length c      = %1.4fcm ',del_c);

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

+

+

diff --git a/1898/CH10/EX10.12/Ex10_12.sce b/1898/CH10/EX10.12/Ex10_12.sce
new file mode 100755
index 000000000..00c449e4f
--- /dev/null
+++ b/1898/CH10/EX10.12/Ex10_12.sce
@@ -0,0 +1,34 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.12 : ")

+

+//Given:

+di = 60/1000; //m

+ri = di/2;

+d0 = 80/1000;//m

+ro = d0/2;

+T = 8000;//Nm

+M = 3500; //Nm

+sigma_y_sqr = 250^2; //MPa

+

+//Calculations:

+c = ro;

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

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

+tou_a = (T*c)/J;

+sigma_a = (M*c)/I;

+

+sigma_avg = (0-sigma_a)/2;

+

+R = sqrt(116.4^2 + sigma_avg^2);

+sigma1 = sigma_avg + R;

+sigma2 = sigma_avg - R;

+

+test = (sigma1^2 - (sigma1*sigma2) + sigma2^2);

+

+

+if(test<sigma_y_sqr)

+    printf("\n\nThe material within the pipe will not yield.");

+end

+

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

diff --git a/1898/CH10/EX10.13/Ex10_13.sce b/1898/CH10/EX10.13/Ex10_13.sce
new file mode 100755
index 000000000..2d4119a06
--- /dev/null
+++ b/1898/CH10/EX10.13/Ex10_13.sce
@@ -0,0 +1,22 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.13 : ")

+

+//Given:

+T = 400; //Nm

+sigma_ult = 150*10^6; //N/m^2

+

+//Calculations:

+

+x = T/(%pi/2);

+r_3 = [x/sigma_ult];

+r = nthroot(r_3, 3);

+r= r*1000; //in mm

+

+//Display:

+

+printf('\n\nThe smallest radius of the solid cast iron shaft     = %1.2fmm ',r);

+

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

+

+

diff --git a/1898/CH10/EX10.14/Ex10_14.sce b/1898/CH10/EX10.14/Ex10_14.sce
new file mode 100755
index 000000000..37aabac96
--- /dev/null
+++ b/1898/CH10/EX10.14/Ex10_14.sce
@@ -0,0 +1,42 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.14 : ")

+

+//Given:

+r = 0.5; //cm

+sigma_yield = 360; //MPa

+T = 3.25; //kN/cm

+A= (%pi*r^2);

+P = 15; //kN

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

+sigma_y_sqr = sigma_yield^2;

+

+//Calculations:

+sigma_x = -(P/A)*10;

+sigma_y = 0;

+tou_xy = (T*r*10)/J;

+

+k = (sigma_x + sigma_y)/2;

+R = sqrt(k^2 + (tou_xy^2));

+

+sigma1 = k+R;

+sigma2 = k-R;

+l = sigma1 - sigma2;

+

+//Maximum Shear Stress Theory:

+test1 = abs(l);

+

+if(test1 >= sigma_yield)

+    

+    printf("\n\nFailure occurs by Maximum Shear Stress Theory.');

+end

+

+

+//Maximum Distortion-Energy Theory:

+test2 = (sigma1^2 - (sigma1*sigma2) + sigma2^2);

+

+

+if(test2<sigma_y_sqr)

+ 

+    printf("\n\nFailure will not occur by Maximum Distortion-Energy Theory.');

+end

diff --git a/1898/CH10/EX10.2/Ex10_2.sce b/1898/CH10/EX10.2/Ex10_2.sce
new file mode 100755
index 000000000..2e207ebca
--- /dev/null
+++ b/1898/CH10/EX10.2/Ex10_2.sce
@@ -0,0 +1,36 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.2 : ")

+

+//Given:

+ep_x = -350;//(*10^-6) Normal Strain

+ep_y = 200; //*(10^-6) Normal Strain

+gamma_xy = 80; //*(10^-6) Shear Strain

+

+//Orientation of the element:

+tan_thetap = (gamma_xy)/(ep_x - ep_y);

+thetap1 = (0.5)*(atan(tan_thetap));

+

+//Principal Strains:

+

+k = (ep_x + ep_y)/2;

+l = (ep_x - ep_y)/2;

+tou = gamma_xy/2;

+R = sqrt( (l)^2 + tou^2);

+ep1 = R + k;

+ep2 = k -R ;

+ep = [ep1 ep2];

+

+ep_x1 = k + l*cos(2*thetap1)+ tou*sin(2*thetap1);

+thetap1 = thetap1*(180/%pi);

+thetap2 = (90 + thetap1);

+thetap =[thetap1 thetap2];

+

+

+//Display:

+

+printf('\n\nThe orientation of the element in the positive counterclockwise direction    = %1.2f degrees, %1.2f  degrees ',thetap);

+printf('\nThe principal strains are                                                    = %1.0f *10^-6 , %1.0f *10^-6 ',ep);

+printf('\nThe principal strain in the new x direction is                               = %1.0f *10^-6 ',ep_x1);

+

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

diff --git a/1898/CH10/EX10.3/Ex10_3.sce b/1898/CH10/EX10.3/Ex10_3.sce
new file mode 100755
index 000000000..8ec926a0a
--- /dev/null
+++ b/1898/CH10/EX10.3/Ex10_3.sce
@@ -0,0 +1,34 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.3 : ")

+

+//Given:

+ep_x = -350;//(*10^-6) Normal Strain

+ep_y = 200; //*(10^-6) Normal Strain

+gamma_xy = 80; //*(10^-6) Shear Strain

+

+//Orientation of the element:

+tan_thetap = -(ep_x - ep_y)/(gamma_xy);

+thetap1 = (0.5)*(atan(tan_thetap));

+

+//Maximum in-plane shear strain:

+

+l = (ep_x - ep_y)/2;

+tou = gamma_xy/2;

+R = sqrt( l^2 + tou^2);

+max_inplane_strain = 2*R;

+

+gamma_xy_1 = (-l*sin(2*thetap1)+ tou*cos(2*thetap1))*2;

+strain_avg = (ep_x + ep_y)/2;

+

+thetap1 = thetap1*(180/%pi);

+thetap2 = (90 + thetap1);

+thetap =[thetap1 thetap2];

+

+//Display:

+

+printf('\n\nThe orientation of the element            = %1.1f degrees, %1.1f  degrees ',thetap);

+printf('\nThe maximum in-plane shear strain         = %1.0f *10^-6 ',max_inplane_strain);

+printf('\nThe average strain                        = %1.0f *10^-6 ',strain_avg);

+

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

diff --git a/1898/CH10/EX10.4/Ex10_4.sce b/1898/CH10/EX10.4/Ex10_4.sce
new file mode 100755
index 000000000..2b7c19def
--- /dev/null
+++ b/1898/CH10/EX10.4/Ex10_4.sce
@@ -0,0 +1,33 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.4 : ")

+

+//Given:

+ep_x = 250;//(*10^-6) Normal Strain

+ep_y = -150; //*(10^-6) Normal Strain

+gamma_xy = 120; //*(10^-6) Shear Strain

+

+//Construction of the circle:

+strain_avg = (ep_x + ep_y)/2;

+tou = gamma_xy/2;

+R = sqrt((ep_x - strain_avg)^2 + (tou^2));

+

+//Principal Strains:

+ep1 = (strain_avg + R);

+ep2 = (strain_avg - R);

+strain = [ep1 ep2];

+

+tan_thetap = (tou)/(ep_x - strain_avg);

+thetap1 = (atan(tan_thetap))/2;

+thetap1 = thetap1*(180/%pi);

+

+//Display:

+

+printf('\n\nThe principal strains are              = %1.0f *10^-6, %1.0f*10^-6 ',strain);

+printf('\nThe orientation of the element         = %1.2f degrees',thetap1);

+printf('\nThe average strain                     = %1.0f *10^-6 ',strain_avg);

+

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

+

+

+

diff --git a/1898/CH10/EX10.5/Ex10_5.sce b/1898/CH10/EX10.5/Ex10_5.sce
new file mode 100755
index 000000000..03e72e13a
--- /dev/null
+++ b/1898/CH10/EX10.5/Ex10_5.sce
@@ -0,0 +1,31 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.5 : ")

+

+//Given:

+ep_x = 250;//(*10^-6) Normal Strain

+ep_y = -150; //*(10^-6) Normal Strain

+gamma_xy = 120; //*(10^-6) Shear Strain

+

+//Orientation of the element:

+thetas = 90 - 2*8.35;

+thetas1 = thetas/2;

+

+//Maximum in-plane shear strain:

+

+l = (ep_x - ep_y)/2;

+tou = gamma_xy/2;

+R = sqrt( l^2 + tou^2);

+max_inplane_strain = 2*R;

+

+

+strain_avg = (ep_x + ep_y)/2;

+

+

+//Display:

+

+printf('\n\nThe orientation of the element            = %1.1f degrees ',thetas1);

+printf('\nThe maximum in-plane shear strain         = %1.0f *10^-6 ',max_inplane_strain);

+printf('\nThe average strain                        = %1.0f *10^-6 ',strain_avg);

+

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

diff --git a/1898/CH10/EX10.6/Ex10_6.sce b/1898/CH10/EX10.6/Ex10_6.sce
new file mode 100755
index 000000000..671c3aa7b
--- /dev/null
+++ b/1898/CH10/EX10.6/Ex10_6.sce
@@ -0,0 +1,41 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.6 : ")

+

+//Given:

+ep_x = -300;//(*10^-6) Normal Strain

+ep_y = -100; //*(10^-6) Normal Strain

+gamma_xy = 100; //*(10^-6) Shear Strain

+theta = 20; //degrees

+

+

+//Construction of the circle:

+strain_avg = (ep_x+ ep_y)/2;

+tou = gamma_xy/2;

+R = sqrt((-ep_x + strain_avg)^2 + tou^2);

+

+//Strains on Inclined Element:

+theta1 = 2*theta;

+

+phi = atan((tou)/(-ep_x +strain_avg));

+phi = phi*(180/%pi);

+psi = theta1 - phi;

+psi = psi*(%pi/180);

+

+ep_x1 = -(-strain_avg+ R*cos(psi));

+gamma_xy1 = -(R*sin(psi))*2;

+

+ep_y1 = -(-strain_avg - R*cos(psi));

+

+//Display:

+

+printf('\n\nThe normal strain in the new x direction       = %1.0f *10^-6 ',ep_x1);

+printf('\nThe normal strain in the new y direction       = %1.1f *10^-6 ',ep_y1);

+printf('\nThe shear strain in the new xy direction       = %1.0f *10^-6 ',gamma_xy1);

+printf('\nThe average strain                             = %1.0f *10^-6 ',strain_avg);

+

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

+

+

+

+

diff --git a/1898/CH10/EX10.7/Ex10_7.sce b/1898/CH10/EX10.7/Ex10_7.sce
new file mode 100755
index 000000000..111e19edc
--- /dev/null
+++ b/1898/CH10/EX10.7/Ex10_7.sce
@@ -0,0 +1,32 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.7 : ")

+

+//Given:

+ep_x = -400;//(*10^-6) Normal Strain

+ep_y = 200; //*(10^-6) Normal Strain

+gamma_xy = 150; //*(10^-6) Shear Strain

+

+//Maximum in-plane Shear Strain:

+strain_avg = (ep_x+ ep_y)/2;

+tou = gamma_xy/2;

+

+R = sqrt((-ep_x + strain_avg)^2 + tou^2); 

+strain_max = strain_avg + R;

+strain_min = strain_avg - R;

+

+max_shear_strain = strain_max - strain_min;

+

+//Absolute Maximum Shear Strain:

+abs_max_shear = max_shear_strain;

+

+//Display:

+

+printf('\n\nThe maximum in-plane principal strain     = %1.0f *10^-6 ',strain_max);

+printf('\nThe minimum in-plane principal strain     = %1.0f *10^-6 ',strain_min);

+printf('\nThe maximum in-plane shear strain         = %1.0f *10^-6 ',max_shear_strain);

+printf('\nThe absolute maximum shear strain         = %1.0f *10^-6 ',abs_max_shear);

+printf('\nThe average strain                        = %1.0f *10^-6 ',strain_avg);

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

+

+

diff --git a/1898/CH10/EX10.8/Ex10_8.sce b/1898/CH10/EX10.8/Ex10_8.sce
new file mode 100755
index 000000000..b8406b7e1
--- /dev/null
+++ b/1898/CH10/EX10.8/Ex10_8.sce
@@ -0,0 +1,63 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.8 : ")

+

+//Given:

+ep_a = 60;//(*10^-6) Normal Strain

+ep_b = 135; //*(10^-6) Normal Strain

+ep_c = 264; //*(10^-6) Normal Strain

+

+theta_a = 0;

+theta_b = 60*(%pi/180);

+theta_c = 120*(%pi/180);

+

+//Using matrices to solve the equations: 

+a1 = (cos(theta_a))^2;

+b1 = (sin(theta_a))^2;

+c1 = cos(theta_a)*sin(theta_a);

+

+a2 = (cos(theta_b))^2;

+b2 = (sin(theta_b))^2;

+c2 = cos(theta_b)*sin(theta_b);

+

+a3 = (cos(theta_c))^2;

+b3 = (sin(theta_c))^2;

+c3 = cos(theta_c)*sin(theta_c);

+

+A = [a1 b1 c1 ; a2 b2 c2; a3 b3 c3 ]

+b = [ep_a ; ep_b ; ep_c];

+strain = A\b;

+

+ep_x = strain(1);

+ep_y = strain(2);

+gamma_xy = strain(3);

+

+strain_avg = (ep_x + ep_y )/2;

+tou = gamma_xy/2;

+

+R = sqrt((-ep_x + strain_avg)^2 + tou^2); 

+

+ep1 = strain_avg + R;

+ep2 = strain_avg - R;

+ep = [ep1 ep2];

+

+tan_thetap =atan(-tou/(-ep_x + strain_avg));

+thetap = tan_thetap/2;

+thetap2 = thetap*(180/%pi);

+

+//Display:

+

+

+printf('\n\nThe maximum in-plane principal strains are     = %1.0f *10^-6 , %1.1f *10^-6',ep);

+printf('\nThe angle of orientation                       = %1.1f degrees',thetap2);

+

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

+

+

+ 

+ 

+

+

+

+

+

diff --git a/1898/CH10/EX10.9/Ex10_9.sce b/1898/CH10/EX10.9/Ex10_9.sce
new file mode 100755
index 000000000..cb3b0f6de
--- /dev/null
+++ b/1898/CH10/EX10.9/Ex10_9.sce
@@ -0,0 +1,25 @@
+clear all; clc;

+

+disp("Scilab Code Ex 10.9 : ")

+

+//Given:

+E_st = 200*10^9; //GPa

+nu_st = 0.3; //Poisson's ratio

+ep1 = 272 *10^-6;

+ep2 = 33.8 *10^-6;

+

+//Solving for sigma using matrices:

+A = [1 -nu_st; -nu_st 1];

+b = [(ep1*E_st) ; (ep2*E_st)];

+sigma = A\b;

+

+sigma1= sigma(1)/(10^6);

+sigma2= sigma(2)/(10^6);

+

+//Display:

+

+printf('\n\nThe principal stresses at point A are    = %1.0f MPa , %1.1f MPa',sigma1, sigma2);

+

+

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

+