Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Механика_сплошной_среды_

.PDF
Скачиваний:
23
Добавлен:
06.03.2016
Размер:
1.41 Mб
Скачать

 

 

 

t

div( v) 0

(3.9)

Уравнение неразрывности

 

Из этого закона можно получить известную формулу вычисления расхода

V vS const

(3.10)

для невязкой и несжимаемой жидкости, текущей в трубопроводе.

Вывод равенства (3.10) с помощью (3.9) дан в Приложении.

3.6. Фильтрация жидкости в скважину

Просачивание, течение сквозь пористую среду называют фильтрацией жидкости.

Пусть имеется вертикальная скважина радиуса rc . Она пронизывает пласт мощности (толщины) h, содержащий жидкость. За счѐт

перепада давления жидкость течѐт в скважину в радиальном направлении. Скорость фильтрации прямо пропорциональна градиенту давления dp / dr (т.е. разности, перепаду давления, приходящего-

ся на единицу длины) и обратно пропорциональна вязкости жидкости:

v

k

 

dp

,

 

(а)

 

dr

 

 

 

 

 

 

где k коэффициент проницаемости пласта,

r

расстояние от оси

скважины до произвольной точки прискважинного пространства.

Во время эксплуатации скважины прилежащие пласты, как правило, загрязняются. Это означает, что коэффициент проницаемости пластов k уменьшается и потому уменьшается приток жидкости в скважину. Найдѐм количественную связь между этими величинами.

Вокруг скважины мысленно опишем цилиндрическую поверхность радиуса r (r rc ) и высоты h. Площадь этой воображаемой по-

верхности S 2 rh.

(б)

Объѐм жидкости, протекающей сквозь S

за единицу времени, ра-

вен V vS. Он не зависит от r. Подставив сюда значения (а) и (б), получим

V k dp 2 rh.dr

Отделим переменные p и r :

41

dp

 

V

 

dr

.

 

 

 

 

2 h

rk

 

 

 

 

 

 

 

 

 

 

Интегрирование по переменной r в интервале от rc до r

даст

 

 

 

 

 

 

 

r

dr

 

 

p(r)

V

 

 

 

 

pс ,

(в)

 

 

 

r rk

 

 

 

2 h

 

 

 

 

 

 

 

 

 

 

c

 

 

 

где pс

p(rс ) давление жидкости у стенки скважины.

 

В свою очередь, пласт, из которого жидкость притекает в скважину, снабжается жидкостью из области, называемой контуром питания. Ради простоты будем считать, что контур питания представ-

ляет собой окружность радиуса rп

с центром на оси скважины. Дав-

ление на расстоянии rп

обозначим pпл

 

и назовѐм пластовым дав-

лением.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, должно выполняться условие

 

p(rп ) pпл .

Подставим его в

(в). Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

V

 

rп dr

p

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r rk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пл

2 h

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 h

 

 

 

 

 

rп

 

dr

1

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

(3.11)

 

 

 

 

 

V

 

