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

Учебное пособие 800292

.pdf
Скачиваний:
7
Добавлен:
01.05.2022
Размер:
1.35 Mб
Скачать

dy

a a2

a a

dy

a

12

a2

a a

 

 

12

12

11 22

;

 

 

 

12

11 22

. (1.8)

 

 

a11

 

 

 

 

a11

 

dx

 

 

 

dx

 

 

 

 

Знак подкоренного выражения определяет тип уравнения (1.2): a122 a11a22 0 гиперболический тип;

a122 a11a22 0 эллиптический тип;

a122 a11a22 0 параболический тип.

Если в некоторой области выражение a122 a11a22 меняет знак, то

уравнение (1.2) называется уравнением смешанного типа в этой области.

Пример 1.3. Определить тип дифференциального уравнения

 

 

 

 

 

1

uy 0.

 

 

 

 

 

 

uxx

yuyy

 

 

 

 

 

Решение. Здесь a

 

2

a2

a a

 

y,

11

1,

a

12

0,

a

22

y и

22

 

 

 

 

 

 

 

12

11

 

т.е. если y 0, то уравнение гиперболического типа,

если

y 0

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

 

 

 

 

 

 

 

 

 

 

 

 

 

Приведем наиболее типичные для каждого типа основные

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

1. Гиперболический тип: волновое уравнение, описывающее малые поперечные колебания струны:

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2u

 

 

2 2u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

utt a uxx

или

 

 

 

a

 

 

 

 

 

 

– одномерное;

 

 

 

 

 

 

 

 

 

t2

 

 

x2

 

 

 

 

u

 

a

2

(u

 

u

 

)

или

2u

a

2

2u

 

 

2u

– двумерное;

 

 

tt

 

xx

yy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

 

 

y2

 

u

tt

a2(u

xx

u

yy

u

zz

)

– трёхмерное,

 

 

где a2

const

(скорость звука).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

2. Эллиптический тип: уравнение Пуассона:

u

 

u

 

f(x,y) или

2u

 

2u

f(x,y).

xx

yy

x2

y2

 

 

 

 

 

3.Если f(x,y) 0, то уравнение называется уравнением Лапласа

иопределяет гармоническую функцию u(x,y):

u

 

u

 

0 или

2u

 

2u

0

– двумерное;

xx

yy

x2

y2

 

 

 

 

 

 

uxx uyy uzz 0 – трёхмерное.

4. Параболический тип: уравнение теплопроводности:

 

 

 

 

 

2

 

 

 

 

 

 

 

 

u

 

 

2

2u

 

 

 

 

 

 

 

 

ut a uxx

 

или

 

a

 

 

 

 

,

 

(a 0)

– одномерное;

 

 

t

 

 

x2

 

 

u

 

a

2

(u

 

 

u

 

 

)

 

 

 

 

u

 

2

2u

2u

t

 

xx

yy

или

 

 

 

 

a

 

 

 

 

 

 

– двумерное;

 

 

t

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

u

t

a

2(u

xx

u

yy

u

zz

)

– трёхмерное,

 

где a2 const .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

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

11

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

(краевых) условий.

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

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

Например, для одномерного однородного волнового уравнения

 

2u

a2

2u

, описывающего малые

 

 

поперечные колебания

 

t2

x2

 

 

 

 

 

 

 

 

 

 

 

 

закрепленной на концах струны, нужно

 

указать отклонение u и

скорость движения

u

 

в начальный момент времени t0 0, то есть

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,t)

 

t 0

u(x,0) f (x);

u

 

(x),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

t 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где f (x), (x) ‒ заданные функции. Так как струна закреплена на концах, то в любой момент времени отклонения в сечениях струны с абсциссами x 0 и x l отсутствуют. Таким образом, граничные условия в этом случае имеют вид

u(0,t) 0,

u(l,t) 0.

12

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

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

Таким образом, при постановке задачи математической физики необходимо:

1)составить уравнение в частных производных;

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

Основные типы задач математической физики

В соответствии с типом дополнительных условий задачи математической физики обычно делят на три основных типа:

