Newer
Older
Oliver Sander
committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
# Number of grid levels
numLevels = 1
# Number of homotopy steps for the Dirichlet boundary conditions
numHomotopySteps = 1
# 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
# 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
# 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.6
# Curvature exponent
q = 2
# Shear correction factor
kappa = 1
[]
### Neumann values, if needed
neumannValues = -1e3 0 0
structuredGrid = false
path = /home/sander/data/shells/wriggers_L_shape/
gridFile = wriggers-L-shape_99.msh