initial commit / add all books

12 files changed, 310 insertions, 0 deletions

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)