# gnuplot: plot the case core files for 8620 carburized

set term qt enhanced font "hack,11"
set key noenhanced
#set colors classic
set colors default
set grid

set pointsize 1.8
set logscale xy

# Strain Life  plot --------------------------------------
set xlabel "2Nf"
set ylabel "Strain Ampl."
set yrange [0.0001:1.0]
set xrange [1:1.0E8]
set title "Compare 8620 Carb.Case files
plot "aisiIt095sae8620_Bhn469.html" u 2:1 w p pt 6, \
"aisiIt090sae8620_Bhn615.html" u 2:1 w p pt 6, \
"aisiIt089sae8620_Bhn688.html" u 2:1 w p pt 6, \
"aisiIt088sae8620_Bhn722.html" u 2:1 w p pt 6, \
"aisiIt087sae8620_Bhn670.html" u 2:1 w p pt 6, \
"aisiIt071sae8620_Bhn615.html" u 2:1 w p pt 6, \
"aisiIt071sae8620caseAll_fitted.html" u 2:1 w l lw 2, \
"aisiIt039sae8620_Bhn430core.html" u 2:1 w p pt 26

#Possible bad material:
#"aisiIt038sae8620_Bhn583.html" u 2:1 w p pt 6, \

#Compare the Core material tests vs fitted:
plot "aisiIt039sae8620_Bhn430core.html" u 2:1 w p pt 26 lc 1, "aisiIt039sae8620_Bhn430_fitted.html"  u 2:1 w lp pt 25 lc 1





# Neuber plot   -----------------------------------------
load "/home/faconle/plotinit"
set grid
set logscale xy
set pointsize 1.8
set xlabel "Reversals to failure, 2Nf"
set ylabel "Neuber Stress Ampl, sqrt(Sts * Str * E),  mpa"
set yrange [100:10000]
plot \
"aisiIt090sae8620_Bhn615.html" u 2:(sqrt($3*$1*202500)) w p pt 9, \
"aisiIt089sae8620_Bhn688.html" u 2:(sqrt($3*$1*206100)) w p pt 9, \
"aisiIt088sae8620_Bhn722.html" u 2:(sqrt($3*$1*201500)) w p pt 9, \
"aisiIt087sae8620_Bhn670.html" u 2:(sqrt($3*$1*198000)) w p pt 9, \
"aisiIt071sae8620_Bhn615.html" u 2:(sqrt($3*$1*207900)) w p pt 9, \
"aisiIt071sae8620caseAll_fitted.html" u 2:(sqrt($3*$1*180000)) w lp pt 33 lc -1

#Probably bad material:
#"aisiIt038sae8620_Bhn583.html" u 2:(sqrt($3*$1*202000)) w p pt 6, \

 set title "Compare AISI Case+Core Composite Iter. 62, 70, 39(Core)"
 plot \
  "aisiIt062sae8620case+core_Bhn615.html" u 2:(sqrt($3*$1*205400)) w p pt 6 lc -1, \
  "aisiIt070sae8620case+core_Bhn615.html" u 2:(sqrt($3*$1*207800)) w p pt 4 lc  1, \
  "aisiIt039sae8620_Bhn430.html" u 2:(sqrt($3*$1*180200)) w p pt 8 lc  7 \
  t "Core aisiIt039sae8620_Bhn430.html"



#!!grep E=  *html
#aisiIt038sae8620_Bhn583.html:#E= 202000
#aisiIt039sae8620_Bhn430core.html:#E= 180000
#aisiIt071sae8620_Bhn615.html:#E= 207900
#aisiIt087sae8620_Bhn670.html:#E= 198000
#aisiIt088sae8620_Bhn722.html:#E= 201500
#aisiIt089sae8620_Bhn688.html:#E= 206100
#aisiIt090sae8620_Bhn615.html:#E= 202500
#aisiIt095sae8620_Bhn469.html:#E= 207800



