[CIG-SHORT] fault meshing, SCEC mesh file

[Birendra jha]

Thank you Brad and Charles. Right now, I am trying Lagrit instead of Trelis, given its flexibility in honoring interior surfaces (faults) of different sizes. I am using the mesh_tet4_1000m.lagrit file from pylith. I understand the basic commands after reading lagrit.lanl.gov/docs/commands. 

Can you please tell me how I can import a file in Lagrit to create a surface? Here is what I do:

1. Import fault.facet file in Trelis with make_elements on so as to get the triangles
2. Create a block from the surface 
3. Export as Abaqus inp file, fault.inp
4. In the lagrit file, use
cmo/create/mo_fault
surface/surf_fault/intrface/sheet/mo_fault/< fault.inp

I get the error Sheet CMO has the wrong dimensionality: mo_fault

Thanks
Birendra 

--------------------------------------------
On Wed, 8/26/15, Charles Williams <willic3 at gmail.com> wrote:

 Subject: Re: [CIG-SHORT] fault meshing, SCEC mesh file
 To: cig-short at geodynamics.org
 Date: Wednesday, August 26, 2015, 3:23 PM
 
 There are two methods you
 can use.  The first one involves creating an interior
 volume.  You first cut the domain volume with the interior
 volume (chop command), then cut the interior volume with the
 fault surface.  The other method is to just cut the domain
 volume with an extended fault surface, but only use a subset
 of the fault nodes to define the fault nodeset.  I hope
 that helps.
 
 Cheers,
 Charles
 
 
 > On 27/08/2015, at 4:11 am, Brad Aagaard
 <baagaard at usgs.gov>
 wrote:
 > 
 > On
 08/26/2015 04:52 AM, Birendra jha wrote:
 >> Hi
 >> 
 >> I want to mesh a domain with a fault
 that does not cut it. Similarly
 >> to
 Fig. 9 in Aagaard et al, JGR 2013. How do I honor the
 fault
 >> surface in Trelis without
 using webcut? Are imprint and merge
 >>
 commands sufficient?
 >> 
 >> I see that Fig. 10 in the above paper
 does this. So I am looking for
 >> the
 mesh file for SCEC benchmark TPV13. I am interested in a 3D
 tet
 >> mesh. If there is a mesh file
 for that, it will help me a lot.
 > 
 > Figure 10 is just a diagram for the
 simulation. It does not show the actual geometry used with
 CUBIT. We are unaware of any methods for creating interior
 surfaces in CUBIT that will be honored by the mesh. If
 someone does know how to do this, please let us know!
 > 
 > In the geometry.jou
 file (https://github.com/geodynamics/pylith_benchmarks/blob/master/dynamic/scecdynrup/tpv210/geometry.jou)
 we create an extended fault surface that spans across the
 volume and use webcut to create the fault surface. We do
 scribe lines on that surface so that we can create a nodeset
 confined to the portion of the surface corresponding to the
 fault and edges of cells line up with discontinuities in the
 simulation parameters.
 > 
 > WARNING: The journal files may or may not
 work exactly as is for your version of CUBIT/Trelis. The ids
 sometimes need to be updated as a result of how entities are
 split by CUBIT. Using Idless journal files does remedy this
 problem. We have tried to migrate the examples bundled with
 PyLith to be independent of the CUBIT/Trelis version, but
 have not migrated all of the benchmark journal files.
 > 
 > Regards,
 > Brad
 > 
 >
 _______________________________________________
 > CIG-SHORT mailing list
 > CIG-SHORT at geodynamics.org
 > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
 
 Charles A. Williams
 Scientist
 GNS Science
 1 Fairway Drive, Avalon
 PO Box
 30368
 Lower Hutt  5040
 New
 Zealand
 ph (office): 0064-4570-4566
 fax (office): 0064-4570-4600
 C.Williams at gns.cri.nz
 
 _______________________________________________
 CIG-SHORT mailing list
 CIG-SHORT at geodynamics.org
 http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short