diff --git a/cosserat-continuum-twisted-strip.parset b/cosserat-continuum-twisted-strip.parset index 3d748989532e8cca885c7ff941dd0c6a7940debf..2aa1856601a018e80ad7c6e9dc308da442ae3a61 100644 --- a/cosserat-continuum-twisted-strip.parset +++ b/cosserat-continuum-twisted-strip.parset @@ -1,8 +1,21 @@ +############################################# +# Grid parameters +############################################# + +structuredGrid = true +lower = 0 0 +upper = 0.1 0.01 +elements = 10 1 + # Number of grid levels numLevels = 1 +############################################# +# Solver parameters +############################################# + # Number of homotopy steps for the Dirichlet boundary conditions -numHomotopySteps = 1 +numHomotopySteps = 12 # Tolerance of the trust region solver tolerance = 1e-8 @@ -38,90 +51,26 @@ baseTolerance = 1e-8 instrumented = 0 ############################ -# Problem specifications +# Material parameters ############################ -# 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 - -# 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 +# Parameters for the twisted 0.1x0.01 strip [materialParameters] # shell thickness -thickness = 0.6 +thickness = 0.002 # 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 = 0.002 # Curvature exponent q = 2 @@ -131,9 +80,12 @@ kappa = 1 [] -### Neumann values, if needed -neumannValues = -1e3 0 0 +############################################# +# Boundary values +############################################# + +### Python predicate specifying all Dirichlet grid vertices +# x is the vertex coordinate +dirichletVerticesPredicate = "x[0] < 0.001 or x[0] > 0.0999" + -structuredGrid = false -path = /home/sander/data/shells/wriggers_L_shape/ -gridFile = wriggers-L-shape_99.msh