From 3c22e19480b3c81a5b1e757060aa3e6c884d7223 Mon Sep 17 00:00:00 2001 From: Oliver Sander <oliver.sander@tu-dresden.de> Date: Thu, 21 Jan 2021 14:15:27 +0100 Subject: [PATCH] Rename problem file 'wriggers-l-shape' to 'identity-deformation' MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Because that is what it is: It implements the identity deformation for three-dimensional continua. The name 'wriggers-l-shape' refers to a simulation that appeared in O. Sander, P. Neff, and M. Bîrsan: Numerical Treatment of a Geometrically Nonlinear Planar Cosserat Shell Model, (2016) It belongs to the dune-gfe module, and only got into dune-elasticity by accident. --- problems/film-on-substrate-deformation.parset | 2 +- .../finite-strain-elasticity-bending.parset | 5 ++--- problems/identity-deformation.py | 8 ++++++++ problems/wriggers-l-shape-dirichlet-values.py | 18 ------------------ 4 files changed, 11 insertions(+), 22 deletions(-) create mode 100644 problems/identity-deformation.py delete mode 100644 problems/wriggers-l-shape-dirichlet-values.py diff --git a/problems/film-on-substrate-deformation.parset b/problems/film-on-substrate-deformation.parset index f4f147e..6626ae5 100644 --- a/problems/film-on-substrate-deformation.parset +++ b/problems/film-on-substrate-deformation.parset @@ -115,7 +115,7 @@ mooneyrivlin_energy = square # log, square or ciarlet; different ways to compute # Boundary values ############################################# -problem = wriggers-l-shape +problem = identity-deformation ### Python predicate specifying all Dirichlet grid vertices # x is the vertex coordinate diff --git a/problems/finite-strain-elasticity-bending.parset b/problems/finite-strain-elasticity-bending.parset index 3eaf945..6fc832f 100644 --- a/problems/finite-strain-elasticity-bending.parset +++ b/problems/finite-strain-elasticity-bending.parset @@ -56,7 +56,6 @@ baseTolerance = 1e-8 energy = stvenantkirchhoff -## For the Wriggers L-shape example [materialParameters] # Lame parameters @@ -79,7 +78,7 @@ mooneyrivlin_12 = 0 mooneyrivlin_03 = 0 mooneyrivlin_k = 1e+6 # volume-preserving parameter -mooneyrivlin_energy = log # log, square or ciarlet; different ways to compute the Mooney-Rivlin-Energy +mooneyrivlin_energy = log # log, square or Ciarlet; different ways to compute the Mooney-Rivlin-Energy [] @@ -87,7 +86,7 @@ mooneyrivlin_energy = log # log, square or ciarlet; different ways to compute th # Boundary values ############################################# -problem = wriggers-l-shape +problem = identity-deformation ### Python predicate specifying all Dirichlet grid vertices # x is the vertex coordinate diff --git a/problems/identity-deformation.py b/problems/identity-deformation.py new file mode 100644 index 0000000..e540d38 --- /dev/null +++ b/problems/identity-deformation.py @@ -0,0 +1,8 @@ +class DirichletValues: + def __init__(self, homotopyParameter): + self.homotopyParameter = homotopyParameter + + # The identity deformation of 3d classical materials + def dirichletValues(self, x): + # Clamp the shape in its reference configuration + return [x[0], x[1], x[2]] diff --git a/problems/wriggers-l-shape-dirichlet-values.py b/problems/wriggers-l-shape-dirichlet-values.py deleted file mode 100644 index eac737a..0000000 --- a/problems/wriggers-l-shape-dirichlet-values.py +++ /dev/null @@ -1,18 +0,0 @@ -class DirichletValues: - def __init__(self, homotopyParameter): - self.homotopyParameter = homotopyParameter - - # Deformation of 3d classical materials - def dirichletValues(self, x): - # Clamp the L-shape in its reference configuration - return [x[0], x[1], x[2]] - - # Deformation of Cosserat shells - def deformation(self, x): - # Clamp the L-shape in its reference configuration - return [x[0], x[1], 0] - - # Orientation of Cosserat materials - def orientation(self, x): - rotation = [[1,0,0], [0, 1, 0], [0, 0, 1]] - return rotation -- GitLab