Tenseurs de contrainte macroscopique, de contrainte de couple et de flux dérivés par équivalence énergétique à partir de volumes représentatifs finis hétérogènes microscopiques continus et discrets

Informations fondamentales

• ID de l'article: 2412.12616
• Titre: Macroscopic stress, couple stress and flux tensors derived through energetic equivalence from microscopic continuous and discrete heterogeneous finite representative volumes
• Auteurs: Jan Eliáš (Université de technologie de Brno), Gianluca Cusatis (Université Northwestern)
• Classification: cs.CE (Ingénierie computationnelle), cond-mat.other
• Date de soumission: Décembre 2024
• Lien de l'article: https://arxiv.org/abs/2412.12616

Résumé

Cet article propose une méthode de dérivation rigoureuse pour évaluer les tenseurs de contrainte macroscopique, les tenseurs de contrainte de couple et les vecteurs de flux dans les systèmes hétérogènes continus et discrets possédant des déplacements et des rotations indépendants. Contrairement à l'homogénéisation classique par développement asymptotique, cet article considère les volumes représentatifs élémentaires (VRE) de taille finie. Les quantités macroscopiques sont d'abord dérivées pour le continuum de Cosserat hétérogène, puis les équations continues sont discrétisées pour fournir les quantités macroscopiques dans les systèmes hétérogènes discrets, et enfin les expressions pour les systèmes discrets sont redérivées directement à partir des propriétés discrètes.

Contexte et motivation de la recherche

1. Problèmes à résoudre

Cet article aborde principalement la question de l'extraction précise des propriétés mécaniques et de transport macroscopiques à partir de systèmes de matériaux hétérogènes à l'échelle microscopique, en particulier:

• Comment dériver le tenseur de contrainte macroscopique à partir d'un volume représentatif élémentaire de taille finie
• Comment traiter le tenseur de contrainte de couple dans les systèmes possédant des degrés de liberté de rotation indépendants
• Comment établir une relation d'équivalence énergétique entre les échelles microscopique et macroscopique

2. Importance du problème

De nombreux matériaux d'ingénierie (tels que le béton, la céramique, les matériaux granulaires, etc.) présentent une hétérogénéité significative aux échelles mésoscopique ou microscopique. La modélisation multi-échelle précise est importante pour:

• Prédire le comportement mécanique global des matériaux
• Développer des modèles numériques fiables
• Établir des relations constitutives phénoménologiques

3. Limitations des méthodes existantes

• Homogénéisation par développement asymptotique: Suppose que le volume représentatif élémentaire tend vers zéro, ignorant les effets de taille finie
• Formule classique de Love-Weber: Principalement applicable au tenseur de contrainte, avec un traitement insuffisant du tenseur de contrainte de couple
• Traitement des conditions aux limites: Les méthodes existantes présentent des insuffisances dans le traitement des effets de bord et du choix du point de référence

4. Motivation de la recherche

Cet article vise à établir un cadre théorique plus complet et rigoureux capable de:

• Traiter les volumes représentatifs élémentaires de taille finie
• Considérer simultanément les formulations en forces internes et externes
• Clarifier l'influence de la position du point de référence dans le tenseur de contrainte de couple

Contributions fondamentales

1. Établissement d'un cadre de dérivation rigoureux basé sur l'équivalence du travail virtuel, applicable aux systèmes hétérogènes continus et discrets
2. Dérivation des expressions de quantités macroscopiques tenant compte de la taille finie du VRE, incluant le tenseur de contrainte, le tenseur de contrainte de couple et le vecteur de flux
3. Proposition de deux variantes formelles: formulations basées sur les forces internes et externes
4. Découverte de termes supplémentaires dans le tenseur de contrainte de couple, expliquant sa dépendance au système de référence et à la position du point macroscopique
5. Fourniture de vérifications complètes, incluant des comparaisons entre solutions analytiques et modèles numériques en conditions stationnaires et transitoires

Détails de la méthode

Définition de la tâche

La tâche fondamentale de cet article est d'établir la relation d'équivalence entre les matériaux hétérogènes microscopiques et le continuum homogène macroscopique, incluant spécifiquement:

• Entrées: Champ de contrainte microscopique σ^mic, champ de contrainte de couple μ^mic, champ de flux a^mic
• Sorties: Tenseur de contrainte macroscopique σ^mac, tenseur de contrainte de couple μ^mac, vecteur de flux macroscopique a^mac
• Contraintes: Satisfaction du principe du travail virtuel et des équations d'équilibre

Cadre théorique

1. Fondements du continuum de Cosserat

Les variables principales du continuum de Cosserat incluent:

