summaryrefslogtreecommitdiff
path: root/3761/CH4/EX4.4/Ex4_4.sce
blob: 025fbc19554d94b4f27141c385724c7cb4a6f119 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
L=6
Df=100
bw=250
b=1000
D=600
d=520
M=200
disp("Example 4.4")
disp("Service Load Moment = 200kNm","Ast=6-28mmdia bar","Effective depth d=520 mm","Depth D=600mm","b=1000mm","bw=250mm","Depth of Flange Df=100mm","Span of Beam L=6m","Given Data:")
disp("Grade of Steel =Fe250","Grade of Concrete=M20")
disp("Verifying for the effective flange width b")
disp("Refering IS456:2000, Clause 23.1.2 c or Eq 4.30(b) from TB")
disp("For T-beams, bf = (lo/((lo/b)+4)+bw")
disp("bw=250mm","lo=6000mm")
lo=6000
bf=((lo)/((lo/b)+4)+bw)
disp("mm",bf,"bf=")

if(bf<b)
       disp(bf,"Since, bf<b,b=")

sigmacbc=7
m=280/(3*sigmacbc)
disp(m,"m=")
dia=28
Ast= %pi*dia*dia*6/4
mAst=m*Ast
disp("Assuming kd</Df, and equating moments of compression and transformed tension areas about the neutral axis,")
disp("bf*(kd)^2/2 = m*Ast*(d-kd)")
b1=(2*mAst/bf)
c=(-mAst*2*d/bf)
a=1
kd1=((-b1+sqrt((b1*b1)-(4*a*c)))/(2*a))
kd2=((-b1-sqrt((b1*b1)-(4*a*c)))/(2*a))
disp("mm",kd1,"therefore, kd=")
disp("kd1>Df, the assumption kd</Df is incorrect")
disp(" For kd>Df, neutral axis located in the web")
disp("Using Equation 4.31 of TB")
disp("(bf-bw)*Df*(kd-Df/2)+ bw*(kd)^2/2 = mAst*(d-kd)")
a=bw/2
B=(bf*Df-bw*Df+mAst)
c=(bw*Df*Df/2 - bf*Df*Df/2-mAst*d)

Dis=(B*B)-(4*a*c)
kd1=((-B+sqrt(Dis))/(2*a))
disp("mm",kd1,"kd=")
disp("Relating the compressive stress fc1 at the flange bottom to fc,")
disp("fc1=fc*((kd-Df)/kd)")
fact=((kd1-Df)/kd1)
disp("*fc",fact,"fc1=")
disp("Compressive Force C= 0.5*fc*bf*(kd)-0.5*fc1*(bf-bw)*(kd-Df)")
disp("Taking moments of forces about the tension steel centriod, using equation 4.34")
fc=((M*10^6)/((0.5*bf*(kd1)*(d-kd1/3)-(0.5*fact*(bf-bw)*(kd1-Df))*(d-Df-((kd1-Df)/3)))))
disp("MPa",fc,"Therefore, on solving we get, fc=")
disp("Applying C= T")
fst= (0.5*fc*bf*(kd1)-0.5*fact*fc*(bf-bw)*(kd1-Df))/Ast
disp("MPa",fst,"Therefore, fst = ")
disp("From the stress distribution diagram: fst = m*fc*((d-kd)/kd)")
fst1=m*fc*((d-kd1)/kd1)
disp("MPa",fst1, "As before")