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

9611

.pdf
Скачиваний:
0
Добавлен:
25.11.2023
Размер:
2.94 Mб
Скачать
Ai 1 Ai
n длина каждого отрезка

если, конечно, при кривой будет

стремиться к нулю.

Обобщение этой вычислительной процедуры, отвлеченное от физического содержания, приводит к понятию криволинейного интеграла1-го рода (по длине дуги). Пусть в пространстве в некоторой области D расположена линия AB и пусть в некоторой окрестности этой линии определена функция f (x, y, z) .Разделим линию AB на nучастков

точками A0 A , A1, A2 , A3 ,..., An B .Длину участка линии от точки

Ai 1 до

точки Ai обозначим li . На участке линии Ai 1 Ai выберем некоторую точку

Pi , ηi , νi ) и сформируем следующую интегральную сумму:

 

 

 

n

 

 

 

Sn f i i , νi ) li .

 

 

 

i 1

 

Если существует предел S

n

при n и li 0 , и он не зависит ни от

 

 

 

способа деления кривой AB на n частей, ни от выбора точек Pi , ηi , νi ) наi- ом участке, то этот предел называется криволинейным интегралом 1-ого

рода (по длине дуги) от функции f (x, y, z)

вдоль кривой AB и обозначается как

f (x, y, z)dl . Таким образом,

 

 

 

 

 

 

AB

 

 

 

 

 

 

 

 

f (x, y, z)dl lim S

 

lim n

f (ξ

, η , ν

) l .

(54.2)

n

n

n

i

i i

i

AB

 

 

i 1

 

 

 

 

Возвращаясь к задаче о нахождении массы неоднородной линии, можно с помощью введённого определения записать, что

M f (x, y, z)dl .

(54.3)

AB

 

Будем говорить, что кривая

AB гладкая, если в каждой точке этой

кривойсуществует касательная и угол наклона касательной непрерывно меняется при движении точки вдоль данной кривой.На рис. 54.2 кривая AB гладкая, акривая LQ кусочно-гладкая, ибо в точках M и P касательная не существует.

Рис. 54.2

Теорема. Если непрерывная кривая AB может быть разбита на конечное число гладких кусков и в некоторой окрестности этой кривой

функция f (x, y, z) непрерывна, то для нее существует криволинейный интеграл (54.2).

Наряду со свойствами, которые имеют все рассмотренные ранее интегралы (постоянное число можно выносить за знак интеграла; интеграл от суммы или разности двух функций равен сумме или разности интегралов от этих функций), отметим еще ряд свойств криволинейного интеграла 1-го рода.

а) Криволинейный интеграл 1-го рода не зависит от направления движения по кривой. Действительно, в формуле (54.2), множитель li равен

длине i-го отрезка кривой и поэтому он не зависит от направления.

б) dl L , где L – длина кривой AB . Действительно, если в (54.2)

AB

f (x, y, z) 1, то интегральная сумма будет равна длине кривой от точки Aдо

точки B .

в) Если точка C находится на кривой AB (см. рис. 54.1), то

f (x, y, z)dl

f (x, y, z)dl f (x, y, z)dl .

AB

AC

CB

54.2. Вычисление криволинейного интеграла 1-го рода. Понятно,

что, записав формулу (54.3), мы еще не дали способа вычисления массы, который бы отличался от вычислительной конструкции, приведенной в начале параграфа. Однако анализ конструкции (54.2) показывает, что при

достаточно общих

предположениях относительно свойств функции

f (x, y, z) и кривой

AB вычисление криволинейного интеграла 1-го рода

сводится к вычислению обычного определенного интеграла. При этом существенную роль играет способ задания кривой AB .

Плоский случай.Явное задание кривой. Пусть кривая AB на плоскости

определена уравнением y y(x),

a x b .

Рассмотрим сначала случай,

когда вдоль кривой f (x, y) 1. Ввиду того, что

 

 

 

n

 

L ,

 

dl L lim

 

l

 

n

i

 

AB

 

 

i 1

 

 

получим формулу для вычисления длины дуги AB .

Рис. 54.3

Если отрезки, на которые делится кривая, достаточно малые, то длина li отрезка кривой Ai 1 Ai приблизительно равна длине хорды Ai 1 Ai . По теореме

Лагранжа приращение yi

функции

y y(x)

