Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Учебные материалы ИКНТ 1-8 семестры / 7 семестр / !Фирсов А.Н. (сост.) Уравнения математической физики (учебное пособие для заочников).pdf
Скачиваний:
50
Добавлен:
21.03.2016
Размер:
419.37 Кб
Скачать

21

 

 

4

p q

 

mpx

 

npy

 

 

 

Amn

=

 

ò òf0 (x, y) sin

 

 

sin

 

 

 

dxdy,

 

p

 

q

 

 

pq 0 0

 

 

 

 

 

 

 

 

 

4

p q

 

mpx

 

 

npy

 

Bmn

=

 

ò òf1 (x, y) sin

sin

dxdy.

 

 

 

 

 

 

akmn pq 0 0

 

 

p

 

 

 

q

Примеры.

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

u(x, y, t)

 

t =0

=

 

 

l

 

sin

px

 

sin

py

. Начальная скорость равна нулю. Вдоль контура

 

 

 

 

 

 

 

 

 

100

 

 

 

l

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мембрана закреплена.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. В рассматриваемом случае

 

 

 

 

 

 

j0 (x, y) =

 

l

sin

px

sin

py

,

j1 (x, y) = 0. Следовательно,

 

 

 

 

 

l

 

 

 

100

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

Bmn = 0, m =1,2,..,

 

n =1,2,...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

l l

l

px

py

mpx

 

npy

 

 

 

 

 

 

 

 

Amn =

 

òò

 

sin

 

sin

 

sin

 

sin

 

dxdy.

 

 

 

 

 

 

 

l 2

100

l

l

l

l

00

Всилу ортогональности тригонометрической системы функций только A11 ¹ 0,

а все остальные

 

 

 

 

 

 

 

 

4

 

æ l

 

 

 

px

 

 

ö2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

= 0.

 

A

=

 

ç

 

sin2

dx÷

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

11

 

100l ç ò

 

 

 

 

l

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è 0

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 æ

1 l

æ

 

 

2px

ö

ö2

 

 

 

1

 

æ

 

l

 

 

 

l

 

 

 

2px

 

l

ö2

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

ç

 

 

ç1

- cos

 

 

 

÷dx ÷

=

 

 

 

 

ç x

 

 

-

 

sin

 

 

 

 

 

÷

=

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100l ç

2 òè

 

 

 

l

ø

÷

 

 

 

 

100l ç

 

0

 

 

 

2p

 

 

 

l

 

 

 

÷

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

0

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

 

0

ø

 

 

 

 

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

u(x, y,t) =

 

 

l

 

cos

ap

 

 

2

 

t sin

px

sin

py

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

l

 

 

 

 

 

 

 

l

 

 

 

 

l

 

 

 

 

 

 

8. Метод разделения переменных (общий случай).

Пусть требуется найти решение уравнения

 

2u

 

æ

u ö

- q(x)u

 

r(x)

 

 

=

 

ç p(x)

 

÷

(8.1)

t

2

 

 

 

 

x è

x ø

 

 

(где r(x), p(x), q(x) - достаточно гладкие функции, причем p(x)>0, r(x)>0,

 

 

 

 

22

 

q(x)³0), удовлетворяющие условиям

 

 

ì

 

u(0, t)

 

 

ïau(0, t) + b

 

 

= 0,

 

x

 

ï

 

 

(8.2)

í

u(l, t)

 

 

ïg u(l, t) + d

= 0.

 

 

 

ï

 

x

 

 

î

 

 

 

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

u(x,0) = j(x),

u(x,0)

=y (x) (0 £ x £ l).

(8.3)

 

 

t

 

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

 

 

 

 

 

 

u(x, t) = T (t) X (x) .

(8.4)

Подставляя (8.4) в уравнение (8.1) получим

 

 

r(x)T"(t) X (x) = T (x)

d

æ

 

 

 

 

dX ö

 

 

 

ç p(x)

 

 

 

÷ - q(x)T (t) X (x) или

 

 

 

 

 

 

 

 

 

 

 

 

dx è

 

 

 

 

