initial commit / add all books

9 files changed, 267 insertions, 0 deletions

diff --git a/3392/CH7/EX7.1/Ex7_1.sce b/3392/CH7/EX7.1/Ex7_1.sce
new file mode 100755
index 000000000..f1e9bdb5e
--- /dev/null
+++ b/3392/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,21 @@
+clc

+// initialization of variables

+clear

+E=200 //G Pa

+Y=250 //M Pa

+SF=1.9

+w=1 //kN/m

+L=3 //m

+S_max=Y

+// Calculations

+E=E*10^9

+Y=Y*10^6

+w=w*10^3

+Mx=-SF*w*L^2/2

+S_max=S_max*10^6

+k=2 // c_max=h/k

+//Formula to be used

+// S_max=abs(Mx)*c_max/Ix

+// Note that c_max=h/2 and Ix=h^4/24

+h=(abs(Mx)*24/(k*S_max))^(1/3)

+printf('h = %.4f m',h)

diff --git a/3392/CH7/EX7.2/Ex7_2.sce b/3392/CH7/EX7.2/Ex7_2.sce
new file mode 100755
index 000000000..838cb60a3
--- /dev/null
+++ b/3392/CH7/EX7.2/Ex7_2.sce
@@ -0,0 +1,32 @@
+clc

+// initialization of variables

+clear

+P1=1.5 //kN

+P2=4.5 //kN

+// part (a)

+A=1000 //mm^2

+A1=500 //mm^2

+A2=500 //mm^2

+//calculation

+A=A*10^-6

+A1=A1*10^-6

+A2=A2*10^-6

+y1=25*10^-3

+y2=55*10^-3

+c1=(A1*y1+A2*y2)/A

+c2=60*10^-3-c1 // c1+c2=60 mm

+y_1=c1-25*10^-3

+y_2=c2-5*10^-3

+b1=50*10^-3

+h1=10*10^-3

+h2=50*10^-3

+b2=10*10^-3

+Ix=1/12*b1*h1^3 + A1*y_1^2 + 1/12*b2*h2^3 + A2*y_2^2

+printf('part (a)')

+R1=2550 //N

+Vy=750 //N

+Mx=975 //N.m

+S_zzT=Mx*c1/Ix

+S_zzC=Mx*(-c2)/Ix

+printf('\n Maximum Tensile stress = %.1f MPa',S_zzT/10^6)

+printf('\n Maximum Compressive stress = %.1f MPa',S_zzC/10^6)

diff --git a/3392/CH7/EX7.3/Ex7_3.sce b/3392/CH7/EX7.3/Ex7_3.sce
new file mode 100755
index 000000000..2f57f1310
--- /dev/null
+++ b/3392/CH7/EX7.3/Ex7_3.sce
@@ -0,0 +1,27 @@
+clc

+// initialization of variables

+clear

+P=12 //kN

+Phi=%pi/3

+// calculations

+L=3 //m

+P=12 //kN

+A=10000 //mm^2

+Ix=39.69*10^6 //mm^4

+yo=82 //mm

+Iy=30.73*10^6 //mm^4

+Ixy=0

+P=P*10^3

+Ix=Ix*10^-12

+Iy=Iy*10^-12

+alpha=atan(-Ix/(Iy*tan(Phi)))

+M=-L*P

+Mx=M*sin(Phi)

+yA=-118*10^-3 //m

+xA=-70*10^-3 //m

+xB=-xA

+yB=82*10^-3 //m

+S_A=Mx*(yA-xA*tan(alpha))/(Ix-Ixy*tan(alpha))

+S_B=Mx*(yB-xB*tan(alpha))/(Ix-Ixy*tan(alpha))

+printf(' Sigma A = %.1f M Pa \n',S_A/10^6)

+printf(' Sigma B = %.1f M Pa',S_B/10^6)

diff --git a/3392/CH7/EX7.4/Ex7_4.sce b/3392/CH7/EX7.4/Ex7_4.sce
new file mode 100755
index 000000000..5f583f839
--- /dev/null
+++ b/3392/CH7/EX7.4/Ex7_4.sce
@@ -0,0 +1,52 @@
+clc

+// initialization of variables

+clear

+P=4 //kN

+L=1.2 //m

+A=1900 //mm^2

+Ix=2.783*10^6 //mm^4

+Iy=1.003*10^6 //mm^4

+Ixy=-0.973*10^6 //mm^4

+P=P*10^3

+Ix=Ix*10^-12

+Iy=Iy*10^-12

+Ixy=Ixy*10^-12

+A=1900 //mm^2

