### Setup terminal ###
reset session
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"



### Setup grille ###
set style line 81 lt -1 lc rgb "gray90" lw 0.5
set grid ytics mytics ls 81
set grid xtics mytics ls 81
set mxtics 1
set mytics 1
set grid x mx y my



### Setup fichiers input ###
array data[2]
    data[1] = "../0.2_Donnees/1-Moy_libre_ChDe.txt"
    data[2] = "Libre_hyper_princ.maple"





### Setup plot 1 ###
# Key & labels
set xrange [0.:]
set xlabel "Déformation [-]"
#set yrange [-5.:25.]
set ylabel "Contrainte [MPa]"

# Plot
set output "2_Figures/1_contrainte_defo.png"
show output
plot data[1] u 2:3 t "Exp" w l lc 1 lw 2, \
    data[2] u ((23.97-$41)/7.97):($3/955) t "num dim vrai" w l lc 2 lw 2



### Setup plot 2 ###
# Key & labels
# Key & labels
set xrange [0.:]
set xlabel "Déformation [-]"
set yrange [-0.5:4.]
set ylabel "Contrainte [MPa]"

# Plot
set output "2_Figures/2_contributions.png"
show output
plot data[1] u 2:3 t "Exp" w l lc 1 lw 2, \
    data[2] u (1-($41-16)/7.97):($3/955) t "num dim vrai" w l lc 2 lw 2, \
	"" u (1-($41-16)/7.97):(-$77*$71) t "Hyper sph" w l lc 4 lw 2, \
	"" u (1-($41-16)/7.97):(-$83*$72) t "Hyper dév" w l lc 6 lw 2, \
	"" u (1-($41-16)/7.97):(-$89*$73) t "Hyste sph" w l lc 7 lw 2, \
	"" u (1-($41-16)/7.97):(-$95*$74) t "Hyste dev" w l lc 8 lw 2



### Setup plot 3 ###
# Key & labels
set xrange [0.:0.03]
set xlabel "Déformation [-]"
set yrange [0.:2.]
set ylabel "Contrainte [MPa]"

# Plot
set output "2_Figures/2.2_contributions.png"
show output
plot data[1] u 2:3 t "Exp" w l lc 1 lw 2, \
    data[2] u (1-($41-16)/7.97):($3/955) t "num dim vrai" w l lc 2 lw 2, \
	"" u (1-($41-16)/7.97):(-$77*$71) t "Hyper sph" w l lc 4 lw 2, \
	"" u (1-($41-16)/7.97):(-$83*$72) t "Hyper dév" w l lc 6 lw 2, \
	"" u (1-($41-16)/7.97):(-$89*$73) t "Hyste sph" w l lc 7 lw 2, \
	"" u (1-($41-16)/7.97):(-$95*$74) t "Hyste dév" w l lc 8 lw 2



### Setup plot 4 ###
# Key & labels
set xrange [1.:0.88]
set xlabel "V_{vol} [-]"
set yrange [-5.:25.]
set ylabel "Contrainte [MPa]"

# Plot
set output "2_Figures/2.3_contributions.png"
show output
plot data[1] u 2:3 t "Exp" w l lc 1 lw 2, \
    data[2] u ($65):($3/955) t "num dim vrai" w l lc 2 lw 2, \
	"" u ($65):(-$77*$71) t "Hyper sph" w l lc 4 lw 2, \
	"" u ($65):(-$83*$72) t "Hyper dév" w l lc 6 lw 2, \
	"" u ($65):(-$89*$73) t "Hyste sph" w l lc 7 lw 2, \
	"" u ($65):(-$95*$74) t "Hyste dév" w l lc 8 lw 2



### Setup plot 5 ###
# Key & labels
set xrange [1.:0.88]
set xlabel "Déformation [-]"
set yrange [-0.05:1.05]
set ylabel "Pondération [-]"
unset xrange


# Plot
set output "2_Figures/3_ponderations.png"
show output
plot data[2] u ($62):($71) t "Hyper sph" w l lc 4 lw 2, \
	"" u 61:($72) t "Hyper dév" w l lc 6 lw 2, \
	"" u ($62):($73) t "Hyste sph" w l lc 7 lw 2, \
	"" u 61:($74) t "Hyste dév" w l lc 8 lw 2



### Setup plot 6 ###
# Key & labels
set xrange [1.:0.88]
set xlabel "Déformation sphérique [-]"
set yrange [0.:1.2]
set ylabel "Pondération [-]"
unset xrange

