Tensores de esfuerzo macroscópico, esfuerzo de pareja y flujo derivados mediante equivalencia energética de volúmenes representativos finitos heterogéneos continuos y discretos microscópicos

Información Básica

• ID del Artículo: 2412.12616
• Título: Macroscopic stress, couple stress and flux tensors derived through energetic equivalence from microscopic continuous and discrete heterogeneous finite representative volumes
• Autores: Jan Eliáš (Universidad Tecnológica de Brno), Gianluca Cusatis (Universidad Northwestern)
• Clasificación: cs.CE (Ingeniería Computacional), cond-mat.other
• Fecha de Presentación: Diciembre de 2024
• Enlace del Artículo: https://arxiv.org/abs/2412.12616

Resumen

Este artículo propone un método de derivación riguroso para evaluar tensores de esfuerzo macroscópico, tensores de esfuerzo de pareja y vectores de flujo en sistemas heterogéneos continuos y discretos con desplazamientos y rotaciones independientes. A diferencia de la homogeneización clásica por expansión asintótica, este trabajo considera volúmenes de elementos representativos (RVE) de tamaño finito. Primero se derivan las cantidades macroscópicas para un continuo de Cosserat heterogéneo, luego se discretizan las ecuaciones continuas para proporcionar cantidades macroscópicas en sistemas heterogéneos discretos, y finalmente se rederivan las expresiones para sistemas discretos directamente a partir de propiedades discretas.

Antecedentes de Investigación y Motivación

1. Problemas a Resolver

Este artículo aborda principalmente cómo extraer con precisión propiedades mecánicas y de transporte macroscópicas de sistemas de materiales heterogéneos a escala microscópica, en particular:

• Cómo derivar tensores de esfuerzo macroscópico a partir de volúmenes de elementos representativos de tamaño finito
• Cómo manejar tensores de esfuerzo de pareja en sistemas con grados de libertad de rotación independientes
• Cómo establecer relaciones de equivalencia energética entre escalas microscópicas y macroscópicas

2. Importancia del Problema

Muchos materiales de ingeniería (como hormigón, cerámica, materiales granulares, etc.) exhiben heterogeneidad significativa a escala mesoscópica o microscópica. El modelado multiesacala preciso es importante para:

• Predecir el comportamiento mecánico general de los materiales
• Desarrollar modelos numéricos confiables
• Establecer relaciones constitutivas fenomenológicas

3. Limitaciones de Métodos Existentes

• Homogeneización por expansión asintótica: Asume que el volumen del elemento representativo tiende a cero, ignorando efectos de tamaño finito
• Fórmula clásica de Love-Weber: Principalmente aplicable a tensores de esfuerzo, con tratamiento incompleto de tensores de esfuerzo de pareja
• Tratamiento de condiciones de frontera: Los métodos existentes son insuficientes en el manejo de efectos de frontera y selección de puntos de referencia

4. Motivación de la Investigación

Este trabajo tiene como objetivo establecer un marco teórico más completo y riguroso que pueda:

• Manejar volúmenes de elementos representativos de tamaño finito
• Considerar simultáneamente formulaciones de fuerzas internas y externas
• Aclarar explícitamente la influencia de la ubicación del punto de referencia en el tensor de esfuerzo de pareja

Contribuciones Principales

1. Establecimiento de un marco de derivación riguroso basado en equivalencia de trabajo virtual, aplicable a sistemas heterogéneos continuos y discretos
2. Derivación de expresiones de cantidades macroscópicas considerando tamaño finito de RVE, incluyendo tensor de esfuerzo, tensor de esfuerzo de pareja y vector de flujo
3. Proposición de dos formas variantes: formulación basada en fuerzas internas y basada en fuerzas externas
4. Descubrimiento de términos adicionales en el tensor de esfuerzo de pareja, explicando su dependencia del sistema de referencia y la ubicación del punto macroscópico
5. Provisión de verificación completa, incluyendo comparación de soluciones analíticas y modelos numéricos bajo condiciones estacionarias y transitorias

Explicación Detallada del Método

Definición de Tareas

La tarea central de este trabajo es establecer relaciones de equivalencia entre materiales heterogéneos microscópicos y cuerpos continuos homogéneos macroscópicos, incluyendo específicamente:

• Entrada: Campo de esfuerzo microscópico σ^mic, campo de esfuerzo de pareja μ^mic, campo de flujo a^mic
• Salida: Tensor de esfuerzo macroscópico σ^mac, tensor de esfuerzo de pareja μ^mac, vector de flujo a^mac
• Restricciones: Satisfacer el principio de equivalencia de trabajo virtual y ecuaciones de equilibrio

Marco Teórico

1. Fundamentos del Continuo de Cosserat

Las variables principales del continuo de Cosserat incluyen:

• Vector de desplazamiento u y vector de rotación independiente θ
• Tensor de deformación γ = ε - E·θ^m y tensor de curvatura κ = ∇⊗θ
• Tensor de esfuerzo σ (no simétrico) y tensor de esfuerzo de pareja μ

