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

Лекция 6. Метод Даламбера

В этой лекции решение задачи Коши для волнового уравнения

Шаг 1. Заменим переменные (x, t) новыми переменными (ξ,η), в которых волновое уравнение примет другой вид: Такая замена выполняется по формулам

После подстановки этих производных в волновое уравнение, получим:

что и требовалось доказать.

Шаг 2. Преобразованное уравнение легко решается двумя последовательными интегрированиями (сначала по переменной η , а затем по ξ):

где C1(η) – произвольная функция от η. Так как C(ξ) – произвольная функция, то и – также произвольная функция.

Окончательно, общее решение U(ξ,η) имеет вид

Шаг 3. Для нахождения общего решения первоначального уравнения подставим в (25) вместо ξ и η выражения (24):

Шаг 4. Определим функции C1 и C2, используя начальные условия из (23). После подстановки первого условия получим

Найдем производную функции U в (26) по переменной t и подставим второе условие:

В результате будем иметь систему уравнений

Если проинтегрировать второе уравнение системы (27) по x в пределах от xo до х , то получим следующую систему:

При сложении этих уравнений получим

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

Подставим теперь полученные функции в общее решение (26):

Поменяем местами пределы интегрирования во втором интеграле, стоящем в скобках в (28). В результате получим решение исходной задачи Коши

Формула (29) называется формулой Даламбера.

Далее мы исследуем решение, определяемое по формуле Даламбера.

Пространственно-временная интерпретация формулы Даламбера

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

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

Такая задача Коши с помощью замены независимой переменной сводится к задаче (23):

Решение преобразованной задачи имеет вид (см. формулу Даламбера (29):

Если теперь в эту формулу вместо τ подставить at, то получится решение исходной задачи

Прежде, чем перейти к физической интерпретации этой формулы, сделаем следующее замечание.

Замечание. Рассмотрим в отдельности функции C1(x-at) и C2(x-at), входящие в общее решение (26) (коэффициент а в них появился потому, что нас сейчас интересует более общее уравнение (30)). Начнем с функции C1(x-at) и построим графики этой функции при возрастающих значениях t: t=to, t=t1, t=t2 и т.д. (см. рис. 8).

Если по очереди проецировать эти картинки на экран (как в мультфильмах), то они «побегут» вправо. Процесс передвижения отклонения по струне называется волной. При этом коэффициент а является скоростью распространения волны. В самом деле, предположим, что параллельно оси х движется наблюдатель со скоростью а. Пусть в некоторый момент to он находился в точке xo. Тогда за промежуток наблюдатель сместится вправо на величину и окажется в точке Если в точке xo наблюдатель видел отклонение струны на величину то в момент t величина отклонения – будет точно такой же! То есть наблюдатель будет видеть форму струны не изменяющейся.

Вторая функция C2(x-at) тоже представляет собой волну, но только она будет распространяться со скоростью а влево. Часто функции C1(x-at) и C2(x-at) называют, соответственно, прямой и обратной волной. Таким образом, общее решение U(x,t) (формула (26)) волнового уравнения является суперпозицией прямой и обратной волны.

Теперь дадим интерпретацию формулы Даламбера для двух частных случаев.

СЛУЧАЙ 1. Предположим, что начальное отклонение отлично от нуля, а начальная скорость равна нулю. Это означает, что начальные условия имеют вид

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

то есть решение U в некоторой точке xo в момент времени to зависит от значений начальной функции φ в двух точках на оси х: в точке (xo — ato) и в точке (xo + ato) (см. рис. 9).

Значение U равно среднему арифметическому значений начальной функции φ в точках (xo — ato) и (xo + ato). На рис. 9 изображена плоскость xOt, которая называется фазовой плоскостью. На оси х указаны точки (xo — ato, 0) и (xo + ato, 0), в которых начальные отклонения струны определяют величину отклонения струны в точке xo в момент времени to. Эти точки являются точками пересечения прямых x — at = xo — ato и x + at = xo + ato с осью х. Указанные прямые называются характеристиками волнового уравнения. Треугольник с вершиной в точке o, to) и основанием, которое получается при пересечении характеристик с осью х (см. рис. 9), называется характеристическим треугольником.

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

