1 files changed, 19 insertions, 0 deletions
diff --git a/3683/CH3/EX3.13/Ex3_13.sce b/3683/CH3/EX3.13/Ex3_13.sce
new file mode 100644
index 000000000..cb4a82b0b
--- /dev/null
+++ b/3683/CH3/EX3.13/Ex3_13.sce
@@ -0,0 +1,19 @@
+Bf=500//width of flange, in mm
+Bw=250//breadth of web, in mm
+Df=100//thickness of flange, in mm
+d=500//effective depth, in mm
+sigma_cbc=5//in MPa
+sigma_st=140//in MPa
+m=18.66//modular ratio
+Ast=2000//in sq mm
+//to find critical depth of neutral axis
+Xc=d/(1+sigma_st/(m*sigma_cbc))//in mm
+//assume x>Df
+x=(m*Ast*d+Bf*Df^2/2)/(m*Ast+Bf*Df)//in mm
+//as x>Xc, beam is over-reinforced
+sigma_cbc_dash=sigma_cbc*(x-Df)/x//in MPa
+//to find lever arm
+z=d-(sigma_cbc+2*sigma_cbc_dash)/(sigma_cbc+sigma_cbc_dash)*Df/3//in mm
+//taking moments about tensile steel
+Mr=Bf*Df*(sigma_cbc+sigma_cbc_dash)*z/2//in N-mm
+mprintf("Moment of resistance of the beam=%f kN-m", Mr/10^6)