Convertir la velocidad de impacto en magnitud de presión

En las pruebas de seguridad y estabilidad de explosivos, se suele utilizar una prueba de caída para determinar la sensibilidad al impacto. En la prueba, se deja caer un impactador de masa conocida. La altura inicial varía a lo largo del experimento. El estado de detonación se registra en el material y se utilizan diversas técnicas para convertir las alturas en métricas utilizadas para caracterizar la estabilidad; por ejemplo, una métrica de uso común es h 50 o la altura del impactador que genera una detonación el 50% de las veces.

El material explosivo que se prueba suele tener forma de varilla o disco cuyo radio es considerablemente más pequeño que el impactador y el yunque. El impacto genera una onda de choque a través de la muestra de material.

Como estoy usando un código Euleriano para simular esto, necesito saber la presión de la onda de choque resultante en el material. ¿Hay alguna manera de determinar esta presión?

Supongamos que conozco todas las propiedades del material (densidad, velocidad del sonido, módulo de Young, etc.), que el impactador es perfectamente rígido y que el impacto está dentro de los límites elásticos de la muestra de material.

El único enfoque que se me ocurre sería tratar el problema como una barra 1D y asignarle una velocidad inicial (digamos en el X dirección) a la varilla de modo que el momento sea el mismo que el del impactador utilizado en el experimento. Luego imponga que la velocidad sea cero en el borde izquierdo de la barra. Esto requeriría algo de trabajo numérico y es menos que ideal. ¿Quizás existe una relación analítica o empírica entre el impulso del impactador y la magnitud de la onda de presión resultante?

Respuestas (1)

Antecedentes
Existe una buena referencia 1 sobre la física del sonido/ondas de choque en sólidos (consulte el Capítulo XI). Encontré lo siguiente (en la página 688) muy interesante y relevante para su pregunta:

En un sólido o líquido, una onda de choque con una fuerza de incluso cien mil atmósferas se considera débil. Tal onda difiere poco de una onda acústica: viaja con una velocidad cercana a la velocidad del sonido, comprime el material solo en un pequeño porcentaje o quizás del orden del diez por ciento, e imparte una velocidad al material detrás del frente que es del orden de una décima parte de la velocidad de la onda misma... entonces una onda de choque fuerte para medios condensados ​​es aquella cuya presión no es inferior a decenas o cientos de millones de atmósferas.

Definamos PAGS como la presión y ε como la energía interna de un material sólido. Estos se pueden dividir en dos partes: un elástico (subíndice C ) y parte térmica. PAGS C y ε C dependen únicamente de la densidad del material, ρ , o el volumen específico, V = 1 / ρ . Estos son iguales a la presión total y la energía interna específica en el cero absoluto o T = 0   k . Supongamos que el volumen específico en T = 0 y PAGS = 0 es dado por V o C , que es solo ~ 1-2% más pequeño que el volumen específico en STP , V o , para la mayoría de los metales.

La curva de energía potencial, o curva que define ε C , es cualitativamente similar a la curva de energía potencial que describe la interacción entre dos átomos en función de la distancia intranuclear, Δ X norte . Cuando V > V o C , las fuerzas de atracción dominan pero disminuyen rápidamente a medida que aumentan las distancias intranucleares (p. ej., a medida que T aumenta). En otras palabras, cuando los átomos se alejan más ε C aumentará asintóticamente a algún valor tu , que es aproximadamente la energía de enlace de los átomos en el cuerpo. De este modo, tu representa la energía requerida para eliminar todos los átomos del objeto hasta el infinito, que es aproximadamente igual al calor de vaporización del material (escribí algunos detalles más sobre el calor de vaporización y proporcioné varios enlaces útiles en esta respuesta ). Por ejemplo, el calor de atomización (similar a la vaporización) para el hierro es aproximadamente 415 k j metro o yo 1 o ~4,3 eV/átomo. De este modo, ε C ( V ) tu como Δ X norte →∼ 2 .

Por el contrario, las fuerzas repulsivas dominan si V < V o C . Podemos definir esto cuantitativamente considerando que el trabajo realizado al comprimir el material será igual al aumento de energía interna. En otras palabras:

PAGS C = ( d ε C d V ) T = 0
lo que equivale a decir que es la ecuación isotérmica / isoentrópica para la compresión en frío. El signo negativo muestra que si se aplicara una fuerza de tracción al cuerpo, las fuerzas de unión entre los átomos actuarían como una fuerza restauradora. la pendiente de la PAGS C curva en PAGS = 0 (o 1 atm) define la compresibilidad del material en condiciones normales (es decir, T = T o 300   k ). Esto viene dado por:
k o = 1 V o ( V PAGS ) T o
Tenga en cuenta que la pendiente de k o define la velocidad de las ondas elásticas dentro del objeto. Así, definamos la velocidad del sonido en el sólido como esta velocidad, dada por:
C o = V 2 ( PAGS V ) S
donde el subíndice S indica una derivada isoentrópica y la derivada parcial será negativa para evitar velocidades imaginarias del sonido.

Aproximación simple
de orden cero Mi suposición instintiva es que el enfoque más simple, dado que asume relaciones de colisión elásticas, es simplemente aproximar el Δ ε C por la energía cinética final de su objeto impactante, asumiendo que el impactado (?) No se mueve después del impacto.

Aproximación de primer orden
[Lo siguiente proviene del Capítulo XI, Secciones 3.14-3.16 en la Referencia 1]

A continuación consideraremos los efectos en una barra cilíndrica (usada por simetría y simplicidad).

Para pequeñas deformaciones, el cambio relativo en volumen, Δ V / V , es dado por:

Δ V V = k   PAGS = PAGS k
dónde k = 1 / k es el módulo volumétrico .

Definamos C 1 como la velocidad de una onda de compresión en el material debido a la aplicación de una presión constante, PAGS , aplicado a un extremo de la varilla en algún momento inicial. El material entre el frente de onda y el extremo de la barra se contrae a una velocidad constante, tu . Bajo estas condiciones, podemos usar la ley de Hooke y demostrar que para cargas y deformaciones pequeñas tenemos:

tu C 1 = PAGS mi
dónde mi es el módulo de Young . Después de algún tiempo, t , la masa de material abarcada por la onda adquirirá un impulso ρ   C 1   t   tu , que debe ser igual a PAGS   t de la ley de Newton, que nos da:
PAGS = ρ   tu   C 1

Respuesta más detallada
Desafortunadamente, no tengo tiempo para realizar la derivación completa, pero sugiero el Capítulo XI en la Referencia 1 y uso la Referencia 2 para obtener información de apoyo. Zel'dovich y Raizer básicamente pasan todo el Capítulo XI discutiendo este tema y analizan todos los matices que se aplicarían a su problema (por ejemplo, la onda de presión inducida por la compresión de la onda de choque). Supongo que gran parte necesitará un análisis numérico, pero hay puntos de partida y aproximaciones que son analíticas que probablemente le ahorrarán mucho tiempo.

Referencias

  1. Zel'dovich, Ya.B. y Yu.P. Raizer (2002) Física de ondas de choque y fenómenos hidrodinámicos de alta temperatura , ed. por WD Hayes y RF Probstein, Mineola, NY, Dover Publications, inc., The Dover Edition; ISBN-13: 978-0486420028.
  2. Whitham, GB (1999), Ondas lineales y no lineales , Nueva York, NY: John Wiley & Sons, Inc.; ISBN: 0-471-35942-4.