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

Учебник по УМФ первая часть

.pdf
Скачиваний:
103
Добавлен:
16.04.2015
Размер:
968.33 Кб
Скачать

Циркуляция и роторвекторного поля

Если в векторном поле, заданном выражением (2), выбрать некоторую линию L и определенное направление на ней, то вектор дифференциала дуги задается выражением

ds = dxi + dyj + dzk .

Циркуляцией вектора А вдоль замкнутого контура называется криволинейный интеграл по этомуконтуруот проекции вектора А на вектор ds, т.е.

As (M ) ds =

A ds = (Axd x + Ayd y + Azd z) .

(9)

L

L

L

 

Рассмотрим теперь предел отношения циркуляции по плоскому контуру L, окружающему точку М, к площади, ограниченной этим контуром, при условии, что контур L стягивается в точку М, оставаясь в одной и той же плоскости:

As (M ) ds

lim

L

 

.

 

S

S0

 

Чтобы вычислить этот предел, следует преобразовать выражение для циркуляции, воспользовавшись формулой Стокса:

 

 

 

 

 

 

 

 

 

 

 

As (M ) ds

= (Axd x + Ay d y + Az d z) =

 

 

 

 

 

L

 

 

L

 

 

 

 

 

 

 

 

 

 

A

 

Ay

 

A

A

Ay

 

A

 

= ∫∫

z

 

 

cosα +

x

z cos β +

 

x cosγ dσ.

 

z

x

D

y

 

 

 

z

x

 

 

y

 

При стягивании контура L в точку P значение подынтегральной функции будет стремиться к её значению в этой точке. Это подынтегральное выражение является

скалярным произведением вектора нормали n = (cos α, cos β, cos γ) кплоскости, в которой лежит контур L, и вектора, проекции которого равны

 

A

Ay

 

 

A

A

 

Ay

 

A

 

z

 

 

,

 

x

z

,

 

 

x .

z

x

 

y

 

 

 

z

x

 

 

 

y

Этот вектор и называется ротором или вихремвектора А(М):

 

A

Ay

 

A

A

Ay

 

A

 

rot A =

z

 

i +

x

z j+

 

x k .

(10)

z

x

 

y

 

 

z

x

 

 

y

 

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

С помощью определения ротора теорему Стокса можно сформулировать в векторной форме:

∫∫rotn A(V )dσ = As (M )ds.

(11)

S

L

 

Если ротор векторного поля равен нулю, то векторное поле является безвихревым или потенциальным. Более подробно это будет обсуждаться в главахVI и VII.

§ 3. Оператор Гамильтона и дифференциальные операторы второго порядка

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

=

i +

j +

k .

(12)

 

x

y

 

z

 

Правила действий с этим вектором таковы:

1.Произведение вектора на скалярную функцию U(x,y,z) дает градиент этой функции

u =

u

u

u

 

x i +

y j+

z k = grad u.

(13)

 

 

 

 

 

2.Скалярное произведение вектора на векторную функцию А(x,y,z) дает дивергенцию этой функции

 

A

Ay

 

A

 

A(x, y, z) =

x +

 

+

z =div A(M ).

(14)

y

 

x

 

z

 

3.Векторное произведение вектора на векторную функцию А(x,y,z) дает ротор этой функции

i j k

∂ ∂ ∂

x y z Ax Ay Az

 

A

Ay

 

A

A

 

Ay

 

A

=

z

 

i +

x

z

j+

 

x k = rot A . (15)

z

x

 

y

 

 

z

x

 

 

 

y

Перейдем теперь к наиболее важным для нас дифференциальным операторам второго порядка:

1. Оператор Лапласа, который образуется путем последовательного применения оператора градиента и оператора дивергенции кскалярной функции U (x,y,z). В результате получаем скалярную функцию

u =

2u

+

2u

+

2u

 

 

 

z2 .

(16)

x2

y2

2. Оператор, представляющий собой последовательное применение оператора ротора и оператора дивергенции кскалярной функции U (x,y,z). Нетрудно убедиться, что в результате мы получим нулевой вектор:

rotgradU = 0 .

(17)

3. Оператор, представляющий собой последовательное применения оператора дивергенции и оператора ротора к векторной функции А(x,y,z). Нетрудно убедиться, что в

этом случае мы в результате получим скаляр:

divrot U = 0 .

(18)

Представленные выше операторы записаны в декартовой системе координат

4. При записи волновых уравнений иногда используют так называемый оператор Даламбера или волновой оператор

u

1

u

(19)

a2 tt

