7 files changed, 111 insertions, 0 deletions
diff --git a/866/CH9/EX9.1/9_1.sce b/866/CH9/EX9.1/9_1.sce
new file mode 100755
index 000000000..ca02d8d2c
--- /dev/null
+++ b/866/CH9/EX9.1/9_1.sce
@@ -0,0 +1,15 @@
+clc
+//initialisation of variables
+w= 10 //KN/m
+L= 6 //m
+sigmaxallowable= 155 //N/mm^2
+Modulusofuniversalbeam= 307600 //mm^3
+Satisfactorybeam= 254*102*28 
+//CALCULATIONS
+Mmaxnormal= (w*L^2)/8
+Zemin= (Mmaxnormal)/sigmaxallowable
+Totalload= w+((28*9.81)/10^3)
+Mmaxload= (Totalload*L^2)/8
+Allowablestress= (Mmaxload*10^3*10^3)/Modulusofuniversalbeam
+//RESULTS
+printf ('Satisafactorybeam=% 2f (254*102*28)',Satisfactorybeam)
diff --git a/866/CH9/EX9.2/9_2.sce b/866/CH9/EX9.2/9_2.sce
new file mode 100755
index 000000000..555bb2e7e
--- /dev/null
+++ b/866/CH9/EX9.2/9_2.sce
@@ -0,0 +1,16 @@
+clc
+//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) 
+
diff --git a/866/CH9/EX9.3/9_3.sce b/866/CH9/EX9.3/9_3.sce
new file mode 100755
index 000000000..05cf97ddc
--- /dev/null
+++ b/866/CH9/EX9.3/9_3.sce
@@ -0,0 +1,13 @@
+clc
+//initialisation of variables
+M= 100*10^6 //Nmm
+BeamB= 300 //mm
+BeamH= 200 //mm
+BeamT= 25 //mm
+BeamT2= 20 //mm
+//CALCULATIONS
+Iy= (2*BeamT2*BeamH^3/12)+((BeamB-2*BeamT2)*BeamT^3/12)
+alphabyz=-M/Iy
+//RESULTS
+printf ('sigmax= %.2f *z',alphabyz)
+
diff --git a/866/CH9/EX9.4/9_4.sce b/866/CH9/EX9.4/9_4.sce
new file mode 100755
index 000000000..c3e7e2746
--- /dev/null
+++ b/866/CH9/EX9.4/9_4.sce
@@ -0,0 +1,19 @@
+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)
diff --git a/866/CH9/EX9.5/9_5.sce b/866/CH9/EX9.5/9_5.sce
new file mode 100755
index 000000000..90d7516b9
--- /dev/null
+++ b/866/CH9/EX9.5/9_5.sce
@@ -0,0 +1,18 @@
+clc
+//initialisation of variables
+BeamB= 200 //mm
+BeamH= 200 //mm
+BeamT= 20 //mm
+BeamT2= 20 //mm
+sigmaxtop= 150 //N/mm^2
+sigmaxbottom= -150 //N/mm^2
+//CALCULATIONS
+y= ((BeamB*BeamT*BeamT/2)+(BeamH*BeamT2*((BeamH/2)+BeamT)))/(BeamB*BeamT+BeamH*BeamT2)
+Iz= ((BeamB*(y^3))/3)-(((BeamB-BeamT2)*(y-BeamT)^3)/3)+(BeamT2*((BeamH+BeamT-y)^3)/3)
+P1= sigmaxtop/((1/(BeamB*BeamT+BeamH*BeamT2))+((y-BeamT/2)*y/Iz))  
+P2= (sigmaxbottom/((1/(BeamB*BeamT+BeamH*BeamT2))+((y-BeamT/2)*(-(BeamH+BeamT-y))/Iz)))
+if(P1>P2)
+    disp(P2)
+else
+    disp(P1) 
+end
diff --git a/866/CH9/EX9.6/9_6.sce b/866/CH9/EX9.6/9_6.sce
new file mode 100755
index 000000000..6343e108b
--- /dev/null
+++ b/866/CH9/EX9.6/9_6.sce
@@ -0,0 +1,10 @@
+clc
+//initialisation of variables
+Height= 7 //m
+Thickness= 0.6 //m
+density= 2000 //Kg/m^3
+//CALCULATIONS
+W= density*9.81*Thickness*Height
+p= 0.1*W*2/(Height*Height)
+//RESULTS
+printf ('p= %.2f N/m^1',p)
diff --git a/866/CH9/EX9.8/9_8.sce b/866/CH9/EX9.8/9_8.sce
new file mode 100755
index 000000000..8b8cb0f3f
--- /dev/null
+++ b/866/CH9/EX9.8/9_8.sce
@@ -0,0 +1,20 @@
+clc
+//initialisation of variables
+M= 1500 //Nm
+x1= 40 //mm
+x2= 80 //mm
+x3= 8 //mm
+y1= 8 //mm
+y2= 80 //mm
+//CALCULATIONS
+Ycentroid= (((x1+x2)*y1*y1/2)+(y2*x3*((y2/2)+y1)))/((x1+x2)*y1+y2*x3)
+Zcentroid= (((y2*x3*y1/2)+((x1+x2)*x3*((y2/2)+y1)/2))/((x1+x2)*y1+y2*x3))
+Iz= (((x1+x2)*(y1)^3)/12)+(x1+x2)*y1*((Ycentroid-(y1/2))^2)+(x3/12)*(y2^3)+y2*x3*((y2/2)-(Ycentroid-y1))^2
+Iy= (y1*(((x1+x2)^3)/12))+(x1+x2)*y1*(y1^2)+((x3^3)/12)*(y2)+y2*x3*((Zcentroid-(x3/2))^2)
+Izy= (x1+x2)*y1*(-y1)*(Ycentroid-(y1/2))+y2*x3*(Zcentroid-(x3/2))*(-(y2/2)+(Ycentroid-y1))
+sigmax= ((-M*10^3*Izy*x3)+(-M*10^3*Iy*(y2-Ycentroid+y1)))/((Iz*Iy)-(Izy)^2)
+//RESULTS
+printf ('Maximum direct stress= %.2f N/mm^2',sigmax)
+
+
+