############################################# # Grid parameters ############################################# structuredGrid = false path = /home/sander/data/shells/wriggers_L_shape/ gridFile = wriggers-L-shape_99.msh # Number of grid levels numLevels = 1 ############################################# # Solver parameters ############################################# # Number of homotopy steps for the Dirichlet boundary conditions numHomotopySteps = 1 # Tolerance of the trust region solver tolerance = 1e-6 # Max number of steps of the trust region solver maxTrustRegionSteps = 500 # 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 ############################ ## For the Wriggers L-shape example [materialParameters] # shell thickness thickness = 0.6e-3 # 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 L_c = 0.6e-3 # Curvature exponent q = 2 # Shear correction factor kappa = 0.01 [] ############################################# # Boundary values ############################################# problem = wriggers-l-shape ### 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" ### Neumann values, if needed #neumannValues = -2e7 0 0 # Initial deformation #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)]"