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

книги / Практические занятия по высшей математике. Численное решение алгебраических и трансцендентных уравнений, матричное исчисление, векторный анализ и интегрирование линейных дифференци

.pdf
Скачиваний:
11
Добавлен:
12.11.2023
Размер:
15.42 Mб
Скачать

Дифференцирование обеих частей каждого из этих равенств по х9 у, г и / дает

дох

32у и

_

3*у 1

до2

3*у

“дГ"“ 373?’

д1

З ^ З /’

“ЗГ “

ЗГЗГ:

 

 

 

Зы* — д*? .

д ° г

_

3*у ,

Зг^ __ З^р

 

 

 

 

37

ЗТЗх ’

37

 

3Тду *

"37

Ъг

Подставляя эти значения в выражения для

Уу%И7, из ус­

ловия задачи,

получим

 

 

 

 

 

 

 

 

^ * “ Й + Й ' ! Е + ^ ' ! 5 + й г 5г 1 “

 

 

- & Ш

+ Й ) ‘ + Й П + Й '

 

 

 

 

 

г 1

3*у

д?

,

3*у

Зу

,

3*у

Зу

у

 

 

3^3/ "*

ду дх

Зх

'

ду%* 35

'

З^Зг *

За

 

 

-Ш(1)’+($’+(йЯ+$

I?

*

=

^ Ь -

3*у Зу

Зау

Зу

3®у З у _

 

 

 

З Г л

Зг Зх дх

дгду Щ

3? Зг

Замечая, что

 

-*д а + (8 ,+йп+а-

( I ) ’ + ( й ) ’ + ( 8 ) ’ - ^ + < + * “ л

предыдущие равенства

перепишем в виде

 

 

 

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

. а отсюда на основании определения градиента функции

заключаем, что

- рай ( ^ »• + §*).

Итак, вектор есть градиент некоторой функции. Это значит, что

он вектор

потенциальный.

Задача

12,4. Дифференциальное уравнение неразрывности для

несжимаемой жидкости имеет вид

 

 

 

до$.

д?у дог

0 .

 

 

 

 

 

 

дх ''

ду '

 

дг

 

 

 

 

Доказать,

 

что

при потенциальном движении жидкости потен­

циал скорости — функция 9

удовлетворяет

уравнению Лапласа

 

 

 

 

 

 

4 .

дг* ~

 

^

-

П

 

 

 

 

дх*ду*^

 

 

 

 

 

Р е ш е н и е .

Если предположить,

что движение

жидкости потен­

циальное, то

 

скорость V= V (х, у, г)

жидкости

является градиен­

том некоторой

функции

9 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V=

§гас1 9 ,

 

 

 

 

 

а отсюда проекции

скорости

на оси

прямоугольной системы коор­

динат равны

частным производным

потенциала скорости 9 по соот­

ветствующим

этим осям

координатам,

т. е.

 

 

 

 

 

 

 

_

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V X~Гх'

 

 

 

 

 

 

 

 

 

 

. Дифференцируя

обе

части

первого равенства по х,

второго —

по4#, третьего — по г, получим

 

 

 

 

 

 

 

 

 

 

 

дух

д*?,

 

 

 

 

.

&>,

дЬ?

 

 

 

 

дх

дх* ’

1&у ** ёу**

дг

дг* *

 

 

Складываем эти равенства почленно:

 

 

 

 

 

 

 

 

д ° х

д ° у

, д у ,

д

\

,

д*ч

, а*у

 

 

 

По условию задачи левая часть этого

равенства есть нуль, сле­

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

 

правая его часть

равна

нулю, т. е.

 

 

 

 

 

Й + 0 + 8 = О-

 

 

 

(12 , 10)

З а к л ю ч е н и е .

Потенциал

скорости

несжимаемой

жидкости

удовлетворяет уравнению Лапласа (12,10).

 

 

 

 

Обычно левая часть этого уравнения

обозначается сокращенно

через Д9 , поэтому

(12,10) записывается так:

 

 

 

Ф ункция <р, удовлетворяющая этому уравнению, называется гармонической. В случае плоского потенциаль­ ного движения несжимаемой жидкости уравнение Лапласа (12,11) имеет вид

 

 

 

 

д12у

,

д\

= 0.

 

 

 

( 12, 12)

 

 

 

 

дх* +

ду>

 

 

 

Задача 12,5 (для самостоятельного решения). Доказать, что функ­

ция

<р =

1п г удовлетворяет

уравнению

Лапласа

(12 , 12)

(г =

= V х? + у*).

 

 

 

 

 

 

 

 

 

 

Задача

1 2 ,6. Дано поле вектора

 

 

 

 

 

 

 

 

 

 

а =» 2х • г у / + г й

 

 

 

 

Найти

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

линий.

 

 

 

 

 

Р е ш е н и е .

Дифференциальные уравнения векторных

линий

имеют вид (12,5). В нашем случае

 

 

 

 

 

 

 

 

 

а , = 2х;

ау « у,

аг -

г

 

 

 

и система

дифференциальных уравнений -(12,5) запишется так:

 

 

 

 

 

42

 

 

 

 

 

или

 

 

 

Тх = - у

= Т

 

 

 

 

 

 

 

 

- Л-

 

 

 

 

 

 

 

 

 

 

1)

2)

*

 

 

 

 

 

 

 

9

у '

> 2х

г '

 

 

 

 

Интегрируя,

получаем семейства векторных линий

 

 

 

 

 

 

1) - у Ь * =

