initial commit / add all books

9 files changed, 175 insertions, 0 deletions

diff --git a/866/CH14/EX14.1/14_1.sce b/866/CH14/EX14.1/14_1.sce
new file mode 100755
index 000000000..7ead5d32e
--- /dev/null
+++ b/866/CH14/EX14.1/14_1.sce
@@ -0,0 +1,19 @@
+clc

+//initialisation of variables

+d= 2 //m

+t= 20 //mm

+p= 1.5 //N/mm^2

+load= 2500 //KN

+alpha= 60 //degrees

+//CALCULATIONS

+Cs= (p*d*10^3)/(2*t)

+Ls= (p*d*10^3)/(4*t)

+Ds= (load*10^3)/(%pi*d*t*10^3)

+Ts= Ls+Ds

+sigman= (Ts*t*(cosd(90-alpha))^2+Cs*t*(cosd(alpha))^2)/t

+tab= Ts*sind(alpha)*cosd(alpha)-Cs*sind(alpha)*cosd(alpha)

+tmax= (Ts-Cs)/2

+//RESULTS

+printf ('direct stress= %.1f N/mm^2',sigman)

+printf (' \n Shear stress=%.1f N/mm^2',tab)

+printf (' \n maximum Shear stress=%.1f N/mm^2',tmax)

diff --git a/866/CH14/EX14.10/14_10.sce b/866/CH14/EX14.10/14_10.sce
new file mode 100755
index 000000000..9358e82f0
--- /dev/null
+++ b/866/CH14/EX14.10/14_10.sce
@@ -0,0 +1,13 @@
+clc

+//initialisation of variables

+BeamL= 250 //mm

+BeamH= 500 //mm

+spanlength= 4 //m

+Failurestress= 1.5 //N/mm^2

+//CALCULATIONS

+W1= (Failurestress*BeamL*BeamH^3)/(10^3*BeamL*12)

+W2= (4*Failurestress*BeamL*BeamH)/3

+if(W1>W2)

+    printf ("W2= %.2e kN",W2)

+else

+    printf ("W2= %.2e kN",W1)

diff --git a/866/CH14/EX14.2/14_2.sce b/866/CH14/EX14.2/14_2.sce
new file mode 100755
index 000000000..12ff51aff
--- /dev/null
+++ b/866/CH14/EX14.2/14_2.sce
@@ -0,0 +1,19 @@
+clc

+//initialisation of variables

+load= 50000 //N

+torque= 1200 //Nm

+d= 60//mm

+t= 1.5 //mm

+alpha= 60 //degrees 

+//CALCULATIONS

+BM= load*t

+axialload= (load*4)/(%pi*d^2)

+bendingmoment= (BM*d*64)/(%pi*d^4*2)

+Ts= axialload+bendingmoment

+shearstress= (torque*10^3*d*32)/(2*%pi*d^4)

+sigman= -Ts*(cosd(alpha-30))^2+shearstress*cosd(alpha-30)*sind(alpha-30)+shearstress*cosd(alpha-30)*sind(alpha-30)

+T= -Ts*sind(alpha)*cosd(alpha)-shearstress*(sind(alpha))^2+shearstress*(cosd(alpha)^2)

+//RESULTS

+printf ('direct stress= %.1f N/mm^2',sigman)

+printf (' \n Shear stress=%.1f N/mm^2',T)

+

diff --git a/866/CH14/EX14.3/14_3.sce b/866/CH14/EX14.3/14_3.sce
new file mode 100755
index 000000000..af9f5eecd
--- /dev/null
+++ b/866/CH14/EX14.3/14_3.sce
@@ -0,0 +1,24 @@
+clc

+//initialisation of variables

+directstress= 80 //N.mm^2

+shearstress= 45 //N/mm^2

+//CALCULATIONS

+Sx= directstress 

+Sy= 0

+Txy= shearstress

+sigma1= ((Sx+Sy)/2)+(sqrt((Sx-Sy)^2+4*Txy^2))/2

+sigma2= ((Sx+Sy)/2)-(sqrt((Sx-Sy)^2+4*Txy^2))/2

+theta= (atand(-((2*Txy)/(Sx-Sy))))/2

+direction1= theta

+direction2= theta-90

+Tmax= (sigma1-sigma2) /2

+theta1= direction1-45

+theta2= direction2-45

+//RESULTS

+printf ('principal stress1= %.1f N/mm^2',sigma1)

+printf (' \n principal stress2=%.1f N/mm^2',sigma2)

+printf (' \n direction of principal stress1=%.1f degrees',direction1)

+printf (' \n.direction of principal stress2=%.1f degrees',direction2)

+printf (' \n.maximum shear stress=%.1f N/mm^2',Tmax)

+printf (' \n.derection of plane1=%.1f degrees',theta1)

+printf (' \n.derection of plane2=%.1f degrees',theta2)

diff --git a/866/CH14/EX14.5/14_5.sce b/866/CH14/EX14.5/14_5.sce
new file mode 100755
index 000000000..d2abc4e3b
--- /dev/null
+++ b/866/CH14/EX14.5/14_5.sce
@@ -0,0 +1,24 @@
+clc

+//initialisation of variables

+sigmax= 60 //N/mm^2

+sigmay= -40 //N/mm^2

+Txy= 50 //N/mm^2

+E= 200000 //N/mm^2

+mu= 0.3

+//CALCULATIONS

+Sx= (sigmax-mu*sigmay)/E

