第二章 数学物理定解问题
什么是定解问题

§1.2 什么是定解问题1. 定解问题定解问题是根据已知物理规律求解特定物理过程的数学条件,它由泛定方程和定解条件两个部分组成,泛定方程也称为数学物理方程。
2. 泛定方程泛定方程是待解物理过程所遵循的物理规律的数学表达式,具体表现为某物理量关于时间和空间变量的偏微分方程,同一类物理过程遵循相同的物理规律,因此泛定方程反映一类物理过程的共性。
方程中物理量对时间变量的偏微分项反映物理过程的因果关联。
方程中物理量对空间变量的偏微分项反映物理过程的内部作用,或内在关联。
例1. 质点运动状态的演化问题在质点动力学问题中常求质点的运动轨迹,一旦求出运动轨迹,则一切与质点运动有关的物理量(如动能、动量、角动量等)都可求出。
质点的运动状态是由质点的位矢和动量完全确定,求质点运动轨迹的方法就是求解质点的运动状态随时间演变的过程,即由前一时刻的位矢和动量推算出下一时刻位矢和动量,从物理上看前后二时刻质点的运动状态的联系为dt t p m t r dt t r t r dt t r )(1)()()()(K K K K K +=+=+, dt t F t p dt t p t p dt t p )()()()()(K K K K K +=+=+ 因此,只要知道质点的受力情况就能由前一时刻的运动状态求出下一时刻的运动状态,这样的推演过程就是求解常微分方程F t r m K K =)(满足初始条件“0000)(,)(v t r r t r K K K K ==”的解。
§1.3 定解条件。
一、初始条件初始条件描述特定物理过程的起因,就t 这个自变数而言,如果泛定方程中物理量u 对t 最高阶偏导数是n 阶偏导数n n tu ∂∂,则要确定具体的定解问题,需要n 个初始条件。
例1:均匀细杆的导热问题满足的泛定方程为02=−xx t u a u ,则要确定具体的导热问题的解只需一个初始条件:)(0x u t ϕ==,即要已知初始温度分布。
第2章定解问题

第2章定解问题第2章定解问题1、何谓数理⽅程?按其描绘的物理过程,它可分为哪⼏类?2、何谓定解问题?它分为哪⼏类?试写出⼀维波动⽅程的Cauchy问题的数学表⽰。
3、何谓定解条件?它包括哪些内容?4、何谓边界条件?它分为哪⼏类?⼀个边界需⽤⼏个边界条件来描述?5、⽤数理⽅程来研究物理问题需要经历哪⼏个步骤?6、在静电场问题中,由介电常数分别为和的两种介质组成的系统的交界⾯S 处的衔接条件有⼏个?应如何表⽰?7、如何导出物理模型的数理⽅程?在推导弦的横振动⽅程时采⽤了哪些近似?由⼩⾓度近似我们得到什么结论?8、热传导⽅程的扩散⽅程有何共同和不同之处?9、在杆的纵振动问题中,若端⾃由,这个边界条件如何写?你能从Hooke定律出发证明吗?10、在杆的导热问题中,若端绝热,这个边界条件该如何写?你能从⼀物理定律出发证明吗?11、在热传导问题中,若热源密度不随时间⽽变化,则热传导⽅程会发⽣怎样的变化?12、在弦的横振动问题中,若弦受到了⼀与速度成正⽐的阻⼒,该阻⼒对于弦的振动问题是否起到了源的作⽤?若受到了⼀与位移成正⽐的回复⼒呢?第3章⾏波法1、⾏波法的解题要领是什么?它适合⽤来求解哪⼀类定解问题?为什么?2、⼀维波动⽅程的通解为什么含有两个任意函数?他们各个有怎样的形式和怎样的物理意义?靠什么确定他们的具体函数形式?3、公式是⽤⾏波法求解弦的横振动问题时推得的,能否⽤公式求解如下定解问题?请说明原因?4、能否⽤公式求解如下定解问题?5、能否⽤⾏波法求解如下定解问题?6、你能否根据直⾓坐标系中的导出球坐标中球对称情况下的的表达式请记住这个结论:7、何谓平均值法?你能通过引⼊球⾯的平均值,将三维的波动⽅程化为关于平均值的⼀维⽅程吗?8、在Poisson 公式中,?若已知9、对于定解问题除了可⽤Poisson 公式求解外?你能否有其他的求解法?10、在弦的横振动⽅程单位质量的弦所受的外⼒,若将则怎样的物理含意?它的量纲是什么?11、冲量原理的精神是什么?12、你能否⽤纯强迫振动的解来求解定解问题13、试述推迟势的物理意义,在推迟势中,若,且局限于⼀单位球内,则其中的体积分该如何计算?14、对于定解问题按下述⽅法进⾏求解是否正确?为什么?令使由公式可求得⽽显然,所满⾜的定解问题的解为所以,原定理问题的解为第4章分离变量法1、分离变量法的物理背景是什么?为什么能将未知函数表⽰为单元函数的乘积?2、分离变量法适于求解哪些定解问题?能⽤分离变量法求解⽆界问题吗?4、分离变量法有哪⼏个求解步骤?其中最关键的是哪⼀步?5、何谓本征值问题?以下两个定解问题是否构成本征值问题?(1)(2)6、仿照上章⽤冲量原理求解⽆界弦的纯迫振动的思想和⽅法,你能否写出⽤冲量原理求有界弦的纯强迫振动的公式?7、在将边界条件齐次化时,为什么通常可选辅助函数为X的⼀次式,⽽当问题的两个端点均有第⼆类边界条件时,必须选辅助数为X的⼆次式?8、在⽤分离变量法求解圆的Dirichlet问题时,需要将边界条件齐次化吗?为什么?9、在⽤分离变量法求解下述问题时,是否需将边界条件齐次化?如何齐次化?10、在柱坐标和极坐标中对分离变量,所得到的的⽅程为…其后为什么要注明…?它是怎样得来的?11、在扇形区域中,⽤分离变量法求Dirichlet问题应选择什么坐标系?所得到的的⽅程仍是…吗?为什么?12、在⽤分离变量法求解定解问题时,应如何选择坐标系?能在直⾓坐标系中求解吗?5章特殊函数>> 1)勒让德多项式1、⽅程是什么⽅程?你能写出它在中的⼀有限解吗?2、试述Legendre⽅程本征值问题的提法,其本征值、本征函数是什么?3、你能证明吗?你能由和之值算出吗?4、Legendre多项式的母函数是什么?何谓母函数法?它有哪些⽤途?5、Legendre多项式的归⼀化因⼦是什么?模是什么?你能得到⼀正交归⼀的Legendre多项式吗?6、积分和之值分别是多少?和7、你能将⽤Legendre多项式表⽰吗?8、你能否⽤关系式导出递推公式9、在球坐标系中,在轴对称的情况下,△u=0的变量分离形式的解是什么?在球内的解是什么?在球外的解呢?10、什么是缔合Legendre函数?它是否⼀定是多项式?为什么?11、试述缔合Legendre⽅程本征值问题的提法,其本征值和本征函数是什么?12、缔合Legendre函数的模和归⼀化因⼦是什么?13、是否等同于?与有何关系?你能否由的正交归⼀性导出的正交归⼀性?15、何谓球函数⽅程?它满⾜下列条件的特解是什么?16、独⽴的l阶球函数共有多少个?17、你能⽤两种不同的形式,写出在球坐标系中,在⾮轴对称的情况下△u=0的解吗?它们对于球内和球外的具体情况,⼜分别是怎样的呢?2)贝塞⽿函数1、⽅程叫什么⽅程?你能写出它的⼀有限解吗?2、何谓Bessel函数的零点?它与Bessel⽅程的何种本征值问题有关?有什么样的关系?3、Bessel函数的母函数是什么?当v不为整数时有⽆母函数?为什么?4、你能利⽤Bessel函数的母函数关系式推导出Bessel函数的递推公式吗?5、Bessel函数有⽆微分表达式?若有,试写出;若⽆,说明为什么?6、什么是三类柱函数?它们是否均满⾜Bessel⽅程?它们互相的关系是怎样的?7、第⼆、三类柱函数是否也满⾜Bessel函数递推公式?为什么?8、9、10、Bessel⽅程的通解是什么?其有限解是什么?11、什么是虚宗量的Bessel⽅程?它经过什么样的代换可变成Bessel⽅程?由此你能推得虚宗量的Bessel ⽅程的⼀个特解吗?12、什么是虚宗量的Bessel函数和虚宗量的Neumann函数?虚宗量Bessel⽅程的通解是什么?13、你能完整地写出在柱坐标中对分离变量后所得到的在柱体内的分离变量形式的解吗?14、⽅程在柱坐标系下分离变量,在什么样的边界条件下会出现虚宗量Bessel⽅程?虚宗量的Bessel⽅程是否会构成本征值问题?15、球Bessel⽅程是什么样的情况下出现的?它与半整数的Bessel⽅程有什么关系?你能理解式给出的⼏个函数是球Bessel ⽅程的特解吗?16、试述球Bessel⽅程本征值问题的提法,其本征值和本征函数是什么?17、你能写出在球坐标系中对所得到的分离变量形式的解吗?第6章积分变换法1、何谓积分变换法?他的解题步骤是怎样的?2、Fourier变换的定义是什么?它的存在条件是什么?你能由周期函数的Fourier级数⽽导出⾮周期函数的Fourier积分从⽽引⼊Fourier变换吗?3、试求函数的Fourier变换(a>0),你能利⽤Fourier变换的某些性质求出和吗?其中,a为常数,t为参变量。
数学物理方程_定解问题

