Newer
Older
#############################################
# Grid parameters
#############################################
structuredGrid = false
path = /home/sander/data/shells/wriggers_L_shape/
gridFile = wriggers-L-shape_99.msh

Oliver Sander
committed
# Number of grid levels
numLevels = 1
#############################################
# Solver parameters
#############################################

Oliver Sander
committed
# Number of homotopy steps for the Dirichlet boundary conditions
numHomotopySteps = 1
# Tolerance of the trust region solver

Oliver Sander
committed
# Max number of steps of the trust region solver

Oliver Sander
committed
# Initial trust-region radius
initialTrustRegionRadius = 0.1
# Number of multigrid iterations per trust-region step
numIt = 200
# Number of presmoothing steps
nu1 = 3
# Number of postsmoothing steps
nu2 = 3
# Number of coarse grid corrections
mu = 1
# Number of base solver iterations
baseIt = 100
# Tolerance of the multigrid solver
mgTolerance = 1e-7
# Tolerance of the base grid solver
baseTolerance = 1e-8
# Measure convergence
instrumented = 0
############################
# Material parameters

Oliver Sander
committed
############################
## For the Wriggers L-shape example
[materialParameters]
# shell thickness
thickness = 0.6e-3

Oliver Sander
committed
# 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
# Cosserat couple modulus
mu_c = 0
# Length scale parameter

Oliver Sander
committed
# Curvature exponent
q = 2
# Shear correction factor

Oliver Sander
committed
[]
#############################################
# Boundary values
#############################################

Oliver Sander
committed
### Python predicate specifying all Dirichlet grid vertices
# x is the vertex coordinate
dirichletVerticesPredicate = "x[0] < 0.001"
### Python predicate specifying all Dirichlet grid vertices
# x is the vertex coordinate
neumannVerticesPredicate = "x[1] < -0.239"

Oliver Sander
committed
### Neumann values, if needed

Oliver Sander
committed
#initialDeformation = "[x[0], x[1], -0.3*x[0]*(0.24+x[1])]"
initialDeformation = "[x[0], x[1], 0 if (x[0] < 0.225 or x[1] <-0.015) else 4*(x[0]-0.225)*(x[1]+0.015)]"