Las ecuaciones de equilibrio son:

-ρü + ∇·σ + f = 0                    (momento lineal)
-ρJ_ρθ̈ + ∇·μ + E:σ + z = 0         (momento angular)

2. Principio de Equivalencia de Trabajo Virtual

La suposición clave es que desplazamientos virtuales y rotaciones virtuales pueden expresarse como series de coordenadas espaciales:

δu_j(x) = α_j^(0) + x_iα_ij^(1) + x_kx_iα_ikj^(2) + ...
δθ_j(x) = β_j^(0) + x_iβ_ij^(1) + ...

A través de la condición de equivalencia de trabajo virtual:

δW_int^mac = δW̄_int^mic = δW̄_ext^mic

3. Derivación del Tensor de Esfuerzo Macroscópico

Para modos de deformación de tracción y corte (términos α^(1)), se obtiene:

σ_ij^mac = ⟨σ_ij^mic⟩ = (1/V)∫_Ω σ_ij^mic dV        (forma de fuerzas internas)
σ_ij^mac = (1/V)[∫_Γ x_it_j^mic dΓ + ∫_Ω x_if_j^mic dV]  (forma de fuerzas externas)

4. Derivación del Tensor de Esfuerzo de Pareja Macroscópico

Para modos de deformación de flexión y torsión (términos β^(1)), se obtiene:

μ_ij^mac = ⟨μ_ij^mic + E_jklx_kσ_il^mic⟩ - E_jklx_k^macσ_il^mac

El hallazgo clave es que el tensor de esfuerzo de pareja contiene un término adicional que depende de la ubicación del punto de referencia macroscópico.

Tratamiento de Sistemas Discretos

1. Discretización de la Formulación de Fuerzas Externas

Discretización directa de expresiones continuas:

σ^mac = (1/V)∑_f x_f ⊗ F_f                    (tensor de esfuerzo)
μ^mac = σ^mac ⊗ x^mac : E + (1/V)∑_z x_z ⊗ Z_z  (tensor de esfuerzo de pareja)
a^mac = -(1/V)∑_q x_qQ_q                       (vector de flujo)

2. Formulación de Fuerzas Internas (Tipo Love-Weber)

A través del equilibrio de fuerzas y momentos, se obtiene la generalización de la famosa fórmula de Love-Weber:

σ^mac ≈ (1/V)∑_e A_el_ee_N^e ⊗ t^e
μ^mac ≈ σ^mac ⊗ x^mac : E + (1/V)∑_e A_el_ee_N^e ⊗ [m^e + E:(x_c^e ⊗ t^e)]
a^mac ≈ (1/V)∑_d S_dh_de_λ^d j_d

3. Término de Brecha Frontera-Radio

La formulación de fuerzas internas ignora la ubicación real de las fuerzas externas, introduciendo un término de error "brecha frontera-radio" que puede despreciarse para sistemas con gran cantidad de partículas.

Configuración Experimental

Estrategia de Verificación

El artículo emplea dos clases de métodos de verificación:

1. Verificación del Problema de Poisson

• Geometría: Dominio cuadrado de 1×1 m²
• Discretización: Aproximadamente 48,000 nodos, 72,000 elementos de conducto
• Casos de Prueba:
  • Prueba de parche de flujo constante: p = 2(2x-1) + 2(2y-1), a = (-4,-4)^T
  • Solución analítica no lineal: p = 16x₁x₂(1-x₁)(1-x₂)

2. Verificación de Problemas Mecánicos

• Geometría: Viga en voladizo, profundidad D=3m, longitud S=12m
• Discretización: 1,034,196 grados de libertad
• Parámetros de Material: E₀=40GPa, α=0.3, β=0~10⁵
• Modelos de Comparación:
  • Modelo discreto
  • Modelo de Cosserat homogeneizado
  • Modelo de Cosserat simplificado
  • Modelo de continuo de Cauchy

Indicadores de Evaluación

• Evaluación de Precisión: Desviación respecto a soluciones analíticas
• Convergencia: Resultados bajo diferentes tamaños de malla
• Consistencia: Comparación de formulaciones de fuerzas internas y externas
• Efectos de Frontera: Impacto del término de brecha frontera-radio

Resultados Experimentales

Resultados Principales

1. Verificación del Problema de Poisson

• Variante I (incluyendo término de brecha frontera-radio): Proporciona siempre soluciones exactas
• Variante II (ignorando término de brecha frontera-radio): Precisión disminuye con reducción del tamaño de malla
• Convergencia: Resultados similares entre variantes con malla de tamaño grande

2. Verificación de Problemas Mecánicos

Análisis Estacionario:

• Buen acuerdo entre modelo discreto y modelo homogeneizado
• Con aumento de rigidez de flexión β, esfuerzo disminuye mientras esfuerzo de pareja aumenta
• Cuando β<10³, efectos de Cosserat pueden despreciarse, puede usarse continuo de Cauchy

