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

УМК

.PDF
Скачиваний:
91
Добавлен:
17.03.2015
Размер:
3.93 Mб
Скачать

преобразуется в уравнение

2 U

 

= 0 , которое имеет общее решение

∂ξ ∂η

 

 

U(ξ, η) = Φ(ξ) + F(η),

где Φ и F - произвольные дважды дифференцируемые функции (см. пример 2.40). Для определения функций Φ и F, т.е. для определения закона колебаний струны, требуется использовать начальные условия, а для некоторых задач и граничные. Если вернуться к старым переменным x и t , то общее решение примет вид

U(x, t) = Φ(x + at) + F(x at).

Здесь F(x at) характеризует прямую волну (кривая F(x)смещается вправо со скоростью a ), а Φ(x + at) - обратную волну (кривая Φ(x) смещает-

ся влево со скоростью a ).

 

 

 

 

 

Если рассматривается задача

Коши

для бесконечной

струны

− ∞ < x < ∞, то по заданным начальным условиям

 

U(x,0) = ϕ(x) , Ut (x,0) = ψ(x)

(2.74)

определяются функции Φ и F, и искомое решение имеет вид

 

U(x, t) = ϕ(x at) + ϕ(x + at) +

1

 

x+a t

 

ψ (z) dz .

(2.75)

 

2

 

2a

x a t

 

Формула (2.75) называется формулой Даламбера. Эта формула доказывает единственность решения задачи Коши.

В частности, когда начальная скорость равна нулю ( ψ (x) = 0), то

U(x, t) = ϕ(x at) + ϕ(x + at) , 2

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

ная форма каждой волны определяется функцией ϕ(x)2, равной половине на-

чального отклонения.

В случае полубесконечной струны, кроме начальных условий (2.74), заданных при 0 x < ∞, необходимо добавить еще граничное условие (конец предполагается в точке x = 0 )

U (0, t) = 0

(2.76)

для закрепленной в точке x = 0 струны,

∂ U (0, t) = 0

(2.77)

x

 

для свободного конца в точке x = 0 ,

 

∂ U (0, t) h U (0, t) = 0 .

 

x

 

для упругого закрепления в точке x = 0 .

В случае однородных граничных условий (2.76) или (2.77) решение задачи о колебании полубесконечной струны сводится к решению задачи о колебании бесконечной струны путем продолжения начальных условий на всю ось

нечетным образом

для условия (2.76), т.е. полагают ϕ(x) = −ϕ(x) ,

ψ(x) = −ψ(x) ,

и четным образом для условия (2.77), т.е. ϕ(x) (x) ,

ψ(x) (x) .

 

 

Примеры решения задач

ПРИМЕР 2.43. Найти форму достаточно длинной струны, определяемой

уравнением

2 U = a 2 2 U , в момент времени t = π 2 , t = π, если заданы

 

 

t2

x2

начальные смещения и скорости:

а)

U(x,0) = sin x,

Ut (x,0) = cos x ;

б)

U(x,0) = 0, Ut (x,0) = cos x ;

в)

U(x,0) = sin x,

Ut (x,0) = 0 .

Решение. По постановке вопроса надо найти решение U (x, t) задачи

Коши (2.73), (2.74) в области: − ∞ < x < ∞, t > 0 . Оно определяется формулой Даламбера (2.75).

СЛУЧАЙ а). Полагая в формуле Даламбера ϕ(x) = sin x , ψ(x) = cos x , найдем смещение U (x, t) в любой точке и любой момент t :

 

1

[sin(x t) + sin(x +t)]+

1

x +t

U(x, t) =

cos z dz =

 

 

2

 

2 xt

= 1 [sin(x t) + sin(x +t)]+ 1 [sin(x +t) sin(x t)] =sin(x +t).

2

2

Откуда определяем форму кривой в указанные моменты времени:

U (x, π2) = cos x , U (x, π) = −sin x .

Кривые изображены на рис. 2.3.

U(x, y)

t = π

t = π 2

 

