2025-11-14T17:43:11.826412

Macroscopic stress, couple stress and flux tensors derived through energetic equivalence from microscopic continuous and discrete heterogeneous finite representative volumes

EliÃ¡Å¡, Cusatis

This paper presents a rigorous derivation of equations to evaluate the macroscopic stress tensor, the couple stress tensor, and the flux vector equivalent to underlying microscopic fields in continuous and discrete heterogeneous systems with independent displacements and rotations. Contrary to the classical asymptotic expansion homogenization, finite size representative volume is considered. First, the macroscopic quantities are derived for a heterogeneous Cosserat continuum. The resulting continuum equations are discretized to provide macroscopic quantities in discrete heterogeneous systems. Finally, the expressions for discrete system are derived once again, this time considering the discrete nature directly. The formulations are presented in two variants, considering either internal or external forces, couples, and fluxes. The derivation is based on the virtual work equivalence and elucidates the fundamental significance of the couple stress tensor in the context of balance equations and admissible virtual deformation modes. Notably, an additional term in the couple stress tensor formula emerges, explaining its dependence on the reference system and position of the macroscopic point. The resulting equations are verified by comparing their predictions with known analytical solutions and results of other numerical models under both steady state and transient conditions.

academic

微視的連続および離散不均質有限代表体積要素からのエネルギー等価性を通じて導出された巨視的応力、偶応力およびフラックステンソル

基本情報

論文ID: 2412.12616
タイトル: Macroscopic stress, couple stress and flux tensors derived through energetic equivalence from microscopic continuous and discrete heterogeneous finite representative volumes
著者: Jan Eliáš (ブルノ工科大学)、Gianluca Cusatis (ノースウェスタン大学)
分類: cs.CE (計算工学)、cond-mat.other
提出時期: 2024年12月
論文リンク: https://arxiv.org/abs/2412.12616

概要

本論文は、独立した変位と回転を有する連続および離散不均質系における巨視的応力テンソル、偶応力テンソルおよびフラックスベクトルを評価するための厳密な導出方法を提案する。古典的な漸近展開均質化とは異なり、本論文は有限サイズの代表体積要素(RVE)を考慮する。まず不均質Cosserat連続体に対して巨視的量を導出し、次に連続体方程式を離散化して離散不均質系における巨視的量を提供し、最後に離散特性から直接離散系の表現式を再導出する。

研究背景と動機

1. 解決すべき問題

本論文は、微視的スケールの不均質材料系から巨視的力学および輸送特性を正確に抽出する方法に関する以下の問題に対処する：

有限サイズの代表体積要素から巨視的応力テンソルを導出する方法
独立した回転自由度を有する系における偶応力テンソルの取り扱い
微視的および巨視的スケール間のエネルギー等価関係の確立

2. 問題の重要性

コンクリート、セラミックス、粒状材料などの多くの工学材料は、細観的または微視的スケールで顕著な不均質性を示す。正確な多スケール建模は以下の点で重要である：

材料の全体的な力学挙動の予測
信頼性の高い数値モデルの開発
現象論的本構関係の確立

3. 既存方法の限界

漸近展開均質化: 代表体積要素がゼロに近づくと仮定し、有限サイズ効果を無視する
古典的Love-Weber公式: 主に応力テンソルに適用可能であり、偶応力テンソルの取り扱いが不十分である
境界条件の処理: 既存方法は境界効果と参照点選択の処理に不足がある

4. 研究の動機

本論文は、以下を可能にする、より完全で厳密な理論的枠組みの確立を目指す：

有限サイズの代表体積要素の処理
内力および外力の両方の表現方式の同時考慮
偶応力テンソルにおける参照点位置の影響の明確化

核心的貢献

仮想仕事等価性に基づく厳密な導出枠組みの確立：連続および離散不均質系に適用可能
有限RVEサイズを考慮した巨視的量の表現式の導出：応力テンソル、偶応力テンソルおよびフラックスベクトルを含む
2つの変体形式の提案：内力ベースおよび外力ベースの表現
偶応力テンソルにおける追加項の発見：参照系および巨視的点位置への依存性を説明
完全な検証の提供：定常および非定常条件下の解析解と数値モデルの比較

方法の詳細

タスク定義

本論文の核心的タスクは、微視的不均質材料と巨視的均質連続体間の等価関係を確立することであり、具体的には以下を含む：

入力: 微視的スケールの応力場σ^mic、偶応力場μ^mic、フラックス場a^mic
出力: 巨視的応力テンソルσ^mac、偶応力テンソルμ^mac、フラックスベクトルa^mac
制約: 仮想仕事等価原理および平衡方程式を満たす

理論的枠組み

1. Cosserat連続体の基礎

Cosserat連続体の主要変数は以下を含む：

変位ベクトル u および 独立回転ベクトル θ
ひずみテンソル γ = ε - E·θ^m および 曲率テンソル κ = ∇⊗θ
応力テンソル σ（非対称）および 偶応力テンソル μ

平衡方程式は以下の通りである：

-ρü + ∇·σ + f = 0                    (線運動量)
-ρJ_ρθ̈ + ∇·μ + E:σ + z = 0         (角運動量)

2. 仮想仕事等価原理

重要な仮定は、仮想変位および仮想回転が空間座標の級数として表現可能であることである：

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

仮想仕事等価条件を通じて：

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

3. 巨視的応力テンソルの導出

引張およびせん断変形モード(α^(1)項)に対して、以下を得る：

σ_ij^mac = ⟨σ_ij^mic⟩ = (1/V)∫_Ω σ_ij^mic dV        (内力形式)
σ_ij^mac = (1/V)[∫_Γ x_it_j^mic dΓ + ∫_Ω x_if_j^mic dV]  (外力形式)

4. 巨視的偶応力テンソルの導出

曲げおよび捩れ変形モード(β^(1)項)に対して、以下を得る：

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

重要な発見は、偶応力テンソルが巨視的参照点位置に依存する追加項を含むことである。

離散系の処理

1. 外力表現の離散化

連続体表現式を直接離散化する：

σ^mac = (1/V)∑_f x_f ⊗ F_f                    (応力テンソル)
μ^mac = σ^mac ⊗ x^mac : E + (1/V)∑_z x_z ⊗ Z_z  (偶応力テンソル)
a^mac = -(1/V)∑_q x_qQ_q                       (フラックスベクトル)

2. 内力表現(Love-Weber型)

力および力矩平衡を通じて、有名な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. 境界-半径間隙項

内力表現は外力の実際の位置を無視し、「境界-半径間隙」誤差項を導入する。大量の粒子系では無視可能である。

実験設定

検証戦略

論文は2つのクラスの検証方法を採用する：

1. Poisson問題の検証

幾何学: 1×1 m²正方形領域
離散化: 約48,000個のノード、72,000個の導管要素
テストケース:
- 定常フラックスパッチテスト：p = 2(2x-1) + 2(2y-1), a = (-4,-4)^T
- 非線形解析解：p = 16x₁x₂(1-x₁)(1-x₂)