ut é ù é ù é ù ê ú ê ú ê ú = ¾¾¾ = =ê ú ê ú ê ú ë û ë û ë û
ut é ù é ù é ù ê ú ê ú ê ú = ¾¾¾ = =ê ú ê ú ê ú ë û ë û ë û
ut é ù é ù é ù ê ú ê ú ê ú = ¾¾¾ = =ê ú ê ú ê ú ë û ë û ë û
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Teoría de Campos - § 1.4 Tensores de<br />
2do. orden<br />
§1.4 ESTUDIO DE LOS TENSORES DE (2)<br />
PANORAMA<br />
•a) Transpuesto de un tensor arbitrario: trasposición en (2)<br />
•b) y c) Determinante y traza de un tensor cualquiera de segundo orden<br />
•d) Inverso de un tensor regular (o a<strong>ut</strong>omorfismo): inversión en A<strong>ut</strong>()<br />
–el grupo lineal de transformaciones GL(n,) := A<strong>ut</strong>()<br />
•e) A<strong>ut</strong>oanálisis de un tensor: a<strong>ut</strong>ovalores y a<strong>ut</strong>ovectores<br />
–Definiciones, cálculo, interpretación geom<strong>é</strong>trica<br />
–Aplicaciones: funciones analíticas de tensores diagonalizables.<br />
•f) Tensores sim<strong>é</strong>tricos y antisim<strong>é</strong>tricos<br />
–Definición y propiedades especiales<br />
–Descomposición espectral de un tensor sim<strong>é</strong>trico: direc. principales del tensor.<br />
–Vector axial de un tensor antisim<strong>é</strong>trico<br />
–Parte sim<strong>é</strong>trica y p. antisimª. de un tensor cualq.:<br />
•g) Tensores ortogonales: rotaciones y simetrías especulares<br />
–Def. y propiedades especiales. el subgrupo ortogonal O(n, ) GL(n,)<br />
•h*) El teorema de descomposición polar*: T : T = D·Q, donde D = (Tt ·T) 1/2<br />
es sim. (parte deformante), y Q = D1 <br />
: S<br />
<br />
: A<br />
t t<br />
(2) T T T T<br />
T : T <br />
2 2<br />
·T es ortogonal (parte reorientadora)<br />
a) Tensor traspuesto en (2) . Trasposición tensorial<br />
•DEFINICIÓN (intrínseca.): "TÎ (2) $!U := T t Î (2) / "xÎ : x · T t = T · x<br />
–Ejemplo: (ab) t = ba , porque: x · ba := (x·b) a = a (b·x) = ab · x<br />
•PROPIEDADES de la TRASPOSICIÓN TENSORIAL:<br />
–linealidad: (S+T) t = S t + T t<br />
–Idempotencia: (T t ) t = T x · T = T t · x (así x·T es una transf. lin. de x, q. coincide con T t )<br />
–traspuesto de un producto contraído: (S·T) t = T t · S t<br />
–determinante* del trasp.: det(T t ) = detT (se repasa determinante de un tensor en apartado c)<br />
–inverso*: (T t ) = (T ) t := T t (se repasa inverso de un tensor regular en apartado b)<br />
•CÁLCULO en compntes. generales: U := T t = (t i j i j ) t = t i j j i := u j i j i = u i j i j<br />
•o sea: u j<br />
i = t j<br />
i matricialmente se dispone del algoritmo de trasponer o trasposición:<br />
t<br />
<strong>é</strong> 1 2 3 1 2 3 1 1 1<br />
u1 u1 u <strong>ù</strong> <strong>é</strong> 1 t 1 t 1 t <strong>ù</strong> <strong>é</strong> 1 t 1 t 2 t <strong>ù</strong><br />
j j<br />
3<br />
1 2 3 1 2 3 2 2 2<br />
ui = t i¾¾¾ <strong>ê</strong>u i ,j 2 u2 u <strong>ú</strong> 2 = <strong>ê</strong>t 2 t 2 t <strong>ú</strong> 2 = <strong>ê</strong>t 1 t 2 t <strong>ú</strong><br />
= f = c<br />
3<br />
<strong>ê</strong> 1 2 3 1 2 3 3 3 3<br />
u u u <strong>ú</strong> <strong>ê</strong>t t t <strong>ú</strong> <strong>ê</strong>t t t <strong>ú</strong><br />
<strong>ë</strong> 3 3 3<strong>û</strong> <strong>ë</strong> 3 3 3<strong>û</strong> <strong>ë</strong> 1 2 3<strong>û</strong><br />
•Análogos: u ij = t ji ; u i j = t j i ; u ij = t ji ; análogo tb. en cartesianas: escr. matricial<br />
•CASOS PARTICULARES: tensores sim<strong>é</strong>tricos y antisim<strong>é</strong>tricos<br />
–tensores sim<strong>é</strong>tricos: T = sim. T t = T (se estudian en particular más adelante)<br />
–tensores antisim<strong>é</strong>tricos: T = antisim. T t = T (se estudian tambi<strong>é</strong>n más adelante)<br />
•Teorema de descomposición: T (2) : T = ½(T+T t ) + ½ (TT t ), llamadas,<br />
respectivamente, parte sim<strong>é</strong>trica y parte antisim<strong>é</strong>trica de T.<br />
•Ejemplo: Descomponer el tensor T de componentes [t i j] = [1 , 3; -2 , 1] siendo G = [1 , 1; 1 , 2].<br />
Curso 2012-13 1
Teoría de Campos - § 1.4 Tensores de<br />
2do. orden<br />
b) Tensor inverso. Inversión tensorial<br />
•DEF. intríns.: "TÎ (2) , regular, $!V := T 1 := apl. lin. inv. / u Î:<br />
w = T·u V·w = u T · V = V · T = 1 (inverso resp. " · ")<br />
–Grupo lineal de transf. o tens. reg.:= Gl(n, ) := (A<strong>ut</strong>(), · ), grupo<br />
multiplicativo con producto "dot" o pr. contraído.<br />
•CÁLCULO en comptes.: → algor. inversa matr. (obs: ¡ da 1 en forma ij !):<br />
–cartes.: V ihT hj = T ihV hj = ij [V ij] = [T ij] 1<br />
–generales mix.: v i ht h j = v i h th j = i j = i j [v i j] = [t i j] 1 .; [v i j ] = [ti j ] 1 .;<br />
–generales puras: v ih t hj = i j ; v ih thj = i j [v ij] = [t ij ] ; [v ij ] = [t ij] <br />
•PROPIEDADES de la inversión tensorial:<br />
–Idempotencia: (T ) = T<br />
–Inverso del m<strong>ú</strong>ltiplo escalar: (A) = (1/) A (la inversión no es lineal)<br />
–Inverso del producto contraído de tens. reg. de (2) : (S·T) = T · S <br />
–determinante: det(A ) = 1/detA<br />
–inverso del traspuesto: (T t ) = (T ) t := T t<br />
•Ejemplos: 1) ¿Existe (a×) —1 ? 2) Si R es una rotac. de rad. alr. de e , identifica el<br />
tensor R . 3) Misma cuestión si S es una simª. resp. {e} . 4) Tens. ortogonal Q 1 = Q t<br />
(más adelante: diapositiva 9)<br />
c y d) determinante. y traza de un tensor<br />
DEF.<br />
Propie<br />
dades<br />
Cálculo<br />
c) determinante d) Traza<br />
T· g1, T· g2, T·<br />
g <br />
3<br />
vol(base<br />
transf.)<br />
g1, g2, g <br />
3<br />
vol(base<br />
orig.)<br />
det T : <br />
Ej.: det1 = 1; det(a×) = 0 (base transformada es<br />
coplanaria.)<br />
•det(T) = n detT (det no es lineal)<br />
•det(S·T) = detS detT<br />
•det(Tt ) = detT<br />
•det(T1 ) = 1/detT (Nota: n = dim)<br />
detT = det[t i j] = det[t i j ] = det[tij]g 1 = det[t ij ]g<br />
•obs: Distinguir entre el algoritmo determ te . de<br />
una matriz (regla para operar), del<br />
determinante de un tensor (coeficiente de<br />
dilatación signado de la definición dada arriba)<br />
Def.: tr(ab) := a·b ;<br />
se extiende por linealidad a T<br />
Ej.: tr1 = 3 (ó 2); tr(a×) = 0<br />
•tr(A) = trA<br />
•tr(A+B) = trA + trB<br />
•tr(A t ) = tr(A)<br />
(la traza sí es lineal)<br />
• trT:= T ijtr(e ie j) = T ij ij = T ii;<br />
• trT= t i j i j = t i i ; trT = t ijg ij<br />
• trT= t i j i j = t i i ; trT = t ij gij<br />
obs.: trT = T··1 (doble contracción,<br />
que se define sobre tensores de (2))<br />
•Ambos son coeficientes del pol. característico y son valores intrínsecos de cada T<br />
•Carácter invariante* por c. de base: todo el polinomio característico es invariante porque: ()<br />
= det([ i j][i j]) = det(C-1 ·[ti j]·C – C1 ·[i j]·C) = … = det([ti j] [i ˆt<br />
j])<br />
Curso 2012-13 2
Teoría de Campos - § 1.4 Tensores de<br />
2do. orden<br />
e) A<strong>ut</strong>ovalores y a<strong>ut</strong>ovectores: 1.- definiciones (!!)<br />
•Definiciones y cálculos<br />
• Î es a<strong>ut</strong>ovalor de T $uÎ, u ≠ 0 / T·u = ·u (T-1)·u = 0 ; u := -a<strong>ut</strong>ovector <br />
T·u | | u (proporcional) el s.v. L({u}) se conserva al actuar T (s.v. unidimensional<br />
invariante por T).<br />
• u = -a<strong>ut</strong>ovector u Î ker(T1) := () = -a<strong>ut</strong>oespacio (subespacio vectorial, invariante<br />
por la acción del tensor, asociado al a<strong>ut</strong>ovalor ). Son lin. ind. los a<strong>ut</strong>ovec. de distinto a<strong>ut</strong>oval.<br />
• Polinomio característico: () := det(T 1) := polinomio característico de T. Es intrínseco<br />
• espectro de T, (T) = {Î / = a<strong>ut</strong>ovalor de T} = 1 (0) = {Î / () = 0} ( los `s<br />
son las raíces reales de la ec. caract.)<br />
•Subesp. vec. invariantes : rectas invariantes () y planos inv.; <strong>é</strong>stos, pueden ser:<br />
de<br />
sólo 2 rectas inv. dos a<strong>ut</strong>oval. simples distintos<br />
(diap.12)<br />
<br />
pl. invar. de<br />
todas las rectas inv. un a<strong>ut</strong>oval. real doble<br />
<br />
sin<br />
rectas invariantes un par de a<strong>ut</strong>oval. compl. conjugados<br />
(diap.14)<br />
• Índice de Jordan de , () : = dim(). Existirá una base de formada por a<strong>ut</strong>ovectores de<br />
T "Î (T) : () = () := orden de multiplicidad de la raíz de () T es<br />
diagonalizable. En general sólo se puede asegurar que () ().<br />
Observaciones: 1) = 0 (T) kerT ≠ {0} T es singular detT = 0 → (0) = kerT.<br />
2) Los coeficientes del polinomio característico, () = 3 + J 1 2 J 2+ J 3, o invariantes de Jordan<br />
de T, resultan tres escalares, que en función de los a<strong>ut</strong>ovalores son: J 1 = tr(T) = 1+ 2+ 3,<br />
J 2 = 1 2+ 2 3+ 1 3 , J 3 = det(T) = 1 2 3, (otro recurso de cálculo para detT y trazaT).<br />
e) A<strong>ut</strong>ovalores y a<strong>ut</strong>ovectores: 2.- ejemplos<br />
•Ejemplos: Aplicando el concepto geom<strong>é</strong>trico (interpretación geom<strong>é</strong>trica o<br />
intrínseca), determinar los a<strong>ut</strong>ovalores y a<strong>ut</strong>ovectores de los siguientes tensores:<br />
–una rotación R(e ; ): plano invariante sin rectas invariantes (pareja de "a<strong>ut</strong>ovalores"<br />
complejos conj.); recta invariante de a<strong>ut</strong>ovalor real = 1: el eje L({e})<br />
–una simetría respecto de un plano {e} : pl. inv. de rectas invariantes (a<strong>ut</strong>ovalor real<br />
doble = 1); recta invariante de a<strong>ut</strong>ovalor = 1).<br />
–el tensor unidad 1 (uso de la definición directamente); una homotecia de razón k, k1<br />
(a<strong>ut</strong>ovalor real triple)<br />
–una simetría respecto de un eje ({e}): ¿existe plano invariante de rectas<br />
invariantes? [N.: en realidad el tensor de simetría axial es una rotación de ángulo ].<br />
–el tensor axial e: analizar sus a<strong>ut</strong>ovalores a partir de interpretar geom<strong>é</strong>tricamente su<br />
acción.<br />
•S.v. invariantes de un tensor TÎ (2) , permiten describirlo intrínsecamente: de esta<br />
forma se puede recuperar la inf. intrínseca de T a partir de cualquier representación<br />
matricial contra-covariante, [t i j], efectuando el a<strong>ut</strong>oanálisis de la matriz.<br />
•Pueden hacerse los ejercicios: PR1.3, 6, 8, 9 (apt. 1 y 2) y 15 (cuestión y apartado 1) y<br />
apartados particulares de otros, que se indicarán cuando se puedan hacer completos).<br />
Curso 2012-13 3
Teoría de Campos - § 1.4 Tensores de<br />
2do. orden<br />
e) A<strong>ut</strong>ovalores y a<strong>ut</strong>ovectores: 3.- aplicaciones (!!)<br />
•Aplicaciones (1): Diagonalización<br />
• Un tensor T es diagonalizable {w i} base de formada por a<strong>ut</strong>ovectores de T<br />
(a<strong>ut</strong>obase o base propia). En una tal base T tiene una matriz mixta diagonal, llamada<br />
una forma diagonal de la representación matricial de T. Formas diagonales<br />
alternativas del mismo tensor se obtienen sólo reordenando la base.<br />
•Son diagonalizables con certeza los tensores con tres a<strong>ut</strong>ovalores reales simples y<br />
los tensores sim<strong>é</strong>tricos (en este caso la a<strong>ut</strong>obase puede conseguirse ortogonal y la<br />
a<strong>ut</strong>obase normalizada se llama base espectral y tambi<strong>é</strong>n direcciones principales)<br />
•No son diagonalizables sobre los tensores con a<strong>ut</strong>ovalores complejos.<br />
•Son dudosos los tensores con a<strong>ut</strong>ovalores reales m<strong>ú</strong>ltiples. Discusión: Se estudiará<br />
si () = (), en cada caso, y si coinciden, el tensor es diagonalizable (si y sólo<br />
si).<br />
•Aplicaciones (2): cálculos que se simplifican<br />
•Potencias: Calcular la matriz A m siendo A = [a,0,0; 0,b,0; 0,0,c] y m, un entero.<br />
Deducir la matriz raíz m-<strong>é</strong>sima de una matriz diagonal: [a,0,0; 0,b,0; 0,0,c] 1/m . Una<br />
vez calculada, se puede regresar a la base original. Ejemplo: PR1.29 (examen Spt11).<br />
•Polinomios: Calcular el polinomio f(x) = 1 + 2x–3x 2 + x 3 aplicado a la matriz<br />
diagonal A del ejemplo anterior.<br />
•Funciones analíticas de un tensor, f(T) : (sigue…)<br />
•Aplicaciones (3 continuación): evaluación de funciones analíticas de<br />
tensores<br />
•Se dice que una función f(A) se conserva por semejanza, si:<br />
f(C 1 ·A·C) = C 1 ·f(A)·C , "A matriz, C matriz regular<br />
•Ejemplos: Si f es una potencia, f(A) = An , o un polinomio o una serie de potencias<br />
(func. analíticas, como exp, Ln, sen, cos, tan,..) entonces f conserva la semejanza.<br />
•Si f conserva la semejanza y D es la forma diagonal de un tensor T del que se tiene<br />
una matriz mixta T / D = C1 ·T·C, entonces:<br />
f(D) =C1 · f(T)·C f(T) =C· f(D)·C1 3 1<br />
•Ejemplo: Calcular la raíz cuadrada de T ~ , dado en una base canónica.<br />
1 3<br />
<br />
<br />
Solución:…<br />
•Funciones analíticas: Se puede calcular exp(A) mediante la serie de potencias x n /n!<br />
•Aplicaciones (4): análisis de la acción de un tensor<br />
•Cálculo gráfico de T·v en t<strong>é</strong>rminos de a<strong>ut</strong>ovalores y a<strong>ut</strong>ovectores de T: la acción de T<br />
se puede descomponer en las direcciones de sus a<strong>ut</strong>ovectores o subespacios invariantes<br />
–por ejemplo, describir la acción del tensor 1 + ee en el espacio 3.<br />
•Pueden hacerse: PR1 n os : 6, 9-3), 11, 12 y 15-2),<br />
Curso 2012-13 4
Teoría de Campos - § 1.4 Tensores de<br />
2do. orden<br />
f) Tens. sim<strong>é</strong>tricos, antisim. y g1) Tens. ortogonales<br />
Df. carctrzn.<br />
Prop. fundamental<br />
Otras propiedades<br />
SIMÉTRICOS ANTISIMÉTRICOS ORTOGONALES<br />
S es sim<strong>é</strong>trico S t = S "v : v·S<br />
= S·v "u, v : u·S·v = v·S·u<br />
descomposición espectral (!!):<br />
S sim<strong>é</strong>trico base ortonl de<br />
a<strong>ut</strong>ovec. {<strong>ê</strong> i} (dir. ppales,) y<br />
a<strong>ut</strong>ovalores reales de modo que: S<br />
= 1<strong>ê</strong> 1<strong>ê</strong> 1+ 2<strong>ê</strong> 2<strong>ê</strong> 2+ 3<strong>ê</strong> 3<strong>ê</strong> 3<br />
clasificación como las formas<br />
cuadráticas x·S·x:<br />
–definidos pos. o neg.<br />
–semidefinidos pos. o ng.<br />
–indefinidos<br />
seg<strong>ú</strong>n los signos de los tres<br />
a<strong>ut</strong>ovalores reales de S.<br />
A es antisim. A t = A "v :<br />
v·A = A·v "u, v : u·A·v = <br />
v·A·u<br />
vector axial:<br />
A = antisim. vector axial<br />
/ A = ×.<br />
El vector axial kerA<br />
pol car. () = (2 +2 )<br />
– plano invariante {} ,<br />
asociado al par a<strong>ut</strong>ovals.<br />
complejos = ±i (sin rectas<br />
invariantes)<br />
– (0) = kerA = {}<br />
ej. Ejercicios*: 1) Si T es regular, probar que T t ·T es sim<strong>é</strong>trico y definido positivo.<br />
Pueden hacerse: PR1 n os . 6, 9, 14, 15<br />
g2) Rotaciones y simetrías<br />
Q es ortogonal Q t = Q 1 Q<br />
es un tensor de c. de b.<br />
ortonormales: {e i} = base<br />
ortonl. : {Q·e i} = nueva base<br />
ortonl.<br />
aplicación conforme:<br />
Q es ortogonal Q conserva el<br />
producto escalar: " u, v :<br />
(Q·u)·(Q·v) = u·v Qes apl.<br />
conforme, o sea conserva<br />
ángulos y módulos<br />
Q ortogonal detQ = ±1; detQ<br />
= 1 Q= rotación y detQ = 1<br />
Q= simetría (respecto de un<br />
plano) o prod. de rotación por<br />
simetría<br />
El grupo lineal (3;) y el<br />
grupo ortogonal (3;)<br />
Expresión intrínseca de una rotación R(e; ) :<br />
R =cos1 +(1cos)ee +sene× → formulario<br />
•Aplicaciones:<br />
•cálculo de componentes de una rotación: basta calcular las de cada sumando<br />
•ejemplo: rotaciones alrededor de los ejes cartesianos<br />
•cálculo del ángulo de giro: es tal que trR = 1+2cos.<br />
•cálculo del eje e : debe cumplir: e×= [R – Rt ]/(2sen); tambi<strong>é</strong>n: está en la dirección del<br />
a<strong>ut</strong>ovector del a<strong>ut</strong>ovalor = 1 que necesariamente tiene la rotación<br />
•Ejemplo: Se dan la rotación R1, alrededor del vector e1e2 y de ángulo 1 = /6, y la rotación<br />
R2, de eje e3 y ángulo 2 = /3, siendo {ei} una base ortonormal de V3. Se pide. Expresar R1 y R2 en la base cartesiana {ei}. Calcular el eje e y el ángulo de la rotación R = R1·R2. (sol.)<br />
•Las rotaciones forman un subgrupo, denotado (3;), del grupo ortogonal (3;). Este es a su<br />
vez es subgrupo del grupo de a<strong>ut</strong>omorfismos o tensores regulares, (3,)<br />
•Ejemplo: Demostrar que el producto de dos rotaciones es una nueva rotación: R1·R2 es un<br />
tensor ortogonal con determinante +1 es una rotación.<br />
Expresión intrínseca de una simetría respecto {e} :<br />
H(e) = 1 –2ee<br />
•Aplicaciones: cálculo de las componentes de H<br />
•Ejemplo 1: Compntes. canónicas de la Simetría respecto del plano vectorial {3x + 2y z = 0}<br />
•Ejemplo 2: El producto de dos simetrías respecto de dos planos 1 y 2 es una rotación:<br />
identificar el eje de la rotación en t<strong>é</strong>rminos de dichos planos.<br />
•PR1: ejercicios nº 8-3), 9, 12, 13-4)<br />
Curso 2012-13 5
Teoría de Campos - § 1.4 Tensores de<br />
2do. orden<br />
Expresión intrínseca de una rotación, R(e;)<br />
P<br />
v<br />
e<br />
Objetivo: Expresar R·v en t<strong>é</strong>rm. de v, y e :<br />
R·v =OP'=OM+MN+NP'<br />
y se tiene:<br />
OM= (v·e)e<br />
MN = |MP'| cos (MP/|MP|) =<br />
=cos [v (v·e)e]<br />
NP' = |MP'| sen (ev)/|ev|=<br />
= |v|sensen (ev)/(|v|sen)=sen (ev)<br />
De donde deducimos:<br />
R·v =[cos1 +(1cos) ee +sene]·v y de ahí, por la arbitrariedad de v:<br />
3v 2<br />
N<br />
M<br />
<br />
<br />
<br />
O<br />
( 2 = 3)<br />
w 2<br />
v 2<br />
R ·v<br />
P'<br />
R =cos 1 +(1 cos)ee +sen e<br />
Descomposición de la acción de un tensor T diagonalizable en 2: observar<br />
cómo construimos gráficamente la imagen de v con los a<strong>ut</strong>ovalores y<br />
a<strong>ut</strong>ovectores asociados.<br />
w 1<br />
Curso 2012-13 6<br />
v<br />
v 1<br />
T·v<br />
2v 1<br />
( 1 = 2)
Teoría de Campos - § 1.4 Tensores de<br />
2do. orden<br />
Subgrupos multiplicativos notables de (2) ( 3)<br />
(2)<br />
(3,) = (3,)<br />
(3,)<br />
(3,)<br />
•Ejemplo de pl. invariante sin rectas invariantes: Tensor de matriz canónica<br />
T·e<br />
2<br />
e<br />
3<br />
e 2<br />
e1<br />
a b<br />
0<br />
<br />
b a 0<br />
<br />
<br />
0 0 2<br />
T·e1<br />
Curso 2012-13 7<br />
a<br />
b