principal strains using eig function

3 次查看(过去 30 天)
BLP
BLP 2015-6-18
回答: Eric Lin 2015-6-19
I am trying to use the eig function to calculate the principal strains from 3D strain tensor. Unfortunately I do get different results than from another tool (<http://www.continuummechanics.org/cm/techforms/Eigen.html>) For
A1 = [4.420000000000000 1.350000000000000 3.509000000000000;
1.350000000000000 3.832000000000000 0;
3.509000000000000 0 3.207000000000000]
Matlab eig function yields
>> [V1, D1] = eig(A1)
V1 =
-0.651153549926562 0.054428335246106 0.756991816825237
0.232116359059687 -0.935349177129471 0.266915553499784
0.722579442189515 0.349513194551449 0.596422230101738
D1 =
0.044859378106921 0 0
0 3.753442970840320 0
0 0 7.660697651052759
while the Web tool gives:
The Values look pretty the same, but their position in vector/matrix is different. As far I understand, the columns of the V1 represent the cosines of the angle between the a principal strain and the Cartesian axis. But which corresponds to which principal strain and how to calculate angles between the principal strains and the Cartesian strain axis?

回答(1 个)

Eric Lin
Eric Lin 2015-6-19
What is happening is the web app is displaying the eigenvectors by row whereas MATLAB is displaying them by column. The relative order of the eigenvectors is not significant.
I am not familiar with the domain but from a quick search it seems that the vectors are related to the strains based on their eigenvalue, i.e. the max principal strain is the eigenvector with the largest eigenvalue. Again, I am not familiar with the domain but angles between vectors can be computed by dot product.

类别

Help CenterFile Exchange 中查找有关 Linear Algebra 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by