( p

пл

p

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

rk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rc

 

 

 

 

 

Рассмотрим два примера применения данной формулы.

1.

В прискважинной зоне [rс , rп ]

коэффициент проницаемости по-

стоянен, k const. Вычислим интеграл, входящий в (3.11):

rп dr 1 rп dr 1

 

rп 1

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

ln r

 

 

ln

 

.

 

 

 

 

 

rk

k

r

k

k

r

 

 

 

 

r

 

 

 

r

 

r

с

 

 

 

 

 

 

 

 

 

 

 

 

с

 

 

 

c

 

с

 

 

 

 

 

 

 

Подставив это значение в (3.11), будем иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

V

2 hk ( pпл pc )

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

rп

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rс

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получилась известная формула Дюпюи1.

2.

В прискважинной зоне [rс , rп ]

коэф-

фициент проницаемости растѐт по линейному закону:

1 Дюпюи Жюль (1804-1866) – французский инженер, механик, экономист.

42

k ar b (a, b 0)

(рис. 3.11). Подставим эту величину в интеграл, входящий в (3.11):

Рис. 3.11

rп

dr

1

 

r

 

 

rп

1

 

 

(arс b) rп

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

ln

 

 

.

r(ar b)

b a

ar b

 

r

b a

(ar b)r

r

 

 

 

 

 

 

 

 

 

 

 

п

 

с

с

 

 

 

 

 

 

 

с

 

 

 

 

 

 

 

 

Равенство (3.11) принимает вид

 

 

 

 

 

 

 

 

 

 

 

 

V

2 h(b a)( pпл pc )

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

(arc

b)rп

 

 

 

 

 

 

 

 

 

 

 

 

 

(ar

b)r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

c

 

 

 

3.7. Закон парности касательных напряжений

Внутри жидкого или твѐрдого тела в какой-либо точке мысленно построим оси координат x, y, z (рис. 3.12). Вблизи начала координат

проведѐм наклонную плоскость, пересекающую оси так, чтобы треугольник, содержащий точку C, был равнобоким (тогда 3dx 3dy,

значит, dx dy). Получится бесконечно малый тетраэдр, ограничен-

ный координатными плоскостями и наклонной плоскостью (на рис. 3.12 тетраэдр показан в увеличенном масштабе). Точки A, B, C, D, лежащие на гранях, будем считать именами этих граней.

43

Рис. 3.12.

Докажем, что если нет ускоренного вращения различных частей этого тела, то будет выполняться закон парности касательных напряжений:

касательные напряжения на взаимно перпендикулярных гранях, примыкающих друг к другу, одинаковы и направлены либо к ребру, либо от ребра.

Малость тетраэдра позволяет считать, что по грани напряжение распределено равномерно (т.к. оно не успевает измениться). Тогда поверхностная сила, действующая на грань, будет приложена к центру тяжести грани, т.е. к точке пересечения медиан.

Вращение тетраэдра вокруг прямой CD возможно только под действием напряжений xy и yx (рис. 3.11). Напряжение xy создаѐт на

грани B силу xy SB , а напряжение yx создаѐт на грани A силу yx S A.

Чтобы вращения не было, одна из этих сил должна стремиться вращать тетраэдр против часовой стрелки (на рис. 3.11 это сила yx S A ),

другая – по часовой стрелке (сила xy SB ). Поэтому суммарный момент этих сил относительно CD

M z DA yx SA DB xy SB

должен быть равен нулю. Тогда

DA yx SA DB xy SB .

44

t произвольный момент вре-

Но DA = DB, S A S B . После сокращений останется

yx xy ,

что и требовалось доказать.

ПРИЛОЖЕНИЕ

1. Скалярное и векторное поля

Пусть M – произвольная точка пространства, мени. Введѐм обозначения:

(M ,t) плотность частицы жидкости или газа в точке M в момент времени t;

скорость частицы жидкости или газа в точке в момент времени t. v(M ,t) M

Эти обозначения говорят, что плотность и скорость – функции от точки пространства и времени. Кроме плотности и скорости существуют другие переменные величины, также зависящие от положения точки M и момента времени t.

Поле – функция от точки и времени.

Так как функция может быть скалярной либо векторной, то и поле бывает скалярным или векторным.

Значит, плотность (M ,t) скалярное поле (скалярная функция, скаляр), а

скорость векторное поле (векторная функция, вектор). v(M ,t)

Поле, не меняющееся во времени, называют стационарным или устано-

вившимся.

 

 

В этом случае для плотности вместо (M ,t) пишем

(M ), а для скорости

 

 

 

вместо v(M ,t) пишем v(M ).

 

Введѐм систему координат Oxyz. Произвольная точка M будет иметь коор-

динаты x, y, z, поэтому вместо (M ,t) и

 

можно писать (x, y, z,t) и

v(M ,t)

 

 

 

v(x, y, z,t).

 

 

45

Далее символом f будем обозначать произвольную скалярную функцию

 

 

 

 

 

 

(скалярное поле, скаляр), а символом F

– произвольную векторную функцию

(векторное поле, вектор).

 

 

 

 

 

 

 

 

 

 

 

Три координаты вектора F

в пространстве

Oxyz будем обозначать P, Q, R

 

 

 

 

 

 

(три функции от x, y, z, t). В этом случае

F

Pi

Qj

Rk или в краткой записи

F (P, Q, R).

2. Вектор площадки. Поток векторного поля

На поверхности (S) мысленно выделим бесконечно малую площадку

(dS)

(рис. 2.1). Построим вектор площадки

 

 

dS – вектор, выходящий из (dS)

и на-

 

 

 

 

правленный перпендикулярно к (dS) так, что из конца dS виден обход границы

площадки против часовой стрелки.

 

 

 

 

нормальный к (dS) (рис. 2.2). Тогда

 

Построим единичный вектор n,

 

 

 

 

(2.1)

 

dS n dS.

Пусть имеется поле скоростей

 

 

 

v. Поместим в этом поле неподвижную пло-

щадку (dS), проницаемую для жидкости.

 

Рис. 2.1

 

Рис. 2.2

Когда площадка (dS)

перпендикулярна векторному полю

 

v (рис. 2.3, при

 

 

 

 

 

этом векторы v

и dS параллельны), то (dS) пересекает наибольшее число век-

 

 

 

сквозь

(dS) максима-

торных линий. В этом случае говорят, что поток поля v

лен. Он равен dV v dS

– объѐму цилиндра высотой v

и площадью основания

dS.

 

 

 

 

 

 

Рис. 2.3

 

 

 

 

 

 

 

 

 

 

Рис. 2.4

На рис.

2.4 показан другой крайний случай, когда площадка параллельна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

векторному полю v (при этом векторы v

и dS перпендикулярны). Здесь векто-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ры v скользят вдоль (dS), не пересекая эту площадку. Поэтому поток вектора v

сквозь (dS)

равен нулю.

 

 

 

 

 

 

 

 

 

 

 

 

Если

 

и

 

располагаются

 

 

 

 

 

 

 

 

 

v

dS

 

 

 

 

 

 

 

 

 

под произвольным углом

(рис.

 

 

 

 

 

 

 

 

 

2.5), то справедлива общая фор-

 

 

 

 

 

 

 

 

 

мула

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.2)

 

 

 

 

 

 

 

 

 

 

