[CIG-SHORT] Output Material Subset

[Brad Aagaard]

Scott,

Extracting slices of state variables like stress and strain from a 3-D 
model is best done inside visualization software. Nearly all 3-D 
visualization tools have this feature. This gives you the most freedom 
to interactively adjust where the slice is and extract either a slice or 
clip the domain at a particular slice.

Regards,
Brad


On 12/01/2013 03:59 PM, Charles Williams wrote:
> 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.
>
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>