Ley de Hooke y tasas de estrés objetivas

A menudo, en artículos que presentan métodos de simulación Lagrangianos actualizados para dinámica de sólidos, se presenta el siguiente procedimiento para actualizar el tensor de tensión (Cauchy):

En primer lugar, el tensor de tensiones de Cauchy se divide en una parte hidrostática y otra desviadora:

σ i j = pag d i j + S i j

La presión se encuentra usando una ecuación de estado. A menudo se utiliza el siguiente enfoque isotérmico:

pag = C 0 2 ( ρ ρ 0 )

C 0 2 siendo la velocidad adiabática del sonido, ρ la densidad y ρ 0 la densidad de referencia. Entonces, se afirma que se asume la ley de Hooke y la parte desviadora del tensor de tensiones evoluciona de la siguiente manera:

d S i j d t = 2 m ( ϵ ˙ i j 1 3 d i j ϵ ˙ k k ) + S i j Ω j k + Ω i k S k j

dónde m es el módulo de cortante,

Ω i j = 1 2 ( d v i d X j d v j d X i )
es el tensor de espín y
ϵ ˙ i j = 1 2 ( d v i d X j + d v j d X i )
es el tensor de tasa de deformación. Ahora, desde:

σ ˙ = σ Δ j + σ Ω + Ω σ

dónde σ Δ j es la tasa de Jaumann que sostiene que:

σ Δ j i j = 2 m ( ϵ ˙ i j 1 3 d i j ϵ ˙ k k )

Ahora vamos a mis preguntas:

  1. ¿Cómo se llega a la ecuación anterior para la tasa de Jaumman? O, en particular, ¿cómo la suposición de la ley de Hooke produce esa ecuación para la tasa de Jaumann?

  2. ¿Esa ecuación para la tasa de Jaumann también es válida para otras tasas de estrés objetivas? Por ejemplo, para la tasa de Truesdell , dando una actualización de estrés de la siguiente manera:

d S d t = 2 m ( ϵ ˙ 1 3 1 T r ( ϵ ˙ ) ) T r ( ϵ ˙ ) S + ϵ ˙ S + S ϵ ˙ T

( T r ( . ) siendo la huella)

Respuestas (1)

Probablemente sepa que, si bien las tensiones de Cauchy son objetivas, su tasa de tensión (derivada material) no lo es. Si q ( t ) es un tensor ortogonal que representa un cambio de marco, la tensión es el nuevo marco es

T = q T q T

Sin embargo, si toma derivados materiales en ambos lados, tiene,

T ˙ = q ˙ T q T + q T ˙ q T + q T q ˙ T ( 1 )

Los términos 1 y 3 anteriores son "extra" -> implican que tomar ingenuamente el material derivado del estrés T no le dará una tasa objetiva.

Surge un problema similar porque la parte de "giro" del tensor de gradiente de velocidad L , (lo llamaré W en vez de Ω ) no es objetivo. es decir, L = W + D , pero mientras D es objetivo, W no es. A partir de los primeros principios, no es difícil demostrar que (agregaré este bit si lo desea)

W = q W q T + q ˙ q T ( 2 )

También porque q es ortogonal ( q T = q 1 ) tienes el lema

q ˙ q T = q q ˙ T

De hecho, podemos usar este segundo problema para construir una solución al primero. eliminaremos q ˙ , q ˙ T de (1) a su vez. Multiplique a la derecha (2) por Q y reorganice; tienes

q ˙ = W q q W
Transponga ambos lados y tendrá una expresión para q ˙ T .
q ˙ T = q T W + W q T
El último resultado utiliza el hecho de que el tensor de espín W es sesgadamente simétrica.

Sustituyamos estas expresiones por q ˙ , q ˙ T de vuelta en (1).

T ˙ = ( W q q W ) T q T + q T ˙ q T + q T ( q T W + W q T ) = W q T q T q W T q T + q T ˙ q T q T q T W + q T W q T ) = W T q W T q T + q T ˙ q T T W + q T W q T (al reorganizar) T ˙ + T W W T = q ( T ˙ + T W W T ) q T

Observe que el lado derecho tiene la forma q R q T , y el lado izquierdo es precisamente la cantidad entre paréntesis en el marco de referencia con estrella. En otras palabras, hemos construido una tasa j = T ˙ + T W W T tal que j ˙ = q j q T , que es objetivo por definición. Ahora puede especializar fácilmente sus ecuaciones en elasticidad.

Lo importante a recordar es que la tasa de Jaumann es simplemente una de un número infinito de tasas de estrés objetivas. En general, las derivadas de Lie de campos tensoriales objetivos son objetivas.

Tenga en cuenta @Matthias: su expresión que relaciona la tasa de Jaumann y σ ˙ tiene un problema de señal.
Gracias por tu extensa respuesta. Su derivación de la tasa de Jaumann es clara y fácil de seguir. Sin embargo, todavía lucho por especializar "fácilmente" las ecuaciones a la elasticidad. ¿Podría comentar mis preguntas 1 y 2 en concreto? Sería muy apreciado. (No puedo votar por falta de reputación)
@Matthias: puede ver la siguiente página (desplácese hacia abajo hasta 3.9.3): solidmechanics.org/text/Chapter3_9/Chapter3_9.htm Tenga en cuenta que hay una aproximación involucrada (es decir, que las tensiones son mucho más pequeñas que los módulos elásticos)
Además, siempre tenga cuidado de usar cantidades conjugadas de trabajo en sus relaciones tensión-deformación.