initial commit / add all books

3 files changed, 109 insertions, 0 deletions

diff --git a/3392/CH3/EX3.1/Ex3_1.sce b/3392/CH3/EX3.1/Ex3_1.sce
new file mode 100755
index 000000000..17d2e6ffc
--- /dev/null
+++ b/3392/CH3/EX3.1/Ex3_1.sce
@@ -0,0 +1,13 @@
+clc

+// initialization of variables

+clear

+// part (a)

+E=72 // G Pa

+v=0.33 // Poisoon's ratio

+h=2 // mm

+R=600 // mm

+//calculations

+sig_cir=E*h/(2*(1-v^2)*R)

+// results

+printf('\n part (a) \n')

+printf(' The maximum circumferential stress is %.d M Pa',sig_cir*10^3)

diff --git a/3392/CH3/EX3.7/Ex3_7.sce b/3392/CH3/EX3.7/Ex3_7.sce
new file mode 100755
index 000000000..7b4339708
--- /dev/null
+++ b/3392/CH3/EX3.7/Ex3_7.sce
@@ -0,0 +1,36 @@
+clc

+// initialization of variables

+clear

+tR=0.02 // t/R ration

+E_A=69 //G Pa

+v_A=0.33  // Poisson's ratio

+alpha_A=21.6*10^-6 // /degree Celcius (Coefficient of expansion)

+E_S=207 // G Pa

+v_S=0.280

+alpha_S=10.8*10^-6 // /degree Celcius (Coefficient of expansion)

+// calculations

+// Sig_LA=a*p+b*delT+c*sig_thS

+// Sig_LS=v_S*Sig_thS+d*delT

+E_S=E_S*10^9

+E_A=E_A*10^9

+a=1/tR*E_A/E_S

+b=-2/3*alpha_S*E_S

+c=-E_A/E_S

+d=-alpha_S*E_S

+// SigthS=e*p+f*delT

+// SigthA=g*p+h*delT

+e=37.16

+f=0.8639*10^6

+g=1/tR-e

+h=-f

+// results

+p=689.4 // kPa

+delT=100 // degree Celcius

+p=p*10^3 // Pa

+SigthA=g*p+h*delT

+SigthS=e*p+f*delT

+Sig_LA=a*p+b*delT+c*SigthS

+Sig_LS=v_S*SigthS+d*delT

+printf('Thus, for p = %.1f k Pa and delT = %.d degree celcius \n',p/10^3,delT)

+printf(' SigthA = %.1f M Pa,   Sig_LA = %.d M Pa \n',SigthA/10^6,Sig_LA/10^6)

+printf(' SigthS = %.1f M Pa,   Sig_LS = %.d M Pa',SigthS/10^6,Sig_LS/10^6)

diff --git a/3392/CH3/EX3.8/Ex3_8.sce b/3392/CH3/EX3.8/Ex3_8.sce
new file mode 100755
index 000000000..e6dc72189
--- /dev/null
+++ b/3392/CH3/EX3.8/Ex3_8.sce
@@ -0,0 +1,60 @@
+clc

+// initialization of variables

+clear

+// Material constants

+Ex=14700 // M Pa

+Ey=1000 // M Pa

+Ez=735 // M Pa

+Gxy=941 // M Pa

+Gxz=1147 // M Pa

+Gyz=103 // M pa

+Vxy=0.292

+Vxz=0.449

+Vyz=0.39

+// Stresses at a point

+Sxx=7 // M pa

+Syy=2.1 // M Pa

+Szz=-2.8 //M Pa

+Sxy=1.4 // M Pa

+Sxz=0 //M Pa

+Syz=0 // M Pa

+// part (a)

+th=1/2*atan(2*Sxy/(Sxx-Syy))*180/%pi

+I1=Sxx+Syy+Szz

+I2=Sxx*Syy-Sxy^2+Szz*Sxx-Sxz^2+Szz*Syy-Syz^2

+M=[Sxx Sxy Sxz

+   Sxy Syy Syz

+   Sxz Syz Szz]

+I3=det(M)   

+p=[1 -I1 I2 -I3]

+S=roots(p)

+// results

+printf('Part (a) \n')

+printf('The maximum principal stress is S1 = %.2f M Pa', S(1))

+printf('\n and occurs in direction th = %.1f degrees',th)

+printf('\n and the intermediate principal stress S2 = %.2f M Pa occurs in the direction th = %.1f degrees \n',S(3),th+90)

+printf(' The minimum principal stress is S3 = Szz = %.1f M Pa', S(2))

+Ex=Ex*10^6

+Ey=Ey*10^6

+Ez=Ez*10^6

+Gxy=Gxy*10^6

+Gxz=Gxz*10^6

+Gyz=Gyz*10^6

+// part (b) is to find strains

+Exx=Sxx/Ex-Vxy*Syy/Ey-Vxz*Szz/Ez

+Eyy=-Vxy*Sxx/Ex+Syy/Ey-Vyz*Szz/Ez

+Ezz=-Vxz*Sxx/Ex-Vyz*Syy/Ey+Szz/Ez

+Exy=Sxy/Gxy

+Exz=Sxz/Gxz

+Eyz=Syz/Gyz

+printf('\n Part (b)')

+printf('\n The srains are')

+printf('\n Exx = %.2e ,   Eyy = %.2e ,  Ezz = %.4e',Exx,Eyy,Ezz)

+printf('\n Exy = %.4e ,  Exz = %.2d ,  Eyz = %.2d',Exy,Exz,Eyz)

+// Wrong Exx value in the textbook

+th=1/2*atan(Exy/(Exx-Eyy))

+th=th*180/%pi

+th2=th+90

+printf('\n part (c)')

+printf('\n theta = %.2f  or theta = %.2f degrees',th,th2)

+// Wrong theta too since Ex given in textbook is wrong