t = 0

 

 

 

 

 

 

 

 

 

π

− π 2 0

π 2

π

3π 2

2π

x

Рис. 2.3

СЛУЧАЙ б) Начальные смещения струны равны нулю, т.к. U(x,0) = 0 . При ϕ(x) = 0 колебательный процесс будет описан по формуле

 

1

 

 

x +t

 

1

 

 

 

xx+tt

 

1

[sin(x +t)sin(x t)]=

U(x, t) =

 

cos z dz =

sin z

 

=

 

 

 

 

 

2

 

 

x t

2

 

 

 

 

2

 

 

=

1

 

2 cos

x + t + x t

sin

x + t x + t

= cos x sin t .

 

 

 

 

2

 

 

2

 

 

 

 

 

 

2

 

В момент времени t = π2 струна имеет форму косинусоиды: U t 2 = cos x ,

а в момент t = π она совпадает с осью абсцисс: U t = 0 .

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

ψ(x)= 0 . Тогда имеем

 

 

 

 

 

 

 

U(x, t) =

sin (x t) + sin (x +t)

=

1

2 sin

x t + x + t

cos

x t x t

=

 

 

 

 

2

2

2

2

 

= sin x cos t .

 

 

 

 

 

 

 

Форма струны в указанные моменты времени определяется уравнениями:

U (x, π2) = 0 , U (x, π) = −sin x .

2.21 РЕШЕНИЕ НЕКОТОРЫХ КРАЕВЫХ ЗАДАЧ МЕТОДОМ ФУРЬЕ. РЕШЕНИЕ УРАВНЕНИЙ КОЛЕБАНИЙ МЕТОДОМ ФУРЬЕ

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

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

задачами ШтурмаЛиувилля.

Искомое решение представляется рядом по про-

изведениям

собственных функций задач

ШтурмаЛиувилля. Приведем

с х е м у э т о г о м е т о д а

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

параболического типов – волнового уравнения и уравнения теплопроводности. Рассмотрим однородные уравнения

 

 

1 2 U

=

2

U

(2.78)

 

 

a 2

t 2

x 2

 

 

 

 

 

 

1 2 U

=

2

U

 

 

 

 

 

t

 

,

(2.79)

 

a 2

 

 

 

x 2

 

для которых граничные условия имеют вид

 

U(0, t)= 0,

U(l, t)= 0 .

(2.80)

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

 

 

 

 

 

U(x, 0)= ϕ(x),

U t (x, 0)= ψ(x)

(2.81)

для (2.78)

 

 

 

 

 

 

 

U(x, 0)= ϕ(x)

(2.82)

для (2.79).

Процесс решения разбивается на два этапа: I – нахождение частных решений; II – нахождение общего решения и удовлетворение начальным условиям.

I. Ищутся всевозможные частные решения (2.78) в виде

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

(2.83)

В результате подстановки функции такого вида в уравнение (2.78) полу-

чаем

1

X(x)T ′′(t)= X ′′(x)T(t) или

T′′(t)

 

=

X′′(x)

 

= −λ (const).

a 2

a 2 T(t)

X(x)

 

 

 

Последнее уравнение распадается на два обыкновенных дифференциальных уравнения относительно функций X(t) и T(t):

T′′(t)+ a 2 λ T(t)= 0 .

(2.84)

X′′(x)+ λ X(x)= 0 .

(2.85)

Как говорят, в уравнении (2.78) переменные разделены.

 

Подставляя (2.83) в граничные условия (2.80), получим

 

X(0)T(t)= 0, X(l)T(t)= 0 , откуда

 

X(0)= 0, X(l)= 0

(2.86)

(т.к. случай T(t)= 0 не представляет интереса, поскольку тогда U 0).

Для определения функции X(x) получена з а д а ч а Ш т у р м а -

Л и у в и л л я: найти такие значения параметра λ, при которых существуют нетривиальные решения уравнения (2.85), удовлетворяющие граничным условиям (2.86). Те значения параметра λ, для которых задача (2.85) – (2.86) имеет нетривиальные решения, называются собственными значениями, а сами решения – собственными функциями. Нетривиальные решения задачи (2.85) – (2.86) возможны лишь при значениях

 

