% Program to formulate global stiffness matrix
clc
clear


%load input.mat

%Soil_type
Type_of_soil = '';

fck = 20;
fy250 = 250;
fy450 = 450;
fy500 = 500;
fy = 500;
strainSteelMax = 10;
gammaConcrete = 1.5;
gammaSteel = 1.15;
ESteel = 200E03;

stressStrainSteel = [0 0];
stressStrainSteel = [stressStrainSteel; [fy250/ESteel*100 fy250]];
stressStrainSteel = [stressStrainSteel; [strainSteelMax fy250]];
stressStrainConcrete = [0:.01:.2]';
stressStrainConcrete(:,2) = 10 * fck * stressStrainConcrete(:,1) -25 * ...
  fck * stressStrainConcrete(:,1) .^2;
stressStrainConcrete =  [stressStrainConcrete; [.35 fck]];
% stressConcrete = [stressConcrete fck];
for ii = 1:rows(stressStrainConcrete)
  stressStrainConcrete(ii,1)
  stressStrainConcrete(ii,2)
end


%% Plot mode shapes

plotHangle = figure('visible', 'off')

plot(stressStrainConcrete(:,1), stressStrainConcrete(:,2),'-r')
hold on
%plot(stressStrainConcrete(:,1), stressStrainConcrete(:,2) * .67,'-g')
%plot(stressStrainConcrete(:,1), stressStrainConcrete(:,2) * .67 / ...
%  gammaConcrete,'-b')
plot(stressStrainSteel(:,1), stressStrainSteel(:,2),'-r')

hold off

% saveas(plotHangle, 'ModeShape.eps','eps')
%print (plotHangle, '-color',  'ModeShape.eps')
saveas(plotHangle, 'ModeShape.png','png')
%saveas(plotHangle, 'ModeShape.pdf')

% End of file
