diff --git a/src/sand-wedge-data/mygeometry.hh b/src/sand-wedge-data/mygeometry.hh index 19e0553330e13e560cc6e906d75fffc56e0cf92a..fdfcdea4e3d156b9bfa192daf4fec4df33e299f0 100644 --- a/src/sand-wedge-data/mygeometry.hh +++ b/src/sand-wedge-data/mygeometry.hh @@ -85,7 +85,7 @@ namespace { double const lengthScale = reference::s; -double const depth = 0.10; +double const depth = 0.60 * lengthScale; LocalVector const A = rotate(reference::A); LocalVector const B = rotate(reference::B); diff --git a/src/sand-wedge-data/parset.cfg b/src/sand-wedge-data/parset.cfg index ab6670d0db9cb2ceb255aac3bbe3854c35a72874..45a42435e108c57c28bdd3f25cd8c707a429f37b 100644 --- a/src/sand-wedge-data/parset.cfg +++ b/src/sand-wedge-data/parset.cfg @@ -9,32 +9,32 @@ writeVTK = false finalTime = 1800 # [s] [body] -bulkModulus = 1e5 # [Pa] -poissonRatio = 0.3 # [1] 0.2 - 0.3 +bulkModulus = 0.5e5 # [Pa] +poissonRatio = 0.3 # [1] [body.elastic] density = 900 # [kg/m^3] -shearViscosity = 1e3 # [Pas] 0 -bulkViscosity = 1e3 # [Pas] 0 +shearViscosity = 1e3 # [Pas] +bulkViscosity = 1e3 # [Pas] [body.viscoelastic] density = 1000 # [kg/m^3] shearViscosity = 1e4 # [Pas] bulkViscosity = 1e4 # [Pas] [boundary.friction] -smallestDiameter= 1e-3 # [m] +smallestDiameter= 2e-3 # [m] C = 10 # [Pa] -mu0 = 0.7 # [1] +mu0 = 0.7 # [ ] V0 = 5e-5 # [m/s] -L = 2e-5 # [m] ? -initialAlpha = 0.916290731874155 # [ ] ? +L = 2.5e-5# [m] +initialAlpha = 0 # [ ] stateModel = AgeingLaw frictionModel = Truncated [boundary.friction.weakening] -a = 0.015 # [1] ? -b = 0.030 # [1] ? +a = 0.002 # [ ] +b = 0.014 # [ ] [boundary.friction.strengthening] -a = 0.030 # [1] ? -b = 0.015 # [1] ? +a = 0.025 # [ ] +b = 0.005 # [ ] [timeSteps] refinementTolerance = 1e-5 @@ -42,22 +42,22 @@ number = 100000 scheme = newmark [u0.solver] -tolerance = 1e-10 +tolerance = 1e-8 maximumIterations = 100000 verbosity = quiet [a0.solver] -tolerance = 1e-10 +tolerance = 1e-8 maximumIterations = 100000 verbosity = quiet [v.solver] -tolerance = 1e-10 +tolerance = 1e-8 maximumIterations = 100000 verbosity = quiet [v.fpi] -tolerance = 1e-10 +tolerance = 1e-8 maximumIterations = 10000 lambda = 0.5