13 files changed, 187 insertions, 0 deletions

diff --git a/3557/CH6/EX6.1/Ex6_1.sce b/3557/CH6/EX6.1/Ex6_1.sce
new file mode 100644
index 000000000..504333ec2
--- /dev/null
+++ b/3557/CH6/EX6.1/Ex6_1.sce
@@ -0,0 +1,14 @@
+//Example 6.1//

+

+s=300*10^6;//Pa //pascal //strain

+a=0.0043;// From the figure

+E1=s/a

+mprintf("E1 = %e GPa= 70GPa ",E1)

+mprintf("  (As G= 10^9)")

+//The 0.2% offset construction gives 

+mprintf("\nY.S. =410MPa")

+//The maximum for the stress stain curve gives

+mprintf("\n T.S = 480MPa")

+ef=0.08;//percent //the strain at fracture

+f=100*ef

+mprintf("\n f = %i percent",f)

diff --git a/3557/CH6/EX6.10/Ex6_10.sce b/3557/CH6/EX6.10/Ex6_10.sce
new file mode 100644
index 000000000..0b5c4fbef
--- /dev/null
+++ b/3557/CH6/EX6.10/Ex6_10.sce
@@ -0,0 +1,12 @@
+//Example 6.10//

+

+ap=5*10^-1;//percent per hour

+Q=2*10^5;//J/mol //activation energy

+R=8.314;//J/mol.K// universal gas constant

+T=1273;//K //Kelvin //absolute temperature

+T1=873;//given //absolute temperature

+C=ap*%e^((Q)/(R*T))

+mprintf("C = %e percent per hour",C)

+//applying this amount to the service temprature yield

+C1=C*%e^-((Q)/(R*T1))

+mprintf("\n C1 = %e percent per hour",C1) 

diff --git a/3557/CH6/EX6.11/Ex6_11.sce b/3557/CH6/EX6.11/Ex6_11.sce
new file mode 100644
index 000000000..1a9a824ff
--- /dev/null
+++ b/3557/CH6/EX6.11/Ex6_11.sce
@@ -0,0 +1,13 @@
+//Example 6.11//

+s=125;//ksi

+s=95;//ksi

+s=65;//ksi

+T=540;//degree C

+T=595;//degree C

+T=650;//degree C

+x=[540 595 650]

+y=[125 95 65 ]

+plot2d(x,y, style=1)

+ylabel("stress (ksi)","fontsize",2 );

+xlabel("T(degree C)")

+mprintf(" T = 585 degree C")

diff --git a/3557/CH6/EX6.12/Ex6_12.sce b/3557/CH6/EX6.12/Ex6_12.sce
new file mode 100644
index 000000000..eef02579e
--- /dev/null
+++ b/3557/CH6/EX6.12/Ex6_12.sce
@@ -0,0 +1,13 @@
+//Example 6.12//

+

+s1=2;//MPa //MegaPascal

+s2=1;//MPa //Megapascal

+a=60;//days //relaxation time for a rubber band at 25 degree C

+t=(a)*log(s1/s2)

+mprintf("t = %f days",t)

+Q=30*10^3;//J/mol //activation energy for the relaxation process

+R=8.314;//J/(mol.K) // universal gas constant

+T1=308;//K //Kelvin //absolute temperature

+T2=298;//K //Kelvin //absolute temperature

+t35=a*exp((Q/R)*((1/T1)-(1/T2)))

+mprintf("\n t35 = %f days",t35)

diff --git a/3557/CH6/EX6.13/Ex6_13.sce b/3557/CH6/EX6.13/Ex6_13.sce
new file mode 100644
index 000000000..e602d2229
--- /dev/null
+++ b/3557/CH6/EX6.13/Ex6_13.sce
@@ -0,0 +1,38 @@
+//Example 6.13//

+a=514;//K //Kelvin //Temperature

+b=273;//K //Kelvin //Temperature

+apt=a+b

+mprintf("apt = %i K for eta = 10^13.4P",apt)

+c=696;//K //Kelvin //Temperature

+spt=c+b//for eta=10^7.6P

+mprintf("\n spt = %i K",spt)

+i=(10^13.4); //P //Pascal //preexponential constant