Замечание. На самом деле начальные отклонения струны не могут быть разрывными в точках х = -1 и х = 1, ведь струна не разрывается. Однако мы не слишком сильно погрешим против истинной картины распространения колебаний, если будем считать их кусочно постоянными. Дело в том, что, во-первых, рассматриваются очень малые колебания струны, и, во-вторых, малые изменения начальных значений незначительно влияют на решение задачи.

На рисунке 10 изображена фазовая плоскость x0t. Решение U(x,t) задачи отлично от нуля только в заштрихованных областях, причем начальное отклонение распространяется с одинаковой скоростью в двух противоположных направлениях – возникает прямая и обратная волны. Границы этих областей – это характеристики волнового уравнения: x — at = -1, x — at = 1, x + at = -1, x + at = 1.

Если рассмотреть процесс колебания некоторой фиксированной точки струны x = xo, то нетрудно заметить, что она колеблется только в конечный промежуток времени: от момента до момента , то есть В остальное время точка xo находится в покое. Говорят, что в момент t1 через точку x = xo проходит передний фронт волны, а в момент t2 — задний фронт волны. Вообще, фронтом волны называется граница между возмущенной (колеблющейся) и невозмущенной областями среды (точками струны). Для прямой волны уравнение переднего фронта x — at = 1, а заднего фронта x — at = -1. Для обратной волны, соответственно, x + at = -1 — уравнение переднего фронта, а x + at = 1 — заднего фронта.

СЛУЧАЙ 2. Пусть начальное отклонение равно нулю, а начальная скорость отлична от нуля. Это означает, что начальные условия имеют вид

В этом случае решение задачи Коши называют волной импульса. Оно имеет вид (см. формулу Даламбера)

то есть решение U в некоторой точке xo в момент времени to зависит от начальных скоростей ψ во всех точках отрезка [xo — ato , xo + ato] (см. рис 11). Значение U равно (интегральному) среднему значению начальной скорости на отрезке [xo — ato , xo + ato], умноженному на промежуток времени t.

На рис. 11 изображена фазовая плоскость x0t. Точки (xo — ato, 0) и (xo + ato, 0) являются точками пересечения характеристик x — at = xo — ato и x + at = xo + ato с осью х. В качестве примера приведем фазовую картину решения следующей задачи:

Рис. 12 описывает процесс колебания струны, которой сообщается начальная единичная скорость на отрезке -1

При вычислении интеграла всегда удобно представить себе характеристический треугольник с вершиной в точке, лежащей в соответствующей области (см. рис 12). Тогда значение U(x,t) будет определяться значениями начальной функции ψ(x) в основании характеристического треугольника.

2. В области 2 функция

3. В области 3 функция

4. В области 4 функция

5. В области 6 функция

Это решение в различные моменты времени можно изобразить на плоскости x0U (см. рис 13). Здесь для простоты положим a=1.

Графики функции U(x,t), изображенные на рис. 13, задают форму струны в различные моменты времени.

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

Основные типы уравнений

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

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

.

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

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

.

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

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

.

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

В выписанных уравнениях искомая функция u зависит от двух переменных t, x или x, y. Рассматриваются также уравнения и для функций с большим числом переменных. Например, волновое уравнение с тремя независимыми переменными имеет вид

,

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

.

Уравнение колебаний струны.

Формулировка краевой задачи

В математической физике струной называют гибкую упругую нить. Пусть струна в начальный момент времени расположена на отрезке 0≤xl оси Ox. Предположим, что ее концы закреплены в точках x=0 и x=l. Если струну отклонить от первоначального положения, а потом предоставить самой себе или придать ее точкам некоторую скорость, то точки струны будут совершать движение. Задача заключается в определении формы струны в любой момент времени и в определении закона движения каждой точки струны в зависимости от времени.

