summaryrefslogtreecommitdiff
path: root/3456/CH12/EX12.6/Ex12_6.sce
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /3456/CH12/EX12.6/Ex12_6.sce
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '3456/CH12/EX12.6/Ex12_6.sce')
-rw-r--r--3456/CH12/EX12.6/Ex12_6.sce33
1 files changed, 33 insertions, 0 deletions
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);