Eigenvalues and Eigenvectors of Right Cauchy-Green Tensor in GetFEM++

General
1

Dear GetFEM community,

I am currently working on a nonlinear solid mechanics problem using GetFEM++, and I would like to compute the eigenvalues and eigenvectors of the right Cauchy-Green deformation tensor , where is the deformation gradient.

Since is a symmetric positive-definite tensor, its spectral decomposition (principal stretches and directions) is particularly useful for my analysis.

My question is:
Is there any built-in functionality in GetFEM++ to compute the eigenvalues and eigenvectors of such a tensor field, either symbolically or numerically (e.g., at Gauss points or element-wise)?

If not directly available, I would appreciate any recommendations on how to extract the tensor at integration points and perform the eigen decomposition using external libraries (e.g., Eigen, NumPy, etc.).

Thank you in advance for your support and suggestions.

2

No, actually this we have not implemented, but we have implemented matrix logarithm function Logm (and also Expm), which can do as much as eigenvalue decomposition.

This is an excerpt from my notes on hyperelasticity, about a Logm/Expm based implementaion of the Ogden material law:

image
image636×270 9.94 KB

image
image636×632 25.1 KB

eigendecomposition (and its derivatives) are as you might know quite ugly in their implementation for the case of multiple eigenvalues. We can implement it if you are interested in, but I would need to find time for it.

3

Thank you very much for your answer, @Konstantinos.Poulios. I actually need it for a specific problem, and I’m quite sure it will also be helpful for the community. In the meantime, I’ll also try to explore some solutions on my side.
Thanks again!