[CIG-SHORT] Stress discontinuity at fault intersection in 3D
Brad Aagaard
baagaard at usgs.gov
Fri Feb 23 12:25:28 PST 2018
Josimar,
You can output the stress tensor in the domain and use post-processing
to compute the stress changes at whatever points you want and then
resolve them into tractions for a specific orientation. We don't have
utilities to do this, so you would need to write some (Python) code to
do this post-processing.
Regards,
Brad
On 2/23/18 11:35 AM, Josimar Alves da Silva wrote:
> Hi Brad,
>
> I have a question regarding fault intersection in 3D. This is a follow
> up question from a previous recent discussion.
>
> When I look at the normal stress at one of my faults I notice a
> discontinuity in the stresses at the fault intersection. See attached
> figure for details.
>
> This is of course a result of the fault_edges definition as you
> explained in a recent post. My questions is:
>
> i) Is it possible to output normal and shear stresses on a plane inside
> the domain without having to define those planes as faults in PyLith ?
> This would perhaps bypass the fault intersection problem for cases where
> one only wants stresses on the fault. The goal is to avoid having to
> define fault buried edges that result in the attached stress discontinuity.
>
> thank you in advance,
> Josimar
>
>
>
More information about the CIG-SHORT
mailing list