傅里叶热传导定律
属土的字-搞笑喜剧
2023年2月21日发(作者：南仁东简介)第八章热传导方程的傅里叶解小结及习题答案
1
第八章热传导方程的傅里叶解
第一节热传导方程和扩散方程的建立
8.1.1热传导方程的建立
推导热传导方程和前面弦振动所用的数学方法完全相用，不同之处在于具体的物理规律
不同。这里用到的是热学方面的两个基本规律，即能量守恒和热传导的傅里叶实验定律。
热传导的傅里叶实验定律：设有一块连续的介质，选定一定的坐标系，并用(,,,)uxyzt
表示介质内空间坐标为的一点在t时刻的温度。若沿x方向有一定的温度差，在x方向也就
一定有热量的传递。从宏观上看，单位时间内通过垂直x方向的单位面积的热量q与温度的
沿x方向的空间变化率成正比，即
x
u
qk
x
（8-1.1）
q称为热流密度，k称为导热系数。公式中的负号表示热流的方向和温度变化的方向正好相
反，即热量由高温流向低温。
研究三维各向同性介质中的热传导，在介质中三个方向上存在温度差，则有
x
u
qk
x
，
y
u
qk
y
，
z
u
qk
z
或
qku
即热流密度矢量q与温度梯度u成正比。
下面以一维均匀细杆为例，根据傅里叶实验定律和能量守恒定律推导介质中的热传导方
程。
第一步，定变量。研究介质x位置处在t时刻的温度(,)uxt。
第二步，取局部。在介质内部隔离出从x到xx一段微元长度，在t到tt时间
内温度的变化(,)(,)uuxttuxt。
x
xx
x
1
Q
2
Q
第三步，立假设。假设均匀介质的横截面积为A，质量密度为，比热为c，热传导系
数为k。
第八章热传导方程的傅里叶解小结及习题答案
2
第四步，找规律。隔离出来的微元长度在t到tt时间内吸收的热量为：
QcmucAxu（8-1.2）
在t到tt时间内，同过x位置处的横截面的热量为：
1xx
x
QqAtkuAt（8-1.3）
在t到tt时间内，同过xx位置处的横截面的热量为：
2xxx
xx
QqAtkuAt
（8-1.4）
如果在微元段内有其他的热源，假设在单位时间单位体积内产生的热量为(,)Fxt，则该热
源在微元内产生的热量为：
(,)
3
QFxttAx（8-1.5）
第五步，列方程。根据能量守恒定律，净流入的热量应该等于介质在此时间内温度升高
所需要的热量。
123
QQQQ
即
(,)
xx
xxx
cAxukuAtkuAtFxttAx
(,)xx
xxx
uu
u
ckFxt
tx
得到：
(,)
txx
kFxt
uu
cc
令
k
a
c
，
(,)
(,)
Fxt
fxt
c
则得到热传导方程为
(,)2
txx
uaufxt
（8-1.6）
当介质内部无其他热源时，热传导方程是齐次的，为
2
txx
uau（8-1.7）
8.1.2扩散方程的建立
第八章热传导方程的傅里叶解小结及习题答案
3
扩散问题研究的是杂质在其他介质中的浓度分布，得到的扩散方程与热传导方程有完全
一样的形式。过程略。
8.1.3热传导问题的定解条件
与弦的振动一样，其定解条件包括边界条件和初始条件。
初始条件为：已知初始时刻细杆上各点的温度分布(,)0ux
其边界条件有三种：
第一边界条件：已知细杆端点的温度(,)0ut或者(,)ult。
第二边界条件：已知通过端点的热量，即已知端点的
x
u。例如：当介质x=0端和外界
绝热，此时
(,)00
x
ut。
第三边界条件：例如，已知端点x=l与某种介质按热传导中的牛顿实验定律进行着热量
交换，已知端点的温度为(,)ult，与其接触的介质的温度为()
1
t，有牛顿实验定律知道：
在单位时间内由端点x=l流入介质的热量为
[(,)()]
1
QhulttA
由傅里叶实验定律可知，在单位时间内，端点x=l流出热量为：
(,)
x
QkultA
由QQ
，就可以得出第三边界条件为
(,)(,)()()
1x
kulthulthtt
其中，k为热传导系数，h为热交换系数。
第二节混合问题的傅里叶解
8.2.1混合问题的解
对于有界杆的热传导问题，我们先考虑齐次方程和齐次边界条件下的混合问题。即：
2
0
0
(0,0),(82.1)
0,0(0),(82.2)
(0).(82.3)
txx
xxl
t
uauxlt
uut
uxxl
第一步,分离变量，将二阶偏微分方程转化为两个常微分方程。
令
第八章热传导方程的傅里叶解小结及习题答案
4
)()(),(tTxXtxu
将此代入泛定方程（8-2.1），得到两个常微分方程：
()()20TtaTt
（8-2.4）
0)()(
xXxX（8-2.5）
第二步，将(,)uxt原来的边界条件转化为()Xx的边界条件。
将此(,)()()uxtXxTt代入边界条件，得()Xx的边界条件：
0)0(X，0)(lX（8-2.6）
第三步,求解本征值问题
通过讨论分析得出只有0时，方程（8-2.5）的解才有意义。因此，0时解（8-2.5）
式得
()cossinXxAxBx.
将这个通解代入边界条件（8-2.6），就有
0;
cossin0.
A
AlBl
即
0;
sin0.
A
Bl
于是
0sinl，即nl,3,2,1n.
得到本征值：
2
l
n
n
,3,2,1n
相应的本征函数是：
x
l
n
xX
n
sin)(
第四步,求特解，并进一步叠加出一般解：
对于每一个本征值
n
，解（8-2.5）式得出相应的)(tT
n
:
2()()
na
t
l
nn
TtCe
.
得到了满足偏微分方程和边界条件的特解:
第八章热传导方程的傅里叶解小结及习题答案
5
2()(,)sin
na
t
l
nn
n
uxtCex
l
,3,2,1n.
得到方程的一般解为
2()
1
(,)sin
na
t
l
n
n
n
uxtCex
l
（8-2.7）
第五步,利用本征函数的正交性确定叠加系数：
现在根据初始条件中的已知函数)(x定出叠加系数
n
C，将上面的一般解代入初始条
件，并利用本征函数sin
m
x
l
的正交性得到系数为
0
2
()sindl
n
nx
Cxx
ll
（8-2.8）
公式（8-2.7）给出了均匀细杆上温度场的分布，表明温度场随时间做指数衰减。
第三节初值问题的傅里叶解
8.3.1利用傅里叶积分求出热传导的初值问题
对于无穷长一维介质上的热传导问题，可以表示为
2
0
(,0),(83.1)
().(83.2)
txx
t
uauxt
uxx
解：令
)()(),(tTxXtxu
代入泛定方程（8-3.1），得到两个常微分方程：
()()20TtaTt
（8-3.3）
0)()(
xXxX（8-3.4）
解式（8-3.3）得到：
2()atTtCe（8-3.5）
由公式（8-3.5）可以看出：当0时，温度随时间的变化将趋于无穷大，这与物理事实
不符，因此，0，令2。（8-3.3）和（8-3.4）的解为与有关系的一系列解，记
为
22()atTte
（8-3.6）
解式（8-3.4）得到：
第八章热传导方程的傅里叶解小结及习题答案
6
()()cos()sinXxAxBx
于是得到热传导的一系列解为
22(,)[()cos()sin]atuxteAxBx
（8-3.7）
由于这里的没有边界条件的限制，所以为任意实数值。则(,)uxt的一般解为公式（8-3.7）
对所有值对应解的叠加，由于为连续实数，因此，(,)uxt的一般解为公式（8-3.7）对
从到进行积分。即
22(,)[()cos()sin]datuxteAxBx
（8-3.8）
把初始条件代入上式得到：
()[()cos()sin]dxAxBx
（8-3.9）
其中傅里叶系数：
()()cos
1
d
2
A
（8-3.10）
()()sin
1
d
2
B
（8-3.11）
把公式（8-3.10）与（8-3.11）带入公式（8.3-9）得到：
221
(,)()[cos()d]d
2
atuxtex
（8-3.12）
利用2dπxex
，得出
2
22
2
()
4
11
cos()d
2
2
x
at
atexe
at
因此，(,)uxt可以写为
2
2
()
4
1
(,)()d
2
x
atuxte
at
（8-3.12）
8.3.2热传导傅里叶解的物理意义
细杆上位置的点热源在整个细杆上引起的温度分布为：
2
2
()
4
1
(,)()d
2
x
atuxte
at
解（8-3.12）式可以看作是由各个瞬时点热源引起的温度分布的叠加。
第八章热传导方程的傅里叶解小结及习题答案
7
第四节一端有界的热传导问题
8.4.1左端有界热传导定解问题的解
2
0
0
(0,0),(84.1)
0(0),(84.2)
(0).(84.3)
txx
x
t
uauxt
ut
uxx
方法1：直接用分离变量法求解。
解：令
)()(),(tTxXtxu
将此代入泛定方程（8-4.1），得到两个常微分方程：
()()20TtaTt
（8-4.4）
0)()(
xXxX（8-4.5）
将此代入边界条件（8-4.2），得到：
()00X（8-4.6）
解式（8-4.4）得到：
2()atTtCe（8-4.7）
由公式（8-4.7）可以看出：当0时，温度随时间的变化将趋于无穷大，这与物理事实
不符，因此，0，令2。（8-4.4）和（8-4.5）的解为与有关系的一系列解，记
为
22()atTte
（8-4.8）
解式（8-4.5）得到：
()()cos()sinXxAxBx
把边界条件（8-4.6）代入上式得到：()0A，因此
()()sinXxBx
于是得到热传导的一系列解为
22(,)[()cos()sin]atuxteAxBx
（8-4.9）
第八章热传导方程的傅里叶解小结及习题答案
8
由于这里的没有边界条件的限制，所以为任意实数值。则(,)uxt的一般解为公式（8-4.9）
对所有值对应解的叠加，由于为连续实数，因此，(,)uxt的一般解为公式（8-4.9）对
从到进行积分。即
22(,)()sindatuxtBex
（8-4.10）
把初始条件代入上式得到：
()()sindxBx
（8-4.11）
得出：
()()sin
0
1
dB
（8-4.12）
把公式（8-4.12）带入公式（8-4.10）得到：
22
0
1
(,)(){[cos()cos()]d}d
2
atuxtexx
（8-4.13）
利用2dπxex
，得出
2
22
2
()
4
11
cos()d
2
2
x
at
atexe
at
2
22
2
()
4
11
cos()d
2
2
x
at
atexe
at
因此，(,)uxt可以写为
22
22
()()
44
0
1
(,)()[-]d
2
xx
atatuxtee
at
（8-4.14）
第八章热传导方程的傅里叶解小结及习题答案
9
方法2：把半无界拓展为无界
如何拓展？
先看无界热传导问题在坐标原点的温度分布具有什么样的特点。
由第三节可知，无界热传导问题的解为
2
2
()
4
1
(,)()d
2
x
atuxte
at
在0x点，有：
2
24
1
(0,)()d
2
atute
at
2
24
2
1
(0,)()d
2
2
at
x
ute
at
at
当()x为奇函数时，(0,)0ut满足第一类齐次边界条件。
当()x为偶函数时，
(0,)0
x
ut满足第二类齐次边界条件。
所以：
（1）当半边界为第一类齐次边界条件时，把半无限问题扩展为无限问题为：
2
0
0
0
(0,0)
0(0)
(0)
(0)
txx
x
t
t
uauxt
ut
uxx
uxx
则其解为
22
22
()()
0
44
0
11
(,)()d()d
22
xx
atatuxtee
atat
把第二项积分变量和区间变为0-
，则
22
22
()()
44
00
11
(,)()d()d
22
xx
atatuxtee
atat
（2）当半边界为第二类齐次边界条件时，把半无限问题扩展为无限问题为：
第八章热传导方程的傅里叶解小结及习题答案
10
2
0
0
0
(0,0)
0(0)
(0)
(0)
txx
x
x
t
t
uauxt
ut
uxx
uxx
则其解为
22
22
()()
0
44
0
11
(,)()d()d
22
xx
atatuxtee
atat
把第二项积分变量和区间变为0-
，则
22
22
()()
44
00
11
(,)()d()d
22
xx
atatuxtee
atat
第八章热传导方程的傅里叶解小结及习题答案
11
非齐次偏微分方程的求解
齐次偏微分方程和齐次边界条件在分离法中起着关键的作用：因为方程和边界条件是齐
次的，分离变量法才得以实现。如果定解问题中的方程和边界条件不是齐次的，还有没有可
能应用分离变量法呢？
在第七章弦的振动问题中针对非齐次边界条件先要进行齐次化处理，才能用分离变量法
已经进行了分析说明。
对于非齐次方程的解法在这里详加分析说明。
例如：强迫振动的定解问题：
2
12
0
00
(,)(0,0),(0.1)
(),()(0),(0.2)
,()(0).(0.3)
ttxx
xxl
t
tt
uaufxtxlt
ututt
uxuxxl
该弦的振动位移可以认为是由三部分干扰引起的：第一部分是由初始位移()x和初始
速度()x引起的振动；第二部分是由边界条件
(),()
12
tt干扰引起的振动；第三部分是由
强迫力(,)fxt干扰引起的振动。因此，求解上述问题强迫振动问题，可以转化为求解下面
三个定解问题：
I:
2
0
00
(0,0),(1.1)
0,0(0),(1.2)
,()(0).(1.3)
ttxx
xxl
t
tt
uauxlt
uut
uxuxxl
II:
2
12
0
00
(0,0),(2.1)
(),()(0),(2.2)
0,0(0).(2.3)
ttxx
xxl
t
tt
uauxlt
ututt
uuxl
III:
2
0
00
(,)(0,0),(3.1)
0,0(0),(3.2)
0,0(0).(3.3)
ttxx
xxl
t
tt
uaufxtxlt
uut
uuxl
设方程I的解为
I
u，方程II的解为
II
u，方程III的解为
III
u，则原定解问题的解为以
上三个定解问题解的和，即
第八章热传导方程的傅里叶解小结及习题答案
12
(,)
IIIIII
uxtuuu
方程I直接用分离变量法求解；方程II为非齐次边界条件，先将边界条件齐次化后用
分离变量法求解。
下面研究方程III的解法。
基本解法一将未知解展开为本征函数法
该方法的前提条件是必须知道对应齐次方程的本征函数，第七章第四节“非齐次方程的
求解”例题用该方法求解，但最后落脚点还是非齐次常微分方程，非齐次常微分方程的解法
用冲量法（基本方法三）或积分变换法（拉普拉斯变换法或傅里叶变换法）。
基本解法二非齐次方程齐次化找出特解(,)vxt
2
0
00
(,)(0,0),(3.1)
0,0(0),(3.2)
0,0(0).(3.3)
ttxx
xxl
t
tt
uaufxtxlt
uut
uuxl
令(,)(,)(,)uxtwxtvxt，保持原有的齐次边界条件不变，使得(,)wxt满足：
2(0,0)
(0,)0,(,)0(0)
(,0)(,0),(,0)(,0)(0)
ttxx
tt
wawxlt
wtwltt
wxvxwxvxxl
则(,)vxt满足常微分方程的边值问题：
2
0
(,)(0,0)
0,0(0)
ttxx
xxl
vavfxtxlt
vvt
该方法的关键在于找出特解(,)vxt，适用于(,)fxt比较简单的情形。
第七章习题第11题、第14题为非齐次方程，其中的自由项比较简单，可以用该方法求
解。
第七章第11题：
2sinh(0,0),(1)
(0,)0,(,)0(0),(2)
(,0)0,(,0)0(0).(3)
ttxx
t
uaubxxlt
utultt
uxuxxl
分析：由于方程（1）的非齐次项知识x的函数，就可以把特解函数也取为只是x的函
数，即令
第八章热传导方程的傅里叶解小结及习题答案
13
(,)(,)()uxtwxtvx
其中()vx满足：
()sinh()
(),()()
2
4
0005
b
vxx
a
vvl
（4）式的解两次积分很容易求出来。求出()vx后，再求(,)wxt的定解问题：
2(0,0)
(0,)0,(,)0(0)
(,0)(,0),(,0)(,0)(0)
ttxx
tt
wawxlt
wtwltt
wxvxwxvxxl
基本方法三冲量定理法
该方法的基本思路是将强迫振动问题转化为无穷多个自由振动的叠加，而每一个自由振
动的初始位移和初始速度为前一个自由振动产生的位移和速度。
2
0
00
(,)(0,0),(3.1)
0,0(0),(3.2)
0,0(0).(3.3)
ttxx
xxl
t
tt
uaufxtxlt
uut
uuxl
思路：将非齐次项(,)fxt（单位质量所受的力或者强迫力引起的加速度）分解成无穷
多个前后相继的瞬时力的叠加。t时刻的力(,)fx不会影响
时刻之前的振动情况，
(,)fx独自干扰产生的弦的振动位移(,)vxt满足：
2()
(0,)0,(,)0
(,)0,(,)(0)
ttxx
t
vavt
vtvlt
vxvfxxl
求解上述(,)vxt的定解问题，然后对(,)vxt按时间叠加后，持续强迫力(,)fxt所产生的振
动位移为：
(,)(,;)
0
d
tuxtvxt
对于热传导问题的非齐次方程，处理方法相同。
第八章热传导方程的傅里叶解小结及习题答案
14
第八章习题解答
第八章热传导方程的傅里叶解小结及习题答案
15
第八章热传导方程的傅里叶解小结及习题答案
16
第八章热传导方程的傅里叶解小结及习题答案
17
第八章热传导方程的傅里叶解小结及习题答案
18
第八章热传导方程的傅里叶解小结及习题答案
19
第八章热传导方程的傅里叶解小结及习题答案
20
第八章热传导方程的傅里叶解小结及习题答案
21
第八章热传导方程的傅里叶解小结及习题答案
22
第八章热传导方程的傅里叶解小结及习题答案
23
第八章热传导方程的傅里叶解小结及习题答案
24
第八章热传导方程的傅里叶解小结及习题答案
25
第八章热传导方程的傅里叶解小结及习题答案
26
第八章热传导方程的傅里叶解小结及习题答案
27
第八章热传导方程的傅里叶解小结及习题答案
28
第八章热传导方程的傅里叶解小结及习题答案
29
第八章热传导方程的傅里叶解小结及习题答案
30
第八章热传导方程的傅里叶解小结及习题答案
31