Solve a problem that is not soft

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-

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