set term x11 font ",12"

elongation_al(eps_almansi) = 1./sqrt(1. - 2.*eps_almansi)
eps_log(eps_almansi) = log(elongation_al(eps_almansi))
###IIeps(EPS_I,EPS_II, NU) = eps_al(EPS_I)**2 + eps_al(EPS_II)**2 + eps_al(-NU/(1.-NU)*(EPS_I+EPS_II))**2
IIeps_EPS33nulle(EPS11,EPS22,EPS12) = EPS11**2 + EPS22**2 + 2.*EPS12**2

set grid
set key left top Left reverse
set xlabel 'temps'
set ylabel 'deformation log'
p   "traction_princ.maple" u 1:(eps_log($5)) every 18 t 'def log(EPS11)' w p lc 0 pt 4
rep "traction_princ.maple" u 1:25 t 'DEF ASSO LOI 11' w l dt 1 lc 0 lw 3
rep "traction_princ.maple" u 1:(eps_log($6)) every 18 t 'def log(EPS22)' w p lc 1 pt 4
rep "traction_princ.maple" u 1:27 t 'DEF ASSO LOI 22' w l dt 1 lc 1 lw 3
rep "traction_princ.maple" u 1:(eps_log($22)) every 18 t 'def log(DEF EPAISSEUR)' w p lc 3 pt 4
rep "traction_princ.maple" u 1:30 t 'DEF ASSO LOI 33' w l dt 1 lc 3 lw 2
rep "traction_princ.maple" u 1:(log($21/0.038)) t 'def log 33 via epaisseur' w l dt 2 lc 3 lw 2
pause -1


set xtics 0.05
set key center center Left reverse
set xlabel 'IIeps (EPS33 = 0)'
set ylabel 'module actuel'
p   "traction_princ.maple" u (IIeps_EPS33nulle($5,$6,$7)):23 t 'K = Kc/3' w l lc 0 dt 1 lw 2
rep "traction_princ.maple" u (IIeps_EPS33nulle($5,$6,$7)):24 t 'mu = 2G' w l lc 1 dt 1 lw 2
pause -1
set xtics autofreq