k π

2

λk

=

 

 

, k = 1,2,3,K

 

 

 

l

 

Этим значениям λk соответствуют собственные функции

 

 

 

 

X k

(t)= sin

k π x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

При λ = λk

общее решение уравнения (2.84) имеет вид

 

T (t)= a

 

cos

k πa t

+ b

 

sin

k πa t

,

 

k

 

k

 

 

k

 

l

 

 

 

l

 

 

 

 

 

 

 

 

 

 

где a k и bk

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

 

Таким образом, произведения функций X k (t) и Tk (t) образуют частные решения уравнения (2.78), удовлетворяющие краевым условиям (2.80)

 

 

 

kπa

 

kπa

 

kπ

U k

(x, t)= a k

cos

 

t + bk

sin

 

t sin

 

x, k = 1,2,K.

 

 

 

 

 

 

l

 

l

 

l

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

U(x, t) = U k (x, t). (2.87)

k=1

Для уравнения (2.78) общее решение имеет вид

 

 

 

 

 

 

 

kπa

 

 

kπa

 

kπ

 

U (x, t) =

a k

cos

 

 

t + bk

sin

 

t sin

 

 

x .

(2.88)

 

 

 

 

 

 

 

 

k=1

 

 

 

l

 

 

 

 

 

l

 

 

l

 

Подставляя решение

(2.88)

в начальные условия (2.81), определяют зна-

чения коэффициентов a k

и bk , пользуясь разложением функций ϕ(x) и ψ(x)

в ряд Фурье (в ряд по системе собственных функций {X n (t)}).

В результате

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

l

 

kπx

 

 

 

 

 

2

 

l

 

 

kπx

 

 

a k =

 

ϕ(x)sin

dx ,

 

bk =

 

 

ψ(x)sin

dx .

(2.89)

 

 

 

kπa

 

 

l

0

 

 

l

 

 

 

 

 

0

 

 

l

 

Решение

уравнения

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

(2.79) получается

применением

этой же схемы с той лишь разницей, что вместо уравнения (2.84) надо решать уравнение

 

 

1

 

T(t)

= −λ

 

 

или

T′ + a 2 λ T = 0 ,

 

 

a 2

 

T(t)

 

 

 

общее решение которого есть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k π a

2

 

 

T

 

(t) = a

 

 

 

t

 

 

 

 

 

 

 

 

k

e

 

l

, k = 0,1,K

 

 

 

k

 

 

 

 

 

 

 

 

 

 

где a k произвольные постоянные.

 

kπx

 

 

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

 

(x, t) = T (t)sin

 

и общее решение (2.87) принимает

k

 

 

 

 

 

k

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вид

U(x, t) = a k

k=1

 

kπa

2

sin

kπx

 

exp

 

 

t

 

.

(2.90)

 

 

 

 

l

 

 

l

 

 

 

 

 

 

 

 

 

 

Потребовав выполнения начального условия (2.82), коэффициенты a k будут найдены как коэффициенты разложения функции ϕ(x) в ряд Фурье:

a k

=

2

l

ϕ(x)sin

kπx

dx .

(2.91)

 

 

 

 

l

0

 

l

 

Приведенные здесь решения для линейных однородных уравнений с од-

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

нии неоднородных и однородных уравнений с однородными и неоднородными

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

 

 

Примеры решения задач

 

ПРИМЕР 2.44. Найти закон колебания однородной струны, закрепленной

на концах x = 0 и x = l. В начальный момент струна оттянута в точке x 0 на

высоту h (рис. 2.4) и затем отпущена без начальной скорости.

 

U

 

 

 

 

A

 

 

 

h

B

 

 

 

 

0

x 0

 

x

 

 

l

 

 

 

Рис. 2.4

 

Решение. Задача сводится к решению уравнения

 

U tt = a 2

U xx

 

 

