data=load('BenchmarkRotationalSilicone_Cum30_2D-DIC_TS_MP(1x32x32_50ov)_Vy_y60-x-t.prm');

[m n]=size(data)

for k=1:2:n-1
    
subplot(2,3,1)
plot(data(:,1),data(:,k+1))

title('Across-strike extension (y = 60 mm)')
xlabel('position x (mm)')
ylabel('cumulative y-displacement (mm)')
xlim([-50 900])
ylim([-5 15])

hold on

end

data=load('BenchmarkRotationalSilicone_Cum30_2D-DIC_TS_MP(1x32x32_50ov)_Vx_y0-y-t.prm');

[m n]=size(data)

for k=1:2:n-1
    
subplot(2,3,3)
plot(data(:,1),data(:,k+1))

title('Along-strike deformation (y = 0 mm)')
xlabel('position x (mm)')
ylabel('cumulative x-displacement (mm)')
xlim([-50 900])
ylim([-20 0])

hold on

end

data=load('BenchmarkRotationalSilicone_Cum30_2D-DIC_TS_MP(1x32x32_50ov)_Vy_x750-y-t.prm');

[m n]=size(data)

for k=1:2:n-1
    
subplot(2,3,5)
plot(data(:,1),data(:,k+1))

title('Across-strike extension (x = 750 mm)')
xlabel('position y (mm)')
ylabel('cumulative y-displacement (mm)')
xlim([-200 200])
ylim([-15 15])

hold on

end



data=load('BenchmarkRotationalSilicone_Cum30_2D-DIC_TS_MP(1x32x32_50ov)_Vy_x50-y-t.prm');

[m n]=size(data)

for k=1:2:n-1
    
subplot(2,3,4)
plot(data(:,1),data(:,k+1))

title('Across-strike extension (x = 50 mm)')
xlabel('position y (mm)')
ylabel('cumulative y-displacement (mm)')
xlim([-200 200])
ylim([-15 15])

hold on

end