Derivada de un tensor (mecánica de medios continuos)

De testwiki
Ir a la navegación Ir a la búsqueda

Las derivadas de escalares, vectores y tensores de segundo orden con respecto a los tensores de segundo orden son de considerable utilidad en mecánica de medios continuos. Estas derivadas se utilizan en las teorías de elasticidad no lineal y plasticidad, particularmente en el diseño de algoritmos de simulación.[1]

La derivada direccional proporciona una forma sistemática de encontrar estas derivadas.[2]

Derivadas con respecto a vectores y tensores de segundo orden

A continuación se dan las definiciones de derivadas direccionales para diversas situaciones. Se supone que las funciones son lo suficientemente suaves como para poder tomar derivadas.

Derivadas de funciones escalares de vectores

Sea f(v) una función con valor real del vector v. Entonces, la derivada de f(v) con respecto a v (o en v) es el vector definido mediante el producto escalar con cualquier vector u

f𝐯𝐮=Df(𝐯)[𝐮]=[ddαf(𝐯+α𝐮)]α=0

para todos los vectores u. El producto escalar anterior produce un escalar, y si u es un vector unitario, da la derivada direccional de f en v, en la dirección u.

Propiedades:

  1. Si f(𝐯)=f1(𝐯)+f2(𝐯) entonces f𝐯𝐮=(f1𝐯+f2𝐯)𝐮
  2. Si f(𝐯)=f1(𝐯)f2(𝐯) entonces f𝐯𝐮=(f1𝐯𝐮)f2(𝐯)+f1(𝐯)(f2𝐯𝐮)
  3. Si f(𝐯)=f1(f2(𝐯)) entonces f𝐯𝐮=f1f2f2𝐯𝐮

Derivadas de funciones vectoriales de vectores

Sea f(v) una función con valor vectorial del vector v. Entonces la derivada de f(v) con respecto a v (o en v) es el tensor de segundo orden definido a través de su producto escalar con cualquier vector u

𝐟𝐯𝐮=D𝐟(𝐯)[𝐮]=[ddα𝐟(𝐯+α𝐮)]α=0

para todos los vectores u. El producto escalar anterior produce un vector, y si u es un vector unitario, da la derivada direccional de f en v, en la dirección u.

Propiedades:

  1. Si 𝐟(𝐯)=𝐟1(𝐯)+𝐟2(𝐯) entonces 𝐟𝐯𝐮=(𝐟1𝐯+𝐟2𝐯)𝐮
  2. Si 𝐟(𝐯)=𝐟1(𝐯)×𝐟2(𝐯) entonces 𝐟𝐯𝐮=(𝐟1𝐯𝐮)×𝐟2(𝐯)+𝐟1(𝐯)×(𝐟2𝐯𝐮)
  3. Si 𝐟(𝐯)=𝐟1(𝐟2(𝐯)) entonces 𝐟𝐯𝐮=𝐟1𝐟2(𝐟2𝐯𝐮)

Derivadas de funciones escalares de tensores de segundo orden

Sea f(𝑺) una función con valor real del tensor de segundo orden 𝑺. Entonces, la derivada de f(𝑺) con respecto a 𝑺 (o en 𝑺) en la dirección 𝑻 es el tensor de segundo orden definido como

f𝑺:𝑻=Df(𝑺)[𝑻]=[ddαf(𝑺+α𝑻)]α=0

para todos los tensores de segundo orden 𝑻.

Propiedades:

  1. Si f(𝑺)=f1(𝑺)+f2(𝑺) entonces f𝑺:𝑻=(f1𝑺+f2𝑺):𝑻
  2. Si f(𝑺)=f1(𝑺)f2(𝑺) entonces f𝑺:𝑻=(f1𝑺:𝑻)f2(𝑺)+f1(𝑺)(f2𝑺:𝑻)
  3. Si f(𝑺)=f1(f2(𝑺)) entonces f𝑺:𝑻=f1f2(f2𝑺:𝑻)

Derivadas de funciones tensoriales de tensores de segundo orden