根据边界条件确定任意函数 f:
令 故
规定,当
时
4、定解问题是一个整体
达朗贝尔公式的求解过程,与大家熟知的常 微分方程的求结果成完全类似。
但遗憾的是,绝大多数偏微分方程很难求出 通解;即是求出通解,用定解条件确定其中待 定函数往往更为困难。这说明,达朗贝尔公式 不适用于普遍的数学物理定解问题的求解?
(7.1.8)
称式(7.1.8)为弦的自由振动方程。
(2) 如果在弦的单位长度上还有横向外力 作用,则式(7.1.8)应该改写为
(7.1.9)
式中
称为力密度,为 时刻作用于
处单位质量上的横向外力
式(7.1.9)称为弦的受迫振动方程.
2、 均匀杆的纵振动
一根杆,只要其中任一小段做纵向移动,必然使 它的邻段压缩或伸长,这邻段的压缩或伸长又使 它自己的邻段压缩或伸长。这样,任一小段的纵 振动必然传播到整个杆,这种振动的传播是纵波.
泊松方程和拉普拉斯方程的定解条件不包含初始条件, 而只有边界条件. 边界条件分为三类:
1、在边界上直接给定未知函数 , 即为第一类边界条件.
2、在边界上给定未知函数导数的值,即为第二类边界条件.
3、在边界上给定未知函数和它的导数的某种线性组合, 即第三类边界条件.
第一、二、三类边界条件可以统一地写成
第二类边界条件 规定了所研究的物理量在边界外法线方向上方向导数 的数值
u n
x0 , y0 ,z0
f (x0, y0, z0,t)
(7.2.3)
第三类边界条件 规定了所研究的物理量及其外法向导数的线性组合在 边界上的数值
(7.2.4)
其中 是时间 的已知函数, 为常系数.
7.2.2 泊松方程和拉普拉斯方程的定解条件
数学物理方程第一章、第二章习题全解

