Tensorial Sinbolos de Cristoffel

UNIVERSIDAD MAYOR DE SAN ANDRES CÓDIGO MATLAB SÍMBOLOS DE Christoffel
FACULTAD DE INGENIERIA
Docente: Msc.Ing. Mario Delgadillo Z LUGAR: FACULTAD DE INGENIERIA FECHA: 26/07/2019

UBICACIÓN[UTM]: 19K X= 592302.91 m E Y= 8175747.11 m S
UNIVERSIDAD MAYOR DE SAN ANDRÉS

FACULTAD DE INGENIERÍA
ANALISIS VECTORIAL Y TENSORIAL
CÓDIGO MATLAB SÍMBOLOS DE

CHRISTOFFEL
✓ PAYE CHIPANA JOSE CI: 6072256 L.P.
La Paz - Bolivia
2019

MATLAB SIMBOLOS DE Christoffel

syms cr cr1 r fron ph H3 H2 Haste Ecu u u1prima u2prima u1 u2 u3 u13
syms u23 u33 u11 u21 u31 u12 u22 u32 u111 u112 u113 u121 u122 u123
syms u131 u132 u133 u211 u212 u213 u221 u222 u223 u231 u232 u233 u311
syms u312 u313 u321 u322 u323 u331 u332 u333 %Definición de los Christoffel
cr(i,j,k)=[cr(1,1,1)=0;
cr(1,2,2)=-r*(sin(ph))^2;
cr(1,3,3)=-r;
cr(2,2,1)=1/r;
cr(2,1,2)=1/r;
cr(3,3,1)=1/r;
cr(3,1,3)=1/r;
cr(2,2,3)=cot(ph);
cr(2,3,2)=cot(ph);
cr(3,2,2)=-sin(ph)*cos(ph);
cr(3,3,3)=0;
%Derivada de los Christoffel cr1(i,j,k,l)=[i,j,k],l
cr1(1,1,1,1)=0;
cr1(1,2,2,1)=-(sin(ph))^2;
cr1(1,2,2,3)=-2*r*sin(ph)*cos(ph);
cr1(1,3,3,1)=-1;
cr1(2,2,1,1)=-1/r^2;
cr1(2,1,2,1)=-1/r^2;
cr1(3,3,1,1)=-1/r^2;
cr1(3,1,3,1)=-1/r^2;
cr1(2,2,3,3)=-(sec(ph))^2;
cr1(2,3,2,3)=-(sec(ph))^2;
cr1(3,2,2,3)=-cos(2*ph);
cr1(3,3,3,3)=0;
%Definición del vector desplazamiento
u(1)=u1;
u(2)=u2;
u(3)=u3;
%Definición de la primera derivada del vector desplazamiento
u1prima(1,1)=0;
u1prima(1,1)=u11;
u1prima(2,1)=u21;
u1prima(3,1)=u31;
u1prima(1,2)=u12;
u1prima(2,2)=u22;
u1prima(3,2)=u32;
u1prima(1,3)=u13;
u1prima(2,3)=u23;
u1prima(3,3)=u33;

%Definición de la segunda derivada del vector desplazamiento

u2prima(1,1,1)=u111;
u2prima(1,1,2)=u112;
u2prima(1,1,3)=u113;
u2prima(1,2,1)=u121;
u2prima(1,2,2)=u122;
u2prima(1,2,3)=u123;
u2prima(1,3,1)=u131;
u2prima(1,3,2)=u132;
u2prima(1,3,3)=u133;
u2prima(2,1,1)=u211;
u2prima(2,1,2)=u212;
u2prima(2,1,3)=u213;
u2prima(2,2,1)=u221;
u2prima(2,2,2)=u222;
u2prima(2,2,3)=u223;
u2prima(2,3,1)=u231;
u2prima(2,3,2)=u232;
u2prima(2,3,3)=u233;
u2prima(3,1,1)=u311;
u2prima(3,1,2)=u312;
u2prima(3,1,3)=u313;
u2prima(3,2,1)=u321;
u2prima(3,2,2)=u322;
u2prima(3,2,3)=u323;
u2prima(3,3,1)=u331;
u2prima(3,3,2)=u332;
u2prima(3,3,3)=u333;
%Definición del tensor efectivo de cuarto orden
syms h1111 h1122 h1133 h2211 h2222 h2233 h3311 h3322 h3333 h1212 h1313 h2323 H
H(1,1,1,1)=h1111;H(1,1,2,2)=h1122;
H(1,1,3,3)=h1133;H(2,2,1,1)=h1122;
H(2,2,2,2)=h1111;H(2,2,3,3)=h1133;
H(3,3,1,1)=h1133;H(3,3,2,2)=h1133;
H(3,3,3,3)=h3333;H(1,2,1,2)=h1212;
H(1,2,2,1)=h1212;H(2,1,2,1)=h1212;
H(2,1,1,2)=h1212;H(1,3,1,3)=h1313;
H(1,3,3,1)=h1313;H(3,1,1,3)=h1313;
H(3,1,3,1)=h1313;H(2,3,2,3)=h1313;
H(2,3,3,2)=h1313;H(3,2,2,3)=h1313;
H(3,2,3,2)=h1313;
%Aplicación de las ecuaciones para tensores de orden inferior
for i=1:3
for k=1:3

for l=1:3
H3(i,k,l)=0;
for m=1:3
for n=1:3
H3(i,k,l)=H3(i,k,l)-cr(l, m, n)*H( i, k, m, n );
end
end
end
end
end for i=1:3
for k=1:3
for l=1:3
for r=1:3
for j=1:3
H3(i,k,l)=H3(i,k,l)+cr(i, j, r)*H( r, j, l, k );
end
end
end
end
end for i=1:3
for k=1:3
for l=1:3
for r=1:3
for j=1:3
H3(i,k,l)=H3(i,k,l)+cr(j, j, r)*H( i, r, l, k );
end
end
end
end
end for i=1:3
for l=1:3
H2(i,l)=0;
for j=1:3
for m=1:3
for n=1:3
H2(i,l)=H2(i,l)-H( i, j, m, n )*cr1( l, m, n, j );
end
end
end
end
end for i=1:3
for l=1:3
for j=1:3
for m=1:3

for n=1:3
for r=1:3
H2(i,l)=H2(i,l)-cr(i,j,r)*H( j, r, m, n )*cr(l,m,n);
end
end
end
end
end
end for i=1:3
for l=1:3
for j=1:3
for m=1:3
for n=1:3
for r=1:3
H2(i,l)=H2(i,l)-cr(j,j,r)*H( r, i, m, n )*cr(l,m,n);
end
end
end
end
end
end for i=1:3
for j=1:3
for l=1:3
Haste(i,j,l)=0;
for m=1:3
for n=1:3
Haste(i,j,l)=Haste(i,j,l)-cr(l, m, n)*H( i, j, m, n );
end
end
end
end
end
%Construcción de la ecuación homogeneizada for i=1:3
Ecu(i,1)=0;
for j=1:3
Ecu(i,1)=Ecu(i,1)+H2(i,l)*u(l)
; end
end for i=1:3
for k=1:3
for l=1:3
Ecu(i,1)=Ecu(i,1)+H3(i,k,l)*u1prima(l,k);
end
end
end for i=1:3

for j=1:3
for l=1:3
for k=1:3
Ecu(i,1)=Ecu(i,1)+H(i,j,l,k)*u2prima(l,k,j);
end
end
end
end %Construcción de la condición de frontera en $\Sigma_2$
syms n1 n2 n3 n
n(1)=n1;
n(2)=n2;
n(3)=n3;
for i=1:3
fron(i,1)=0;
for j=1:3
for k=1:3
for l=1:3
fron(i,1)=fron(i,1)+H(i,j,l,k)*u1prima(l,k)*n(j);
end
end
end
end for i=1:3
for j=1:3
for l=1:3
fron(i,1)=fron(i,1)+Haste(i,j,l)*u(l)*n(j);
end
end
end

Tensorial Sinbolos de Cristoffel

Cargado por

Copyright:

Formatos disponibles

Tensorial Sinbolos de Cristoffel

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tensorial Sinbolos de Cristoffel

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD MAYOR DE SAN ANDRES CÓDIGO MATLAB SÍMBOLOS DE Christoffel

Docente: Msc.Ing. Mario Delgadillo Z LUGAR: FACULTAD DE INGENIERIA FECHA: 26/07/2019

UNIVERSIDAD MAYOR DE SAN ANDRÉS

CÓDIGO MATLAB SÍMBOLOS DE

✓ PAYE CHIPANA JOSE CI: 6072256 L.P.

Docente: Msc.Ing. Mario Delgadillo Z LUGAR: FACULTAD DE INGENIERIA FECHA: 26/07/2019

MATLAB SIMBOLOS DE Christoffel

Docente: Msc.Ing. Mario Delgadillo Z LUGAR: FACULTAD DE INGENIERIA FECHA: 26/07/2019

%Definición de la segunda derivada del vector desplazamiento

Docente: Msc.Ing. Mario Delgadillo Z LUGAR: FACULTAD DE INGENIERIA FECHA: 26/07/2019

Docente: Msc.Ing. Mario Delgadillo Z LUGAR: FACULTAD DE INGENIERIA FECHA: 26/07/2019

Docente: Msc.Ing. Mario Delgadillo Z LUGAR: FACULTAD DE INGENIERIA FECHA: 26/07/2019

También podría gustarte