[CIG-SHORT] Output Material Subset

[Charles Williams]

You can definitely output the stresses at a desired set of points, but it won’t be resolved onto a surface because there’s no orientation information.  For a surface consisting of a nodeset, there also isn’t any orientation information.  You only get that if the surface is a fault.  I would say that your options are:

1.  Output the stresses to a nodeset and then do your own post processing to get the normal and shear components.
2.  Make the surface for which you want the stresses into a fault (it can have zero slip), and the stresses will be resolved onto it.

Charles


On 2/12/2013, at 11:42 am, Matthew Knepley <knepley at mcs.anl.gov> wrote:

> On Sun, Dec 1, 2013 at 4:37 PM, Scott Henderson <sth54 at cornell.edu> wrote:
> In Pylith is it possible to output stresses resolved onto an arbitrary surface (not a defined fault) for a 3D model?
> 
> I assume you can do this with OutputSolnPoints, just like its done in the GreenFunc examples.
> 
>    Matt
>  
> I’m imagining something like OutputSolnSubset for materials where you specify a Cubit nodeset based on a material surface. For example, this would be useful for plotting stress contours in a profile through a 3D domain:
> 
> #---------------
> [pylithapp.problem.materials.elastic]
> output =  pylith.meshio.OutputSolnSubset
> 
> [pylithapp.problem.materials.elastic.output]
> label = face_y
> cell_filter = pylith.meshio.CellFilterAvgSubMesh
> cell_data_fields = [total-strain, stress]
> writer = pylith.meshio.DataWriterHDF5SubMesh
> writer.filename = output/step01/NorthSouthMaterialCrossSection.h5
> #———————
> 
> This is similar to the question of getting stresses at specified points (http://www.geodynamics.org/pipermail/cig-short/2012-October/001165.html). So I know it’s possible to do this as post-processing in Python or Matlab. Maybe there is also an easy way to export ’slices' of 3D data from Paraview too? Any suggestions would be appreciated,
> 
> Thanks,
> Scott
> 
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
> 
> 
> 
> -- 
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
> -- Norbert Wiener
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://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

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


Notice: This email and any attachments are confidential.
If received in error please destroy and immediately notify us.
Do not copy or disclose the contents.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/cig-short/attachments/20131202/cdb5d52d/attachment-0001.html>