18
数学物理方程与特殊函数导教·导学·导考
2δρ ut ( x , 0 ) = k ( c - δ≤ x ≤ c + δ) 在这个小段外,初速度仍为零, 我们想得到的是 x = c 处受到冲 击的初速度 , 所 以 最后 还 要 令 δ→ 0。此 外 , 弦是 没 有 初 位 移的 , 即 u( x, 0) = 0 , 于是初始条件为
3. 有一均匀杆 , 只要杆中任一小段有纵向位移或速度 , 必导致 邻段的压缩或伸长, 这种伸缩传开去, 就有纵波沿着杆传播, 试推导 杆的纵振动方程。
解 如图 1 9 所示, 取杆
长方向为 x 轴正向, 垂直于杆长
方向的 各截 面 均 用 它 的 平 衡 位 置 x 标记 , 在时刻 t, 此截面相对
u( x, 0) = 0 0,
ut ( x , 0 ) = δkρ,
| x - c| >δ | x - c | ≤ δ (δ→ 0)
所以定解问题为
utt - a2 uxx = 0
u(0 , t) = u( l, t) = 0 u( x, 0) = 0 , ut ( x , 0 ) =
0, | x - c| > δ δkρ, | x - c | ≤ δ (δ→ 0 )
16
数学物理方程与特殊函数导教·导学·导考
第一章 课后习题全解
1 .4 习题全解
1. 长为 l 的均匀杆 , 侧面绝缘 , 一端温度为零 , 另一端有恒定热
流 q进入 ( 即单位时间内通过单位截面积流入的热量为 q) , 杆的初始
温度分布是 x( l 2
x) ,试写出相应的定解问题。
解 见图 1 8, 该问题是一维热传导方程, 初始条件题中已给
u x
定解问题复习