+j=(10^7.6);//P // Pascal //preexponential constant

+f=8.314;//J/(mol K) //universal gas constant

+a1=log(i/j); //(Taking antilog of i and j to remove exponential term)

+//mprintf("\na1 = %f ",a1)

+b1=(1/apt)-(1/spt);//(subtracting the temperature)

+//mprintf("\nb1 = %e ",b1)

+Q=(a1/b1)*f

+mprintf("\nQ = %e J/mol  = 465kJ (As 1K = 10^3)",Q)

+eta0=i*%e^-((Q)/(f*apt))

+mprintf("\n eta0 = %e P",eta0)

+h=10^4;//given

+//for eta=10^4 P and eta=10^8 P

+//for eta = 10^4

+T=Q/((f)*log(h/eta0))

+mprintf("\n T = %i K = 858 degree C",T)

+//for eta=10^8P

+h1=10^8;//P //Pascal 

+T1=Q/((f)*log(h1/eta0))

+mprintf("\n T1 = %i K = 680 degree C",T1)

+//Therefore working range = 680 to 858 degree C

+

+//For melting range eta=50 to 500 P

+ eta=50;//P //Pascal

+T2=Q/((f)*log(eta/eta0))

+mprintf("\n T2 = %i K = 993 degree C",T2)

+ eta1 = 500;// P //Pascal

+T3=Q/((f)*log(eta1/eta0))

+mprintf("\n T3 = %i K = 931 degree C",T3)

+mprintf("\n melting range = 931 to 993 degree C")

diff --git a/3557/CH6/EX6.2/Ex6_2.sce b/3557/CH6/EX6.2/Ex6_2.sce
new file mode 100644
index 000000000..7adbb0732
--- /dev/null
+++ b/3557/CH6/EX6.2/Ex6_2.sce
@@ -0,0 +1,12 @@
+//Example 6.2//

+p=50000;//N //tensile load

+A0=5*10^-3;//m //area of the sample parallel to the applied load

+s=p/(%pi*A0^2)

+mprintf("s = %e N/m^2  637 MPa",s)

+mprintf("  (As M= 10^6)")

+s1=637*10^6;//Pa //Pascal //modulus of elasticity

+E=200*10^9;//Pa // Pascal //Youngs Modulus

+E1=s1/E

+mprintf("\n E1 = %e ",E1)

+

+

diff --git a/3557/CH6/EX6.3/Ex6_3.sce b/3557/CH6/EX6.3/Ex6_3.sce
new file mode 100644
index 000000000..2877cf472
--- /dev/null
+++ b/3557/CH6/EX6.3/Ex6_3.sce
@@ -0,0 +1,23 @@
+//Example 6.3//

+

+P=6*10^3//N //Newton // load on the sample

+A0=(10/2)*10^-3;//N/m^2

+s=P/(%pi*A0^2)

+mprintf("s = %e N/m^2 = 76.4 MPa",s)

+mprintf("  (As M= 10^6)")

+s1=76.4;//MPa //Megapascal //modulus od elasticity

+E=70*10^3;//MPa//Megapascal //Young's Modulus

+e=s1/E

+mprintf("\n e = %e",e)

+//the strain of diameter is calculated as

+v=0.33;//given

+ed=-v*e

+mprintf("\n ed = %e ",ed)

+//resulting diameter

+d0=10;//mm 

+df=d0*(ed+1)

+mprintf("\n df = %f mm",df)

+//compressive stress

+ed1=+3.60*10^-4;//the diameter strain will be of equal magnitude but opposite sign

+df1=d0*(ed1+1);

+mprintf("\n df1 = %f mm",df1)

diff --git a/3557/CH6/EX6.4/Ex6_4.sce b/3557/CH6/EX6.4/Ex6_4.sce
new file mode 100644
index 000000000..bea1a2bd2
--- /dev/null
+++ b/3557/CH6/EX6.4/Ex6_4.sce
@@ -0,0 +1,11 @@
+//Example 6.4//

+

+rO=0.132;//nm //Ionic radius of Oxygen (From appendix 2)

+p=2*rO

+mprintf("p = %f nm",p)