dx ø

 

 

d æ

dX ö

- q(x) X (x)

 

 

 

 

 

 

 

 

 

 

 

 

ç p(x)

 

÷

 

 

T"(t)

 

 

 

 

 

 

 

 

 

 

 

 

dx è

dx ø

 

 

 

 

 

 

 

=

= -l,

где l - некоторая постоянная. Отсюда

 

 

r(x) X (x)

 

 

 

 

 

T (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d æ

 

 

 

dX ö

 

 

 

 

 

 

 

 

 

 

 

 

ç p(x)

 

 

 

 

÷ + (lr(x)

- q(x)) X = 0, (8.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx è

 

 

 

 

dx ø

 

 

 

 

 

 

 

 

 

 

 

T "(t) + lT (t) = 0

 

 

(8.6)

Так как T (t) ¹ 0, то для того чтобы функция (8.4) удовлетворяла краевым условиям (8.2), необходимо и достаточно выполнение условий:

ìaX (o) + bX '(0) = 0,

í

(8.7)

îgX (l) + dX '(l) = 0.

 

Таким образом, для определения функции X(x) мы пришли к следующей краевой задачи для обыкновенного дифференциального уравнения:

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

Свойства собственных значений и собственных функций краевой задачи: 1. Существует счетное множество собственных значений l1 < l2 < ... < ln < ...,

которым соответствуют собственные функции X1 (x), X 2 (x),...

 

23

2. При q(x) ³ 0 и (p(x) X n

(x) X n' (x))

 

x =l £ 0 все собственные значения ln поло-

 

жительны.

 

 

x =0

 

 

 

 

 

3.Собственные функции на отрезке [0, l] образуют ортогональную с весомr(x) и нормированную систему:

l

 

ì0,

при

m ¹ n,

ò r(x) X n

(x) X m

(x) =í

при

(8.8)

0

 

î1,

m = n.

4. (Теорема Стеклова). Всякая функция f (x) , удовлетворяющая краевым условиям (8.7) и имеющая непрерывную первую производную и кусочно-непрерывную вторую производную, разлагается в абсолютно и равномерно сходящийся ряд по соб-

 

 

 

 

 

¥

 

 

 

 

 

l

ственным функциям X n (x) :

f (x) = åan X n (x),

an = ò r(x) X n (x) f (x)dx.

 

 

 

 

 

 

 

n =1

 

 

 

0

 

Далее, для каждого собственного значения ln решаем уравнение (8.6). Общее ре-

шение уравнения (8.6) при l = ln (обозначим его через Tn (t) ) имеет вид

Tn (t) = an cos

 

t + bn sin

 

t, где an

и bn -произвольные постоянные.

ln

ln

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

(8.1) вида un (x, t) = Tn (t) X n (x) = (an cos

 

 

t + bn sin

 

t) X n (x).

 

ln

ln

Чтобы удовлетворить начальным условиям (8.3), составим ряд

¥

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x, t) = å(an cos

ln

t + bn sin

ln

t) X n (x)

(8.9)

n =1

 

 

 

 

 

 

 

 

 

 

Если этот ряд сходится равномерно, так же как и ряды, получающиеся из него двукратным почленным дифференцированием поx и t, то сумма его будет удовлетворять уравнению (8.1) и краевым условиям (8.2).

Для выполнения начальных условий(8.3) найдем постоянные an и bn как коэффициенты разложения функций j и y в обобщенные ряды Фурье по ортонор-

мированной (с весом) системе функций (X n ).

Теперь можно сделать некоторые общие замечания относительно области применения метода разделения переменных.

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

A(x)

2u

+ B(t)

2u

+ C(x)

u

+ D(t)

u

+ (F (x) + F (t))u = 0 или приводиться

 

 

 

 

 

x2

 

t 2

 

x

 

t

1

2

 

 

 

 

 

 

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

24

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

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

9. Вывод уравнения теплопроводности для стержня.

Рассмотрим однородный теплоизолированный с боков стержень конечной длины l, имеющий постоянную по длине толщину, и настолько тонкий, чтобы в лю-

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

Выберем ось x (направив ее по оси стержня) так, чтобы стержень совпадал с отрезком [0;l] оси x.

 

 

 

 

 

 

o

 

 

 

 

x

x

 

 

 

x + Dx

Рисунок 5.

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

Выделим элемент стержня [x,x+Dx] и составим для него уравнение теплового баланса, согласно которому скорость изменения количества тепла в рассматриваемом объеме (изменение количества тепла в единицу времени), обусловленная теплоемкостью материала, равна количеству тепла, поступившему в этот объем в единицу времени вследствие теплопроводности.

Скорость изменения тепла в выделенном элементе стержня равна

x + Dx

crs

u ( x, t )

dx ,

 

ò

где с - теплоемкость материала стержня; r - плотность;

 

x

 

t

 

s - площадь поперечного сечения. Применяя к этому интегралу теорему о среднем,

x +Dx

u(x, t)

 

u(x +q1Dx,t)

 

 

 

получим ò cr s

dx = cr s

Dx, где 0

<q1

<1.

t

t

x

 

 

 

 

мая функция u(x,t)
u(x, t)

25

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

времени,

 

равно - k

 

u(x, t)

s ,

где k - коэффициент теплопроводности, а s - пло-

 

 

 

щадь сечения.

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

u(x, t)

 

æ

 

 

 

 

u(x + Dx, t)

 

ö

æ

u(x + Dx, t)

 

u(x, t) ö

 

- ks

 

 

 

 

s

- ç

- k

 

 

 

s ÷

= ksç

 

 

 

 

-

 

 

÷

=

 

x

 

x

 

 

x

 

 

 

x

 

 

 

è

 

 

 

 

 

 

 

ø

è

 

 

 

 

 

ø

 

= ks

2u(x +q

2Dx, t)

Dx,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где 0 <q2 <1. ( Здесь применяется формула конечных приращений Лагранжа к

функции

 

u(x, t)

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Составим уравнение теплового баланса

 

2u(x +q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

crs

u(x +q

Dx, t)

Dx = ks

 

2

Dx, t)

Dx.

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разделим обе части этого уравнения на sDx (объем выделенного элемента стержня) и устремим Dx к нулю (стягивая выделенный элемент стержня к сечению). Получим

u(x,t)

 

2 2u(x, t)

æ

2

 

k

ö

 

 

= a

 

 

 

ça

 

=

 

÷.

(9.1)

 

 

 

2

 

 

t

 

 

x

ç

 

 

 

÷

 

 

 

 

è

 

 

cr ø

 

Это уравнение называется уравнением теплопроводности для однородного стержня.

k

Величина a = называется коэффициентом температуропроводности. Иско-

cr

должна удовлетворять уравнению (9.1), начальному условию

t =0 = u(x,0) = f (x) (0 £ x £ l), где f(x)- заданная функция от x (это усло-

вие выражает закон распределения температуры по длине стержня в начальный момент времени t =0), и граничным условиям

u(x, t)

 

 

= u(0,t) = j (t),ü

 

 

 

x =0

 

1

ï

(0 £ t £ +¥), где j1 (t) и j2 (t) - заданные функции

 

 

 

u(x, t)

 

 

= u(l, t) = j

 

 

ý

 

 

 

(t). ï

 

 

 

x =l

 

2

 

þ

 

 

 

 

 

 

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

26

10. Первая краевая задача для уравнения теплопроводности.

Пусть W - конечная область трехмерного пространства. Обозначим через Q

в пространстве переменных (x,y,z,t) цилиндр, основание которого есть область W и образующие которого параллельны оси Ot. Пусть QT -часть этого цилиндра, ограниченного снизу плоскостью t=0 и сверху плоскостью t=T(T>0).Часть границы

цилиндра

QT , состоящую из его нижнего основания (t=0) и боковой поверхности,

обозначим через Г .

 

 

 

 

 

в цилиндре QT решение уравнения

 

 

Рассмотрим

 

 

следующую задачу: найти

теплопроводности

 

 

 

 

 

 

 

 

 

 

u

 

æ

2

u

 

 

2

u

 

2

u

ö

 

 

= a

2 ç

 

 

 

 

 

÷

 

 

 

 

 

 

2

+

 

 

 

 

 

 

2 +

 

 

2

(10.1)

 

t

ç

x

 

y

z

÷,

 

 

è

 

 

 

 

 

 

ø

 

удовлетворяющее начальному условию

 

 

 

 

 

 

 

 

 

u

 

t =0 = j(x, y, z)

(10.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и граничному условию

 

 

 

 

(t Î[0,T ]),

 

 

 

 

 

 

 

 

 

u

 

s

=y (P, t)

(10.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р - точка поверхности S. Функцииj и y непрерывны,

где S- граница области W,

причем значения y при t=0 совпадают со значениями j на границе S.

Задача нахождения решения уравнения (10.1) при условиях (10.2), (10.3) называется первой краевой задачей для уравнения теплопроводности.

Теорема. Функция u(x,y,z,t), удовлетворяющая однородному уравнению теплопроводности (10.1) внутри цилиндра QT и непрерывная вплоть до его границы, принимает наибольшее и наименьшее значения наГ, т.е. или при t=0, или на боковой поверхности цилиндра QT .

Из этой теоремы непосредственно вытекает, что:

1)Решение первой краевой задачи (10.1)-(10.3) единственно.

2)Решение первой краевой задачи(10.1)-(10.3) непрерывно зависит от правых частей начального и граничного условий.

3)Если функция и, непрерывная на замыкании QT и удовлетворяющая первой краевой задаче (10.1)-(10.3), равна нулю на Г, то она тождественно равна нулю в замыкании QT .

11. Интегрирование уравнения распространения тепла в ограниченном стержне методом Фурье.

Задача о распространении тепла в теплоизолированном с боков стержне длины l приводится к нахождению решения уравнения

 

 

 

27

u

= a2

2u

(11.1)

t

x2

 

 

в области 0£ x £ l, 0£ t <+¥, удовлетворяющего начальному условию

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

(11.2)

u(0, t) = j1

(t),ü

(0

£ t < +¥).

и граничным условиям

ý

u(l, t) = j2 (t) þ

 

 

Ограничимся рассмотрением случая, когда на концах стержня поддерживается постоянная температура, т.е. когда граничные условия имеют вид:

u(0, t) = u0 = const,ü

(0 £ t < +¥). (11.3)

u(l, t) = u1 = const

ý

þ

 

Не умаляя общности можно считать, что u0 = 0, u1 = 0, ибо в противном случае этого всегда можно добиться при помощи замены искомой функции u(x, t) по формуле

 

 

 

 

 

 

 

 

v(x, t) = u(x, t) - u0 -

u1 - u0

x,

(11.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

где v - новая неизвестная функция. Действительно, так как

 

 

v

=

u

,

2v

=

2v

, то функция v удовлетворяет тому же уравнению, что и

 

t

t

x2

x2

 

 

 

 

 

 

 

 

 

 

 

функция и:

 

v

= a2

2v

. Далее из (11.4) и (11.3) следует, что

 

 

x2

 

 

 

 

 

 

t

 

 

 

v(0, t) = 0,ü

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(0

£ t < +¥).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ý

 

 

 

 

 

 

 

 

 

 

 

 

 

v(l, t) = 0 þ

 

 

 

Таким образом, достаточно найти решение уравнения (11.1), удовлетворяющее начальному условию (11.2) и граничным условиям

u(0, t) = 0,ü

(0 £ t < +¥).

u(l,t) = 0

ý

þ

 

Как и в случае волнового уравнения, будем искать решение уравнения (11.1) в виде произведения двух функций

u = X (x)T (t),

(11.5)

одна из которых зависит только от x , а другая -только от t; причем X (x) ¹ 0 и

T (t) ¹ 0 , ибо в противном случае u(x,t)º 0, что невозможно: функция иº 0 не удо-

влетворяет начальному условию (11.2), поскольку предполагается, что f(x)¹ 0.