initial commit / add all books

6 files changed, 200 insertions, 0 deletions

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)