на участке (xi 1 , xi ) равно

значению производной в некоторой точке

ξ i

этого участка, умноженному

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на приращение xi , т.е. y

y (

ξ

) x . Таким образом, для плоского случая

 

 

 

i

 

 

 

i

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(рис. 54.3) будем иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

A A

 

x2

y2 1 y 2 (

ξi

) x

 

 

i

 

 

i 1

i

 

 

 

 

i

 

 

i

 

 

 

 

 

i

и, следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

L dl lim

1 y 2 (ξi ) xi .

 

 

 

 

 

 

 

 

 

AB

 

 

n

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Li 0

 

 

 

 

 

 

 

 

 

 

 

Нетрудно видеть, что под знаком предела стоит интегральная сумма для функции ( ) = √1 + 2( ) на отрезке (a, b) , и потому последний предел

равен определенному интегралу от этой функции на этом отрезке. Таким образом,

L dl b 1 y 2 (x)dx

AB a

Пусть теперь f (x, y) – произвольная непрерывная функция, определенная в точках близких к кривой AB . Так как, согласно определению, точка Pi на

отрезке кривой Ai 1 Ai может выбираться произвольным образом, то выберем ее так, чтобы она имела координаты (ξi , y'(ξi)) , где ξi есть то значение

аргумента x на отрезке (xi 1 , xi ) , при котором yi y'(ξ i) xi . Отсюда

 

 

n

 

 

 

 

 

 

 

 

 

 

f (x, y)dl lim

f (ξi , yi )) 1 y 2 (ξi ) xi .

 

 

n

i 1

 

 

AB

Li 0

 

 

 

 

 

 

 

 

 

 

 

 

Приглядевшись к выражению под знаком предела, мы видим, что оно

 

 

 

представляет интегральную сумму для функции g(x) f (x, y(x))

1 y'2 (x)

на отрезке (a, b) и потому предел равен определенному интегралу от этой функции на данном отрезке. В результате мы приходим к формуле

 

f (x, y)dl b

 

 

 

 

 

 

f (x, y(x))

1 y 2 (x)dx.

(54.5)

AB

a

 

 

 

 

 

Если кривая AB на плоскости определена уравнение м

= ( ), ≤ ≤ ,

то

 

f (x, y)dl d

 

 

 

 

 

 

f (x( y), y)

1 x 2 (y)dy .

(54.6)

AB

c

 

 

 

 

 

 

 

Пример. Найти массу части кривой y=2-x2 соединяющей

точки A(0, 2) и B(1,1) ,

 

если

плотность распределения массы вдоль

кривой задана функцией

ρ(x, y) 2x (рис.53.4).

Рис. 54.4

Решение. Так как y 2x , то по формуле (54.3) искомая масса будет равна

 

 

 

 

 

 

1

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

2x 1 4x2 dx

1 4x2 d(1 4x2 )

M

 

ρ(x, y)d l

 

 

4

 

 

 

AB

 

 

0

 

 

 

 

 

0

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

((1 4x2 )3 / 2

1

 

 

 

 

 

 

 

 

 

 

 

(5 5 1) 1, 7.

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

0

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Плоский случай. Параметрическое задание кривой.Пусть кривая AB

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

x x(t); y y(t); α t β,

где x(t) и y(t) дифференцируемые функции, производные которых

непрерывны, причем значению t α соответствует точка A, а значению t β соответствует точка B . В этом случае криволинейный интеграл

1-го рода вычисляется по формуле:

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

f (x, y)dl f (x(t), y(t)) x

(t) y

(t)dt.

(54.7)

 

 

AB

 

 

 

 

 

 

 

Для её вывода достаточно в формуле (54.7) осуществить замену переменных

x x(t); a x(α); b x(β), dx xt dt и вспомнить, что yx yt / xtt .

Плоский случай. Уравнение кривой определено в полярных координатах. Попробуем вывести формулу для вычисления

криволинейного интеграла 1-го рода

в случае, когда кривая AB на

плоскости задана уравнением в полярных координатах:

r r( )

 

(54.8)

 

1

2

Вспомним, что связь между декартовыми и полярными координатами определяется формулами x r cos , y r sin . Если в них вместо r

подставить его выражение из уравнения (54.8), то получится параметрическое уравнение кривой AB :

x r( )cos , y r( )sin .(54.9)

вкотором в качестве параметра выступает полярный угол , изменяющийся

впределах от 1 до 2 . При этом

x 2 ( ) y 2 ( ) (r ( )cos r( )sin )2 (r ( )sin r( )cos )2

r2 ( ) r 2 ( ).

Подставляя это выражение в соотношение (54.7), получим формулу для вычисления криволинейных интегралов 1-го рода в случае, когда кривая АВ на плоскости задана в полярных координатах:

 

 

 

 

 

 

 

 

 

f (x, y)dl 2 f (r cos ,r sin )

r2 ( ) r 2 ( )d .

(54.10)

AB

1

 

 

 

Пример. Найти массу половины кардиоиды

 

r (1 cos ),

0 π,

 

 

 

если плотность

ρ r в каждой её точке.

 

Рис. 54.5

Пользуемся формулами (54.3) и (54.10):

M dl r r2 ( ) r 2 ( )d

AB 0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 cos )2

