Dear Cummunity,
I want to do some post processing in paraview.
In my exodus file I have 9 stress components stored both as the vectors and scalers. As the vector I have variables named as stress_x
, stress_y
, and stress_z
. But as the scaler I have stress_x_X
, stress_x_Y
, stress_x_Z
, stress_y_X
, stress_y_Y
, stress_y_Z
, stress_z_X
, stress_z_Y
, and stress_z_Z
. On one of the (2D) surfaces (named as fault
) in domain I want to calculate the normal stress. The first thing is that I should find the three components of the normal on the surface (fault
). I followed this discussion but I could find a solution for it. In each triangle on the surface I need to calculate the three coefficients whoen in the following figures as n1
, n2
and n3
. Then, the following matrix multiplication will allow to run a matrix calculation:
The t_vector
can be finally used to calculate the normal stress value as a scaler (by the dot production of t_vector
and normal coefficients):
I have also uploaded a zipped exodus file here.
Thanks a lot for any support regarding my issue for normal stress calculation. At the end I should have one normal stress value for each triangle of the surface.