Scribd
Upload
86 views4 pages

MatLab Code For Response of Pre-Stressed Members Subjected To Axial Loading

The document analyzes the strain variation and axial load response of four prestressed concrete members (A, B, C, D) subjected to axial loading over time. It defines variables, calculates the stress in concrete and reinforcement/strands, and determines the axial load for each member. Graphs of the load-deformation curves are plotted and compared to analyze the long-term response under varying strain conditions.

Uploaded by

Syed Gilani
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
86 views4 pages

MatLab Code For Response of Pre-Stressed Members Subjected To Axial Loading

The document analyzes the strain variation and axial load response of four prestressed concrete members (A, B, C, D) subjected to axial loading over time. It defines variables, calculates the stress in concrete and reinforcement/strands, and determines the axial load for each member. Graphs of the load-deformation curves are plotted and compared to analyze the long-term response under varying strain conditions.

Uploaded by

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

%%%____Response of Pre-Stressed Members Subjected to Axial Loading___%%%

% Strain Variation from -0.003 to +0.003


fprintf("Assignment - 2 \n Syed Kaleem Gilani \n Roll # 326957" + ...
" \n Prestressed Concrete Structures \n Long Term Response\n\n ")
syms fcp fc ecf ecp es ec Es As Ac fcr Ap fy fpy delta_ep
disp("fc = ")
fc = fcp * (2* (ecf/ecp)-(ecf/ecp)^2);
disp(fc)
fcp = -5; % Design compressive stress (Ksi)
Es = 30000; % Modulus of Elasticity of steel (Ksi)
Ep = Es; % strands modulus of elasticity (Ksi)
Ep_eff = 0.95*Ep; % Effective Modulus of elasticity (5% loss of pre stress)
ecp = -7.4*10^-3; % Maximum strain
esh = -0.0004; % Shrinkage strain in concrete
Delta_T = 50; %Temperature Drop (°F)
alpha_csp = -6*10^-6; % Co-efficient of thermal variation in concrete,steel,strands.
fcr = 0.3; % Tensile Maximum stress (Ksi)
fy = 60; % Yield stress reinforcement (Ksi)
fpy = 240; % Yield stress strands (Ksi)
%%%_________setting variable strain matrix_____________%%%
s = 7;
ecc = zeros(s,1);
J = (-0.003:0.001:0.003);
j = transpose(J);
for i = 1:s
ecc(i,1) = j(i,1);
end
e_sh = ones(s,1) .* esh;
e_cth = ones(s,1) .* (Delta_T*alpha_csp);
ecf = ecc- e_sh- e_cth
%%%___________________ MEMBER - A __________________________%%%
disp("For Member A - Non-Prestressed - 8 #6 Bars (10in x 10in)")
As = 3.53 % Area of reinforcement (in^2)
Ac = 100 - As % Net area of concrete (in^2)
%%%__________ STRESS IN CONCRETE - Member - A ____________%%%
disp("Stress in concrete fc - Member - A = ")
for i = 1:s
fc(i,1) = fcp* (2*(ecf(i,1)/ecp) - (ecf(i,1)/ecp) ^2);
if fc(i,1) > fcr
fc(i,1) = 0;
end
end
disp(double(fc))
%%%____________Stress in Reinforcements - Member - A _________%%%
disp("Stress in Reinforcements - Member - A = ")
fs = zeros(7,1);
es = ecc;
esf = es - (Delta_T*alpha_csp);
for i = 1:s
fs(i,1) = Es* esf(i,1);
if fs(i,1) > fy
fs(i,1) = fy;
elseif fs(i,1) < -fy

1
fs(i,1) = -fy;
end
end
disp(fs)
%%%____________________Axial Load, N, (Kips)- Member - A ___________%%%
NA = zeros(s,1);
for i = 1:s
NA(i,1) = fc(i,1)* Ac + fs(i,1)* As;
end
disp("Axial Load - Member A, N(Kips) = ")
double(NA);
disp(double(NA))
%%%_________________________END____________________________%%%
%%_________________________ MEMBER - B ______________________________%%
disp("For Member B - Fully-Prestressed - 6# 1/2in Bars (10in x 10in)")
Ap = 0.918 % Area of strands (in^2)
Ac = 100 - Ap % Net area of concrete (in^2)
fpbed = 180 % Prestressing (Ksi)
delta_ep = fpbed/ Ep
%%%__________ STRESS IN CONCRETE - Member - B ____________%%%
disp("Stress in concrete fc - Member - B = ")
for i = 1:s
fc(i,1) = fcp* (2*(ecf(i,1)/ecp) - (ecf(i,1)/ecp) ^2);
if fc(i,1) > fcr
fc(i,1) = 0;
end
end
disp(double(fc))
%%%____________Stress in prestressed strands - Member - B _________%%%
disp("Stress in strands - Member - B = ")
fp = zeros(7,1);
ep = ecc + delta_ep;
epf = ep- (Delta_T* alpha_csp);
for i = 1:s
fp(i,1) = Ep_eff* epf(i,1);
if fp(i,1) > fpy
fp(i,1) = fpy;
elseif fp(i,1) < -fpy
fp(i,1) = -fpy;
end
end
disp(fp)
%%%____________________Axial Load, N, (Kips) - Member - B ___________%%%
NB = zeros(s,1);
for i = 1:s
NB(i,1) = fc(i,1)* Ac + fp(i,1)* Ap;
end
disp("Axial Load - Member - B, N(Kips) = ")
double(NB);
disp(double(NB))
%%%_________________________END____________________________%%%

