set terminal pngcairo size 1080,1080

set encoding utf8
set key left top
set tics font ",22"
set xlabel font ",22"
set ylabel font ",22"
set title font ",22"
set key font ",22"

set style line 81 lt 0 lc rgb "#808080" lw 2
set grid ytics mytics
set grid xtics mytics
set mxtics 2 
set mytics 2 
set grid

set xrange [0.:]
set xlabel "Déformation [-]"
set yrange [-5.:25.]
set ylabel "Contrainte [MPa]"

set output "2_Figures/1_contrainte_defo.png"
plot "1_Donnees/Essais_CD_GEM_essai_2.txt" u 5:6 t "Exp" w l lc 1 lw 2, \
    "Essais_contact_princ.maple" u (1-($48-16)/8.05):($6/955) t "num dim vrai" w l lc 2 lw 2, \
    "1_Donnees/Quadra_princ.maple" u ((24. - $111)/8):($12/955) t "num dim theo" w l lc 3 lw 2



set xrange [0.:]
set xlabel "Déformation [-]"
set yrange [-5.:25.]
set ylabel "Contrainte [MPa]"

set output "2_Figures/2_contributions.png"
plot "1_Donnees/Essais_CD_GEM_essai_2.txt" u 5:6 t "Exp" w l lc 1 lw 2, \
    "Essais_contact_princ.maple" u (1-($48-16)/8.05):($6/955) t "num dim vrai" w l lc 2 lw 2, \
	"" u (1-($48-16)/8.05):(-$74*$68) t "Hyper sph" w l lc 4 lw 2, \
	"" u (1-($48-16)/8.05):(-$80*$69) t "Hyper dév" w l lc 6 lw 2, \
	"" u (1-($48-16)/8.05):(-$86*$70) t "Hyste sph" w l lc 7 lw 2, \
	"" u (1-($48-16)/8.05):(-$92*$71) t "Hyste dev" w l lc 8 lw 2



set xrange [1.:0.88]
set xlabel "V vol [-]"
set yrange [0.:4.]
set ylabel "Pondération [-]"

set output "2_Figures/3_ponderations.png"
plot "Essais_contact_princ.maple" u 66:($68) t "Hyper sph" w l lc 4 lw 2, \
	"" u 66:($69) t "Hyper dév" w l lc 6 lw 2, \
	"" u 66:($70) t "Hyste sph" w l lc 7 lw 2, \
	"" u 66:($71) t "Hyste dev" w l lc 8 lw 2
	



set key center top
set xrange [0.:]
set xlabel "Déformation bague radiale [-]"
set yrange [-5.:25.]
set ylabel "Contrainte [MPa]"
unset xrange
unset yrange
set format x "%.0s*10^{%T}"

set output "2_Figures/4_radial.png"
plot "Essais_contact_princ.maple" u ($42/22.9999-1):($6/955) t "exterieur" w lp lc 1 lw 1 ps 1, \
    "" u ($52/15.0001-1):($6/955) t "interieur" w lp lc 2 lw 1 ps 1




set xrange [0.:]
set xlabel "Déformation bague radiale [-]"
set yrange [-5.:25.]
set ylabel "Contrainte [MPa]"
unset xrange
unset yrange
unset format x

set output "2_Figures/5_radial.png"
plot "Essais_contact_princ.maple" u ($42):($6/955) t "exterieur" w lp lc 1 lw 1 ps 1, \
    #"" u ($52/15.0001-1):($6/955) t "interieur" w lp lc 2 lw 1 ps 1
