clear 
clc

% This script plots the comparisons for Assignment #3

res_case1=importdata('elastic_linGeom.txt');
res_case2=importdata('elastic_nonlinGeom.txt');
res_case3a=importdata('springs_linGeom.txt');
res_case3b=importdata('fiber_linGeom.txt');
res_case4a=importdata('springs_nonlinGeom.txt');
res_case4b=importdata('fiber_nonlinGeom.txt');

h2=figure;
plot(res_case1(:,1),res_case1(:,2)/1e3)
hold on
plot(res_case2(:,1),res_case2(:,2)/1e3)
plot(res_case3a(:,1),res_case3a(:,2)/1e3)
plot(res_case3b(:,1),res_case3b(:,2)/1e3)
plot(res_case4a(:,1),res_case4a(:,2)/1e3)
plot(res_case4b(:,1),res_case4b(:,2)/1e3)
legend('Case 1','Case 2','Case 3a','Case 3b','Case 4a','Case 4b','Location','best')
xlabel('u_b [mm]')
ylabel('V_{tot} [kN]')
xlim([0,1000])
ylim([0,2000])
grid on

plot_settings_ASCE(h2)