initial commit / add all books

6 files changed, 232 insertions, 0 deletions

diff --git a/3289/CH4/EX4.1/Ex4_1.sce b/3289/CH4/EX4.1/Ex4_1.sce
new file mode 100755
index 000000000..40f90ee61
--- /dev/null
+++ b/3289/CH4/EX4.1/Ex4_1.sce
@@ -0,0 +1,47 @@
+clc;

+

+sigmayp=350 //MPa

+sigma3=0

+M=8 //kN

+Mt=24 //kNm

+N=2

+v=0.3

+

+// sigma= My/I ==32M/%pid^3

+// tau= Mt*r/J ==16Mt/%pid^3

+//sigma1=(16*(M+sqrt(M^2+Mt^2)))/(%pi*d^3)

+//sigma2=(16*(M-sqrt(M^2+Mt^2)))/(%pi*d^3)

+

+//solution a: max principal stress theory

+//(16*(M+sqrt(M^2+Mt^2)))/(%pi*d^3)=sigmayp/N

+

+a=(16*(M+sqrt(M^2+Mt^2)))/%pi

+disp(a)

+b=sigmayp*10^6/N

+disp(b)

+d=(a/b)^(1/3)

+disp(d,"diameter of the bar in meter is= ")

+

+//solution b:max shearing stress theory

+

+c=(32*sqrt(M^2+Mt^2))/%pi

+disp(c)

+d=(c/b)^(1/3)

+disp(d,"diameter of the bar in meter is= ")

+

+//solution c:max principal strain theory

+//epsilon1=(sigma1-v(sigma2+sigma3))/E=epsilonyp/N=sigmayp/EN    Or

+//sigma1-v(sigma2+sigma3)=b    given

+//sigma1=b+v(sigma2+sigma3)    substituting the values of the sigma 1,2 and 3 we get

+//(16*(M+sqrt(M^2+Mt^2)-vM-v*sqrt(M^2+Mt^2)))/(%pi*d^3)=b

+e=(16*(M+sqrt(M^2+Mt^2)-v*M-v*sqrt(M^2+Mt^2)))/%pi

+disp(e)

+d=(e/b)^(1/3)

+disp(d,"diameter of the bar in meter is= ")

+

+//solution d:max energy of distortion theory

+

+f=(16*sqrt(4*M^2+3*Mt^2))/%pi

+disp(f)

+d=(f/b)^(1/3)

+disp(d,"diameter of the bar in meter is= ")

diff --git a/3289/CH4/EX4.3/EX4_3.sce b/3289/CH4/EX4.3/EX4_3.sce
new file mode 100755
index 000000000..de12c5daf
--- /dev/null
+++ b/3289/CH4/EX4.3/EX4_3.sce
@@ -0,0 +1,43 @@
+clc

+

+sigmau1=300  //MPa

+sigmau2=700  //MPa

+b=0.105 //m outer diameter

+a=0.100 //m inner diameter

+

+//sigma1==(-sigma2)==tau

+

+sigma3=0

+

+//Mt=J*tau/r= (%pi*(b^2-a^2))/2*b

+//Mt=((%pi*(b^4-a^4))/(2*b))*tau         equation a

+q=(%pi*(b^4-a^4))/(2*b)

+

+//solution a: max principal stress theory

+tau=sigmau1

+Mt=q*tau*10^6

+disp(Mt,"max principal stress in Nm is= ")

+

+//solution b:max shearing stress theory

+// |sigma1-sigma2|=sigmau1

+// 2*sigma1==sigmau1==2*tau  or

+

+Mt=q*tau*10^6

+disp(Mt,"max shearing stress in Nm is= ")

+

+//solution c:Coulomb- mohr theory

+//(tau/sigmau1)-(-tau/sigmau2)=1

+tau=1*((sigmau1*sigmau2)/(sigmau1+sigmau2))

+disp(tau,"tau in MPa is= ")

+Mt=q*tau*10^6

+disp(Mt,"Coulomb- mohr in Nm is= ")

+

+

+

+

+

+

+

+

+

+

diff --git a/3289/CH4/EX4.4/Ex4_4.sce b/3289/CH4/EX4.4/Ex4_4.sce
new file mode 100755
index 000000000..0df8e6719
--- /dev/null
+++ b/3289/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,22 @@
+clc

+

+a=0.05 // m

+Fm=90 //kN

+sigmacr=210 // MPa

+sigmayp=280 // MPa

+

+//sigmaa=Ma*c/I         equation 1

+//Ma=0.025*Fa

+c=0.025

+I=(a^4)/12

+disp(I)

+//sigmaa=((0.025*Fa)*c)/I  substituting the values

+

+

+//sigmam=Fm/A           equation 2

+sigmam=Fm/(a*a)

+disp(sigmam,"in kilo Pa is= ")

+

+//(((1200*Fa)/sigmacr)+(sigmam/sigmayp))=1

+Fa=(1-(sigmam/sigmayp))*(sigmacr/1200)

+disp(Fa,"load Fa in N is= ")                         //wrong ans in textbook

diff --git a/3289/CH4/EX4.5/Ex4_5.sce b/3289/CH4/EX4.5/Ex4_5.sce
new file mode 100755
index 000000000..f785f9700
--- /dev/null
+++ b/3289/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,27 @@
+clc

+

+r=0.04 //m

+t=5 //mm

+sigmae=250 //MPa

+sigmay=300 //MPa

+

+//sigmathetamax=(p*r)/t =8*p        max values of tangential stresses