(-sin )2 d

 

 

 

 

 

 

(1 cos )

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

2cos2

2cos

d 8 (1 sin2 ) d (sin ) 8

 

sin

sin

3

16

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

0

2

2

 

 

0

2

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

x x(t); y y(t); z z(t) α t β,

причем параметруt α соответствует точка A, а параметруt β соответствует точка B . Предполагая, что производные x (t), y (t), z (t) непрерывны при α t β, а функция f(x,y,z) непрерывна в некоторой

области D, окружающей кривую AB, и, проводя рассуждения, подобные приведенным для плоского случая, придем к следующей формуле:

 

β

 

 

 

 

 

 

 

 

 

2

(t) y

2

(t) z

2

(t)dx .

(54.11)

 

f (x, y, z)dl f (x(t), y(t), z(t)) x

 

 

AB

α

 

 

 

 

 

 

 

54.3. Некоторые приложения криволинейного интеграла 1-ого рода. Определение площади цилиндрической поверхности. На рис. 54.6

изображено ограждение переменной высоты некоторой горизонтально расположенной территории, ограниченной кривой L. Вопрос о расходе краски, необходимой для окраски этого ограждения, сведется к нахождению его площади. Если ввести систему координат так, чтобы территория находилась в плоскости XOY, то можно считать, что данное ограждение

ρ(x, y, z)

представляет собой часть цилиндрической поверхности с направляющей L и образующей, параллельной оси OZ, причем высота ограждения задается функцией z f (x, y), определенной в каждой точке кривой L(рис.54.6).

Если рассмотреть часть поверхности на небольшом участке кривой Ai 1 Ai , то

можно считать, что высота этой части постоянна и равна значению функции z f (x, y)в некоторой точке Pi (xi , yi ) . В таком случае площадь

этого участка поверхности приблизительно равна Si z(xi , yi ) L i .

Рис. 54.6

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

поверхности ограждения будет определяться по формуле S z(x, y)dl .

L

Вычисление массы, координат центра тяжести и моментов инерции материальной кривой.Часто математическую идеализацию пространственного материального объекта можно представить в виде пространственной кривой AB, вдоль которой распределена масса с линейной плотностью (рис.54.1). (Канаты, на которых

подвешен мост через пролив Босфор, соединяющий Черное и Средиземные моря, можно рассматривать в качестве примера такого объекта. В качестве математической модели канатов можно рассматривать их осевые линии с постоянной линейной плотностью, равной массе тонкого слоя поперечного сечения). В таком случае основные механические характеристики подобных объектов могут быть вычислены с помощью криволинейного интеграла 1- ого рода.

а)Масса материальной кривой. Как мы уже знаем, масса М такой кривой может быть вычислена по формуле (54.3).

б) Моменты инерции I x, I y, I z относительно осей OX , OY , OZ .

Перечисленные моменты вычисляются по следующим формулам

I

x

 

 

(y2 z2 )ρ (x, y, z)d l ,

I

y

 

 

(x2 z2 )ρ (x, y, z)d l ,

I

z

 

 

(x2

y2 )ρ (x, y, z)d l.

(54.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

AB

 

 

 

 

AB

 

 

 

 

AB

 

 

 

Для их вывода, как обычно в таких ситуациях, разделяем кривую на мелкие участки и считаем, что вся масса отрезка кривой (приблизительно равная

ρ(ξ

,η , ν

) L ) сосредоточена

в некоторой ее точке

P

, η

, ν

) . Момент

