2set duplicate block elements on
5webcut volume 1 with plane yplane imprint merge
7# Sete block with volume with gel material
13# Set Dirichelt boundary condition for solvent concentration
15block 6 name 'CHEMICAL_LOAD_1'
16block 6 attribute count 1
17block 6 attribute index 1 1
19# Set solvent flux, i.e. Neumann boundary condition
20block 5 surface 8 15 5 14 10 3
21block 5 name 'FLUX_CHEMICAL_LOAD'
22block 5 attribute count 1
23block 5 attribute index 1 0
25# Kinematic boundary condition for mechanical field
26create displacement on surface 9 16 dof 1 fix 0
27create displacement on vertex 7 dof 1 dof 2 dof 3 fix 0
28create displacement on vertex 4 dof 2 fix 0
29create displacement on vertex 8 dof 3 fix 0
31# Static boundary condition, i.e. pressure
32create pressure on surface 6 magnitude 1
35volume all scheme Tetmesh
36volume all size auto factor 8
39# Set block 3 and set 10 node tetrahedrons in that block
41block 3 element type TETRA10