Sea 𝑭(𝑺) una función tensorial de segundo orden del tensor de segundo orden 𝑺. Entonces la derivada de 𝑭(𝑺) con respecto a 𝑺 (o en 𝑺) en la dirección 𝑻 es el tensor de cuarto orden definido como

𝑭𝑺:𝑻=D𝑭(𝑺)[𝑻]=[ddα𝑭(𝑺+α𝑻)]α=0

para todos los tensores de segundo orden 𝑻.

Propiedades:

  1. Si 𝑭(𝑺)=𝑭1(𝑺)+𝑭2(𝑺) entonces 𝑭𝑺:𝑻=(𝑭1𝑺+𝑭2𝑺):𝑻
  2. Si 𝑭(𝑺)=𝑭1(𝑺)𝑭2(𝑺) entonces 𝑭𝑺:𝑻=(𝑭1𝑺:𝑻)𝑭2(𝑺)+𝑭1(𝑺)(𝑭2𝑺:𝑻)
  3. Si 𝑭(𝑺)=𝑭1(𝑭2(𝑺)) entonces 𝑭𝑺:𝑻=𝑭1𝑭2:(𝑭2𝑺:𝑻)
  4. Si f(𝑺)=f1(𝑭2(𝑺)) entonces f𝑺:𝑻=f1𝑭2:(𝑭2𝑺:𝑻)

Gradiente de un campo tensorial

El gradiente, 𝑻, de un campo tensorial 𝑻(𝐱) en la dirección de un vector constante arbitrario c se define como:

𝑻𝐜=limα0ddα𝑻(𝐱+α𝐜)

El gradiente de un campo tensorial de orden n es un campo tensorial de orden n+1.

Coordenadas cartesianas

Si 𝐞1,𝐞2,𝐞3 son los vectores base en un sistema de coordenadas cartesianas, con las coordenadas de los puntos indicadas por (x1,x2,x3), entonces el gradiente del campo tensorial 𝑻 viene dado por

𝑻=𝑻xi𝐞i

Plantilla:Demostración

Dado que los vectores de la base no varían en un sistema de coordenadas cartesiano, tenemos las siguientes relaciones para los gradientes de un campo escalar ϕ, un campo vectorial v' y un campo tensorial de segundo orden 𝑺.

ϕ=ϕxi𝐞i=ϕ,i𝐞i𝐯=(vj𝐞j)xi𝐞i=vjxi𝐞j𝐞i=vj,i𝐞j𝐞i𝑺=(Sjk𝐞j𝐞k)xi𝐞i=Sjkxi𝐞j𝐞k𝐞i=Sjk,i𝐞j𝐞k𝐞i

Coordenadas curvilíneas

Plantilla:AP

Si 𝐠1,𝐠2,𝐠3 son los vectores de una base contravariante en un sistema de coordenadas curvilíneas, con las coordenadas de los puntos indicadas por (ξ1,ξ2,ξ3), entonces el gradiente del campo tensorial 𝑻 viene dado por (consulte[3] para obtener una demostración).

𝑻=𝑻ξi𝐠i

De esta definición se obtienen las siguientes relaciones para los gradientes de un campo escalar ϕ, un campo vectorial v y un campo tensorial de segundo orden 𝑺.

ϕ=ϕξi𝐠i𝐯=(vj𝐠j)ξi𝐠i=(vjξi+vkΓikj)𝐠j𝐠i=(vjξivkΓijk)𝐠j𝐠i𝑺=(Sjk𝐠j𝐠k)ξi𝐠i=(SjkξiSlkΓijlSjlΓikl)𝐠j𝐠k𝐠i

donde los símbolos de Christoffel Γijk se definen usando

Γijk𝐠k=𝐠iξjΓijk=𝐠iξj𝐠k=𝐠i𝐠kξj

Coordenadas polares cilíndricas

En coordenadas cilíndricas, el gradiente viene dado por

