mofem/html/horizontal__crack_8jou_source.html

reset


#Geometry

brick x 5 y 10 z 1

create cylinder height 1 radius 1

webcut volume 2 with plane yplane imprint merge

chop volume 1  with volume 2

chop volume 5 with volume 3

imprint all

merge all

color Volume all win2kgray


#Material properties

block 1 volume all

block 1 name "MAT_ELASTIC1"

block 1 attribute count 10

block 1 attribute index 1 1000

block 1 attribute index 2 0.3


#Create nodesets and sidesets for fracture constraints

sideset 100 curve 2 11 8 12 1 5 7 3 4 6 9 10

nodeset 101 vertex 8 5 2 3 4 7 6 1

sideset 102 surface all

sideset 102 surface 16 18 12 remove


sideset 200 surface 18

sideset 201 curve 27 29


#Boundary conditions and loads

create pressure on surface 5 magnitude -1 bottom

create pressure on surface 3 magnitude -1 bottom

create displacement on vertex 4 dof 1 dof 2 dof 3 fix

create displacement on vertex 1 dof 2 dof 3 fix

create displacement on vertex 3 dof 3 fix


#Meshing

volume all scheme tetmesh

volume all size auto factor 6

mesh volume 4

mesh volume 6

mesh volume 7