при граничных условиях U(0, t) = U(0, x) = 0 и начальных условиях

 

 

 

h

 

 

 

,

 

 

 

 

 

x 0

U(0, x) =

h (l − x)

 

 

 

,

 

 

 

l

x 0

 

 

0 < x x 0

 

U t (x,0) = 0,

0 < x ≤ l ,

x 0 < x ≤ l

 

где

h

x = U - уравнение прямой ОА,

h (l − x)

= U – уравнение прямой

 

 

 

x 0

 

l − x 0

АВ (оба записываются как уравнение прямой с угловым коэффициентом). Решение поставленной задачи определяется рядом (2.88)

 

aπ

 

aπ

 

 

kπ

U(x, t) = a k

Cos

 

t + bk

Sin

 

t

Sin

 

x ,

 

 

 

k =1

 

l

 

l

 

 

l

где a k , bk - коэффициенты Фурье для функций, которые вычисляются по формулам (2.89).

a k

=

2

l ϕ(x) Sin

kπ

 

x dx ,

где ϕ(x) = U(x,0) .

 

 

 

 

 

 

 

 

l

0

 

 

 

l

 

 

 

 

 

 

 

В нашем случае

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xo

 

 

 

 

 

 

 

 

l

 

 

a k

=

2

 

h

x Sin

kπ

x dx +

l

h( x)

Sin

kπ

x dx .

 

x o

 

 

 

 

 

l

 

0

 

 

 

l

 

xo

l - x o

l

 

Вычислим первый интеграл

 

h

 

xo

x Sin

kπ

x dx =

 

x = u,

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

o

0

 

 

 

 

l

 

dx = du,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

xl

 

 

 

kπ

 

xo

 

l xo

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

Cos

 

x

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

kπ xo

 

 

 

x o

kπ

 

 

 

l

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kπ

 

 

Sin

 

x dx = dv

 

 

 

 

 

 

l

 

 

 

 

=

 

 

 

 

 

 

 

 

v = −

 

l

Cos

kπ

x

 

kπ

 

 

 

 

 

 

l

 

 

 

kπ

 

 

=

 

 

Cos

x dx

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hl

 

 

 

kπ

 

 

 

 

 

 

 

hl

 

 

 

l

 

 

 

 

 

kπ

 

 

xo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

Cos

x 0

+

 

 

 

 

 

 

Sin

x

0

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kπ

 

 

 

l

 

 

 

 

x o kπ

 

 

kπ

 

 

 

l

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

 

hl

Cos

kπ

+

 

 

 

 

hl2

 

 

 

Sin

kπ

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

x

0 .

 

 

 

 

 

kπ

 

 

xo (kπ)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

Вычислим второй интеграл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

h(l − x)

Sin

kπx

dx =

 

 

 

 

h

 

l (l − x) Sin

kπx

dx =

 

 

 

 

 

 

 

 

 

 

 

 

xo

 

 

l - x o

 

 

 

 

l

 

 

 

 

 

 

l - x o xo

 

 

 

 

l

 

l − x = u,

 

Sin

kπx

dx = dv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

l

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kπx

 

 

 

 

 

 

 

 

 

 

du = − dx,

 

v = −

 

 

l

 

Cos

 

 

 

 

 

 

 

 

 

 

 

kπ

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

h

l (l − x)

Cos

kπx

 

 

 

 

 

 

 

 

 

l - x o kπ

 

l

 

 

 

 

 

l

 

l l

kπx

 

 

 

Cos

 

dx =

 

 

 

x

o

kπ xo

l

 

 

 

 

 

 

=

hl

Cos

kπx 0

hl

 

l

Sin

kπx

kπ

 

 

 

 

 

 

l

kπ (l − x o ) kπ

 

l

l

=

xo

 

 

 

 

 

=

hl

 

Cos

kπx

0

+

 

 

 

hl2

 

 

 

 

 

 

Sin

kπx

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kπ

 

 

 

 

l

 

 

 

(

 

 

)2

(

l

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kπ

 

 

 

xo )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

