RZAMR uses an octree type refinement applied to an existing hexahedral mesh to all nodes in a specified region. No additional elements are kept, the intention is that the resulting node distribution will be passed to connect to generate a tetrahedra mesh.
rzamr/region_name/number_of_levels
region_name is the name of region to refine. If blank, all regions will be refined. An element will be refined if any node of
the element is in the specified region.
number_of_levels is the number of times the refinement will be performed. After each level, the code will determine which of
the new nodes are in the specified region and will refine the associated elements. Default is 1.
rzamr
refine the entire mesh
rzamr /r1/3
refine elements with nodes in the region r1 three times
Create a hex mesh, refine with rzamr, connect into a tet mesh.
cmo/create/cmo///hex
* define geometry
surface/inside/reflect/box/0,0,0/1,1,1
surface/diag/intrface/plane/0,0,0/1,0,1/1,1,1
region/lin/ le inside and ge diag /
region/rin/ le inside and lt diag /
mregion/mlin/ le in=side and gt diag /
mregion/mrin/ le inside and lt diag /
* distribute nodes
quadxyz/2,2,2/0.,0.,0./1.,0.,0./1.,1.,0./0.,1.,0./ &
0.,0.,1./1.,0.,1./1.,1.,1./0.,1.,1./
* set node types and materials
setpts
* connect up the hex mesh
rzbrick/xyz/2,2,2/1,0,0/connect
* refine the hex mesh
rzamr/lin/1
rzamr/rin/3
* create the tet mesh
cmo/create/cmot///tet
* define geometry again for tet mesh
surface/inside/reflect/box/0,0,0/1,1,1
surface/diag/intrface/plane/0,0,0/1,0,1/1,1,1
region/lin/ le inside and ge diag /
region/rin/ le inside and lt diag /
mregion/mlin/ le inside and gt diag /
mregion/mrin/ le inside and lt diag /
* copy in the nodes from the hex mesh to the tet mesh
copypts/cmot/cmo
cmo/select/cmot
cmo/release/cmo
* reset node types and materials so that setpts will use
* geometry to figure out the correct values
cmo/setatt/cmot/itp/1,0,0/0
cmo/setatt/cmot/imt/1,0,0/0
* set node types and materials
setpts
* connect up the tet mesh
connect
* set element materials
* and create parent/child nodes on interfaces
settets
dump/gmv/gmvtet
finish