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

Введение в математическую физику. Учебно-методическое пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
932 Кб
Скачать

Легко убедиться, чтопри y < 0 – этоуравнениегиперболического типа; приy > 0 – эллиптическоготипа; а при y = 0 – параболического типа.

При y < 0 уравнения характеристик принимают вид

1 .

Поэтому кривые

 

 

 

и

 

 

 

являются в

 

 

данном случаехарактеристиками уравненияТрикоми. Преобразование

3

 

 

 

3

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

виду:

приводит уравнение Трикоми кξ,каноническому

η

 

1

 

 

 

.

 

 

∂ξ η

6 ξ η

 

ξ

η

 

Если y > 0, то решение системы уравнений (2.2) можно записать в следующем виде:

 

3

 

x

2 iy 2

C .

 

3

 

Выбрав в качестве новых переменных

ξ, η ,

получим:

.

Задачи

В задачах 2.1–2.5 определите тип уравнения 2-го порядка и приведите его к каноническому виду.

2.1.

2.2.

2

 

2

2

 

3

2

2

 

 

6

 

0.

2

2

 

 

2

2

 

 

 

2

 

 

 

 

 

2

 

4

 

 

 

5

 

 

2

 

 

 

2

 

 

 

0.

 

 

 

 

 

 

 

 

 

 

11

2.3.

2.4.

2.5.

2

 

6

 

2

 

2

 

 

 

 

0

.

 

 

 

 

 

 

2

2

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

5

 

.

 

2

2

 

2

 

 

 

 

2

 

 

10

225

 

 

 

0

.

 

2

2

 

 

 

 

2

 

α

 

 

 

β

 

 

0

 

 

 

 

 

 

 

 

 

2.6. Найдите область гиперболичности, параболичности и эллиптичности уравнения

22

2

2

 

2

2

0.

 

 

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

2.7.

2.8.

2

 

 

 

 

2

 

 

.

 

22

 

 

 

22

 

 

2

 

 

2

0

0.

 

 

 

Вопросы

1.Запишите общий вид линейного дифференциального уравнения второго порядка.

2.Какие типы линейных дифференциальных уравнений с частными производными второго порядка вам известны?

3.

К какому типу принадлежит уравнение, если

0?

Каков ка-

нонический вид уравнения параболического типа?

 

5.

Каков канонический вид уравнения гиперболического типа?

12

Практическое занятие 3

ОСНОВНЫЕ УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ И ОБЛАСТИ ИХ ПРИМЕНЕНИЯ

Теоретическая часть

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

I.Волновое уравнение:

. (3.1)

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

II.Уравнение теплопроводности, или уравнение Фурье:

. (3.2)

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

III. Уравнение Лапласа:

0. (3.3)

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

Вуравнениях (3.1) и (3.2) искомая функция u зависит от четырех,

ав уравнении (3.3) от трех независимых переменных. Однако в ряде задач о колебаниях струны или мембраны число переменных сокра-

щается. Так, волновоеуравнение, описывающееколебанияструны в плоскости, принимает вид

13

.

(3.1’)

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

.

(3.2’)

Уравнение теплопроводности в стационарном случае переходит в уравнение Лапласа. Это уравнение в случае трех независимых переменных решается аналитически для случая сферической симметрии и то только в некоторых частных случаях. В свою очередь, некоторые трехмерные задачи можно привести к задачам двумерным, описывающимся уравнением Лапласа с двумя переменными:

0. (3.3’)

Такие уравнения решаются относительно просто методом разделения переменных.

Задачи

3.1.Запишите волновое уравнение для случая трех независимых переменных.

3.2.ЗапишитеуравнениеФурье для случая трех независимых пе-

ременных.

3.3. Запишитеи решитеуравнениеЛапласа дляодномерного случая.

Вопросы

1.Какие уравнения относятся к уравнениям математической физики?

2.Прирешении какихфизических задачиспользуютсяуравнения гиперболического типа?

3.Прирешении какихфизических задачиспользуютсяуравнения параболического типа?

4.Прирешении какихфизических задачиспользуютсяуравнения

эллиптического типа?

14

Практическое занятие 4

ВЫВОД И РЕШЕНИЕ УРАВНЕНИЙ ГИПЕРБОЛИЧЕСКОГО ТИПА, ПРИМЕРЫ ЗАДАЧ

Теоретическая часть

4.1.Вывод и решение уравнения колебания струны

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