hl

Cos

kπx

o

 

 

 

hl2

 

 

Sin

kπx

o

 

 

 

 

 

 

 

 

 

 

 

a k

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

x o (kπ)2

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

hl

 

 

Cos

 

kπx

o

+

 

 

 

 

 

 

 

hl

2

 

 

 

 

 

 

 

Sin

kπx

o

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kπ

 

 

 

 

 

 

 

 

 

 

 

 

(kπ)2 (l − x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

o

)

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 hl2

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

kπx o

 

 

 

 

 

 

 

 

 

 

 

 

2 hl2

 

 

 

 

 

kπx o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

2 Sin

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sin

 

 

 

 

 

 

x o (l − x o ) (kπ)

 

 

 

l (kπ)

 

 

 

 

x o

 

 

l − x o

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

l

 

Определим bk .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bk =

 

 

2

 

 

 

 

 

l

ψ(x) Sin

kπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ(x) =

∂U(x,0) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x dx ,

где

 

 

 

 

 

 

 

 

 

 

a kπ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

В нашем случае

ψ(x) = 0 ,

тогда bk

 

= 0 и решение

 

U(x, t)

имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 hl

2

 

 

 

 

 

 

 

 

 

1

 

 

Sin

kπx o

Cos

aπ

t Sin

kπ

 

 

 

U(x, t) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x .

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

l

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

π x o (l − x o )

k=1 k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученный ряд описывает колебательный процесс: смещение

U(x, t)

точки x струны в любой момент времени t . Чтобы определить форму струны в момент t 0 , надо протабулировать функцию U(x, t 0 ), ограничившись несколькими.

2.22УРАВНЕНИЕ ТЕПЛОПРОВОДНОСТИ. ПОСТАНОВКА КРАЕВЫХ ЗАДАЧ

Вп. 2.18 отмечалось, что нестационарные процессы теплопроводности, диффузии и другие описываются уравнениями параболического типа (2.63) – (2.66). Ограничимся рассмотрением процесса распространения тепла в одномерном теле – стержне, достаточно тонком, чтобы в любой момент времени температуру U во всех точках поперечного сечения можно считать одинаковой.

U = a 2 2 U + f (x, t).

(2.92)

t

x 2

 

Примеры решения задач

ПРИМЕР 2.45. Поставить задачу о распределении температуры в тонком однородном стержне 0 x ≤ l с теплоизолированной боковой поверхностью. Внутри стрежня происходит свободный теплообмен, т.е. отсутствуют источники и поглотители тепла. Начальная температура есть некоторая функция от x , а на концах поддерживается следующий температурный режим:

а) на концах задается температура a′ − постоянная, a′′ − переменная;

б) на концы стержня подаются извне заданные тепловые потоки; в) концы стрежня теплоизолированы;

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

Решение. Задача о распределении температуры в однородном стержне без

тепловых

источников

сводится к решению уравнения (2.92), в котором

f (x, t) 0 , т.е.

 

 

 

2 U = a 2 2 U , o < x < l, t > 0,

 

 

 

 

t

x 2

где a =

k

 

коэффициент температуропроводности ( k коэффициент теп-

 

 

 

c ρ

 

лопроводности материала стержня, c удельная теплоемкость, ρ − плотность

стержня). Начальное условие запишется в виде

U(x,0) = ϕ(x), 0 x ≤ l,

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

СЛУЧАЙ

а) поскольку концы стержня ( x = 0,

x = l ) поддерживаются

при заданных постоянных температурах, задача

aформулируется

математически следующим образом: найти решение уравнения

U = a 2 2 U ,

0 < x < l , t > 0

 

t

x 2

 

 

при начальном условии

 

U(x,0) = ϕ(x) ,

0 x ≤ l ,

 

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

 

U(0, t) = u1 ,

U(l, t) = u 2 , t > 0, где u1

= u 2 = const .

В случае a′′, когда концы стержня находятся при заданных переменных температурах, формулировка задачи следующая: найти решение уравнения

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