9. Поверхностные интегралы

Как считать площадь кривой поверхности, массу тонкой оболочки и поток векторного поля через поверхность. Всё сводится к двойному интегралу по проекции, главное — правильно найти нормаль и элемент площади $dS$.

пример КР2 в РПД, задание 5ДЗ4, задача 1РГР2, задание 5aэкзамен: вопросы 27–30, 40, 41

Откуда известно, что это нужно

По РПД курса (рабочая программа):

Теоремы Стокса и Остроградского-Гаусса, дивергенция и ротор — в главе 10. Здесь они упоминаются только как проверка.

Нужно уверенно уметь: двойной интеграл и переход в полярные координаты (глава 6), нормаль к поверхности через градиент (глава 4). Полезна аналогия с криволинейными интегралами: там $dl$ — элемент длины кривой, здесь $dS$ — элемент площади поверхности; интеграл 1 рода не зависит от направления, 2 рода — зависит.

1. Поверхности, нормали, стороны

вопрос экзамена №27вопрос экзамена №28

Три способа задать поверхность

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

ЗаданиеПримерВектор нормали $\vec N$ (не единичный)
Явно: $z = f(x,y)$, $(x,y)\in D$параболоид $z = x^2+y^2$$\vec N = (-f'_x,\ -f'_y,\ 1)$ — смотрит вверх ($\cos\gamma > 0$); противоположный $(f'_x, f'_y, -1)$ — вниз
Неявно: $F(x,y,z) = 0$сфера $x^2+y^2+z^2 = R^2$$\vec N = \grad F = (F'_x,\ F'_y,\ F'_z)$ или $-\grad F$
Параметрически: $\vec r(u,v) = (x(u,v),\ y(u,v),\ z(u,v))$, $(u,v)\in\Omega$цилиндр $(R\cos\varphi,\ R\sin\varphi,\ z)$$\vec N = \vec r'_u\times\vec r'_v$ или $\vec r'_v\times\vec r'_u$

Явное задание — частный случай параметрического: $\vec r(x,y) = (x,\ y,\ f(x,y))$, тогда $\vec r'_x = (1, 0, f'_x)$, $\vec r'_y = (0, 1, f'_y)$ и

$$\vec r'_x\times\vec r'_y = \begin{vmatrix}\vec i & \vec j & \vec k\\ 1 & 0 & f'_x\\ 0 & 1 & f'_y\end{vmatrix} = (-f'_x,\ -f'_y,\ 1).$$

Явное задание — и частный случай неявного: $F = f(x,y) - z$, $\grad F = (f'_x, f'_y, -1)$, тот же вектор с минусом. Так что формулы не надо зубрить по отдельности.

Определение (гладкая поверхность)

Поверхность $\vec r(u,v)$, $(u,v)\in\Omega$, называется гладкой, если $\vec r'_u$, $\vec r'_v$ непрерывны и $\vec r'_u\times\vec r'_v \ne \vec 0$ (касательные векторы не параллельны). Тогда в каждой точке есть касательная плоскость и нормаль, меняющаяся непрерывно. Кусочно-гладкая поверхность склеена из конечного числа гладких кусков (поверхность куба, конуса с дном, цилиндра с крышками).

Векторное произведение: что нужно помнить

Для $\vec a = (a_1,a_2,a_3)$, $\vec b = (b_1,b_2,b_3)$:

$$\vec a\times\vec b = \begin{vmatrix}\vec i & \vec j & \vec k\\ a_1 & a_2 & a_3\\ b_1 & b_2 & b_3\end{vmatrix} = (a_2b_3 - a_3b_2,\ \ a_3b_1 - a_1b_3,\ \ a_1b_2 - a_2b_1).$$

Треугольник $ABC$ — половина параллелограмма, натянутого на $\overrightarrow{AB}$ и $\overrightarrow{AC}$ (диагональ $BC$ делит его на два равных треугольника), поэтому $S_{ABC} = \frac12|\overrightarrow{AB}\times\overrightarrow{AC}|$.

Пример 1 (явное задание)

Найти единичную нормаль к параболоиду $z = x^2+y^2$ в точке $(1,1,2)$, направленную вверх (верхняя сторона).

Решение

$f'_x = 2x = 2$, $f'_y = 2y = 2$. Вектор, смотрящий вверх: $\vec N = (-f'_x, -f'_y, 1) = (-2, -2, 1)$, $|\vec N| = \sqrt{4+4+1} = 3$.

$\vec n = \left(-\frac23,\ -\frac23,\ \frac13\right)$. Проверка «вверх»: третья координата $\frac13 > 0$, то есть угол с $Oz$ острый.

Нижняя сторона в той же точке: $\vec n = \left(\frac23, \frac23, -\frac13\right)$.

Пример 2 (неявное задание)

Найти внешнюю единичную нормаль к сфере $x^2+y^2+z^2 = 9$ в точке $(1,2,2)$.

Решение

$F = x^2+y^2+z^2-9$, $\grad F = (2x, 2y, 2z) = (2,4,4)$ в точке. Градиент смотрит в сторону роста $F$, то есть наружу (снаружи сферы $F > 0$). $|\grad F| = 6$.

$\vec n = \left(\frac13, \frac23, \frac23\right)$. Для сферы с центром в начале координат всегда $\vec n_{\text{внеш}} = \frac{(x,y,z)}{R}$ — радиус-вектор, делённый на радиус.

Пример 3 (параметрическое задание: цилиндр)

Цилиндр $x^2+y^2 = R^2$ задан как $\vec r(\varphi, z) = (R\cos\varphi,\ R\sin\varphi,\ z)$. Найти $\vec r'_\varphi\times\vec r'_z$ и его длину.

Решение

$\vec r'_\varphi = (-R\sin\varphi,\ R\cos\varphi,\ 0)$, $\vec r'_z = (0, 0, 1)$.

$$\vec r'_\varphi\times\vec r'_z = \begin{vmatrix}\vec i & \vec j & \vec k\\ -R\sin\varphi & R\cos\varphi & 0\\ 0 & 0 & 1\end{vmatrix} = (R\cos\varphi,\ R\sin\varphi,\ 0).$$

Это $(x, y, 0)$ — горизонтальный вектор от оси наружу: внешняя нормаль. Длина $R$. Третья координата нулевая: у цилиндра с вертикальными образующими нормаль горизонтальна, это понадобится в разделе 4.

Пример 4 (параметрическое задание: сфера)

Сфера $\vec r(\theta,\varphi) = (R\sin\theta\cos\varphi,\ R\sin\theta\sin\varphi,\ R\cos\theta)$, $\theta\in[0,\pi]$, $\varphi\in[0,2\pi]$ (сферические координаты при $\rho = R$). Найти $\vec r'_\theta\times\vec r'_\varphi$.

Решение

$\vec r'_\theta = (R\cos\theta\cos\varphi,\ R\cos\theta\sin\varphi,\ -R\sin\theta)$, $\vec r'_\varphi = (-R\sin\theta\sin\varphi,\ R\sin\theta\cos\varphi,\ 0)$.

Компоненты: $x$: $R\cos\theta\sin\varphi\cdot 0 - (-R\sin\theta)\cdot R\sin\theta\cos\varphi = R^2\sin^2\theta\cos\varphi$; $y$: $(-R\sin\theta)(-R\sin\theta\sin\varphi) - R\cos\theta\cos\varphi\cdot 0 = R^2\sin^2\theta\sin\varphi$; $z$: $R\cos\theta\cos\varphi\cdot R\sin\theta\cos\varphi + R\cos\theta\sin\varphi\cdot R\sin\theta\sin\varphi = R^2\sin\theta\cos\theta$.

$$\vec r'_\theta\times\vec r'_\varphi = R\sin\theta\cdot(R\sin\theta\cos\varphi,\ R\sin\theta\sin\varphi,\ R\cos\theta) = R\sin\theta\cdot\vec r,\qquad |\vec r'_\theta\times\vec r'_\varphi| = R^2\sin\theta.$$

При $0 < \theta < \pi$ множитель $R\sin\theta > 0$, так что это внешняя нормаль. На полюсах ($\theta = 0, \pi$) произведение нулевое — параметризация там вырождается, но это две точки, на интеграл они не влияют.

Стороны поверхности и ориентация

У листа бумаги две стороны: можно покрасить одну в красный, другую в синий, и краски не встретятся. Математически «сторона» — это выбор одного из двух направлений нормали в каждой точке, причём выбор должен меняться непрерывно.

Определение

