Толстых-Уравнматем-физики
.pdf8.5 |
2 |
l |
cx l x |
|
|
|
|
nx |
|
dx udv uv |
|
vdu |
|
|
|
|
|||||||||||||||||||||
b |
|
sin |
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
n |
l |
0 |
|
l2 |
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
u |
|
|
|
|
dv |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4c |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 n 1, |
n 1,2,... |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3n3 |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
b2k 0, |
b2k 1 |
|
|
|
|
|
8c |
|
|
|
, |
k 1,2,... |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
3 2k 1 3 |
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Ответ: искомый температурный режим |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
8c |
|
|
1 |
|
|
|
|
|
|
a 2k 1 2 |
|
|
2k 1 x |
|
|
|
|||||||||||
|
|
|
|
u x,t |
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
e |
|
l |
|
|
|
sin |
. |
|
|||||||||||||||||||
|
|
|
3 |
|
|
|
|
|
|
3 |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
k 1 |
2k 1 |
|
|
|
|
|
|
|
|
l |
|
|
|||||||||||||
Задача 8.2. Найти решение уравнения теплопроводности |
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
u |
2 |
2u |
|
|
|
|
0 x 1, t 0 |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|||
при краевых условиях u x,0 x |
|
|
|
x |
, |
|
u 0,t |
0, |
u 1,t |
|
. |
||||||||||||||||||||||||||
2 |
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
Решение. Ищем решение в виде u x,t v x,t u* x , где функция
u* x удовлетворяет граничным условиям, например, |
u* 0,5x, а |
v x,t |
|||||
удовлетворяет задаче |
|
|
|
||||
|
v |
2 |
2v |
, v x,0 u x,0 u* x x2, v 0,t v 1,t 0, |
|
||
|
t |
x2 |
|
||||
|
|
|
|
|
|
||
|
|
|
|
8.8 |
|
|
|
|
|
|
|
u x,t |
0,5x bne 2 2n2t sin nx. |
|
|
n 1
Предлагаем читателям самостоятельно получить коэффициенты, интегрируя по частям:
bn |
8.5 |
1 |
x x2 |
sin nxdx udv uv |
|
vdu |
41 1 |
n |
|
|||||
2 |
|
|
||||||||||||
|
3 |
3 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
0 |
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
u |
|
|
dx |
|
|
|
|
|
|
|
61
b2k 0, |
b2k 1 |
|
8 |
, |
k 1,2,... |
||||
3 |
2k 1 3 |
||||||||
|
|
|
|
|
|
||||
Ответ: решение данной краевой задачи |
|
|
|||||||
u x,t 0,5x |
8 |
|
1 |
|
e 2 2 2k 1 2t sin 2k 1 x |
||||
|
|
||||||||
3 |
|
|
3 |
||||||
|
k 1 |
2k 1 |
|
|
|||||
(ср. задачу 8.1 при a2 2, |
c l |
1). |
|
|
Задача 8.3. Найти распределение температуры в стержне 0 x l с теплоизолированной боковой поверхностью, если один конец теплоизолирован, а второй поддерживается при нулевой температуре, в начальный момент времени температура во всех точках стержня f(x).
Решение. а) Пусть левый конец теплоизолирован, а правый поддерживается при нулевой температуре.
Необходимо решить краевую задачу
|
u |
a2 |
2u |
, |
u x,0 f x , |
u |
|
|
|
0, |
u l,t 0. |
||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
t |
|
|
|
|
x2 |
|
|
|
|
|
|
|
x |
|
x 0 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Решение ищем в виде u x,t X x T t : |
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
TX |
|
|
T |
|
|
X |
|
2 |
const 0 ; |
||||
|
|
|
|
|
|
|
|
a2T |
|
X |
|||||||||||||
|
|
|
|
|
|
T X a |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
T a2 2T T e a2 2t, |
|
||||||||||||||||
|
|
|
|
|
|
|
X 2 X 0 X Acos x Bsin x. |
||||||||||||||||
Из условия |
u |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
x |
|
|
|
следует, что X 0 0, |
|
||||||||||||||||||
|
|
|
|
x 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X A sin x B cos x B 0.
Из условия u l,t 0 следует X l 0 Acos l
cos l 0 A 0 l 2n 1)
2
n 2n 1 , n 0,1,2,...
2l
62
Итак, температурный режим стержня
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a 2n 1 2 |
|
|
|
2n 1 x |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
и х,t ane |
|
2l |
|
cos |
. |
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
n 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2l |
|
|||||||
Из начального условия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2n 1 x |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
и х,0 f x an cos |
, |
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n 0 |
|
|
|
|
2l |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 l |
|
|
|
|
|
2n 1 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
an |
|
|
f x cos |
|
|
|
|
|
|
|
dx, |
n 0,1,2,... |
|
||||||||||||||||||
|
|
l |
|
|
|
2l |
|
|
||||||||||||||||||||||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Например, при f x и0 |
|
const |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
2 l |
|
2n 1 x |
|
|
|
|
|
4u0 |
|
|
|
2n 1 x |
|
l |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
an |
|
u0 cos |
|
|
|
|
|
|
|
dx |
|
|
|
sin |
|
|
|
|
|
|
, |
|||||||||||||
l |
|
|
2l |
|
|
|
2n 1 |
|
2l |
|
|
|
||||||||||||||||||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
||||||
а |
|
|
1 n 4u0 |
, n 0,1,2,.... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
n |
2n 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Ответ: температурный режим стержня |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
4u0 |
|
|
|
n |
|
a 2n 1 2 |
|
|
2n 1 x |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
u x,t |
|
|
1 |
|
e |
2l |
|
|
|
cos |
. |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
n 0 2n 1 |
|
|
|
|
|
|
|
|
2l |
|
б) Аналогично решается задача, если правый конец теплоизолирован, а левый поддерживается при нулевой температуре:
|
|
|
|
|
u 0,t ) 0, |
u |
|
|
0. |
|
|
|
|
|
|
|
x |
|
|||||
Как и выше, T e a2 2t , |
|
|
x l |
|
|||||||
|
|
|
|||||||||
X Acos x Bsin x. |
|
||||||||||
u 0,t 0 X 0 A 0, |
|
X Bsin x, |
B 0. |
||||||||
|
u |
|
|
|
0 X l Bcos l 0, cos l 0 |
||||||
|
|
||||||||||
|
x |
|
x l |
||||||||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
63
l |
2n 1 |
n |
2n 1 |
, |
n 0,1,2,... |
||||||||||||||
|
|
|
|
|
|
|
|||||||||||||
2 |
|
|
|
|
|
|
|
2l |
|
|
|
|
|
|
|
||||
Температурный режим |
|
a 2n 1 2 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
2n 1 x |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|||||
|
u x,t bne |
2l |
|
|
|
sin |
|
|
|
|
|
. |
|
||||||
|
|
|
|
|
|
2l |
|
||||||||||||
|
|
|
|
n 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Из начального условия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
2n 1 x |
|
|
|
||||
|
|
u x,0 f x bn sin |
, |
|
|
||||||||||||||
|
|
|
|
||||||||||||||||
откуда |
|
|
n 0 |
|
|
|
|
|
|
|
|
2l |
|
|
|||||
2n 1 x |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
2 l |
|
|
|
|
|
|
|
|
|
|
|
||||||
bn |
|
f (x)sin |
|
|
|
dx, |
n 0,1,2,... |
|
|||||||||||
l |
|
2l |
|
||||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8.2. Интеграл Пуассона
Для бесконечного стержня (достаточно длинного) нужно учитывать только начальное условие, то есть необходимо решить задачу (8.1), (8.2). Решение получается с помощью интеграла Фурье.
Действительно, частное решение имеет вид (как и выше):
u x,t e a2 2t A cos x B sin x R.
Интегрируя по , получим решение в виде:
u x,t e a2 2t A cos x B sin x d .
Из начального условия (8.2)
u(x,0) f (x) A cos x B sin x d ,
|
|
|
|
||
откуда |
|
|
|
||
|
1 |
|
|
||
|
|
||||
A |
|
f z cos z dz, |
|||
2 |
|||||
|
|
||||
|
1 |
|
|
||
|
|
||||
B |
|
|
f z sin z dz. |
||
|
2 |
||||
|
|
|
(8.10)
(8.11)
64
Подставляя (8.11) в (8.10), после преобразований получим:
|
1 |
|
|
|
|
|
|
|
|
|
|
||
u x,t |
|
f z dz e a2 2t cos z x d . |
(8.12) |
||||||||||
|
|||||||||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|||
Если воспользоваться формулой |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
2 |
|
|||
|
|
1 |
|
|
|
|
e |
|
|||||
e a2 2t cos d |
|
|
|
4a2t |
, |
|
|||||||
2a |
|
|
t |
|
|||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
то решение (8.12) можно переписать в виде:
|
1 |
|
|
|
|
z x 2 |
|
|
|
u x,t |
|
|
f z e |
4a |
2t dz. |
(8.13) |
|||
|
|
|
|||||||
|
|
|
|||||||
2a |
t |
|
|
||||||
|
|
|
|
|
|
|
|
Эта формула носит название интеграла Пуассона и представляет собой решение задачи о распространении тепла в неограниченном стержне.
Установим физический смысл формулы (8.13). Рассмотрим функцию
|
|
|
|
0, |
|
|
x x |
0 |
, x x |
0 |
|
h, |
|
|
|
|||
|
|
f* x |
|
|
|
|
|
|
|
|
|
|
|
|
||||
Тогда функция |
|
|
f x , x0 x x0 h. |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
|
|
|
|
|
z x 2 |
|
|
|
1 |
|
|
|
x0 h |
|
z x 2 |
|
u* x,t |
|
|
f * z e |
|
4a2t dz |
|
|
|
f (z)e 4a2t dz |
|||||||||
|
|
|
|
|
|
|
|
|||||||||||
2a |
t |
|
|
2a |
t |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
есть решение уравнения (8.1), |
принимающее при t 0 |
значение f * x . |
Применив теорему о среднем к последнему интегралу, получим
|
|
x,t |
f ( )h |
|
|
|
x 2 |
|
||
u |
* |
|
e |
2 |
|
x0 x0 h. |
||||
|
|
|
|
4a t , |
||||||
|
|
|
|
|||||||
|
2a t |
|||||||||
|
|
|
|
|
|
|
|
|
||
Эта формула дает значение температуры в любой точке стержня в |
||||||||||
любой момент времени, если при |
t 0 |
температура u* 0 вне отрезка |
||||||||
x0; x0 h , где она равна |
f x . Если линейная плотность стержня, с – |
65
теплоемкость материала, то количество тепла в элементе x0; x0 h при t 0 будет Q f h c. Тогда значение
1 |
|
|
e |
x 2 |
|
u* x,t |
|
u* x,t c |
|
|
|
4a2t |
|
|
|||||
|
|
|
f h |
Q |
|||||
2a t |
|||||||||
|
|
|
|
|
определяет температуру в любой точке стержня в любой момент времени, если при t 0 в сечении (предельный случай при h 0) был мгновенный источник тепла Q c .
Задача 8.4. Найти распределение температуры в бесконечном стержне, если начальное распределение температуры определяется
функцией |
|
|
|
u |
|
|
, |
|
|
|
x x x |
|
, |
|||
|
|
|
|
|
|
|
|
|
|
|
||||||
u x,0 f x |
0 |
|
|
|
|
1 |
|
|
2 |
|
||||||
|
|
|
|
|
0, |
|
|
|
|
x x1, |
|
x x2. |
||||
Решение. Так как стержень бесконечный, то решение уравнения |
||||||||||||||||
теплопроводности |
u |
a2 |
2u |
|
запишется в |
виде интеграла Пуассона |
||||||||||
t |
x2 |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
(8.13). Для данной функции |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
u x,t |
u0 |
|
|
|
x |
2 |
|
z x 2 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
e 4a2t dz. |
|
||||||||||
|
|
|
|
|
|
|
||||||||||
|
|
t |
|
|||||||||||||
|
|
|
|
|
2a |
|
|
x |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
Сделаем замену в интеграле:
z x |
|
y, |
|
z |
|
|
|
|
x1 |
|
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
a 2t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
dz a |
|
|
|
y |
y1 |
|
x1 |
x |
|
|
y2 |
x2 |
x |
|
||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||
2tdy, |
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
a |
|
2t |
|
|
|
2t |
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
||||||
|
|
|
|
|
|
u x,t |
u0 |
|
y |
2e |
y2 |
dy u |
|
Ф y |
|
|
Ф y , |
|
||||||||||||
|
|
|
|
|
|
|
|
|
2 |
0 |
2 |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
2 y |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ф x |
|
1 |
|
x |
2 /2dy – функция Лапласа (табулирована). |
|
|
|
e y |
||||
|
|
|
||||
2 |
||||||
|
|
0 |
Ф 0 0, Ф x Ф x . |
|||
|
|
|
|
|
66
Замечание. Если сделать замену z x y, то
|
|
|
|
|
|
|
2a t |
||||
u x,t |
u0 |
Ф y |
|
Ф y , |
y |
i |
|
xi |
x |
|
, Ф y |
|
|
|
|
|
|||||||
2 |
|
2 |
1 |
|
|
2a t |
y
2e y2 dy.
0
Задача 8.5. По формуле Пуассона найти распределение температуры в бесконечном стержне, если начальное распределение температуры
u x,0 f x e 4x2 4x .
Решение. Решение ищем в виде интеграла Пуассона.
8.13 |
1 |
|
|
|
|
|
|
|
z x 2 |
|
|
|
|
1 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
u x,t |
|
|
|
|
e 4z2 4z e |
|
4a2t dz |
|
|
|
|
|
e A z dz , |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
2a t |
2a |
t |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
A z 4z |
2 |
|
|
|
z x 2 |
|
|
|
1 |
2 |
|
|
|
|
|
x |
x2 |
||||||||
|
4z |
|
|
|
|
4 |
|
|
z |
|
2 |
2 |
|
|
|
|
|
z |
|
. |
|||||
|
|
4a2t |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
4a2t |
|
|
|
|
|
4a2t |
4a2t |
В последнем выражении выделяем полный квадрат:
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
8a2t x |
2 |
|
|
|
|
x2 |
|
|
|
|
|
|
8a2t x 2 |
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
A z 4 |
|
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
16a |
2 |
|
1 4a |
2 |
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
4a |
|
|
|
|
|
16a |
t |
1 |
|
|
|
|
4a |
t |
t |
t |
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
|
16a2t 1 |
|
|
|
8a2t x |
2 |
|
|
16a2t 4x 4x2 |
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
. |
|
|
|
||||||||||||
|
|
|
|
4a |
t |
|
|
|
|
z |
|
|
|
1 |
|
|
|
|
|
|
|
16a |
t 1 |
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
16a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
2 |
t 1 |
|
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
16a |
|
t 4x 4x |
|
|
|
|
|
|
|
16a |
|
|
|
8a t x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
Тогда u x,t |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
t |
|
z |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
e |
|
|
16a2t 1 |
|
|
|
|
e |
|
|
|
4a |
|
|
|
16a t 1 dz. |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
2a |
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Сделаем замену |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
8a2t x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
16a2t 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2a |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
y |
, dz |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
2a t |
|
|
|
|
|
|
16a |
t 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
16a2t 1 |
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
67
|
|
|
|
|
16a2t 4x 4x2 |
|
|
|
|
|
|
|
|
|
u x,t |
1 |
|
|
|
|
|
|
2a t |
|
|
||||
|
|
e |
16a2t 1 |
|
e y2 |
|
|
|
dy, |
|||||
|
|
|
|
|
|
|||||||||
|
|
|
||||||||||||
|
2a t |
|
|
|
|
16a2t 1 |
Учитывая, что e y2 dy , получаем распределение температуры
u x,t |
1 |
|
|
16a2t 4x 4x2 |
|
||
|
|
e |
16a2t 1 |
. |
|||
|
|
|
|
||||
|
16at |
2 1 |
|
Задачи для самостоятельного решения
1. Методом Фурье найти решение уравнения теплопроводности
u a2 2u при следующих условиях.
t x2
а) |
|
|
|
a 1, |
|
|
|
|
u(x,0) |
x, |
0 x 3/2, |
||
|
|
|
|
u 0,t u 3,t 0 , |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
3 x, 3/2 x 3. |
|||||
б) |
|
|
|
a 1, |
|
|
|
|
u x,0 |
x, |
0 x l / 2, |
||
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
u 0,t u l,t 0 , |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
l x, l/ 2 x l. |
|||||
в) |
u |
|
|
x 0 |
0, |
u l,t u0, |
|
|
u x,0 f x . |
||||
|
|
|
|
||||||||||
|
|||||||||||||
|
x |
|
|
v x,t . |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||||||
Указание. u u0 |
|
|
|
|
|
||||||||
г) a 1, |
u 0,t 0, |
u |
|
|
0, |
u x,0 5x. |
|||||||
|
|
||||||||||||
x |
|||||||||||||
|
|
|
|
|
|
|
|
x 2 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||
Указание. Воспользоваться задачей 8.3. |
|
||||||||||||
д) |
a 2, |
u 0,t 0, |
u(1,t) 1, |
|
u x,0 x2. |
68
2. Найти |
распределение u x,t температуры в |
однородном стержне |
0 x l |
с теплоизолированной поверхностью, |
если на концах его |
поддерживается постоянная температура: u1 – на левом конце, u2 – на правом конце, а начальное распределение температуры имеет вид
u x,0 Аx, A const.
3.Начальное распределение температуры бесконечного алюминиевого стержня имеет вид:
u |
|
50 , |
x 0,1м х 0,2м х |
|
, |
u x,0 |
0 |
|
1 |
2 |
|
0, |
|
х 0,1 ; 0,2 . |
|
|
|
|
|
|
|
|
|
Найти температуру стержня в точке x 0 через час, если
k196ккал/м ч град, 2700 кг/м3 ,
с0,214ккал/кг град, a2 R 0,34 м2/ч.
c
4. Найти закон распределения температуры внутри бесконечного однородного стержня, если в начальный момент времени температура на
участке l x l была 1000 , а в остальных точках 0 .
Указание. В задачах 3 и 4 использовать задачу 8.4.
5. На концах стержня поддерживается постоянная температура
u 0,t 0, u l,t 100 .
Найти распределение температуры по стержню в произвольный момент времени, если в начальный момент
|
200x |
, |
0 x |
l |
, |
|||
|
|
|
|
|||||
|
|
|
||||||
|
l |
|
|
|
|
2 |
|
|
u x,0 |
|
|
l |
|
|
|
||
100, |
|
|
x l. |
|
||||
|
|
|
|
|||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
69
9. Уравнение Лапласа и задачи, к нему приводящие
Уравнение Лапласа имеет вид:
u |
2u |
|
2u |
|
2u |
0. |
(9.1) |
||
x2 |
y2 |
z |
2 |
||||||
|
|
|
|
|
Функция u, дважды непрерывно дифференцируемая в некоторой области V, удовлетворяющая уравнению Лапласа, называется
гармонической в этой области.
Рассмотрим некоторые задачи, приводящие к уравнению Лапласа.
9.1. Потенциал стационарного электрического поля
Пусть в однородной среде объема V проходит электрический ток,
плотность которого задана вектор-функцией J x,y,z Jxi Jy j Jzk.
Предположим, что плотность тока не зависит от времени t и в области V
нет источников тока. Тогда поток вектора |
J через любую замкнутую |
||||||
поверхность S V равен нулю: |
|
|
|
|
|
то есть из формулы Гуасса- |
|
|
J |
nds, |
|||||
S |
|
||||||
Остроградского дивергенция |
|
||||||
|
|
0. |
(9.2) |
||||
divJ |
На основании обобщенного закона Ома электрическая сила E J / или J E , где проводимость среды (считаем ее постоянной). Из общих уравнений электрического поля следует, что если процесс
стационарный, то поле E безвихревое, то есть ротор |
rotE 0. Это |
||||||||||||||
означает, что электрическая сила |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
grad |
|
|
|
|
|
|
|
|
|
|
|
, |
(9.3) |
|
E |
i |
j |
k |
|||||||||||
|
x |
y |
|
||||||||||||
|
|
|
|
|
|
|
|
|
z |
|
где потенциал поля.
Окончательно с учетом (9.2), (9.3) |
div grad 0, то есть потенциал |
|||||||||
стационарного электрического поля удовлетворяет уравнению: |
|
|||||||||
|
2 |
|
2 |
|
2 |
|
0. |
(9.4) |
||
x2 |
y |
2 |
z2 |
|||||||
|
|
|
|
|
70