+a=7.0*10^9;//Pa //The theoretical strength of the defect free glass

+p1=0.264*10^-9//m  

+c=1*10^-6;//m //crack length

+s=(1/2)*a*(p1/c)^(1/2)

+mprintf("\ns = %e Mpa  = 57MPa (As M =10^6)",s)

+

diff --git a/3557/CH6/EX6.5/Ex6_5.sce b/3557/CH6/EX6.5/Ex6_5.sce
new file mode 100644
index 000000000..8244baad5
--- /dev/null
+++ b/3557/CH6/EX6.5/Ex6_5.sce
@@ -0,0 +1,8 @@
+//Example 6.5//

+L=50*10^-3;//m //Distance between support

+m=404*10^3;//N/m //Initial slope of load-deflection curve

+b=13*10^-3;//m //test piece geometry

+h=7*10^-3;//m //test piece geometry

+E=((L^3)*m)/(4*b*h^3)

+mprintf("E = %e N/m^2  =2830MPa  (As M= 10^6)",E)

+

diff --git a/3557/CH6/EX6.6/Ex6_6.sce b/3557/CH6/EX6.6/Ex6_6.sce
new file mode 100644
index 000000000..79144c479
--- /dev/null
+++ b/3557/CH6/EX6.6/Ex6_6.sce
@@ -0,0 +1,15 @@
+//Example 6.6//

+

+E=830;//MPa //Megapascal //Young's Modulus

+s=1;//MPa //MegaPascal //modulus of elasticity

+e=s/E

+mprintf("e = %e",e)

+//(b)

+E1=1.3;//MPa//Megapascal //Young's Modulus

+e1=s/E1

+mprintf("\n e1 = %f",e1)

+//(c)E=200 GPa= 2*10^5 Mpa

+E2=2*10^5;//MPa//Megapascal //Young's Modulus

+e2=s/E2//Mpa

+mprintf("\n e2 = %e",e2)

+

diff --git a/3557/CH6/EX6.7/Ex6_7.sce b/3557/CH6/EX6.7/Ex6_7.sce
new file mode 100644
index 000000000..8d018d6aa
--- /dev/null
+++ b/3557/CH6/EX6.7/Ex6_7.sce
@@ -0,0 +1,11 @@
+//Example 6.7//

+

+//From hooke law'

+s=0.2489;//nm //nanometer // modulus of elasticity

+s1=0.2480;// nm //nanometer // modulus of elasticity

+e=(s-s1)/s1

+printf("e = %f ",e)

+s2=1000;//Mpa //MegaPascal //sigma

+E=s2/e

+mprintf("\n E = %e",E)

+mprintf("  275   GPa  (As G=10^9) (Answer calculated in the textbook is wrong)")

diff --git a/3557/CH6/EX6.8/Ex6_8.sce b/3557/CH6/EX6.8/Ex6_8.sce
new file mode 100644
index 000000000..9e87eaa2b
--- /dev/null
+++ b/3557/CH6/EX6.8/Ex6_8.sce
@@ -0,0 +1,9 @@
+//Example 6.8//

+si=0.690;//MPa //Megapascal //tensile stress

+a=cosd(40);//degree

+b=cosd(60);//degree

+torque=si*a*b

+mprintf("torque = %f MPa (38.3psi)",torque)

+t=0.94;//MPa //MegaPascal //torque

+sig=t/(a*b)

+mprintf("\n sig = %f Mpa (356psi)",sig)

diff --git a/3557/CH6/EX6.9/Ex6_9.sce b/3557/CH6/EX6.9/Ex6_9.sce
new file mode 100644
index 000000000..e6de54b79
--- /dev/null
+++ b/3557/CH6/EX6.9/Ex6_9.sce
@@ -0,0 +1,8 @@
+//Example6.9//

+P=3000;//kg //load

+D=10;//mm//diamter sphere of tungsten carbide

+d=3.91;//mm //diameter impression in the iron surface

+BHN=(2*P)/((%pi*D)*(D-sqrt(D^2-d^2)))

+mprintf("BHN = %i",BHN)

+//From the Figure 6.28b

+printf("\n(TS)BHN=240 = 800 Mpa")