From 84975a57bc81d4d3664e8531a75a3cc732d9b134 Mon Sep 17 00:00:00 2001 From: Elias Pipping <elias.pipping@fu-berlin.de> Date: Mon, 3 Sep 2012 18:50:40 +0200 Subject: [PATCH] Solve a problem that is not soft --- src/one-body-sample.parset | 21 +++++++-------------- src/one-body-sample.py | 2 +- 2 files changed, 8 insertions(+), 15 deletions(-) diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset index 1f8449c4..b20c4226 100644 --- a/src/one-body-sample.parset +++ b/src/one-body-sample.parset @@ -17,12 +17,8 @@ enable_timer = false refinements = 4 [body] -E = 2e7 -# "Most steels and rigid polymers when used within their design limits -# (before yield) exhibit values of about 0.3" -# -- http://en.wikipedia.org/wiki/Poisson's_ratio -nu = 0.3 # Needs to be strictly below 0.5; the closer we get, the - # tougher the problem. But we do have convergence! +E = 5e7 +nu = 0.3 # The close we get to 0.5, the more wiggly everything gets [solver] tolerance = 1e-11 @@ -51,31 +47,28 @@ acceptFactor = 1.0 requiredResidual = 1e-12 [boundary.friction] -normalstress = 10000000 -# "Most dry materials in combination have friction coefficient values -# between 0.3 and 0.6" -# -- http://en.wikipedia.org/wiki/Friction#Coefficient_of_friction +normalstress = 50000 # "mu_0 is the nominal coefficient of friction that has values near 0.6" # -- James H. Dieterich and Brian Kilgore: Implications of fault # constitutive properties for earthquake prediction # http://earthquake.usgs.gov/research/physics/lab/prediction.pdf -mu = 0.5 +mu = 0.6 eta = 1 model = Exponential [boundary.friction.state] # log(alpha(0)) -initial = -2.23805 +initial = 0 model = Dieterich # Ruina [boundary.friction.ruina] # "For rocks, typical values of A and B range from 0.005 to 0.015" # -- Ronaldo I. Borja and Craig D. Foster: # Continuum mathematical modeling of slip weakening in geological systems -a = 0.0075 +a = 0.010 b = 0.015 # If L is too large, reaching a steady state takes too long -L = 1e-3 +L = 1e-5 # The second-order a and b terms in [..] are very important in # determining whether sliding is stable or unstable and (in appli- diff --git a/src/one-body-sample.py b/src/one-body-sample.py index db967172..e6db78f3 100644 --- a/src/one-body-sample.py +++ b/src/one-body-sample.py @@ -13,7 +13,7 @@ class neumannCondition: class dirichletCondition: def __call__(self, x): - return 2 + return .005 # return 0 fst = 3e-1 snd = 5e-1 -- GitLab