New problem (with velocity weakening)

src/one-body-sample.parset

@@ -74,16 +74,16 @@ model = Ruina
 [boundary.friction.state]
 evolve = true
 # log(alpha(0))
-initial = 0
+initial = -2.42037
 
 [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.015
-b = 0.005
+a = 0.0075
+b = 0.015
 # If L is too large, reaching a steady state takes too long
-L = 1e-4
+L = 1e-3
 
 # 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

@@ -14,9 +14,9 @@ class neumannCondition:
 class dirichletCondition:
     def __call__(self, x):
         # return 0
-        fst = 5e-4
-        snd = 2e-4
-        trd = 5e-4
+        fst = 3e-4
+        snd = 5e-4
+        trd = 3e-4
         if x < 1.0/3:
             return fst
         elif x < 2.0/3: