manus_continuum_granular1_exp

commit 16a5c0dc2c937b07e389530ace60f8162d3b2728
parent c3a4389a4834fb680012f552ad7b8dd89a89319b
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 19 Nov 2019 10:47:06 +0100

Play with parameter choice for strain distribution

Diffstat:
M
fig-strain_distribution/Makefile
 | 
5
++++-
M
fig-strain_distribution/fig.gp
 | 
8
++++----

2 files changed, 8 insertions(+), 5 deletions(-)
diff --git a/fig-strain_distribution/Makefile b/fig-strain_distribution/Makefile
@@ -10,7 +10,10 @@ strain_distribution_N10kPa.txt: $(BIN) Makefile
 		--friction-coefficient 0.4 \
 		--set-shear-velocity 0.0369 \
 		--origo 0.05 --length 0.64 --normal-stress $${P}e3 \
-		--nonlocal-amplitude 0.48 \
+		--porosity 0.275 \
+		--grain-size 0.04 \
+		-A 0.2 \
+		-b 0.2 \
 		--normalize > strain_distribution_N$${P}kPa.txt; \
 	done'
 
diff --git a/fig-strain_distribution/fig.gp b/fig-strain_distribution/fig.gp
@@ -16,18 +16,18 @@ set yrange [0.0:0.73]
 #set xrange [-0.13:1.13]
 set xrange [-0.1:0.8]
 
-xdisp=0.737
+xdisp=0.74
 
-set key bottom right #samplen 0.9
+set key bottom right font ",10" #samplen 0.9
 
 plot "strain_distribution_N10kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#ee2430" title "{/Symbol s}_n' 10 kPa", \
 	"damsgaard2013-sheardisp/sheardisp10-binned.dat" w p lw 2 lc "#ee2430" title "DEM, {/Symbol s}_n' 10 kPa", \
 	"strain_distribution_N40kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#2f54ad" title "{/Symbol s}_n' 40 kPa", \
 	"damsgaard2013-sheardisp/sheardisp40-binned.dat" w p lw 2 lc "#2f54ad" title "DEM, {/Symbol s}_n' 40 kPa", \
-	"strain_distribution_N80kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#6ccfd9" title "{/Symbol s}_n' 80 kPa", \
-	"damsgaard2013-sheardisp/sheardisp80-binned.dat" w p lw 2 lc "#6ccfd9" title "DEM, {/Symbol s}_n' 80 kPa", \
 	"strain_distribution_N120kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#231f20" title "{/Symbol s}_n' 120 kPa", \
 	"damsgaard2013-sheardisp/sheardisp120-binned.dat" w p lw 2 lc "#231f20" title "DEM, {/Symbol s}_n' 120 kPa"
 
+	#"strain_distribution_N80kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#6ccfd9" title "{/Symbol s}_n' 80 kPa", \
+	#"damsgaard2013-sheardisp/sheardisp80-binned.dat" w p lw 2 lc "#6ccfd9" title "DEM, {/Symbol s}_n' 80 kPa", \
 	#"strain_distribution_N20kPa.txt" u 2:1 w l lw 2 lc "#53bf3a" title "20 kPa", \
 	#"damsgaard2013-sheardisp/sheardisp20-binned.dat" w p lw 2 lc "#53bf3a" title "DEM, 20 kPa", \