Возьмём точку $M$ на поверхности и выберем в ней одно из двух направлений единичной нормали. Будем двигать точку по замкнутому контуру на поверхности (не пересекая её края), непрерывно перенося нормаль. Если для любого такого контура нормаль возвращается в исходное положение, поверхность двусторонняя. Если есть контур, после обхода которого нормаль оказывается противоположной, поверхность односторонняя.

Сторона двусторонней поверхности — это непрерывное поле единичных нормалей $\vec n(M)$ на ней. Выбрать сторону = ориентировать поверхность. У двусторонней связной поверхности ровно две стороны.

Лист Мёбиуса

Стандартный пример односторонней поверхности: полоску бумаги перекручивают на пол-оборота и склеивают концы. Параметризация:

$$\vec r(u,v) = \left((1 + v\cos\tfrac u2)\cos u,\ \ (1 + v\cos\tfrac u2)\sin u,\ \ v\sin\tfrac u2\right),\quad u\in[0,2\pi],\ v\in[-0.4,\ 0.4].$$

При $v = 0$ получается окружность $x^2+y^2 = 1$ — средняя линия ленты. Параметр $u$ — угол обхода, $v$ — смещение поперёк ленты. Пока $u$ проходит полный круг, отрезок поперёк ленты поворачивается на угол $\frac u2$, то есть на пол-оборота.

Пример 5

Проверить вычислением, что лист Мёбиуса односторонний: найти нормаль $\vec r'_u\times\vec r'_v$ вдоль средней линии $v = 0$ и сравнить её при $u = 0$ и $u = 2\pi$.

Решение

При $v = 0$: $\vec r'_u = (-\sin u,\ \cos u,\ 0)$, $\vec r'_v = (\cos\frac u2\cos u,\ \cos\frac u2\sin u,\ \sin\frac u2)$.

$$\vec r'_u\times\vec r'_v = \left(\cos u\sin\tfrac u2,\ \ \sin u\sin\tfrac u2,\ \ -\cos\tfrac u2\right).$$

Это непрерывная функция $u$, и её длина всегда 1 ($\cos^2 u\sin^2\frac u2 + \sin^2u\sin^2\frac u2 + \cos^2\frac u2 = 1$), то есть мы честно переносим нормаль вдоль средней линии.

При $u = 0$: точка $(1,0,0)$, нормаль $(0, 0, -1)$. При $u = 2\pi$: та же точка $(1,0,0)$, нормаль $(0, 0, -\cos\pi) = (0, 0, 1)$. Обошли контур — нормаль перевернулась. Значит, лист Мёбиуса односторонний.

Ещё одно наблюдение: $\vec r(2\pi, v) = (1 - v,\ 0,\ 0) = \vec r(0, -v)$ — края полоски склеены «наоборот», это и есть перекрут.

Зачем это нужно

Поверхностный интеграл 2 рода (поток, раздел 4) определяется для выбранной стороны: поток «через поверхность в сторону нормали». У листа Мёбиуса стороны не выбрать, поэтому поток через него не определён. Интегралы 1 рода (площадь, масса) от стороны не зависят, их можно считать по любой гладкой поверхности.

Пример Шварца: почему площадь не определяют через вписанные многогранники

Длину кривой определяют как предел длин вписанных ломаных (глава 8). Хочется так же определить площадь поверхности: вписать многогранник с мелкими гранями и перейти к пределу. Шварц показал, что так нельзя: даже для обычного цилиндра предел зависит от того, как именно измельчать грани, и может быть бесконечным.

Конструкция («сапог Шварца», «фонарь Шварца»). Цилиндр радиуса $R$ и высоты $H$ режем горизонтальными плоскостями на $m$ слоёв. Каждую из $m+1$ окружностей делим на $n$ равных частей, причём на соседних окружностях точки сдвинуты на полшага. Соединяем точки треугольниками: две соседние точки одной окружности и точку между ними на соседней окружности. Получается вписанный многогранник из $2nm$ одинаковых треугольников, все его вершины лежат на цилиндре.

Площадь. Основание треугольника — хорда $2R\sin\frac{\pi}{n}$. Вершина лежит на соседней окружности напротив середины хорды: по вертикали до неё $\frac Hm$, по горизонтали $R - R\cos\frac\pi n$ (середина хорды ближе к оси, чем окружность). Высота треугольника $\sqrt{\left(\frac Hm\right)^2 + R^2\left(1-\cos\frac\pi n\right)^2}$, и всего

$$S(n,m) = 2nm\cdot\frac12\cdot 2R\sin\frac\pi n\cdot\sqrt{\frac{H^2}{m^2} + R^2\Bigl(1-\cos\frac\pi n\Bigr)^2} = 2nR\sin\frac\pi n\cdot\sqrt{H^2 + m^2R^2\Bigl(1-\cos\frac\pi n\Bigr)^2}.$$

При $n\to\infty$ множитель $2nR\sin\frac\pi n\to 2\pi R$, а $1-\cos\frac\pi n\approx\frac{\pi^2}{2n^2}$, так что под корнем стоит $H^2 + \frac{\pi^4R^2}{4}\cdot\frac{m^2}{n^4}$. Всё решает, как быстро растёт $m$ по сравнению с $n^2$:

Как измельчаем$\frac{m}{n^2}$$\lim S(n,m)$
$m = n$$\to 0$$2\pi RH$ — правильная площадь
$m = n^2$$= 1$$2\pi R\sqrt{H^2 + \frac{\pi^4R^2}{4}}$ — больше правильной
$m = n^3$$\to\infty$$+\infty$

При этом в любом случае все грани становятся сколь угодно мелкими и все вершины лежат на цилиндре.

Как понимать

Если слоёв гораздо больше, чем точек на окружности, слои тонкие, а треугольники всё равно уходят внутрь на $R(1-\cos\frac\pi n)$. Многогранник превращается в «гармошку»: грани почти перпендикулярны поверхности цилиндра, и их площадь не похожа на площадь куска цилиндра под ними. Для кривых такого не бывает: маленькая хорда всегда почти параллельна касательной. Для поверхностей маленькая грань может стоять под большим углом к касательной плоскости.

Вывод: площадь поверхности определяют иначе — через касательные плоскости, то есть формулой с $|\vec r'_u\times\vec r'_v|$ (раздел 2). Можно и через вписанные многогранники, но тогда надо требовать, чтобы грани приближали касательные плоскости.

У какой поверхности нельзя выбрать сторону?
Лист Мёбиуса односторонний: при обходе средней линии нормаль переворачивается (пример 5). У сферы и цилиндра есть внешняя и внутренняя стороны, у графика — верхняя и нижняя.

2. Площадь поверхности

вопрос экзамена №29пример КР2 в РПД, задание 5

Откуда формула

Разобьём область параметров $\Omega$ на маленькие прямоугольники $du\times dv$. Каждый переходит на поверхности в маленький криволинейный четырёхугольник. При малых $du$, $dv$ он почти совпадает с параллелограммом, натянутым на векторы $\vec r'_u\,du$ и $\vec r'_v\,dv$ (они лежат в касательной плоскости). Площадь параллелограмма — модуль векторного произведения:

$$dS = |\vec r'_u\,du\times\vec r'_v\,dv| = |\vec r'_u\times\vec r'_v|\,du\,dv.$$ du dv плоскость (u, v) r(u,v) r′u du r′v dv dS ≈ |r′u × r′v| du dv
Площадь поверхности

Параметрически $\vec r(u,v)$, $(u,v)\in\Omega$:

$$S = \iint\limits_\Omega |\vec r'_u\times\vec r'_v|\,du\,dv.$$

Явно $z = f(x,y)$, $(x,y)\in D$ ($D$ — проекция поверхности на плоскость $xOy$):

$$S = \iint\limits_D \sqrt{1 + (f'_x)^2 + (f'_y)^2}\,dx\,dy.$$

Аналогично для $x = g(y,z)$: $\sqrt{1 + (g'_y)^2 + (g'_z)^2}\,dy\,dz$ по проекции на $yOz$; для $y = h(x,z)$ — по проекции на $xOz$.