#Compare Strain life,  Axial vs bending ----------------------------------------------
#Compare 8620 carb. core Axial  vs  8615 carb core Bending
plot "aisiIt039sae8620_Bhn430core.html" u 2:1 w p pt 26 lc 1, "aisiIt039sae8620_Bhn430_fitted.html"  u 2:1 w lp pt 25 lc 1, "aisiIt156-sae8615simCoreBending.txt" u 2:1 w p pt 33


#Compare 8620 carb. case Axial vs  8615 carb case Bending
#  mv SAE_Iter163.txt aisiIt163-sae8620+carb-Bend-CA.txt
set xlabel "2Nf"
set ylabel "Strain Ampl."
set yrange [0.0001:1.0]
set xrange [1:1.0E8]
plot "aisiIt095sae8620_Bhn469.html" u 2:1 w p pt 6, \
"aisiIt090sae8620_Bhn615.html" u 2:1 w p pt 6, \
"aisiIt089sae8620_Bhn688.html" u 2:1 w p pt 6, \
"aisiIt088sae8620_Bhn722.html" u 2:1 w p pt 6, \
"aisiIt087sae8620_Bhn670.html" u 2:1 w p pt 6, \
"aisiIt071sae8620_Bhn615.html" u 2:1 w p pt 6, \
"aisiIt071sae8620caseAll_fitted.html" u 2:1 w l lw 2 lc -1, \
"aisiIt159-sae8615+carb-Bend-CA.txt" u 2:1 w p pt 33 lc 3, \
"aisiIt163-sae8615deepCarb-Bend-CA.txt" u 2:1 w p pt 26 lc 3, \
"aisiIt039sae8620_Bhn430core.html" u 2:1 w p pt 26 lc 1, "aisiIt039sae8620_Bhn430_fitted.html"  u 2:1 w lp pt 25 lc 1, \
"sim8615coreSmax=800.txt" u 2:1 w lp pt 9 lc 1



#Compare 8615  Deep case,  bend vs axial

set term qt enhanced font "hack,11"
set key noenhanced
#set colors classic
set colors default
set grid
set pointsize 1.8
set logscale xy

# Strain Life  plot --------------------------------------
set xlabel "2Nf"
set ylabel "Strain Ampl."
set yrange [0.0001:1.0]
set xrange [1:1.0E8]
set title "Compare 8615 Deep Case Bending vs Axial files
plot \
  "aisiIt155sae8615simCoreAxial.html"            u 2:1 w p pt 7 lc -1, \
  "aisiIt156-sae8615simCoreBending.html"         u 2:1 w p pt 4 lc -1 , \
  "aisiIt161sae8615deepCase+CoreAxialCA.html"    u 2:1 w p pt 8 lc  1, \
  "aisiIt163-sae8615deepCarb-Bend-CA.txt"        u 2:1 w p pt 4 lc  1

#  "aisiIt162sae8615deepCase+CoreAxialPOL.html"   u 2:1 w p pt 6 lc  1, \
#  "aisiIt164-sae8615deepCarb-Bend-POL.txt"       u 2:1 w p pt 4 lc  1

plot \
  "aisiIt155sae8615simCoreAxial.html"            u 2:1 w p pt 7 lc -1, \
  "aisiIt156-sae8615simCoreBending.html"         u 2:1 w p pt 4 lc -1 , \
  "aisiIt163-sae8615deepCarb-Bend-CA.txt"        u 2:1 w p pt 4 lc  1

plot \
  "aisiIt161sae8615deepCase+CoreAxialCA.html"    u 2:1 w p pt 8 lc  1, \
  "aisiIt163-sae8615deepCarb-Bend-CA.txt"        u 2:1 w p pt 4 lc  1

# 8620   Iter 62+70 +39(Core)
 set title "Compare AISI Case+Core Composite Iter. 62, 70, 39(Core)"
 plot \
  "aisiIt062sae8620case+core_Bhn615.html" u 2:1 w p pt 6 lc -1, \
  "aisiIt070sae8620case+core_Bhn615.html" u 2:1 w p pt 4 lc  1, \
  "aisiIt039sae8620_Bhn430.html" u 2:1 w p pt 8 lc  7 \
  t "Core aisiIt039sae8620_Bhn430.html"