Макроскопические тензоры напряжений, парных напряжений и потоков, полученные через энергетическую эквивалентность из микроскопических непрерывных и дискретных неоднородных конечных представительных объемов

Основная информация

• 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
• Ссылка на статью: https://arxiv.org/abs/2412.12616

Аннотация

В статье предложен строгий метод вывода для оценки макроскопических тензоров напряжений, парных напряжений и векторов потоков в непрерывных и дискретных неоднородных системах с независимыми перемещениями и поворотами. В отличие от классического асимптотического разложения гомогенизации, работа рассматривает представительные объемные элементы (РОЭ) конечного размера. Сначала выводятся макроскопические величины для неоднородного континуума Коссера, затем уравнения континуума дискретизируются для получения макроскопических величин в дискретных неоднородных системах, и, наконец, выражения для дискретной системы переводятся заново непосредственно из дискретных свойств.

Исследовательский контекст и мотивация

1. Решаемые проблемы

Статья решает проблему точного извлечения макроскопических механических и транспортных свойств из микроскопических неоднородных материальных систем, в частности:

• Как вывести макроскопический тензор напряжений из РОЭ конечного размера
• Как обрабатывать тензоры парных напряжений в системах с независимыми степенями свободы вращения
• Как установить энергетическую эквивалентность между микроскопической и макроскопической шкалами

2. Значимость проблемы

Многие инженерные материалы (такие как бетон, керамика, гранулированные материалы и т.д.) демонстрируют значительную неоднородность в мезо- или микромасштабе. Точное многомасштабное моделирование важно для:

• Прогнозирования общего механического поведения материалов
• Разработки надежных численных моделей
• Установления феноменологических определяющих соотношений

3. Ограничения существующих методов

• Асимптотическое разложение гомогенизации: предполагает, что представительный объемный элемент стремится к нулю, игнорируя эффекты конечного размера
• Классическая формула Лава-Вебера: в основном применима к тензору напряжений, недостаточно развита для обработки тензоров парных напряжений
• Обработка граничных условий: существующие методы имеют недостатки в обработке граничных эффектов и выборе опорной точки

4. Исследовательская мотивация

Статья направлена на установление более полной и строгой теоретической базы, способной:

• Обрабатывать представительные объемные элементы конечного размера
• Одновременно рассматривать внутренние и внешние силовые формулировки
• Явно определить влияние положения опорной точки в тензоре парных напряжений

Основные вклады

1. Установлена строгая базовая схема вывода на основе принципа виртуальной работы, применимая к непрерывным и дискретным неоднородным системам
2. Выведены выражения макроскопических величин с учетом конечного размера РОЭ, включая тензоры напряжений, парных напряжений и векторы потоков
3. Предложены две вариантные формулировки: на основе внутренних сил и на основе внешних сил
4. Обнаружены дополнительные члены в тензоре парных напряжений, объясняющие его зависимость от системы отсчета и положения макроскопической точки
5. Предоставлена полная верификация, включая сравнение аналитических решений и численных моделей в стационарных и нестационарных условиях

Подробное описание методологии

Определение задачи

Основная задача работы состоит в установлении отношения эквивалентности между микроскопическим неоднородным материалом и макроскопическим однородным континуумом, конкретно включая:

• Входные данные: микроскопические поля напряжений σ^mic, поля парных напряжений μ^mic, поля потоков a^mic
• Выходные данные: макроскопический тензор напряжений σ^mac, тензор парных напряжений μ^mac, вектор потока a^mac
• Ограничения: удовлетворение принципу виртуальной работы и уравнениям равновесия

Теоретическая база

1. Основы континуума Коссера

Основные переменные континуума Коссера включают:

• Вектор перемещения 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. Формулировка внутренних сил (тип Лава-Вебера)

Через равновесие сил и моментов получается обобщение знаменитой формулы Лава-Вебера:

σ^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. Верификация задачи Пуассона

• Геометрия: квадратная область 1×1 м²
• Дискретизация: примерно 48 000 узлов, 72 000 трубчатых элементов
• Тестовые случаи:
  • Тест постоянного потока: p = 2(2x-1) + 2(2y-1), a = (-4,-4)^T
  • Нелинейное аналитическое решение: p = 16x₁x₂(1-x₁)(1-x₂)

2. Верификация механической задачи

• Геометрия: консольная балка, глубина D=3м, длина S=12м
• Дискретизация: 1 034 196 степеней свободы
• Параметры материала: E₀=40ГПа, α=0.3, β=0~10⁵
• Сравниваемые модели:
  • Дискретная модель
  • Гомогенизированная модель Коссера
  • Упрощенная модель Коссера
  • Континуум Коши

Показатели оценки

• Оценка точности: отклонение от аналитического решения
• Сходимость: результаты при различных размерах сетки
• Согласованность: сравнение формулировок внутренних и внешних сил
• Граничные эффекты: влияние граничного радиусного зазора

Результаты экспериментов

Основные результаты