5. Оператор u – вектор:

u = grad divu rot rot u =

6.Оператор, представляющий собой последовательное применение оператора ротора и оператора дивергенции кскалярной функции U (x,y,z). Нетрудно убедиться, что в

результате мы получим вектор, равный нулю, а именно:

rot grad U = 0

(20)

7. Оператор, представляющий собой последовательное применения оператора дивергенции и оператора ротора квекторной функции А (x,y,z). Нетрудно убедиться, что в этом случае мы в результате получим скалярную величину, равную нулю:

divrot U = 0

(21)

В заключение заметим, что при записи уравнений в частныхпроизводныхта часть, в

которую входят частные производные, часто рассматривают как дифференциальный оператор. При этом используют запись L[u]или Lu, например

 

 

u

 

L[u] =

 

k(x)

 

q(x)u

 

 

x

x

 

Гл а в а II. Одномерное волновое уравнение

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

распространении звука). В отсутствии внешнихсил волновое уравнениеимеет вид:

2u

= c2

 

2u

+

2u

+

2u

,

(1)

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

2

 

2

t

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где с является постоянной скоростью распространения волны.

Мы начнем знакомство с волновым уравнением с одного из самыхпростыхслучаев, а именно с задачи о малыхпоперечныхколебанияхструны без учета затухания.

§1. Уравнение малых поперечных колебаний струны

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

движения струны с этими свойствами с учетом действия внешней силы. Строго говоря, заранее мы не можем говорить, что это будет уравнение колебаний, поскольку только

решение уравнения покажет нам, каким будет движение.

Расположим ось x вдоль линии натяжения бесконечной струны в состоянии равновесия, тогда перпендикулярная ей ось u будет параллельна поперечным смещениям

этой струны (Рис. 5). Сила натяжения Т для всех точекструны будет одинакова по величине и направленав каждой из этихточекпо касательной кструне.

Рис. 5. Расположение струныи силынатяженияпри отклоненияот оси x.

Пусть ρ (x) – линейная плотность струны, а p(x,t) – удельная сила (в данном случае приходящаяся на единицу длины), приложенная перпендикулярно струне в ее исходном положении. Будем рассматривать поперечные смещения струны и считать их настолько

малыми, что можно будет отбрасывать величины второго порядка малости по сравнению с ux . В этом случае элемент струны dx после отклонения от положения равновесия будет иметьдлину

 

 

 

 

 

 

 

 

u 2

 

 

ds = (dx)

2

+(du)

2

= dx

1

dx .

(2)

 

 

+

 

 

 

 

 

 

 

 

 

x

 

 

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

Какая же сила будет стремиться возвратить каждую точку струны в положение равновесия? Дело в том, что при отклонении от положения равновесия в каждой точке будет появляться составляющая вектора Т по оси u, абсолютная величина которой будет

равна будет равна T sin α . Это и будет сила, под действием которой точки струны будут стремиться вернуться в положение равновесия.

Рис.6.К выводу уравненияколебаний струны

Посчитаем теперь суммупроекций всехсил, приложенных кэлементу струны ds, на ось u (Рис. 6). Проекция силы Т на ось u будет равна:

T sinα =T du

=T

 

 

 

du

 

 

T

u .

 

 

 

 

 

 

ds

 

dx

1

 

u

2

 

x

 

 

+

 

 

 

 

 

 

 

 

 

 

x

 

 

 

А равнодействующая сил, приложенныхна концахэлемента ds, составит:

 

u

 

u

 

2u

 

T

 

T

 

=T

 

 

 

dx .

x

2

 

x x+dx

 

x x

 

 

x

 

 

 

 

 

 

 

 

ср

 

Сила, приложенная кэлементуds, будет равна p(x,t)dx. По второмузаконуНьютона сумма приведенных выше сил должна быть равна массе элемента ds= ρ(x)dx, умноженной на его

ускорение 2u , т.е.

t2

T 2ux2 xср

dx + p(x,t)dx = ρ(x) 2u dx .

t2

Разделив все члены этого уравнения на ρ и на dx, устремляя одновременно dx к нулю, мы и получим уравнение колебания струны:

T 2u

2u

= −

1

p(x,t).

(3)

 

 

 

 

ρ x2

t2

ρ(x)

 

 

 

 

В случае, если ρ – величина постоянная, то можно ввести величину a = Tρ , тогда уравнение (3) примет вид:

a

2 2u

2u

= −

p(x,t)

,

(4)

 

x2

t2

ρ

 

 

 

 

 

 

или:

u

a2u

