tchange parameters to lab values - cngf-pf-exp2 - experiments for second paper with continuum granular model
(HTM) git clone git://src.adamsgaard.dk/cngf-pf-exp2
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) Submodules
(DIR) LICENSE
---
(DIR) commit d1fc82817e52edc0718033ee34282be9507004a7
(DIR) parent 62b50dfc3030a206a816cdd6bba0fff53de8adf2
(HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Wed, 4 Aug 2021 17:58:30 +0200
change parameters to lab values
Diffstat:
M fig-strain_distribution/Makefile | 31 +++++++++++++++----------------
M fig-strain_distribution/fig.gp | 39 +++++++++++++++++--------------
2 files changed, 36 insertions(+), 34 deletions(-)
---
(DIR) diff --git a/fig-strain_distribution/Makefile b/fig-strain_distribution/Makefile
t@@ -3,23 +3,22 @@ FIG != basename $(PWD)
default: ../$(FIG).pdf
-strain_distribution_N10kPa.txt: $(BIN) Makefile
- /bin/sh -c '\
- for N in 10 20 40 60 80 120; do \
- ./$(BIN) \
- -m 0.4 \
- -s 0.0369 \
- -o 0.05 \
- -L 0.64 \
- -n $${N}e3 \
- -p 0.275 \
- -d 0.04 \
- -A 0.5 \
- -b 0.022 \
- -N > strain_distribution_N$${N}kPa.txt; \
- done'
+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 \
+ ./$(BIN) \
+ -s "$$(awk -v s="$$s" 'BEGIN {print s/1000.0/60.0}')" \
+ -e "$$(awk -v s="$$s" 'BEGIN {print 0.1/(s/1000.0/60.0)}')" \
+ -n $${N}e3 \
+ -p 0.275 \
+ -d 0.004 \
+ -L 0.09 \
+ -o -0.09 \
+ > strain_distribution_s$${s}m-s_N$${N}kPa.txt; \
+ done; \
+ done
-../$(FIG).pdf: fig.gp strain_distribution_N10kPa.txt
+../$(FIG).pdf: fig.gp strain_distribution_s7.0m-s_N50kPa.txt
gnuplot fig.gp > $@
clean:
(DIR) diff --git a/fig-strain_distribution/fig.gp b/fig-strain_distribution/fig.gp
t@@ -2,30 +2,33 @@
reset
-#set terminal pdfcairo enhanced color size 4.45 cm, 4.45 cm font ",10"
-set terminal pdfcairo enhanced color size 7 cm, 4.5 cm font ",10"
-
-set margins 8,2,3.2,1
+set terminal pdfcairo enhanced color size 8 cm, 16 cm font ",10"
+set multiplot layout 3,1
+#set margins 8,2,3.2,1
set colorsequence podo
-set xlabel "{/:Normal Shear displacement {/Symbol D}x [m]}"
-set ylabel "{/:Normal Vertical position z [m]}" offset 1.0
-
-set yrange [0.0:0.73]
-set xrange [-0.05:0.8]
-set xtics 0.2
-set ytics 0.2
+set xlabel "{/:Normal Shear displacement {/Symbol D}x [mm]}"
+set ylabel "{/:Normal Vertical position z [mm]}" offset 1.0
-xdisp=0.74
+xdisp=100.0
+#set yrange [0.0:0.73]
+#set xrange [0.0:xdisp]
+#set xtics 0.2
+#set ytics 0.2
-set label "d" at screen 0.02,0.95
+#set label "d" at screen 0.02,0.95
set key bottom right font ",9" samplen 1.0 maxrows 4
-pointsize=0.6
+#pointsize=0.6
+
+set title "N = 6 kPa"
+plot "strain_distribution_s0.07m-s_N6kPa.txt" u ($2*xdisp):1 w l t "{/:Bold CNGF-PF} 0.07 mm/min", \
+ "strain_distribution_s0.7m-s_N6kPa.txt" u ($2*xdisp):1 w l t "{/:Bold CNGF-PF} 0.7 mm/min", \
+ "strain_distribution_s7.0m-s_N6kPa.txt" u ($2*xdisp):1 w l t "{/:Bold CNGF-PF} 7.0 mm/min"
-plot "strain_distribution_N10kPa.txt" u ($2*xdisp):1 w l lw 2 lc 1 title "{/:Bold CNGF-PF} 10 kPa", \
- "strain_distribution_N20kPa.txt" u ($2*xdisp):1 w l lw 2 lc 2 title "{/:Bold CNGF-PF} 20 kPa", \
- "strain_distribution_N40kPa.txt" u ($2*xdisp):1 w l lw 2 lc 3 title "{/:Bold CNGF-PF} 40 kPa", \
- "strain_distribution_N80kPa.txt" u ($2*xdisp):1 w l lw 2 lc 4 title "{/:Bold CNGF-PF} 80 kPa", \
+set title "N = 50 kPa"
+plot "strain_distribution_s0.07m-s_N50kPa.txt" u ($2*xdisp):1 w l t "{/:Bold CNGF-PF} 0.07 mm/min", \
+ "strain_distribution_s0.7m-s_N50kPa.txt" u ($2*xdisp):1 w l t "{/:Bold CNGF-PF} 0.7 mm/min", \
+ "strain_distribution_s7.0m-s_N50kPa.txt" u ($2*xdisp):1 w l t "{/:Bold CNGF-PF} 7.0 mm/min"