+xo=19.74 //mm

+yo=39.74 //mm

+Phi=2*%pi/3

+Nr=Ixy-Ix/tan(Phi)

+Dr=Iy-Ixy/tan(Phi)

+alpha=atan(Nr/Dr)

+M=L*P

+Mx=M*sin(Phi)

+yA=39.74*10^-3 //m

+xA=-60.26*10^-3 //m

+xB=19.74*10^-3

+yB=-80.26*10^-3 //m

+S_A=Mx*(yA-xA*tan(alpha))/(Ix-Ixy*tan(alpha))

+S_B=Mx*(yB-xB*tan(alpha))/(Ix-Ixy*tan(alpha))

+printf('part (a)')

+printf('\n Sigma A = %.1f M Pa \n',S_A/10^6)

+printf(' Sigma B = %.1f M Pa',S_B/10^6)

+

+// part (b)

+th=1/2*atan(-2*Ixy/(Ix-Iy))

+th1=0.415 //rad

+th2=-1.156 //rad

+IX=Ix*(cos(th1))^2+Iy*(sin(th1))^2-2*Ixy*sin(th1)*cos(th1)

+IY=Ix+Iy-IX

+Phi=2*%pi/3-th1

+alphA=-IX/(IY*tan(Phi))

+alpha=alphA+th1

+XA=xA*cos(th1)+yA*sin(th1)

+YA=yA*cos(th1)-xA*sin(th1)

+XB=xB*cos(th1)+yB*sin(th1)

+YB=yB*cos(th1)-xB*sin(th1)

+MX=M*sin(Phi)

+MY=-M*cos(Phi)

+S_A=MX*YA/IX-MY*XA/IY

+S_B=MX*YB/IX-MY*XB/IY

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

+printf('\n Sigma A = %.1f M Pa \n',S_A/10^6)

+printf(' Sigma B = %.1f M Pa',S_B/10^6)

diff --git a/3392/CH7/EX7.5/Ex7_5.sce b/3392/CH7/EX7.5/Ex7_5.sce
new file mode 100755
index 000000000..1512877ab
--- /dev/null
+++ b/3392/CH7/EX7.5/Ex7_5.sce
@@ -0,0 +1,31 @@
+clc

+// initialization of variables

+clear

+A=3085.9 //mm^2

+Ix=29.94e-6 //m^4

+Iy=4.167e-6 //m^4

+Ixy=0

+ybar=207.64 //mm

+tau_max=165e6 //Pa

+//calculations

+A=A*1e-6

+ybar=ybar*1e-3

+Mxk=-6.1*cos(%pi/6) // Mx=Mxk*P

+Myk=-6.1*sin(%pi/6) //My=Myk*P

+// Equation to be followed

+// S_zz=Mx*y/Ix-My*x/Iy

+// At A x=100 mm  y=-92.36 mm 

+x=100e-3

+y=-92.36e-3

+S_zzA=Mxk*y/Ix-Myk*x/Iy //Sigma_zz=S_zz*P

+// At B x=-100 mm  y=-92.36 mm 

+x=-100e-3

+y=-92.36e-3

+S_zzB=Mxk*y/Ix-Myk*x/Iy //Sigma_zz=S_zz*P

+// At C x=-3.125 mm  y=207.64 mm 

+x=-3.125e-3

+y=207.64e-3

+S_zzC=Mxk*y/Ix-Myk*x/Iy //Sigma_zz=S_zz*P

+// To find P

+P=2*tau_max/max(S_zzA,S_zzB,S_zzC)

+printf('P = %.2f kN',P/10^3)

diff --git a/3392/CH7/EX7.6/Ex7_6.sce b/3392/CH7/EX7.6/Ex7_6.sce
new file mode 100755
index 000000000..df799b495
--- /dev/null
+++ b/3392/CH7/EX7.6/Ex7_6.sce
@@ -0,0 +1,29 @@
+clc

+// initialization of variables

+clear

+P=35 //kN

+Phi=5*%pi/9

+E=72e9 //Pa

+L=3 //m

+Ix=39.69*10^6 //mm^4

+Iy=30.73*10^6 //mm^4

+Ixy=0

+//calculations

+P=P*1e3

+Ix=Ix*10^-12

+Iy=Iy*10^-12

+alpha=atan(-Ix/(Iy*tan(Phi)))

+M=P*L/4

+Mx=M*sin(Phi)

+yA=-118*10^-3 //m

+xA=70*10^-3 //m

+xB=-xA

+yB=82*10^-3 //m

+S_comp=Mx*(yA-xA*tan(alpha))/(Ix-Ixy*tan(alpha))

