From 73b81a548334903e6e439244c8f7cedd59754cfa Mon Sep 17 00:00:00 2001 From: Oliver Sander <oliver.sander@tu-dresden.de> Date: Fri, 2 Aug 2024 21:10:58 +0200 Subject: [PATCH] Update one parameter file to work with cosserat-continuum.cc again --- problems/cosserat-continuum-twisted-strip.parset | 13 ++++++++++--- 1 file changed, 10 insertions(+), 3 deletions(-) diff --git a/problems/cosserat-continuum-twisted-strip.parset b/problems/cosserat-continuum-twisted-strip.parset index a10d3197..84cad73c 100644 --- a/problems/cosserat-continuum-twisted-strip.parset +++ b/problems/cosserat-continuum-twisted-strip.parset @@ -2,7 +2,7 @@ # Grid parameters ############################################# -structuredGrid = true +structuredGrid = cube lower = 0 -0.005 upper = 0.1 0.005 elements = 10 1 @@ -21,7 +21,7 @@ numHomotopySteps = 24 tolerance = 1e-8 # Max number of steps of the trust region solver -maxTrustRegionSteps = 200 +maxSolverSteps = 200 solverScaling = 1 1 1 1 1 1 @@ -80,6 +80,12 @@ q = 2 # Shear correction factor kappa = 1 +# TODO: These three parameters are not actually used, +# but the current implementation requires them nevertheless. +b1 = 1 +b2 = 1 +b3 = 1 + [] ############################################# @@ -90,7 +96,8 @@ problem = twisted-strip ### Python predicate specifying all Dirichlet grid vertices # x is the vertex coordinate -dirichletVerticesPredicate = "x[0] < 0.001 or x[0] > 0.0999" +dirichletVerticesPredicate = "[x[0] < 0.001 or x[0] > 0.0999, x[0] < 0.001 or x[0] > 0.0999, x[0] < 0.001 or x[0] > 0.0999]" +dirichletRotationVerticesPredicate = "x[0] < 0.001 or x[0] > 0.0999" # Initial deformation initialDeformation = "[x[0], x[1], 0]" -- GitLab