mofem/html/layered__gel__example_8jou_source.html

reset

set duplicate block elements on

brick x 0.002 y 0.00025 z 0.00125

brick x 0.002 y 0.00025 z 0.00125

move volume 1 y 0.00025

merge surface 3 with surface 11

webcut volume all with plane zplane noimprint nomerge

delete Volume 3 4

# Set block with volume with gel material

block 1 volume 1

block 1 name 'GEL_1'

block 2 volume 2

block 2 name 'GEL_2'

# Make a mesh

volume all scheme Tetmesh

volume all size auto factor 7

mesh volume all

#surface 3 scheme trimesh

#surface 3 size auto factor 2

#mesh surface 3

#

#surface 5 9 scheme trimesh

#surface 5 9 size auto factor 2

#mesh surface 5 9

#

#volume all scheme Tetmesh

#volume all sizing function type skeleton scale 7 time_accuracy_level 2 min_size auto max_size auto max_gradient 1.2 min_depth 2 max_depth 7 min_num_layers_3d 1 min_num_layers_2d 1 min_num_layers_1d 1 max_span_ang_surf 45 max_span_ang_curve 45

#volume all sizing function type skeleton

#mesh volume all

#volume all scheme Tetmesh

#volume 1 sizing function type skeleton scale 7 time_accuracy_level 2 min_size auto max_size auto max_gradient 1.2 min_depth 2 max_depth 7 min_num_layers_3d 1 min_num_layers_2d 1 min_num_layers_1d 1 max_span_ang_surf 45 max_span_ang_curve 45

#volume 1 sizing function type skeleton

#mesh volume 1

# Set Dirichelt boundary condition for solvent concentration

block 3 surface 17 27 1 7 14 26

block 3 name 'CHEMICAL_LOAD_1'

block 3 attribute count 1

block 3 attribute index 1 1.5e06

# Set solvent flux, i.e.

block 5 surface 13 23 15 25

block 5 name 'FLUX_CHEMICAL_LOAD_2'

block 5 attribute count 1

block 5 attribute index 1 0

# Kinematic boundary condition for mechanical field

create displacement on surface 15 dof 1 fix 0

# Kinematic boundary condition for mechanical field

create displacement on surface 25 dof 1 fix 0

# Kinematic boundary condition for mechanical field

create displacement on curve 29 dof 2 fix 0

# Kinematic boundary condition for mechanical field

create displacement on surface 1 dof 3 fix 0

# Kinematic boundary condition for mechanical field

create displacement on surface 7 dof 3 fix 0

# Kinematic boundary condition for mechanical field

create displacement on surface 13 dof 3 fix 0

# Kinematic boundary condition for mechanical field

create displacement on surface 23 dof 3 fix 0

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

block 4 tet all

block 4 element type TETRA10