Scribd
Upload
22 views3 pages

Exp 10

This MATLAB program designs a singly reinforced concrete beam using the limit state method. It calculates the effective depth, bending moment, shear force, tension reinforcement area and number of bars needed, and stirrup spacing. Input values like beam dimensions, grades of concrete and steel, and loads are given. The program outputs the effective depth, reinforcement details, and checks that design requirements are met.

Uploaded by

Shravan kumarssk
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
22 views3 pages

Exp 10

This MATLAB program designs a singly reinforced concrete beam using the limit state method. It calculates the effective depth, bending moment, shear force, tension reinforcement area and number of bars needed, and stirrup spacing. Input values like beam dimensions, grades of concrete and steel, and loads are given. The program outputs the effective depth, reinforcement details, and checks that design requirements are met.

Uploaded by

Shravan kumarssk
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 3

NAME:M.

Srinith
21311A0123
MATLAB EXP 10
% Exp 10- Program for design of a Singly Reinforced concrete
beam by % limit State Method
clear all
InpFile=fopen('D:\matlabdata\r123inp10.txt','r');
OutFile=fopen('D:\matlabdata\r123exp10OUT.txt','w');
size1=[1];
L=fscanf(InpFile,'%f',size1); %Clear span of beam (m)
b=fscanf(InpFile,'%f',size1); %Width of beam (mm)
fck=fscanf(InpFile,'%f',size1); %Grade of concrete (20 for
M20)
fy=fscanf(InpFile,'%f',size1); %Grade of main steel bar
fys=fscanf(InpFile,'%f',size1); %Grade of steel for stirrups
IL=fscanf(InpFile,'%f',size1); %Imposed load (kN/m)
Ww=fscanf(InpFile,'%f',size1); % Width of the wall (mm)
dm=fscanf(InpFile,'%f',size1); % Diameter of main bar (mm)
ds=fscanf(InpFile,'%f',size1); % Diameter of stirrups in
(mm)
%calculation of effective depth
d=L*1000/20;
%calculation of effective length
EL1=L+(d/1000);
EL2=L+(Ww/1000);
if EL1<EL2
EL=EL1;
else
EL=EL2;
end
fprintf(OutFile,'Effective length of beam is %.1f\r\n',EL);
%calculation of bending moment and shear force
Dok=0;
while (Dok == 0)
DL=(b/10000)*(d/1000)*24; %Dead Load in kN/m
Wu=(1.5*DL)+(1.5*IL); %Design Load in kN/m
Mu=(Wu*EL^2)/8; %Max Bending Moment (kN-m)
Vu=(Wu*EL)/2; %Max Shear force (kN)
fprintf(OutFile,'The design bending moment(kNm) is %.1f\n
and design shear force (kN) is %.1f\r\n',Mu,Vu);
% check for Depth requirement
Dmin=(Mu*10^6/(.138*fck*b))^0.5
if Dmin>d %revise d
d=ceil(1.1*Dmin); Dok=0
else Dok=1
end
NAME:M.Srinith
21311A0123
MATLAB EXP 10
end
fprintf(OutFile,'Effective depth of beam is: %.1f\r\n',d);
fprintf(OutFile,'Overall depth of beam is: %.1f\r\n',d+20);
%Calculation of area of steel in tension
syms Ast;
Eqn=0.87*fy*Ast*(d-(0.42*0.87*fy*Ast)/
(0.36*b*fck))==Mu*10^6;
Ast=solve(Eqn,Ast);
At=min(double(Ast)) %considering the minimum value of two
solutions for Ast
if (isreal(At)==0)
disp('NO REAL VALUE solution. REVISE the data'\n);
fprintf(OutFile, 'NO real value solution. Revise the data'\
n);
end
bars=At/((pi*dm^2)/4);
No_bars=ceil(bars);
fprintf(OutFile, 'Tension reinforcement is: %.1f\r\n and No.
of bars (%d dia) are %d\n',At,dm,No_bars);
%Calculation of shear reinforcement
tv=Vu*1000/(b*d);
%Astp=fscanf(InpFile,'%f',size1); %Area of steel provided in
tension:
Psteel=At*100/(b*d);
tc= 0.28+Psteel*(0.82-0.28)/(3.0-0.15); %Approximate fitting
table from IS456
tcmax=fscanf(InpFile,'%f',size1); %tcmax from T20 of IS456
if tv<=tcmax
if tv>tc
Vc=tc*b*d;
Vus=Vu*1000-Vc; %Balance shear load that needs reinforcement
Sv1=0.87*fys*(pi/4)*(ds^2)*d/Vus;
S=[Sv1,300,0.75*d];
Sv=min(S);
fprintf(OutFile,'c/c Spacing of stirrups is: %.1f\r\n',Sv);
else
fprintf(OutFile,'Concrete can take care of shear force,
provide minimum shear reinforcement');
end
else
fprintf(OutFile,'Redesign the section by changing the
dimensions');
end
NAME:M.Srinith
21311A0123
MATLAB EXP 10
fclose(InpFile);
fclose(OutFile);

INPUT:

6.5

250

20

415

250

10

250

13

2.8

OUTPUT:

Effective length of beam is 6.8

The design bending moment(kNm) is 87.1

and design shear force (kN) is 51.6

The design bending moment(kNm) is 87.4

and design shear force (kN) is 51.8

Effective depth of beam is: 391.0

Overall depth of beam is: 411.0

Tension reinforcement is: 736.1

and No. of bars (13 dia) are 6

c/c Spacing of stirrups is: 293.3

You might also like