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