1.Задача Коши (задача начальных условий). В ней требуется отыскать решение данного ДУЧП, удовлетворяющее начальным условиям, задающим значения функции (и её производной по времени, если требуется) в начальный момент, т.е. при t 0.

2.Граничная задача. В ней требуется в некоторой области отыскать решение данного ДУЧП, удовлетворяющее заданным условиям на границе этой области.

3.Краевая или смешанная задача. В ней требуется в некоторой области отыскать решение данного ДУЧП, удовлетворяющее некоторым начальным и некоторым граничным условиям.

13

Замечание. Также часто граничную задачу называют краевой. В этом случае краевую задачу обязательно называют смешанной или начально-краевой.

ГЛАВА 2. ОДНОМЕРНОЕ ВОЛНОВОЕ УРАВНЕНИЕ

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

Допустим, что начальное положение струны длины l совпадает с осью OX, левый конец ее помещен в начало координат и колебания проходят в вертикальной плоскости Oxu. Пусть в силу тех или иных причин струна выведена из положения равновесия (рис. 2.1, а). Струна при этом изменит свою форму: каждая ее точка испытывает некоторое смещение. Затем струна начнет колебаться. Для определения положения струны в момент времени t 0 нужно в этот момент указать величину смещения u(x,t)каждой ее точки x 0,l (рис. 2.1, б). Итак, весь процесс колебания струны полностью описывается вещественной функцией двух переменных u u t,x .

u u

u(x,0)

u(x,t)

 

 

M

M

0

 

x

 

l

x

0

x

x+

x l

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

 

 

 

 

 

а)

 

 

 

 

 

Рис. 2.1. Поперечные колебания струны

Будем представлять струну как гибкую упругую нить, на которую не действуют внешние силы. С математической точки зрения гибкость заключается в том, что напряжения, возникающие в струне, всегда направлены по касательной к ее мгновенному профилю (рис. 2.1, б). Кроме того, станем рассматривать малые

14

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

x

случае, как известно из механики, величина натяжения T в каждой точке M струны одинакова и не зависит от времени, причем

проекция

Tx и

 

Tu натяжения T на оси OX и Ou равны

соответственно:

 

 

 

u

 

 

T

x

T cos T ,

T T sin T tg T

,

 

 

 

 

 

u

x

 

 

 

 

 

 

 

где α – угол касательной к мгновенному профилю струны в точке M

(рис. 2.1, б).

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

движения участка MM по оси Ou в момент времени t равна

x x u ,t d ,

x

t

а изменение количества движения за промежуток времени Δt есть

x x

 

u

 

u

 

 

 

 

 

 

,t t

 

,t

d .

 

 

 

t

 

t

 

 

x

 

 

 

 

 

Приравниваем изменение количества движения за промежуток времени Δt к импульсу действующих сил, складывающихся из натяжения

T u x x,t T u x,t

x

x

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

x x

u

 

 

 

 

,t

t

x

 

 

u

t ,t d

t

t t

t

 

u

 

u

 

 

T

 

 

x x,

 

x,

d . (2.1)

 

 

 

x

 

x

 

 

 

 

 

 

 

Станем рассматривать только такие значения u поставленной задачи, которые представляют собой функции, непрерывные вместе со своими частными производными на полосе 0,l 0, и имеющие внутри этой полосы вторые частные производные. Тогда после применения теоремы о среднем формула (2.1) примет вид:

15

 

 

 

 

u

 

 

 

 

u

 

 

 

 

 

u

 

 

 

 

 

u

 

 

 

(x

)

 

 

x

,t t

 

 

 

 

x ,t x

 

 

 

 

x x,t

1

 

 

x ,t

t,

 

 

 

 

 

 

 

 

 

 

 

1

 

t

1

 

 

 

t

1

 

 

 

x

 

 

 

 

x

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

x1 x,x x ,

t1 t, t t ,

а

 

 

после применения

теоремы

Лагранжа:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x

)

2u

x ,t

 

t x

2u

 

x

 

 

,t

 

x t,

 

 

(2.2)

 

 

 

 

t2

 

x2

 

 

 

 

 

 

 

 

 

 

1

 

1

 

2

 

 

 

 

 

2

 

2

 

 

 

 

 

 

где x2 x,x x ,

 

