 |
tAdd experiment with varying velocity - hansen-zoet-exp - comparison of lab and… |
 |
git clone git://src.adamsgaard.dk/hansen-zoet-exp |
 |
Log |
|
Files |
|
Refs |
|
Submodules |
|
LICENSE |
|
--- |
|
commit 8c8b651fafb09a5d73c64d9523566e1ac3238950 |
|
parent 56b744fa7d25b9f14cab23f2cdb764147208a31a |
|
Author: Anders Damsgaard <[email protected]> |
|
Date: Tue, 29 Oct 2019 11:22:30 +0100 |
|
|
|
Add experiment with varying velocity |
|
|
|
Diffstat: |
|
M fig-effective_stress/Makefile | 7 ++++--- |
|
A fig-velocity/Makefile | 42 +++++++++++++++++++++++++++++… |
|
A fig-velocity/fig-strain_distributi… | 31 +++++++++++++++++++++++++++… |
|
A fig-velocity/fig-velocity_stress.gp | 23 +++++++++++++++++++++++ |
|
|
|
4 files changed, 100 insertions(+), 3 deletions(-) |
|
--- |
|
diff --git a/fig-effective_stress/Makefile b/fig-effective_stress/Makefile |
|
t@@ -1,7 +1,7 @@ |
|
BIN = ../1d_fd_simple_shear/1d_fd_simple_shear |
|
FIG = $(shell basename $(CURDIR)) |
|
|
|
-default: ../$(FIG).pdf ../$(FIG)-mohr_coulomb.pdf |
|
+default: ../$(FIG)-strain_distribution.pdf ../$(FIG)-mohr_coulomb.pdf |
|
|
|
strain_distribution_N51kPa.txt: $(BIN) Makefile |
|
sh -c '\ |
|
t@@ -28,7 +28,7 @@ mohr_coulomb.txt: strain_distribution_N51kPa.txt |
|
tail -n 1 strain_distribution_N$${P}kPa.txt | cut -f3- >> $@; \ |
|
done' |
|
|
|
-../$(FIG).pdf: fig-strain_distribution.gp strain_distribution_N51kPa.txt |
|
+../$(FIG)-strain_distribution.pdf: fig-strain_distribution.gp strain_distribut… |
|
gnuplot fig-strain_distribution.gp > $@ |
|
|
|
../$(FIG)-mohr_coulomb.pdf: fig-mohr_coulomb.gp mohr_coulomb.txt |
|
t@@ -36,6 +36,7 @@ mohr_coulomb.txt: strain_distribution_N51kPa.txt |
|
|
|
clean: |
|
rm -f *.txt |
|
- rm -f ../$(FIG).pdf |
|
+ rm -f ../$(FIG)-strain_distribution.pdf |
|
+ rm -f ../$(FIG)-mohr_coulomb.pdf |
|
|
|
.PHONY: default clean |
|
diff --git a/fig-velocity/Makefile b/fig-velocity/Makefile |
|
t@@ -0,0 +1,42 @@ |
|
+BIN = ../1d_fd_simple_shear/1d_fd_simple_shear |
|
+FIG = $(shell basename $(CURDIR)) |
|
+ |
|
+default: ../$(FIG)-strain_distribution.pdf ../$(FIG)-velocity_stress.pdf |
|
+ |
|
+strain_distribution_vx3.17e-6ms.txt: $(BIN) Makefile |
|
+ sh -c '\ |
|
+ for vx in 3.17e-6 9.51e-6 1.90e-5 2.85e-5 3.81e-5; do \ |
|
+ ./$(BIN) \ |
|
+ --grain-size 1e-3 \ |
|
+ --friction-coefficient 0.661 \ |
|
+ --cohesion 30e3 \ |
|
+ --porosity 0.25 \ |
|
+ --density 1600 \ |
|
+ --set-shear-velocity $$vx \ |
|
+ --origo 0.0 \ |
|
+ --length 0.12 \ |
|
+ --normal-stress 101e3 \ |
|
+ --nonlocal-amplitude 0.48 \ |
|
+ --rate-dependence 0.9377 \ |
|
+ > strain_distribution_vx$${vx}ms.txt; \ |
|
+ done' |
|
+ |
|
+velocity_stress.txt: strain_distribution_vx3.17e-6ms.txt |
|
+ sh -c '\ |
|
+ rm -f $@; \ |
|
+ for vx in 3.17e-6 9.51e-6 1.90e-5 2.85e-5 3.81e-5; do \ |
|
+ tail -n 1 strain_distribution_vx$${vx}ms.txt | cut -f2- >> $@;… |
|
+ done' |
|
+ |
|
+../$(FIG)-strain_distribution.pdf: fig-strain_distribution.gp strain_distribut… |
|
+ gnuplot fig-strain_distribution.gp > $@ |
|
+ |
|
+../$(FIG)-velocity_stress.pdf: fig-velocity_stress.gp velocity_stress.txt |
|
+ gnuplot fig-velocity_stress.gp > $@ |
|
+ |
|
+clean: |
|
+ rm -f *.txt |
|
+ rm -f ../$(FIG)-strain_distribution.pdf |
|
+ rm -f ../$(FIG)-velocity_stress.pdf |
|
+ |
|
+.PHONY: default clean |
|
diff --git a/fig-velocity/fig-strain_distribution.gp b/fig-velocity/fig-strain_… |
|
t@@ -0,0 +1,31 @@ |
|
+#!/usr/bin/env gnuplot |
|
+ |
|
+reset |
|
+ |
|
+set terminal pdfcairo enhanced color size 7.5 cm, 10 cm |
|
+set multiplot layout 2,1 |
|
+ |
|
+#set lmargin 7.0 |
|
+#set bmargin 3.5 |
|
+#set rmargin 2.0 |
|
+#set tmargin 2.0 |
|
+ |
|
+#set xlabel "Normalized horizontal velocity, v_x [-]" |
|
+set xlabel "Horizontal velocity, v_x [m/a]" |
|
+set ylabel "Vertical position, z [m]" |
|
+set key bottom right font ",10" #samplen 0.9 |
|
+plot "strain_distribution_vx3.17e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title… |
|
+ "strain_distribution_vx9.51e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title… |
|
+ "strain_distribution_vx1.90e-5ms.txt" u ($2*365*24*3600):1 w l lw 2 title… |
|
+ "strain_distribution_vx2.85e-5ms.txt" u ($2*365*24*3600):1 w l lw 2 title… |
|
+ "strain_distribution_vx3.81e-5ms.txt" u ($2*365*24*3600):1 w l lw 2 title… |
|
+ |
|
+#set xlabel "Normalized horizontal velocity, v_x [-]" |
|
+set xlabel "Shear-strain rate, dv_x/dz [1/a]" |
|
+set ylabel "Vertical position, z [m]" |
|
+set key bottom right font ",10" #samplen 0.9 |
|
+plot "strain_distribution_vx3.17e-6ms.txt" u ($5*365*24*3600):1 w l lw 2 title… |
|
+ "strain_distribution_vx9.51e-6ms.txt" u ($5*365*24*3600):1 w l lw 2 title… |
|
+ "strain_distribution_vx1.90e-5ms.txt" u ($5*365*24*3600):1 w l lw 2 title… |
|
+ "strain_distribution_vx2.85e-5ms.txt" u ($5*365*24*3600):1 w l lw 2 title… |
|
+ "strain_distribution_vx3.81e-5ms.txt" u ($5*365*24*3600):1 w l lw 2 title… |
|
diff --git a/fig-velocity/fig-velocity_stress.gp b/fig-velocity/fig-velocity_st… |
|
t@@ -0,0 +1,23 @@ |
|
+#!/usr/bin/env gnuplot |
|
+ |
|
+reset |
|
+ |
|
+set terminal pdfcairo enhanced color size 7.5 cm, 8 cm |
|
+set multiplot layout 2,1 |
|
+ |
|
+#set lmargin 7.0 |
|
+#set bmargin 3.5 |
|
+#set rmargin 2.0 |
|
+#set tmargin 2.0 |
|
+ |
|
+#set yrange [0.0:0.73] |
|
+#set xrange [-0.13:1.13] |
|
+#set key bottom right #samplen 0.9 |
|
+ |
|
+#set xlabel "Shear velocity, v_x' [m/a]" |
|
+set ylabel "Shear stress, {/Symbol t} [kPa]" |
|
+plot "velocity_stress.txt" u ($1*365*24*3600):($3*$2/1e3) w lp title "" |
|
+ |
|
+set xlabel "Shear velocity, v_x [m/a]" |
|
+set ylabel "Shear friction, {/Symbol m} [-]" |
|
+plot "velocity_stress.txt" u ($1*365*24*3600):3 w lp title "" |
