Principal Tensor directions error
3 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
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 comentarios
Respuestas (0)
Ver también
Categorías
Más información sobre Condensed Matter & Materials Physics en Help Center y File Exchange.
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!