t2 t,t t . Сократив обе части равенства (2.2) на

t x

и переходя к пределу при t 0,

x 0, получим:

 

 

 

 

 

 

 

 

 

(x)

2u

x,t

2u

x,t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

 

 

 

 

 

 

 

то есть искомая функция u должна быть решением однородного линейного уравнения с частными производными 2-го порядка:

(x)

2u

 

2u

.

(2.3)

t2

 

 

 

x2

 

В случае постоянной плотности ρ этому уравнению обычно придают вид

2u

2

2u

,

где

 

T

 

.

t2

 

x2

 

 

 

 

Если на струну действует сила, плотность которой, рассчитанная на единицу массы, равна f, то процесс колебаний описывается неоднородным линейным уравнением

2u

2

2u

f x,t .

(2.4)

t2

x2

 

 

 

Полученные уравнения (2.3) и (2.4) называются простейшими уравнениями колебаний струны. Они, очевидно, принадлежат гиперболическому типу в каждой точке множества 0,l 0, .

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

16

2u

 

2u

2u

 

f x, y,t ,

(2.5)

 

 

2

 

 

 

 

 

 

 

 

2

 

2

 

2

t

 

 

x

 

y

 

 

 

 

 

 

 

 

 

 

 

 

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

2u

 

2u

2u

2u

 

f x, y,z,t .

(2.6)

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

2

 

2

t

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

2.1. Решение задачи Коши для однородного одномерного волнового уравнения методом д’Аламбера

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

Дифференциальное уравнение свободных колебаний струны имеет вид однородного одномерного волнового уравнения

utt(x, t) a2uxx(x, t)

( x ,

t 0).

(2.7)

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

17

u(x, 0) f(x),

ut(x, 0) F(x),

(2.8)

где f(x) и F(x) – заданные при всех значениях x функции. Математически эта задача формулируется так: найти решение

уравнения (2.7), удовлетворяющее начальным условиям (2.8). Это типичная задача начальных условий, или задача Коши.

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

x at,

x at.

(2.9)

Преобразуем (2.7) к новым переменным и . По правилу дифференцирования сложных функций, считая, что

u x,t u (x,t), (x,t) ,

найдём:

ux u x u x u (x at)x u (x at)x

u u ;

ut u t

u t u (x at)t u (x at)t a (u u );

 

ux x (ux )x (u u )x (u )x (u )x

 

(u ) x (u ) x (u ) x (u ) x u 2u u ;

ut t (ut )t

a (u u )t a (u )t -a (u )t

 

a (u ) t

a (u ) t a (u ) t

a (u ) t

 

 

 

=a2 (u 2u u ).

 

 

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

utt a2uxx 4a2u ,

и тогда уравнение (2.7) сводится к уравнению u 0. Переписав

18

его в виде (u ) 0, мы

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

зависит от

 

и,

значит, является функцией только от :

 

 

 

 

 

u ( ).

 

 

(2.10)

Далее, интегрируя (2.10), получаем

 

 

 

 

u ( ) d 1( ) 1( ) 2( ),

 

(2.11)

где 2( ) ( ) d , а

1( ) – произвольная

функция

от

,

играющая при этом интегрировании роль произвольной постоянной. Переходя в (2.11) к старым переменным x и t , получим:

 

u(x, t) 1(x at) 2(x at),

(2.12)

где 1

и 2 – произвольные функции своих аргументов.

 

Полученное решение уравнения свободных колебаний струны содержит две произвольные функции. Его называют общим решением уравнения (2.7).

Для решения поставленной задачи Коши для уравнения (2.7) в найденном общем решении (2.12) определим функции 1 и 2 так,

чтобы это решение удовлетворяло начальным условиям (2.8). Продифференцировав (2.12) по переменной t , получим:

ut(x, t) a ( 1(x at))x at a ( 2(x at))x at .

(2.13)

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

 

u(x, 0) 1(x) 2(x) f(x);

(2.14)

ut (x, 0) a 1(x) a 2(x) F(x).

(2.15)

Проинтегрируем (2.15) по отрезку [0, x ]:

 

 

1

x

 

1(x) 2(x)

F(z) dz C .

(2.16)

 

a

0

 

 

 

 

19