+//sigmathetamin=((-p/4)*r)/t =-2*p   min values of tangential stresses

+//sigmazamax=(p*r)/2*t =4*p       axial principal stresses

+//sigmazmin=((-p/4)*r)/2*t =-p

+

+//sigmathetaa=(sigmathetamax-sigmathetamin)/2= 5p    alternating and mean stresses

+//sigmathetam=(sigmathetamax+sigmathetamin)/2= 3p

+//sigmaza=(sigmazamax-sigmazmin)/2 =2.5p

+//sigmazm=(sigmazamax+sigmazmin)/2 =1.5p

+

+//sqrt(sigmathetaa^2-sigmathetaa*sigmaza+sigmaza^2)=sigmaea

+//sqrt(sigmathetam^2-sigmathetam*sigmazm+sigmazm^2)=sigmaem

+

+//sqrt(25p^2-12.3p^2+6.25p^2)=sigmaea

+//sqrt(9p^2-4.5p^2+2.25p^2)=sigmaem           solving this equation we get

+sigmaea=4.33 //p

+sigmaem=2.60 //p

+

+p=1/((sigmaea/sigmae)+(sigmaem/sigmay))

+disp(p,"the value of p in MPa is= ")

diff --git a/3289/CH4/EX4.6/Ex4_6.sce b/3289/CH4/EX4.6/Ex4_6.sce
new file mode 100755
index 000000000..bf4bdc443
--- /dev/null
+++ b/3289/CH4/EX4.6/Ex4_6.sce
@@ -0,0 +1,50 @@
+clc

+

+a=[700 14 0; 14 -350 0; 0 0 -350]

+disp(a)

+c=[-660 -7 0; -7 -350 0; 0 0 -350]

+disp(c)

+sigmau=2400 //MPa

+K=1

+sigmae=800 //MPa

+Nf=1 //cycles for SAE

+Nff=10^3 //cycles for Gerber

+Ne=10^8 //cycles

+

+sigmaxa=(700+660)/2

+disp(sigmaxa,"alternating and mean values of stresses in MPa is= ")

+sigmaxm=(700-660)/2

+disp(sigmaxm,"alternating and mean values of stresses in MPa is= ")

+sigmaya=(-350+350)/2

+disp(sigmaya,"alternating and mean values of stresses in MPa is= ")

+sigmaym=(-350-350)/2

+disp(sigmaym,"alternating and mean values of stresses in MPa is= ")

+sigmazm=(-350-350)/2

+disp(sigmazm,"alternating and mean values of stresses in MPa is= ")

+tauxya=(14+7)/2

+disp(tauxya,"alternating and mean values of stresses in MPa is= ")

+tauxym=(14-7)/2

+disp(tauxym,"alternating and mean values of stresses in MPa is= ")

+

+sigmaea=sqrt(((sigmaxa-sigmaya)^2+(sigmaya-sigmaxa)^2+6*(tauxya)^2)/2)

+disp(sigmaea,"in MPa is =")

+sigmaem=sqrt(((sigmaxm-sigmaym)^2+(sigmaym-sigmaxm)^2+6*(tauxym)^2)/2)

+disp(sigmaem,"in MPa is= ")

+

+//solution a: 

+sigmacr=sigmaea/(1-(sigmaem/2400))

+disp(sigmacr)

+b=log(sigmau/sigmae)/log(1/Ne)

+disp(b)

+

+Ncr=1*(sigmacr/2400)^(1/b)

+disp(Ncr,"in cycles is= ")

+

+//solution b:

+sigmacr=sigmaea/(1-(sigmaem/sigmau)^2)

+disp(sigmacr,"in MPa is= ")

+b=log(0.9*2400/sigmae)/log(Nff/Ne)

+disp(b)

+

+Ncr=Nff*(sigmacr/(0.9*2400))^(-11.587)

+disp(Ncr,"in cycles is= ")

diff --git a/3289/CH4/EX4.7/Ex4_7.sce b/3289/CH4/EX4.7/Ex4_7.sce
new file mode 100755
index 000000000..1be3ce0ba
--- /dev/null
+++ b/3289/CH4/EX4.7/Ex4_7.sce
@@ -0,0 +1,43 @@
+clc

+

+W=180 //N

+h=0.1 //m

+L=1.16 //m

+w=0.025 //m

+d=0.075 //m

+E=200 //GPa

+k=180 //kN/m

+

+I=w*d^3

+disp(I)

+//deltast=(W*L^3)/(48*E*I)       equation 1

+deltast=(W*L^3*12)/(48*E*10^9*I)

+disp(deltast,"deflection of a point in meter is= ")

+

+//deltastmax=Mc/I                equation 2

+deltastmax=(W*L*12*0.0375)/(4*I)

+disp(deltastmax,"deflection of a point in Pa is= ")

+

+//solution a:

+a=1+sqrt(1+((2*h)/deltast))

+disp(a,"imapct factor is= ")

+deltamax=deltast*a

+disp(deltamax,"in meter is =")

+sigmamax=deltastmax*a

+disp(sigmamax,"in Pa is= ")

+

+//solution b:

+deltast=deltast+(90/180000)

+disp(deltast,"static deflection of the beam in meter is= ")

+a=1+sqrt(1+((2*h)/deltast))

+disp(a,"imapct factor is= ")

+deltamax=deltast*a

+disp(deltamax,"in meter is =")

+sigmamax=deltastmax*a

+disp(sigmamax,"in Pa is= ")

+

+

+

+

+

+