[CIG-SHORT] [EXTERNAL] Re: defining fault_edge nodes

[Charles Williams]

Dear Niloufar,

If I’m understanding your question correctly, you need the ‘geometry/vertices’ field, which gives you the vertex coordinates, and you need the ‘topology/cells’ field, which gives you the vertices associated with each cell.  If you are filtering your results to the cell center, the cell center associated with each cell will just be the average of the vertex coordinates for the cell.  Remember that all indices are zero-based (rather than one-based as in Matlab).

Cheers,
Charles


> On 30/06/2018, at 10:29 AM, Niloufar Abolfathian <niloufar.abolfathian at gmail.com> wrote:
> 
> Hi,
> 
> I had another question regarding reading the output files. I am using matlab for reading the h5 files and further plottings.  
> For reading stress and strain related to each material I do have the stress/strain regarding each gird cell plus the topology information (only number of vertices not the geometrical values). I could not figure out how to relate the topology cells with the geometrical vertices.
> 
> Thanks,
> Niloufar
> 
> 
> On Thu, Jun 21, 2018 at 7:14 PM, Brad Aagaard <baagaard at usgs.gov <mailto:baagaard at usgs.gov>> wrote:
> On 6/21/18 5:50 PM, Niloufar Abolfathian wrote:
> Thanks, Brad. By my second question, I mean in cases where the fault is extended to the external boundaries, and we want to set boundary conditions on these surfaces(external boundaries), fault nodes should be all excluded. Is this statement right?
> 
> Yes, that is correct.
> 
> 
> 
> 
> On Thu, Jun 21, 2018 at 5:40 PM, Brad Aagaard <baagaard at usgs.gov <mailto:baagaard at usgs.gov> <mailto:baagaard at usgs.gov <mailto:baagaard at usgs.gov>>> wrote:
> 
>     On 6/21/18 12:38 PM, Niloufar Abolfathian wrote:
> 
>         Hi,
> 
>         I am modeling a 3D vertical fault in an upper crust, lying on a
>         lower crust. It is as your 3d/hex example instead of that the
>         fault is only located in the upper crust.  I set up the nodes
>         related to the fault and fault_edge. From my meshing and node
>         selection the fault should only be located in the upper crust,
>         but after running the model the fault is extended below the
>         upper crust. I guess it may be due to the selection of fault and
>         fault_edge nodes. Hope you can help me with this issue.
> 
>         Question: Should the fault_edge nodes be also included in fault
>         nodes? (I suppose the answer is yes, due to the figure 6.4 in
>         the manual.)
>         I tried excluding fault_edge nodes from the fault but then the
>         model will not be converged.
> 
> 
>     Yes. The fault edge nodes are included in the fault nodes.
> 
>     The edge nodes are any buried edges of the fault. That is, any nodes
>     on edges that are not on the external domain boundary should be
>     included in the fault edges.
> 
>         I have another question. In the example 3d/hex model, if we want
>         to define a boundary condition on y_neg/y_pos , the fault nodes
>         be excluded az z_neg, right?
> 
> 
>     I don't understand this question.
> 
>     Regards,
>     Brad
>     _______________________________________________
>     CIG-SHORT mailing list
>     CIG-SHORT at geodynamics.org <mailto:CIG-SHORT at geodynamics.org> <mailto:CIG-SHORT at geodynamics.org <mailto:CIG-SHORT at geodynamics.org>>
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short>
>     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short>>
> 
> 
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org <mailto:CIG-SHORT at geodynamics.org>
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short

Charles Williams I Geodynamic Modeler
GNS Science I Te Pῡ Ao
1 Fairway Drive, Avalon 5010, PO Box 30368, Lower Hutt 5040, New Zealand
Ph 0064-4-570-4566 I Mob 0064-22-350-7326 I Fax 0064-4-570-4600
http://www.gns.cri.nz/ <http://www.gns.cri.nz/> I Email: C.Williams at gns.cri.nz <mailto:your.email at gns.cri.nz>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20180701/5ad5400e/attachment.html>