[CIG-SHORT] Pylith stress output

Sun Apr 2 16:59:40 PDT 2017

Dear Teng,

A quick look shows that your elements appear to be ordered starting in the middle, row-by-row on the right, and then row-by-row on the left (see attached image).  You could probably modify the script to take this into account, but I believe that any script that relies on elements being numbered in a particular way is going to be very fragile.  I think that a more robust method would be to reorder the cells yourself before plotting.  This would be easy to do, since you also have access to the vertex and cell information.

Cheers,
Charles

p.s.  How was the mesh generated?

> On 3/04/2017, at 11:31 AM, Li, Teng <tengli2 at illinois.edu> wrote:
> 
> Hi Charles,
> I just used the HDF5 writer for the stress output. Now I am using Matlab to post process stress in h5 file. And I find that the size of the stress is (3,1600,205), in which 3 are stress components, 1600 are number of elements, 205 is the total time steps.
> I try to plot stress in 2d domain, however, I do not know the element order. So I can not determine the position of element. 
> Here I attach the file in case you will run Pylith and use the script Plot.m to open material.h5 file to see the stress. 
> 
> 
> 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>
> 
> 
> From: CIG-SHORT [cig-short-bounces at geodynamics.org <mailto:cig-short-bounces at geodynamics.org>] on behalf of Charles Williams [willic3 at gmail.com <mailto:willic3 at gmail.com>]
> Sent: Sunday, April 02, 2017 5:41 PM
> To: cig-short at geodynamics.org <mailto:cig-short at geodynamics.org>
> Subject: Re: [CIG-SHORT] Pylith stress output
> 
> Dear Teng,
> 
> The components correspond to xx, yy, and xy — see Table 5.2 in the manual.  Also, it is probably more reasonable to use HDF5 output rather than VTK output.
> 
> Cheers,
> Charles
> 
> 
>> On 2/04/2017, at 12:04 PM, Li, Teng <tengli2 at illinois.edu <mailto:tengli2 at illinois.edu>> wrote:
>> 
>> Hi,
>> In a 2d plane strain problem, I tried to output stress in VTK video form. And I find there are three components in the Paraview which are x, y, z. So what is the meaning of these three components? Three normal stress in x, y, z directions? However I think we should also have shear stress in the plane. So how to output all the stress components?
>>  
>> 
>> 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 <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>
> <fault_stress.exo><batch.pbs><Makefile.am><matprops.spatialdb><secondary_fault.cfg><traction.spatialdb><traction1.spatialdb><material.h5><Plot.m>_______________________________________________
> 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>
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/20170403/40e5bfeb/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cell_order.png
Type: image/png
Size: 13874 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170403/40e5bfeb/attachment-0001.png>