[CIG-SHORT] Cubit tet mesh not being read into pylith - journal files
Brad Aagaard
baagaard at usgs.gov
Thu Oct 22 15:54:08 PDT 2009
Oliver-
You never cut out the small box from your big box. You need to add the
following 2 lines to your geometry.jou file after deleting surface 20.
subtract volume 2 4 from volume 6 keep
delete volume 6
Your revised journal file should look like what I have below.
NOTE TO ALL USERS: There are many different ways to create the solid
model geometry in LaGriT and CUBIT. Regardless of how you create your
geometry, you need to make sure the volumes do *not* overlap.
Brad
# ----------------------------------------------------------------------
# Create block
# ----------------------------------------------------------------------
# Block is 1000 km x 1000 km x 100 km
# -500 km <= x <= 500 km
# -500 km <= y <= 500 km
# -100 km <= z <= 0 km
reset
brick x 1000000 y 1000000 z 100000
volume 1 move x 0 y 0 z -50000
# ----------------------------------------------------------------------
# Create fault surface and subdivide fault volume
# ----------------------------------------------------------------------
create vertex 0 -65e3 -20e3
create vertex 0 65e3 -20e3
create vertex 0 65e3 -40e3
create vertex 0 -65e3 -40e3
create surface vertex 9 to 12
sweep surface 7 vector 1 0 0 distance 100e3
sweep surface 12 vector -1 0 0 distance 100e3
create surface vertex 9 to 12
webcut volume 2 with plane surface 17
delete vertex 9 to 12
delete surface 17
# ----------------------------------------------------------------------
# Create material surface and subdivide material volume
# ----------------------------------------------------------------------
create planar surface with plane zplane offset -20000
webcut volume 1 with plane surface 20
delete surface 20
subtract volume 2 4 from volume 6 keep
delete volume 6
# ----------------------------------------------------------------------
# Name volumes and surfaces
# ----------------------------------------------------------------------
volume 1 name "elastic"
volume 2 4 6 name "visco"
surface 18 name "fault_surface"
# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
imprint all with volume all
merge all
More information about the CIG-SHORT
mailing list