ϕ=ϕr𝐞r+1rϕθ𝐞θ+ϕz𝐞z𝐯=vrr𝐞r𝐞r+1r(vrθvθ)𝐞r𝐞θ+vrz𝐞r𝐞z+vθr𝐞θ𝐞r+1r(vθθ+vr)𝐞θ𝐞θ+vθz𝐞θ𝐞z+vzr𝐞z𝐞r+1rvzθ𝐞z𝐞θ+vzz𝐞z𝐞z𝑺=Srrr𝐞r𝐞r𝐞r+Srrz𝐞r𝐞r𝐞z+1r[Srrθ(Sθr+Srθ)]𝐞r𝐞r𝐞θ+Srθr𝐞r𝐞θ𝐞r+Srθz𝐞r𝐞θ𝐞z+1r[Srθθ+(SrrSθθ)]𝐞r𝐞θ𝐞θ+Srzr𝐞r𝐞z𝐞r+Srzz𝐞r𝐞z𝐞z+1r[SrzθSθz]𝐞r𝐞z𝐞θ+Sθrr𝐞θ𝐞r𝐞r+Sθrz𝐞θ𝐞r𝐞z+1r[Sθrθ+(SrrSθθ)]𝐞θ𝐞r𝐞θ+Sθθr𝐞θ𝐞θ𝐞r+Sθθz𝐞θ𝐞θ𝐞z+1r[Sθθθ+(Srθ+Sθr)]𝐞θ𝐞θ𝐞θ+Sθzr𝐞θ𝐞z𝐞r+Sθzz𝐞θ𝐞z𝐞z+1r[Sθzθ+Srz]𝐞θ𝐞z𝐞θ+Szrr𝐞z𝐞r𝐞r+Szrz𝐞z𝐞r𝐞z+1r[SzrθSzθ]𝐞z𝐞r𝐞θ+Szθr𝐞z𝐞θ𝐞r+Szθz𝐞z𝐞θ𝐞z+1r[Szθθ+Szr]𝐞z𝐞θ𝐞θ+Szzr𝐞z𝐞z𝐞r+Szzz𝐞z𝐞z𝐞z+1rSzzθ𝐞z𝐞z𝐞θ

Divergencia de un campo tensorial

La divergencia de un campo tensorial 𝑻(𝐱) se define usando la relación recursiva

(𝑻)𝐜=(𝐜𝑻T);𝐯=tr(𝐯)

donde c es un vector constante arbitrario y v es un campo vectorial. Si 𝑻 es un campo tensorial de orden n > 1, entonces la divergencia del campo es un tensor de orden n- 1.

Coordenadas cartesianas

En un sistema de coordenadas cartesiano se tienen las siguientes relaciones para un campo vectorial v' y un campo tensorial de segundo orden 𝑺

𝐯=vixi=vi,i𝑺=Sikxi𝐞k=Sik,i𝐞k

donde con la notación tensorial indexada para derivadas parciales se utiliza en las expresiones situadas más a la derecha. Tenga en cuenta que

𝑺𝑺T.

Para un tensor simétrico de segundo orden, la divergencia también suele escribirse como[4]

𝑺=Skixi𝐞k=Ski,i𝐞k

La expresión anterior se utiliza a veces como definición de 𝑺 en forma de componente cartesiano (a menudo también escrito como div𝑺). Téngase en cuenta que dicha definición no es coherente con el resto de este artículo (consúltese la sección sobre coordenadas curvilíneas).

La diferencia surge de si la diferenciación se realiza respecto de las filas o columnas de 𝑺, y es convencional. Esto se demuestra con un ejemplo. En un sistema de coordenadas cartesiano, el tensor (matriz) de segundo orden 𝐒 es el gradiente de una función vectorial 𝐯

(𝐯)=(vi,j𝐞i𝐞j)=vi,ji𝐞i𝐞i𝐞j=(𝐯),j𝐞j=(𝐯)[(𝐯)T]=(vj,i𝐞i𝐞j)=vj,ii𝐞i𝐞i𝐞j=2vj𝐞j=2𝐯

La última ecuación es equivalente a la definición/interpretación alternativa[4]

