cngf-pf

commit f32aea8be407e1d445ddc7fff296966d7b010605
parent 948b628088cd1e756595eb0e9e9b423e4e6b6d65
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 11 Apr 2019 13:21:04 +0200

Add stress vs. strain rate analysis

Diffstat:
M
1d_fd_simple_shear.png
 | 
0

A
1d_fd_simple_shear_rheology.gp
 | 
24
++++++++++++++++++++++++
A
1d_fd_simple_shear_rheology.png
 | 
0

M
Makefile
 | 
16
+++++++++++++++-

4 files changed, 39 insertions(+), 1 deletion(-)
diff --git a/1d_fd_simple_shear.png b/1d_fd_simple_shear.png
Binary files differ.
diff --git a/1d_fd_simple_shear_rheology.gp b/1d_fd_simple_shear_rheology.gp
@@ -0,0 +1,24 @@
+#!/usr/bin/env gnuplot
+
+set terminal pngcairo color size 18.6 cm, 17.6 cm
+set output "1d_fd_simple_shear_rheology.png"
+set xlabel "Applied friction, mu = tau/N [-]"
+set ylabel "Shear velocity, v [m/s]" offset 2
+
+set xrange [0.0:0.7]
+set yrange [0.0:0.01]
+set key top left #samplen 0.9
+
+set style line 1 linetype 1 linewidth 3 pointtype 1 pointsize 1
+
+plot "1d_fd_simple_shear_rheology_b0.20.txt" u 1:2 w lp lw 2 t "b = 0.2", \
+         "1d_fd_simple_shear_rheology_b0.30.txt" u 1:2 w lp lw 2 t "b = 0.3", \
+         "1d_fd_simple_shear_rheology_b0.40.txt" u 1:2 w lp lw 2 t "b = 0.4", \
+         "1d_fd_simple_shear_rheology_b0.50.txt" u 1:2 w lp lw 2 t "b = 0.5", \
+         "1d_fd_simple_shear_rheology_b0.60.txt" u 1:2 w lp lw 2 t "b = 0.6", \
+         "1d_fd_simple_shear_rheology_b0.70.txt" u 1:2 w lp lw 2 t "b = 0.7", \
+         "1d_fd_simple_shear_rheology_b0.80.txt" u 1:2 w lp lw 2 t "b = 0.8", \
+         "1d_fd_simple_shear_rheology_b0.90.txt" u 1:2 w lp lw 2 t "b = 0.9", \
+         "1d_fd_simple_shear_rheology_b1.00.txt" u 1:2 w lp lw 2 t "b = 1.0"
+
+set xtics norotate # Restore defaults
diff --git a/1d_fd_simple_shear_rheology.png b/1d_fd_simple_shear_rheology.png
Binary files differ.
diff --git a/Makefile b/Makefile
@@ -4,7 +4,8 @@ SRC=$(wildcard *.c)
 OBJ=$(patsubst %.c,%.o,$(SRC))
 HDR=$(wildcard *.h)
 
-default: 1d_fd_simple_shear.png
+default: 1d_fd_simple_shear.png \
+	1d_fd_simple_shear_rheology.png
 
 1d_fd_simple_shear: $(OBJ) $(HDR)
 	$(CC) $(LDFLAGS) $(OBJ) -o $@
@@ -18,6 +19,18 @@ default: 1d_fd_simple_shear.png
 	./$< -P 120e3 -N > $<_P120kPa.txt
 	gnuplot $<.gp > $@
 
+1d_fd_simple_shear_rheology.png: 1d_fd_simple_shear 1d_fd_simple_shear_rheology.gp
+	/bin/bash -c '\
+	for b in $$(seq 0.20 0.10 1.00); do \
+		out="$<_rheology_b$$b.txt"; \
+		rm "$$out"; \
+		for t in $$(seq 0.01 0.01 0.8); do \
+			printf "$$t\t" >> "$$out"; \
+			./$< -P 20e3 --stress-ratio $$t -b $$b | \
+			tail -n 1 | cut -f2 >> "$$out"; \
+	done; done'
+	gnuplot $<_rheology.gp > $@
+
 .PHONY: watch
 watch:
 	echo $(SRC) $(HDR) | tr ' ' '\n' | entr -s 'make && ./1d_fd_simple_shear'
@@ -27,3 +40,4 @@ clean:
 	$(RM) *.o
 	$(RM) 1d_fd_simple_shear
 	$(RM) 1d_fd_simple_shear.png
+	$(RM) 1d_fd_simple_shear_rheology.png