set xlabel 'ecrasement (%)'
set ylabel 'contrainte nominale (MPa)'

U_max = 0.02283875
ecrasement(t) = 100*U_max*t

r = 1.
S_0 = pi*r**2

set key left top
set grid
set title 'influence du frottement (reaction en bas) => oscillation avec coef=0.4'
p   "calcul_1x1elts.maple.coef_0p0"   u (ecrasement($1)):(abs($3/S_0)) t '1 elt frot=0'       w lp lc 0 dt 1 pt 4
rep "calcul_10x10elts.maple.coef_0p0" u (ecrasement($1)):(abs($3/S_0)) t '100 elts frot=0'    w lp lc 1 dt 1 pt 2
rep "calcul_10x10elts.maple.coef_0p1" u (ecrasement($1)):(abs($3/S_0)) t '100 elts frot=0.1'  w lp lc 2 dt 1 pt 2
rep "calcul_10x10elts.maple.coef_0p2" u (ecrasement($1)):(abs($3/S_0)) t '100 elts frot=0.2'  w lp lc 3 dt 1 pt 2
rep "calcul_10x10elts.maple.coef_0p3" u (ecrasement($1)):(abs($3/S_0)) t '100 elts frot=0.3'  w lp lc 4 dt 1 pt 2
rep "calcul_10x10elts.maple.coef_0p4" u (ecrasement($1)):(abs($3/S_0)) t '100 elts frot=0.4'  w lp lc 5 dt 1 pt 2
pause -1

set key left top
set grid
set title 'difference de reaction entre le haut et le bas du maillage'
p   "calcul_1x1elts.maple.coef_0p0"   u (ecrasement($1)):(abs($3/S_0)) t '1 elt frot=0'    w lp lc 0 dt 1 pt 4
rep "calcul_10x10elts.maple.coef_0p2" u (ecrasement($1)):(abs($6/S_0)) t '100 elts frot=0.2 (Nord)'w lp lc 1 dt 1 pt 2 lw 2
rep "calcul_10x10elts.maple.coef_0p2" u (ecrasement($1)):(abs($3/S_0)) t '100 elts frot=0.2 (Sud)'  w lp lc 1 dt 2 pt 2 lw 2
rep "calcul_10x10elts.maple.coef_0p4" u (ecrasement($1)):(abs($6/S_0)) t '100 elts frot=0.4 (Nord)' w lp lc 3 dt 1 pt 2 lw 2
rep "calcul_10x10elts.maple.coef_0p4" u (ecrasement($1)):(abs($3/S_0)) t '100 elts frot=0.4 (Sud)' w lp lc 3 dt 2 pt 2 lw 2
pause -1