dV v dS.

 

 

 

 

 

 

 

 

 

 

 

Скалярное

произведение

в

 

 

 

 

 

 

 

 

 

правой

части

можно

записать

 

 

 

 

 

 

 

 

 

подробней:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV v dS cos .

 

 

 

 

 

 

 

 

 

 

 

 

Проинтегрировав (2.2) по по-

 

 

 

 

 

 

 

 

 

верхности (S ), получим

 

 

 

 

 

 

 

 

 

 

Рис. 2.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V v

dS

(2.3)

 

 

 

 

 

 

 

 

 

 

 

(S )

 

 

 

 

 

 

 

 

 

 

 

сквозь поверхность (S ).

 

поток векторного поля v

 

Если

 

 

то

 

 

v dS.

Если при этом и значение v

во всех точках (S)

v (S),

v

dS

одинаково, то из (2.3) будем иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V v S.

(2.4)

Если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v – скорость течения жидкости, то формулы (2.3) и (2.4) дают расход

жидкости: объѐм жидкости, протекающей сквозь (S) за единицу времени.

 

 

 

 

 

3. Инвариантные определения

 

 

 

 

 

градиента, дивергенции, ротора

Пусть M произвольная точка пространства.

 

Градиент скалярной функции f f (M )

определяется по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dS f

 

 

 

 

 

 

 

 

 

grad f

 

(dS )

 

 

.

