Вариационные задачи с ограничениями и связями

Ограниченные системы в физике и механике

В реальном мире движение тел почти всегда ограничено: маятник движется по дуге окружности, автомобиль катится без проскальзывания, жидкость течёт по трубе. Такие ограничения называются связями (constraints). Вариационное исчисление с ограничениями — это теория, позволяющая системно работать с такими задачами. Разница между голономными (конечными) и неголономными (дифференциальными) связями принципиальна и имеет глубокие физические следствия.

Голономные связи

Голономная связь задаёт условие вида φ(x, y(x)) = 0 — это функциональное уравнение на кривые. Называется так от греческого «holos» (целый) — ограничение полностью определяет конфигурацию.

Пример: математический маятник в декартовых координатах. Частица движется в ℝ², но прикреплена к точке нитью длиной l. Голономная связь: x² + y² = l². Система имеет 2 − 1 = 1 степень свободы. В обобщённой координате θ (угол): x = l sin θ, y = −l cos θ, связь исчезает!

Метод решения: либо явная параметризация (как с маятником), либо метод множителей Лагранжа-функций.

Неголономные связи

Неголономная связь задаёт дифференциальное условие: g(x, y, y') = 0. Её нельзя «интегрировать» до конечного соотношения — она накладывает ограничение на скорости, а не на конфигурации.

Классический пример: качение монеты без скольжения. Монета катится по плоскости без проскальзывания. Скорость точки контакта = 0: ẋ − r θ̇ cos φ = 0, ẏ − r θ̇ sin φ = 0. Это два дифференциальных условия на 4 переменных. Нельзя свести к меньшему числу конфигурационных переменных!

Для задачи min J = ∫F dx при g(x,y,y') = 0 вводим функцию-множитель Лагранжа λ(x):

∫[F + λ(x)g(x,y,y')] dx

и уравнение ЭЛ для вариации по y даёт систему, определяющую y(x) и λ(x).

Теорема Каратеодори о полных вариационных задачах

Общая постановка с n переменными yᵢ(x) и m ограничениями Gⱼ(x,y,y') = 0:

Необходимые условия оптимальности:

1. Уравнения ЭЛ (модифицированные): Fᵧᵢ − d/dx F_{yᵢ'} + Σⱼ λⱼ(x)(Gⱼ){yᵢ} − d/dx Σⱼ λⱼ(Gⱼ){yᵢ'} = 0

2. Ограничения: Gⱼ(x, y, y') = 0 вдоль траектории

3. Условие поперечности на свободных концах

4. Условия Лежандра-Якоби для достаточности

Вариационное неравенство

В задачах с односторонними ограничениями (y ≥ 0 — «препятствие» снизу) допустимое множество — не линейное пространство, а выпуклый конус.

Условие экстремума: вариационное неравенство:

δJ[y*; η − y*] ≥ 0 для всех допустимых η

Если y* в interior → δJ = 0 (обычное уравнение ЭЛ). Если y* на границе (y* = 0) → неравенство строгое.

Задача с препятствием (obstacle problem): найти наименьшую поверхность, лежащую выше заданного препятствия ψ(x,y). Решение удовлетворяет:

−∆u ≥ 0 в Ω, u ≥ ψ в Ω, (−∆u)(u − ψ) = 0

Это система нелинейных условий — сильное неравенство. Применяется в финансовой математике (цена американского опциона) и теории упругопластичности.

Полный разбор: маятник через метод Лагранжа

Задача: частица массы m прикреплена нитью длиной l к точке O = (0,0). Движется в силе тяжести. Найти уравнения движения.

Переменные: x(t), y(t) — координаты. Связь: φ = x² + y² − l² = 0 (нить нерастяжима).

Лагранжиан: L = T − U = m(ẋ² + ẏ²)/2 − mgy. Действие: S = ∫L dt.

Уравнения с множителем: уравнение Лагранжа с реакцией связи:

mẍ = 2λx, mÿ = 2λy − mg, x² + y² = l²

Из связи дифференцированием дважды: xẍ + yÿ + ẋ² + ẏ² = 0.

Подставляя: x(2λx/m) + y(2λy/m − g) + v² = 0 → 2λl²/m − gy + v² = 0 → λ = m(gy − v²)/(2l²).

Реакция нити: T = −2λ = m(v² − gy)/l² (центростремительное ускорение + вес по нормали).

Переход к обобщённой координате: θ (угол). x = l sin θ, y = −l cos θ. L = ml²θ̇²/2 + mgl cos θ. ЭЛ: ml²θ̈ = −mgl sin θ → θ̈ + (g/l) sin θ = 0 — уравнение маятника. Связь исчезла, λ исключилась!

Связи Лагранжа vs неголономные связи

В вариационных задачах с ограничениями различают два типа связей:

• Голономные: φ(x, y) = 0 — связь только на координатах. Можно явно разрешить и подставить, понизив размерность.
• Неголономные: φ(x, y, y') = 0 — связь включает производные и в общем случае не интегрируется. Классический пример: качение шара без проскальзывания (5 голономных + 2 неголономных связи). Эти связи требуют принципиально других методов: множители Лагранжа применимы, но требуют осторожности.

Принцип Гамильтона-Понтрягина

Для задач с ограничениями на управление u ∈ U Понтрягин (1956) сформулировал принцип максимума: оптимальное управление в каждый момент времени максимизирует «гамильтониан» H = pᵀf − F. Это ослабленная версия уравнения Эйлера-Лагранжа, допускающая разрывные управления (бэнг-бэнг). Принцип максимума — стандартный инструмент в оптимальном управлении, заменяющий вариационное исчисление при наличии ограничений на управление.

Численные методы

Задачи с ограничениями требуют специализированных численных методов:

• Прямое преобразование в NLP: дискретизация состояния и управления, добавление ограничений как нелинейных условий, решение через IPOPT, SNOPT
• Метод проекции: после каждого шага градиентного спуска проектировать решение на множество допустимых функций
• Augmented Lagrangian: штрафная функция плюс множители — устойчивое решение задач с активными связями

Применения

• Робототехника: движение манипулятора с механическими связями (шарниры, контакты)
• Транспортное планирование: учёт ограничений на ускорение, рывок (jerk), угол поворота
• Аэрокосмическая отрасль: оптимизация траектории при ограничениях на тягу, нагрев, перегрузки
• Биомеханика: моделирование походки человека как оптимизация энергии при анатомических ограничениях суставов
• Управление химическими процессами: ограничения на температуру, давление, концентрацию реагентов — стандартная постановка вариационной задачи с активными связями