1пу + 1п С

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) у 1п х 1п 2 + 1п Са,

 

 

 

 

а отсюда уже

 

 

 

 

 

 

 

 

 

 

1 )

\п{У~ху)ш* 1пС; У х у - С ;

ху*^С Л; У = %,

г д е С ,= С г.

2)

1п Ух = 1п С2г;

У х = Сгг\

х = С\гг

или

гг

С3х,

где

С= —

*Сп­

итак, имеем следующие уравнения семейства векторных линий:

= ^ и г* = СЛх.

Надо иметь в виду, что начало координат является здесь особой точкой: если х-»-0 , то у-*- оо при ( ^ # 0.

Задача 12 ,7. Движение жидкости задано проекциями скорости

У* —Ьу: Уу — кх,

где к — постоянная величина. Требуется найти линии тока и направ­ ление движения.'

Решен| Ге . Если в уравнениях

(12,5) ах, а„,

аг — проекции

скорости частицы жидкости, то эти

уравнения при

установившемся

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

В нашем случае эти дифференциальные уравнения (12,5) запи­

шутся так:

Лх = 0у —ку кх*

Отсюда

кхйх —ку йу.

Сокращая на к и интегрируя, получаем

х' + у* = С

— линии тока — семейство концентрических окружностей с

цент­

ром в начале координат.

 

 

нас Ух =

—ку;

Теперь определим направление движения. У

Уу ■» кх, поэтому

 

 

 

 

У= VУ1+Уу =V *У+кгх*;

V =*к\/хг+ у*;

 

О Т (К '> “ Т “ , - 7 5 = Ь >

» - ■

у Я + р '

 

 

 

 

х

 

 

 

 

х*+ у*

 

Из равенства сох (У, х) =

___ 2— . следует,

что если

х > 0

и у> 0, то соз (Р, х) < 0 и скорость V с положительным направ­ лением оси Ох образует тупой угол. Это говорит о том, что дви­ жение происходит против хода часовой стрелки

Задача 12,7а. Движение жидкости задано проекциями скорости

У, = —* + * У „~у + и

Определить:

1. Уравнение семейства линий тока, а также линию тока, про­

ходящую через точку А (—2, —3) в момент времени 1 = 0.

 

2. Траекторию

частицы

 

жидкости,

которая

в момент времени

( = 0

находилась в точке

А (—2,

—3).

 

 

 

 

 

Р е ш е н и е ,

К

Дифференциальное

уравнение семейства

линий

тока

на основании (12,5) имеет вид

 

 

 

 

 

 

 

 

 

йх

йу

 

 

 

 

 

 

 

 

- х

+ 1 ~

у + Г

 

 

 

 

Считая I — фиксированным и интегрируя, получим

 

откуда

— 1п (—х

 

1) = 1п +

0 — 1пС,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( * - * ) 0 / +

0 “

С.

 

 

(А)

Линиями тока в каждый момент времени

является семейство

гипербол. При

1 = 0 это семейство имеет уравнение

 

 

 

 

 

 

—ху С.

 

 

 

 

 

Подставляя сюда

на основании условия

задачи

координаты точки

А (—2, —3), получим С =

6. Искомым уравнением линии тока,

соответствующим

условию задачи,

будет

 

 

 

 

 

 

 

 

 

ху —6.

 

 

 

 

 

2.

Чтобы ответить на второй вопрос задачи, надо знать следу­

ющее:

 

 

 

 

 

 

 

 

 

 

1. Движение

жидкости'

называется

неустановившимся,

если

проекции скорости Ух, Уу

и Уг являются

функциями не

только

координат, но и времени, т. е. если

 

/);

Уг —У,(*. У>2, 0-

К = Ух{х, у, 2, /); Vу = Уу (х, у, 2,

2. Дифферейциальные

уравнения траекторий жидкой .частицы

в этом случае

имеют вид

 

 

 

 

 

 

'<

 

 

____ йх____________ йд____________ йх

_ ..

 

 

(х, у, г, /) — (ж, у, г. I) ~

Уг (х, у,

/)

 

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

ЗГ " ух (*• У> *. 0;

= Уу(х, у, г, Г>; *

ат в у г (*• у>г>0 -

Каждое

из них — линейное неоднородное уравнение первого по­

рядка с

постоянными коэффициентами.

Перепишем их в виде

Интегрируем эти уравнения по правилам 'интегрирования ли­ нейных дифференциальных уравнений:

х = Схе-< + (— 1; у ^ С ^ — 1— 1.

(В)

Подставляем в эти уравнения / = О, х — —2, у = 3.

Получится,

что Сх = 1 , С* — —2 и уравнения траектории примут вид

х= ~е~' + ( - 1;

уа ~2е1 — 1— 1.

За м е ч а н и е . Сравнивая уравнения (А) и (В), видим, что при

неустановившемся движении жидкости линии тока не совпадают с траекториями жидкой частицы (при установившемся движении жидкости линии тока являются одновременно и траекториями частиц).

Задача 12 ,8. (для самостоятельного решения). Поле скоростей жидкости задано следующим образом:

 

4* ' (** + у * + г*) у X* + у» +

г*'

где (2 — постоянная величина.

 

 

Определить линии

тока.

 

 

О т в е т .

 

 

 

через ось 0г)\

у штСхх

(семейство

плоскостей,

проходящих

г = С2х

(семейство

плоскостей,

проходящих

через ось Оу).

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

другого семейства — прямые, проходящие через

начало координат.

Задача 12,9. (для самостоятельного решения). Напряженность

электростатического поля определяется вектором

г ,

где

е— по­

ложительный электрический заряд, г расстояние

точки

поля до

заряда,

г° — единичный вектор

вектора г,

соединяющего

точку

поля с

зарядом.

 

 

 

 

 

 

 

 

 

Определить векторные линии

поля.

 

 

 

 

 

 

От в е |. Лучи,

выходящие из заряда.

 

 

 

 

 

Задача

1 2 ,10 .

Твердое тело

вращается

с

постоянной

угловой

скоростью о) вокруг некоторой оси. Известно, что

во вращательном

движении

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

 

 

 

0 =

0» х

г,

 

 

 

 

 

 

где г — радиус-вектор точки А (х, у,

г)

тела. Найти

го1 о.

Р е ш е н и е . Направим вектор и»

по

оси

вращения,

которую

примем за ось Ог, в сторону, откуда вращение представляется про­ исходящим против движения часовой стрелки. Проекции векторов

о» и л на координатные оси равны:

и>х = 0;

% = 0;

<ог = о»;

 

гх — *;

Гу = у,

г1 = г.

(А)

Так как вектор о равен векторному

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

ш и

г, то

 

 

 

 

 

1

/

к

 

V= (1>X г =

%

Гг

 

 

г х

ГУ

 

Подставляя значения проекций векторов и» и г из равенства (А), получаем

1 / к

V = 0 0 (1) = —<*>у7+л>*/ ,

X У г

откуда следует, что проекции вектора V равны:

ох = —<оу; ьу — и>х; Уг = 0.

Чтобы воспользоваться формулой (12,6) для определения ротора вектора, найдем входящие в эту формулу частные производные

от проекции вектора о:

до

Ч

дх и’ ду

до

д ? ~ ш' & 0;

ч

до

до.

0.

а ? 1 ду *“ ЗГ-* =

Подставляем эти значения в формулу (12,6):

го(Ъ=

Ь = [«>— (—«>)] • Ь — 2шк.

Итак,

го! V2<л>к.

В связи с тем что вектор ш — приходим к заключению: го( V =» 2ш,

а отсюда следует: ротор линейной скорости точек вращающегося

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

Задача 12,11 (для самостоятельного решения). Най,ти векторные линии поля вектора

О — —0)1/1 + а>ДГ/ + кЛ,

где со и А —* величины постоянные.

Ук а з а н и я .

1.Уравнения векторных линий запишутся так!

— со^ а>*

Л

 

или

 

 

а г - —

з? = Л-

(А)

2. Из первых двух уравнений получаем, дифференцируя обе части каждого из них по I,

г - — * а - - *

3. Общими решениями втих линейных однородных уравнений с постоянными коэффициентами будут:

х «

Ахсое(о/ + Вхз!п<о/;

(В)

у *

созсо/ + з!п(о/.

 

4. Учесть уравнения (А) н показать, что между произвольными постоянными существуют соотношения

Вх —Аъ Вг Ах

итогда уравнения (В) перепишутся так:

х= Л, соз (о/ — Аг51п ш/;

у— АгС05 <о( 4- Ах51П (I)/.

5. Если Л, *= Я соза, а Л = К $ т а , т. е. Я = УлГЙГл!,

у |

= т * т0 уравнения (С) запишутся так:

х = Л соз (а>/ + «);

уЯ $1п ((■>/ 4 - а).

Третье уравнение системы (А) дает

г = М 4 - А,.

Таким образом, введены только три про­ извольные Я, л и А3, как и должно быть.

Уравнения (О) и (Е) определяют винтовые линии, расположенные на ци­

линдре радиуса Я, ось которого совпа­ дает с осью Ог.

Задача 12,12. Найти линейный ин­

теграл вектора а — дЛ — г/4/ вдоль пер­ вой четверти окружности

(С)

а —

О»

(Е)

х = Я соз(; у = Я$1п(.

Р е ш е н и е . Чтобы

найти значение

вектора а

на

заданной ок­

ружности,

заменим в его выражении

х на Ясоз*,

а у — на

Тогда

 

Я1СОЗ* I I

Я4зш4

 

 

Отсюда

а =

 

 

ах «* Я4соз4*; ау = —Я1 зш41.

 

 

 

 

(А)

Чтобы

воспользоваться формулой

(12,3), найдем

йх и йу из

уравнения окружности

х Ясоз 1\

 

 

 

 

 

 

 

 

Отсюда

 

у = ДзШ*.

 

 

 

 

йх — —Я з т / Л ;

 

 

 

 

 

(В)

 

 

йу Я соз I й(.

 

 

Подставляя значения (А) и (В) в формулу (12,3), получаем, учи-

Т

и = { ах6х + ау6у = (* Л4со$4/ (—Я $ ш/) Л —

Задача

12,13

(для

самостоятельного

решения). Найти линейный интеграл

вектора

а — Л — у3] вдоль

первой чет­

верти окружности X— /? С05 (; у = /? 51П и

О т в е т , а *= — у # 4.

 

Задача

12,14.

Найти

циркуляцию

* вектора

а » дя — у)

вдоль

замкнутой

кривой,

образованной

осями координат

 

и

первой

четвертью астроиды

х =

Р е ш е н и е . Линия I

— адсо$3/;

у =

а,ы п8<.

АВ

в

формуле

(12,2) состоит из дуги

астроиды и отрезков АО и ОВ координатных осей. Поэтому

 

ах6х + ау6у = ( ах6х + ау6у +

{ ах6х + ау6у +

 

 

6 а

 

 

Х о

 

+

Сах 6х + ау6у.

 

 

 

ов

 

 

 

 

Каждый из этих интегралов вычислим отдельно: при вычислении первого интеграла по дуге астроиды АВ следует учесть, что так как

а = XI — у/, то ах = х; ау — —у.

На астроиде д: ~ а, со$3 /; у = ад $1П8/, поэтому

 

ах =■ х аусо?8/; ау ~ у = —а 1

$ т 3 1.

Из уравнений астроиды определим и

6.x = — 30| со$* /$1П 161; ■= 3ад $1па I со$ 161.

зоо

Соседние файлы в папке книги