длины L в начальный момент направлена по отрезку оси 0X от 0 до L (рис. 4.1). Предположим, что концы струны закреплены в точках X = 0 и X = L. Если струну отклонить от еепервоначального положения, а потом предоставить самой себе, то струна начнет колебаться, т.е. точки струны будут совершать движения. Задача заключается в определении формы струны в любой момент времени и в определении закона движения каждой точки струны в зависимости от времени.

Мы будемрассматриватьтолькомалыеотклоненияточекструны от начального положения в одной плоскости (плоскость задана системой координат – осями 0X и 0U). В силу этого можно предположить:

1) движение точек струны происходит перпендикулярно оси 0X, при этом процесс колебания струны описывается одной функцией u(x, t), которая дает величину перемещения точки струны с абсциссой x в момент времени t (см. рис. 4.1);

2) длина элемента струны MM’ равняется ее проекции на ось 0X,

т.е. MM’ = x.

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

Рассмотрим элемент струны MM’. На концах этого элемента по касательным к струне действуют силы Т (см. рис. 4.1). Пусть каса-

тельныеобразуют с осью 0X углы и + . Тогда проекция на ось 0U сил, действующих на элемент MM’, будет равна sin φ

15

∆φ

sin.

φ

. Так как угол мал, можно предположить, что

sinφ

tg φ

Вследствие чего, будем иметь*:

 

 

 

 

sinφ

tg φ

∆φ .

tg φ

 

 

sin φ ∆φ

(4.1)

 

 

 

 

 

 

∆ ,

 

,

 

 

 

* Последнее равенство вытекает из геометрического смысла производной. Как известно, u’(x) при данном значении аргумента x равняется тангенсу угла, образованного касательной к графику функции u(x) в соответствующей точке M(x, y) с положительным направлением оси 0X. Таким образом, tg( ) = u’(x).

Применив формулуЛагранжа к выражению, стоящемув скобках, получим:

 

∆ ,

,

0 θ

,

,

∆ .

(4.2)

 

Согласно теореме Лагранжа

1.

 

 

,

(4.3)

где

, т. е. приращение функции

наξ

интервале равно’ ξ

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

ξ

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

U

 

 

 

M

T

 

 

 

M

 

+

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x, t)

 

 

 

 

 

 

 

 

 

0

 

x

x+ x

 

L

X

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

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

– линейнаяплотность струны. Тогда масса элемента струны будет

16

x. Ускорение элемента равно

 

 

 

. Следовательно, по принципу

 

 

 

Даламбера будем иметь:

 

 

 

 

 

 

 

 

 

 

 

 

 

Сокращая на

и

 

ρ∆

 

 

= T

 

 

 

.

(4.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, получаем уравнениедвиже-

ния:

 

обозначая

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

.

 

 

 

(4.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, уравнениеколебанияструны есть волновоеуравнение. Однако для полного определения движения струны одного уравнения (4.5) недостаточно. Искомая функция u(x, t) должна удовлетворять ее граничным условиям, указывающим, что делается на концах струны (x = 0 и x = L), иначальным условиям, описывающим состояниеструны в начальный момент (t = 0). Совокупность граничных и начальных условий называется краевыми условиями.

Так как по условию задачи концы струны жестко закреплены в точках X = 0 и X = L (и потому неподвижны), при любом t должны выполняться равенства

0,

0.;

(4.6)

(4.6’)

граничными условиями нашей задачи.

Эти равенства и являются,

0

 

Вначальный моментt = 0 струнаимеет определеннуюформу, которую ей придали. Пусть этаформа определяетсяфункцией f (x). Та-

ким образом, должно быть

 

.

(4.7)

 

 

 

В начальный момент

должна быть задана скорость

в каждой

,0

|

 

точке струны, которая определяется функцией (x). Таким образом, должно быть

 

 

.

(4.8)

Условия (4.7) и (4.8) являются

начальными условиями.

φ

 

есть постоянная. Ее

Для однородной струны величина

 

 

значениезависит от Т и ρ (т.е. от материала струны и еегеометрических параметров).

17

Будемискатьрешениеэтогоуравненияввиде

).

Подставляя это выражение в уравнение (4.5), получим,

:

или*

’’

 

’’

),

 

 

’’

’’

.

 

(4.9)

* В левой части этого равенства стоит функция, которая не зависит от х, а в правой части – функция, не зависящая от t. Равенство (4.9) возможно только в том случае, когдалеваяиправаячастинезависят ниотx, ниотt, т.е. равны постоянному числу:

 

’’

 

 

