Tensori di stress macroscopico, stress di coppia e flusso derivati attraverso equivalenza energetica da volumi rappresentativi eterogenei finiti microscopici continui e discreti

Informazioni Fondamentali

• ID Articolo: 2412.12616
• Titolo: Macroscopic stress, couple stress and flux tensors derived through energetic equivalence from microscopic continuous and discrete heterogeneous finite representative volumes
• Autori: Jan Eliáš (Brno University of Technology), Gianluca Cusatis (Northwestern University)
• Classificazione: cs.CE (Computational Engineering), cond-mat.other
• Data di Presentazione: Dicembre 2024
• Link Articolo: https://arxiv.org/abs/2412.12616

Riassunto

Il presente articolo propone un metodo di derivazione rigoroso per valutare i tensori di stress macroscopico, i tensori di stress di coppia e i vettori di flusso in sistemi eterogenei continui e discreti con spostamenti e rotazioni indipendenti. A differenza dell'omogeneizzazione classica mediante espansione asintotica, il presente lavoro considera volumi rappresentativi (RVE) di dimensioni finite. Inizialmente vengono derivate le quantità macroscopiche per il continuo di Cosserat eterogeneo, successivamente le equazioni del continuo vengono discretizzate per fornire quantità macroscopiche in sistemi eterogenei discreti, e infine le espressioni per il sistema discreto vengono riderivate direttamente dalle proprietà discrete.

Contesto di Ricerca e Motivazione

1. Problemi da Risolvere

L'articolo affronta principalmente il problema di come estrarre accuratamente le proprietà meccaniche e di trasporto macroscopiche da sistemi di materiali eterogenei a scala microscopica, in particolare:

• Come derivare il tensore di stress macroscopico da volumi rappresentativi di dimensioni finite
• Come gestire il tensore di stress di coppia in sistemi con gradi di libertà rotazionali indipendenti
• Come stabilire relazioni di equivalenza energetica tra scale microscopiche e macroscopiche

2. Importanza del Problema

Molti materiali ingegneristici (come calcestruzzo, ceramiche, materiali granulari, ecc.) mostrano eterogeneità significativa a scala mesoscopica o microscopica. La modellazione multiscala accurata è essenziale per:

• Prevedere il comportamento meccanico complessivo dei materiali
• Sviluppare modelli numerici affidabili
• Stabilire relazioni costitutive fenomenologiche

3. Limitazioni dei Metodi Esistenti

• Omogeneizzazione mediante espansione asintotica: Assume che il volume rappresentativo tenda a zero, trascurando gli effetti di dimensione finita
• Formula classica di Love-Weber: Principalmente applicabile al tensore di stress, con trattamento incompleto del tensore di stress di coppia
• Gestione delle condizioni al contorno: I metodi esistenti presentano insufficienze nel trattamento degli effetti di bordo e nella scelta del punto di riferimento

4. Motivazione della Ricerca

L'articolo mira a stabilire un quadro teorico più completo e rigoroso, capace di:

• Gestire volumi rappresentativi di dimensioni finite
• Considerare simultaneamente formulazioni basate su forze interne ed esterne
• Chiarire l'influenza della posizione del punto di riferimento nel tensore di stress di coppia

Contributi Principali

1. Stabilimento di un quadro di derivazione rigoroso basato sull'equivalenza del lavoro virtuale, applicabile a sistemi eterogenei continui e discreti
2. Derivazione di espressioni per quantità macroscopiche considerando la dimensione finita dell'RVE, inclusi tensore di stress, tensore di stress di coppia e vettore di flusso
3. Proposizione di due formulazioni varianti: basate su forze interne e su forze esterne
4. Scoperta di termini aggiuntivi nel tensore di stress di coppia, che spiegano la sua dipendenza dal sistema di riferimento e dalla posizione del punto macroscopico
5. Fornitura di verifiche complete, incluso il confronto tra soluzioni analitiche in condizioni stazionarie e transitorie e modelli numerici

Dettagli Metodologici

Definizione del Compito

Il compito centrale dell'articolo è stabilire relazioni di equivalenza tra materiali eterogenei microscopici e continui omogenei macroscopici, includendo specificamente:

• Input: Campo di stress microscopico σ^mic, campo di stress di coppia μ^mic, campo di flusso a^mic
• Output: Tensore di stress macroscopico σ^mac, tensore di stress di coppia μ^mac, vettore di flusso a^mac
• Vincoli: Soddisfare il principio del lavoro virtuale e le equazioni di equilibrio

Quadro Teorico

1. Fondamenti del Continuo di Cosserat

Le variabili principali del continuo di Cosserat includono:

• Vettore di spostamento u e vettore di rotazione indipendente θ
• Tensore di deformazione γ = ε - E·θ^m e tensore di curvatura κ = ∇⊗θ
• Tensore di stress σ (non simmetrico) e tensore di stress di coppia μ

Le equazioni di equilibrio sono:

-ρü + ∇·σ + f = 0                    (quantità di moto lineare)
-ρJ_ρθ̈ + ∇·μ + E:σ + z = 0         (quantità di moto angolare)

2. Principio del Lavoro Virtuale Equivalente

L'ipotesi chiave è che gli spostamenti virtuali e le rotazioni virtuali possono essere espressi come serie di coordinate spaziali:

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

Attraverso la condizione di equivalenza del lavoro virtuale:

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

3. Derivazione del Tensore di Stress Macroscopico

Per i modi di deformazione di trazione e taglio (termini α^(1)), si ottiene:

σ_ij^mac = ⟨σ_ij^mic⟩ = (1/V)∫_Ω σ_ij^mic dV        (formulazione forze interne)
σ_ij^mac = (1/V)[∫_Γ x_it_j^mic dΓ + ∫_Ω x_if_j^mic dV]  (formulazione forze esterne)

4. Derivazione del Tensore di Stress di Coppia Macroscopico

Per i modi di deformazione di flessione e torsione (termini β^(1)), si ottiene:

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

La scoperta chiave è che il tensore di stress di coppia contiene un termine aggiuntivo che dipende dalla posizione del punto di riferimento macroscopico.

Trattamento dei Sistemi Discreti

1. Discretizzazione della Formulazione Basata su Forze Esterne

Discretizzazione diretta dell'espressione del continuo:

σ^mac = (1/V)∑_f x_f ⊗ F_f                    (tensore di stress)
μ^mac = σ^mac ⊗ x^mac : E + (1/V)∑_z x_z ⊗ Z_z  (tensore di stress di coppia)
a^mac = -(1/V)∑_q x_qQ_q                       (vettore di flusso)

2. Formulazione Basata su Forze Interne (Tipo Love-Weber)

Attraverso l'equilibrio di forze e momenti, si ottiene l'estensione della celebre formula di 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. Termine di Lacuna Bordo-Raggio

La formulazione basata su forze interne trascura la posizione reale delle forze esterne, introducendo un termine di errore "lacuna bordo-raggio", che può essere trascurato per sistemi con un gran numero di particelle.

Configurazione Sperimentale

Strategia di Verifica

L'articolo impiega due classi di metodi di verifica:

1. Verifica del Problema di Poisson

• Geometria: Dominio quadrato 1×1 m²
• Discretizzazione: Circa 48.000 nodi, 72.000 elementi condotto
• Casi di prova:
  • Prova di patch a flusso costante: p = 2(2x-1) + 2(2y-1), a = (-4,-4)^T
  • Soluzione analitica non lineare: p = 16x₁x₂(1-x₁)(1-x₂)

2. Verifica del Problema Meccanico

• Geometria: Trave a sbalzo, profondità D=3m, lunghezza S=12m
• Discretizzazione: 1.034.196 gradi di libertà
• Parametri materiali: E₀=40GPa, α=0,3, β=0~10⁵
• Modelli di confronto:
  • Modello discreto
  • Modello Cosserat omogeneizzato
  • Modello Cosserat semplificato
  • Modello continuo di Cauchy

Indicatori di Valutazione

• Valutazione della precisione: Deviazione dalle soluzioni analitiche
• Convergenza: Risultati con diverse dimensioni di mesh
• Coerenza: Confronto tra formulazioni basate su forze interne ed esterne
• Effetti di bordo: Influenza del termine di lacuna bordo-raggio

Risultati Sperimentali

Risultati Principali

1. Verifica del Problema di Poisson

• Variante I (con termine di lacuna bordo-raggio): Fornisce sempre soluzioni esatte
• Variante II (omettendo il termine di lacuna bordo-raggio): La precisione diminuisce con la riduzione della dimensione della mesh
• Convergenza: I risultati delle due varianti sono simili con mesh di grandi dimensioni

2. Verifica del Problema Meccanico

Analisi stazionaria:

• Il modello discreto concorda bene con il modello omogeneizzato
• Con l'aumento della rigidità di flessione β, lo stress diminuisce mentre lo stress di coppia aumenta
• Quando β<10³, gli effetti di Cosserat sono trascurabili e può essere utilizzato il continuo di Cauchy

Analisi transitoria:

• La risposta di spostamento di tutti i modelli è altamente coerente
• L'evoluzione temporale dello stress e dello stress di coppia è conforme alle aspettative
• Verifica l'efficacia delle formule in condizioni transitorie

Scoperte Chiave

1. Effetto della Lacuna Bordo-Raggio

L'esperimento conferma l'importanza del termine di lacuna bordo-raggio:

