cylinder8.jou

#######                 #####  
#How to run this file       #
#playback '~/cylinder8.jou' #
######                 #####
reset
#1) Create a Cylinder: volume 1
create Cylinder height 10 radius 1 
#2) Create another cylinder : volume 2
create Cylinder height 10 radius 0.5 
#3) Subtract smaller radius cylinder from the larger: volume 2 from volume 1
subtract volume 2  from volume 1
#4) Make the entire holed-cylinder into a block.
block 1 volume 1
#5) Choose element type hex8 for the block created
block 1 element type hex8 
#6) Refine mesh size if needed
#    factor:  7    6    5    4.5    4      3.1     2.5
#  elem num: 110  368  1176  1635  4180  12400  17670
volume 1 size auto factor 2
#7) Mesh the volume created
mesh volume 1
# auto factor 3.5 + numsplit 1 = 54K element
# auto factor 3   + numsplit 1 = 99K element
# auto factor 2   + numsplit 1 = 206K element
# auto factor 2.7 + numsplit 2 = 399K element
# auto factor 1   + numsplit 1 = 250K element
# auto factor 2.5 + numsplit 2 = 477K element
#7.1) Mesh refinement
refine volume 1 numsplit 1 bias 1.0 depth 1 smooth
#Choose nodesets and sidesets from Boundary conditions
#8) Create a sidesets. Call it Dirch for Dirichlet (inner wall of the cylinder)
sideset 999 surface 9 #left
sideset 998 surface 8 #right
sideset 997 surface 1  #outer
sideset 996 surface 7  #inner
#nodeset 999 surface 9 #left
#nodeset 998 surface 8 #right
#nodeset 997 surface 1  #outer
#nodeset 996 surface 7  #inner
#9) Create a nodeset.
#sideset 1 surface 1
#10) Turn on element and Node maps on
#set idmap on
#print(cubit.get_quad_count())
#11) export the mesh
export mesh "cylinder8_e_4ss_us.exo" dimension 3 overwrite