read minrelativetau param

src/foam/foam.cc

@@ -513,7 +513,7 @@ int main(int argc, char *argv[]) {
         stepBase(parset, contactNetwork, totalDirichletNodes, globalFriction, frictionNodes,
                  externalForces, stateEnergyNorms);
 
-    const auto minTau = parset.get<double>("timeSteps.minRelativeTau");
+    const auto minTau = parset.get<double>("initialTime.minRelativeTau");
     AdaptiveTimeStepper<NonlinearFactory, std::decay_t<decltype(contactNetwork)>, Updaters, std::decay_t<decltype(stateEnergyNorms)>>
         timeStepper(stepBase, contactNetwork, current,
                             programState.relativeTime, programState.relativeTau, minTau,

src/foam/foam.cfg

 # -*- mode:conf -*-
 [general]
-outPath = pipping-2013-newmark-double-1e5 # output written to ./output/outPath
+outPath = pipping-2013-newmark-double-adaptive  # output written to ./output/outPath
 gravity         = 9.81     # [m/s^2]
 
 [body0]
-length          = 6.0      # [m]
+length          = 5.0      # [m]
 height          = 1.0     # [m]
 bulkModulus     = 4.12e7 #4.12e9   # [Pa] #2190
 poissonRatio    = 0.3     # [1]  #0.11
@@ -38,7 +38,7 @@ mu0             = 0.6      # [ ]
 V0              = 1e-6     # [m/s]
 L               = 1e-5  # [m]
 initialAlpha    = -10.0  # [ ]
-stateModel      = AgeingLaw
+stateModel      = AgeingLaw #AgeingLaw
 frictionModel   = Truncated #Tresca #None #Truncated #Regularised
 [boundary.friction.weakening]
 a               = 0.010    # [ ]
@@ -60,7 +60,7 @@ printProgress   = true
 restarts.first  = 0
 restarts.spacing= 50
 restarts.write  = true #true
-vtk.write       = true
+vtk.write       = false
 
 [problem]
 finalTime       = 50     # [s] #1000
@@ -69,11 +69,12 @@ bodyCount       = 2
 [initialTime]
 timeStep = 0
 relativeTime = 0.0
-relativeTau = 1e-5 # 1e-6
+relativeTau = 1e-3 # 1e-6
+minRelativeTau = 1e-6
 
 [timeSteps]
 scheme = newmark # newmark, backwardEuler
-timeSteps = 1e5
+timeSteps = 1e6
 
 [u0.solver]
 maximumIterations = 100