summaryrefslogtreecommitdiff
path: root/3683/CH3/EX3.11/Ex3_11.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3683/CH3/EX3.11/Ex3_11.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3683/CH3/EX3.11/Ex3_11.sce')
-rw-r--r--3683/CH3/EX3.11/Ex3_11.sce27
1 files changed, 27 insertions, 0 deletions
diff --git a/3683/CH3/EX3.11/Ex3_11.sce b/3683/CH3/EX3.11/Ex3_11.sce
new file mode 100644
index 000000000..34b03a4d7
--- /dev/null
+++ b/3683/CH3/EX3.11/Ex3_11.sce
@@ -0,0 +1,27 @@
+Bf=1500//width of flange, in mm
+Df=150//thickness of flange, in mm
+d=600//effective depth, in mm
+sigma_cbc=7//in MPa
+sigma_st=230//in MPa
+m=13.33//modular ratio
+Ast=1964//in sq mm
+Asc=1140//in sq mm
+top_cover=50//in mm
+//to find critical depth of neutral axis
+Xc=d/(1+sigma_st/(m*sigma_cbc))//in mm
+//assume x>Df; equating moments of area on compression and tension sides about N.A.
+x=(m*Ast*d+Bf*Df^2/2+(1.5*m-1)*Asc*top_cover)/(m*Ast+Bf*Df+(1.5*m-1)*Asc)// in mm
+//we find that x<Df, hence our assumption that x>Df is wrong
+//to find x using Bf(x^2)/2 + (1.5m-1)Asc(x-d')=mAst(d-x), which becomes of the form px^2+qx+r=0
+p=Bf/2
+q=m*Ast+(1.5*m-1)*Asc
+r=-(m*Ast*d+(1.5*m-1)*Asc*top_cover)
+//solving quadratic equation
+x=(-q+sqrt(q^2-4*p*r))/(2*p)//in mm
+//as x<Xc, beam is under-reinforced
+sigma_cbc=sigma_st/m*x/(d-x)//in MPa
+sigma_cbc_dash=sigma_cbc*(x-top_cover)/x//stress in concrete at level of compression steel, in MPa
+//taking moments about tensile steel
+Mr=Bf*x*sigma_cbc*(d-x/3)/2+(1.5*m-1)*Asc*sigma_cbc_dash*(d-top_cover)//in N-mm
+mprintf("Moment of resistance of the beam=%f kN-m", Mr/10^6)
+//answer given in textbook is incorrect