# Plot
set output "2_Figures/3.1_ponderations.png"
show output
plot data[2] u ($62):($71) t "Hyper sph" w l lc 4 lw 2, \
	"" u ($62):($73) t "Hyste sph" w l lc 7 lw 2



### Setup plot 7 ###
# Key & labels
set xrange [1.:0.88]
set xlabel "Déformation déviatorique [-]"
set yrange [0.:1.2]
set ylabel "Pondération [-]"
unset xrange

# Plot
set output "2_Figures/3.2_ponderations.png"
show output
plot data[2]  u 61:($72) t "Hyper dév" w l lc 6 lw 2, \
	"" u 61:($74) t "Hyste dév" w l lc 8 lw 2



### Setup plot 8 ###
# Key & labels
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}"

# Plot
set output "2_Figures/4_radial.png"
show output
plot data[2] u ($35/22.945-1):($3/955) t "exterieur" w lp lc 1 lw 1 ps 1, \
    "" u ($45/15.06-1):($3/955) t "interieur" w lp lc 2 lw 1 ps 1



### Setup plot 9 ###
# Key & labels
set xrange [0.:]
set xlabel "Déformation bague radiale [-]"
set yrange [-5.:25.]
set ylabel "Contrainte [MPa]"
unset xrange
unset yrange
unset format x

# Plot
set output "2_Figures/5_radial.png"
show output
plot data[2] u ($35):($3/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



### Setup plot 10 ###
# Key & labels
set xrange [0.:]
set xlabel "Déformation [-]"
set yrange [-5.:25.]
set ylabel "Contrainte [MPa]"

# Plot
set output "2_Figures/6_contributions_brut.png"
show output
plot data[1] u 2:3 t "Exp" w l lc 1 lw 2, \
    data[2] u (1-($41-16)/7.97):($3/955) t "num dim vrai" w l lc 2 lw 2, \
	"" u (1-($41-16)/7.97):(-$77) t "Hyper sph" w l lc 4 lw 2, \
	"" u (1-($41-16)/7.97):(-$83) t "Hyper dév" w l lc 6 lw 2



### Setup plot 11 ###
# Key & labels
set xrange [0.:25.]
set xlabel "Contrainte [MPa]"
set yrange [0.:4.]
set ylabel "Pondération [-]"

# Plot
set output "2_Figures/7_ponderations_sigma.png"
show output
plot data[2] u ($3/955):($71) t "Hyper sph" w l lc 4 lw 2, \
	"" u ($3/955):($72) t "Hyper dév" w l lc 6 lw 2



### Setup plot 2 ###
# Key & labels
set xrange [0.:0.12]
set xlabel "Déformation [-]"
set yrange [0.:4.]
set ylabel "Pondération [-]"

# Plot
set output "2_Figures/8_ponderations_eps.png"
show output
plot data[2] u (1-($41-16)/7.97):($71) t "Hyper sph" w l lc 4 lw 2, \
	"" u (1-($41-16)/7.97):($72) t "Hyper dév" w l lc 6 lw 2



### Setup plot 13 ###
# Key & labels
unset xrange
set xlabel "Temps [-]"
unset yrange
set ylabel "V vol [-]"

# Plot
set output "2_Figures/9_Vvol.png"
show output
plot data[2] u 1:($65) t "V vol" w l lc 1 lw 2



### Setup plot 14 ###
# Key & labels
set xlabel "Temps [-]"
set ylabel "V forme [-]"

# Plot
set output "2_Figures/10_Forme.png"
show output
plot data[2] u 1:(($41-16)/(($35-15.06)-($45-15.06))) t "V forme" w l lc 1 lw 2



### Setup plot 15 ###
# Key & labels
set arrow from 0.87, 0.87 to 1, 1 nohead dt 2
set xlabel "V forme [-]"
set xrange [1.002:0.86]
set ylabel "V vol [-]"
set yrange [1.002:0.86]

F0 = (7.97)/(22.93-15.06)

# Plot
set output "2_Figures/11_evolution.png"
show output
plot data[2] u ((($41-16)/(($35-15.06)-($45-15.06))) / F0):($69) t "V vol vs V forme" w l lc 1 lw 2
unset arrow



### Setup plot 16 ###
# Key & labels
set xlabel "Deformation Mises [-]"
set xrange [1.025:0.88]
set ylabel "Contrainte Mises [MPa]"
set yrange [0.88:1.005]
unset xrange
unset yrange

# Plot
set output "2_Figures/12_mises.png"
show output
plot data[2] u ($61):($59/3) t "V vol vs V forme" w l lc 1 lw 2




### Fin fichier ###
unset output
