162 lines
3.1 KiB
Markdown
Executable File
162 lines
3.1 KiB
Markdown
Executable File
|
|
* create a 3-D hexahedral mesh object and name it 3dmesh
|
|
|
|
cmo/create/3dmesh///hex/
|
|
|
|
* create a cylinder inside a cylinder
|
|
|
|
* inner cylinder will be replaced by tet grid
|
|
|
|
surface/cyl1/reflect/cylinder/0.0,0.0,0.0/0.0,0.0,10.0/5.0/
|
|
|
|
surface/top/reflect/plane/0.,0.,10./1.,0.,10./1.,1.,10. /
|
|
|
|
surface/bot/reflect/plane/0.,0.,0./1.,0.,0./1.,1.,0. /
|
|
|
|
surface/cyl2/intrface/cylinder/0.0,0.0,1.0/0.0,0.0,9.0/4.0/
|
|
|
|
surface/cutpl1/intrface/plane/0.,0.,1./1.,0.,1./1.,1.,1./
|
|
|
|
surface/cutpl2/intrface/plane/0.,0.,9./1.,0.,9./1.,1.,9./
|
|
|
|
*define geometric regions
|
|
|
|
region/inner/ lt cyl2 and lt cutpl2 and gt cutpl1 /
|
|
|
|
region/outer/ le cyl1 and le top and ge bot and &
|
|
|
|
( ge cyl2 or ge cutpl2 or le cutpl1 )
|
|
|
|
* define material regions
|
|
|
|
mregion/inner/ lt cyl2 and lt cutpl2 and gt cutpl1 /
|
|
|
|
mregion/outer/ le cyl1 and le top and ge bot and &
|
|
|
|
( gt cyl2 or gt cutpl2 or lt cutpl1 )
|
|
|
|
* distribute points in cylindrical arrangement
|
|
|
|
rz/rtz/11,61,21/0.,0.,0./5.,360.,10./1,1,1/
|
|
|
|
* assign material colors to the points
|
|
|
|
* identify points that are on material interfaces
|
|
|
|
* identify constrained points
|
|
|
|
* connect hex mesh
|
|
|
|
rzbrick/xyz/61,11,21/1,0,0/connect/
|
|
|
|
zq/itp/1,0,0/0
|
|
|
|
zq/imt/1,0,0/0
|
|
|
|
zq/icr/1,0,0/0
|
|
|
|
setpts
|
|
|
|
dump/x3d/x3d1
|
|
|
|
settets
|
|
|
|
zq/itp/100,130,1
|
|
|
|
dump/gmv/gmv1
|
|
|
|
* remove inner region
|
|
|
|
rmmat/inner
|
|
|
|
rmpoint/compress
|
|
|
|
geniee
|
|
|
|
*zq/xic/1,0,0
|
|
|
|
*zq/itp/1,0,0
|
|
|
|
dump/gmv/gmv2
|
|
|
|
* grid tet part
|
|
|
|
cmo/create/tetgrid///tet
|
|
|
|
rz/rtz/9,61,17/0.,0.,1./1.,360.,9./1,1,1/
|
|
|
|
filter/1,0,0/
|
|
|
|
search
|
|
|
|
zq/imt/1,0,0/1
|
|
|
|
zq/xic/1,0,0
|
|
|
|
zq/itp/1,0,0
|
|
|
|
* add meshes
|
|
|
|
dump/x3d/x3d1
|
|
|
|
addmesh/pyramid/hybrid/tetgrid/3dmesh/
|
|
|
|
surface/cyl1/reflect/cylinder/0.0,0.0,0.0/0.0,0.0,6.0/2.0/
|
|
|
|
surface/top/reflect/plane/0.,0.,6./1.,0.,6./1.,1.,6. /
|
|
|
|
surface/bot/reflect/plane/0.,0.,0./1.,0.,0./1.,1.,0. /
|
|
|
|
surface/cyl2/intrface/cylinder/0.0,0.0,2.0/0.0,0.0,4.0/1.0/
|
|
|
|
surface/cutpl1/intrface/plane/0.,0.,2./1.,0.,2./1.,1.,2./
|
|
|
|
surface/cutpl2/intrface/plane/0.,0.,4./1.,0.,4./1.,1.,4./
|
|
|
|
*define geometric regions
|
|
|
|
region/inner/ le cyl2 and le cutpl2 and ge cutpl1 /
|
|
|
|
region/outer/ le cyl1 and le top and ge bot and &
|
|
|
|
( gt cyl2 or gt cutpl2 or lt cutpl1 )
|
|
|
|
* define material regions
|
|
|
|
mregion/inner/ lt cyl2 and lt cutpl2 and gt cutpl1 /
|
|
|
|
mregion/outer/ le cyl1 and le top and ge bot and &
|
|
|
|
( gt cyl2 or gt cutpl2 or lt cutpl1 )
|
|
|
|
filter/1,0,0/
|
|
|
|
dump/x3d/x3d2
|
|
|
|
geniee
|
|
|
|
zq/imt/1,0,0/0/
|
|
|
|
zq/itp/1,0,0/0/
|
|
|
|
zq/icr/1,0,0/0/
|
|
|
|
setpts
|
|
|
|
settets
|
|
|
|
* dump mesh to some output form
|
|
|
|
dump/gmv/gmv.hybrid/hybrid
|
|
|
|
* terminate processing
|
|
|
|
finish
|
|
|
|
|
|
|
|
[Return to LaGriT Home Page](index.md)
|
|
|
|
|
|
|