xx

=

p(x,t)

.

(4.a)

 

tt

 

 

ρ

 

 

 

 

 

 

Обозначив p(x,t) через F(x,t) можем переписать это уравнение следующим образом:

ρ

u a2u

xx

= F(x,t).

(4.б)

tt

 

 

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

При отсутствии внешней силы это уравнение будет иметь вид:

u tt = a2u x x ,

(5)

Это уравнение является уравнением свободных колебаний струны и отражает

простой факт пропорциональности вертикального ускорения в каждой точке струны utt её вогнутости в этой точке uxx с коэффициентом пропорциональности a2.

Начальные условия

В случае бесконечной струны нужно сформулировать только начальные условия в следующем виде:

u(x,0) =ϕ1(x), ut (x,0) =ϕ 2(x) .

(6)

С физической точки зрения первое из этих условий задает смещение каждой точки струны, а второе – скорость каждой точки струны в начальный момент времени. Задача

отыскания решения уравнения (4) или (5) при выполнении условий (6) называют задачей Коши.

§ 2. Случай ограниченной струны.

Рассмотрим теперь струну длины l. Для ограниченной струны наряду с начальными условиями надо сформулировать и граничные условия. Существует три основных типа

граничныхусловий.

Граничные условия I рода. В этом случае задаются перемещения на концах струны как функции времени (Рис. 7)

u(0,t) =η1(t)

(7)

u(l,t) =η2 (t)

 

Если на концах струна неподвижно закреплена, то в условиях (6) η1 (t) =η 2 (t) = 0 . Такие условия называют однородными, а само неподвижное закрепление – защемлением.

u

Рис.7.Граничные условияпервого рода

Рис. 8. Граничные условиявторого рода

Граничные условия II рода. В этом случае задаются проекции силы на ось u на концах струны как функции времени. Как мы выяснили, каждая из них равна Tu x , поэтомуна концах задаются производные ux как функции времени (Рис. 8)

ux (0,t) =φ1

(t)

(8)

u x(l,t) =φ 2(t)

 

Если концы струны свободны, хотя и перемещаются только по вертикали, то в условиях(7) g1(t) и g2(t) будут равны нулю.

Граничные условия III рода. В этом случае лучше начать с механического смысла этих условий. Сначала предположим, что концы струны связаны с осью x (см.рис.13), тогда со стороны каждой из пружин на концы струны по закону Гука будут действовать

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

условия на концахследующим образом:

ux (0,t) = k1 u(0,t)

 

 

T

,

(9)

 

k2

u x(l,t) =

u(l,t)

 

T

 

 

 

 

где k1 и k2 – жесткости пружин, прикрепленныхкструне на концахструны. Такое граничное условие часто называют упругим закреплением струны в граничныхточках

Теперь рассмотрим общий случай, когда один конец пружины прикреплен к струне,

а второй движется по вертикали по заданномузаконувремени θ (t) (Рис. 9). Тогда условия (9) изменятся следующим образом:

 

u

x

(0,t) = k1

[u(0,t) θ (t)]

 

 

 

T

1

(10)

 

 

 

 

 

u x(l,t) = k2

 

 

[u(l,t) θ2 (t)]

 

 

 

 

T

 

 

После

переноса слагаемых и переобозначений мы получим граничные условия III

рода в традиционном виде

 

 

 

 

 

ux (0,t) γ1u(0,t) = g1(t)

(11)

 

 

 

u x(l,t) γ2u(0,t) = g2 (t)

 

 

 

 

где γi = ki ,

gi (t) = − ki

θi (t) (i =1,2)

 

T

T

 

 

 

 

Рис. 9.Граничные условиятретьего рода (частный случай).

Рис. 10.Граничные условиятретьего рода (общий случай).

К перечисленным видам граничных условий нужно добавить граничное условие, которое описывает действие на конце струны вертикальной силы, пропорциональной

скорости и направленной в противоположном направлении. Такую силу обычно называют

силой вязкого трения. Она может создаваться гидравлическим демпфером, соединяющим конец пружины с неподвижной точкой на оси x. Условие такого типа записывается (для

левого конца) следующим образом:

 

Tu x (0,t) = −but (0,t)

 

 

(12)

или

u x (0,t) + βut (0,t) = 0

, где β =

b

(13)

T

 

 

 

 

В некоторых физических задачах приходится прибегать к использованию на разных

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

условие защемления, а на правом упругое закрепление.

§ 3. Решение задачи Коши

Формула Даламбера

