From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 3705/CH8/EX8.1/Ex8_1.sce   | 29 +++++++++++++++++++++++++++++
 3705/CH8/EX8.10/Ex8_10.sce | 22 ++++++++++++++++++++++
 3705/CH8/EX8.11/Ex8_11.sce | 32 ++++++++++++++++++++++++++++++++
 3705/CH8/EX8.12/Ex8_12.sce | 28 ++++++++++++++++++++++++++++
 3705/CH8/EX8.13/Ex8_13.sce | 24 ++++++++++++++++++++++++
 3705/CH8/EX8.14/Ex8_14.sce | 29 +++++++++++++++++++++++++++++
 3705/CH8/EX8.3/Ex8_3.sce   | 28 ++++++++++++++++++++++++++++
 3705/CH8/EX8.4/Ex8_4.sce   | 17 +++++++++++++++++
 3705/CH8/EX8.5/Ex8_5.sce   | 27 +++++++++++++++++++++++++++
 3705/CH8/EX8.6/Ex8_6.sce   | 23 +++++++++++++++++++++++
 3705/CH8/EX8.7/Ex8_7.sce   | 32 ++++++++++++++++++++++++++++++++
 3705/CH8/EX8.9/Ex8_9.sce   | 19 +++++++++++++++++++
 12 files changed, 310 insertions(+)
 create mode 100644 3705/CH8/EX8.1/Ex8_1.sce
 create mode 100644 3705/CH8/EX8.10/Ex8_10.sce
 create mode 100644 3705/CH8/EX8.11/Ex8_11.sce
 create mode 100644 3705/CH8/EX8.12/Ex8_12.sce
 create mode 100644 3705/CH8/EX8.13/Ex8_13.sce
 create mode 100644 3705/CH8/EX8.14/Ex8_14.sce
 create mode 100644 3705/CH8/EX8.3/Ex8_3.sce
 create mode 100644 3705/CH8/EX8.4/Ex8_4.sce
 create mode 100644 3705/CH8/EX8.5/Ex8_5.sce
 create mode 100644 3705/CH8/EX8.6/Ex8_6.sce
 create mode 100644 3705/CH8/EX8.7/Ex8_7.sce
 create mode 100644 3705/CH8/EX8.9/Ex8_9.sce

(limited to '3705/CH8')

