diff options
Diffstat (limited to '3557/CH8')
-rw-r--r-- | 3557/CH8/EX8.2/Ex8_2.sce | 9 | ||||
-rw-r--r-- | 3557/CH8/EX8.3/Ex8_3.sce | 11 | ||||
-rw-r--r-- | 3557/CH8/EX8.4/Ex8_4.sce | 6 | ||||
-rw-r--r-- | 3557/CH8/EX8.5/Ex8_5.sce | 15 | ||||
-rw-r--r-- | 3557/CH8/EX8.6/Ex8_6.sce | 6 | ||||
-rw-r--r-- | 3557/CH8/EX8.7/Ex8_7.sce | 16 |
6 files changed, 63 insertions, 0 deletions
diff --git a/3557/CH8/EX8.2/Ex8_2.sce b/3557/CH8/EX8.2/Ex8_2.sce new file mode 100644 index 000000000..46d7adf89 --- /dev/null +++ b/3557/CH8/EX8.2/Ex8_2.sce @@ -0,0 +1,9 @@ +//Example 8.2//
+
+Y=1;// dimensionless geometry factor
+YS=1460//MPa //MegaPascal // overall stress applied at failure
+b=0.5;//Y.S //given
+Kic=98;//MPa sqrt(m) //fracture toughness
+a=(Kic^2)/((%pi)*(b*YS)^2)
+mprintf("a = %e m = 5.74 mm (As 1milli = 10^-3 )",a)
+
diff --git a/3557/CH8/EX8.3/Ex8_3.sce b/3557/CH8/EX8.3/Ex8_3.sce new file mode 100644 index 000000000..36b09377f --- /dev/null +++ b/3557/CH8/EX8.3/Ex8_3.sce @@ -0,0 +1,11 @@ +//Example 8.3//
+
+a=25*10^-6;// m // length of surface crack
+// (a) For Sic,
+b=3;//MPa sqrt(m) //fracture toughness
+s1=b/(sqrt(%pi*a))
+mprintf("s1 = %i MPa",s1)
+// (b) For PSZ,
+c=9;//MPa sqrt(m)// fracture toughness
+s2=c/(sqrt(%pi*a))
+mprintf("\n s2 = %i MPa (Answer calculated in textbook is wrong)",s2)
diff --git a/3557/CH8/EX8.4/Ex8_4.sce b/3557/CH8/EX8.4/Ex8_4.sce new file mode 100644 index 000000000..e2566e4b4 --- /dev/null +++ b/3557/CH8/EX8.4/Ex8_4.sce @@ -0,0 +1,6 @@ +//Example 8.4//
+T.S=800;//MPa
+F.S=T.S/4
+mprintf("F.S = %i MPa",F.S)
+ss=F.S/2
+mprintf("\n ss = %i Mpa",ss)
diff --git a/3557/CH8/EX8.5/Ex8_5.sce b/3557/CH8/EX8.5/Ex8_5.sce new file mode 100644 index 000000000..51098a279 --- /dev/null +++ b/3557/CH8/EX8.5/Ex8_5.sce @@ -0,0 +1,15 @@ +//Example8.5//
+Q=78.6*10^3;//J/mol //Activation energy
+R=8.314;//J/mol //universal gas constant
+T=323;//K //Kelvin //absolute temperature
+T1=223;//K //Kelvin //absolute temperature
+C=1/(%e^-((Q)/(R*T)))
+mprintf("C = %e s^-1",C)
+t50=C*(%e^-(Q/(R*T1)))
+mprintf("\n t50 = %e s^-1",t50)
+t=5.0*10^5;//s //seconds
+a=1;//h //hour
+b=3.6*10^3;//s //seconds
+t1=t*(a/b)
+mprintf("\n t1 = %i h =5days, 20h (Answer calculated in the textbook is wrong)",t1)
+
diff --git a/3557/CH8/EX8.6/Ex8_6.sce b/3557/CH8/EX8.6/Ex8_6.sce new file mode 100644 index 000000000..fd599a804 --- /dev/null +++ b/3557/CH8/EX8.6/Ex8_6.sce @@ -0,0 +1,6 @@ +//Example 8.6//
+
+u=0.293;//mm^-1 //linear absorption coefficient for the material
+x=10;//mm //x-ray beam intensity transmitted
+I=%e^-(u*x)
+disp(I)
diff --git a/3557/CH8/EX8.7/Ex8_7.sce b/3557/CH8/EX8.7/Ex8_7.sce new file mode 100644 index 000000000..f72128898 --- /dev/null +++ b/3557/CH8/EX8.7/Ex8_7.sce @@ -0,0 +1,16 @@ +//Example 8.7//
+//(a)
+a=2.70;//Mg/m^3 //Density of aluminium (From appendix 1)
+b=6320;//m/s //velocity of sound
+ZAl=a*b
+mprintf("ZAl = %e Mg/(m^2s)",ZAl)
+a1=7.85;//Mg/m^3 //Density of Manganese (From Appendix 1)
+b1=5760;//m/s //Velocity of sound
+Zst=a1*b1
+mprintf("\n Zst = %e Mg/(m^2s)",Zst)
+Ir=[(Zst-ZAl)/(Zst+ZAl)]^2
+mprintf("\n Ir = %f ",Ir)
+
+//(b) For the reverse direction of ultrasonic-pulse travel
+Ir1=[(ZAl-Zst)/(ZAl+Zst)]^2
+mprintf("\n Ir1 = %f ",Ir1)
|