Чтобы выяснить, как будет вести себя бесконечная струна в отсутствии внешних сил, необходимо решить задачу Коши для уравнения (5) с начальными условиями (6), которые мы перепишем в следующем виде:

u tt = a2u x x ,

 

(14)

u(x,0) = f (x),

,

(15)

 

= g(x).

ut (x,0)

 

 

Одно из условий (15) может быть нулевым, но не оба, ибо в этом случае струна будет оставаться в исходном, недеформированном состоянии.

Решение задачи Коши осуществим в три этапа:

1-й шаг (приведение уравнения(7)кканоническомувиду): Введем новые переменные ξи ηследующим образом:

ξ = x + at

(16)

η = x at

 

и запишем равнение (14) в новыхпеременных. Для этого получим выражения для uxx и utt:

ux = uξ ξx +uη ηx = uξ +uη ,

uxx = uξξ + 2uξη +uηη ,

u

= u

ξ

+u

 

η

= c(u

u

),

t

ξ

t

 

η

t

 

ξ

η

 

u

= a2 (u

 

2u

+u

).

 

 

tt

 

ξξ

 

 

ξη

ηη

 

 

 

После подстановки полученныхвыражений висходное уравнение (14) получим:

uξη = 0.

(17)

Уравнение (17) является каноническим видом уравнения (14), а переменные ξ и η

каноническими переменными.

2-й шаг (решение преобразованного уравнения):

Проинтегрируем полученное уравнение по ξ, а затем по η. После интегрирования по ξ, получим:

uη (ξ,η) = F(η) ,

где F(η) – произвольная функцияот η. После интегрирования по η получим

u(ξ,η) = F(η)dη +ψ (ξ) =ϕ (η) +ψ (ξ),

где φ(η) и ψ(ξ) – произвольные функции своихаргументов.

3-й шаг:

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

u (x,t) =ϕ (x at) +ψ (x + at) .

(18)

Физический смысл полученного решения состоит в том, что оно представляет собой

суперпозицию двух волн, распространяющихся вдоль оси x в противоположных направлениях со скоростью a, как это показано на рис.11, причем формы этих волн,

определяемые функциями ϕ и ψ, являются произвольными.

4-й шаг (Использование начальныхусловий):

Для конкретизации формы воспользуемся начальными условиями. Для этого подставим полученное решение и его производную по времени в начальные условия (15).

Рис. 11.Распространяющиеся волны

В результате получим два уравнения относительно неизвестных функций φ(x) и

ψ(x):

u (x,0) =ϕ (x) +ψ (x) = f (x),

ut (x,0) = −aϕ '(x) + aψ '(x) = g(x).

Проинтегрировав второе уравнение от x0 до x, получим систему уравнений относительно искомых функций:

ϕ (x) +ψ (x) = f (x),

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

aϕ(x) + aψ (x) = g(ς) dς +C,

 

 

 

 

 

 

x0

 

 

 

 

где С – произвольная константа. Решая этусистему, получим:

 

1

 

 

1

 

x

 

C

 

ψ (x) =

f (x) +

 

 

g(ς) dς +

 

,

 

 

 

 

 

 

 

2

 

 

2a x

 

2a

 

 

 

0

 

 

 

 

 

1

 

 

1

 

x

 

C

 

 

ϕ(x) =

f (x)

 

g(ς) dς

 

.

 

 

 

 

 

 

2

 

2a x

2a

 

 

 

0

 

 

 

 

Записывая полученные решения для любого времени t, будем иметь

 

1

 

1

 

x+at

 

 

C

 

ψ (x + at) =

f (x + at) +

 

g(ς) dς +

 

,

2

 

2a

 

2a

 

 

 

x

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

1

 

1

 

xat

 

 

C

 

 

ϕ(x at) =

f (x at)

 

g(ς) dς

 

.

 

 

 

 

 

 

 

2

 

2a

x

 

2a

 

 

 

 

 

 

 

0

 

 

 

 

 

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

 

1

[ f (x at) + f (x + at)]+

1

x+at

 

u (x,t) =

g(ς) dς.

(19)

2

2a

 

 

xat

 

 

 

 

 

 

Эта формула носит название формулы Даламбера. Она была получена в предположении, что решение поставленной задачи Коши существует, а её наличие

доказывает и единственность решения. Можно показать и непрерывную зависимость решения от начальных данных. В самом деле, для любого ε > 0 можно указать такое δ > 0,

что если заменить f и g наf1 и g1,так, что

f (x) f1(x) <δ, g(x) g1(x) <δ (−∞ < x < ∞) ,