[CIG-SHORT] Fault output normal traction

[Li, Teng]

Dear Charles,

I have a 2d rectangular model and I use the slip weakening model.

Please find the attachments which are the fault.h5 file.

And I use the following Matlab code to find the output:


hinfoh=hdf5info('/Users/tengli/pylith/Local/fault.h5');

coords=hdf5read(hinfoh.GroupHierarchy.Groups(1).Datasets(1));

shear=hdf5read(hinfoh.GroupHierarchy.Groups(3).Datasets(3));


And I find there 2 rows in shear data. So I guess one is the shear traction and another is normal traction.


Best,

Teng


Teng Li

Master Candidate in Structures

Department of Civil and Environmental Engineering

University of Illinois at Urbana-Champaign

205 North Mathews Ave, Urbana, IL. 61801

Phone:(217)8196210, Email: tengli2 at illinois.edu



________________________________
From: CIG-SHORT [cig-short-bounces at geodynamics.org] on behalf of Charles Williams [willic3 at gmail.com]
Sent: Thursday, August 03, 2017 10:59 PM
To: cig-short at geodynamics.org
Subject: Re: [CIG-SHORT] Fault output normal traction

Dear Teng,

There should be more than one component of traction.  What is the dimension of your problem, and what slip model are you using?

Cheers,
Charles



On 4/08/2017, at 2:09 PM, Li, Teng <tengli2 at illinois.edu<mailto:tengli2 at illinois.edu>> wrote:

Hi!

I would like to output normal traction of the fault in my model. So I use:

vertex_data_fields = [slip,slip_rate,traction]
writer = pylith.meshio.DataWriterHDF5
writer.filename = output/fault.h5

Finally I got the fault.h5 file.
However, in the fault.h5 file, I find the traction is only the shear traction along the fault. So I was wondering how I can output the normal traction along the fault?

Best,
Teng



Teng Li

Master Candidate in Structures

Department of Civil and Environmental Engineering

University of Illinois at Urbana-Champaign

205 North Mathews Ave, Urbana, IL. 61801

Phone:(217)8196210, Email: tengli2 at illinois.edu<mailto:tengli2 at illinois.edu>



_______________________________________________
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

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/ 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/20170804/c4c2ff5d/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fault.h5
Type: application/octet-stream
Size: 7880720 bytes
Desc: fault.h5
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170804/c4c2ff5d/attachment-0001.obj>