()alt(𝐯)=()alt(vi,j𝐞i𝐞j)=vi,jj𝐞i𝐞j𝐞j=2vi𝐞i=2𝐯

Coordenadas curvilíneas

Plantilla:AP

En coordenadas curvilíneas, las divergencias de un campo vectorial v' y de un campo tensorial de segundo orden 𝑺 son

𝐯=(viξi+vkΓiki)𝑺=(SikξiSlkΓiilSilΓikl)𝐠k

Más generalmente,

𝑺=[SijqkΓkilSljΓkjlSil]gik𝐛j=[Sijqi+ΓiliSlj+ΓiljSil]𝐛j=[Sjiqi+ΓiliSjlΓijlSli]𝐛j=[SijqkΓiklSlj+ΓkljSil]gik𝐛j

Coordenadas polares cilíndricas

En coordenadas polares cilíndricas

𝐯=vrr+1r(vθθ+vr)+vzz𝑺=Srrr𝐞r+Srθr𝐞θ+Srzr𝐞z+1r[Sθrθ+(SrrSθθ)]𝐞r+1r[Sθθθ+(Srθ+Sθr)]𝐞θ+1r[Sθzθ+Srz]𝐞z+Szrz𝐞r+Szθz𝐞θ+Szzz𝐞z

Rotacional de un campo tensorial

El rotacional de un campo tensorial de orden n > 1 𝑻(𝐱) también se define usando la relación recursiva

(×𝑻)𝐜=×(𝐜𝑻);(×𝐯)𝐜=(𝐯×𝐜)

donde c es un vector constante arbitrario y v es un campo vectorial.

Rotacional de un campo tensorial (vectorial) de primer orden

Considérese un campo vectorial v y un vector constante arbitrario c. En notación indexada, el producto cruzado viene dado por

𝐯×𝐜=εijkvjck𝐞i

donde εijk es el símbolo de permutación, también conocido como símbolo de Levi-Civita. Entonces,

(𝐯×𝐜)=εijkvj,ick=(εijkvj,i𝐞k)𝐜=(×𝐯)𝐜

Por lo tanto,

×𝐯=εijkvj,i𝐞k

Rotacional de un campo tensorial de segundo orden

Para un tensor de segundo orden 𝑺

𝐜𝑺=cmSmj𝐞j

Por tanto, utilizando la definición de la curvatura de un campo tensorial de primer orden,

×(𝐜𝑺)=εijkcmSmj,i𝐞k=(εijkSmj,i𝐞k𝐞m)𝐜=(×𝑺)𝐜

Por lo tanto, se tiene que

×𝑺=εijkSmj,i𝐞k𝐞m

Identidades que involucran la curvatura de un campo tensorial

La identidad más comúnmente utilizada que involucra la curvatura de un campo tensorial, 𝑻, es

×(𝑻)=0

Esta identidad es válida para campos tensoriales de todos los órdenes. Para el caso importante de un tensor de segundo orden, 𝑺, esta identidad implica que

×(𝑺)=0Smi,jSmj,i=0

Derivada del determinante de un tensor de segundo orden

La derivada del determinante de un tensor de segundo orden 𝑨 viene dada por

𝑨det(𝑨)=det(𝑨)[𝑨1]T.

En términos ortonormales, las componentes de 𝑨 se pueden escribir como una matriz A. En ese caso, el lado derecho corresponde a los cofactores de la matriz.

Plantilla:Demostración

Derivadas de los invariantes de un tensor de segundo orden

Los principales invariantes de un tensor de segundo orden son

I1(𝑨)=tr𝑨I2(𝑨)=12[(tr𝑨)2tr𝑨2]I3(𝑨)=det(𝑨)

Las derivadas de estos tres invariantes con respecto a 𝑨 son

I1𝑨=1I2𝑨=I11𝑨TI3𝑨=det(𝑨)[𝑨1]T=I21𝑨T(I11𝑨T)=(𝑨2I1𝑨+I21)T

Plantilla:Demostración

Derivada del tensor de identidad de segundo orden

