Criterio de estabilidad para ondas en sólidos anisotrópicos

8

Las ecuaciones de movimiento para un sólido elástico están dadas por

σ+f=ρu¨σ=Cεε=12(u+[u]T)

o en notación de índice

σij,j+fi=ρui¨σij=Cijklεklε=12(ui,j+uj,i)

u es el vector de desplazamiento,f es la fuerza del cuerpo (término fuente),σ es el tensor de tensión,ε es el tensor de tensión yC es el tensor de rigidez. En el caso de los sólidos isotrópicos, el tensor de rigidez se escribe en términos dedos constantes diferentes, para dominios no acotados, la ecuación admite dos tipos de ondas que están desacopladas y el criterio de estabilidad viene dado por el peor de los dos casos diferentes (es decir, , el de mayor velocidad).

Para los materiales isotrópicos transversales hay 5 parámetros independientes que definen el tensor y 3 tipos de ondas (2 de ellas están acopladas). En el caso más general, el número de parámetros es 21 y la onda está acoplada.

Pregunta: ¿Cómo encuentra el criterio de estabilidad en un algoritmo de marcha temporal para el caso general?

nicoguaro
fuente
1
¿Cómo define "estabilidad" y qué tipo de criterios está buscando?
Wolfgang Bangerth
Estoy buscando estabilidad en un esquema numérico explícito. El equivalente de la condición CFL.
nicoguaro

Respuestas:

6

Las ecuaciones de onda como esta pueden reescribirse como un sistema hiperbólico de leyes de conservación de primer orden:

qt+F(q)=0.

F

dFdF

Una buena referencia para esta teoría es el texto de LeVeque sobre métodos de volumen finito . La elasticidad se trata en detalle en el Capítulo 22.

Se aplican todas las advertencias habituales sobre la condición de CFL: es una condición necesaria pero no suficiente para la estabilidad. Pero la condición suficiente para la estabilidad de una discretización dada generalmente viene dada por la condición CFL multiplicada por alguna constante. Para encontrar su criterio de estabilidad, necesita conocer tanto la velocidad máxima de la onda (según las ecuaciones que está resolviendo) como esa constante (según la discretización que está utilizando).

David Ketcheson
fuente
ABC
2

[ρc2δijCijklnjnl][uk]=0
ρcδnjj

ρc2

Por lo tanto, en el caso anisotrópico tridimensional, todavía tenemos tres velocidades de fase diferentes para una dirección particular de propagación, la mayor de las cuales debe usarse para el análisis de CFL, de manera similar a la forma en que se usa la velocidad longitudinal en un problema isotrópico

DanielRch
fuente
1
maxθ,ϕmaxiλi(θ,ϕ)
λi
1

Expandiré la respuesta proporcionada por @DavidKetcheson. Primero, las ecuaciones se reescriben como un sistema hiperbólico de leyes de conservación de primer orden:

qt+F(q)=0

o

qt+Aqx+Bqy+Cqz=0

Donde es un vector de estado formado con los componentes del tensor de tensión y componentes del vector de velocidad .q(σ11,σ22,σ33,σ12,σ23,σ13)(u,v,w)

q=(σ11σ22σ33σ12σ23σ13uvw),

A=(000000c11c16c15000000c12c26c25000000c13c36c35000000c14c46c45000000c15c56c55000000c16c66c561ρ000000000001ρ00000000001ρ000),

B=(000000c16c12c14000000c26c22c24000000c36c23c34000000c46c24c44000000c56c25c45000000c66c26c460001ρ0000001ρ000000000001ρ0000),

C=(000000c15c14c13000000c25c24c23000000c35c34c33000000c45c44c34000000c55c45c35000000c56c46c36000001ρ00000001ρ0000001ρ000000).

Para calcular las velocidades del problema (como se describió anteriormente) necesitamos formar la matriz , donde es un vector unitario y determina la dirección de propagación. Para encontrar la condición CFL es necesario resolverA^(n1,n2,n3)=n1A+n2B+n3Cn=(n1,n2,n3)

max(θ,ϕ)maxiγi(θ,ϕ)

donde son ángulos esféricos y son los valores propios de la matriz .(θ,ϕ)γiA^(θ,ϕ)

En base a esto, y la respuesta proporcionada por @DavidKetcheson, es más sencillo calcular los valores propios de la ecuación de Christoffel y resolver el problema de optimización

max(θ,ϕ)maxiλi(θ,ϕ)

con valores propios de la ecuación de Christoffel. Y la velocidad es solo . c = λic=λi/ρ

nicoguaro
fuente