From 0bfdb4d6a24aa20e5f5e7274d27fe2c836b3b519 Mon Sep 17 00:00:00 2001
From: Romain Boman <romain.boman@gmail.com>
Date: Mon, 2 Dec 2024 12:44:38 +0100
Subject: [PATCH] a narch (nb of archives of metafor) as a parameter
---
models/Panthera_pardus/Panthera_pardus_RC_30.py | 1 +
models/bonemodel.py | 5 +++--
2 files changed, 4 insertions(+), 2 deletions(-)
diff --git a/models/Panthera_pardus/Panthera_pardus_RC_30.py b/models/Panthera_pardus/Panthera_pardus_RC_30.py
index 0726ed9..ede6e27 100644
--- a/models/Panthera_pardus/Panthera_pardus_RC_30.py
+++ b/models/Panthera_pardus/Panthera_pardus_RC_30.py
@@ -70,6 +70,7 @@ def parms(d={}):
p['Poisson'] = 0.3 # [-]
# p['use_gmshOld'] = True
+ p['narch'] = 2 # number of archives (Metafor)
p.update(d)
return p
diff --git a/models/bonemodel.py b/models/bonemodel.py
index 6193e0b..b95a1bb 100644
--- a/models/bonemodel.py
+++ b/models/bonemodel.py
@@ -81,7 +81,8 @@ def parms(d={}):
# numerical parameters
p['tolNR'] = 1e-6 # [-] equilibrium tolerance
p['dt0'] = 1.0 # [s] time step size
-
+ p['narch'] = 1 # number of archives
+
# gmsh toolbox
p['use_gmshOld'] = False # use old gmsh interface
@@ -232,7 +233,7 @@ def getMetafor(p={}):
tsm = metafor.getTimeStepManager()
tsm.setInitialTime(0.0, p['dt0'])
- tsm.setNextTime(1.0, 1, 1.0)
+ tsm.setNextTime(1.0, p['narch'], 1.0)
mim = metafor.getMechanicalIterationManager()
mim.setResidualTolerance(p['tolNR'])
--
GitLab