Hi Brad,

Using PyLith, can we directly extract/output the surface traction vectors (shear-horizontal, shear-vertical, normal) stresses for any arbitrary surface defined using a subdomain?

Or do we necessarily need post processing to get them.

Thank you

Viven

PyLith does not currently support extracting surface traction vectors on an arbitrary surface. You need to either

- Output the stress tensor and compute the surface tractions in post processing, or
- Output the tractions on a fault surface (the slip can be set to zero).

Thank you for the reply.

Can we impose Dirichlet BC on a surface and at the same time define it as an interface to extract traction vector?

Viven

No. Dirichlet BC can only be applied to external boundaries. Faults can only be defined for interior interfaces.