From ddd0687b19f8910cb643566fe47f726cb0ce7e61 Mon Sep 17 00:00:00 2001
From: podlesny <podlesny@zedat.fu-berlin.de>
Date: Mon, 20 Sep 2021 15:12:32 +0200
Subject: [PATCH] read minrelativetau param
---
src/foam/foam.cc | 2 +-
src/foam/foam.cfg | 13 +++++++------
2 files changed, 8 insertions(+), 7 deletions(-)
diff --git a/src/foam/foam.cc b/src/foam/foam.cc
index 75a93bf5..22ac0604 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 572b7afc..56c4e935 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
--
GitLab