2set duplicate block elements on
3brick x 0.0008 y 0.0016 z 0.00015
4brick x 0.0008 y 0.0016 z 0.00015
5rotate volume 2 angle 15 about z include_merged
6move Volume 2 y 0.0015 include_merged
7move Volume 2 x -0.0004 include_merged
8brick x 0.0008 y 0.0016 z 0.00015
9rotate volume 3 angle 345 about z include_merged
10move Volume 3 y 0.0015 include_merged
11move Volume 3 x 0.0004 include_merged
12subtract volume 2 from volume 1
13subtract volume 3 from volume 1
14brick x 0.0008 y 0.0016 z 0.00015
15rotate volume 4 angle 345 about z include_merged
16move Volume 4 y -0.0015 include_merged
17move Volume 4 x -0.0004 include_merged
18brick x 0.0008 y 0.0016 z 0.00015
19rotate volume 5 angle 15 about z include_merged
20move Volume 5 y -0.0015 include_merged
21move Volume 5 x 0.0004 include_merged
22subtract volume 4 from volume 1
23subtract volume 5 from volume 1
24brick x 0.001 y 0.0016 z 0.00015
25move Volume 5 y -0.0015 include_merged
28move Volume 1 y 0.0011342 include_merged
29move Volume 7 y -0.0011342 include_merged
30subtract volume 1 from volume 6
31subtract volume 7 from volume 6
32brick x 0.001 y 0.0016 z 0.00015
33rotate volume 9 angle 15 about z include_merged
34move Volume 9 y 0.00112 include_merged
35brick x 0.001 y 0.0016 z 0.00015
36rotate volume 10 angle 345 about z include_merged
37move Volume 10 y -0.00112 include_merged
38subtract volume 9 from volume 8
39subtract volume 10 from volume 8
40brick x 0.001 y 0.0016 z 0.00015
41move Volume 11 x -0.0005 include_merged
43rotate volume 12 angle 180 about z include_merged
44move Volume 8 x 0.0001 include_merged
45move Volume 12 x -0.0001 include_merged
47move Volume 13 x 0.001 include_merged
48subtract volume 11 from volume 8
49subtract volume 13 from volume 12
50move Volume 12 x -0.00005 include_merged
51move Volume 8 x 0.00005 include_merged
52subtract volume 12 from volume 6
53subtract volume 8 from volume 6
54brick x 0.001 y 0.0016 z 0.00015
56move Volume 14 x -0.00095 include_merged
57move Volume 15 x 0.00095 include_merged
58subtract volume 14 from volume 6
59subtract volume 15 from volume 6
60brick x 0.001 y 0.0016 z 0.00015
62move Volume 16 y 0.001527405 include_merged
63move Volume 17 y -0.001527405 include_merged
64subtract volume 16 from volume 6
65subtract volume 17 from volume 6
67# Set block with volume with gel material
71volume all scheme Tetmesh
72volume all size auto factor 4
74# Set Dirichelt boundary condition for solvent concentration
75block 3 surface 81 84 83 82 138 137 139 75 74 76 73 131 130 132 182 153 184 174 146 177 209 206 201 214
76block 3 name 'CHEMICAL_LOAD_1'
77block 3 attribute count 1
78block 3 attribute index 1 3.98950543760577e06
79# Set solvent flux, i.e.
80block 5 surface 179 180 172 171 173 186 178 181 185 183 210 204 202 212 208 207 199 200 213 211
81block 5 name 'FLUX_CHEMICAL_LOAD_2'
82block 5 attribute count 1
83block 5 attribute index 1 0
84# Kinematic boundary condition for mechanical field
85create displacement on surface 208 207 199 200 dof 2 fix 0
86# Kinematic boundary condition for mechanical field
87create displacement on surface 210 204 202 212 dof 1 fix 0
88# Kinematic boundary condition for mechanical field
89create displacement on surface 211 dof 3 fix 0
90# Set block 4 and set 10 node tetrahedrons in that block
92block 4 element type TETRA10
