Skip to content
Snippets Groups Projects
Commit f634ad1f authored by Lisa Julia Nebel's avatar Lisa Julia Nebel
Browse files

Adjust and cleanup film-on-substrate.parset

parent 02a76392
No related branches found
No related tags found
1 merge request!45Feature/proximal newton solver
...@@ -9,13 +9,13 @@ structuredGrid = true ...@@ -9,13 +9,13 @@ structuredGrid = true
lower = 0 0 0 lower = 0 0 0
upper = 200 100 200 upper = 200 100 200
elements = 10 5 5 elements = 4 2 2
# Number of grid levels, all elements containing surfaceshell grid vertices will get adaptively refined # Number of grid levels, all elements containing surfaceshell grid vertices will get adaptively refined
numLevels = 2 numLevels = 1
# When starting from a file, the stress-free configuration of the surfaceShell is read from a file, this file needs to match the *finest* grid level! # When starting from a file, the stress-free configuration of the surfaceShell is read from a file, this file needs to match the *finest* grid level!
startFromFile = true startFromFile = false
pathToGridDeformationFile = ./ pathToGridDeformationFile = ./
gridDeformationFile = deformation gridDeformationFile = deformation
...@@ -26,7 +26,7 @@ gridDeformation="[1.3*x[0], x[1], x[2]]" ...@@ -26,7 +26,7 @@ gridDeformation="[1.3*x[0], x[1], x[2]]"
# Boundary values # Boundary values
############################################# #############################################
problem = cantilever dirichletValues = identity-dirichlet-values
### Python predicate specifying all Dirichlet grid vertices ### Python predicate specifying all Dirichlet grid vertices
# x is the vertex coordinate # x is the vertex coordinate
...@@ -34,7 +34,7 @@ dirichletVerticesPredicate = "x[0] < 0.01" ...@@ -34,7 +34,7 @@ dirichletVerticesPredicate = "x[0] < 0.01"
### Python predicate specifying all surfaceshell grid vertices, elements conataining these vertices will get adaptively refined ### Python predicate specifying all surfaceshell grid vertices, elements conataining these vertices will get adaptively refined
# x is the vertex coordinate # x is the vertex coordinate
surfaceShellVerticesPredicate = "x[2] > 199.99" surfaceShellVerticesPredicate = "x[2] > 199.99 and x[0] > 49.99 and x[0] < 150.01"
### Python predicate specifying all Neumann grid vertices ### Python predicate specifying all Neumann grid vertices
# x is the vertex coordinate # x is the vertex coordinate
...@@ -60,7 +60,7 @@ numHomotopySteps = 1 ...@@ -60,7 +60,7 @@ numHomotopySteps = 1
tolerance = 1e-3 tolerance = 1e-3
# Max number of solver steps # Max number of solver steps
maxSolverSteps = 300 maxSolverSteps = 1
# Measure convergence # Measure convergence
instrumented = 0 instrumented = 0
...@@ -70,7 +70,7 @@ instrumented = 0 ...@@ -70,7 +70,7 @@ instrumented = 0
############################################# #############################################
# initial regularization parameter # initial regularization parameter
initialRegularization = 10000000000 initialRegularization = 1000000
############################################# #############################################
# Solver parameters specific for trust-region solver using multigrid solver # Solver parameters specific for trust-region solver using multigrid solver
...@@ -102,14 +102,11 @@ mgTolerance = 1e-5 ...@@ -102,14 +102,11 @@ mgTolerance = 1e-5
# Tolerance of the base grid solver # Tolerance of the base grid solver
baseTolerance = 1e-8 baseTolerance = 1e-8
# Measure convergence
instrumented = 0
############################ ############################
# Material parameters # Material parameters
############################ ############################
energy = stvenantkirchhoff energy = mooneyrivlin
## For the Wriggers L-shape example ## For the Wriggers L-shape example
[materialParameters] [materialParameters]
...@@ -125,13 +122,7 @@ surfaceShellParameters = surface-shell-parameters-1-3 ...@@ -125,13 +122,7 @@ surfaceShellParameters = surface-shell-parameters-1-3
mu_c = 0 mu_c = 0
# Length scale parameter # Length scale parameter
L_c = 0.2 L_c = 0.2
# Curvature exponent
q = 2
# Shear correction factor
kappa = 1
b1 = 1 b1 = 1
b2 = 1 b2 = 1
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment