Origen de la propiedad de simetría mayor del tensor de elasticidad

En la teoría de la elasticidad lineal, el tensor de tensión σ está relacionado con el tensor de deformación ϵ a través del tensor elástico C . Específicamente

σ i j = C i j k yo ϵ k yo

Porque σ y ϵ ambos son tensores simétricos de segundo rango, σ i j = σ j i y ϵ k yo = ϵ yo k de modo que C tiene las llamadas "simetrías menores":

C i j k yo = C j i k yo = C i j yo k

Puedo convencerme de eso, pero con lo que tengo problemas es con la "simetría principal":

C i j k yo = C k yo i j

Esto supuestamente proviene de la simetría de la densidad de energía de la tensión. ψ ( ϵ ) ,

ψ = 1 2 C i j k yo ϵ i j ϵ k yo

o alguna propiedad de sus segundas derivadas parciales

C i j k yo = 2 ψ ϵ i j ϵ k yo = 2 ψ ϵ k yo ϵ i j = C k yo j i

Para mí, parece que solo cambia o vuelve a etiquetar los índices en lugar de la transposición. ¿Alguna idea de lo que me estoy perdiendo?

** AMPLIADO DESPUÉS DE LA PUBLICACIÓN DE Fénix **

Suponiendo un sólido 2D en tamaño reducido ( I j ) notación la densidad de deformación es

ψ = ϵ I C I j ϵ j

ψ = [ mi F ] [ C 11 C 12 C 21 C 22 ] [ a b ]

ψ = mi [ C 11 a + C 12 b ] + F [ C 21 a + C 22 b ]

y las segundas derivadas son

2 ψ mi b = C 12
y
2 ψ b mi = C 12

Incluso en la situación especial de mi = a y F = b esto todavía sale como

2 ψ a b = C 12 + C 21
y
2 ψ b a = C 12 + C 21

es decir, nada allí requiere C ser simétrico.

Su edición no está relacionada con la respuesta de Phoenix, en 2D, el índice múltiple es I = {11}, {12} o {22}, por lo que debe tener una matriz de 3x3. El multiíndice es irrelevante para la idea, que es que si arreglas i , j , k , y yo entonces C i j k yo = 2 ψ ϵ i j ϵ k yo = 2 ψ ϵ k yo ϵ i j = C k yo i j .
Tuve la misma pregunta y entendí la respuesta muy fácilmente de este video, la respuesta está en la parte: Reducción 2. youtube.com/watch?v=IhHZ3V4uBsI
@doobs ¡De hecho, la prueba anterior basada en la densidad de energía de tensión es evidentemente incorrecta! En su definición, la posible parte antisimétrica de C (que evidentemente se admite en la primera fórmula que escribiste) se cancela. Entonces no es posible recuperar C de ψ ! En su lugar, se necesita algún argumento físico.

Respuestas (3)

Debes comenzar con la densidad de energía de deformación. ψ , luego define:

C i j k yo = 2 ψ ϵ i j ϵ k yo ,

y luego definir

σ i j = C i j k yo ϵ k yo

El resto de mi respuesta se tratará de explicar por qué tienes que hacerlo de esa manera. En primer lugar, es físico, realmente hay energía asociada con la tensión, y si no la hubiera, no habría estrés. En segundo lugar, es lineal exactamente porque estamos considerando los cambios de energía debidos a pequeñas deformaciones.

Pero volvamos a las simetrías menores. Nosotros necesitamos C i j k yo = C j i k yo porque de otra manera σ i j σ j i (y luego obtenemos velocidades angulares arbitrariamente grandes, por lo tanto no físicas, para regiones cada vez más pequeñas). Pero la otra simetría menor no es necesaria. Si alguien te entregó un tensor aleatorio de rango cuatro, llamémoslo B i j k yo , y lo llamamos tensor elástico y no tenía la segunda simetría menor, puedes definir C i j k yo = ( B i j k yo + B i j yo k ) / 2 D i j k yo = ( B i j k yo B i j yo k ) / 2 y luego B = C + D pero cuando D se contrae con un tensor simétrico de rango dos (como el tensor de deformación) da cero. Entonces, la parte del tensor de rango cuatro sin esa segunda simetría menor simplemente no contribuye, como actor no hace nada (cuando crees que todo lo que hace la elasticidad es darte estrés por tensión). Entonces, también podría asumir que su tensor tiene ambas simetrías menores porque actúa como si tuviera la segunda ( B y C actuar de la misma en tensores simétricos) y tiene que tener el primero.

¿Lo mencioné para ser pedante? No, lo mencioné porque sucede lo mismo si contraes el tensor de elasticidad con una combinación simétrica de rango cuatro del tensor de deformación. La parte del tensor de elasticidad sin mayor simetría no contribuye a la densidad de energía de deformación. Entonces, un tensor aleatorio necesita que la primera simetría menor sea física. Pero también podría asumir que tiene la segunda simetría menor, ya que no afecta la relación tensión-deformación. Y también podría asumir que tiene la mayor simetría porque la parte que no la tiene no contribuirá a la densidad de energía de tensión.

Pero es la densidad de energía de tensión lo que es físico, y cómo cambia es lo que es la elasticidad. Entonces, en realidad no está derivando estas simetrías tanto como diciendo que solo las simétricas generan las cosas físicas que desea, energía cuando se le aplica tensión. Y una derivación real debería comenzar con la densidad de energía de deformación y la deformación, y luego simplemente definir la elasticidad a partir de eso.

tienes razón en que la parte antisimétrica no juega un papel en la densidad de energía, pero lo que le impide jugar un papel en la relación entre σ y ϵ ? No hay razón a priori por la que la parte antisimétrica deba desaparecer también de allí.
@ Phoenix87 La densidad de energía es lo físico principal. El estrés se define correctamente en términos de energía, no al revés. Una vez que define la densidad de energía (en función de la deformación), obtiene la elasticidad en términos de la densidad de energía y luego la tensión en términos de la elasticidad (y la deformación). Así que digo que la razón a priori es que la energía es fundamental y que la elasticidad no es ni más ni menos que cómo la energía depende de la tensión. El hecho de que el estrés dependa de la tensión es, para mí, simplemente un hecho derivado.
El comentario sobre la segunda simetría menor es muy útil y explica lo que a menudo se oculta: en todas partes se lee que la simetría "viene de la simetría de ϵ "; como demuestras, esta es una elección (no restrictiva). ¡Gracias!
@Timaeus (Refiriéndose a las dos fórmulas iniciales) ¿Por qué el tensor C en la última identidad es el mismo que en la primera? Quiero decir, esto último es un hecho físico, no una definición. Con su enfoque, está imponiendo desde cero que C es simétrico (simetría mayor) en una ley física sin motivaciones, en mi opinión.
Debería ser más honesto declarar que la simetría mayor es un hecho físico. @ Poenix87 tiene razón en mi opinión.
Sin ofender, pero en mi opinión, no tiene sentido ocultar las dificultades físicas detrás del formalismo matemático. Las matemáticas son una mera herramienta. Todos los hechos físicos deben ser declarados. La segunda identidad es una ley física ya que en ambos lados, σ y ϵ pueden medirse, en principio, por separado. C puede resultar simétrico o no según la física no según una definición.

Desde ϵ es un tensor simétrico, tiene 6 componentes independientes que lo determinan. Por lo tanto, use un índice múltiple I { ( i , j ) | 1 i j 3 } para denotarlos. La densidad de energía de deformación entonces se convierte en (quizás hay que tener cuidado con los términos "diagonales" aquí para obtener los coeficientes correctos)

ψ = C I j ϵ I ϵ j
donde la suma está destinada a múltiples índices repetidos. Está claro que, dado que supone que el orden de las segundas derivadas es el mismo ϵ I s no importa, se deduce que
C I j = C j I ,
para cualquier par de multiíndices I y j . Entonces, al final, este es el mismo argumento que muestra que una matriz hessiana es simétrica si el orden de las segundas derivadas no importa.

Supongo que esto supone que ϵ I y ϵ j están en la misma base, mientras que tuve la sensación de que podría haber una especie de permutación o rotación entre sus ejes de referencia para que C ! = C T en general.
Cuando calculo las segundas derivadas, simplemente sale como C I j = C I j independientemente del orden de las derivadas parciales.
Echa un vistazo a en.wikipedia.org/wiki/Hooke%27s_law#Anisotropic_materials . El argumento que necesitas está ahí.

Lo siguiente es básicamente lo que dice Ashcroft/Mermin al respecto.

La idea es la siguiente:

en aproximación armónica un desplazamiento relativo tu da como resultado una energía

tu = 1 4 ( tu ( R ) tu ( R ) ) D ( tu ( R ) tu ( R ) )

el tensor D ya tiene simetría natural D m v = D v m , lo que significa que la interacción es simétrica. aproximando

tu ( R ) tu ( R ) + ( R R ) tu

uno obtiene

tu = 1 2 tu m X σ R σ D m v R τ tu v X τ = 1 2 tu m X σ mi σ m τ v tu v X τ

Entonces se afirma que la energía armónica no debe verse afectada por una rotación. Esto requiere que la energía dependa de las derivadas de la forma

ϵ σ m = 1 2 ( tu m X σ + tu σ X m )

Esto es algo así como un invariante de Lifschitz, supongo. No debería ser demasiado difícil probar eso.

Entonces se afirma que la energía se puede escribir en la forma

22.78) tu = 1 2 ϵ σ m C σ m τ v ϵ τ v

con

22.79) C σ m τ v = 1 8 ( mi σ m τ v + mi m σ τ v + mi σ m v τ + mi m σ v τ )

Mirando la definición de mi σ m τ v esto proporciona la simetría de interés. Entonces se origina en el requisito de que la energía sea invariable con la rotación.

Sin embargo, en este momento tengo algunas dificultades con el último paso (que es el paso 22.78 a 22.79 en Ashcroft/Mermin, de ahí los números) ya que creo que produce algunos derivados que no están presentes en la forma original. Si alguien puede explicar fácilmente por qué el formulario 22.78 requiere C para ser de la forma 22.79, me encantaría que me lo explicaran.