Load and plot length-tension data
Three load and unload cycles were done on a 19.5cm restlength giantband. Graph shows that there is hysteresis. Note that the force data takes between 30 seconds and a minute to stabilise.
l0=19.5; % unstretched length of giantband. Appears to have a tollerance of 5mm T = readtable('GiantbandLTdata.xlsx'); offset=l0-T.length(1); leng=T.length+offset; % avoid using length as a variable as it is also a matlab function plot(leng,T.force) hold on indx1=[12 40 70 95 125 150]; x=leng(indx1); y=T.force(indx1); u=diff([leng(indx1) leng(indx1+5)]'); v=diff([T.force(indx1) T.force(indx1+5)]'); quiver(x,y,u',v') hold off xlabel('extension in cm') ylabel('force in Newtons') title('Three cycles of a giant band.') text(15,16,'Arrows show direction of loading and unloading');grid on
disp(sprintf('Publish source is %s.m on %s',mfilename,date));
Publish source is bandhysteresis.m on 07-Feb-2021