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. 기존 방법의 한계

점근 전개 균질화: 대표 체적 단위가 0으로 수렴한다고 가정하여 유한 크기 효과를 무시함
고전적 Love-Weber 공식: 주로 응력 텐서에 적용되며, 우력 응력 텐서의 처리가 불완전함
경계 조건 처리: 기존 방법은 경계 효과 및 기준점 선택 처리에 부족함

4. 연구 동기

본 논문은 다음을 수행할 수 있는 더욱 완전하고 엄밀한 이론적 틀을 확립하는 것을 목표로 한다:

유한 크기의 대표 체적 단위 처리
내력 및 외력 표현 방식 동시 고려
우력 응력 텐서에서 기준점 위치의 영향 명확화

핵심 기여

가상 일 등가성을 기반으로 한 엄밀한 유도 틀 확립, 연속 및 이산 이질 시스템에 적용 가능
유한 RVE 크기를 고려한 거시 량 표현식 유도, 응력 텐서, 우력 응력 텐서 및 플럭스 벡터 포함
두 가지 변형 형태 제시: 내력 기반 및 외력 기반 표현
우력 응력 텐서의 추가 항 발견, 기준계 및 거시 점 위치에 대한 의존성 설명
완전한 검증 제공, 정상 및 과도 조건에서의 해석해 및 수치 모델 비교

방법 상세 설명

작업 정의

본 논문의 핵심 작업은 미시 규모 이질 재료와 거시 규모 균질 연속체 간의 등가 관계를 확립하는 것이며, 구체적으로 다음을 포함한다:

입력: 미시 규모의 응력장 σ^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₂)