initial commit / add all books

6 files changed, 119 insertions, 0 deletions

diff --git a/3456/CH12/EX12.1/Ex12_1.sce b/3456/CH12/EX12.1/Ex12_1.sce
new file mode 100644
index 000000000..6ef4d4644
--- /dev/null
+++ b/3456/CH12/EX12.1/Ex12_1.sce
@@ -0,0 +1,17 @@
+//Example 12.1

+//Mean Stress

+//Page No. 387

+clc;clear;close;

+

+sigma_u=158;             // in ksi

+sigma0=147;             // in ksi

+sigma_e=75;             // in ksi

+l_max=75;             // in ksi

+l_min=-25;             // in ksi

+sf=2.5;                //no unit

+sigma_m=(l_max+l_min)/2;

+sigma_a=(l_max-l_min)/2;

+sigma_e=sigma_e/sf;

+A=sigma_a/sigma_e+sigma_m/sigma_u;

+D=sqrt(4*A/%pi);

+printf('\nBar Diameter = %g in',D);

diff --git a/3456/CH12/EX12.2/Ex12_2.sce b/3456/CH12/EX12.2/Ex12_2.sce
new file mode 100644
index 000000000..8566ca212
--- /dev/null
+++ b/3456/CH12/EX12.2/Ex12_2.sce
@@ -0,0 +1,14 @@
+//Example 12.2

+//Low Cycle Fatigue

+//Page No. 391

+clc;clear;close;

+

+sigma_b=75;               //in MPa

+e_b=0.000645;              //no unit

+e_f=0.3;                   //no unit

+E=22*10^4;                 //in MPa

+c=-0.6;                    //no unit

+d_e_e=2*sigma_b/E;

+d_e_p=2*e_b-d_e_e;

+N=(d_e_p/(2*e_f))^(1/c)/2;

+printf('\nd_e_e = %g\nd_e_p = %g\nNumber of Cycles = %g cycles',d_e_e,d_e_p,N);

diff --git a/3456/CH12/EX12.3/Ex12_3.sce b/3456/CH12/EX12.3/Ex12_3.sce
new file mode 100644
index 000000000..650365760
--- /dev/null
+++ b/3456/CH12/EX12.3/Ex12_3.sce
@@ -0,0 +1,15 @@
+//Example 12.3

+//Fatigue Crack Proportion

+//Page No. 401

+clc;clear;close;

+

+ai=0.5;                  //in mm

+ai=ai*10^-3;           //conversion to m

+sigma_max=180;             //in MPa

+Kc=100;                  //MPam^(1/2)

+alpha=1.12;                //no unit

+p=3;                       //no unit

+A=6.9*10^-12;              //in MPam^(1/2)

+af=(Kc/(sigma_max*alpha))^2/%pi;

+Nf=(af^(1-(p/2))-ai^(1-(p/2)))/((1-p/2)*A*sigma_max^3*%pi^(p/2)*alpha^p);

+printf('Fatigue Cycles = %g cycles',Nf);

diff --git a/3456/CH12/EX12.4/Ex12_4.sce b/3456/CH12/EX12.4/Ex12_4.sce
new file mode 100644
index 000000000..f2d2e5a70
--- /dev/null
+++ b/3456/CH12/EX12.4/Ex12_4.sce
@@ -0,0 +1,25 @@
+//Example 12.4

+//Stress Concentration of Fatigue

+//Page No. 404

+clc;clear;close;

+

+rho=0.0004;            //no unit  

+S_u=190;            //in ksi

+S_u=S_u*1000;           //conversion to psi

+M=200;               //in inches-lb

+Pm=5000;            //in lb

+D=0.5;                //in inches

+dh=0.05;               //in inches

+r=dh/2;

+Kt=2.2;               //no unit

+Kf=1+(Kt-1)/(1+sqrt(rho/r));

+q=(Kf-1)/(Kt-1);

+A=%pi/4*D^2;

+sigma_m=Pm/A;

+I=%pi/64*D^4;

+sigma_a=Kf*((M*D)/(2*I));

+sigma_max=sigma_a+sigma_m;

+sigma_min=sigma_a-sigma_m;

+sigma_e=S_u/2;

+sigma_a1=sigma_e/Kf*(1-sigma_m/S_u);

+printf('\nMean Stress = %g psi\nFluctuating Bending Stress = %g psi\nEffective Maximum Stress = %g psi\nEffective Minimum Stress = %g psi\nsigma_a = %g psi\n\n\nNote: Calculation Errors in the book',sigma_m,sigma_a,sigma_max,sigma_min,sigma_a1);

diff --git a/3456/CH12/EX12.5/Ex12_5.sce b/3456/CH12/EX12.5/Ex12_5.sce
new file mode 100644
index 000000000..e906edbb8
--- /dev/null
+++ b/3456/CH12/EX12.5/Ex12_5.sce
@@ -0,0 +1,15 @@
+//Example 12.5

+//Infinite Life Design

+//Page No. 422

+clc;clear;close;

+

+Kt=1.68;               //no unit

+q=0.9;                 //no unit

+sigma_ed=42000;         //in psi

+Cs=0.9;                 //no unit

+Cf=0.75;               //no unit

+Cz=0.81;                //no unit

+Kf=q*(Kt-1)+1;

+sigma_e=sigma_ed*Cs*Cf*Cz;

+sigma_en=sigma_e/Kf;

+printf('\nFatigue Limit = %g psi',sigma_en);

diff --git a/3456/CH12/EX12.6/Ex12_6.sce b/3456/CH12/EX12.6/Ex12_6.sce
new file mode 100644
index 000000000..153defcf9
--- /dev/null
+++ b/3456/CH12/EX12.6/Ex12_6.sce
@@ -0,0 +1,33 @@
+//Example 12.6

+//Local Strain method

+//Page No. 424

+clc;clear;close;

+

+funcprot(0);

+K=189;              //in ksi

+n=0.12;             //no unit

+ef=1.06;             //no unit

+sigma_f=190;          //in ksi

+b=-0.08;              //no unit

+c=-0.66;           //no unit

+E=30*10^6;           //in psi

+E=E/1000;        //conversion to ksi\

+s=200;            //in ksi

+sigma_m=167;         //in ksi

+sigma_a=17;            //in ksi

+se=s^2/E;

+deff('y=f(ds)','y=(ds^2)/(2*E)+(ds^((1+n)/n))/(2*K)^(1/n)-se/2');

+[ds,v,info]=fsolve(0,f);

+de=se/ds;

+deff('y=f1(N2)','y=N2^-b*(sigma_f/E)+ef*N2^-c-de/2');

+[N2,v,info]=fsolve(0,f1);

+N2=1/N2;

+N_1=N2/2;

+de_e2=sigma_a/E;

+deff('y=f2(N2)','y=N2^-b*((sigma_f-sigma_m)/E)+ef*N2^-c-de_e2');

+[N2,v,info]=fsolve(0,f2);

+N2=1/N2;

+N_2=N2/2;

+C_pd=2*60*60*8;

+f=N_2/C_pd;

+printf('\nNumber of cycles = %g cycles\nFatigue damage per cycle = %g\nNumber of cycles with correction of mean stress= %g cycles\nFatigue damage per cycle  with correction of mean stress= %g damage per year\nShaft will fail in %g days',N_1,1/N_1,N_2,1/N_2,f);