1 files changed, 16 insertions, 19 deletions
diff --git a/866/CH9/EX9.4/9_4.sce b/866/CH9/EX9.4/9_4.sce
index c3e7e2746..7899cb205 100755
--- a/866/CH9/EX9.4/9_4.sce
+++ b/866/CH9/EX9.4/9_4.sce
@@ -1,19 +1,16 @@
-clc
-//initialisation of variables
-BeamB= 300 //mm
-BeamL= 200 //mm
-BeamT= 25 //mm
-BeamT2= 20 //mm
-alpha= 30 //degrees
-Mz= 100*10^6*cosd(alpha)//Nmm
-My= 100*10^6*sind(alpha)//Nmm
-//CALCULATIONS
-Iz= ((BeamL*BeamB^3)/12)-((BeamL-BeamT)*(BeamB-2*BeamT2)^3)/12
-Iy= (2*BeamT2*BeamL^3/12)+((BeamB-2*BeamT2)*BeamT^3/12)
-sigmaxtl= (+Mz/Iz)*(BeamB/2)+(My/Iy)*(BeamH/2)
-sigmaxtr= (+Mz/Iz)*(BeamB/2)+(My/Iy)*(-BeamH/2)
-inclination= atand((My*Iz)/(Mz*Iy))
-//RESULTS   
-printf ('Stress at top left of the beam= %.2fN/mm^2(Tension)',sigmaxtl) 
-printf ('\n Stress at top right of the beam= %.2fN/mm^2(compression)',sigmaxtr) 
-printf ('\n Incination= %.2fdegrees',inclination)
+clc
+clear
+//initialisation of variables
+M= 100*10^6 //Nmm
+BeamB= 300 //mm
+BeamL= 200 //mm
+BeamT= 25 //mm
+BeamT2= 20 //mm
+//CALCULATIONS
+Iz= ((BeamL*BeamB^3)/12)-((BeamL-BeamT)*(BeamB-2*BeamT2)^3)/12
+sigmaxbyY= -M/Iz
+SB= sigmaxbyY*(BeamB/2)
+ST= sigmaxbyY*(-BeamB/2)
+//RESULTS   
+printf ('Stress at top of the beam= %.2fN/mm^2(Tension)',ST) 
+printf ('\n Stress at bottom of the beam= %.2fN/mm^2(compression)',SB)
+\ No newline at end of file