Тройные интегралы

Как ракетные инженеры вычисляют центр масс топливо-загруженного Falcon 9 с точностью 0.5%? Тройным интегралом по сложной 3D-области.

• SpaceX Falcon 9: центр масс и момент инерции топливо-загруженных ступеней считаются как тройные интегралы
• GE Aviation: интегралы напряжений, температуры, воздушного потока по секциям компрессора реактивного двигателя — везде цилиндрические координаты
• Гравитационное моделирование NASA: интеграл плотности массы по нерегулярным астероидам для прогноза траекторий спутников
• ANSYS Fluent CFD: суммарные концентрации в камере сгорания считаются как тройные интегралы по миллионам тетраэдральных ячеек

Предварительные знания

• Двойные интегралы и теорема Фубини
• Полярные координаты и якобиан
• Визуализация 3D-тел

• Двойные интегралы

Тройные интегралы в декартовых координатах

ANSYS Fluent разбивает 3D-камеру сгорания на 10 миллионов тетраэдральных ячеек, вычисляет концентрацию топлива в каждой и суммирует, чтобы получить полную массу — тройной интеграл. Та же конструкция работает для момента инерции лопатки турбины, заряда внутри конденсатора, кинетической энергии вращающегося маховика.

Координаты центра масс: x_bar = (1/m) * тройной интеграл от x rho dV, аналогично для y_bar, z_bar. Моменты инерции: I_z = тройной интеграл от (x^2 + y^2) rho dV. Все стандартные величины механики твёрдого тела — это тройные интегралы.

Нарисуйте тело (или хотя бы его проекцию на xy) перед записью пределов. Самая внутренняя переменная ограничена двумя поверхностями; средняя — двумя кривыми на проекции; внешняя — двумя числами.

Цилиндрические координаты

GE Aviation моделирует компрессор реактивного двигателя как 17 уложенных в стопку цилиндрических секций. Интегралы напряжений, температуры, воздушного потока по каждой секции записаны в цилиндрических координатах, потому что геометрия по природе вращательно-симметрична. Декартовы увеличили бы сложность кода на порядок.

Используйте цилиндрические координаты, когда тело имеет ось симметрии вдоль z: цилиндры, конусы, параболоиды, тороидальные формы. Множитель r в объёмном элементе часто хорошо ладит с подынтегральными функциями вида f(x^2 + y^2, z).

Типичная ошибка: забыть якобиан r. Без него вы вычисляете другой (неверный) интеграл. Правило: каждый раз при переходе в цилиндрические умножайте на r.

Сферические координаты

NASA считает гравитационное притяжение астероида на спутник, интегрируя плотность массы по объёму астероида в сферических координатах. Команда SpaceX Starlink интегрирует усиление антенны по небесной сфере для оптимизации бюджета линии связи спутник-земля. Всё со сферической симметрией — гравитация, электрические поля, радар — просит сферических.

Обозначения разнятся. Многие учебники физики меняют ролями phi и theta. Принятая здесь математическая конвенция: phi — полярный (от оси z), theta — азимутальный. Всегда проверяйте конвенцию в любом тексте.

Если подынтегральная функция зависит от x^2 + y^2 + z^2 = rho^2 или область это шар/конус/сферический слой — переходите в сферические сразу. Подынтегральная функция упрощается, а объёмный элемент впитывает геометрический множитель.

Масса, центр масс, моменты инерции

Когда SpaceX проектирует вторую ступень Falcon 9, ей нужны точный центр масс и момент инерции каждой топливо-загруженной конфигурации. Промах в 2% по центру — ракета теряет управление по тангажу во время горения. Каждое из этих чисел — тройной интеграл по 3D-области, считаемый по данным высокой точности.

Теорема Паппа (I-IV вв.): объём тела вращения равен площади образующей области, умноженной на расстояние, пройденное её центроидом. Тройной интеграл сводится к однострочному вычислению, если геометрия симметрична.

Перед выбором координат проверьте симметрию. Вращательная вокруг z — цилиндрические. Сферическая — сферические. Нет симметрии — декартовы. Правильный выбор сокращает работу вдвое (или удваивает её).

Куда ведут тройные интегралы

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

• Замена переменных — Общая теория якобианов объединяет цилиндрические, сферические и произвольные координатные отображения
• Теорема о дивергенции — Превращает тройной интеграл от div F по телу в поверхностный по границе
• Механика сплошных сред — Напряжения, деформации интегрируются по 3D-телам через тройные интегралы
• Вероятность в R^3 — Совместная плотность трёх случайных величин интегрируется для получения вероятностей

Итоги

• Тройной интеграл = предел 3D-сумм Римана: тройной интеграл от f dV = lim sum f(x*, y*, z*) Delta x Delta y Delta z
• Декартова повторная форма: сначала по z между двумя поверхностями, потом по y между двумя кривыми, потом по x между двумя числами
• Цилиндрические: x = r cos theta, y = r sin theta, z = z; dV = r dr d theta dz — для осевой симметрии
• Сферические: x = rho sin phi cos theta, y = rho sin phi sin theta, z = rho cos phi; dV = rho^2 sin phi d rho d phi d theta — для радиальной симметрии
• Приложения: объём = тройной интеграл от dV; масса = тройной интеграл от rho dV; центр масс = (1/m) * тройной интеграл от x rho dV; момент инерции I_z = тройной интеграл от (x^2 + y^2) rho dV

Связанные уроки

• calc-21-double-integrals — Тройные интегралы расширяют конструкцию двойного на одну размерность
• calc-23-coordinate-changes — Цилиндрические и сферические координаты — рабочие лошадки замены переменных в 3D
• calc-26-divergence-theorem — Теорема о дивергенции превращает тройной интеграл по телу в поверхностный по его границе