mofem/html/brick__of__gel_8jou_source.html

reset

set duplicate block elements on


brick x 10 y 1 z 1

webcut volume 1 with plane yplane imprint merge


# Sete block with volume with gel material

block 1 volume 1

block 1 name 'GEL_1'

block 2 volume 2

block 2 name 'GEL_2'


# Set Dirichelt boundary condition for solvent concentration

block 6 surface 11 13

block 6 name 'CHEMICAL_LOAD_1'

block 6 attribute count 1

block 6 attribute index 1 1


# Set solvent flux, i.e. Neumann boundary condition

block 5 surface 8 15 5 14 10 3

block 5 name 'FLUX_CHEMICAL_LOAD'

block 5 attribute count 1

block 5 attribute index 1 0


# Kinematic boundary condition for mechanical field

create displacement on surface 9 16  dof 1 fix 0

create displacement on vertex 7 dof 1 dof 2 dof 3 fix 0

create displacement on vertex 4 dof 2 fix 0

create displacement on vertex 8 dof 3 fix 0


# Static boundary condition, i.e. pressure

create pressure  on surface 6 magnitude 1


# Make a mesh

volume all scheme Tetmesh

volume all size auto factor 8

mesh volume all


# Set block 3 and set 10 node tetrahedrons in that block

block 3 tet all

block 3 element type TETRA10