Неявно $F(x,y,z) = 0$, если поверхность однозначно проектируется на $xOy$: $dS = \dfrac{|\grad F|}{|F'_z|}\,dx\,dy$.

Смысл корня

$\sqrt{1 + f_x'^2 + f_y'^2} = |\vec N|$ для $\vec N = (-f'_x, -f'_y, 1)$, а $\cos\gamma = \frac{1}{|\vec N|}$ — косинус угла между нормалью и осью $Oz$. Значит, $dS = \frac{dx\,dy}{|\cos\gamma|}$: кусочек поверхности, наклонённый на угол $\gamma$, проектируется на $xOy$ с уменьшением площади в $|\cos\gamma|$ раз. Горизонтальная плоскость ($f'_x = f'_y = 0$) — корень равен 1, площадь равна площади проекции. Чем круче поверхность, тем больше корень.

Если поверхность где-то вертикальна ($\cos\gamma = 0$), корень бесконечен, а проекция вырождается в линию — такую поверхность на $xOy$ не проектируют.

Алгоритм (явное задание)
  1. Нарисовать поверхность и найти её проекцию $D$ на удобную координатную плоскость. Поверхность должна проектироваться однозначно (каждая точка $D$ — одна точка поверхности), иначе разбить на куски.
  2. Выразить поверхность как функцию двух координат этой плоскости, найти частные производные.
  3. Посчитать $\sqrt{1 + f_x'^2 + f_y'^2}$ и упростить с учётом уравнения поверхности (для сферы и конуса упрощается сильно).
  4. Посчитать двойной интеграл по $D$; для круглых $D$ — в полярных координатах.
Пример 6

Найти площадь сферы радиуса $R$ через параметризацию.

Решение

Из примера 4: $|\vec r'_\theta\times\vec r'_\varphi| = R^2\sin\theta$.

$$S = \int_0^{2\pi}d\varphi\int_0^\pi R^2\sin\theta\,d\theta = 2\pi\cdot R^2\cdot\bigl[-\cos\theta\bigr]_0^\pi = 2\pi R^2\cdot 2 = 4\pi R^2.$$
Пример 7 (часть сферы внутри цилиндра)

Найти площадь части сферы $x^2+y^2+z^2 = 4$, лежащей внутри цилиндра $x^2+y^2 = 1$ при $z \ge 0$.

Решение

Верхняя полусфера: $z = \sqrt{4 - x^2 - y^2}$. Проекция на $xOy$: круг $D$: $x^2+y^2\le 1$.

$z'_x = -\frac{x}{\sqrt{4-x^2-y^2}} = -\frac xz$, $z'_y = -\frac yz$.

$$1 + z_x'^2 + z_y'^2 = \frac{z^2 + x^2 + y^2}{z^2} = \frac{4}{z^2}\ \Rightarrow\ \sqrt{\ldots} = \frac{2}{\sqrt{4-x^2-y^2}}.$$

В полярных координатах ($dx\,dy = r\,dr\,d\varphi$):

$$S = \int_0^{2\pi}d\varphi\int_0^1\frac{2r\,dr}{\sqrt{4-r^2}} = 2\pi\cdot 2\Bigl[-\sqrt{4-r^2}\Bigr]_0^1 = 4\pi\,(2-\sqrt3)\approx 3.37.$$

Проверка формулой площади сферического сегмента $2\pi Rh$: высота сегмента $h = 2 - \sqrt{4-1} = 2-\sqrt3$, $2\pi\cdot2\cdot(2-\sqrt3)$ — совпало.

Пример 8 (часть параболоида)

Найти площадь части параболоида $z = x^2+y^2$, вырезанной цилиндром $x^2+y^2 = 2$.

Решение

$z'_x = 2x$, $z'_y = 2y$, $\sqrt{1 + 4x^2 + 4y^2} = \sqrt{1+4r^2}$. Проекция — круг $r\le\sqrt2$.

$$S = \int_0^{2\pi}d\varphi\int_0^{\sqrt2}\sqrt{1+4r^2}\,r\,dr = 2\pi\cdot\frac{1}{12}(1+4r^2)^{3/2}\Big|_0^{\sqrt2} = \frac{\pi}{6}(27 - 1) = \frac{13\pi}{3}.$$

Внутренний интеграл — подстановка $t = 1 + 4r^2$, $dt = 8r\,dr$: $\int\sqrt t\,\frac{dt}{8} = \frac{1}{12}t^{3/2}$.

Пример 9 (часть конуса)

Найти площадь части конуса $z = \sqrt{x^2+y^2}$, лежащей внутри цилиндра $x^2+y^2 = 2x$.

Решение

$z'_x = \frac{x}{\sqrt{x^2+y^2}}$, $z'_y = \frac{y}{\sqrt{x^2+y^2}}$, $z_x'^2 + z_y'^2 = \frac{x^2+y^2}{x^2+y^2} = 1$, поэтому $\sqrt{1+1} = \sqrt2$ — постоянный. У конуса с углом $45°$ всегда так.

Проекция: $x^2 + y^2 \le 2x \iff (x-1)^2 + y^2\le 1$ — круг радиуса 1, его площадь $\pi$.

$$S = \iint_D\sqrt2\,dx\,dy = \sqrt2\cdot\pi.$$

Двойной интеграл даже не пришлось считать: при постоянном корне площадь поверхности = корень × площадь проекции.

Пример 10 из примера КР2 в РПД

Треугольник $ABC$ образован плоскостью $3x - 2y - z + 6 = 0$, отсечённой координатными плоскостями. При помощи поверхностного интеграла найти его площадь. Проверить результат по формуле $S = \frac12\left|\overrightarrow{AB}\times\overrightarrow{AC}\right|$.

Решение

Вершины. Точки пересечения с осями: $y = z = 0 \Rightarrow 3x + 6 = 0$, $A(-2, 0, 0)$; $x = z = 0 \Rightarrow -2y + 6 = 0$, $B(0, 3, 0)$; $x = y = 0 \Rightarrow -z + 6 = 0$, $C(0, 0, 6)$.

Площадь как поверхностный интеграл. $S = \iint_{ABC} dS$. Выразим $z$: $z = 3x - 2y + 6$. $z'_x = 3$, $z'_y = -2$:

$$dS = \sqrt{1 + 9 + 4}\,dx\,dy = \sqrt{14}\,dx\,dy.$$

Проекция треугольника на $xOy$ — треугольник $D$ с вершинами $A'(-2,0)$, $B'(0,3)$, $O(0,0)$ (точка $C$ проектируется в начало координат). Он ограничен осями и прямой $3x - 2y + 6 = 0$, то есть $y = \frac32x + 3$.

$$S = \sqrt{14}\int_{-2}^0dx\int_0^{\frac32x+3}dy = \sqrt{14}\int_{-2}^0\Bigl(\frac32x + 3\Bigr)dx = \sqrt{14}\Bigl[\frac34x^2 + 3x\Bigr]_{-2}^0 = \sqrt{14}\,(0 - (3 - 6)) = 3\sqrt{14}.$$

(Интеграл от 1 — площадь $D$: $\frac12\cdot2\cdot3 = 3$, так что можно сразу $S = 3\sqrt{14}$.)

Проверка. $\overrightarrow{AB} = (2, 3, 0)$, $\overrightarrow{AC} = (2, 0, 6)$.

$$\overrightarrow{AB}\times\overrightarrow{AC} = \begin{vmatrix}\vec i & \vec j & \vec k\\ 2 & 3 & 0\\ 2 & 0 & 6\end{vmatrix} = \vec i\,(18 - 0) - \vec j\,(12 - 0) + \vec k\,(0 - 6) = (18, -12, -6).$$

$|(18,-12,-6)| = \sqrt{324 + 144 + 36} = \sqrt{504} = 6\sqrt{14}$, $S = \frac12\cdot6\sqrt{14} = 3\sqrt{14}\approx 11.2$. Совпало.

Заметь: $(18,-12,-6) = 6\cdot(3,-2,-1)$ — векторное произведение параллельно нормали плоскости, как и должно быть.

Пример 11 (параметрическая поверхность: геликоид)

Найти площадь куска винтовой поверхности $\vec r(u,v) = (u\cos v,\ u\sin v,\ v)$, $0\le u\le1$, $0\le v\le 2\pi$.

Решение

$\vec r'_u = (\cos v,\ \sin v,\ 0)$, $\vec r'_v = (-u\sin v,\ u\cos v,\ 1)$.

$$\vec r'_u\times\vec r'_v = (\sin v\cdot1 - 0,\ \ 0 - \cos v\cdot 1,\ \ u\cos^2v + u\sin^2v) = (\sin v,\ -\cos v,\ u),\qquad |\ldots| = \sqrt{1+u^2}.$$ $$S = \int_0^{2\pi}dv\int_0^1\sqrt{1+u^2}\,du = 2\pi\cdot\frac12\Bigl[u\sqrt{1+u^2} + \ln\bigl(u + \sqrt{1+u^2}\bigr)\Bigr]_0^1 = \pi\bigl(\sqrt2 + \ln(1+\sqrt2)\bigr)\approx 7.21.$$

Явно $z = f(x,y)$ эту поверхность не записать: над каждой точкой круга несколько витков. Поэтому параметризация здесь — единственный удобный путь.

3. Поверхностный интеграл 1 рода

вопрос экзамена №30ДЗ4, задача 1

Идея: масса оболочки

Тонкая оболочка (кусок жести, мыльная плёнка) имеет поверхностную плотность $\rho(x,y,z)$ — массу на единицу площади. Масса маленького кусочка площади $\Delta S$ около точки $M$ примерно $\rho(M)\Delta S$; складываем по всем кусочкам и измельчаем — получаем интеграл. Это прямой аналог криволинейного интеграла 1 рода (масса проволоки), только вместо $dl$ стоит $dS$.

Определение

Пусть на поверхности $S$ задана функция $f(x,y,z)$. Разобьём $S$ на кусочки $S_1,\ldots,S_n$ с площадями $\Delta S_k$, выберем в каждом точку $M_k$ и составим интегральную сумму $\sum_{k=1}^n f(M_k)\,\Delta S_k$. Если при стремлении наибольшего диаметра кусочков к нулю суммы имеют предел, не зависящий от разбиения и выбора точек, он называется поверхностным интегралом 1 рода:

$$\iint\limits_S f(x,y,z)\,dS = \lim_{\lambda\to0}\sum_{k=1}^n f(M_k)\,\Delta S_k.$$

Сторона поверхности в определении не участвует: $\Delta S_k > 0$ всегда.

Сведение к двойному интегралу

Если $S$: $z = z(x,y)$, $(x,y)\in D$, функция $z$ гладкая, $f$ непрерывна на $S$, то

$$\iint\limits_S f(x,y,z)\,dS = \iint\limits_D f\bigl(x,\ y,\ z(x,y)\bigr)\sqrt{1 + z_x'^2 + z_y'^2}\,dx\,dy.$$

Если $S$: $\vec r(u,v)$, $(u,v)\in\Omega$:

$$\iint\limits_S f\,dS = \iint\limits_\Omega f\bigl(\vec r(u,v)\bigr)\,|\vec r'_u\times\vec r'_v|\,du\,dv.$$

То есть: подставить уравнение поверхности в $f$, заменить $dS$ на $\sqrt{\ldots}\,dx\,dy$ и интегрировать по проекции.

Свойства — как у любого интеграла: линейность, аддитивность по кускам поверхности, $\iint_S 1\,dS = S$ (площадь), оценка $|\iint f\,dS| \le \max|f|\cdot S$, теорема о среднем.

Приложения

ВеличинаФормула ($\rho$ — поверхностная плотность)
Площадь$S = \iint_S dS$
Масса$M = \iint_S\rho\,dS$
Статические моменты относительно координатных плоскостей$M_{yz} = \iint_S x\rho\,dS$, $M_{xz} = \iint_S y\rho\,dS$, $M_{xy} = \iint_S z\rho\,dS$
Центр масс$x_c = \frac{M_{yz}}{M}$, $y_c = \frac{M_{xz}}{M}$, $z_c = \frac{M_{xy}}{M}$
Моменты инерции относительно осей$I_x = \iint_S(y^2+z^2)\rho\,dS$, $I_y = \iint_S(x^2+z^2)\rho\,dS$, $I_z = \iint_S(x^2+y^2)\rho\,dS$
Момент инерции относительно начала координат$I_O = \iint_S(x^2+y^2+z^2)\rho\,dS$

В моменте инерции стоит квадрат расстояния до оси: до оси $Oz$ расстояние $\sqrt{x^2+y^2}$. Индекс у статического момента — плоскость, от которой меряем: расстояние до $yOz$ равно $x$.

Алгоритм
  1. Нарисовать $S$ и выбрать плоскость, на которую она проектируется однозначно. Найти проекцию $D$.
  2. Выразить третью координату через две другие, найти $dS = \sqrt{1 + z_x'^2 + z_y'^2}\,dx\,dy$ (или по аналогии для других плоскостей).
  3. Подставить уравнение поверхности в подынтегральную функцию: на $S$ координата $z$ — не переменная, а $z(x,y)$.
  4. Посчитать двойной интеграл по $D$ (часто в полярных координатах). Поискать симметрию: нечётные по $x$ или $y$ слагаемые на симметричной области дают 0.
Пример 12 (масса пластинки)

Найти массу части плоскости $x + y + z = 1$ в первом октанте с плотностью $\rho = z$.

Решение

$z = 1 - x - y$, $z'_x = z'_y = -1$, $dS = \sqrt3\,dx\,dy$. Проекция — треугольник $D$: $x\ge0$, $y\ge0$, $x + y\le1$.

$$M = \iint_D(1 - x - y)\sqrt3\,dx\,dy = \sqrt3\int_0^1dx\int_0^{1-x}(1-x-y)\,dy = \sqrt3\int_0^1\frac{(1-x)^2}{2}dx = \sqrt3\cdot\frac16 = \frac{\sqrt3}{6}.$$

Внутренний интеграл: $\int_0^{a}(a - y)\,dy = \frac{a^2}{2}$ при $a = 1 - x$.

Пример 13

Вычислить $\iint_S z\,dS$, где $S$ — верхняя полусфера $x^2 + y^2 + z^2 = R^2$, $z\ge0$.

Решение

Как в примере 7: $dS = \frac{R}{z}\,dx\,dy$, где $z = \sqrt{R^2 - x^2 - y^2}$. Тогда $z\,dS = R\,dx\,dy$ — корень сократился.

$$\iint_S z\,dS = \iint_{x^2+y^2\le R^2}R\,dx\,dy = R\cdot\pi R^2 = \pi R^3.$$
Пример 14 (центр масс)

Найти центр масс однородной полусферы $x^2+y^2+z^2 = R^2$, $z\ge0$.

Решение

Однородная — $\rho = \text{const}$, она сократится. По симметрии $x_c = y_c = 0$ (полусфера симметрична относительно плоскостей $xOz$ и $yOz$, интегралы от $x$ и от $y$ нулевые).

$M = \rho\cdot 2\pi R^2$ (половина площади сферы), $M_{xy} = \rho\iint_S z\,dS = \rho\pi R^3$ (пример 13).

$$z_c = \frac{\rho\pi R^3}{2\pi\rho R^2} = \frac R2.$$

Центр масс полусферы-оболочки — ровно посередине высоты. Для сплошного полушара было бы $\frac{3R}{8}$: у оболочки больше массы «наверху».

Пример 15 (момент инерции)

Найти момент инерции однородной сферы $x^2+y^2+z^2 = R^2$ массы $M$ относительно оси $Oz$.

Решение

$I_z = \rho\iint_S(x^2+y^2)\,dS$. Удобный приём через симметрию: сфера одинаково выглядит относительно всех осей, значит $\iint x^2\,dS = \iint y^2\,dS = \iint z^2\,dS$. Сумма трёх равна $\iint(x^2+y^2+z^2)\,dS = \iint R^2\,dS = R^2\cdot4\pi R^2$ (на сфере $x^2+y^2+z^2 = R^2$). Каждый — $\frac{4\pi R^4}{3}$.

$$\iint_S(x^2+y^2)\,dS = \frac{8\pi R^4}{3},\qquad I_z = \rho\cdot\frac{8\pi R^4}{3} = \frac{M}{4\pi R^2}\cdot\frac{8\pi R^4}{3} = \frac23MR^2.$$

Напрямую через проекцию тоже можно: две полусферы, $dS = \frac{R\,dx\,dy}{\sqrt{R^2-r^2}}$, $\iint(x^2+y^2)\,dS = 2\int_0^{2\pi}d\varphi\int_0^R\frac{R\,r^3\,dr}{\sqrt{R^2-r^2}} = \frac{8\pi R^4}{3}$ (проверено).

Пример 16 (масса оболочки параболоида)

Найти массу части параболоида $z = x^2 + y^2$, $z\le1$, с плотностью $\rho = \sqrt{1 + 4z}$.

Решение

$dS = \sqrt{1 + 4x^2 + 4y^2}\,dx\,dy$, на поверхности $\rho = \sqrt{1 + 4(x^2 + y^2)}$ — тот же корень. Проекция — круг $r\le1$.

$$M = \iint_D(1 + 4r^2)\,r\,dr\,d\varphi = 2\pi\int_0^1(r + 4r^3)\,dr = 2\pi\Bigl(\frac12 + 1\Bigr) = 3\pi.$$
Пример 17 (боковая поверхность цилиндра)

Вычислить $\iint_S z\,dS$, где $S$ — боковая поверхность цилиндра $x^2 + y^2 = 1$, $0\le z\le2$.

Решение

Ловушка. На $xOy$ цилиндр проектируется в окружность — у неё нулевая площадь, формула с $\sqrt{1 + z_x'^2 + z_y'^2}$ не работает ($z$ не функция от $x$, $y$).

Способ 1: параметризация. $\vec r(\varphi, z) = (\cos\varphi, \sin\varphi, z)$, из примера 3 $|\vec r'_\varphi\times\vec r'_z| = 1$, то есть $dS = d\varphi\,dz$.

$$\iint_S z\,dS = \int_0^{2\pi}d\varphi\int_0^2 z\,dz = 2\pi\cdot2 = 4\pi.$$

Способ 2: проекция на $xOz$. Цилиндр состоит из двух половин $y = \pm\sqrt{1-x^2}$, каждая однозначно проектируется на прямоугольник $-1\le x\le1$, $0\le z\le2$. $y'_x = \mp\frac{x}{\sqrt{1-x^2}}$, $y'_z = 0$, $\sqrt{1 + y_x'^2} = \frac{1}{\sqrt{1-x^2}}$.

$$\iint_S z\,dS = 2\int_{-1}^1\frac{dx}{\sqrt{1-x^2}}\int_0^2z\,dz = 2\cdot\pi\cdot2 = 4\pi.$$

Совпало. Параметризация короче, и делить на половинки не нужно.

Пример 18 из ДЗ4 в РПД

Вычислить поверхностный интеграл первого рода $\iint_S(xy + yz + zx)\,dS$, где $S$ — часть конической поверхности $z = \sqrt{x^2+y^2}$, расположенная внутри цилиндра $x^2+y^2 = 2x$.

Решение

Поверхность и проекция — как в примере 9 (картинка там же): $dS = \sqrt2\,dx\,dy$, $D$: $(x-1)^2 + y^2\le1$.

Полярные координаты. Граница $x^2+y^2 = 2x$ в полярных: $r^2 = 2r\cos\varphi$, то есть $r = 2\cos\varphi$; круг лежит справа от оси $Oy$, поэтому $\varphi\in\left[-\frac\pi2, \frac\pi2\right]$, $0\le r\le2\cos\varphi$.

Подынтегральная функция на конусе ($z = r$): $xy + z(x + y) = r^2\cos\varphi\sin\varphi + r\cdot r(\cos\varphi + \sin\varphi) = r^2(\sin\varphi\cos\varphi + \cos\varphi + \sin\varphi)$.

$$I = \sqrt2\int_{-\pi/2}^{\pi/2}(\sin\varphi\cos\varphi + \cos\varphi + \sin\varphi)\,d\varphi\int_0^{2\cos\varphi}r^2\cdot r\,dr = \sqrt2\int_{-\pi/2}^{\pi/2}(\sin\varphi\cos\varphi + \cos\varphi + \sin\varphi)\cdot4\cos^4\varphi\,d\varphi.$$

Симметрия. Отрезок симметричен относительно 0, а слагаемые $\sin\varphi\cos^5\varphi$ и $\sin\varphi\cos^4\varphi$ нечётные — их интегралы равны 0. Остаётся

$$I = 4\sqrt2\int_{-\pi/2}^{\pi/2}\cos^5\varphi\,d\varphi = 8\sqrt2\int_0^{\pi/2}\cos^5\varphi\,d\varphi = 8\sqrt2\cdot\frac{4\cdot2}{5\cdot3} = \frac{64\sqrt2}{15}\approx 6.03.$$

$\int_0^{\pi/2}\cos^5\varphi\,d\varphi$ можно и подстановкой $t = \sin\varphi$: $\int_0^1(1-t^2)^2dt = 1 - \frac23 + \frac15 = \frac{8}{15}$.

Ответ проверен в sympy.

4. Поверхностный интеграл 2 рода (поток)

вопрос экзамена №40вопрос экзамена №41

Идея: сколько жидкости протекает через поверхность

Жидкость течёт со скоростью $\vec F(x,y,z)$. Сколько её проходит за единицу времени через маленькую площадку $\Delta S$ с единичной нормалью $\vec n$? За время $1$ частицы, прошедшие через площадку, заполняют наклонный цилиндр с основанием $\Delta S$ и образующей $\vec F$. Его объём = площадь основания × высота, а высота — проекция $\vec F$ на нормаль: $(\vec F,\vec n)$. Итого через площадку проходит $(\vec F,\vec n)\,\Delta S$.

Знак зависит от того, какую нормаль выбрали, то есть от стороны поверхности. В этом главное отличие от интеграла 1 рода.

Определение (поток)

Пусть $S$ — двусторонняя гладкая поверхность с выбранной стороной (полем единичных нормалей $\vec n = (\cos\alpha, \cos\beta, \cos\gamma)$), $\vec F = (P, Q, R)$ — векторное поле. Потоком $\vec F$ через $S$ в сторону $\vec n$ называется

$$\Pi = \iint\limits_S(\vec F,\vec n)\,dS = \iint\limits_S(P\cos\alpha + Q\cos\beta + R\cos\gamma)\,dS.$$

Это и есть поверхностный интеграл 2 рода. Его записывают ещё так:

$$\iint\limits_S P\,dy\,dz + Q\,dz\,dx + R\,dx\,dy,\qquad dy\,dz = \cos\alpha\,dS,\ \ dz\,dx = \cos\beta\,dS,\ \ dx\,dy = \cos\gamma\,dS.$$

Здесь $dx\,dy$ — не обычный элемент площади, а площадь проекции кусочка поверхности на $xOy$ со знаком: «$+$», если нормаль смотрит вверх ($\cos\gamma > 0$), «$-$», если вниз.

Связь 1 и 2 рода, свойства $$\iint\limits_S P\,dy\,dz + Q\,dz\,dx + R\,dx\,dy = \iint\limits_S(P\cos\alpha + Q\cos\beta + R\cos\gamma)\,dS,$$

где $\cos\alpha, \cos\beta, \cos\gamma$ — направляющие косинусы нормали выбранной стороны. Отсюда:

Как выбрать знак нормали

Способ 1: через одну проекцию

Поверхность $z = f(x,y)$, $(x,y)\in D$. Нормаль верхней стороны $\vec n = \frac{(-f'_x, -f'_y, 1)}{\sqrt{1 + f_x'^2 + f_y'^2}}$, а $dS = \sqrt{1 + f_x'^2 + f_y'^2}\,dx\,dy$. Корни сокращаются, и $\vec n\,dS = (-f'_x, -f'_y, 1)\,dx\,dy$:

Поток через график $$\iint\limits_S(\vec F,\vec n)\,dS = \pm\iint\limits_D\Bigl(-P\,f'_x - Q\,f'_y + R\Bigr)\Big|_{z = f(x,y)}dx\,dy,$$

«$+$» для верхней стороны, «$-$» для нижней. Для параметрической поверхности аналогично:

$$\iint\limits_S(\vec F,\vec n)\,dS = \pm\iint\limits_\Omega\bigl(\vec F(\vec r(u,v)),\ \vec r'_u\times\vec r'_v\bigr)\,du\,dv,$$

«$+$», если $\vec r'_u\times\vec r'_v$ смотрит в выбранную сторону.

Алгоритм (одна проекция)
  1. Нарисовать $S$, понять, какая сторона нужна.
  2. Выбрать плоскость, на которую $S$ проектируется однозначно, найти проекцию $D$.
  3. Записать неединичную нормаль $\vec N$ (например $(-f'_x, -f'_y, 1)$) и проверить, что она смотрит в нужную сторону; если нет — сменить знак.
  4. Посчитать $(\vec F,\vec N)$ и подставить уравнение поверхности.
  5. Посчитать двойной интеграл по $D$. Никаких корней: они сократились.

Способ 2: через три проекции

Интеграл $\iint P\,dy\,dz + Q\,dz\,dx + R\,dx\,dy$ считают как сумму трёх, каждый — своей проекцией:

$$\iint\limits_S R\,dx\,dy = \pm\iint\limits_{D_{xy}}R\bigl(x, y, z(x,y)\bigr)\,dx\,dy,\qquad(+\ \text{при}\ \cos\gamma > 0),$$ $$\iint\limits_S P\,dy\,dz = \pm\iint\limits_{D_{yz}}P\bigl(x(y,z), y, z\bigr)\,dy\,dz,\qquad(+\ \text{при}\ \cos\alpha > 0),$$ $$\iint\limits_S Q\,dz\,dx = \pm\iint\limits_{D_{xz}}Q\bigl(x, y(x,z), z\bigr)\,dz\,dx,\qquad(+\ \text{при}\ \cos\beta > 0).$$

Каждое слагаемое считают со своим знаком: у плоскости $x + 2y - z = 0$ нормаль $(1, 2, -1)$, и если выбрать её, то два слагаемых берутся с «$+$», а третье с «$-$». Если поверхность проектируется на плоскость неоднозначно, её режут на куски (пример 22). Если поверхность перпендикулярна плоскости, соответствующее слагаемое равно 0. Способ длиннее, зато иногда он единственный, если интеграл уже записан в форме с $dy\,dz$, $dz\,dx$, $dx\,dy$.

Пример 19 (поток через часть плоскости)

Найти поток поля $\vec F = (x, y, z)$ через часть плоскости $x + y + z = 1$ в первом октанте, в сторону верхней нормали.

Решение

$z = 1 - x - y$, $\vec N = (-f'_x, -f'_y, 1) = (1, 1, 1)$ — третья координата $> 0$, смотрит вверх, подходит.

$(\vec F,\vec N) = x + y + z = 1$ на плоскости. Проекция $D$ — треугольник площадью $\frac12$.

$$\Pi = \iint_D 1\,dx\,dy = \frac12.$$

Через единичную нормаль то же самое: $\vec n = \frac{(1,1,1)}{\sqrt3}$, $(\vec F,\vec n) = \frac{1}{\sqrt3}$, $dS = \sqrt3\,dx\,dy$, произведение $= dx\,dy$.

Пример 20 (нижняя сторона)

Вычислить $\iint_S z\,dx\,dy$ по нижней стороне поверхности $z = x^2+y^2$ над треугольником $D$: $x\ge0$, $y\ge0$, $x + y\le1$.

Решение

Это интеграл вида $\iint R\,dx\,dy$ с $R = z$. Нижняя сторона: $\cos\gamma < 0$, знак «$-$».

$$\iint_S z\,dx\,dy = -\iint_D(x^2+y^2)\,dx\,dy = -\int_0^1dx\int_0^{1-x}(x^2+y^2)\,dy = -\int_0^1\Bigl(x^2(1-x) + \frac{(1-x)^3}{3}\Bigr)dx = -\Bigl(\frac1{12} + \frac1{12}\Bigr) = -\frac16.$$
Пример 21 (три проекции против одной)

Вычислить $\iint_S x^2\,dy\,dz + y^2\,dz\,dx + z^2\,dx\,dy$, где $S$ — верхняя сторона части плоскости $x + y + z = 1$ в первом октанте.

Решение

Нормаль верхней стороны $\frac{(1,1,1)}{\sqrt3}$: все три косинуса положительны, все знаки «$+$».

Три проекции. $\iint_S x^2\,dy\,dz$: на $yOz$ проекция — треугольник $y, z\ge0$, $y + z\le1$, $x = 1 - y - z$:

$$\iint_{D_{yz}}(1 - y - z)^2\,dy\,dz = \int_0^1dy\int_0^{1-y}(1-y-z)^2dz = \int_0^1\frac{(1-y)^3}{3}dy = \frac{1}{12}.$$

Два других слагаемых по симметрии такие же. Ответ: $\frac{3}{12} = \frac14$.

Одна проекция. $\vec N = (1,1,1)$, $(\vec F,\vec N) = x^2 + y^2 + z^2$ с $z = 1 - x - y$:

$$\iint_D\bigl(x^2 + y^2 + (1-x-y)^2\bigr)dx\,dy = \frac1{12} + \frac1{12} + \frac1{12} = \frac14.$$

Совпало.

Пример 22 (поток через боковую поверхность цилиндра)

Найти поток поля $\vec F = (x, y, z)$ через боковую поверхность цилиндра $x^2 + y^2 = R^2$, $0\le z\le H$, наружу.

Решение

Через единичную нормаль. Внешняя нормаль цилиндра $\vec n = \frac{(x, y, 0)}{R}$ (пример 3). $(\vec F,\vec n) = \frac{x^2+y^2}{R} = \frac{R^2}{R} = R$ — постоянна на поверхности.

$$\Pi = \iint_S R\,dS = R\cdot 2\pi RH = 2\pi R^2H.$$

Через три проекции (для тренировки). $\Pi = \iint x\,dy\,dz + y\,dz\,dx + z\,dx\,dy$.

  • $\iint z\,dx\,dy = 0$: цилиндр вертикальный, $\cos\gamma = 0$.
  • $\iint x\,dy\,dz$: на $yOz$ цилиндр проектируется неоднозначно (каждая точка прямоугольника $|y|\le R$, $0\le z\le H$ — две точки цилиндра), режем на половины. Передняя $x = \sqrt{R^2-y^2}$: внешняя нормаль смотрит в сторону $+x$, знак «$+$», вклад $\iint\sqrt{R^2-y^2}\,dy\,dz$. Задняя $x = -\sqrt{R^2-y^2}$: нормаль в сторону $-x$, знак «$-$», вклад $-\iint(-\sqrt{R^2-y^2})\,dy\,dz$ — то же самое. Итого $2H\int_{-R}^R\sqrt{R^2-y^2}\,dy = 2H\cdot\frac{\pi R^2}{2} = \pi R^2H$.
  • $\iint y\,dz\,dx = \pi R^2H$ так же.

Сумма $2\pi R^2H$. Совпало, но первый способ в разы короче.

Пример 23 (поток через сферу)

Найти поток поля $\vec F = (x, y, z)$ через сферу $x^2+y^2+z^2 = R^2$ наружу.

Решение

$\vec n = \frac{(x,y,z)}{R}$, $(\vec F,\vec n) = \frac{x^2+y^2+z^2}{R} = R$ на сфере.

$$\Pi = R\cdot4\pi R^2 = 4\pi R^3.$$

Сравнение с Остроградским-Гауссом (глава 10): поток наружу через замкнутую поверхность равен $\iiint_V\div\vec F\,dV$, где $\div\vec F = P'_x + Q'_y + R'_z$. Здесь $\div\vec F = 1 + 1 + 1 = 3$, объём шара $\frac43\pi R^3$, и $3\cdot\frac43\pi R^3 = 4\pi R^3$. Совпало.

Через три проекции было бы так: $\iint z\,dx\,dy$ по внешней стороне = (верхняя полусфера, знак «$+$») + (нижняя, знак «$-$», $z = -\sqrt{\ldots}$) $= 2\iint_D\sqrt{R^2-x^2-y^2}\,dx\,dy = \frac43\pi R^3$, и три таких слагаемых дают $4\pi R^3$.

Пример 24 (поток через параболоид и проверка Гауссом)

Найти поток поля $\vec F = (x, y, z)$ через верхнюю сторону параболоида $z = 1 - x^2 - y^2$, $z\ge0$.

Решение

$f'_x = -2x$, $f'_y = -2y$, $\vec N = (-f'_x, -f'_y, 1) = (2x, 2y, 1)$ — смотрит вверх, подходит.

$$(\vec F,\vec N) = 2x^2 + 2y^2 + z = 2x^2 + 2y^2 + 1 - x^2 - y^2 = 1 + x^2 + y^2.$$

Проекция — круг $r\le1$ ($z = 0$ при $x^2+y^2 = 1$).

$$\Pi = \int_0^{2\pi}d\varphi\int_0^1(1 + r^2)\,r\,dr = 2\pi\Bigl(\frac12 + \frac14\Bigr) = \frac{3\pi}{2}.$$

Проверка через Остроградского-Гаусса (глава 10). Замкнём поверхность кругом $z = 0$, $x^2+y^2\le1$. Внешняя нормаль у круга $(0,0,-1)$, $(\vec F,\vec n) = -z = 0$ на нём, поток через дно 0. По теореме весь поток наружу $= \iiint\div\vec F\,dV = 3V$, а $V = \iint_D(1 - r^2)\,dx\,dy = 2\pi\cdot\frac14 = \frac\pi2$. Получаем $\frac{3\pi}{2} + 0 = \frac{3\pi}{2}$. Совпало.

Поток $\vec F$ через верхнюю сторону поверхности равен 5. Чему равен поток через нижнюю сторону и интеграл $\iint_S|\vec F|\,dS$ — меняется ли он?
Смена стороны меняет знаки всех направляющих косинусов, поэтому поток меняет знак. В интеграле 1 рода сторона не участвует, $dS > 0$.
Чему равен $\iint_S(x^2 + e^y)\,dx\,dy$ по внешней стороне боковой поверхности цилиндра $x^2+y^2 = 4$, $0\le z\le 3$?
Образующие цилиндра параллельны $Oz$, $\cos\gamma = 0$, проекция на $xOy$ — окружность нулевой площади. Любой интеграл вида $\iint R\,dx\,dy$ по такой поверхности равен 0.

5. Как это выглядит в РГР

РГР2, задание 5a

В типовом расчёте (Попов, «Типовой расчёт, 3 семестр») поверхностные интегралы есть в задании 5. Его пункты:

$K$, $L$, $M$ — точки пересечения $\sigma$ с осями $Ox$, $Oy$, $Oz$. Требования к оформлению из РПД: отдельная тетрадь, полное условие перед каждым заданием, чертежи на миллиметровке.

Опечатка в образце решения

В образце в начале пособия про направляющие косинусы написано, что перед радикалом берут знак «−» и в одном случае, и в другом. Правильно: знак «+», если нормаль образует острый угол с $Oz$ ($\cos\gamma > 0$), и «−» в противном случае. Из самого образца это видно: там $\cos\gamma > 0$ и выбран «+».

Пример 25 РГР, задание 5a, вариант 1

Дано векторное поле $\vec a = (y - x)\,\vec i + (3y - z)\,\vec k$ и плоскость $\sigma$: $2x + y + z = 4$. Найти поток $Q$ поля $\vec a$ через часть $S$ плоскости $\sigma$, вырезанную координатными плоскостями, в сторону нормали $\vec n$, направленной от начала координат $O(0,0,0)$.

Решение

1. Чертёж. Точки пересечения с осями: $K(2,0,0)$ ($y = z = 0$, $2x = 4$), $L(0,4,0)$, $M(0,0,4)$. Часть $S$ — треугольник $KLM$ в первом октанте (рисунок ниже).

2. Формула. $Q = \iint_S(\vec a,\vec n)\,dS$, где $\vec n = (\cos\alpha, \cos\beta, \cos\gamma)$. Компоненты поля: $a_x = y - x$, $a_y = 0$, $a_z = 3y - z$.

3. Нормаль. Из уравнения плоскости $z = 4 - 2x - y$, $p = z'_x = -2$, $q = z'_y = -1$, $\sqrt{1 + p^2 + q^2} = \sqrt{1 + 4 + 1} = \sqrt6$.

$$\cos\alpha = \frac{-p}{\pm\sqrt6},\quad\cos\beta = \frac{-q}{\pm\sqrt6},\quad\cos\gamma = \frac{1}{\pm\sqrt6}.$$

Треугольник лежит над плоскостью $xOy$, и нормаль «от начала координат» смотрит вверх: $\cos\gamma > 0$, берём «+». (Другая проверка: в начале координат $2x + y + z = 0 < 4$, значит от начала смотрит вектор $(2,1,1)$ — в сторону роста $2x + y + z$.)

$$\vec n = \Bigl(\frac{2}{\sqrt6},\ \frac1{\sqrt6},\ \frac1{\sqrt6}\Bigr).$$

4. Скалярное произведение на плоскости.

$$(\vec a,\vec n) = \frac{2(y - x) + 0\cdot1 + (3y - z)\cdot1}{\sqrt6}.$$

Подставляем $z = 4 - 2x - y$: $3y - z = 3y - 4 + 2x + y = 2x + 4y - 4$, и числитель $2y - 2x + 2x + 4y - 4 = 6y - 4$:

$$(\vec a,\vec n) = \frac{6y - 4}{\sqrt6}.$$

5. Сведение к двойному. $dS = \sqrt6\,dx\,dy$, корни сокращаются. Проекция $D_{xy}$ — треугольник $OKL$: $0\le x\le2$, $0\le y\le4 - 2x$.

$$Q = \iint_{D_{xy}}(6y - 4)\,dx\,dy = \int_0^2dx\int_0^{4-2x}(6y - 4)\,dy.$$

6. Вычисление. Внутренний: $\bigl[3y^2 - 4y\bigr]_0^{4-2x} = (4 - 2x)\bigl(3(4 - 2x) - 4\bigr) = (4-2x)(8-6x) = 12x^2 - 40x + 32$.

$$Q = \int_0^2(12x^2 - 40x + 32)\,dx = \Bigl[4x^3 - 20x^2 + 32x\Bigr]_0^2 = 32 - 80 + 64 = 16.$$

Ответ: $Q = 16$.

Проверка через три проекции. $Q = \iint_S(y - x)\,dy\,dz + 0\,dz\,dx + (3y - z)\,dx\,dy$, все косинусы положительны.

  • На $yOz$: $x = \frac{4 - y - z}{2}$, $D_{yz}$: $y, z\ge0$, $y + z\le4$ (площадь 8). $\iint_{D_{yz}}\Bigl(y - \frac{4-y-z}{2}\Bigr)dy\,dz = \frac12\iint_{D_{yz}}(3y + z - 4)\,dy\,dz = \frac12\Bigl(3\cdot\frac{32}{3} + \frac{32}{3} - 32\Bigr) = \frac{16}{3}$ (здесь $\iint_{D_{yz}} y\,dy\,dz = \iint_{D_{yz}} z\,dy\,dz = \frac{4^3}{6} = \frac{32}{3}$).
  • На $xOy$: $\iint_{D_{xy}}(2x + 4y - 4)\,dx\,dy = 2\cdot\frac83 + 4\cdot\frac{16}{3} - 4\cdot4 = \frac{32}{3}$ (здесь $\iint_{D_{xy}}x = \frac83$, $\iint_{D_{xy}}y = \frac{16}{3}$, площадь 4).

$\frac{16}{3} + \frac{32}{3} = 16$. Совпало.

Связь с пунктом b) (сам пункт — в главе 10). Поток через три грани тетраэдра на координатных плоскостях наружу: грань $x = 0$ (нормаль $-\vec i$): $-\iint_{D_{yz}}y = -\frac{32}{3}$; грань $y = 0$: $a_y = 0$, поток 0; грань $z = 0$ (нормаль $-\vec k$): $-\iint_{D_{xy}}3y = -16$. Полный поток наружу: $16 - \frac{32}{3} + 0 - 16 = -\frac{32}{3}$. По Остроградскому-Гауссу: $\div\vec a = -1 + 0 - 1 = -2$, объём тетраэдра $\frac16\cdot2\cdot4\cdot4 = \frac{16}{3}$, $-2\cdot\frac{16}{3} = -\frac{32}{3}$. Всё сходится.

Шаблон оформления задания 5a
  1. Найти $K$, $L$, $M$, сделать чертёж треугольника и проекции.
  2. Выразить $z$ из уравнения плоскости, найти $p$, $q$ и $\sqrt{1 + p^2 + q^2}$.
  3. Выбрать знак перед корнем по условию «от начала координат» и выписать $\cos\alpha, \cos\beta, \cos\gamma$.
  4. Посчитать $(\vec a,\vec n)$, подставить $z$.
  5. $dS = \sqrt{1 + p^2 + q^2}\,dx\,dy$, расставить пределы по треугольнику $OKL$, посчитать.

Если плоскость пересекает $Oz$ ниже нуля ($M$ под плоскостью $xOy$), нормаль от начала смотрит вниз, знак «−». Если $\sigma$ параллельна $Oz$, на $xOy$ проектировать нельзя — проектируют на другую плоскость.

6. Тренировка

Сначала реши сам, потом открывай решение. Все ответы проверены в sympy.

Задача 1

Найти площадь части плоскости $\frac x2 + \frac y3 + \frac z4 = 1$ в первом октанте двумя способами: поверхностным интегралом и через векторное произведение.

Решение

$z = 4 - 2x - \frac43y$, $\sqrt{1 + 4 + \frac{16}{9}} = \frac{\sqrt{61}}{3}$. Проекция — треугольник с катетами 2 и 3, площадь 3. $S = 3\cdot\frac{\sqrt{61}}{3} = \sqrt{61}$.

Вершины $A(2,0,0)$, $B(0,3,0)$, $C(0,0,4)$; $\overrightarrow{AB} = (-2,3,0)$, $\overrightarrow{AC} = (-2,0,4)$, $\overrightarrow{AB}\times\overrightarrow{AC} = (12, 8, 6)$, $|\ldots| = \sqrt{244} = 2\sqrt{61}$, $S = \sqrt{61}$.

Задача 2

Найти площадь части поверхности $z = xy$, лежащей внутри цилиндра $x^2+y^2 = 1$.

Решение

$z'_x = y$, $z'_y = x$, $\sqrt{1 + x^2 + y^2} = \sqrt{1 + r^2}$.

$$S = 2\pi\int_0^1\sqrt{1+r^2}\,r\,dr = 2\pi\cdot\frac13(1+r^2)^{3/2}\Big|_0^1 = \frac{2\pi}{3}\bigl(2\sqrt2 - 1\bigr).$$
Задача 3

Найти площадь части цилиндра $x^2 + y^2 = R^2$, $x\ge0$, между плоскостями $z = 0$ и $z = x$.

Решение

Параметризация $(R\cos\varphi, R\sin\varphi, z)$, $dS = R\,d\varphi\,dz$. $x\ge0$: $\varphi\in\left[-\frac\pi2,\frac\pi2\right]$; $0\le z\le x = R\cos\varphi$.

$$S = \int_{-\pi/2}^{\pi/2}d\varphi\int_0^{R\cos\varphi}R\,dz = R^2\int_{-\pi/2}^{\pi/2}\cos\varphi\,d\varphi = 2R^2.$$

На $xOy$ проектировать нельзя (цилиндр вертикален), только параметризация или проекция на $xOz$/$yOz$.

Задача 4

Вычислить $\iint_S(x + y + z)\,dS$, где $S$ — часть плоскости $x + y + z = a$ ($a > 0$) в первом октанте.

Решение

На плоскости $x + y + z = a$ — подынтегральная функция постоянна. $dS = \sqrt3\,dx\,dy$, проекция — треугольник площадью $\frac{a^2}{2}$.

$$\iint_S a\,dS = a\cdot\sqrt3\cdot\frac{a^2}{2} = \frac{\sqrt3\,a^3}{2}.$$
Задача 5

Вычислить $\iint_S\frac{dS}{(1 + x + y)^2}$, где $S$ — часть плоскости $x + y + z = 1$ в первом октанте.

Решение
$$\sqrt3\int_0^1dx\int_0^{1-x}\frac{dy}{(1+x+y)^2} = \sqrt3\int_0^1\Bigl(\frac{1}{1+x} - \frac12\Bigr)dx = \sqrt3\Bigl(\ln2 - \frac12\Bigr).$$

Внутренний: $\Bigl[-\frac{1}{1+x+y}\Bigr]_0^{1-x} = -\frac12 + \frac1{1+x}$.

Задача 6

Найти массу полусферы $z = \sqrt{R^2 - x^2 - y^2}$ с плотностью $\rho = x^2 + y^2$.

Решение

По примеру 15 по всей сфере $\iint(x^2+y^2)\,dS = \frac{8\pi R^4}{3}$, полусфера — половина: $M = \frac{4\pi R^4}{3}$.

Напрямую: $M = \int_0^{2\pi}d\varphi\int_0^R\frac{r^2\cdot R\,r\,dr}{\sqrt{R^2-r^2}}$. Подстановка $t = R^2 - r^2$: $\int_0^R\frac{r^3dr}{\sqrt{R^2-r^2}} = \frac12\int_0^{R^2}\frac{R^2 - t}{\sqrt t}dt = \frac12\Bigl(2R^3 - \frac23R^3\Bigr) = \frac23R^3$. $M = 2\pi R\cdot\frac23R^3 = \frac{4\pi R^4}{3}$.

Задача 7

Найти поток поля $\vec F = (0, 0, z)$ через верхнюю сторону полусферы $z = \sqrt{R^2 - x^2 - y^2}$.

Решение

Это $\iint_S z\,dx\,dy$ по верхней стороне: знак «+», $z$ заменяем на $\sqrt{R^2 - x^2 - y^2}$.

$$\Pi = \iint_{x^2+y^2\le R^2}\sqrt{R^2 - x^2 - y^2}\,dx\,dy = 2\pi\int_0^R\sqrt{R^2 - r^2}\,r\,dr = 2\pi\cdot\frac{R^3}{3} = \frac{2\pi R^3}{3}.$$

Это объём полушара: поле $(0,0,z)$ «считает» объём под поверхностью.

Задача 8

Найти поток поля $\vec F = (x, y, z)$ наружу через полную поверхность цилиндра $x^2+y^2\le1$, $0\le z\le1$ (боковая поверхность и две крышки). Проверить по Остроградскому-Гауссу.

Решение
  • Боковая: по примеру 22 при $R = H = 1$ поток $2\pi$.
  • Верхняя крышка $z = 1$, нормаль $\vec k$: $(\vec F,\vec n) = z = 1$, поток $= 1\cdot\pi = \pi$.
  • Нижняя крышка $z = 0$, нормаль $-\vec k$: $(\vec F,\vec n) = -z = 0$, поток 0.

Итого $3\pi$. По Гауссу: $\div\vec F = 3$, объём $\pi$, $3\pi$. Совпало.

Задача 9

Вычислить $\iint_S z^2\,dx\,dy$ по внешней стороне сферы $x^2 + y^2 + z^2 = R^2$.

Решение

Сфера проектируется на $xOy$ неоднозначно, режем на полусферы. Верхняя ($\cos\gamma > 0$, «+»): $\iint_D(R^2 - x^2 - y^2)\,dx\,dy$. Нижняя ($\cos\gamma < 0$, «−»): $-\iint_D(R^2 - x^2 - y^2)\,dx\,dy$ ($z^2$ на ней такой же). Сумма 0.

По Гауссу то же: поле $(0, 0, z^2)$, $\div = 2z$, интеграл от $2z$ по шару 0 по симметрии.

Задача 10

Найти поток поля $\vec F = (y, -x, z)$ через часть конуса $z = \sqrt{x^2+y^2}$, $z\le1$, в сторону нормали с $\cos\gamma > 0$.

Решение

$f'_x = \frac xr$, $f'_y = \frac yr$ ($r = \sqrt{x^2+y^2}$), $\vec N = \bigl(-\frac xr, -\frac yr, 1\bigr)$ — третья координата 1, подходит.

$$(\vec F,\vec N) = -\frac{xy}{r} + \frac{xy}{r} + z = z = r.$$ $$\Pi = \int_0^{2\pi}d\varphi\int_0^1r\cdot r\,dr = \frac{2\pi}{3}.$$

Вращательная часть поля $(y, -x, 0)$ касается конуса и ничего не даёт в поток.

7. Типичные ошибки

Что ломается чаще всего

Мини-шпаргалка главы

Нормаль:

  • $z = f(x,y)$: $\vec N = (-f'_x, -f'_y, 1)$ (вверх);
  • $F = 0$: $\vec N = \grad F$;
  • $\vec r(u,v)$: $\vec N = \vec r'_u\times\vec r'_v$;
  • сфера: $\frac{(x,y,z)}{R}$, цилиндр: $\frac{(x,y,0)}{R}$ (наружу).

Элемент площади:

$$dS = \sqrt{1 + f_x'^2 + f_y'^2}\,dx\,dy = \frac{dx\,dy}{|\cos\gamma|},$$ $$dS = |\vec r'_u\times\vec r'_v|\,du\,dv.$$

Сфера: $dS = \frac{R}{|z|}dx\,dy = R^2\sin\theta\,d\theta\,d\varphi$; конус $z = \sqrt{x^2+y^2}$: $\sqrt2\,dx\,dy$; цилиндр: $R\,d\varphi\,dz$.

1 род (площадь, масса, центр масс, моменты):

$$\iint_S f\,dS = \iint_D f\bigl(x,y,z(x,y)\bigr)\sqrt{1 + z_x'^2 + z_y'^2}\,dx\,dy.$$

$M = \iint\rho\,dS$, $z_c = \frac{1}{M}\iint z\rho\,dS$, $I_z = \iint(x^2+y^2)\rho\,dS$. От стороны не зависит.

2 род (поток):

$$\iint_S P\,dy\,dz + Q\,dz\,dx + R\,dx\,dy = \iint_S(\vec F,\vec n)\,dS.$$

Через одну проекцию: $\pm\iint_D(-Pf'_x - Qf'_y + R)\,dx\,dy$, «+» для верхней стороны.

Через три: $\iint R\,dx\,dy = \pm\iint_{D_{xy}}R(x,y,z(x,y))\,dx\,dy$, знак = знак $\cos\gamma$; аналогично для $P$ ($\cos\alpha$) и $Q$ ($\cos\beta$).

Свойства потока: смена стороны — смена знака; вертикальный цилиндр — $\iint R\,dx\,dy = 0$; замкнутая поверхность — можно проверить Остроградским-Гауссом $\oiint_S(\vec F,\vec n)\,dS = \iiint_V\div\vec F\,dV$ (глава 10).

Треугольник: $S_{ABC} = \frac12|\overrightarrow{AB}\times\overrightarrow{AC}|$. Лист Мёбиуса — односторонний, поток через него не определён. Пример Шварца: площади вписанных многогранников могут не стремиться к площади поверхности.