Если предположить, что движение точек струны происходит перпендикулярно оси Ox и в одной плоскости, то процесс колебания струны описывается одной функцией u(x,t), которая определяет величину перемещения точки струны с абсциссой x в момент t.

Доказано, что при отсутствии внешней силы функция u(x,t) должна удовлетворять дифференциальному уравнению в частных производных второго порядка

.

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

Пусть, например, концы струны при x=0 и x=l неподвижны. Тогда при любом t должны выполняться равенства

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

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

.

Эти два условия называются начальными условиями.

Колебания бесконечной струны.

Формула Даламбера решения задачи Коши

для волнового уравнения

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

Рассматривая свободные колебания, мы должны решить однородное уравнение

при начальных условиях

, ,

где функции f(x) и g(x) заданы на всей числовой оси. Такая задача называется задачей с начальными условиями или задачей Коши.

Преобразуем волновое уравнение к каноническому виду, содержащему смешанную производную. Уравнение характеристик

распадается на два уравнения:

интегралами которых служат прямые

Введем новые переменные ξ=xat, η=x + at и запишем волновое уравнение для переменных ξ и η.

, ,

,

,

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

.

Интегрируя полученное равенство по η при фиксированном ξ, придем к равенству . Интегрируя это равенство по ξ при фиксированном η, получим

,

где φ и ψ являются функциями только переменных ξ и η соответственно. Следовательно, общим решением исходного уравнения является функция

. (8)

Найдем функции φ и ψ так, чтобы удовлетворялись начальные условия:

.

,

.

Интегрируя последнее равенство, получим:

,

где х0 и С – постоянные. Из системы уравнений

Таким образом, мы определили функции φ и ψ через заданные функции f и g, причем полученные равенства должны иметь место для любого значения аргумента. Подставляя в (8) найденные значения φ и ψ, будем иметь

.

Найденное решение называется формулой Даламбера решения задачи Коши для волнового уравнения

Пример. Решить уравнение при начальных условиях , .

Используя формулу Даламбера, сразу получаем

.

Решение волнового уравнения

методом разделения переменных

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

, (9)

удовлетворяющее краевым условиям

u(x,0)=f(x), . (12),(13)

Частное решение уравнения (9), удовлетворяющее граничным условиям (10) и (11), ищут в виде произведения двух функций:

Подставляя функцию u(x,t) в уравнение (9) и преобразовывая его, получим

.

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

, где λ>0. (14)

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

и . (15)

Общее решение этих уравнений

,

,

где A, B, C, D – произвольные постоянные.

Постоянные A и B подбирают так, чтобы выполнялись условия (10) и (11), из которых следует, что X(0)=X(l)=0, так как T(t)≠0 (в противном случае u(x,t)=0). Учитывая полученные равенства, находим

А=0 и .

Так как B≠0 (иначе, было бы X=0 и u=0, что противоречит условию), то должно выполняться равенство

,

.

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

Заметим, что, если в равенстве (14) вместо – λ взять число λ (λ>0), то первое из уравнений (15) будет иметь решение в виде

.

Отличное от нуля решение в такой форме не может удовлетворять граничным условиям (10) и (11).

Зная , можем записать

.

Для каждого n получаем решение уравнения (9)

.

Так как исходное уравнение (9) линейное и однородное, то сумма решений также является решением, и потому функция

(16)

будет решением дифференциального уравнения (9), удовлетворяющим граничным условиям (10) и (11).

Найденное частное решение должно еще удовлетворять начальным условиям (12) и (13). Из условия (12) получим

.

Далее, дифференцируя члены ряда (16) по переменной t, из условия (13) будем иметь

.

Правые части двух последних равенств есть ряды Фурье для функций f(x) и φ(x), разложенных по синусам на интервале (0, l). Поэтому

. (17)

Итак, ряд (16), для которого коэффициенты Cn и Dn определяются по выписанным формулам, если он допускает двукратное почленное дифференцирование, представляет решение уравнения (9), удовлетворяющее граничным и начальным условиям.

Пример. Найти решение краевой задачи для волнового уравнения

, 0


источники:

http://pandia.ru/text/79/052/35879.php