tImprove comments - cngf-pf - continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
HTML git clone git://src.adamsgaard.dk/cngf-pf
DIR Log
DIR Files
DIR Refs
DIR README
DIR LICENSE
---
DIR commit c1aa4dce15eba37cc67913aa232337e9205feb70
DIR parent 7a535a7e52e53c1d5967573ab4b97a062198c120
HTML Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Fri, 5 Apr 2019 08:55:34 +0200
Improve comments
Diffstat:
M 1d_fd_simple_shear.jl | 42 +++++++++++++++++++------------
1 file changed, 26 insertions(+), 16 deletions(-)
---
DIR diff --git a/1d_fd_simple_shear.jl b/1d_fd_simple_shear.jl
t@@ -12,43 +12,53 @@ G = 9.81
## Wall parameters
+### Effective normal stress exerted by top wall
# A larger normal stress deepens the deformational depth
P_wall_ = [10, 20, 40, 60, 80, 120] .* 1e3 # normal stress [Pa]
-v_x_bot = 0.0 # bottom velocity along x [m/s]
-μ_wall = 0.40 # stress ratio at top wall
+### bottom velocity along x [m/s]
+v_x_bot = 0.0
-## Nodes along z
+# stress ratio at top wall
+μ_wall = 0.40
+
+### Nodes along z
nz = 100
## Material properties
+### nonlocal amplitude [-]
# lower values of A mean that the velocity curve can have sharper curves, e.g.
# at the transition from μ ≈ μ_s
-A = 0.48 # nonlocal amplitude [-]
-
+A = 0.48 # Henann and Kamrin 2016
+### rate dependence beyond yield [-]
# lower values of b mean larger shear velocity for a given stress ratio above
# yield
-b = 0.9377 # rate dependence beyond yield [-]
+b = 0.9377 # Henann and Kamrin 2016
-μ_s = 0.3819 # static yield friction coefficient [-]
+### bulk and critical state static yield friction coefficient [-]
+#μ_s = 0.3819 # Henann and Kamrin 2016
μ_s = atan(deg2rad(22.0)) # Damsgaard et al 2013
-ϕ = 0.38 # porosity [-]
+### porosity [-]
+#ϕ = 0.38 # Henann and Kamrin 2016
+ϕ = 0.28 # Damsgaard et al 2013
+
+# representative grain size [m]
# lower values of d mean that the shear velocity curve can have sharper curves,
# e.g. at the transition from μ ≈ μ_s
-#const d = 1e-3 # representative grain size [m]
-d = 0.04 # representative grain size [m] - Damsgaard et al 2013
+#d = 1e-3 # Henann and Kamrin 2016
+d = 0.04 # Damsgaard et al 2013
-ρ_s = 2.485e3 # grain material density [kg/m^3]
+### grain material density [kg/m^3]
+#ρ_s = 2.485e3 # Henann and Kamrin 2016
+ρ_s = 2.6e3 # Damsgaard et al 2013
-## Domain size
+## Spatial settings
origo_z = 0.0
-#const L_z = 20.0*d
-L_z = 0.7
-
-## Spatial coordiantes
+#L_z = 20.0*d # Henann and Kamrin 2016
+L_z = 0.7 # Damsgaard et al 2013
z = collect(range(origo_z, L_z, length=nz))
Δz = z[2] - z[1]