From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1898/CH10/EX10.1/Ex10_1.sce   | 29 ++++++++++++++++++++
 1898/CH10/EX10.10/Ex10_10.sce | 38 ++++++++++++++++++++++++++
 1898/CH10/EX10.11/Ex10_11.sce | 34 +++++++++++++++++++++++
 1898/CH10/EX10.12/Ex10_12.sce | 34 +++++++++++++++++++++++
 1898/CH10/EX10.13/Ex10_13.sce | 22 +++++++++++++++
 1898/CH10/EX10.14/Ex10_14.sce | 42 +++++++++++++++++++++++++++++
 1898/CH10/EX10.2/Ex10_2.sce   | 36 +++++++++++++++++++++++++
 1898/CH10/EX10.3/Ex10_3.sce   | 34 +++++++++++++++++++++++
 1898/CH10/EX10.4/Ex10_4.sce   | 33 +++++++++++++++++++++++
 1898/CH10/EX10.5/Ex10_5.sce   | 31 +++++++++++++++++++++
 1898/CH10/EX10.6/Ex10_6.sce   | 41 ++++++++++++++++++++++++++++
 1898/CH10/EX10.7/Ex10_7.sce   | 32 ++++++++++++++++++++++
 1898/CH10/EX10.8/Ex10_8.sce   | 63 +++++++++++++++++++++++++++++++++++++++++++
 1898/CH10/EX10.9/Ex10_9.sce   | 25 +++++++++++++++++
 14 files changed, 494 insertions(+)
 create mode 100755 1898/CH10/EX10.1/Ex10_1.sce
 create mode 100755 1898/CH10/EX10.10/Ex10_10.sce
 create mode 100755 1898/CH10/EX10.11/Ex10_11.sce
 create mode 100755 1898/CH10/EX10.12/Ex10_12.sce
 create mode 100755 1898/CH10/EX10.13/Ex10_13.sce
 create mode 100755 1898/CH10/EX10.14/Ex10_14.sce
 create mode 100755 1898/CH10/EX10.2/Ex10_2.sce
 create mode 100755 1898/CH10/EX10.3/Ex10_3.sce
 create mode 100755 1898/CH10/EX10.4/Ex10_4.sce
 create mode 100755 1898/CH10/EX10.5/Ex10_5.sce
 create mode 100755 1898/CH10/EX10.6/Ex10_6.sce
 create mode 100755 1898/CH10/EX10.7/Ex10_7.sce
 create mode 100755 1898/CH10/EX10.8/Ex10_8.sce
 create mode 100755 1898/CH10/EX10.9/Ex10_9.sce

(limited to '1898/CH10')

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--------------------------------------------------------------------------------------
+
-- 
cgit