% post traitt plaxis loi de comportement

% resultats plaxis
M=dlmread('50b.vlt','\t',239,0);
cs='cd'; %cd, ou cuu ou oedo ou iso
e1=-M(:,3)*100;
e3=-M(:,4)*100;
s1=-M(:,9);
s3=-M(:,10);

p=(s1+2*s3)/3;
q=s1-s3;
ptot=q/3+p(1);
uw=ptot-p;
ev=e1+2*e3;

if strcmp(cs,'cd')==1 || strcmp(cs,'cuu')==1
    figure(1)
    hold on
    plot(e1,q);
    xlabel('\epsilon_1 (%)')
    ylabel('q (kP_a)')
    grid on
end

if strcmp(cs,'cd')==1
    figure(2)
    hold on
    plot(e1,ev);
    xlabel('\epsilon_1 (%)')
    ylabel('\epsilon_v (%)')
    grid on
end

if strcmp(cs,'cuu')==1
    figure(2)
    hold on
    plot(e1,uw);
    xlabel('\epsilon_1 (%)')
    ylabel('uw (kP_a)')
    grid on
    
    figure(3)
    hold on
    plot(p,q);
    xlabel('p (kP_a)')
    ylabel('q (kP_a)')
    grid on
end