cngf-pf-exp2

commit 37b9d310a582f98f857826e2c2b35355346074d8
parent 041625c73a95fbd7014bc2a45cc58d1a3d011308
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 17 Aug 2021 15:21:03 +0200

test new parameter values and plot porosity evolution

Diffstat:
M
fig-strain_distribution/Makefile
 | 
32
+++++++++++++++++++++++++-------
A
fig-strain_distribution/porosity-evolution.gp
 | 
35
+++++++++++++++++++++++++++++++++++
R
fig-strain_distribution/fig.gp -> fig-strain_distribution/strain-distrib.gp
 | 
0

3 files changed, 60 insertions(+), 7 deletions(-)
diff --git a/fig-strain_distribution/Makefile b/fig-strain_distribution/Makefile
@@ -1,29 +1,47 @@
 BIN = ../cngf-pf/cngf-pf
 FIG != basename $(PWD)
 
-default: ../$(FIG).pdf
+default: \
+	../$(FIG).pdf \
+	../$(FIG)-porosity-evolution.pdf \
 
 strain_distribution_s7.0m-s_N50kPa.txt: $(BIN) Makefile
 	for N in 6 50; do \
 		for s in 0.07 0.7 7.0; do \
+			id="s$${s}m-s_N$${N}kPa"; \
 			./$(BIN) \
 			-s "$$(awk -v s="$$s" 'BEGIN {print s/1000.0/60.0}')" \
+			-I "$$(awk -v s="$$s" 'BEGIN {print 0.1/(s/1000.0/60.0)/100.0}')" \
 			-e "$$(awk -v s="$$s" 'BEGIN {print 0.1/(s/1000.0/60.0)}')" \
-			-n $${N}e3 \
-			-p 0.275 \
+			-n "$${N}e3" \
+			-p 0.235 \
 			-d 0.004 \
 			-o -0.08 \
 			-L 0.07 \
-			> strain_distribution_s$${s}m-s_N$${N}kPa.txt; \
+			-T \
+			"$${id}" \
+			> "strain_distribution_s$${s}m-s_N$${N}kPa.txt"; \
+			> "$${id}-timeseries-mean.tsv"; \
+			> "$${id}-timeseries-min.tsv"; \
+			> "$${id}-timeseries-max.tsv"; \
+			for f in "$${id}".output*.txt; do \
+				mean < "$${f}" >> "$${id}-timeseries-mean.tsv"; \
+				min < "$${f}" >> "$${id}-timeseries-min.tsv"; \
+				max < "$${f}" >> "$${id}-timeseries-max.tsv"; \
+			done; \
 		done; \
 	done
 
-../$(FIG).pdf: fig.gp strain_distribution_s7.0m-s_N50kPa.txt
-	gnuplot fig.gp > $@
+../$(FIG).pdf: strain-distrib.gp strain_distribution_s7.0m-s_N50kPa.txt
+	gnuplot strain-distrib.gp > $@
+
+../$(FIG)-porosity-evolution.pdf: porosity-evolution.gp strain_distribution_s7.0m-s_N50kPa.txt
+	gnuplot porosity-evolution.gp > $@
 
 clean:
 	rm -f *.txt
-	rm -f ../$(FIG).pdf
+	rm -f *.tsv
+	rm -f ../$(FIG).pdf ../$(FIG)-porosity-evolution.pdf
 
 .PHONY: default clean
 
diff --git a/fig-strain_distribution/porosity-evolution.gp b/fig-strain_distribution/porosity-evolution.gp
@@ -0,0 +1,35 @@
+#!/usr/bin/env gnuplot
+
+reset
+
+set terminal pdfcairo enhanced color size 8 cm, 10 cm font ",10"
+set multiplot layout 2,1
+
+#set margins 8,2,3.2,1
+#set colorsequence podo
+
+set xlabel "{/:Normal Shear displacement [mm]}"
+set ylabel "{/:Normal Mean porosity [-]}" offset 1.0
+
+xscale=1
+yscale=xscale
+set xrange [-5:105]
+#set yrange [-90:0]
+#set xtics 0.2
+#set ytics 0.2
+
+#set label "d" at screen 0.02,0.95
+
+set key bottom left font ",9" samplen 1.0 maxrows 4
+
+#pointsize=0.6
+
+set title "N = 6 kPa"
+plot "s0.07m-s_N6kPa-timeseries-mean.tsv" u ($0*xscale):($7*yscale) w l t "{/:Bold CNGF-PF} 0.07 mm/min", \
+	"s0.7m-s_N6kPa-timeseries-mean.tsv" u ($0*xscale):($7*yscale) w l t "{/:Bold CNGF-PF} 0.7 mm/min", \
+	"s7.0m-s_N6kPa-timeseries-mean.tsv" u ($0*xscale):($7*yscale) w l t "{/:Bold CNGF-PF} 7.0 mm/min"
+
+set title "N = 50 kPa"
+plot "s0.07m-s_N50kPa-timeseries-mean.tsv" u ($0*xscale):($7*yscale) w l t "{/:Bold CNGF-PF} 0.07 mm/min", \
+	"s0.7m-s_N50kPa-timeseries-mean.tsv" u ($0*xscale):($7*yscale) w l t "{/:Bold CNGF-PF} 0.7 mm/min", \
+	"s7.0m-s_N50kPa-timeseries-mean.tsv" u ($0*xscale):($7*yscale) w l t "{/:Bold CNGF-PF} 7.0 mm/min"
diff --git a/fig-strain_distribution/fig.gp b/fig-strain_distribution/strain-distrib.gp