From bcd82778083a5aa13751a35057a6fbc4d03988b8 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 27 Jun 2014 10:02:14 +0000
Subject: [PATCH] Cleanup and remove all parameters not belonging to the
 Wong/Pellegrino example

And add a little quirk: I don't know how to do 'import math' in dunepython yet.
Therefore I have to replace the cosine in the initial configuration by a
series expansion...

[[Imported from SVN: r9804]]
---
 ...continuum-wong-pellegrino-wrinkling.parset | 105 +++++-------------
 1 file changed, 29 insertions(+), 76 deletions(-)

diff --git a/cosserat-continuum-wong-pellegrino-wrinkling.parset b/cosserat-continuum-wong-pellegrino-wrinkling.parset
index fba98912..6e911f02 100644
--- a/cosserat-continuum-wong-pellegrino-wrinkling.parset
+++ b/cosserat-continuum-wong-pellegrino-wrinkling.parset
@@ -1,6 +1,19 @@
+#############################################
+#  Grid parameters
+#############################################
+
+structuredGrid = true
+lower = 0 0
+upper = 0.38 0.128
+elements = 15 5
+
 # Number of grid levels
 numLevels = 1
 
+#############################################
+#  Solver parameters
+#############################################
+
 # Number of homotopy steps for the Dirichlet boundary conditions
 numHomotopySteps = 1
 
@@ -38,90 +51,27 @@ baseTolerance = 1e-8
 instrumented = 0
 
 ############################
-#   Problem specifications
+#   Material parameters
 ############################
 
 
-# Parameters for the shearing/wrinkling example from Taylor, Bertoldi, Steigmann
+# Parameters for the shearing/wrinkling example from Wong/Pellegrino 2006
 # We use 'meters' as the length unit
-#[materialParameters]
-
-# shell thickness
-#thickness = 2.5e-5
-
-# Lame parameters
-# # corresponds to E = 3.5GPa, nu=0.31
-#mu = 5.6452e+09
-#lambda = 2.1796e+09
-
-# Cosserat couple modulus
-#mu_c = 0
-
-# Length scale parameter
-# h * sqrt(1/12)
-#L_c = 7.2169e-06
-
-# Curvature exponent
-#q = 2
-
-# Shear correction factor
-#kappa = 1
-
-[]
-
-#structuredGrid = true
-#lower = 0 0
-#upper = 0.38 0.128
-#elements = 15 5
-
-
-## Parameters for the twisted 10x1 strip
-#[materialParameters]
-
-# shell thickness
-#thickness = 0.002
-
-# Lame parameters
-# # corresponds to E = 3.5GPa, nu=0.31
-#mu = 5.6452e+09
-#lambda = 2.1796e+09
-
-# Cosserat couple modulus
-mu_c = 0
-
-# Length scale parameter
-#L_c = 0.002
-
-# Curvature exponent
-#q = 2
-
-# Shear correction factor
-#kappa = 1
-
-#[]
-
-#structuredGrid = true
-#lower = 0 0
-#upper = 0.1 0.01
-#elements = 10 1
-
-## For the Wriggers L-shape example
 [materialParameters]
 
 # shell thickness
-thickness = 0.6
+thickness = 2.5e-5
 
 # Lame parameters
-# corresponds to E = 71240 N/mm^2, nu=0.31
-# However, we use units N/m^2
-mu = 2.7191e+10
-lambda = 4.4364e+10
+# corresponds to E = 3.5GPa, nu=0.31
+mu = 5.6452e+09
+lambda = 2.1796e+09
 
 # Cosserat couple modulus
 mu_c = 0
 
 # Length scale parameter
-L_c = 0.6
+L_c = 2.5e-5
 
 # Curvature exponent
 q = 2
@@ -131,12 +81,15 @@ kappa = 1
 
 []
 
-###  Neumann values, if needed
-neumannValues =  -1e3 0 0
+#############################################
+#  Boundary values
+#############################################
 
-structuredGrid = false
-path = /home/sander/data/shells/wriggers_L_shape/
-gridFile = wriggers-L-shape_99.msh
+###  Python predicate specifying all Dirichlet grid vertices
+# x is the vertex coordinate
+dirichletVerticesPredicate = "x[1] < 0.0001 or x[1] > 0.128 - 0.0001"
 
 # Initial deformation
-initialDeformation = "[x[0], x[1], 0.002*cos(1e4*x[0])]"
+#initialDeformation = "[x[0], x[1], 0.002*cos(1e4*x[0])]"
+# Series expansion of the cosine, because I haven't figured out how to use math stuff in dunepython yet
+initialDeformation = "[x[0], x[1], 1e-15*(1- 0.5*(1e4*x[0])*(1e4*x[0]) + (1.0/24)*(1e4*x[0])*(1e4*x[0])*(1e2*x[0])*(1e4*x[0]))]"
-- 
GitLab