Análisis Transitorio:

• Respuesta de desplazamiento altamente consistente entre todos los modelos
• Evolución temporal de esfuerzo y esfuerzo de pareja conforme a lo esperado
• Verifica validez de fórmulas bajo condiciones transitorias

Hallazgos Clave

1. Efecto de Brecha Frontera-Radio

Experimentos confirman importancia del término de brecha frontera-radio:

• Para volúmenes de control de pequeño tamaño, este término no puede ignorarse
• En sistemas con gran cantidad de partículas, impacto de este término es menor
• Proporciona evaluación cuantitativa de errores

2. Influencia de la Ubicación del Punto de Referencia

Verifica predicciones teóricas del término de ubicación del punto de referencia en tensor de esfuerzo de pareja:

• Cuando punto macroscópico está en centroide del dominio, término adicional desaparece
• Otras ubicaciones producen contribuciones adicionales de esfuerzo de pareja

3. Influencia de Parámetros de Material

• Parámetro β afecta significativamente longitud característica de Cosserat ℓ
• Cuando β=0, ℓ≈4.6mm; cuando β=10⁵, ℓ≈0.7m
• Parámetros de Lamé λ,μ relativamente estables

Trabajos Relacionados

Fundamentos Teóricos Clásicos

• Fórmula de Love-Weber: Formulación clásica de tensor de esfuerzo para sistemas discretos
• Homogeneización por expansión asintótica: Teoría de RVE infinitesimal de Forest et al.
• Principio de trabajo virtual: Método de equivalencia energética de Chang y Kuhn

Innovaciones de Este Trabajo

1. RVE de tamaño finito: Considera efectos de tamaño en comparación con teoría asintótica
2. Tratamiento completo de esfuerzo de pareja: Incluye contribución de esfuerzo microscópico sobre distancia
3. Marco unificado: Maneja simultáneamente sistemas continuos y discretos
4. Dos formulaciones: Correspondencia completa entre fuerzas internas y externas

Conclusiones y Discusión

Conclusiones Principales

1. Contribución Teórica: Establece teoría de homogeneización rigurosa para RVE finito, generaliza fórmula clásica de Love-Weber
2. Innovación Metodológica: Propone método completo para manejar dependencia de punto de referencia en tensor de esfuerzo de pareja
3. Valor Práctico: Proporciona herramientas más precisas para modelado multiesacala

Limitaciones

1. Selección de Modos de Deformación Virtual: Ciertos equilibrios de trabajo virtual de orden superior se ignoran
2. Restricciones de Condiciones de Frontera: Requiere que todas las fronteras sean condiciones de frontera estática
3. Complejidad Computacional: Término de brecha frontera-radio aumenta carga computacional

Direcciones Futuras

1. Extensión No Lineal: Generalización a relaciones constitutivas no lineales
2. Problemas Acoplados: Extensión a problemas acoplados de mecánica-transporte
3. Teoría de Orden Superior: Consideración de modos de deformación virtual de orden superior

Evaluación Profunda

Ventajas

1. Rigor Teórico: Derivación matemática rigurosa basada en principio de trabajo virtual
2. Practicidad Fuerte: Directamente aplicable a modelos discretos y continuos existentes
3. Verificación Suficiente: Incluye soluciones analíticas, comparaciones numéricas y estudios paramétricos
4. Escritura Clara: Lógica clara, expresión matemática precisa

Insuficiencias

1. Rango de Aplicabilidad: Principalmente enfocado en problemas de elasticidad lineal, extensión no lineal requiere investigación adicional
2. Eficiencia Computacional: Costo computacional de fórmulas completas relativamente alto
3. Sensibilidad Paramétrica: Ciertos resultados muy sensibles al parámetro de material β

Impacto

1. Valor Académico: Proporciona base teórica importante para mecánica multiesacala
2. Aplicación de Ingeniería: Puede usarse directamente en modelado de hormigón, geotecnia y otros materiales heterogéneos
3. Implementación de Software: Código abierto disponible, facilita difusión y aplicación

Escenarios Aplicables

• Modelado multiesacala de materiales granulares
• Homogeneización de materiales compuestos heterogéneos
• Sistemas de materiales con rotación microscópica
• Problemas de ingeniería que requieren consideración de efectos de tamaño

Referencias Bibliográficas

El artículo cita 63 referencias relacionadas, incluyendo principalmente:

• Literatura clásica de teoría de continuo de Cosserat
• Método de elementos discretos y teoría de homogeneización
• Trabajos relacionados con modelado multiesacala
• Literatura de métodos numéricos utilizados en verificación

Evaluación General: Este es un artículo excelente con teoría rigurosa y alto valor práctico, proporcionando herramientas teóricas y métodos numéricos importantes para modelado multiesacala de materiales heterogéneos. El proceso de derivación es claro, la verificación es suficiente, y posee importante valor académico y aplicado para campos relacionados.