cngf-pf-exp2

commit bcaea4a5b55a046416cad83424e21390535bc0ea
parent c5fa926d91065d4a893101511e77e32c9c80a193
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Wed,  4 Aug 2021 13:26:16 +0200

add root Makefile and strain distribution figure

Diffstat:
A
Makefile
 | 
24
++++++++++++++++++++++++
A
fig-strain_distribution/Makefile
 | 
31
+++++++++++++++++++++++++++++++
A
fig-strain_distribution/fig.gp
 | 
31
+++++++++++++++++++++++++++++++

3 files changed, 86 insertions(+), 0 deletions(-)
diff --git a/Makefile b/Makefile
@@ -0,0 +1,24 @@
+REPO = cngf-pf
+FIGDIRS != find . -maxdepth 1 -type d -name 'fig-*' | sed 's/\.\///'
+FIGURES = $(FIGDIRS:=.pdf)
+
+default: $(FIGURES)
+
+fig-rate_dependence.pdf: $(REPO)/$(REPO)
+	(cd $(@:.pdf=) && make)
+
+fig-strain_distribution.pdf: $(REPO)/$(REPO)
+	(cd $(@:.pdf=) && make)
+
+$(REPO)/$(REPO): $(REPO)/Makefile
+	make -C $(REPO)
+
+$(REPO)/Makefile:
+	git submodule update --init --recursive
+
+clean:
+	make -C $(REPO)/ clean
+	rm -f fig*.pdf
+	for d in $(FIGDIRS); do (cd $$d && make clean); done
+
+.PHONY: default figures clean
diff --git a/fig-strain_distribution/Makefile b/fig-strain_distribution/Makefile
@@ -0,0 +1,31 @@
+BIN = ../cngf-pf/cngf-pf
+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'
+
+../$(FIG).pdf: fig.gp strain_distribution_N10kPa.txt
+	gnuplot fig.gp > $@
+
+clean:
+	rm -f *.txt
+	rm -f ../$(FIG).pdf
+
+.PHONY: default clean
+
+
diff --git a/fig-strain_distribution/fig.gp b/fig-strain_distribution/fig.gp
@@ -0,0 +1,31 @@
+#!/usr/bin/env gnuplot
+
+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 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
+
+xdisp=0.74
+
+set label "d" at screen 0.02,0.95
+
+set key bottom right font ",9" samplen 1.0 maxrows 4
+
+pointsize=0.6
+
+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", \