Вариационные методы в механике сплошных сред

Вариационные принципы как язык механики

Механика сплошных сред — теория деформируемых тел: балок, пластин, жидкостей, резин. Её уравнения (равновесия, движения) выводятся из вариационных принципов. Это не просто математическое удобство: вариационная формулировка напрямую порождает метод конечных элементов — главный инструмент инженерных расчётов. Каждый раз, когда инженер запускает ANSYS или COMSOL, он решает вариационную задачу.

Принцип виртуальных перемещений

Постановка: тело Ω в равновесии под действием объёмных сил f и поверхностных нагрузок t на части ∂Ω_t.

Принцип виртуальных перемещений: тело в равновесии тогда и только тогда, когда суммарная виртуальная работа равна нулю для любого допустимого виртуального перемещения δu:

δW = ∫_Ω σᵢⱼ δεᵢⱼ dV − ∫Ω fᵢ δuᵢ dV − ∫{∂Ω_t} tᵢ δuᵢ dS = 0

Расшифровка терминов:

σᵢⱼ — тензор напряжений (6 независимых компонент для симметричного σ)
δεᵢⱼ = (δuᵢ,ⱼ + δuⱼ,ᵢ)/2 — виртуальная деформация (симметризованный градиент)
fᵢ — компоненты объёмных сил (например, гравитация, fi = ρgᵢ)
tᵢ = σᵢⱼnⱼ — компоненты поверхностных нагрузок

Это «слабая» формулировка уравнений равновесия. Из него следует «сильная» через интегрирование по частям: ∂ⱼσᵢⱼ + fᵢ = 0 в Ω, σᵢⱼnⱼ = tᵢ на ∂Ω_t.

Принцип минимума потенциальной энергии

Для упругих тел (не диссипативных):

Полная потенциальная энергия: Π[u] = ∫_Ω W(ε(u)) dV − ∫Ω f·u dV − ∫{∂Ω_t} t·u dS

где W(ε) — плотность упругой потенциальной энергии (Гуков материал: W = λ(trε)²/2 + μ tr(ε²)).

Принцип: из всех допустимых полей перемещений (удовлетворяющих кинематическим граничным условиям) равновесное минимизирует Π.

Уравнение ЭЛ для Π: ∂ⱼσᵢⱼ + fᵢ = 0 — уравнения равновесия!

Таким образом, уравнения равновесия механики — это условия первого порядка для минимума потенциальной энергии. Вариационное исчисление «выводит» законы механики.

МКЭ: метод Ритца с конечными элементами

МКЭ — численная реализация метода Ритца для вариационных задач механики.

Идея: разбить Ω на конечные элементы Ωₑ (треугольники, тетраэдры). В каждом элементе аппроксимировать u полиномиальными функциями форм Nᵢ(x): u ≈ Σᵢ uᵢ Nᵢ(x), где uᵢ — узловые перемещения.

Подстановка в Π: Π[u] = Π(u₁,...,uN) — функция от N вещественных чисел.

∂Π/∂uᵢ = 0 → KU = F (система линейных уравнений)

где K = ∫_Ω Bᵀ CB dV — матрица жёсткости (B — матрица деформаций, C — тензор упругости), F — вектор нагрузок.

Размерность: N = 3 × (число узлов). Для модели самолёта — миллионы степеней свободы. ANSYS решает такие системы за часы.

Смешанные принципы

Классический принцип Лагранжа — это «однопольный» (только u). Смешанные принципы вводят несколько полей.

Принцип Хеллингера-Рейсснера (двупольный, u и σ):

Π_HR[u, σ] = ∫_Ω [σ:ε(u) − W*(σ)] dV − ∫Ω f·u dV − ∫{∂Ω_t} t·u dS

где W*(σ) — сопряжённая энергия (преобразование Лежандра от W(ε)).

Уравнения ЭЛ → одновременно уравнения совместности деформаций и равновесия. Удобно при независимой аппроксимации σ и u.

Принцип Кастильяно: для конструкций с известными усилиями — перемещение в точке приложения силы = производная дополнительной энергии по этой силе. Прост для ферм и балок.

