mesh__for__prims_8jou_source.html

reset


#GEOMETRY

create brick x 15 y 15 z 5

create brick x 15 y 15 z 0.1

#create brick x 6.4 y 6.4 z 5.1

create brick x 5.6 y 5.6 z 5.1

move surface 8 location surface 1 include_merged

align volume 3 surface 14 with surface 2

imprint volume all

merge volume all

delete volume 3

move surface 25 location 0 0 0 include_merged


#BOUNDARY CONDITIONS

sideset 101 surface 21

#sideset 202 surface 25

create displacement on surface 22 20 dof all fix 0


#BLOCKS OF MATERIALS

set duplicate block elements on

block  1 volume 1

block 1 name "MAT_ELASTIC_1"

block 1 attribute count 2

#block  2 volume 2

#block 2 name "MAT_ELASTIC_2"

#block 2 attribute count 2


#POLYMER PARAMETERS

#block 2 attribute index 1 0.25

#block 2 attribute index 2 0.49


#HYDROGEL PARAMETERS

block 1 attribute index 1 3.6e-3

block 1 attribute index 2 0.4


delete volume 2


#{mesh_factor=1}

surface 21 size auto factor {mesh_factor}

surface 21 scheme TriMesh

mesh surface 21

#surface 21 scheme copy source surface 25 source curve 25 target curve 40 nosmoothing

#mesh surface 21

volume all size auto factor 7

volume all scheme tetmesh


#volume 1 sizing function type skeleton

mesh volume all


save as "/Users/karollewandowski/moFEM/users_modules/cell_engineering/mesh_cell_force_map.cub" overwrite

list tet