i

i i

i

 

i

i

i

 

инерции этой материальной точки относительно, например, оси OX, будет равен ее массе, умноженной на квадрат расстояния от оси, т.е. ρ(ξi i , νi ) Li i 2 i 2 ) . Суммируя по всем отрезкам кривой и переходя к

пределу, получаем искомые выражения.

в)Kоординаты центра тяжести. Координаты центра тяжести C(xc , yc , zc )

материальной пространственной кривой могут быть высчитаны по формулам:

 

M

 

 

 

 

 

M

 

 

 

M

 

x

1

 

xρ (x, y, z)d l ,

y

 

 

1

 

yρ (x, y, z)d l , z

 

 

1

 

(z ρ (x, y, z)d l , (54.13)

 

c

 

c

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

AB

 

 

 

 

 

 

AB

 

 

 

 

AB

 

где M – масса кривой.

 

 

 

 

 

 

 

 

 

 

 

Пример.

Найдем

механические

характеристики верхней

полуокружности радиуса R с центром в начале координат, вдолькоторой равномерно распределена масса с плотностью ρ . В данном примере мы имеем дело с плоской кривой и потому в соответствующих формулах мы должны опустить переменную z. Очевидно, что масса M равна длине полуокружности, умноженной на ρ , т.е. M πρR .

Рис. 54.7

Для нахождения моментов инерции Ix и I y сначала запишем уравнение полуокружности в параметрическом виде

x Rcost,

y Rsint,

0 t π,

а затем воспользуемся формулой (54.11) (для плоского случая) и (54.8). Будем иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ix

 

y2ρdl ρy2 (t)

 

x 2 (t) y 2 (t)dt

 

 

 

AB

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ R2 sin2 t

 

R2 sin2 t R2 cos2 tdt

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

ρ R

3

 

 

 

 

 

 

 

 

 

 

 

 

ρ R3 sin2 tdt

(

 

t

sin 2t)

 

 

.

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

4

2

 

 

 

0

 

 

 

 

 

 

 

 

 

Нетрудно понять, что I y имеет такое же значение.

Что касается координат центра тяжести C(xc , yc ) , то в силу симметрии xc 0 , а для вычисления yc опять используем параметрическое уравнение полуокружности:

 

1

 

 

 

 

 

1

 

 

 

 

 

 

yc

 

 

 

yρdl

 

ρy(t) x 2 (t) y 2 (t)dt

 

 

 

 

 

 

 

M

 

 

M

 

 

 

 

 

 

 

 

 

 

 

AB

 

 

0

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2R

 

 

 

 

ρRsin t

 

R

2 sin2 t R2 cos2 t

dt

.

R

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Лекция 55. Криволинейные интегралы 2-го рода (по координатам)

В данной лекции мы рассмотрим еще одну интегральную конструкцию, которая, как и все рассмотренные ранее, изначально возникла при стремлении решить некоторую прикладную задачу (именно в физике),

а потом оказалась применимой при решении других теоретических и прикладных задач. Речь идет о криволинейных интегралах 2-го рода.

55.1. Определение и обозначения. Физическая задача, в которой подобная конструкция возникает, формулируется следующим образом. Пусть в некоторой области D пространства определено силовое поле, т.е. на материальную точку единичной массы, помещенную в область D, действуетзависящая от ее местонахождения вектор-сила F . И пусть в этой области материальная точка перемещается по кривой L из положения A в положение B . Требуется определить работу W сил данного поля при таком перемещении.

Рис. 55.1

На рисунке 55.1 изображено плоское поле сил. Ясно, что при перемещении точки по кривой L по участку кривойACBсилы поля оказывают положительное воздействие, т.е. совершают положительную работу, а при движении по участку кривой DEFони оказывают отрицательное воздействие, т.е. совершают отрицательную работу. Точно также при спуске с горы мы ощущаем положительное воздействие силы тяжести, а при подъеме в гору нам приходится преодолевать отрицательное воздействие этой силы.

Если в пространстве введена декартова система координат, то сила, действующая в точке M (x, y, z) может задаваться посредством ее разложения

по единичным векторам i , j,k

F(x, y, z) P(x, y, z)i Q(x, y, z) j R(x, y, z)k ,

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]