diff --git a/src/foam/foam.cc b/src/foam/foam.cc index 75a93bf5c095d8a6aa700482af9a764e1e37ac9a..22ac060428c8aa641974462e38b7daa76343b36a 100644 --- a/src/foam/foam.cc +++ b/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, diff --git a/src/foam/foam.cfg b/src/foam/foam.cfg index 572b7afc11e7d0d89e5580a814662605ac3e0d03..56c4e93594d6bc86471156d1bf37309ce955e55e 100644 --- a/src/foam/foam.cfg +++ b/src/foam/foam.cfg @@ -1,10 +1,10 @@ # -*- 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