+Sy= (sigmay-mu*sigmax)/E

+G= E/(2*(1+mu))

+gamaxy= Txy/G

+e1= ((Sx+Sy)/2+(sqrt((Sx-Sy)^2+(gamaxy)^2))/2)

+e2= ((Sx+Sy)-(sqrt((Sx-Sy)^2+gamaxy^2)))/2

+theta= atand(-gamaxy/(Sx-Sy))/2

+theta1= -90+theta

+//RESULTS

+printf ('direct strain in x direction= %.2e ',Sx)

+printf (' \n direct strain in y direction=%.2e ',Sy)

+printf (' \n Gamaxy=%.2e ',gamaxy)

+printf (' \n.principal strain1=%.9e ',e1)

+printf (' \n.principal strain2=%.9e ',e2)

+printf (' \n.direction of plane1=%.1f degrees',theta)

+printf (' \n.direction of plane2=%.1f degrees',theta1)

diff --git a/866/CH14/EX14.6/14_6.sce b/866/CH14/EX14.6/14_6.sce
new file mode 100755
index 000000000..22f8b2afc
--- /dev/null
+++ b/866/CH14/EX14.6/14_6.sce
@@ -0,0 +1,21 @@
+clc

+//initialisation of variables

+ea= 1000*10^-6

+eb= -200*10^-6

+ec= -300*10^-6

+sigmay= 0 //N/mm^2

+E= 70000 //N/mm^2

+mu= 0.3

+d= 50 //mm

+//CALCULATIONS

+e1= ((ea+ec)/2)+sqrt(((ea-eb)^2+(ec-eb)^2)/2)

+e2= ((ea+ec)/2)-sqrt(((ea-eb)^2+(ec-eb)^2)/2)

+sigma1= (E*(e1+mu*e2))/(1-mu^2)

+sigma2= (E*(e2+mu*e1))/(1-mu^2)

+sigmax= sigma1+sigma2

+Txy= sqrt((((sigma1*2)-(sigmax))^2-(sigmax^2)))/2

+P= sigmax*((%pi/4)*(d^2))

+T= Txy*((%pi/32)*(50^4)/25)

+//RESULTS

+printf ('T= %.2fNm',T) 

+printf ('\n P= %.2fN',P) 

diff --git a/866/CH14/EX14.7/14_7.sce b/866/CH14/EX14.7/14_7.sce
new file mode 100755
index 000000000..32d194fe7
--- /dev/null
+++ b/866/CH14/EX14.7/14_7.sce
@@ -0,0 +1,21 @@
+clc

+//initialisation of variables

+sigmax= 140 //N/mm^2

+sigmay= -70 //N/mm^2

+Txy= 60 //N/mm^2

+T= 225 //N/mm^2

+//CALCULATIONS

+sigma1= (sigmax+sigmay)/2+(sqrt((sigmax-sigmay)^2+4*(Txy^2)))/2

+sigma2= (sigmax+sigmay)/2-(sqrt((sigmax-sigmay)^2+4*(Txy^2)))/2

+sigmaY= sigma1-sigma2

+sigmaY1= sqrt(sigma1^2+sigma2^2-sigma1*sigma2)

+if(sigmaY>T)

+    disp("Tresca theory failure has ocurred")

+else

+    disp("Tresca theory failed")

+end

+if(sigmaY1<T)

+    disp("According to Von Mises theory material has not failed")

+else

+    disp("According to Von Mises theory material has failed")

+end

diff --git a/866/CH14/EX14.8/14_8.sce b/866/CH14/EX14.8/14_8.sce
new file mode 100755
index 000000000..388468b65
--- /dev/null
+++ b/866/CH14/EX14.8/14_8.sce
@@ -0,0 +1,19 @@
+clc

+//initialisation of variables

+M= 250 //KNm

+Torquemax= 200 //KNm

+allowablestress= 180 //N/mm^2

+L=250 //mm

+B= 500 //mm

+t= 10 //mm

+t1= 12 //mm 

+//CALCULATIONS

+Stressmax= (Torquemax*10^6)/(2*B*L*t)

+I= (2*t1*L*L^2)+((2*t*B^3)/12)

+sigma= (M*10^6*B)/(2*I)

+Stressallowable= sqrt(sigma^2+3*(Stressmax^2))

+if(Stressallowable<allowablestress)

+    disp("Box grider section is satisfactory")

+else

+    disp("Box grider section is  not satisfactory")

+end

diff --git a/866/CH14/EX14.9/14_9.sce b/866/CH14/EX14.9/14_9.sce
new file mode 100755
index 000000000..c1977dd17
--- /dev/null
+++ b/866/CH14/EX14.9/14_9.sce
@@ -0,0 +1,15 @@
+clc

+//initialisation of variables

+BeamL= 60 //mm

+BeamB= 100 //mm

+Load= 60000 //N

+Stress=  150 //N/mm^2

+//CALCULATIONS

+sigmax= Load/(BeamL*BeamB)

+stressmax1= sqrt((Stress/2)^2-(sigmax^2/4))

+Sy1= (2*stressmax1*BeamL*BeamB)/3

+stressmax2= sqrt((Stress^2-(sigmax^2/4))/3-(sigmax^2/4))

+Sy2= (2*stressmax2*BeamL*BeamB)/3  

+//RESULTS

+printf ('Shear force using Tresca= %.2e N',Sy1)

+printf (' \n Shear force using Von Misses=%.2e N',Sy2)