1 files changed, 30 insertions, 0 deletions
diff --git a/278/CH5/EX5.14/ex_5_14.sce b/278/CH5/EX5.14/ex_5_14.sce
new file mode 100755
index 000000000..4547aa4db
--- /dev/null
+++ b/278/CH5/EX5.14/ex_5_14.sce
@@ -0,0 +1,30 @@
+//cal stress at A and B
+clc
+//solution
+//given
+W=3000//N
+T=10^6//N-mm
+P=15000//N
+d=50//mm
+x=250//mm
+pi=3.14
+A=(pi/4)*d^2//mm^2//area of shaft
+f1=P/A//tensile stress at A and B
+M=W*x//N-mm
+Z=(pi/32)*d^3//mm^3
+f2=M/Z//N/mm^2
+fa=f1+f2//N/mm^2
+fb=f2-f1//N/mm^2//tensile stress at B
+fs=16*T/(pi*d^3)//N/mm^3
+Fama=(fa/2)+0.5*sqrt((fa)^2+4*(fs)^2)//max stress at A
+Fami=(fa/2)-0.5*sqrt((fa)^2+4*(fs)^2)//min stress at A
+Tama=0.5*[sqrt(((fa)^2)+(4*(fs)^2))]// max shear stress at A
+Fbma=(fb/2)+0.5*sqrt((fb)^2+4*(fs)^2)//max stress at B
+Fbmi=(fb/2)-0.5*sqrt((fb)^2+4*(fs)^2)//min stress at B
+Tbma=0.5*[sqrt(((fb)^2)+(4*(fs)^2))]//max shear stress at B
+printf("the max stress at A is,% f N/mm^2\n",Fama)
+printf("the min stress at A is,% f N/mm^2\n",Fami)
+printf("the max stress at B is,% f N/mm^2\n",Fbma)
+printf("the max stress at B is,% f N/mm^2\n",Fbmi)
+printf("the max shear stress at A is,%f N/mm^2\n",Tama)
+printf("the max shear stress at B is,%f N/mm^2",Tbma)