• Vecteur de déplacement u et vecteur de rotation indépendant θ
• Tenseur de déformation γ = ε - E·θ^m et tenseur de courbure κ = ∇⊗θ
• Tenseur de contrainte σ (non-symétrique) et tenseur de contrainte de couple μ

Les équations d'équilibre sont:

-ρü + ∇·σ + f = 0                    (moment linéaire)
-ρJ_ρθ̈ + ∇·μ + E:σ + z = 0         (moment angulaire)

2. Principe du travail virtuel équivalent

L'hypothèse clé est que les déplacements virtuels et les rotations virtuelles peuvent être exprimés comme des séries de coordonnées spatiales:

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

Par la condition d'équivalence du travail virtuel:

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

3. Dérivation du tenseur de contrainte macroscopique

Pour les modes de déformation en traction et cisaillement (termes α^(1)), on obtient:

σ_ij^mac = ⟨σ_ij^mic⟩ = (1/V)∫_Ω σ_ij^mic dV        (formulation en forces internes)
σ_ij^mac = (1/V)[∫_Γ x_it_j^mic dΓ + ∫_Ω x_if_j^mic dV]  (formulation en forces externes)

4. Dérivation du tenseur de contrainte de couple macroscopique

Pour les modes de déformation en flexion et torsion (termes β^(1)), on obtient:

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

La découverte clé est que le tenseur de contrainte de couple contient un terme supplémentaire dépendant de la position du point de référence macroscopique.

Traitement des systèmes discrets

1. Discrétisation de la formulation en forces externes

Discrétisation directe de l'expression continue:

σ^mac = (1/V)∑_f x_f ⊗ F_f                    (tenseur de contrainte)
μ^mac = σ^mac ⊗ x^mac : E + (1/V)∑_z x_z ⊗ Z_z  (tenseur de contrainte de couple)
a^mac = -(1/V)∑_q x_qQ_q                       (vecteur de flux)

2. Formulation en forces internes (type Love-Weber)

Par l'équilibre des forces et des moments, on obtient une généralisation de la célèbre formule 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. Terme d'écart limite-rayon

La formulation en forces internes ignore la position réelle des forces externes, introduisant un terme d'erreur "d'écart limite-rayon", qui peut être négligé pour les systèmes contenant un grand nombre de particules.

Configuration expérimentale

Stratégie de vérification

L'article emploie deux classes de méthodes de vérification:

1. Vérification du problème de Poisson

• Géométrie: Domaine carré 1×1 m²
• Discrétisation: Environ 48 000 nœuds, 72 000 éléments de conduit
• Cas de test:
  • Test de patch à flux constant: p = 2(2x-1) + 2(2y-1), a = (-4,-4)^T
  • Solution analytique non-linéaire: p = 16x₁x₂(1-x₁)(1-x₂)

2. Vérification du problème mécanique

• Géométrie: Poutre en porte-à-faux, profondeur D=3m, longueur S=12m
• Discrétisation: 1 034 196 degrés de liberté
• Paramètres matériaux: E₀=40GPa, α=0,3, β=0~10⁵
• Modèles de comparaison:
  • Modèle discret
  • Modèle de Cosserat homogénéisé
  • Modèle de Cosserat simplifié
  • Modèle de continuum de Cauchy

Indicateurs d'évaluation

• Évaluation de la précision: Écart par rapport à la solution analytique
• Convergence: Résultats pour différentes tailles de maille
• Cohérence: Comparaison des formulations en forces internes et externes
• Effets de bord: Impact du terme d'écart limite-rayon

Résultats expérimentaux

Résultats principaux

1. Vérification du problème de Poisson

