elongation_al(eps_almansi) = 1./sqrt(1. - 2.*eps_almansi)
eps_log(eps_almansi) = log(elongation_al(eps_almansi))

set title "traction uniaxiale X"
set xlabel "eps log 11 (-)"
set ylabel "SIG11 (MPa)"
set y2label "eps log 22 (-)"
set y2tics
set key outside
p   "calcul_HZ_princ.maple" u (eps_log($5)):11 t 'ref SIG11' w lp lc 0 pt 2 ps 1.2 dt 1 lw 2
rep "calcul_HZ_princ.maple" u (eps_log($5)):(eps_log($6)) axis x1y2 t 'ref eps log 22' w lp lc 0 pt 2 ps 1.2 dt 3 lw 2
rep "calcul_Umat_princ.maple" u (eps_log($5)):11 t 'Umat SIG11' w lp lc 1 pt 4 dt 1
rep "calcul_Umat_princ.maple" u (eps_log($5)):(eps_log($6)) axis x1y2 t 'Umat eps log 22' w lp lc 1 pt 4 dt 3
pause -1
set xlabel "temps (s)"
p   "calcul_HZ_princ.maple" u 1:11 t 'ref SIG11' w lp lc 0 pt 2 ps 1.2 dt 1 lw 2
rep "calcul_HZ_princ.maple" u 1:(eps_log($6)) axis x1y2 t 'ref eps log 22' w lp lc 0 pt 2 ps 1.2 dt 3 lw 2
rep "calcul_Umat_princ.maple" u 1:11 t 'Umat SIG11' w lp lc 1 pt 4 dt 1
rep "calcul_Umat_princ.maple" u 1:(eps_log($6)) axis x1y2 t 'Umat eps log 22' w lp lc 1 pt 4 dt 3
pause -1