+S_tens=Mx*(yB-xB*tan(alpha))/(Ix-Ixy*tan(alpha))

+printf(' Tensile strength = %.1f M Pa \n',S_tens/10^6)

+printf(' Compressive Strength = %.1f M Pa',S_comp/10^6)

+v=P*L^3*sin(Phi)/(48*E*Ix)

+u=-v*tan(alpha)

+delta=sqrt(u^2+v^2)

+printf('\n The total deflection is %.2f mm',delta*10^3)

diff --git a/3392/CH7/EX7.7/Ex7_7.sce b/3392/CH7/EX7.7/Ex7_7.sce
new file mode 100755
index 000000000..40dc9e304
--- /dev/null
+++ b/3392/CH7/EX7.7/Ex7_7.sce
@@ -0,0 +1,29 @@
+clc

+// initialization of variables

+clear

+L=3 //m

+Ix=56.43e6 //mm^4

+Iy=18.11e6 //mm^4

+Ixy=22.72e6 //mm^4

+Phi=%pi/3

+E=200e9 //Pa

+Y=300e6 //Pa

+//calculations

+Ix=Ix*10^-12

+Iy=Iy*10^-12

+Ixy=Ixy*10^-12

+yA=-120*10^-3 //m

+xA=-91*10^-3 //m

+Nr=Ixy-Ix/tan(Phi)

+Dr=Iy-Ixy/tan(Phi)

+alpha=atan(Nr/Dr)

+// M=-L*P To know P we do the following

+Mxk=-L*sin(Phi)//Mx=Mxk*P

+P=Y*(Ix-Ixy*tan(alpha))/(Mxk*(yA-xA*tan(alpha)))

+printf('P = %.2f kN \n',P/10^3)

+v=P*L^3*sin(Phi)/(3*E*(Ix-Ixy*tan(alpha)))

+u=-v*tan(alpha)

+delta=sqrt(u^2+v^2)

+printf(' deflection = %.2f mm',delta*10^3)

+// Wrong calculation starting from v in Textbook

+

diff --git a/3392/CH7/EX7.8/Ex7_8.sce b/3392/CH7/EX7.8/Ex7_8.sce
new file mode 100755
index 000000000..3edb5283b
--- /dev/null
+++ b/3392/CH7/EX7.8/Ex7_8.sce
@@ -0,0 +1,26 @@
+clc

+// initialization of variables

+clear

+Ix=937e+06 //mm^4

+Iy=18.7e+6 //mm^4

+Ixy=0

+yA=305 //mm

+xA=90.5 //mm

+Phi=1.5533 //rad

+//calculations

+Ix=Ix*10^-12

+Iy=Iy*10^-12

+Ixy=Ixy*10^-12

+yA=yA*10^-3 //m

+xA=xA*10^-3 //m

+alpha=atan(-Ix/(Iy*tan(Phi)))

+Mxk=sin(Phi) // Mx=Mxk*M

+Sigma_Ak1=Mxk*(yA-xA*tan(alpha))/(Ix-Ixy*tan(alpha)) 

+//Sigma_A=Aigma_Ak*M

+// When the plane of the loads coincide with the y axes

+Sigma_Ak2=yA/Ix

+ratio=Sigma_Ak1/Sigma_Ak2

+percent=(ratio-1)*100

+printf('alpha = %.3f rad',alpha)

+printf('\n The maximum stress in the beam is increased %.1f percent when the plane of the loads is merely 1 degre from the symmetrical vertical plane',percent)

+// Wrong alpha given in the textbook

diff --git a/3392/CH7/EX7.9/Ex7_9.sce b/3392/CH7/EX7.9/Ex7_9.sce
new file mode 100755
index 000000000..b9d20c2fd
--- /dev/null
+++ b/3392/CH7/EX7.9/Ex7_9.sce
@@ -0,0 +1,20 @@
+clc

+// initialization of variables

+clear

+Y=280 //MPa

+AB=40 //mm

+BC=60 //mm

+//calculations

+Y=Y*10^6

+alpha=atan(BC/AB)

+C11=20/3 //mm

+C12=-10 //mm

+C21=-20/3 //mm

+C22=10 //mm

+Beta=atan((C11-C21)/(C22-C12))

+Phi=%pi/2+Beta

+d=sqrt((AB/2-C11)^2+(BC/2-C22)^2)

+d=d*10^-3 //m

+At=1/2*AB*BC/2*10^-6

+Mp=At*Y*d

+printf('Mp = %.3f kN.m',Mp/10^3)