cosserat-continuum.parset

# Number of grid levels
numLevels = 3

# Number of homotopy steps for the Dirichlet boundary conditions
numHomotopySteps = 60

# Tolerance of the trust region solver
tolerance = 1e-8

# Max number of steps of the trust region solver
maxTrustRegionSteps = 200

# 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

############################
#   Problem specifications
############################


# Parameters for the shearing/wrinkling example from Taylor, Bertoldi, Steigmann
# 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
#mu_c = 3.8462e+05

# 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
#
## Lame parameters
## corresponds to E = 71240 N/mm^2, nu=0.31
#mu = 2.7191e+04
#lambda = 4.4364e+04
#
## Cosserat couple modulus
#mu_c = 0
##mu_c = 2.7191e+04
#
## Length scale parameter
#L_c = 10
#
## Curvature exponent
#q = 2
#
## Shear correction factor
#kappa = 1
#
#[]

###  Neumann values, if needed
neumannValues =  -1e3 0 0



# Wriggers L-shape
path = /home/sander/data/shells/wriggers_L_shape/
gridFile = wriggers-L-shape.msh