选择合 适的坐 标系
非齐次方程, 齐次边界条件
特解法
非齐次方程, 齐次定解条件 特征函数法
齐次方程,齐 次边界条件 分离变量法
齐次边界非齐次初始条件下非齐 次方程的解法:
齐次定解条件非齐次方程的解:
方程类型 波动方程 Tn 通解
输运方程
齐次边界非齐次初始条件非齐次方 程的解:
方程类型 Tn
波动方程
解:
定解问题
令 U (x, t) = v (x , t) + w (x, t) ,代入定解问题
即 定解问题
定解问题
定解问题的特解为:
例题
有一长为 l ,侧面绝热而初始温度为零度的均匀细杆, 它的一端保持温度始终为零度,而另一端温度随时间直 线上升,求杆的温度分布。
设杆长方向为 x 轴,x = l 端保持温度始终为零度, x = 0 端 温度随时间直线上升,比例系数为常数 c ,则定解问题为:
2
3
研究细杆导热问题,初始时刻杆的一端温度为零度, 另一端 跟外界绝热,杆上初始温度为 ,试求无热源时细杆上 温度的变化。
长为 l ,两端固定的弦,在单位长度上受横向力 g(x) sinwx 的 作用下做小振动,已知弦的初始位移 和 速度分别为j (x) 和 f (x) ,求其横振动的规律。 有一长为 l ,侧面绝热而初始温度为零度的均匀细杆,它的 一端保持温度始终为零度,而另一端温度随时间直线上升, 求杆的温度分布。
解
令 U (x, t) = v (x , t) + w (x, t) ,代入定解问题
视 v(x, t) 为原方程的特解,考虑到非齐次边界条件,取
将 v(x, t) 代入原定解问题的边界条件,得
v x 0 B ct v x l Al B 0 A ct l
数学物理方程-福州大学-江飞-2.1热传导方程及其定解问题的导出

