Двойной интеграл является мощным инструментом для определения объёмов тел, ограниченных поверхностями. Введение объясняет связь между двойным интегралом и объёмом, формулирует задачи статьи: демонстрацию методов вычисления объёмов и их практическое применение для различных типов областей.
Объём тела, ограниченного сверху и снизу поверхностями z = f _1( x , y ), z = f _2( x , y ), (x,y)∈D, с боков - цилиндрической поверхностью с образующими, параллельными оси Oz , равен v=∬_D[f_1(x,y)−f_2(x,y)]dxdy; эта формула очевидно следует из геометрического смысла двойного интеграла.
Основной вопрос, который надо решить - на какую координатную плоскость проектировать тело, чтобы выкладки были наиболее простыми.
Пример 1
Найти объём тела
V:y=0,;z=0,x+y+z=4,;2x+z=4.
Решение:
Тело изображено на рисунке. Перебором возможностей убеждаемся, что проще всего описать это тело, если отправляться от его проекции на ось Oxz :
V:[(x,z)∈D,0⩽y⩽4−x−z.
Область D - треугольник, ограниченный прямыми x = 0, z = 0, 2 x + z = 4, поэтому
Найти объём области, ограниченной поверхностями x 2+ y 2+ z 2= R 2,
( x 2+ y 2)3= R 2( x 2+ y 2).
Решение:
Первая поверхность - сфера, вторая - цилиндрическая - с образующими, параллельными оси Oz { в уравнении нет z в явной форме). Построить в плоскости Oxy кривую шестого порядка, заданную уравнением ( x 2+ y 2)3= R 2( x 2+ y 2), в декартовой системе координат невозможно, можно только сказать, что она симметрична относительно осей { чётные степени } и точка О(0,0) принадлежит этой кривой. Пробуем перейти к полярным координатам. r6=R2r4(cos4φ+sin4φ);r2=R2((cos2φ+sin2φ)2−2cos2φsin2φ)=R2(1−2sin22φ)=
=R2(1−41−cos4φ)=R243+cos4φ;r=R23+cos4φ. Эту кривую построить уже можно. r(φ) максимально, когда cos4φ=1;(φ=0,42π=2π,44π=π,46π=23π), минимально, когда
cos4φ=−1;(φ=4π,43π,45π,47π), и гладко меняется между этими пределами { точка О(0,0) не принадлежит этой кривой, где мы её потеряли? } .
Пользуясь симметрией, получаем V=16∬_DR2−x2−y2dxdy=16∬_DR2−r2rdrdφ==16∫_04πdφ∫_0R23+cos4φR2−r2rdr==−8∫_04πdφ∫_0R23+cos4φR2−r2d(R2−r2)=−832∫_04π(R2−r2)23_0R23+cos4φdφ=−316R3∫_04π[(2sin22φ)23−1]dφ= и т.д.
Пример 3
Найти объем тела в первом октанте, ограниченного плоскостями (y = 0,) (z = 0,) (z = x,) (z + x = 4.)
Решение:
Данное тело показано на рисунке.
Из рисунка видно, что основание (R) является квадратом. Для заданных (x, y) значение (z) изменяется от (z = x) до (z = 4 - x.) Тогда объем равен V=∬_R[(4−x)−x]dxdy=∫_02[∫_02(4−2x)dy]dx=∫_02[(4y−2xy)∣_y=02]dx=∫_02(8−4x)dx=(8x−2x2)_02=16−8=8.
Пример 4
Описать тело, объем которого определяется интегралом (V = \int\limits\_0^1 { dx } \int\limits\_0^ { 1 - x } { \left( { { x^2 } + { y^2 } }\right)dy } .)
Решение:
Данное тело расположено над треугольной областью (R,) ограниченной координатными осями (Ox,) (Oy) и прямой (y = 1 - x) ниже параболической поверхности (z = { x^2 } + { y^2 } .) Объем тела равен V=∫_01dx∫_01−x(x2+y2)dy=∫_01[(x2y+3y3)_y=01−x]dx=∫_01[x2(1−x)+3(1−x)3]dx==∫_01(x2−x3+31−3x+3x2−x3)dx=∫_01(2x2−34x3−x+31)dx=(32x3−34⋅4x4−2x2+3x)_01=32−31−21+31=61.
Пример 5
Вычислить объем тела, ограниченного поверхностями (z = xy,) (x + y = a,) (z = 0.)
Решение:
Данное тело лежит над треугольником (R) в плоскости (Oxy) ниже поверхности (z = xy.) Объем тела равен V=∬_Rxydxdy=∫_0a[∫_0a−xxydy]dx=∫_0a[(2xy2)_y=0a−x]dx=21∫_0ax(a−x)2dx=21∫_0ax(a2−2ax+x2)dx==21∫_0a(a2x−2ax2+x3)dx=21(a2⋅2x2−2a⋅3x3+4x4)_0a=21(2a2−32a4+4a4)=24a4.
Пример 6
Найти объем тела, ограниченного поверхностями (z = 0,) (x + y = 1,) ( { x^2 } + { y^2 } = 1,) (z = 1 - x.)
Решение:
Как видно из рисунков, в области интегрирования (R) при (0 \le x \le 1) значения (y) изменяются от (1 - x) до (\sqrt { 1 - { x^2 } } .)
Сверху тело ограничено плоскостью (z = 1 - x.) Следовательно, объем данного тела равен V=∬_R(1−x)dxdy=∫_01[∫_1−x1−x2(1−x)dy]dx=∫_01[(1−x)y∣_1−x1−x2]dx=∫_01(1−x)(1−x2−1+x)dx==∫_01(1−x2−x1−x2−1+2x−x2)dx=∫_011−x2dx−∫_01x1−x2dx−∫_01(1+2x−x2)dx.
Вычислим полученные три интеграла отдельно. I_1=∫_011−x2dx. Сделаем замену: (x = \sin t.) Тогда (dx = \cos tdt.) Видно, что (t = 0) при (x = 0) и (t = \large\frac { \pi } { 2 } \normalsize) при (x = 1.) Следовательно, I_1=∫_011−x2dx=∫_02π1−sin2tcostdt=∫_02πcos2tdt=∫_02π21+cos2tdt=21∫_02π(1+cos2t)dt=21(t+2sin2t)_02π=21(2π+2sinπ)=4π. { Сравните с площадью сектора единичного круга в первом квадранте).
Вычислим второй интеграл ( { I\_2 } = \int\limits\_0^1 { x\sqrt { 1 - { x^2 } } dx } ,) используя замену переменной. Полагаем (1 - { x^2 } = w.) Тогда (-2xdx = dw) или (xdx = \large\frac { { - dw } } { 2 } \normalsize.) Находим, что (w = 1) при (x = 0) и, наоборот, (w = 0) при (x = 1.) Интеграл равен I_2=∫_01x1−x2dx=∫_10w(−2dw)=−21∫_10wdw=21∫_01wdw=21∫_01w21dw=21(32w23)_01=31. Наконец, вычислим третий интеграл. I_3=∫_01(1−2x+x2)dx=(x−x2+3x3)_01=1−1+31=31. Таким образом, объем тела равен V=I_1−I_2−I_3=4π−31−31=4π−32≈0,12.
Пример 7
Вычислить объем единичного шара.
Решение:
Уравнение сферы радиусом (1) имеет вид ( { x^2 } + { y^2 } + { z^2 } = 1). В силу симметрии, ограничимся нахождением объема верхнего полушара и затем результат умножим на (2.) Уравнение верхней полусферы записывается как z=1−(x2+y2). Преобразуя это уравнение в полярные координаты, получаем z(r,θ)=1−r2. В полярных координатах область интегрирования (R) описывается множеством
(R = \left[{ \left( { r,\theta }\right)|;0 \le r \le 1,0 \le \theta \le 2\pi }\right].) Следовательно, объем верхнего полушара выражается формулой V_21=∬_R1−r2rdrdθ=∫_02πdθ∫_011−r2rdr=2π∫_011−r2rdr. Сделаем замену переменной для оценки последнего интеграла. Пусть (1 - { r^2 } = t.) Тогда (-2rdr = dt) или (rdr = - \large\frac { { dt } } { 2 } \normalsize.) Уточним пределы интегрирования: (t = 1) при (r = 0) и, наоборот, (t = 0) при (r = 1.) Получаем V_21=2π∫_011−r2rdr=2π∫_10t(−2dt)=−π∫_10tdt=π∫_01t21dt=π(23t23)_01=32π. Таким образом, объем единичного шара равен V=2V_21=34π.
Пример 8
Используя полярные координаты, найти объем конуса высотой (H) и радиусом основания (R).
Решение:
Сначала получим уравнение поверхности конуса. Используя подобные треугольники, можно записать Rr=HH−z,;;где;;r=x2+y2. Следовательно, H−z=RHr;;или;;z(x,y)=H−RHr=RH(R−r)=RH(R−x2+y2). Тогда объем конуса равен V=∬_Rz(x,y)dxdy=∬_RRH(R−x2+y2)dxdy=RH∬_R(R−r)rdrdθ=RH∫_02π[∫_0R(R−r)drd]dθ=RH∫_02πdθ∫_0R(Rr−r2)dr=R2πH∫_0R(Rr−r2)dr==R2πH(2Rr2−3r3)_r=0R=R2πH(2R3−3R3)=R2πH⋅6R3=3πR2H.