1 files changed, 74 insertions, 0 deletions
diff --git a/698/CH4/EX4.2/P2_offset_bar.sce b/698/CH4/EX4.2/P2_offset_bar.sce
new file mode 100644
index 000000000..220bd7f0c
--- /dev/null
+++ b/698/CH4/EX4.2/P2_offset_bar.sce
@@ -0,0 +1,74 @@
+clc
+//Example 4.2
+//Offset bar
+
+//------------------------------------------------------------------------------
+
+//Given data
+//Couple causing force
+F=900 // N
+//Distance between forces
+d=0.14 //m
+//Dimensions of bar
+h=0.025 // m
+R=0.09 // m
+t=0.05 // m
+
+res2=mopen(TMPDIR+'2_offset_bar.txt','wt')
+
+//Bending moment
+M=d*F
+mfprintf(res2,'(a) The bending moment at every section is=%d Nm\n\n',M)
+
+mfprintf(res2,'(b) Refer to figure 4.3[A. Hall; Schaum''s Outline on Machine Design]\n')
+mfprintf(res2,'Tension occurs in upper fibres and compression in lower fibres.\n')
+mfprintf(res2,'Maximum normal stress occurs at R=90mm and is compressive\n')
+mfprintf(res2,'Location of maximum tensile stress is to be determined\n\n')
+
+//Radii of curvature of fibres
+ri=R- h/2
+ro=R+ h/2
+rn=h/ log(ro/ri)
+
+//Distance between centre of gravity and neutral axis
+e=R-rn
+mfprintf(res2,'(c)For section R=90mm\n\tDistance between centre of gravity and neutral axis e=%0.6f m\n\n',e)
+
+//Checking for maximum tensile stress at R=90 mm
+Si=(M*0.01192)/(h*t*e*ri)
+mfprintf(res2,'(d)For section with R=90mm, stess at inner fibre\n\t')
+mfprintf(res2,'Si=M*hi/ A*e*ri\n')
+mfprintf(res2,'\t=%0.2f MN/m^2 compression\n',Si* 10^-6)
+s=Si/2
+So=(M*0.0131)/(h*t*e*ro)
+mfprintf(res2,'Stess at outer fibre\n\t')
+mfprintf(res2,'So=M*ho/ A*e*ro\n')
+mfprintf(res2,'\t=%0.2f MN/m^2 tension\n\n',So* 10^-6)
+
+//Checking for maximum tensile stress at R=115 mm
+R=0.115
+ri=R- h/2
+ro=R+ h/2
+rn=h/ log(ro/ri)
+e=R-rn
+e=ceil(e* 10^4)
+e=e* 10^-4
+mfprintf(res2,'For section R=115mm\n\tDistance between centre of gravity and neutral axis e=%0.4f m\n\n',e)
+Si=(M*0.012)/(h*t*e*ri)
+mfprintf(res2,'(e)For section with R=115mm, stess at inner fibre\n\t')
+mfprintf(res2,'Si=M*hi/ A*e*ri\n')
+mfprintf(res2,'\t=%0.2f MN/m^2 tension\n',Si* 10^-6)
+So=(M*0.013)/(h*t*e*ro)
+mfprintf(res2,'Stess at outer fibre\n\t')
+mfprintf(res2,'So=M*ho/ A*e*ro\n')
+mfprintf(res2,'\t=%0.2f MN/m^2 compression\n\n',So* 10^-6)
+
+//Maximum stresses
+mfprintf(res2,'(f) Maximum tension occurs at inner fibre of section R=115mm\n')
+mfprintf(res2,'Maximum compression occurs at inner fibre of section R=90mm\n')
+mfprintf(res2,'Maximum shear stress=half of greatest difference between any two extremes=%0.1f MN/m^2\n',s* 10^-6)
+
+mclose(res2)
+editor(TMPDIR+'2_offset_bar.txt')
+//------------------------------------------------------------------------------
+//-----------------------------End of program-----------------------------------
+\ No newline at end of file