Principal Tensor directions error
5 次查看(过去 30 天)
显示 更早的评论
Hi !
I'm currently working in a solver for the principal direction of a stress tensor, but while solving the system for the directions a get a fail error. Please if someone can help me with any information it would be great.
Solving manually i got this
My current work is here:
clear all clc
%Input
disp('Informe as componentes do tensor ')
T11=input(' Elemento 1;1 ');
T12=input(' Elemento 1;2 ');
T13=input(' Elemento 1;3 ');
T21=input(' Elemento 2;1 ');
T22=input(' Elemento 2;2 ');
T23=input(' Elemento 2;3 ');
T33=input(' Elemento 3;3 ');
%Tensor informado
T=[T11 T12 T13;T21 T22 T23;T13 T23 T33]
%Solução Tensões/Deformações Principais syms x TL=[T11-x T12 T13;T21 T22-x T23;T13 T23 T33-x]; deter=det(TL); princ=solve(deter,x)
%Solução Direções Principais %syms N1 N2 N3 %A=[(T11-princ(1)) T12 T13;T21 (T22-princ(1)) T23;T13 T23 (T33-princ(1))]; %B=[0;0;0]; %N=[N1;N2;N3]; %EQ=A*N; %Dire_1=vpasolve(A*N==B,[N1 N2 N3])
A=[(T11-princ(1)) T12 T13;T21 (T22-princ(1)) T23;T13 T23 (T33-princ(1))]; B=[0;0;0]; Direcao_1=inv(A)*B
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Condensed Matter & Materials Physics 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!