%% ploting strain-stress-time diagrams
set(0,'DefaultAxesFontSize',20,'DefaultAxesFontName','helvetica','DefaultAxesTickLength',[.02 .02],'DefaultAxesLineWidth',2)
close all
clc
%%
figure
subplot(3,5,1)
plot(EtnaBasalt_EB2(:,1),EtnaBasalt_EB2(:,2),'LineWidth',3)
ylabel('strain, %')
yyaxis right
plot(EtnaBasalt_EB2(:,1),EtnaBasalt_EB2(:,3),'LineWidth',3)

ylim([200 400])
title('EB-2')
xlabel('time (min)')

subplot(3,5,2)
plot(EtnaBasalt_EB3(:,1),EtnaBasalt_EB3(:,2),'LineWidth',3)
ylabel('strain, %')
yyaxis right
plot(EtnaBasalt_EB3(:,1),EtnaBasalt_EB3(:,3),'LineWidth',3)

ylim([200 400])
title('EB-3')
xlabel('time (min)')

subplot(3,5,3)
plot(EtnaBasalt_EB4(:,1),EtnaBasalt_EB4(:,2),'LineWidth',3)
ylabel('strain, %')
yyaxis right
plot(EtnaBasalt_EB4(:,1),EtnaBasalt_EB4(:,3),'LineWidth',3)
ylabel('stress')
ylim([200 400])
title('EB-4')
xlabel('time (min)')

subplot(3,5,4)
plot(IcelandBasalt_OR2_M(:,1),IcelandBasalt_OR2_M(:,2),'LineWidth',3)
ylabel('strain, %')
yyaxis right
plot(IcelandBasalt_OR2_M(:,1),IcelandBasalt_OR2_M(:,3),'LineWidth',3)
ylabel('stress')
ylim([50 100])
title('OR2-M')
xlabel('time (min)')

subplot(3,5,5)
plot(Inada_granite(:,1),Inada_granite(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Inada Granite')
xlabel('time (min)')

subplot(3,5,6)
plot(Thala_limestone(:,1),Thala_limestone(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Thala limestone')
xlabel('time (min)')

subplot(3,5,7)
plot(Darley_Dale_84(:,1),Darley_Dale_84(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Darley Dale, 84%')
xlabel('time (min)')

subplot(3,5,8)
plot(Darley_Dale_88(:,1),Darley_Dale_88(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Darley Dale, 88%')
xlabel('time (min)')

subplot(3,5,9)
plot(Darley_Dale_93(:,1),Darley_Dale_93(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Darley Dale, 93%')
xlabel('time (min)')

subplot(3,5,10)
plot(Etna_basalt_80(:,1),Etna_basalt_80(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Etna Basalt, 80%')
xlabel('time (min)')

subplot(3,5,11)
plot(Etna_basalt_86(:,1),Etna_basalt_86(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Etna Basalt, 86%')
xlabel('time (min)')

subplot(3,5,12)
plot(Etna_basalt_93(:,1),Etna_basalt_93(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Etna Basalt, 93%')
xlabel('time (min)')

subplot(3,5,13)
plot(Etna_basalt_97(:,1),Etna_basalt_97(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Etna Basalt, 97%')
xlabel('time (min)')

subplot(3,5,14)
plot(Glass_1(:,1),Glass_1(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Glass 1')
xlabel('time (min)')

subplot(3,5,15)
plot(Glass_3(:,1),Glass_3(:,2),'LineWidth',3)
ylabel('strain, %')
%ylim([50 100])
title('Glass 3')
xlabel('time (min)')