n
k
u n
k1 u
u1 uΒιβλιοθήκη k k1u nu1
一般形式:u
u n (x,y,z)
g(x, y, z,t)
或
u
u n
g
泛定方程:u t
a2
2u x2
2u y2
2u z2
f
柯 西
初始条件 u
a2u f
问 题
t0
初 边
u g
值边
热管道 1D : ut a2uxx f
t1
则有热源的热传导方程为 ut a2u f a2u F / c .
2. 扩散方程的导出
扩散物从浓度高流向浓度低
* Nerst扩散定律
在该点的扩散系数
扩散物在无穷小时段dt内沿法线方向流过一个无穷
小面积dS的质量dm与扩散物浓度沿曲面dS法线
方向的方向导数N 成正比,即
n
t1,t2
由能量守恒:Q流入 Q吸收
t2 k(x, y, z) udSdt
t1
n
N-L公式及交换下积分次
c(x, y, z)(x, y, z)[u(x, y, z,t2) u(x, y序, z,t1)]dxdydz
t2 t1
ctudxdydzdt
利用高维N-L积分公式,
左端 t2 k(x, y, z) udSdt
dm D(x, y, z) N dSdt
n
因此类似热方程推导:
t2 D(x, y, z) NdSdt
t1
n
(N(x, y, z,t2) N(x, y, z,t1))dxdydz
tN(x, y, z,t) x DxN x DyN x DzN
第二章定解问题
T (x,t) T (x) (t)
对x受力分析,由牛顿第二定律 得
T2 cos 2 T1 cos
T2 sin 2 T1 sin 1 F(x 1x,t)x (x)utt (x 2x, t) ( 2 1)
注意到
sin
tan tan2
ux 1 ux2
ux
sin 1 ux (x, t)
§2.1 引言
一、数学物理方程简介:
数学物理方程是指从物理问题中导 出的反映客观物理量在各个空间、时刻 之间相互制约关系的一些偏微分方程。 方程可以分为线性和非线性方程。
偏微分方程的基本概念:
u u u
mu
F(x1, x2,
, xn,u, x1 , x2 ,
,, xn
, x1m1x2m2
xnmn ) 0
(2)物理问题的数学抽象:
1)由于弦是“细长”的,所以 (x,t) t
忽略重力
2)由于弦“绷紧”于AB两点,这说明弦中各相邻部分之间有 拉力即“张力”作用;由于弦是“柔软”的,所以相邻小段张 力总是弦线的切线方向;
3)由于弦作“微小”的横向振动,故相邻点沿振动方向位移的 差别很小,即
u | ux || x | 1 无穷小量
若 f 0
称为弦的自由振动,振动过程中不受外力。
utt a2uxx
齐次波动方程
事实上,除了以上一维波动方程,像薄膜振动(二维),电 磁场方程(三维)等,均属于波动方程:
utt a2u f (x, y, t)
uxx
uyy
2u x2
2u y 2
utt a22u a2 (uxx uyy uzz )
§2.2 三类数理方程的导出
一、弦的横振动方程(波动方程的建立)
数学物理方程 2-3章课后部分习题答案 李明奇主编 电子科技大学出版社
数学物理方程 李明奇主编 电子科技大学出版社2-3章部分习题答案习题2.14.一根长为L 、截面面积为1的均匀细杆,其x=0端固定,以槌水平击其x=L 端,使之获得冲量I 。
试写出定解问题。
解:由Newton 定律: tt x x Sdxu t x YSu t dx x SYu ρ=-+),(),(,其中,Y 为杨氏模量,S 为均匀细杆的横截面积,x u 为相对伸长率。
化简之后,可以得到定解问题为:⎪⎪⎩⎪⎪⎨⎧-==========)(|,0|0|,0|)/(0002L x Iu u u u u a u Y u t t t L x x x xx xx tt δρρ。
习题2.23.设物体表面的绝对温度为u ,它向外辐射出去的热量,按斯特凡-波尔兹曼定律正比于4u ,即dSdt ku dQ 4=,设物体与周围介质之间,只有热辐射而无热传导,周围介质的绝对温度为已知函数),,,(t z y x ϕ,。
试写出边界条件。
解:由Fourier 热传导实验定律dSdt nuk dQ ∂∂-=1,其中1k 称为热传导系数。
可得dSdt u k dSdt nuk )(441ϕ-=∂∂-,即可得边界条件:)(441ϕ--=∂∂u k k nus。
习题2.34.由静电场Gauss 定理⎰⎰⎰⎰⎰⋅=⋅VsdV dS E ρε01,求证:0ερ=⋅∇E ,并由此导出静电势u 所满足的Poisson 方程。
证明:⎰⎰⎰⎰⎰⎰⎰⎰⋅=⋅=⋅VVsdV dV divE dS E ρε01,所以可以得到:0ερ=divE 。
由E divE ⋅∇=与u E -∇=,可得静电势u 所满足的Poisson 方程:2ερ-=∇u 。
习题2.42.求下列方程的通解:(2):;032=-+yy xy xx u u u (5):;031616=++yy xy xx u u u解:(2):特征方程:03)(2)(2=--dx dy dx dy解得:1-=dx dy 和3=dxdy。
数学物理方程课后参考答案第二章
第 二 章 热 传 导 方 程§1 热传导方程及其定解问题的提1. 一均匀细杆直径为l ,假设它在同一截面上的温度是相同的,杆的表面和周围介质发生热交换,服从于规律dsdt u u k dQ )(11-= 又假设杆的密度为ρ,比热为c ,热传导系数为k ,试导出此时温度u 满足的方程。
解:引坐标系:以杆的对称轴为x 轴,此时杆为温度),(t x u u =。
记杆的截面面积42l π为S 。
由假设,在任意时刻t 到t t ∆+内流入截面坐标为x 到x x ∆+一小段细杆的热量为t x s xuk t s x u k t s x u k dQ x x x x ∆∆∂∂=∆∂∂-∆∂∂=∆+221 杆表面和周围介质发生热交换,可看作一个“被动”的热源。
由假设,在时刻t 到t t ∆+在截面为x 到x x ∆+一小段中产生的热量为()()t x s u u lkt x l u u k dQ ∆∆--=∆∆--=111124π又在时刻t 到t t ∆+在截面为x 到x x ∆+这一小段内由于温度变化所需的热量为()()[]t x s tuc x s t x u t t x u c dQ t ∆∆∂∂=∆-∆+=ρρ,,3由热量守恒原理得:()t x s u u lk t x s x uk t x s t u c x t ∆∆--∆∆∂∂=∆∆∂∂11224ρ消去t x s ∆∆,再令0→∆x ,0→∆t 得精确的关系:()11224u u l kxu k t u c --∂∂=∂∂ρ或 ()()11222112244u u l c k xu a u u l c k x u c k t u --∂∂=--∂∂=∂∂ρρρ 其中 ρc k a =22. 试直接推导扩散过程所满足的微分方程。
解:在扩散介质中任取一闭曲面s ,其包围的区域 为Ω,则从时刻1t 到2t 流入此闭曲面的溶质,由dsdt nuDdM ∂∂-=,其中D 为扩散系数,得 ⎰⎰⎰∂∂=21t t sdsdt nuDM 浓度由u 变到2u 所需之溶质为()()[]⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰ΩΩΩ∂∂=∂∂=-=2121121,,,,,,t t tt dvdt t uC dtdv t u C dxdydz t z y x u t z y x u C M两者应该相等,由奥、高公式得:⎰⎰⎰⎰⎰⎰⎰⎰ΩΩ∂∂==⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎭⎫ ⎝⎛∂∂∂∂=21211t t t t dvdt t uC M dvdt z uD z y u D y x u D x M 其中C 叫做孔积系数=孔隙体积。
数学物理方法第2章
§3.1 (1+1)维齐次方程的分离变量法
一、有界弦的自由振动 由第2章的讨论可知,讨论两端固定弦的自由振动规律问题可以归结 为求解下列定解问题:
2u 2u a 2 2 (0 x l , t 0), 2 x t u x 0 0, u x l 0 t 0 , u ( x), u x 0 x l . t 0 t t 0
§3.4
非齐次边界条件的处理
前面所讨论的定解问题,无论方程是齐次的还是 非齐次的,边界条件都是齐次的。如果遇到非齐次边界条件 的情况,应该如何处理?总的原则是设法将边界条件化成齐 次的。具体地说,就是取一个适当的未知函数之间的代换, 使对新的未知函数,边界条件是齐次的。现在仍以一维波动 方程的定解问题为例,说明选取代换的方法。
物理学、力学、工程科学甚至经济和社会科学中等许 多问题都可以归结为偏微分方程的定解问题。第二章中我 们讨论了怎样将一个物理问题表达为定解问题,这一章以 及以下几章的任务是怎样去求解这些定解问题,也就是说 在已经列出方程和定解条件之后,怎样去求既满足方程又 满足定解条件的解. 从微积分学得知,在计算诸如多元函数的微分和积分 (重积分等)时总是把它们转化为单元函数的相应问题来 解决,与此类似,求解偏微分方程的定解问题也可以设法 把它们转化为常微分方程的定解问题来求解。分离变量法 就是这样一种常用的转化方法。在这一章中,我们将通过 一些实例,讨论分离变量法及其应用。
2.1.1 2.1.2
(2.1.3)
这个问题的特点是,偏微分方程是线性齐次的,边界条件也是齐次的。 求解这样的问题可以运用叠加原理。我们知道,在求解常系数齐次常微分方 程的初值问题时,是在先不考虑初始条件的情况下,求出满足方程的足够多 的特解,再利用叠加原理做出这些特解的线性组合,构成方程的通解,然后 利用初始条件来确定通解中的任意常数,得到初值问题的特解。这就启发我 们要求解定解问题(2.1.1)——(2.1.3),须首先寻求齐次方程(2.1.1) 满足边界条件(2.1.2)的足够多的具有简单形式(变量被分离的形式)的 特解,再利用它们做线性组合,得到方程满足边界条件的一般解,再使这个 一般解满足初始条件(2.1.3)。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
图2-2(a)
图2-2(b)
(1)由x=c点的初始位移,求出其它点的初始位移。设。由于[0,c]
段是直线,设,由x=0时=0和x=c时=h定出A=h/c和B=0。同样,
[c,l]段也是直线,设,由x=c时,=h和x=l时,=0,定出和,综上得
另外,由题意,显然。 (2)先求初始位移,设拉力与张力T平衡时点x=c得位移为h,由 ,所以
, 而x点的应变则是
, 图2-1 弹性杆的微小纵振动 由于振动是微小的(不超过杆的弹性限度),由胡克定律有
。 设杆的横截面为S(设为常数),则由牛顿第二定律,段的运动方程是
其中常数满足,并利用了胡克定律式,而且将函数在处展开为泰勒级数 并取了前两项。以除上式的两端后, 令取极限,得到
记
则方程最后变为
2、长为的均匀杆,侧面绝缘一端温度为零,另一端有恒定热流进入 (即单位时间内通过单位截面积流入的热量为),杆的初始温度分布 是,试写出相应的定解问题。
3、一均匀杆的原长为,一端固定,另一端沿杆的轴线方向拉长而静 止,突然放手任其振动,试建立振动方程与定解条件。
4、把下列方程化为标准型。 5、 若是两个任意二次连续可微函数,验证 满足方程
的热量守恒问题(设杆的侧面是绝热的),有
其中S为杆的截面积,常数满足,将(1)代入,并令取极限,得到 (3)
类似地,对于小段,有 , 满足。将(1)代入,并令取极限,得到 (4)
(3) 所谓牛顿冷却定律是指,在介质边界面S上,热流强度地外法向分量是 同介质温
与周围介质温度之差成正比,即 (5)
其中的比例常数b>0称为热交换系数。在x=l处,的方向就是x的正方向,所 以(5)成
计算得原方程为
即
(3) 由 ,所以特征方程为
(Ⅰ)若,则,该方程为双曲型的。
特征线为和,或写为及。
令
(1)
(2)
原方程化为
(3)
将减式得, 所以 ,代入得标准型为
(Ⅱ)若,则,这时方程为椭圆型的。 其特征方程为 特征线为 及。 于是令 计算可得原方程化为
或
即
。
(4) 因为,所以该方程是椭圆型的,其特征方程为
x=l端
为绝热(即热流为零);(2)两端均有热流流入;(3)两端以牛顿冷却定
律与周围介质(其
温度恒为)进行热交换。
解:(1) 显然。设热流强度为,由热学中的Fourier定律,热流强
与介质的温度的关系是
,
(1)
是x方向的单位矢量),将(1)代入,记得
ห้องสมุดไป่ตู้
(2)
(2)设x=0处的热流强度为,x=i处的热流强度为。考虑杆的小段
性的,故可设任何方向上单位长度上所受的张力T也是常数。在非边界的膜 上任取一小块
到,现在来分析此小块横向上的受力和运动情况: 沿平行于r的方向,设张力与平衡位置的夹角为,有 在r边上受力 ,
边上受力 , 沿平行于的方向,设张力与平衡位置的夹角为,有
边上受力 , 边上受力 , 在小振动近似下, 由此并根据牛顿第二定律可
,
(4)
,
(5)
。
(6)
由方程(5)、(6)消去p后所得得方程再与方程(4)消去,即得
,
(7)
,由方程(3)、(7)和,得到
,
(8)
由于流体是无旋的,即 , 必存在速度势,使得
(9)
这样,求三个未知函数归结为求一个未知函数u。方程(5)和(6)消去p
后,将
(9)代入,有
(10)
方程(4)与(10)消去s后,将(9)代入,得到
(3)
其中k和是常数,它们由流体的性质决定。
原则上五个方程(1)-(3)((1)是矢量方程,它相当于三个标
量方程)可以确定五个未知函数和p,然而它们是非线性的,其求解已
超出本课程的范围。这里我们设法将它们线性化。
设平衡时流体的密度为,引入稠密度以代替。因为声波在传播过程
中,都是小量,略去高阶小量后,方程(1)-(3)成为
故定解问题为
2 解:杆的初始温度分布是,即有初始条件:
由杆的一端温度为零,有边界条件:
杆的另一端有恒定热流,即:
故定解问题为:
3 解:以表示小段的质心位移,设为杆的横截面积,为杆的密 度,是在点处所受的力。
由牛顿第二运动定律,有 当时, ,,又因为 ,故有
令 ,可得振动方程为 由题意放手时即是振动的初始时刻,此时杆振动的速度为零,即 而端拉离平衡位置使整个杆伸长了,故整个系统的初始位移为 再看边界条件,一端固定即该端没有位移,有 另一端由于放手任其振动时未受外力,有 故振动方程与定解条件为:
4 解 (1) 因为,所以该方程是椭圆型的,其特征方程为
特征线为 和。故可令,
为计算方便,又令
原方程化为 改变自变量的记号得:
(2) 由于 (Ⅰ)若,则,该方程为双曲型。 其特征方程为 ,和 特征线为 和。故可令,, 所以原方程化为
(Ⅱ) 若,则,该方程为椭圆型。 其特征方程为 ,和 特征线为 和。故可令,,为计算方便,又令
。在x=0段,的方向与x的正方向相反,所以有 。
例6 试将方程
化成标准型。 解:因为,所以该方程是双曲型的,其特征方程为
特征线为 和。令, 。 例7 讨论方程 的类型,并化为标准型. 解:所以
故所给方程为双曲型的,其特征方程为 即 解之有 特征曲线为
令 则原方程化为
1、长为的弦两端固定,开始时在受冲量的作用,试写出相应的定解问 题。
, 由此得,利用(1)小题的结果,得到
同样有。 (3)由题意有。在点x=c,由冲量定理,有,即是
弦的质量密度);在其它点处。因此。 说明: 也可以用极限的形式写出。将冲量I看成均匀分布在小段
为小量)内,则在此小段上单位长度受到的冲量是,于是,由冲 量定理,得
在计算结果最后取。
导出长为l的杆的一维热传导问题下的边界条件:(1)x=0端保持恒温,
第二章 数学物理定解问题
2.2.1 例题分析
2.2 基础训练
例1 导出均匀弹性杆的微小纵振动方程,设杆的弹性模量(杆伸长 单位长度所需要的力)为E,质量密度为,作用于杆的外力密度为 F(x,t)
解:取x轴沿杆的轴线方向,以u(x,t)表示x点,t时刻的纵向位移。 使用微元法,考虑杆上的一小段的运动情况。以记杆上x点,t时刻的应 力(杆在伸缩过程中各点相互之间单位截面上的作用力),其方向沿x 轴,现在求杆上x点,t时刻的应变(相对伸长)。如图(2。1)所示, A’B’表AB段(平衡位置)在t时刻所处的位置,则AB段的相对伸长是
特征线为 和。故令, 原方程变为
, 或
即: 5 验证:由题意有 及 及
可得 6证明:由题意只有验证 是方程的解即可。将代入方程(1)左边
式中,并利用条件是方程(1)的解,可得 证毕。
6、验证线性奇次方程的叠加原理。即若 均是线性二阶齐次方程 (1)
的解,其中都只是的函数,而且级数收敛,其中为任意常数,并对可以 逐次微分两次,求证仍是原方程的解。
2.2.3 解答与提示 1 解:设弦的两端为:,由题意有 弦的振动方程为
定解条件为:
在点,由冲量定理有 ,即 ;在其它处,或写成
注意:也可用极限的形式写出:将冲量看成均匀分布在小段 (是无穷 小量),则在此小段上单位长度受到的冲量是,于是,由冲量定理得
这就是杆的纵振动方程,也是一维波动方程。
例2 导出无旋流体中声波的传播方程。
解:运用规律法,直接从现有的理想流体动力学方程和热力学的物
态方程出发导出声波方程。
设流体密度为,空间中一点(x,y,z)的流体元的速度为, 压强为,而连续
性方程是
;
(1)
理想流体的动力学方程是
;
(2)
此外,声波传播过程是绝热过程,它的物态方程是
(11)
方程(7)、(8)和(11)都是声波的传播方程。如果都不
随t变化,则方程(7)、(8)和(11)都变成Laplace方程:
,
(12)
,
(13)
,
(14)
这些方程给出的是稳定场方程。
例3 均匀、各向同性的弹性圆膜,沿圆周固定,试列出膜的横振动方程和
边界条件。
解:设圆膜的质量面密度为为常数(因为圆模是均匀的),又因为圆膜是 弹性各向同
除方程两边,并且令取极限,得 , ,
。 边界条件为: 。 对弦的横振动问题导出下列情况下的初始条件:
(1)弦的两端点x=0和x=l固定, 用手将弦上的点x=c (0<c<l)拉 开使之与平衡位置的偏离为h(图2。2(a),并设h<<l),然后放手;
(2)弦的两端点x=0和x=l固定,用横向力拉弦上的点x=c (0<c<l), 达到平衡后放手(图2。2(b)); (3)弦的两端点x=0和x=l固定,以槌击弦上的点x=c (0<c<l) 使之获得冲量