• Per volumi di controllo di piccole dimensioni, il termine non può essere trascurato
• In sistemi con un gran numero di particelle, l'influenza del termine è minore
• Fornisce una valutazione quantitativa dell'errore

2. Influenza della Posizione del Punto di Riferimento

Verifica le previsioni teoriche del termine di posizione del punto di riferimento nel tensore di stress di coppia:

• Quando il punto macroscopico si trova al baricentro del dominio, il termine aggiuntivo scompare
• Altre posizioni producono contributi di stress di coppia aggiuntivi

3. Influenza dei Parametri Materiali

• Il parametro β influenza significativamente la lunghezza caratteristica di Cosserat ℓ
• Quando β=0, ℓ≈4,6mm; quando β=10⁵, ℓ≈0,7m
• I parametri di Lamé λ,μ rimangono relativamente stabili

Lavori Correlati

Fondamenti Teorici Classici

• Formula di Love-Weber: Formulazione classica del tensore di stress per sistemi discreti
• Omogeneizzazione mediante espansione asintotica: Teoria dell'RVE infinitesimale di Forest e altri
• Principio del lavoro virtuale: Metodo di equivalenza energetica di Chang e Kuhn

Innovazioni del Presente Articolo

1. RVE di dimensioni finite: Considera gli effetti di dimensione rispetto alla teoria asintotica
2. Trattamento completo dello stress di coppia: Include il contributo dello stress microscopico sulla distanza
3. Quadro unificato: Gestisce simultaneamente sistemi continui e discreti
4. Due formulazioni: Corrispondenza completa tra forze interne ed esterne

Conclusioni e Discussione

Conclusioni Principali

1. Contributo teorico: Stabilimento di una teoria di omogeneizzazione rigorosa per RVE di dimensioni finite, estensione della classica formula di Love-Weber
2. Innovazione metodologica: Proposizione di un metodo completo per gestire la dipendenza dal punto di riferimento del tensore di stress di coppia
3. Valore pratico: Fornitura di strumenti più precisi per la modellazione multiscala

Limitazioni

1. Scelta dei modi di deformazione virtuali: Alcuni equilibri di lavoro virtuale di ordine superiore vengono omessi
2. Restrizioni sulle condizioni al contorno: Richiede che tutti i bordi siano condizioni al contorno statiche
3. Complessità computazionale: Il termine di lacuna bordo-raggio aumenta l'onere computazionale

Direzioni Future

1. Estensione non lineare: Generalizzazione a relazioni costitutive non lineari
2. Problemi accoppiati: Estensione a problemi accoppiati meccanica-trasporto
3. Teoria di ordine superiore: Considerazione di modi di deformazione virtuale di ordine superiore

Valutazione Approfondita

Vantaggi

1. Rigore teorico: Derivazione matematica rigorosa basata sul principio del lavoro virtuale
2. Forte praticità: Direttamente applicabile a modelli discreti e continui esistenti
3. Verifica sufficiente: Incluso il confronto tra soluzioni analitiche, modelli numerici e studi parametrici
4. Scrittura chiara: Logica trasparente, espressione matematica accurata

Insufficienze

1. Ambito di applicabilità: Principalmente focalizzato su problemi di elasticità lineare, l'estensione non lineare richiede ulteriori ricerche
2. Efficienza computazionale: Il costo computazionale della formula completa è relativamente elevato
3. Sensibilità ai parametri: Alcuni risultati sono molto sensibili al parametro materiale β

Impatto

1. Valore accademico: Fornisce fondamenti teorici importanti per la meccanica multiscala
2. Applicazione ingegneristica: Direttamente utilizzabile per la modellazione di materiali eterogenei come calcestruzzo e geotecnici
3. Implementazione software: Codice open source fornito, facilitando la diffusione e l'applicazione

Scenari Applicabili

• Modellazione multiscala di materiali granulari
• Omogeneizzazione di materiali compositi eterogenei
• Sistemi materiali con rotazioni microscopiche
• Problemi ingegneristici che richiedono considerazione degli effetti di dimensione

Bibliografia

L'articolo cita 63 riferimenti correlati, principalmente includenti:

• Letteratura classica sulla teoria del continuo di Cosserat
• Metodo degli elementi discreti e teoria dell'omogeneizzazione
• Lavori correlati sulla modellazione multiscala
• Letteratura sui metodi numerici utilizzati per la verifica

Valutazione Complessiva: Questo è un articolo eccellente con teoria rigorosa e alto valore pratico, fornendo strumenti teorici e metodi numerici importanti per la modellazione multiscala di materiali eterogenei. Il processo di derivazione dell'articolo è chiaro, la verifica è sufficiente, e possiede importante valore accademico e applicativo per i campi correlati.