Sea 1 el tensor de identidad de segundo orden. Entonces la derivada de este tensor con respecto a un tensor de segundo orden 𝑨 viene dada por

1𝑨:𝑻=𝟢:𝑻=0

Esto se debe a que 1 es independiente de 𝑨.

Derivada de un tensor de segundo orden con respecto a sí mismo

Sea 𝑨 un tensor de segundo orden. Entonces

𝑨𝑨:𝑻=[α(𝑨+α𝑻)]α=0=𝑻=𝖨:𝑻

Por lo tanto,

𝑨𝑨=𝖨

Aquí 𝖨 es el tensor de identidad de cuarto orden. En notación indexada con respecto a una base ortonormal

𝖨=δikδjl𝐞i𝐞j𝐞k𝐞l

Este resultado implica que

𝑨T𝑨:𝑻=𝖨T:𝑻=𝑻T

donde

𝖨T=δjkδil𝐞i𝐞j𝐞k𝐞l

Por lo tanto, si el tensor 𝑨 es simétrico, entonces la derivada también es simétrica y se obtiene

𝑨𝑨=𝖨(s)=12(𝖨+𝖨T)

donde el tensor de identidad simétrico de cuarto orden es

𝖨(s)=12(δikδjl+δilδjk)𝐞i𝐞j𝐞k𝐞l

Derivada del inverso de un tensor de segundo orden

Sean 𝑨 y 𝑻 dos tensores de segundo orden, entonces

𝑨(𝑨1):𝑻=𝑨1𝑻𝑨1

En notación indexada con respecto a una base ortonormal

Aij1AklTkl=Aik1TklAlj1Aij1Akl=Aik1Alj1

También se tiene que

𝑨(𝑨T):𝑻=𝑨T𝑻T𝑨T

En notación indexada

Aji1AklTkl=Ajk1TlkAli1Aji1Akl=Ali1Ajk1

Si el tensor 𝑨 es simétrico entonces

Aij1Akl=12(Aik1Ajl1+Ail1Ajk1)

Plantilla:Demostración

Integración por partes

Dominio Ω, su frontera Γ y el vector normal unitario exterior 𝐧

Otra operación importante relacionada con las derivadas tensoriales en la mecánica continua es la integración por partes. La fórmula de integración por partes se puede escribir como

Ω𝑭𝑮dΩ=Γ𝐧(𝑭𝑮)dΓΩ𝑮𝑭dΩ

donde 𝑭 y 𝑮 son campos tensoriales diferenciables de orden arbitrario, 𝐧 es la unidad normal hacia afuera con respecto al dominio sobre el cual se definen los campos tensoriales, representa un operador del producto tensorial generalizado y es un operador de gradiente generalizado. Cuando 𝑭 es igual al tensor de identidad, se obtiene el teorema de la divergencia

Ω𝑮dΩ=Γ𝐧𝑮dΓ.

Se puede expresar la fórmula de integración por partes en coordenadas cartesianas con notación indexada como

ΩFijk....Glmn...,pdΩ=ΓnpFijk...Glmn...dΓΩGlmn...Fijk...,pdΩ.

Para el caso especial donde la operación del producto tensorial es una contracción de un índice y la operación del gradiente es una divergencia, y tanto 𝑭 como 𝑮 son tensores de segundo orden, se tiene que

Ω𝑭(𝑮)dΩ=Γ𝐧(𝑮𝑭T)dΓΩ(𝑭):𝑮TdΩ.

En notación indexada,

ΩFijGpj,pdΩ=ΓnpFijGpjdΓΩGpjFij,pdΩ.

Véase también

Referencias

Plantilla:Listaref

Plantilla:Control de autoridades

  1. J. C. Simo and T. J. R. Hughes, 1998, Computational Inelasticity, Springer
  2. J. E. Marsden and T. J. R. Hughes, 2000, Mathematical Foundations of Elasticity, Dover.
  3. R. W. Ogden, 2000, Nonlinear Elastic Deformations, Dover.
  4. 4,0 4,1 Plantilla:Cite book