mofem/html/radwanska__benchmark2_8jou_source.html

reset


create Cylinder height 10 radius 83.35

create Cylinder height 10 radius 83.25

brick x 200 y 200 z 20


move Volume 3 x 0 y -100 z 0 include_merged


subtract volume 2 from volume 1

Volume 3 copy rotate +86.561 about z

Volume 3 copy rotate -86.561 about z


delete volume 3

subtract volume 4 from volume 1

subtract volume 5 from volume 1


move Volume 1 location 0 0 0 include_merged

brick x 10 y 10 z 10

move Volume 6 x 5 y 0 z 5 include_merged

intersect volume all


block 1 volume 1

block 1 name "MAT_ELASTIC_1"

block 1 attribute count 2

block 1 attribute index 1 3.4

block 1 attribute index 2 0.2


nodeset 1 vertex 57


create Displacement  on surface 49  dof 1 fix 0

create Displacement  on surface 50  dof 3 fix 0

create Displacement  on curve 78  dof 2 fix 0

create Displacement  on curve 77  dof 2 fix 0


create pressure  on surface 48  magnitude 1

create pressure  on surface 46  magnitude 1


create force  on vertex 57  vector 0 0.5e-6 0 0 0 0

nodeset 5 name 'Imperfection'


volume all scheme Tetmesh

volume all size auto factor 10

mesh volume all


block 2 tet all

block 2 element type tetra10


save as "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/radwanska_benchmark2.cub" overwrite