’’

const

.

(4.10)

 

 

 

 

 

 

 

2

2

> 0.

Обозначим константу в уравнении (4.10) через

, где

 

Тогда:

 

’’

 

λ ;.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

’’

λ

 

 

 

 

Изэтих равенств получаем два линейных однородных дифферен-

циальных уравнения второго порядка

0; .

(4.12)

’’

λ

 

Общиерешения этих’’

уравненийλ

 

 

 

(4.11)

будут

0выглядеть следующим об-

разом*:

cosλ

 

sinλ ;.

(4.13)

 

 

 

 

(4.14)

 

cosaλ

 

 

sinaλ

 

* Решение однородного линейного дифференциального уравнения второго порядка заключается в составлении характеристического уравнения и определении его корней. Характеристические уравнения для дифференциальных уравнений (4.11) и (4.12) запишутся соответственно:

k2+ 2 = 0;

cosλ

sinλ 2 + 2а2 = 0.

Легко убедиться, что корни полученных квадратных уравнений:

, λ,

,λ.

18

Так каквобоих случаях корнихарактеристическогоуравнения– чистомнимые, решением дифференциального уравнения (4.11) будет функция v(x) вида (4.13), а решением уравнения (4.12) – функция w (t) вида (4.14).

Следовательно, общее решение уравнения (4.5) принимает вид

При этом

 

cosλ

,

)

 

 

sinaλ

.

(4.15)

 

sinλ

cosaλ

 

 

 

 

функция u(x) должна удовлетворят условиям (4.6) и

(4.6’). Иными словами, необходимо, чтобы

 

 

 

 

0

. Под-

ставляя значения x = 0 и x = L в равенство

(4.13), получаем:

 

 

 

0

0,

 

 

 

0

С

 

1 =

С

;

0

 

cosλ

 

sinλ

;

 

 

0

 

 

 

 

 

 

C1

0;

 

 

 

 

sinλ

0

.

 

 

 

C2 0, так какв противном случаебыло бы

 

 

, чтопро-

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

 

0

, откуда

 

 

 

 

 

λ

,

 

 

 

sinλ

0 (4.16)

где (n = 1, 2,…)*.

* Значение n = 0 не берется, так как в этом случае

0.

Найденные значения λ называются собственными значениями для данной краевой задачи. Соответствующие им функции называ-

ются собственными функциями. Подставив для каждого значения n и учитывая, что С1 = 0, получаем набор решений уравнения (4.5), удовлетворяющих граничным условиям (4.6) и (4.6’) *:

,

sin

cos

 

sin

 

.

(4.17)

 

 

* Для каждого значения n необходимо брать свои С2, С3 иС4. В формуле (4.17) постоянная С2n включена вС3n и C4n.

Так как уравнение (4.5) является линейным и однородным, то и сумма егорешений является решением, поэтомуфункция, представленная рядом

;

;

sin

π

cos

π

 

 

sin

,

(4.18)

также будет решением дифференциального уравнения (4.5).

19

Решение (4.18) должно еще удовлетворять начальным условиям (4.7) и (4.8). Этого можно добиться подбором постоянных С3n и C4n. Подставив в равенство (4.18) t = 0, получим

sin

.

(4.19)

Определим коэффициенты Фурье для ряда (4.19):

sin

.

(4.20)

Далее необходимо продифференцировать члены равенства (4.18) по t и подставить t = 0. Из условия (4.8) получается равенство:

φ

 

sin .

(4.21)

 

С4n определяются как коэффициенты Фурье для ряда (4.21):

φ

sin

.

(4.22)

Таким образом, ряд (4.18), гдекоэффициенты С3n иC4n определяютсяпоформулам(4.20) и(4.22), являетсярешениемуравнения(4.5) и удовлетворяет граничным и начальным условиям.

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

Поставить краевую задачу, соответствующую данной физической задаче, – значит выбрать функцию, характеризующую физический процесс, а затем:

вывести дифференциальное уравнение для этой функции;

установить для нее граничные условия;

сформулировать начальные условия.

Рассмотрим распространение продольных волн в однородном стержнедлиной L сплотностью . Стержень имеет цилиндрическую форму с площадью поперечного сечения S (рис. 4.2).

Согласно закону Гука на элемент стержня х действуют силы F1 = S (х) слева и F2 = S (х+ х) справа, где (х) – нормальноеупругое напряжение. Таким образом, на элемент х действует результирующая сила

F = S (х + х) – S (х).

(4.23)

20

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