(П.3.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дивергенция векторной функции

 

F F (M ) определяется по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dS

F

 

 

 

 

 

 

 

 

 

 

 

(dS )

 

 

 

 

 

 

 

 

 

 

 

 

 

div

F

 

 

 

 

 

.

(П.3.2)

 

 

 

 

 

 

 

 

 

 

dV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ротор векторной функции

F

F (M ) определяется по формуле

47

dS F

 

 

(dS )

 

 

rot F

 

.

(П.3.3)

dV

 

 

 

 

В этих формулах:

dV объѐм, содержащийся внутри бесконечно малой поверхности (dS);

dS вектор площадки поверхности (dS);

точка ( ) – знак скалярного умножения; косой крест ( ) – знак векторного умножения.

Правые части данных формул не содержат координат переменной точки M . Значит, эти формулы справедливы в любой системе координат, они не зависят от выбора системы координат. Иными словами, они инвариантны.

Правые части этих формул можно записать, используя знак предела. Например, (П.3.3) можно записать так:

 

 

 

 

 

dS

F

 

( S )

 

 

 

rot F (M ) lim

 

.

(П.3.4)

 

 

 

 

( V ) M

V

 

Здесь (V ) M означает, что и область (V ), и замкнутая поверхность (S)

сжимаются в точку M. Поэтому величины (П.3.1)-(П.3.3) не зависят от формы поверхности (dS) или (S ). Если в формуле (П.3.4) убрать знак предела, полу-

чится приближѐнное равенство

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dS

F

 

 

rot

 

( S )

 

 

 

 

F

 

 

,

 

V

 

 

 

 

 

 

или равенство

 

 

 

 

 

 

 

 

 

 

(rot F ) V dS

F,

(П.3.5)

 

 

 

( S )

 

 

которое тем точнее, чем меньше (V ).

 

 

 

 

4. ТЕОРЕМА О ГРАДИЕНТЕ, ДИВЕРГЕНЦИИ, РОТОРЕ

 

Теорема. Если область (V ) ограничена поверхностью (S ), то

 

(grad

 

 

 

 

 

f ) dV dS f ,

(П.4.1)

(V )

 

 

 

(S )

 

 

(div

 

 

 

 

F )dV dS

F ,

(П.4.2)

(V )

 

 

(S )

 

 

(rot

 

 

 

 

 

 

F )dV dS

F.

(П.4.3)

(V )

 

 

(S )

 

 

Формально эти равенства получаются сразу из определений (П.3.1)-(П.3.3).

Равенство (П.4.2) называется формулой Остроградского-Гаусса.

Строгое доказательство этих формул опирается на тот факт, что поверхностные интегралы в правых частях изменяют свой знак при изменении ориентации поверхности.

48

5. Градиент, дивергенция и ротор

вдекартовых координатах

Вдекартовой системе координат градиент, дивергенция и ротор

определяются по формулам

grad f

f

 

 

f

 

 

f

 

x

i

y

j

z

k ,

 

 

 

 

 

 

div F P

Q

R ,

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

z

 

(П.5.1)

(П.5.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot F

 

 

 

 

,

(П.5.3)

 

 

 

 

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

Q

 

R

 

 

 

 

 

 

 

 

где

f f (x, y, z) скаляр (скалярная функция, скалярное поле),

 

 

 

 

 

 

вектор (векторная функция, векторное поле)

F

F (x, y, z) Pi

Qj

Rk

( P, Q, R – функции от x, y, z ).

В П.3 мы выяснили, что формулы (П.3.1)-(П.3.3) не зависят от формы замкнутой поверхности (dS), ограничивающей фигуру (dV ). Поэтому ради удобства

возьмѐм фигуру (dV ) – бесконечно малый параллелепипед со сторонами dx, dy, dz. Его объѐм dV dxdydz. Направим координатные оси вдоль рѐбер (рис. 5.1).

Докажем формулу (П.5.1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(П.3.1) (grad f )dV dS

f

 

 

 

 

 

 

 

 

 

 

(dS )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n dS) f

= Поверхность (dS)

 

 

 

 

 

(dS )

 

 

 

 

 

 

 

 

 

 

 

 

состоит из шести граней =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n f

dS

n f dS

 

 

 

 

 

 

 

 

 

( ЛИРУЛ )

 

 

( АСОКА)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n f

dS

n f

dS

 

 

 

 

 

 

 

(СОРИС)

 

( АКУЛА)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n f dS

 

n f dS

 

 

 

 

 

 

 

 

 

( КУРОК)

 

( АЛИСА)

 

 

 

 

 

 

 

 

 

= На грани (ЛИРУЛ) имеем

 

 

 

 

 

внеш-

 

 

 

 

 

 

 

 

 

 

 

 

ний вектор

 

 

 

 

 

 

 

Рис. 5.1

 

n

i , на противополож-

 

 

 

ной грани (АСОКА) будет

 

 

и т.д. =

 

 

 

 

n

i

 

 

 

 

 

 

f

 

 

f

 

 

 

 

 

 

 

f dS

 

i

dS i

dS j

f dS j

n f

dS k

f dS k

 

( ЛИРУЛ )

( АСОКА)

 

(СОРИС)

( АКУЛА)

 

( КУРОК)

( АЛИСА)

49

= Грань (ЛИРУЛ)

мала, поэтому в пределах этой грани функция f не успе-

вает измениться. Значит, на ней

f f (Л ) const. Аналогично, на грани (АСО-

КА) будет

f f ( А) const. И т.д. =

 

 

dS

 

 

 

f dS

i f (Л )

i f ( А)

dS jf (С)

dS jf ( А)

( ЛИРУЛ )

( АСОКА)

(СОРИС)

( АКУЛА)

kf (К )

(

 

f dS

dS kf ( А)

КУРОК)

( АЛИСА)

= Площадь грани (ЛИРУЛ) равна dy dz и равна площади грани (АСОКА). И т.д. =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i f (Л )dydz i f ( А)dydz

jf (С)dxdz jf ( А)dxdz

kf (К )dxdy kf ( А)dxdy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f ( А)]dxdy

 

 

 

 

 

i [ f (Л ) f ( А)]dydz

j[ f (С) f ( А)]dxdz k[ f (К )

= При перемещении от точки А к Л изменяется лишь координата x,

поэтому

 

 

 

 

 

 

 

 

 

 

 

 

f (Л ) f ( А)

x f

f

dx. И т.д. =

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

dxdydz

f

 

f

 

f

f

f

 

 

 

 

i

 

j

y

dydxdz k

z

dzdxdy i

x

j

y

k

dV .

 

 

 

 

 

x

 

 

 

 

 

 

 

z

 

 

 

 

 

 

dV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученное

равенство

 

 

 

f

f

f

даѐт

нам

формулу

(grad f )dV i

x

j

y

k

dV

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

(П.5.1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Остальные две формулы можно вывести таким же путѐм, но мы поступим иначе. Выпишем полученные равенства

 

 

 

 

f

f

f

(grad f )dV

dS

f i

x

j

y

k

dV .

 

 

 

(dS )

 

 

 

z

 

 

 

 

 

 

 

 

Заметим, что в (а) вектор

заменяется в правой части на векторы

dS

Посмотрим на определение дивергенции:

(а)

 

 

 

i ,

j , k .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(div

F )dV (П.3.2)

dS

F.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(dS )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В отличие от выражения (а) здесь присутствует вектор F

вместо скаляра f . По-

этому в (а) заменим f

 

на

 

и обычное умножение заменим на скалярное, обо-

 

F

значаемое точкой. Будем иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(div

 

 

 

 

F

 

F

F

 

 

 

F )dV

dS

F

i

 

j

 

 

k

 

dV ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(dS )

 

 

 

 

x

 

 

y

 

z

 

 

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

F

 

F

 

(i F)

 

( j F)

 

(k

F)

.

div F

i

x

j

y

k

z

 

x

 

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q,

Но i

F

i

(Pi

Qj

Rk ) P,

аналогично,

j

F

 

 

 

 

 

 

div F P

Q R ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

z

k F R, поэтому

50