Newer
Older
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
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
# Number of grid levels
numLevels = 1
# Type of the domain decomposition algorithm
# Currently supported:
# - FixedPointIteration
# - RichardsonIteration
ddType = RichardsonIteration
# What is the preconditioner, if we are doing a Richardson iteration?
# Possible values:
# - DirichletNeumann
# - NeumannDirichlet
# - NeumannNeumann
# - RobinRobin
preconditioner = DirichletNeumann
# Tolerance of the Dirichlet-Neumann solver
ddTolerance = 1e-9
# Max number of Dirichlet-Neumann steps
maxDirichletNeumannSteps = 10
# Tolerance of the trust-region solver for the rod problem
trTolerance = 1e-12
# Max number of steps of the trust region solver
maxTrustRegionSteps = 20
# Verbosity of the trust-region solver
trVerbosity = 0
# Number of multigrid iterations per time step
numIt = 20
# 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 = 1000
# Tolerance of the multigrid solver
mgTolerance = 1e-13
# Tolerance of the base grid solver
baseTolerance = 1e-13
# Initial trust-region radius
initialTrustRegionRadius = 1
# Damping
damping = 0.6
############################
# Problem specifications
############################
#path = /home/sander/data/multicoupling/simplecoupling/
#gridFile = hexarod.grid
#dirichletNodes = hexarod.dn
#dirichletValues = hexarod.nodv
#interfaceNodes = hexarod.ifn
#numRodBaseElements = 4
## Cross-section area
#rodA = 0.0625
## Geometric moments (here: square of edge length 0.25)
#rodJ1 = 0.0013021
#rodJ2 = 0.0013021
## Material parameters
#rodE = 1e6
#rodNu = 0.3
#E = 1e6
#nu = 0.3
#rodRestEndPoint0X = 0.125
#rodRestEndPoint0Y = 0.125
#rodRestEndPoint0Z = 1
#rodRestEndPoint1X = 0.125
#rodRestEndPoint1Y = 0.125
#rodRestEndPoint1Z = 2
## Dirichlet values
#dirichletValueX = 0.125
#dirichletValueY = 0.125
#dirichletValueZ = 2.5
#dirichletAxisX = 0
#dirichletAxisY = 0
#dirichletAxisZ = 1
#dirichletAngle = 0
########################################################
#path = /home/sander/data/multicoupling/simplecoupling/
#gridFile = cube_5x5x5.grid
#dirichletNodes = cube_5x5x5.dn
#dirichletValues = cube_5x5x5.nodv
#interfaceNodes = cube_5x5x5.ifn
#numRodBaseElements = 5
# Cross-section area
#rodA = 1
# Geometric moments (here: square of edge length one)
#rod J1 = 0.0833333
#rod J2 = 0.0833333
# Material parameters
#rodE = 2.5e5
#rodNu = 0.3
## Dirichlet values
#dirichletValueX = 0.5
#dirichletValueY = 1.5
#dirichletValueZ = 10
#dirichletAxisX = 0
#dirichletAxisY = 0
#dirichletAxisZ = 1
#dirichletAngle = 90
########################################################
path = /home/sander/data/multicoupling/simplecoupling/
gridFile = cube_4x4x4.grid
dirichletNodes = cube_4x4x4.dn
dirichletValues = cube_4x4x4.nodv
interfaceNodes = cube_4x4x4.ifn
numRodBaseElements = 4
# Cross-section area
rodA = 0.0625
# Geometric moments (here: square of edge length 0.25)
#rodJ1 = 0.0013021
#rodJ2 = 0.0013021
rodJ1 = 0.00032552
rodJ2 = 0.00032552
# Material parameters
rodE = 1e6
rodNu = 0.3
E = 1e6
nu = 0.3
rodRestEndPoint0X = 0.625
rodRestEndPoint0Y = 0.625
rodRestEndPoint0Z = 1
rodRestEndPoint1X = 0.625
rodRestEndPoint1Y = 0.625
rodRestEndPoint1Z = 2
#rodRestEndPoint1X = 0.625
#rodRestEndPoint1Y = -0.082
#rodRestEndPoint1Z = 1.707
# Dirichlet values
dirichletValueX = 0.625
dirichletValueY = 0.625
dirichletValueZ = 2.5
dirichletAxisX = 0
dirichletAxisY = 0
dirichletAxisZ = 1
dirichletAngle = 0