top of page
If using this toolbox for research or industrial purposes, please cite:
Advances in Engineering Software. Vol 105. March 2017. Pages 9-16. (2017)
Abaqus2Matlab v.3.0
A new way to post-process FEA
bottom of page
Dear Yadav,
The orientation associated with each principal stress is given in terms of the three components of the vector that is parallel to the principal stress considered. The three components of this vector are simply its projections along the global X-, Y- and Z- axes.
You can find the direction cosine angles (in radians) of the orientation vector of each principal stress by the following code, as a continuation of the code of my previous post:
% Direction cosine angles of the orientation of the first % principal stress (rad)
phi11=acos(or1(1)/norm(or1));
phi12=acos(or1(2)/norm(or1));
phi13=acos(or1(3)/norm(or1));
% Direction cosine angles of the orientation of the second
% principal stress (rad)
phi21=acos(or2(1)/norm(or2));
phi22=acos(or2(2)/norm(or2));
phi23=acos(or2(3)/norm(or2));
% Direction cosine angles of the orientation of the third
% principal stress (rad)
phi31=acos(or3(1)/norm(or3));
phi32=acos(or3(2)/norm(or3));
phi33=acos(or3(3)/norm(or3));