bt400.m

This is MATLAB program

% Chuan bi so lieu cho Anh
fcA=40;
gammamA=1.5;
bienA=2.44*sqrt(fcA/gammamA)/10000;
khacA=0.67*fcA/gammamA;
x1=0:0.0001:0.0035;
paramsA(1)=bienA;
paramsA(2)=fcA;
paramsA(3)=gammamA;
paramsA(4)=khacA;
% Chuan bi so lieu cho Phap
fcP=40;
gammamP=1.5;
bienP=0.002;
khacP=0.85*fcP/gammamP;
paramsP(1)=bienP;
paramsP(2)=fcP;
paramsP(3)=gammamP;
paramsP(4)=khacP;
% Chuan bi so lieu cho Nhat
fcN=40;
gammamN=1.3;
bienN=0.002;
khacN=0.85*fcN/gammamN;
paramsN(1)=bienN;
paramsN(2)=fcN;
paramsN(3)=gammamN;
paramsN(4)=khacN;
% Chuan bi so lieu cho My
xMy=0:0.0001:0.0035;
fcM=41.3685; %quy doi 6000 psi thanh N/mm2
gammamM=1.1;
bienM=1.71*fcM/(57000*sqrt(6000)*0.006894757293);%quy doi psi - mm2
khacM=0.90*fcM;
paramsM(1)=bienM;
paramsM(2)=fcM;
paramsM(3)=gammamM;
paramsM(4)=khacM;
% Ve do thi cua Anh
plot(x1,Anh(x1,paramsA),'o-r','LineWidth',1,'MarkerSize',4);
text(0.0004,Anh(0.0004,paramsA),'Anh\rightarrow',...
	'HorizontalAlignment','right');
hold on
% Ve do thi cua Phap
plot(x1,Phap(x1,paramsP),'g*-','LineWidth',1,'MarkerSize',4);
text(0.0007,Phap(0.0007,paramsP),'\leftarrowPhap');
%text(0.0003,Phap(0.0003,paramsP),'\bullet\leftarrow\fontname...
%{times}fcu/gama{\ite}^{0.67{\itt}}at{\itt}=0.0003);
hold on
% Ve do thi cua Nhat
plot(x1,Nhat(x1,paramsN),'square-k','LineWidth',1,'MarkerSize',4);
text(0.001,Nhat(0.001,paramsN),'Nhat\rightarrow',...
   'HorizontalAlignment','right');
hold on
% Ve do thi cua My
plot(xMy,My1(xMy,paramsM),'x-r','LineWidth',1,'MarkerSize',4);
text(0.0008,My1(0.0008,paramsM),'My\rightarrow',...
 	'HorizontalAlignment','right');
axis('auto');
hold off;
xlabel('Strain');
ylabel('Stress - N/mm2');
title('Short term design stress - strain curve for concrete C40');
%grid off;
%zoom on;