Полный разбор: изгиб балки Эйлера-Бернулли

Задача: балка длиной L, закреплена в точке x=0 (защемление: u(0)=u'(0)=0), свободный конец нагружен силой P при x=L. Найти форму изогнутой балки.

Функционал: Π[u] = ∫₀ᴸ (EI/2)(u'')² dx − P·u(L)

где EI — жёсткость балки.

Уравнение ЭЛ: ∂Π/∂(δu) = 0 → EI u'''' = 0 в (0,L).

Граничные условия: u(0) = u'(0) = 0 (защемление), EI u''(L) = 0, EI u'''(L) = P.

Решение: u'''' = 0 → u = ax³ + bx² + cx + d. Из BC: d = 0, c = 0, 6aL + 2b = 0, 6aEI = P.

a = P/(6EI), b = −PL/(2EI). Прогиб: u(x) = Px²(3L−x)/(6EI).

Прогиб на конце: u(L) = PL³/(3EI) — стандартная формула строительной механики!

Нелинейная упругость и гиперупругость

Для больших деформаций (резина, биоткани) линейный закон Гука не работает.

Запас упругости: W = W(F), где F = I + ∇u — градиент деформации. Задача: min ∫_Ω W(F) dV − ∫ f·u dV.

Уравнение ЭЛ: Div P + f = 0, P = ∂W/∂F — Пиола-Кирхгофф.

Модели: неогуковский (W = μ(I₁−3)/2), Муни-Ривлин (W = C₁(I₁−3) + C₂(I₂−3)). Используются в биомеханике (моделирование сердца, сосудов, кожи), резиноизделиях, мягкой робототехнике.

Принцип минимума потенциальной энергии

В статике твёрдого деформируемого тела фундаментальный принцип: равновесная конфигурация минимизирует полную потенциальную энергию Π = U_внутр − W_внешн, где U_внутр — упругая энергия деформации, W_внешн — работа внешних сил. Это вариационная задача: минимум Π[u] по полю смещений u(x), удовлетворяющему граничным условиям. Уравнение Эйлера-Лагранжа этого функционала — уравнения равновесия упругого тела.

Уравнения Эйлера для жидкости и Навье-Стокса

Движение идеальной (невязкой) жидкости описывается уравнениями Эйлера: ∂v/∂t + (v·∇)v = −∇p/ρ. Они выводятся из принципа Гамильтона с лагранжианом L = (ρ/2)|v|². Для вязкой жидкости — уравнения Навье-Стокса с дополнительным членом νΔv. Существование и гладкость решений Навье-Стокса — одна из задач тысячелетия (1 миллион долларов от Института Клэя).

Метод конечных элементов (МКЭ)

МКЭ — практический инструмент решения вариационных задач механики сплошных сред:

Разбить тело на конечные элементы (треугольники, тетраэдры, гексаэдры)
Аппроксимировать поле смещений в каждом элементе линейной/квадратичной функцией базиса
Свести вариационную задачу к системе линейных уравнений KU = F (K — матрица жёсткости, F — вектор внешних сил)
Решить и постобработать (вычислить напряжения, деформации)

Промышленные пакеты МКЭ: ANSYS, ABAQUS, COMSOL, NASTRAN. Современные расширения: изогеометрический анализ (NURBS-базис), мешфри-методы (SPH, MPM для пластических деформаций).

Принцип виртуальных работ

Альтернативная формулировка: при равновесии работа всех внешних и внутренних сил на любом виртуальном перемещении равна нулю. Это даёт «слабую формулировку» вариационной задачи — основу МКЭ.

Применения

Гражданское строительство: расчёт мостов, небоскрёбов, плотин (Бурдж-Халифа, мост Акаси-Кайкё)
Авиация: расчёт несущих конструкций самолёта (Boeing 787, Airbus A350)
Автомобилестроение: краш-тесты численным моделированием
Биомеханика: моделирование костей, имплантатов, кровотока
Геофизика: моделирование тектоники, землетрясений