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

---
 3289/CH8/EX8.1/Ex8_1.sce | 29 +++++++++++++++++++++++++++++
 3289/CH8/EX8.2/Ex8_2.sce | 43 +++++++++++++++++++++++++++++++++++++++++++
 3289/CH8/EX8.3/Ex8_3.sce | 23 +++++++++++++++++++++++
 3289/CH8/EX8.4/Ex8_4.sce | 34 ++++++++++++++++++++++++++++++++++
 3289/CH8/EX8.5/Ex8_5.sce | 48 ++++++++++++++++++++++++++++++++++++++++++++++++
 3289/CH8/EX8.6/Ex8_6.sce | 23 +++++++++++++++++++++++
 6 files changed, 200 insertions(+)
 create mode 100755 3289/CH8/EX8.1/Ex8_1.sce
 create mode 100755 3289/CH8/EX8.2/Ex8_2.sce
 create mode 100755 3289/CH8/EX8.3/Ex8_3.sce
 create mode 100755 3289/CH8/EX8.4/Ex8_4.sce
 create mode 100755 3289/CH8/EX8.5/Ex8_5.sce
 create mode 100755 3289/CH8/EX8.6/Ex8_6.sce

(limited to '3289/CH8')

diff --git a/3289/CH8/EX8.1/Ex8_1.sce b/3289/CH8/EX8.1/Ex8_1.sce
new file mode 100755
index 000000000..c5999edcf
--- /dev/null
+++ b/3289/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,29 @@
+clc
+
+di=0.3 //m
+de=0.4 //m
+v=0.3
+sigmathetamax=250*10^6 //Pa
+p0=0
+pi=0
+
+//solution a:
+a=0.15
+b=0.2
+r=a
+//sigmathetamax=pi*((b^2+a^2)/(b^2-a^2))
+pi=sigmathetamax*((b^2-a^2)/(b^2+a^2))
+disp(pi,"in Pa is= ")
+
+//solution b:
+r=a
+//sigmathetamax=-2*p0*(b^2/(b^2-a^2))
+p0=-(-sigmathetamax)*((b^2-a^2)/(2*b^2))
+disp(p0,"in Pa is= ")
+
+//solution c:
+u=((a^3*pi)/(b^2-a^2))*(0.7+1.3*(b^2/a^2))
+disp(u,"in per E meter is= ")
+sigmaz=(pi*a^2-p0*b^2)/(b^2-a^2)
+disp(sigmaz,"for longitudinal stress is")
+
diff --git a/3289/CH8/EX8.2/Ex8_2.sce b/3289/CH8/EX8.2/Ex8_2.sce
new file mode 100755
index 000000000..62973ed4b
--- /dev/null
+++ b/3289/CH8/EX8.2/Ex8_2.sce
@@ -0,0 +1,43 @@
+clc
+
+sigmayp=340 //MPa
+tauyp=sigmayp/2 //MPa
+disp(tauyp,"in MPa is=")
+a=0.1 //m
+b=0.15 //m
+v=0.3 
+//pi=4*p0
+//sigmatheta=(pi*(a^2+b^2)-2*p0*b^2)/(b^2-a^2)
+//sigmatheta=1.7*pi
+
+//sloution a: maxi principal stress theory
+sigmatheta=1.7
+pi=sigmayp/sigmatheta
+disp(pi,"in MPa is= ")
+
+//sloution b: maxi shearing stress theory
+//(sigmatheta-sigmar)/2=1.35*pi
+pi=tauyp/1.35
+disp(pi,"in MPa is= ")
+
+//solution c: energy of distortion theory
+sigmar=-1
+sigmayp1=sqrt(sigmatheta^2+sigmar^2-sigmatheta*sigmar)//*pi
+disp(sigmayp1)
+pi=sigmayp/sigmayp1
+disp(pi,"in MPa is=")
+
+//solution d: maxi principal strain theory
+//(sigmatheta-v*sigmar)/E=sigmayp/E
+pi=sigmayp/(sigmatheta-v*sigmar)
+disp(pi,"in MPa is= ")
+
+//solution e: octahedral shearing stress theory:
+pi=(sqrt(2)*sigmayp)/sqrt((sigmatheta-sigmar)^2+sigmar^2+(-sigmatheta)^2)
+disp(pi,"in MPa is= ")
+
+
+
+
+
+
diff --git a/3289/CH8/EX8.3/Ex8_3.sce b/3289/CH8/EX8.3/Ex8_3.sce
new file mode 100755
index 000000000..cc50373ac
--- /dev/null
+++ b/3289/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,23 @@
+clc
+
+a=0.15 //m
+b=0.2 //m
+c=0.25 //m
+E=200*10^9 //Pa
+delta=0.0001 //m
+140 //MPa
+
+p=((E*delta)/8)*(((b^2-a^2)*(c^2-b^2))/(2*(b^2)*(c^2-a^2)))
+disp(p,"the contact pressure in Pa is= ") // textbook ans is wrong
+
+p=12.3*10^6
+sigmatheta=p*((b^2+c^2)/(c^2-b^2)) // where r=0.2
+disp(sigmatheta,"tangential stresses in the outer cylinder in Pa is= ")
+sigmatheta1=(2*p*b^2)/(c^2-b^2) // where r=0.25
+disp(sigmatheta1,"tangential stresses in the outer cylinder in Pa is= ")
+sigmatheta3=-(2*p*b^2)/(b^2-a^2) // where r=0.15
+disp(sigmatheta3,"tangential stresses in the inner cylinder in Pa is= ")
+sigmatheta4=-p*((b^2+a^2)/(b^2-a^2)) // where r=0.2
+disp(sigmatheta4,"tangential stresses in the inner cylinder in Pa is= ")
+
+
diff --git a/3289/CH8/EX8.4/Ex8_4.sce b/3289/CH8/EX8.4/Ex8_4.sce
new file mode 100755
index 000000000..dbc16b730
--- /dev/null
+++ b/3289/CH8/EX8.4/Ex8_4.sce
@@ -0,0 +1,34 @@
+clc
+
+dn=0.1 //m
+do=0.5 //m
+t=0.08 //m
+w=6900*(2*%pi/60) //rpm
+row=7.8*10^3//Ns^2/m^4
+E=200*10^9 //Pa
+v=0.3
+b=0.05
+c=0.25
+
+
+//solution a:
+//ud=((0.05*3.3*0.7)*(0.0025+0.0625-(1.3/3.3)*0.0025+(1.3/0.7)*0.0625)*row*w^2)/(8*E)
+ud=((0.05*3.3*0.7)*(b^2+c^2-(1.3/3.3)*b^2+(1.3/0.7)*c^2))/(8)
+disp(ud,"radial displacement of the disk in meter is= ")
+
+//us=((0.05*0.7)*(3.3*0.0025-1.3*0.0025)*row*w^2)/(8*E)
+us=((0.05*0.7)*(3.3*b^2-1.3*b^2))/(8)
+disp(us,"radial displacement of the shaft in meter is= ")
+delta=(ud-us)*row*w^2/E
+disp(delta)
+
+//solution b:
+//p=E*delta*(c^2-b^2)/(2*b*c^2)
+p=E*delta*(c^2-b^2)/(2*b*c^2)
+disp(p,"in Pa is= ")
+sigmathetamax=p*(c^2+b^2)/(c^2-b^2)
+disp(sigmathetamax,"in Pa is= ")
+
+//solution c:
+sigmathetamax=3.3*(b^2+c^2-(1.9/3.3)*b^2+c^2)*row*w^2/8
+disp(sigmathetamax,"in Pa is= ")
diff --git a/3289/CH8/EX8.5/Ex8_5.sce b/3289/CH8/EX8.5/Ex8_5.sce
new file mode 100755
index 000000000..73f126afb
--- /dev/null
+++ b/3289/CH8/EX8.5/Ex8_5.sce
@@ -0,0 +1,48 @@
+clc
+
+ti=0.075 //m
+to=0.015//m
+a=0.05//m
+b=0.25//m
+delta=0.05 //mm
+w=6900*(2*%pi/60) //rpm
+s=1
+row=7.8*10^3//Ns^2/m^4
+E=200 //GPa
+
+//solution a:
+t1=ti*a^s
+disp(t1,"t1 is=")
+t1=to*b^2
+disp(t1,"t1 is=")
+//(ti/to)=(t1*a^-s)/(t1*b^-s)=(b/a)^s
+c=(b/a)^s
+
+(ti/to)==c
+disp(ti/to,"ti/t0 is=")
+m1=-0.5+sqrt((0.5)^2+(1+0.3*1))
+disp(m1,"m1 is=")
+m2=-0.5-sqrt((0.5)^2+(1+0.3*1))
+disp(m2,"m2 is=")
+
+//sigmar=0=(c1/t1)*(0.05)^m1+(c2/t1)*(0.05)^(m2)-0.00176*row*w^2     // r=0.05
+//sigmar=0=(c1/t1)*(0.25)^m1+(c2/t1)*(0.25)^(m2)-0.0439*row*w^2     // r=0.25
+
+c1=t1*0.12529*row*w^2
+disp(c1,"c1 is=")
+c2=t1*-6.272*10^-5*row*w^2
+disp(c2,"c2 is=")
+
+r=0.05
+sigmar=(0.12529*r^0.745-6.272*10^-5*r^(-1.745)-0.70*r^2)//*row*w^2
+disp(sigmar,"sigmar is= ")
+
+sigmatheta=(0.09334*r^0.745+1.095*10^-4*r^(-1.745)-0.40*r^2)//*row*w^2
+disp(sigmatheta,"sigmatheta is= ")
+
+//solution b:
+r=0.05
+//ur=(r*sigmatheta)/E
+ur=(r*sigmatheta)
+disp(ur,"ur is= ")
+
diff --git a/3289/CH8/EX8.6/Ex8_6.sce b/3289/CH8/EX8.6/Ex8_6.sce
new file mode 100755
index 000000000..0fcdb75cd
--- /dev/null
+++ b/3289/CH8/EX8.6/Ex8_6.sce
@@ -0,0 +1,23 @@
+clc
+
+b=0.25 //m
+w=6900*(2*%pi/60) //rpm
+t1=0.075 //m
+t2=0.015 //m
+row=7.8*10^3//Ns^2/m^4
+c1=t1
+
+x=t2/t1
+disp(x)
+
+//(t2/t1)==(c1*exp(-(row*w^2/2*sigma)*b^2))/c1
+//exp(-(row*w^2/2*sigma)*b^2)=x
+//log(x)=-(row*w^2*b^2/2*sigma)
+y=2*log(x)
+disp(y)
+sigma=-(row*w^2*b^2)/y
+disp(sigma,"in Pa is= ")
+
+//t=c1*exp(-row*(w^2/2*sigma)*r^2)
+z=row*(w^2/(2*sigma))
+disp(z)
-- 
cgit