2
%%------------------------------MEMBER - C ----------------------------------&&
disp("For Member C - Partially-Prestressed - " + ...
"2# 1/2in and 4 #6 bars (10in x 10in)")
As = 1.24 % Area of reinforcement (in^2)
Ap = 0.306 % strands area (in^2)
Ac = 100 - As - Ap % Net area of concrete (in^2)
fpbed = 180 % prestressing (Ksi)
delta_ep = fpbed/ Ep
%%%__________ STRESS IN CONCRETE - Member - C ____________%%%
for i = 1:s
fc(i,1) = fcp* (2*(ecf(i,1)/ecp) - (ecf(i,1)/ecp) ^2);
if fc(i,1) > fcr
fc(i,1) = 0;
end
end
disp(double(fc))
%%%____________Stress in prestressed strands - Member - C _________%%%
disp("Stress in strands - Member - C = ")
fp = zeros(7,1);
ep = ecc + delta_ep;
epf = ep- (Delta_T* alpha_csp);
for i = 1:s
fp(i,1) = Ep_eff* epf(i,1);
if fp(i,1) > fpy
fp(i,1) = fpy;
elseif fp(i,1) < -fpy
fp(i,1) = -fpy;
end
end
disp(fp)
%%%____________Stress in Reinforcements - Member - C _________%%%
disp("Stress in Reinforcements - Member - C = ")
fs = zeros(7,1);
es = ecc;
esf = es - (Delta_T*alpha_csp);
for i = 1:s
fs(i,1) = Es* esf(i,1);
if fs(i,1) > fy
fs(i,1) = fy;
elseif fs(i,1) < -fy
fs(i,1) = -fy;
end
end
disp(fs)
%%%____________________Axial Load, N, (Kips)- Member C___________%%%
NC = zeros(s,1);
for i = 1:s
NC(i,1) = fc(i,1)* Ac + fp(i,1)* Ap + fs(i,1)* As;
end
disp("Axial Load - Member - C, N(Kips) = ")
double(NC);
disp(double(NC))
%%%_________________________END____________________________%%%

3
%%>>>>>>>>>>>>>>>>>>>>>>>>>>>>> MEMBER - D <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<%%
disp("For Member D - Non-Prestressed strands - 4# 1/2in bars (10in x 10in)")
Ap = 0.612 % strands area (in^2)
Ac = 100 - Ap % Net area of concrete (in^2)
fpbed = 0 % Prestressing (Ksi)
delta_ep = fpbed/Ep

%%%__________ STRESS IN CONCRETE - Member - D ____________%%%


disp("Stress in concrete fc - Member - D = ")
for i = 1:s
fc(i,1) = fcp* (2*(ecf(i,1)/ecp) - (ecf(i,1)/ecp) ^2);
if fc(i,1) > fcr
fc(i,1) = 0;
end
end
disp(double(fc))
%%%____________Stress in Non-prestressed strands - Member - D _________%%%
disp("Stress in strands - Member - D = ")
fp = zeros(7,1);
ep = ecc + delta_ep;
epf = ep- (Delta_T* alpha_csp);
for i = 1:s
fp(i,1) = Ep_eff* epf(i,1);
if fp(i,1) > fpy
fp(i,1) = fpy;
elseif fp(i,1) < -fpy
fp(i,1) = -fpy;
end
end
disp(fp)
%%%____________________Axial Load, N, (Kips) - Member - D ___________%%%
ND = zeros(s,1);
for i = 1:s
ND(i,1) = fc(i,1)* Ac + fp(i,1)* Ap;
end
disp("Axial Load - Member - D, N(Kips) = ")
double(ND);
disp(double(ND))
%%_______________________PLOTTING THE LOAD DEFORMATION CURVES____________%%
plot(ecc,NA,"-+",ecc,NB,"-o",ecc,NC,"-*",ecc,ND,"--")
title("Load-Deformation Curve-Long Term")
xlabel("Strain ec")
ylabel("Load (Kips)")
legend("A", "B", "C", "D","location", "southeast")
set(gca,'XGrid','on','YGrid','off')

%%%_________________________END OF Assignment ____________________________%%%

You might also like