• Variante I (incluant le terme d'écart limite-rayon): Fournit toujours la solution exacte
• Variante II (négligeant le terme d'écart limite-rayon): La précision diminue avec la réduction de la taille de maille
• Convergence: Les deux variantes donnent des résultats similaires pour les mailles de grande taille

2. Vérification du problème mécanique

Analyse stationnaire:

• Le modèle discret s'accorde bien avec le modèle homogénéisé
• Avec l'augmentation de la rigidité en flexion β, la contrainte diminue tandis que la contrainte de couple augmente
• Pour β<10³, les effets de Cosserat peuvent être négligés, permettant l'utilisation du continuum de Cauchy

Analyse transitoire:

• La réponse en déplacement de tous les modèles est hautement cohérente
• L'évolution temporelle de la contrainte et de la contrainte de couple est conforme aux attentes
• Valide l'efficacité des formules en conditions transitoires

Découvertes clés

1. Effet d'écart limite-rayon

L'expérience confirme l'importance du terme d'écart limite-rayon:

• Pour les petits volumes de contrôle, ce terme ne peut pas être négligé
• Dans les systèmes contenant un grand nombre de particules, l'impact de ce terme est plus faible
• Fournit une évaluation d'erreur quantitative

2. Influence de la position du point de référence

Valide les prédictions théoriques du terme de position du point de référence dans le tenseur de contrainte de couple:

• Lorsque le point macroscopique est situé au centre de masse du domaine, le terme supplémentaire disparaît
• D'autres positions produisent une contribution supplémentaire de contrainte de couple

3. Influence des paramètres matériaux

• Le paramètre β affecte significativement la longueur caractéristique de Cosserat ℓ
• Pour β=0, ℓ≈4,6mm; pour β=10⁵, ℓ≈0,7m
• Les paramètres de Lamé λ,μ restent relativement stables

Travaux connexes

Fondements théoriques classiques

• Formule de Love-Weber: Formulation classique du tenseur de contrainte pour les systèmes discrets
• Homogénéisation par développement asymptotique: Théorie du VRE infinitésimal de Forest et al.
• Principe du travail virtuel: Méthode d'équivalence énergétique de Chang et Kuhn

Points novateurs de cet article

1. VRE de taille finie: Considère les effets de taille par rapport à la théorie asymptotique
2. Traitement complet de la contrainte de couple: Inclut la contribution de la contrainte microscopique sur la distance
3. Cadre unifié: Traite simultanément les systèmes continus et discrets
4. Deux formulations: Correspondance complète entre les forces internes et externes

Conclusions et discussion

Conclusions principales

1. Contribution théorique: Établissement d'une théorie d'homogénéisation rigoureuse pour les VRE de taille finie, généralisant la formule classique de Love-Weber
2. Innovation méthodologique: Proposition d'une méthode complète pour traiter la dépendance du point de référence du tenseur de contrainte de couple
3. Valeur pratique: Fournit des outils plus précis pour la modélisation multi-échelle

Limitations

1. Choix des modes de déformation virtuelle: Certains équilibres de travail virtuel d'ordre supérieur sont négligés
2. Restrictions sur les conditions aux limites: Exige que toutes les limites soient des conditions aux limites statiques
3. Complexité computationnelle: Le terme d'écart limite-rayon augmente la charge de calcul

Directions futures

1. Extension non-linéaire: Généralisation aux relations constitutives non-linéaires
2. Problèmes couplés: Extension aux problèmes couplés mécanique-transport
3. Théories d'ordre supérieur: Considération de modes de déformation virtuelle d'ordre plus élevé

Évaluation approfondie

Avantages

1. Rigueur théorique: Dérivation mathématique rigoureuse basée sur le principe du travail virtuel
2. Forte praticabilité: Directement applicable aux modèles discrets et continus existants
3. Vérification suffisante: Inclut des solutions analytiques, des comparaisons numériques et des études paramétriques
4. Rédaction claire: Logique claire et expression mathématique précise

Insuffisances

1. Portée d'application: Principalement orientée vers les problèmes d'élasticité linéaire; l'extension non-linéaire nécessite des recherches supplémentaires
2. Efficacité computationnelle: Le coût de calcul des formules complètes est relativement élevé
3. Sensibilité aux paramètres: Certains résultats sont très sensibles au paramètre matériau β

Impact

1. Valeur académique: Fournit une base théorique importante pour la mécanique multi-échelle
2. Application en ingénierie: Peut être directement utilisée pour la modélisation de matériaux hétérogènes tels que le béton et les géomatériaux
3. Implémentation logicielle: Code open-source fourni, facilitant la promotion et l'application

Scénarios applicables

• Modélisation multi-échelle de matériaux granulaires
• Homogénéisation de matériaux composites hétérogènes
• Systèmes de matériaux possédant des rotations microscopiques
• Problèmes d'ingénierie nécessitant de considérer les effets de taille

Références bibliographiques

L'article cite 63 références pertinentes, incluant principalement:

• Littérature classique sur la théorie du continuum de Cosserat
• Méthode des éléments discrets et théorie de l'homogénéisation
• Travaux connexes en modélisation multi-échelle
• Littérature sur les méthodes numériques utilisées pour la vérification

Évaluation globale: Cet article est une excellente contribution théoriquement rigoureuse et pratiquement utile, fournissant des outils théoriques et des méthodes numériques importants pour la modélisation multi-échelle de matériaux hétérogènes. Le processus de dérivation est clair, la vérification est suffisante, et l'article possède une valeur académique et une utilité pratique significatives pour le domaine connexe.