1. Верификация задачи Пуассона

• Вариант I (с граничным радиусным зазором): всегда дает точное решение
• Вариант II (без граничного радиусного зазора): точность снижается с уменьшением размера сетки
• Сходимость: при больших размерах сетки результаты обоих вариантов близки

2. Верификация механической задачи

Стационарный анализ:

• Дискретная модель хорошо согласуется с гомогенизированной моделью
• С увеличением жесткости на изгиб β напряжения уменьшаются, а парные напряжения увеличиваются
• При β<10³ эффекты Коссера можно пренебречь, можно использовать континуум Коши

Нестационарный анализ:

• Отклик перемещения всех моделей высоко согласован
• Временная эволюция напряжений и парных напряжений соответствует ожиданиям
• Подтверждена справедливость формул в нестационарных условиях

Ключевые открытия

1. Эффект граничного радиусного зазора

Эксперименты подтвердили важность граничного радиусного зазора:

• Для контрольных объемов малого размера этот член нельзя пренебречь
• В системах с большим количеством частиц влияние этого члена меньше
• Обеспечена количественная оценка ошибок

2. Влияние положения опорной точки

Подтверждены теоретические предсказания члена положения опорной точки в тензоре парных напряжений:

• Когда макроскопическая точка находится в центре масс области, дополнительный член исчезает
• В других положениях возникает дополнительный вклад парных напряжений

3. Влияние параметров материала

• Параметр β значительно влияет на характеристическую длину Коссера ℓ
• При β=0 ℓ≈4.6мм, при β=10⁵ ℓ≈0.7м
• Параметры Ламе λ,μ относительно стабильны

Связанные работы

Классические теоретические основы

• Формула Лава-Вебера: классическое представление тензора напряжений в дискретных системах
• Асимптотическое разложение гомогенизации: теория бесконечно малых РОЭ Фореста и др.
• Принцип виртуальной работы: метод энергетической эквивалентности Чанга и Куна

Инновации в данной работе

1. РОЭ конечного размера: рассмотрение эффектов размера в отличие от асимптотической теории
2. Полная обработка парных напряжений: включение вклада микроскопических напряжений на расстояниях
3. Унифицированная база: одновременная обработка непрерывных и дискретных систем
4. Две формулировки: полное соответствие между внутренними и внешними силами

Выводы и обсуждение

Основные выводы

1. Теоретический вклад: установлена строгая теория гомогенизации для РОЭ конечного размера, обобщена классическая формула Лава-Вебера
2. Методологическая инновация: предложен полный метод обработки зависимости тензора парных напряжений от опорной точки
3. Практическая ценность: предоставлены более точные инструменты для многомасштабного моделирования

Ограничения

1. Выбор виртуальных режимов деформации: некоторые высокопорядковые балансы виртуальной работы игнорируются
2. Ограничения граничных условий: требуются статические граничные условия на всех границах
3. Вычислительная сложность: граничный радиусный зазор увеличивает вычислительную нагрузку

Направления будущих исследований

1. Нелинейное расширение: обобщение на нелинейные определяющие соотношения
2. Связанные задачи: расширение на механико-транспортные связанные задачи
3. Высокопорядковая теория: рассмотрение более высокопорядковых виртуальных режимов деформации

Глубокая оценка

Преимущества

1. Теоретическая строгость: строгий математический вывод на основе принципа виртуальной работы
2. Высокая практичность: прямое применение к существующим дискретным и непрерывным моделям
3. Полная верификация: включает аналитические решения, численные сравнения и исследования параметров
4. Ясное изложение: логичная структура, точная математическая нотация

Недостатки

1. Область применения: в основном ориентирована на линейно-упругие задачи, нелинейное расширение требует дальнейших исследований
2. Вычислительная эффективность: полные формулы имеют высокую вычислительную стоимость
3. Чувствительность параметров: некоторые результаты очень чувствительны к параметру материала β

Влияние

1. Академическая ценность: предоставляет важную теоретическую базу для многомасштабной механики
2. Инженерные приложения: может быть непосредственно применена к моделированию бетона, грунтов и других неоднородных материалов
3. Реализация программного обеспечения: предоставлен открытый исходный код, облегчающий распространение

Применимые сценарии

• Многомасштабное моделирование гранулированных материалов
• Гомогенизация неоднородных композитных материалов
• Материальные системы с микроскопическими поворотами
• Инженерные задачи, требующие учета эффектов размера

Библиография

Статья цитирует 63 связанные работы, в основном включая:

• Классические работы по теории континуума Коссера
• Метод дискретных элементов и теория гомогенизации
• Связанные работы по многомасштабному моделированию
• Литература по численным методам, используемым для верификации

Общая оценка: Это отличная статья с строгой теорией и высокой практической ценностью, предоставляющая важные теоретические инструменты и численные методы для многомасштабного моделирования неоднородных материалов. Процесс вывода ясен, верификация полна, работа имеет важную академическую и прикладную ценность для соответствующих областей.