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

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

基本信息

  • 论文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áš (Brno University of Technology), Gianluca Cusatis (Northwestern University)
  • 分类: cs.CE (Computational Engineering), cond-mat.other
  • 提交时间: 2024年12月
  • 论文链接: https://arxiv.org/abs/2412.12616

摘要

本文提出了一种严格的推导方法，用于评估连续和离散异质系统中宏观应力张量、偶应力张量和通量矢量，这些系统具有独立的位移和旋转。与经典的渐近展开均匀化不同，本文考虑了有限尺寸的代表性体积单元(RVE)。首先为异质Cosserat连续体推导宏观量，然后将连续体方程离散化以提供离散异质系统中的宏观量，最后直接从离散性质出发重新推导离散系统的表达式。

研究背景与动机

1. 要解决的问题

本文主要解决如何从微观尺度的异质材料系统中准确提取宏观力学和传输特性的问题，特别是：

  • 如何从有限尺寸的代表性体积单元中推导宏观应力张量
  • 如何处理具有独立旋转自由度的系统中的偶应力张量
  • 如何建立微观和宏观尺度之间的能量等价关系

2. 问题的重要性

许多工程材料（如混凝土、陶瓷、颗粒材料等）在细观或微观尺度上表现出显著的异质性。准确的多尺度建模对于：

  • 预测材料的整体力学行为
  • 开发可靠的数值模型
  • 建立现象学本构关系 具有重要意义。

3. 现有方法的局限性

  • 渐近展开均匀化：假设代表性体积单元趋于零，忽略了有限尺寸效应
  • 经典Love-Weber公式：主要适用于应力张量，对偶应力张量的处理不够完善
  • 边界条件处理：现有方法在处理边界效应和参考点选择方面存在不足

4. 研究动机

本文旨在建立一个更加完整和严格的理论框架，能够：

  • 处理有限尺寸的代表性体积单元
  • 同时考虑内力和外力两种表述方式
  • 明确偶应力张量中参考点位置的影响

核心贡献

  1. 建立了基于虚功等价的严格推导框架，适用于连续和离散异质系统
  2. 推导了考虑有限RVE尺寸的宏观量表达式，包括应力张量、偶应力张量和通量矢量
  3. 提出了两种变体形式：基于内力和基于外力的表述
  4. 发现了偶应力张量中的附加项，解释了其对参考系和宏观点位置的依赖性
  5. 提供了完整的验证，包括稳态和瞬态条件下的解析解和数值模型对比

方法详解

任务定义

本文的核心任务是建立微观异质材料与宏观均匀连续体之间的等价关系，具体包括：

  • 输入：微观尺度的应力场σ^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. 边界-半径间隙项

内力表述忽略了外力的真实位置，引入"边界-半径间隙"误差项，对于大量粒子系统可以忽略。

实验设置

验证策略

论文采用两类验证方法：

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₂)

2. 力学问题验证

  • 几何: 悬臂梁，深度D=3m，长度S=12m
  • 离散化: 1,034,196个自由度
  • 材料参数: E₀=40GPa, α=0.3, β=0~10⁵
  • 对比模型:
    • 离散模型
    • 均匀化Cosserat模型
    • 简化Cosserat模型
    • Cauchy连续体模型

评价指标

  • 精度评估: 与解析解的偏差
  • 收敛性: 不同网格尺寸下的结果
  • 一致性: 内力和外力表述的对比
  • 边界效应: 边界-半径间隙项的影响

实验结果

主要结果

1. Poisson问题验证

  • 变体I（含边界-半径间隙项）: 始终提供精确解
  • 变体II（忽略边界-半径间隙项）: 随着网格尺寸减小精度下降
  • 收敛性: 大尺寸网格时两种变体结果相近

2. 力学问题验证

稳态分析:

  • 离散模型与均匀化模型吻合良好
  • 随着弯曲刚度β增加，应力减小而偶应力增大
  • β<10³时Cosserat效应可忽略，可用Cauchy连续体

瞬态分析:

  • 所有模型的位移响应高度一致
  • 应力和偶应力的时间演化符合预期
  • 验证了公式在瞬态条件下的有效性

关键发现

1. 边界-半径间隙效应

实验证实了边界-半径间隙项的重要性：

  • 对于小尺寸控制体积，该项不可忽略
  • 大量粒子系统中该项影响较小
  • 提供了定量的误差评估

2. 参考点位置的影响

验证了偶应力张量中参考点位置项的理论预测：

  • 当宏观点位于域质心时，附加项消失
  • 其他位置会产生额外的偶应力贡献

3. 材料参数的影响

  • β参数显著影响Cosserat特征长度ℓ
  • β=0时ℓ≈4.6mm，β=10⁵时ℓ≈0.7m
  • Lamé参数λ,μ相对稳定

相关工作

经典理论基础

  • Love-Weber公式: 离散系统应力张量的经典表述
  • 渐近展开均匀化: Forest等人的无限小RVE理论
  • 虚功原理: Chang和Kuhn的能量等价方法

本文的创新点

  1. 有限尺寸RVE: 相比渐近理论考虑了尺寸效应
  2. 完整的偶应力处理: 包含了微观应力在距离上的贡献
  3. 统一框架: 同时处理连续和离散系统
  4. 两种表述: 内力和外力的完整对应关系

结论与讨论

主要结论

  1. 理论贡献: 建立了有限RVE的严格均匀化理论，推广了经典Love-Weber公式
  2. 方法创新: 提出了处理偶应力张量参考点依赖性的完整方法
  3. 实用价值: 为多尺度建模提供了更精确的工具

局限性

  1. 虚变形模式的选择: 某些高阶虚功平衡被忽略
  2. 边界条件限制: 要求所有边界为静力边界条件
  3. 计算复杂性: 边界-半径间隙项增加了计算负担

未来方向

  1. 非线性扩展: 推广到非线性本构关系
  2. 耦合问题: 扩展到力学-传输耦合问题
  3. 高阶理论: 考虑更高阶的虚变形模式

深度评价

优点

  1. 理论严谨性: 基于虚功原理的严格数学推导
  2. 实用性强: 直接适用于现有的离散和连续模型
  3. 验证充分: 包含解析解、数值对比和参数研究
  4. 写作清晰: 逻辑清楚，数学表述准确

不足

  1. 适用范围: 主要针对线弹性问题，非线性扩展需要进一步研究
  2. 计算效率: 完整公式的计算成本较高
  3. 参数敏感性: 某些结果对材料参数β非常敏感

影响力

  1. 学术价值: 为多尺度力学提供了重要的理论基础
  2. 工程应用: 可直接用于混凝土、岩土等异质材料建模
  3. 软件实现: 已提供开源代码，便于推广应用

适用场景

  • 颗粒材料的多尺度建模
  • 异质复合材料的均匀化
  • 具有微观旋转的材料系统
  • 需要考虑尺寸效应的工程问题

参考文献

论文引用了63篇相关文献，主要包括：

  • Cosserat连续体理论的经典文献
  • 离散元方法和均匀化理论
  • 多尺度建模的相关工作
  • 验证所用的数值方法文献

总体评价: 这是一篇理论严谨、实用价值高的优秀论文，为异质材料的多尺度建模提供了重要的理论工具和数值方法。论文的推导过程清晰，验证充分，对相关领域具有重要的学术和应用价值。