clear all
close all

e0 = [0:.01:.25];
T0 = .225*1550*(exp(10.5*e0)-1);



T1 = [300 500:250:4000];
dL1 = [ 0 2 7 15 18 24 30 34 39 43 47 50 54 59 63 66];
P1 = [0 0 0 0 0 5 10 22 30 40 47 52 60 70 75 82];


T2 = [300 500:250:3500];
dL2 = [0 3 7 13 19 27 34 41 47 53 58 64 68 72];
dV2 = [0 50 70 175 500 600 800 900 900 1100 1000 1200 1000 1500];
V2 = [0];
for i = 2:length(dV2)
  V2(i) = V2(i-1)+dV2(i);
end

plot(e0,T0,dL1./(8*39),T1);  %,dL2./(8*39),T2);
[legend_h,object_h,plot_h,text_strings] = legend('T = 350*(exp(10.5*e)-1)','Test w/ Internal Volume Sealed',0);%,'Test w/ Internal Volume Vented',2);
set(object_h(1),'FontSize',14)

xlabel('e','FontSize',14);
ylabel('Tension (lbs)','FontSize',14);

print -dpsc2 StressStrain

figure
[AX,H1,H2] = plotyy(T2,V2/1000,T1,P1);
xlabel('Tension (lbs)');
set(get(AX(1),'Ylabel'),'String','Internal Volume Reduction at zero internal pressure (Liters)');
set(get(AX(2),'Ylabel'),'String','Internal Pressure at zero internal volume change (psi)');

print -dpsc2 PVVsTension