diff --git a/3705/CH8/EX8.1/Ex8_1.sce b/3705/CH8/EX8.1/Ex8_1.sce
new file mode 100644
index 000000000..d62c4b00a
--- /dev/null
+++ b/3705/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,29 @@
+
+clear//
+
+//Variable Declaration
+p=125 //Pressure in psi
+r=24 //Radius of the vessel in inches
+t=0.25 //Thickness of the vessel in inches
+E=29*10**6 //Modulus of Elasticity in psi
+v=0.28 //poisson ratio
+
+//Calcualtions
+//Part 1
+sigma_c=p*r*t**-1 //Circumferential Stress in psi
+sigma_l=sigma_c/2 //Longitudinat Stress in psi
+e_c=E**-1*(sigma_c-(v*sigma_l)) //Circumferential Strain using biaxial Hooke's Law 
+delta_r=e_c*r //Change in the radius in inches
+
+//Part 2
+sigma=(p*r)*(2*t)**-1 //Stress in psi
+e=E**-1*(sigma-(v*sigma)) //Strain using biaxial Hooke's Law
+delta_R=e*r //Change inradius of end-cap in inches
+
+//Result
+printf("\n Part 1 Answers")
+printf("\n Stresses are sigma_c= %0.0f psi and sigma_l= %0.0f psi",sigma_c,sigma_l)
+printf("\n Change of radius of cylinder= %0.5f in",delta_r)
+printf("\n Part 2 Answers")
+printf("\n Stresses are sigma= %0.0f psi",sigma)
+printf("\n Change in radius of end cap= %0.5f in",delta_R)
diff --git a/3705/CH8/EX8.10/Ex8_10.sce b/3705/CH8/EX8.10/Ex8_10.sce
new file mode 100644
index 000000000..f0042a5f9
--- /dev/null
+++ b/3705/CH8/EX8.10/Ex8_10.sce
@@ -0,0 +1,22 @@
+
+clear//
+
+//Variable Declaration
+t=0.01 //Thickness of the shaft in m
+p=2 //Internal Pressure in MPa
+r=0.45 //Mean radius of the vessel in m
+tw=50 //Working shear stress in MPa
+
+//Calculation
+sigma_x=(p*r)/(2*t) //Stress in MPa
+sigma_y=(p*r)/t //Stress in MPa
+
+R=100-67.5 //From the diagram in MPa
+tau_xy=sqrt((R**2-(sigma_y-67.5)**2)) //Stress in MPa
+
+J=2*%pi*r**3*t //Polar Moment of inertia in mm^4
+
+T=1000*(tau_xy*J)/r //Maximum allowable Torque in kN.m
+
+//Result
+printf("\n The largest allowable Torque is %0.0f kN.m",T)
diff --git a/3705/CH8/EX8.11/Ex8_11.sce b/3705/CH8/EX8.11/Ex8_11.sce
new file mode 100644
index 000000000..d763804e0
--- /dev/null
+++ b/3705/CH8/EX8.11/Ex8_11.sce
@@ -0,0 +1,32 @@
+
+clear//
+
+//Variable Declaration
+L=15 //Length of the shaft in inches
+r=3.0/8.001 //Radius of the shaft in inches
+T=540 //Torque applied in lb.in
+
+//Calculations
+V=30 //Transverse Shear Force in lb
+M=15*V //Bending Moment in lb.in
+I=(%pi*r**4)/4.0 //Moment of Inertia in in^4
+J=2*I //Polar Moment Of Inertia in in^4
+
+//Part 1
+sigma=(M*r)/I //Bending Stress in psi
+tau_t=10**-3*(T*r)/J //Shear Stress in ksi
+
+sigma_max1=13.92 //From the Mohr Circle in ksi
+
+//Part 2
+Q=(2*r**3)/3.0 //First Moment in in^3
+b=2*r // in
+
+tau_V=10**-3*(V*Q)/(I*b) //Shear Stress in ksi
+tau=tau_t+tau_V //Total Shear in ksi
+
+sigma_max2=tau //Maximum stress in ksi
+
+//Result
+printf("\n The maximum normal stress in case 1 is %0.3f ksi",sigma_max1)
+printf("\n The Maximum normal stress in case 2 is %0.2f ksi",sigma_max2)
diff --git a/3705/CH8/EX8.12/Ex8_12.sce b/3705/CH8/EX8.12/Ex8_12.sce
new file mode 100644
index 000000000..f4a0adade
--- /dev/null
+++ b/3705/CH8/EX8.12/Ex8_12.sce
@@ -0,0 +1,28 @@
+
+clear//
+
+//Variable Declaration
+ex=800*10**-6 //Strain in x (no units)
+ey=200*10**-6 //Strain in y(no units)
+y_xy=-600*10**-6 //Strain in xy(no units)
+
+//Calculations
+e_bar=(ex+ey)*0.5 //Strain 
+R=sqrt(2*300**2)*10**-6 //Resultant 
+
+//Part 1
+e1=e_bar+R //Strain in Principal Axis(no units)
+e2=e_bar-R //Strain in Principal Axis(no units)
+
+//Part 2
+alpha=15*180**-1*%pi //From the Mohr Circle in degrees
+e_xbar=e_bar-(R*cos(alpha)) //Strain in x (no units)
+e_ybar=e_bar+(R*cos(alpha)) //Strain in y(no units)
+
+shear_strain=2*R*sin(alpha) //Shear follows 
+
+//Result
+printf("\n The principal Strains are")
+printf("\n e1= %0.6f e2= %0.6f ",e1,e2)
+printf("\n The follows components are")
+printf("\n ex= %0.6f ey= %0.6f y_xy= %0.6f ",e_xbar,e_ybar,shear_strain)
diff --git a/3705/CH8/EX8.13/Ex8_13.sce b/3705/CH8/EX8.13/Ex8_13.sce
new file mode 100644
index 000000000..29735881b
--- /dev/null
+++ b/3705/CH8/EX8.13/Ex8_13.sce
@@ -0,0 +1,24 @@
+
+clear//
+
+//Variable Declaration
+e_x=800*10**-6 //Strain in x
+e_y=200*10**-6 //Strain in y
+y_xy=-600*10**-6 //Strain in xy
+v=0.30 //Poissons Ratio
+E=200 //Youngs Modulus in GPa
+R_e=424.3*10**-6 //Strain
+e_bar=500*10**-6 //Strain
+
+//Calculations
+//Part 1
+R_sigma=10**-6*R_e*(E*10**9/(1+v)) //Stress in MPa
+sigma_bar=10**-6*e_bar*(E*10**9/(1-v)) //Stress in MPa
+
+//Part 2
+sigma1=sigma_bar+R_sigma //Principal Stress in MPa
+sigma2=sigma_bar-R_sigma //Principal Stress in MPa
+
+//Result
+printf("\n The principal Stresses are as follows")
+printf("\n Sigma1= %0.0f MPa and Sigma2= %0.1f MPa",sigma1,sigma2)
diff --git a/3705/CH8/EX8.14/Ex8_14.sce b/3705/CH8/EX8.14/Ex8_14.sce
new file mode 100644
index 000000000..5383e59be
--- /dev/null
+++ b/3705/CH8/EX8.14/Ex8_14.sce
@@ -0,0 +1,29 @@
+
+clear//
+
+//Variable Declaration
+e_a=100*10**-6 //Strain 
+e_b=300*10**-6 //Strain
+e_c=-200*10**-6 //Strain
+E=180 //Youngs Modulus in GPa
+v=0.28 //Poissons Ratio 
+
+//Calculations
+y_xy=(e_b-(e_a+e_c)*0.5) //Strain in xy
+e_bar=(e_a+e_c)*0.5 //Strain 
+R_e=sqrt(y_xy**2+(150*10**-6)**2) //Resultant Strain
+
+//Corresponding Parameters from Mohrs Diagram
+sigma_bar=(E/(1-v))*e_bar*10**3 //Stress in MPa
+R_sigma=(E/(1+v))*R_e*10**3 //Resultant Stress in MPa
+
+//Principal Stresses
+sigma1=sigma_bar+R_sigma //MPa
+sigma2=sigma_bar-R_sigma //MPa
+
+theta=atan(y_xy/(150*10**-6))*180*%pi**-1*0.5 //Degrees
+
+//Result
+printf("\n The Principal Stresses are as follows")
+printf("\n Sigma1= %0.1f MPa and Sigma2= %0.2f MPa",sigma1,sigma2)
+printf("\n The plane orientation is %0.2f degrees",theta)
diff --git a/3705/CH8/EX8.3/Ex8_3.sce b/3705/CH8/EX8.3/Ex8_3.sce
new file mode 100644
index 000000000..5e4560e43
--- /dev/null
+++ b/3705/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,28 @@
+
+clear//
+
+//Variable Declaration
+b=6 //Width in inches
+h=10 //Depth in inches
+P1=6000 //Force in lb
+P2=3000 //Force in lb
+L=4 //Length in ft
+P=-13400 //Load in lb
+M=6000 //Moment in lb.ft
+y=5 //Depth in inches
+P2=-9800 //Load in lb
+M2=-12000 //Moment in lb.ft
+
+//Calculations
+A=b*h //Area in in^2
+I=b*h**3*12**-1 //Moment of inertia in in^4
+T=(P1*L+P2*L*3)*(6)**-1 //Tension in the cable in lb
+
+//Computation of largest stress
+sigma_B=(P*A**-1)-(M*y*12*I**-1) //Maximum Compressive Stress caused by +ve BM in psi
+sigma_C=(P2*A**-1)-(M2*-y*12*I**-1) //Maximum Compressive Stress caused by -ve BM in psi
+
+sigma_max=max(-sigma_B,-sigma_C) //Maximum Compressive Stress in psi
+
+//Result
+printf("\n The maximum Stress is %0.0f psi",sigma_max)
diff --git a/3705/CH8/EX8.4/Ex8_4.sce b/3705/CH8/EX8.4/Ex8_4.sce
new file mode 100644
index 000000000..ec1510039
--- /dev/null
+++ b/3705/CH8/EX8.4/Ex8_4.sce
@@ -0,0 +1,17 @@
+
+clear//
+
+//Variable Declaration
+theta=(60*%pi)/180 //Angle in radians (Twice as declared)
+sigma_x=30 // Stress in x in MPa
+sigma_y=60 //Stress in y in MPa
+tau_xy=40 //Stress in MPa
+
+//Calcualtions
+sigma_xdash=0.5*(sigma_x+sigma_y)+0.5*(sigma_x-sigma_y)*cos(theta)+tau_xy*sin(theta) //Stress at x' axis in MPa
+sigma_ydash=0.5*(sigma_x+sigma_y)-0.5*(sigma_x-sigma_y)*cos(theta)-tau_xy*sin(theta) //Stress at y' axis in MPa
+tau_x_y=-0.5*(sigma_x-sigma_y)*sin(theta)+tau_xy*cos(theta) //Stress at x'y' in shear in MPa
+//Result
+printf("\n The new stresses at new axes are as follows")
+printf("\n sigma_x= %0.1f MPa sigma_y= %0.1f MPa",sigma_xdash,sigma_ydash)
+printf("\n And tau_xy= %0.0f MPa",tau_x_y)
diff --git a/3705/CH8/EX8.5/Ex8_5.sce b/3705/CH8/EX8.5/Ex8_5.sce
new file mode 100644
index 000000000..39cad6526
--- /dev/null
+++ b/3705/CH8/EX8.5/Ex8_5.sce
@@ -0,0 +1,27 @@
+
+clear//
+
+//Variable Declaration
+sigma_x=8000 //Stress in x in psi
+sigma_y=4000 //Stress in y in psi
+tau_xy=3000 //Stress in xy in psi
+
+//Calculations
+R=sqrt(((sigma_x-sigma_y)*0.5)**2+tau_xy**2) //Resultant Stress in psi
+
+//Principal Stresses
+sigma1=(sigma_x+sigma_y)*0.5+R //Principal Stress in psi
+sigma2=(sigma_x+sigma_y)*0.5-R //Principal Stress in psi
+
+//Principal Direction
+theta1=atan(2*tau_xy*(sigma_x-sigma_y)**-1)*0.5*180*%pi**-1 //Principal direction in degrees
+theta2=theta1+90 //Second pricnipal direction in degrees
+
+//Normal Stress
+sigma_xdash=0.5*(sigma_x+sigma_y)+0.5*(sigma_x-sigma_y)*cos(2*theta1*%pi*180**-1)+tau_xy*sin(2*theta1*%pi*180**-1)
+
+//Result
+printf("\n The principal stresses are as follows")
+printf("\n sigma1= %0.0f psi and sigma2= %0.0f psi",sigma1,sigma2)
+printf("\n The corresponding directions are")
+printf("\n Theta1= %0.1f degrees and Theta2= %0.1f degrees",theta1,theta2)
diff --git a/3705/CH8/EX8.6/Ex8_6.sce b/3705/CH8/EX8.6/Ex8_6.sce
new file mode 100644
index 000000000..66c8e4b54
--- /dev/null
+++ b/3705/CH8/EX8.6/Ex8_6.sce
@@ -0,0 +1,23 @@
+
+clear//
+
+//Variable Declaration
+sigma_x=40 //Stress in x in MPa
+sigma_y=-100 //Stress in y in MPa
+tau_xy=-50 //Shear stress in MPa
+
+//Calculations
+tau_max=sqrt(((sigma_x-sigma_y)*0.5)**2+tau_xy**2) //Maximum in-plane shear in MPa
+
+//Orientation of Plane
+theta1=atan(-((sigma_x-sigma_y)*(2*tau_xy)**-1))*180*%pi**-1*0.5 //Angle in Degrees
+theta2=theta1+90 //Angle in degrees
+
+//Plane of max in-plane shear
+tau_x_y=-0.5*(sigma_x-sigma_y)*sin(2*theta1*%pi*180**-1)+tau_xy*cos(2*theta1*%pi*180**-1) 
+
+//Normal Stress
+sigma=(sigma_x+sigma_y)*0.5 //Stress in MPa
+
+//Result
+printf("\n The maximum in-plane Shear is %0.0f MPa",tau_x_y)
diff --git a/3705/CH8/EX8.7/Ex8_7.sce b/3705/CH8/EX8.7/Ex8_7.sce
new file mode 100644
index 000000000..fefb67c2f
--- /dev/null
+++ b/3705/CH8/EX8.7/Ex8_7.sce
@@ -0,0 +1,32 @@
+
+clear//
+
+//Vairable Declaration
+sigma_x=40 //Stress in x in MPa
+sigma_y=20 //Stress in y in MPa
+tau_xy=16 //Shear in xy in MPa
+
+//Calculations
+sigma=(sigma_x+sigma_y)*0.5 //Normal Stress in MPa
+R=sqrt(((sigma_x-sigma_y)*0.5)**2+tau_xy**2) //Resultant Stress in MPa
+
+//Part 1
+sigma1=sigma+R //Principal Stress in MPa
+sigma2=sigma-R //Principal Stress in MPa
+theta=atan(tau_xy*((sigma_x-sigma_y)*0.5)**-1)*180*%pi**-1*0.5 //Orientation in degrees
+
+//Part 2
+tau_max=18.87 //From figure in MPa
+
+//Part 3
+sigma_xdash=sigma+tau_max*cos((100-theta*2)*%pi*180**-1) //Stress in MPa
+sigma_ydash=sigma-tau_max*cos((100-theta*2)*%pi*180**-1) //Stress in MPa
+tau_x_y=tau_max*sin((100-2*theta)*%pi*180**-1) //Shear stress in MPa
+
+//Result
+printf("\n The principal Stresses are")
+printf("\n Sigma1= %0.1f MPa and Sigma2= %0.1f MPa",sigma1,sigma2)
+printf("\n The Principal Plane is at %0.0f degrees",theta)
+printf("\n The  Maximum Shear Stress is %0.3f MPa",tau_max)
+printf("\n Sigma_x= %0.0f MPa and Sigma_y= %0.2f MPa",sigma_xdash,sigma_ydash)
+printf("\n Tau_xy= %0.2f MPa",tau_x_y)
diff --git a/3705/CH8/EX8.9/Ex8_9.sce b/3705/CH8/EX8.9/Ex8_9.sce
new file mode 100644
index 000000000..503b2ef6f
--- /dev/null
+++ b/3705/CH8/EX8.9/Ex8_9.sce
@@ -0,0 +1,19 @@
+
+clear//
+
+//Variable Declaration
+sigma_w=120 //Working Stress in MPa
+tau_w=70 //Working Shear in MPa
+
+//Calcualtions
+//Section a-a
+M=3750 //Applied moment at section a-a in N.m
+T=1500 //Applied Torque at section a-a in N.m
+
+//After carrying out the variable based computation we compute d
+d1=((124.62)/(sigma_w*10**3*%pi))**0.3333 //Diameter of the shaft in m
+d2=((65.6)/(tau_w*10**3*%pi))**0.3333 //Diameter of the shaft in m
+d=max(d1,d2) //Diameter of the shaft to be selected in m
+
+//Result
+printf("\n The diameter of the shaft to be selected is %0.1f mm",d*1000)
-- 
cgit