Hi,
I am building a 2d/3d elasticity model. By following the tutorial (e.g., Beam_3d.edp), I can solve the stresses {Sxx, Syy, Sxy} (in 2d) or {Sxx, Syy, Szz, Sxy, Syz, Sxz} (in 3d).
Can I further calculate the principal stresses for all elements? Is using EigenValue function a right direction?
Like solving:
( A - \lambda I )x = 0
where A is the stress tensor:
A = [Sxx, Sxy, Sxz,
Sxy, Syy, Syz,
Sxz, Syz, Szz]
\lamba
is 3 eigenvalues
Thank you very much!
Yan