9 files changed, 200 insertions, 0 deletions
diff --git a/3763/CH12/EX12.1/Ex12_1.sce b/3763/CH12/EX12.1/Ex12_1.sce
new file mode 100644
index 000000000..7f5cc48d1
--- /dev/null
+++ b/3763/CH12/EX12.1/Ex12_1.sce
@@ -0,0 +1,15 @@
+clear
+//
+//
+//
+
+//Variable declaration
+sigma0=8.55
+K=2.45      
+sigma=10**-3      //steel size(mm)
+
+//Calculation
+sigma=sigma0+(K/sqrt(sigma))      //yield strength
+
+//Result
+printf("\n yield strength is %0.3f  kg/mm**2",sigma)
diff --git a/3763/CH12/EX12.2/Ex12_2.sce b/3763/CH12/EX12.2/Ex12_2.sce
new file mode 100644
index 000000000..f5aa02dfb
--- /dev/null
+++ b/3763/CH12/EX12.2/Ex12_2.sce
@@ -0,0 +1,15 @@
+clear
+//
+//
+//
+
+//Variable declaration
+E=70*10**9     //young's modulus(Pa)
+gama=1     //surface energy(joule/m**2)
+C=1*10**-6    //depth(m)
+
+//Calculation
+sigma_f=sqrt(2*E*gama/(%pi*C))      //fracture strength(GPa)
+
+//Result
+printf("\n fracture strength is %0.3f  GPa",sigma_f/10**9)
diff --git a/3763/CH12/EX12.3/Ex12_3.sce b/3763/CH12/EX12.3/Ex12_3.sce
new file mode 100644
index 000000000..2c2e8d605
--- /dev/null
+++ b/3763/CH12/EX12.3/Ex12_3.sce
@@ -0,0 +1,21 @@
+clear
+//
+//
+//
+
+//Variable declaration
+ml=57800     //load(N)
+d=10*10**-3    //diameter(m)
+D=5      //diameter after fracture(mm)
+l=50    //guage length(mm)
+L=55    //length after fracture(mm)
+
+//Calculation
+ts=ml/(%pi*(d/2)**2)    //ultimate tensile strength(MPa)
+de=(L-l)*100/l          //ductility % of elongation(%)
+dr=((2*D)**2-D**2)*100/(2*D)**2      //ductility % of reduction(%)
+t=(2/3)*ts*de/100        //modulus of toughness(Pa)
+
+//Result
+printf("\n ultimate tensile strength is %0.0f  MPa",ts/10**6)
+printf("\n modulus of toughness is %0.0f  *10**6 Pa",t/10**6)
diff --git a/3763/CH12/EX12.4/Ex12_4.sce b/3763/CH12/EX12.4/Ex12_4.sce
new file mode 100644
index 000000000..6bda7fc58
--- /dev/null
+++ b/3763/CH12/EX12.4/Ex12_4.sce
@@ -0,0 +1,26 @@
+clear
+//
+//
+//
+
+//Variable declaration
+pl1=206850*10**3     //proportional limit(Pa)
+pl2=310275*10**3     //proportional limit(Pa)
+pl3=413700*10**3     //proportional limit(Pa)
+s2=0.0615      //strain
+s3=0.2020      //strain
+Y=2.0685*10**11    //young's modulus(Pa)
+
+//Calculation
+e1=pl1/Y     //elastic strain in 1st case
+e2=pl2/Y     //elastic strain in 2nd case
+p2=s2-e2     //plastic strain in 2nd case
+r2=e2*100/p2   //ratio of elastic and plastic strain in 2nd case
+e3=pl3/Y     //elastic strain in 2nd case   
+p3=s3-e3     //plastic strain in 2nd case 
+r3=e3*100/p3   //ratio of elastic and plastic strain in 3rd case
+
+//Result
+printf("\n elastic strain in 1st case is %0.3f",e1)
+printf("\n ratio of elastic and plastic strain in 2nd case is %0.3f",r2)
+printf("\n ratio of elastic and plastic strain in 3rd case is %0.3f ",r3)
diff --git a/3763/CH12/EX12.5/Ex12_5.sce b/3763/CH12/EX12.5/Ex12_5.sce
new file mode 100644
index 000000000..dc9ea51e1
--- /dev/null
+++ b/3763/CH12/EX12.5/Ex12_5.sce
@@ -0,0 +1,22 @@
+clear
+//
+//
+//
+
+//Variable declaration
+s=12411*10**3      //stress(Pa)
+t=0.0168     //tension
+e=0.127     //elongation(cm)
+l=15.24     //length(cm)
+g=9.8
+L=68.04     //load(kg)
+
+//Calculation
+E_eff=s/t     //effective modulus(Pa)
+S=e/l      
+W=E_eff*S
+A=L*g/W       //cross sectional area(m**2)
+
+//Result
+printf("\n effective modulus is %0.3f  *10**3 Pa",E_eff/10**3)
+printf("\n cross sectional area is %0.4f  *10**-4 m**2",A*10**4)
diff --git a/3763/CH12/EX12.6/Ex12_6.sce b/3763/CH12/EX12.6/Ex12_6.sce
new file mode 100644
index 000000000..6826d3552
--- /dev/null
+++ b/3763/CH12/EX12.6/Ex12_6.sce
@@ -0,0 +1,19 @@
+clear
+//
+//
+//
+
+//Variable declaration
+E=35*10**10     //youngs modulus(Pa)
+gama=2     //specific surface energy(J/m**2)
+C=2*10**-6    //length(m)
+x=17700    
+y=2.1
+z=31.25
+
+//Calculation
+sigma_f=sqrt(2*E*gama/(%pi*C))     //fracture stress(Pa)
+T=x/((sigma_f/(9.8*10**6))-y+z)    //transition temperature(K)
+
+//Result
+printf("\n transition temperature is %0.0f  K",T)
diff --git a/3763/CH12/EX12.7/Ex12_7.sce b/3763/CH12/EX12.7/Ex12_7.sce
new file mode 100644
index 000000000..1e7623821
--- /dev/null
+++ b/3763/CH12/EX12.7/Ex12_7.sce
@@ -0,0 +1,26 @@
+clear
+//
+//
+//
+
+//Variable declaration
+h1=1
+h2=1
+k1=1
+k2=1
+l1=1
+l2=1
+l3=0
+s=3.5*10**6      //stress(Pa)
+
+//Calculation
+x=sqrt(h1**2+k1**2+l1**2)
+y=sqrt(h2**2+k2**2+l2**2)
+z=sqrt(h2**2+k2**2+l3**2)
+cos_phi=((h1*h2)-(k1*k2)+(l1*l2))/(x*y)
+sin_phi=sqrt(1-(cos_phi)**2)
+cos_theta=((h1*h2)+(k1*k2)+(l1*l3))/(x*z)
+ss=s*cos_theta*cos_phi*sin_phi        //critical resolved shear stress(Pa)
+
+//Result
+printf("\n critical resolved shear stress is %0.3f  MPa",ss/10**6)
diff --git a/3763/CH12/EX12.8/Ex12_8.sce b/3763/CH12/EX12.8/Ex12_8.sce
new file mode 100644
index 000000000..68515b1cd
--- /dev/null
+++ b/3763/CH12/EX12.8/Ex12_8.sce
@@ -0,0 +1,40 @@
+clear
+//
+//
+//
+
+//Variable declaration
+dz1=4*10**-18    //diffusivity(m**2/s)
+dz2=5*10**-13    //diffusivity(m**2/s)
+T1=773    //temperature(K)
+T2=1273   //temperature(K)
+T3=573   //temperature(K)
+T4=973   //temperature(K)
+
+//Calculation
+x1=(log(dz1))
+
+y1=(log(dz2))
+
+x2=(-1/(8.314*T1))
+
+y2=(-1/(8.314*T2))
+
+x=((x1-y1))
+
+y=((x2-y2))
+
+Q=x/y      //activation energy(J/mol)
+z=(y1-(y2*Q))
+
+D0=exp(z)          //diffusion coefficient(m**2/Vs)
+D1=D0*exp(-Q/(8.314*T3))    //diffusivity at 300 C(m**2/s)
+D2=D0*exp(-Q/(8.314*T4))    //diffusivity at 700 C(m**2/s)
+
+//Result
+printf("\n activation energy is %0.3f  kJ/mol",Q/10**3)
+printf("\n answer varies due to rounding off errors")
+printf("\n diffusion coefficient is %0.3f  *10**-4 m**2/s",D0*10**4)
+printf("\n diffusivity at 300 C is %0.2f  *10**-23 m**2/s",D1*10**23)
+printf("\n diffusivity at 700 C is %0.3f  *10**-15 m**2/s",D2*10**15)
+printf("\n answer given in the book is wrong")
diff --git a/3763/CH12/EX12.9/Ex12_9.sce b/3763/CH12/EX12.9/Ex12_9.sce
new file mode 100644
index 000000000..0b93d8550
--- /dev/null
+++ b/3763/CH12/EX12.9/Ex12_9.sce
@@ -0,0 +1,16 @@
+clear
+//
+//
+//
+
+//Variable declaration
+D0=0.73*10**-4     //diffusion coefficient(m**2/s)
+Q=170*10**3     //activation energy(J/mol)
+R=8.314 
+T=873    //temperature(K)
+
+//Calculation
+D=D0*exp(-Q/(R*T))     //diffusion(m**2/s)
+
+//Result
+printf("\n diffusion is %0.1f  *10**-15 m**2/s",D*10**15)