diff --git a/problems/film-on-substrate-deformation.parset b/problems/film-on-substrate-deformation.parset index f4f147e5d2ad9143fcb239c7a21d4319c5df4046..6626ae5c000bcc6505429a73fbf983f35cf94599 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 3eaf945fc5d06b37bbd7433adf63a813374489ec..6fc832fc7e0af15be0708414404a71493c8ddce1 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 0000000000000000000000000000000000000000..e540d38b4f0df2974fcc039c484792967b9c6054 --- /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 eac737a0f584c379922fdf950115dd5654daac28..0000000000000000000000000000000000000000 --- 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