yyf§7.1 数学物理方程的导出

合集下载

数学物理方程和定解条件的导出

数学物理方程和定解条件的导出
dx
dA
u u 在 dt ) x dAdt k ( ) x dtdA 。 n x u u k ( ) x dx dtdA k ( ) x dx dtdA 。 n x
,在 dt 时间内电流流过该段导线所产生的热量
j 2 dtdAdx
为: Q I 2 Rdt j 2 (dA)2 由热平衡方程式
第六章
数学物理方程和定解条件的导出
6.1 波动问题
1. 一长为 l 的均匀细杆, x 0 端固定,另一端沿杆的轴线方向被拉长
b 后静止(在弹性限度内) ,突然放手任其振动,写出振动方程与定
解条件。 解: (1) 方程:
2u [ ( x dx) ( x)]s t 2 2u u ( x dx, t ) u ( x, t ) 2u dx 2 Y [ ] Y 2 dx t x x x Y utt u xx a 2u xx
u u 2u ( )dx gdx h dx dx 2 x x t t 2 2 u u u T 2 g h 2 x t t
因为 g 这项很小,可以忽略不计
2u T 2u h u 所以 2 0 t x 2 t
亦即 Tx dx Tx T 且 sin 1 tan 1
u T x T
u u ,sin 2 tan 2 x x x
x dx
u u 2u gdx h dx dx 2 x x t t x dx
2)初速度,在 x c 段,由动量定理:ΔP t Fdt I ,而动量的变化
1
t2
为ΔP mut ( x, 0) 2 ut ( x, 0) ,将两式联立,有 ut ( x, 0) 在 x c 段,没有受到外界作用,故 ut ( x, 0) 0,

梁昆淼_数学物理方法第7章

梁昆淼_数学物理方法第7章
考虑小振动
T1
x
x+x
x
ds (dx) (dy ) dx
2 2
sin 2 tg 2 u x
x x
sin 1 tg1 u x
x
(Tu x ) x dx (Tu x ) x dxutt
T2 T1 T
T (u x
T (u x
x dx
2 2 2 2 2 2 x y z


2 2 2 2 2 2 x y z
2 2 2 2 2 2 x y z

u utt 2 t
2
u ut t
u xx
2u x 2

有时记
2 2 2 2 2 x y
2 2 2 3 2 2 2 x y z
(二)、数学物理方程的导出
1、弦的横振动
弦的横向位移为 u(x,t)
dm ds
T2 cos 2 T1 cos1 0
u u( x ,t )
1
M1
M2
T2
2
T2 sin 2 T1 sin 1 dmutt
0
x
a
则热流强度与杆端 u|x=a 和周围介质温度有差关系
qx
xa
u k n
u x a k x
xa
h(u xa )
(u Hu x ) xa
H k /h
x=0 处
0
x
a
x 0
qx
x 0
u u k x 0 k n ( x) u h(u x0 ) k x 0 x
dV
E / 0

数学物理方程第三版第一章答案(全)

数学物理方程第三版第一章答案(全)

数学物理方程第三版答案第一章. 波动方程§1 方程的导出。

定解条件1.细杆(或弹簧)受某种外界原因而产生纵向振动,以u(x,t)表示静止时在x 点处的点在时刻t 离开原来位置的偏移,假设振动过程发生的张力服从虎克定律,试证明),(t x u 满足方程()⎪⎭⎫ ⎝⎛∂∂∂∂=⎪⎭⎫ ⎝⎛∂∂∂∂x u E x t u x t ρ 其中ρ为杆的密度,E 为杨氏模量。

证:在杆上任取一段,其中两端于静止时的坐标分别为 x 与+x x ∆。

现在计算这段杆在时刻t 的相对伸长。

在时刻t 这段杆两端的坐标分别为:),();,(t x x u x x t x u x ∆++∆++其相对伸长等于 ),()],([)],([t x x u xxt x u x t x x u x x x ∆+=∆∆-+-∆++∆+θ令0→∆x ,取极限得在点x 的相对伸长为x u ),(t x 。

由虎克定律,张力),(t x T 等于),()(),(t x u x E t x T x =其中)(x E 是在点x 的杨氏模量。

设杆的横截面面积为),(x S 则作用在杆段),(x x x ∆+两端的力分别为x u x S x E )()(x u x x S x x E t x )()();,(∆+∆+).,(t x x ∆+于是得运动方程tt u x x s x ⋅∆⋅)()(ρx ESu t x =),(x x x x x ESu x x |)(|)(-∆+∆+利用微分中值定理,消去x ∆,再令0→∆x 得tt u x s x )()(ρx∂∂=x ESu () 若=)(x s 常量,则得22)(tu x ∂∂ρ=))((x u x E x ∂∂∂∂即得所证。

2.在杆纵向振动时,假设(1)端点固定,(2)端点自由,(3)端点固定在弹性支承上,试分别导出这三种情况下所对应的边界条件。

解:(1)杆的两端被固定在l x x ==,0两点则相应的边界条件为 .0),(,0),0(==t l u t u(2)若l x =为自由端,则杆在l x =的张力x ux E t l T ∂∂=)(),(|lx =等于零,因此相应的边界条件为x u∂∂|lx ==0 同理,若0=x 为自由端,则相应的边界条件为xu∂∂∣00==x (3)若l x =端固定在弹性支承上,而弹性支承固定于某点,且该点离开原来位置的偏移由函数)(t v 给出,则在l x =端支承的伸长为)(),(t v t l u -。

第七章第一节数学物理方程的导出

第七章第一节数学物理方程的导出

q x x
dxdydz
D u dxdydz x x
❖ 同理沿y, z方向净流入量分别为
y
D
u y
dxdydz,
z
D
u z
dxdydz
❖ 根据粒子数守恒定律,如果平行六面体中没有源和汇. 则单 位时间内增加的粒子数=单位时间内净流入的粒子数. 即
u t
dxdydz
x
(1)要研究的物理量是什么? 杆沿纵向的位移
(2)被研究的物理量遵循哪些物理定 理?牛顿第二定律,Hooke定律
(3)按物理定理写出数学物 理方程(即建立泛定方程)
取杆长方向为x方向,垂直于 杆长方向的各截面均用它的平衡 位置x标记.
Px, t S
Px dx,tS
x x dx
在任一时刻t,此截面相对于平 衡位置的位移为u(x,t).
杆的受迫振动方程跟弦的受迫振动方程完全一样,只是其中
f (x, t)应是杆的单位长度上单位横截面积所受纵向外力。
(五) 流体力学与声学问题 流体力学中研究的物理量是流体的流动速度v、压强p和密度ρ。 对于声波在空气中的传播,相应地要研究空气质点在平衡位置 附近的振动速度v、空气的压强p和密度ρ。物体的振动引起周 围空气压强和密度的变化,使空气中形成疏密相间的状态,这 种疏密相间的状态向周围的传播形成声波。
第七章 数学物理定解问题
§7.1 数学物理方程的导出 §7.2 定解条件
在科学技术和生产实际中常常要求研究某个物理量(电场强度、 电势、磁感应强度、声压、杂质浓度)在空间的某个区域的分布 情况,以及它们怎样随着时间而变化。这些问题中的自变数不 仅有时间,而且还有空间坐标。
如波动微分方程
2y 1 2y x 2 u 2 t 2

数学物理方程第3版答案

数学物理方程第3版答案

2u u g [(l x) ] 。 2 x x t
5. 验证
u( x, y, t )
1 t x y
2 2 2
在锥 t x y >0 中都满足波动方程
2 2 2
2u 2u 2u 1 2 2 证:函数 u( x, y, t ) 在锥 t 2 x 2 y 2 >0 内对变量 2 2 2 2 t x y t x y
令 x+at=0 得 ( x) =F(2x)+G(0) 所以 F(x)= ( ) -G(0). G(x)= ( ) -F(0). 且 所以 F(0)+G(0)= (0) (0). u(x,t)= (

x 是 常 量 , Ex E也 是 常 量 .令a 2
2 2u u 2 u b a . t t 2 x 2
E

,则 得 方 程
1.
§ 2 达朗贝尔公式、 波的传抪 证明方程
数学物理方程答案
2 2 x u 1 x 2 u h 0常数 1 2 1 2 x h x h a t
3
u (t 2 x 2 y 2 ) 2 x x
数学物理方程答案
2u x 2
t
t

2
3 5 2 2 2 2 2 2 2 x y 3t x y x 2



同理
5 2u 2 2 2 2 2 t x y t x 2 2 y 2 2 y
F x
1 h x x 1 h d c 2 2a x 2
o
x
1 1 h d c Gx h x x 2 2a x 2

第一大节:课程介绍与三类典型方程的导出

第一大节:课程介绍与三类典型方程的导出
2u 2 a 3 u f (t , X ) 2 t
声波在三维弹性介质中 传播,介质的运动类似:
弹性介质受到微小扰动后的运动方程(波动方程)一般可表述为:
三个典型方程的导出-弦的横振动方程
提问一: 如果弦在粘稠的液体中运动,受到一与速度成正比的阻尼,其 运动方程会有什么变化?
u b t
其中拉格朗日余项
f ( n 1) ( ) Rn ( x ) ( x x0 ) n 1 (n 1)!
(在x0 与x之间)
三个典型方程的导出-弦的横振动方程
2 2u u 2 a f (t , x) 2 2 t x
令a
T

f (t , x)
g (t , x)
三个典型方程的导出-弦的横振动方程
(3)讨论任意一段 x 在
t
时刻的受力情况
外力
三个典型方程的导出-弦的横振动方程
(4)根据牛顿第二定律写出在
t 时刻运动方程
(5)消掉未知量并化简,考虑利用胡克定律
x x
M 1M 2

x
1 (u x ) 2
由于位移很小,相对位移(相邻两点Hale Waihona Puke 移之差与两点距离的比)也很小,所以:
长江大学地物学院教学课件
《数学物理方程》
《Mathematical Equations for Physics》
课程介绍
主讲教师:王婧慈 电子邮箱:851211wjc@
课程介绍
本课程的研究对象 物理问题中提出的数学方程,本课程中讨论的主要是偏微分 方程(Partial Differential Equation):含有多元未知函数的 偏导数的方程 。 它们反映了未知函数关于时间的导数和关于空间变量的导数 之间的制约关系,同时刻画了物理现象和物理过程的基本规律。 数学物理方程的概念 数学物理方程是指从物理、工程问题中,导出的反映客观物 理量在各个地点、时刻之间相互制约关系的一些偏微分方程。

方程的导出、定解条件

方程的导出、定解条件
2
ut a2 (uxx uyy uzz ) f ( x, y, z, t ),
uxx uyy uzz f ( x, y, z, t ).
(2)非线性方程:
(a)拟线性方程:方程对未知函数的最高阶导数总 体来说是线性的方程。比如
ut buux uxxx ,
ut uxx f (u).
t

x x
x
[ utt ( x, t ) Tuxx ( x, t ) F ( x, t )]dxdt 0.
表示单位质量在每点处所受的外力
仍有 x, t 的任意性,知
utt auxx ( x, t ) f ( x, t ) : F ( x, t ) / .
进一步推广到高维情况:
*目标函数: 弦上质点相对于平衡位置的位移图
OL l
x
L
*物理守恒律(转化方程等式):
牛顿第二定律:F ma.
或冲量定理:F t p mv.
教材采用
下面我们推导弦振动方程,先考虑无外力情况:
u
O
t 时刻示意图
OL l
x
x x x
*初始条件:设弦在初始时刻 t 0 时的位置和速度为
u( x,0) ( x), ut ( x,0) ( x), (0 x l ).
*边界条件(注意该提法的物理背景): (a)第一类边界条件(狄利克雷边界条件): 在前面的推导中,弦的两端被固定在 x 0 和 x l 两点,即 u(0, t ) u(l , t ) 0. (b)第二类边界条件(诺伊曼边界条件): 设弦的一端 x 0 处于自由状态,即可以在垂直于 x 轴 的直线上自由滑动,且未受到垂直方向的外力。由于在 边界右端的张力的垂直方向分量是 Tu x ,于是边界处应 有 ux x0 0 也可考虑更一般非零情况。

21-22讲 数理方程的导出

21-22讲 数理方程的导出

(4)由数理方程发展为定解问题。
第七章 数学物理定解问题
在研究某个物理量时,不仅要研究其随时间的变化,还要考虑 其所处的位置,因此这种方程为偏微分方程。 同一类物理现象满足同一物理规律(即同一类数学方程),这 称为共性;而对于具体问题,又有各自的个性。在解决这些具体 问题时,不仅要考虑其所处的环境(即边界条件),还要考虑其 所在的特定历史情况(即初始条件)。边境条件和初始条件称为 定解条件。 表达物理规律的偏微分方程称为数学物理方程,其解与具体条 件无关,在数学上称其为泛定方程。
第二篇 数学物理方程
自然界的许多现象都可以用一种数学模型来抽象,而解决这种 抽象的方程便是数学物理方程,该方程通常为偏微分方程。 数理方程是在微积分产生后形成的。18世纪,Taylor和Bernoulli 研究弦的横振动、Fourier研究热传导、Euler和Lagrange研究流体 力学、Laplace研究位函数时,均给出了相应的数学方程,并给出 经典解。到19世纪中叶,逐渐形成了经典偏微分方程理论。 到20世纪,又提出了大量的数学物理问题,也推动了数学的发 展,出现了许多数学分支:泛函分析、拓扑学、群论、微分几何 等。而这些又为解决数理方程提供了强有力的手段,使数理方程 得到很大发展,主要表现为: (1)由齐次偏微分方程推广到非齐次; (2)偏微分方程的类型出现交叉,如双曲-抛物型; (3)数理方程已深入到生物、化学、经济等领域;
Y
声波沿方向传播时,微元受到的合力为(p为压强):
p dFz [ p( x, y, z, t ) p( x, y, z dz, t )]dxdy dzdxdy. z
同理,声波沿X、Y方向传播时,微元受到的合力分别为:
dFx p dxdydz, x

7.1数学物理方程的导出

7.1数学物理方程的导出

§7.1 数学物理方程的导出1、拿图7-7的B 段弦作为代表,推导弦振动方程. 解答:如图图7-7B 段弦的动力学方程:112222112222cos cos 0sin sin ()()T T d u d u T T dm ds dt dt ααααρ-=⎧⎪⎨-==⎪⎩,考虑微小振动,有10α≈、20α≈、1cos 1α≈、2cos 1α≈、11sin tan x u x αα∂≈=-∂、22sin tan x dxu xαα+∂≈=-∂,∴12T T T ==,22x x dxu uu d uT TT dx ds x xx x dtρ+∂∂∂∂⎛⎫-+== ⎪∂∂∂∂⎝⎭,又2222()()11tan dy ds dx dy dx dx dx dx α⎛⎫=+=+=+≈ ⎪⎝⎭,得22220d u u T dt x ρ∂-=∂,若令2Ta ρ=,则222220d u u a dt x ∂-=∂. 3、弦在阻尼介质中振动,单位长度的弦所受阻力t F Ru =-(比例常数R 叫作阻力系数),试推导弦在这阻力介质中的振动方程.解答:如图,对于纵振动弦,阻力F 作用在纵向,B 段弦的动力学方程:221122221122cos cos 0sin sin ()()T T d u d u T T F dm ds dt dt ααααρ-=⎧⎪⎨-+==⎪⎩,考虑微小振动,有10α≈、20α≈、1cos 1α≈、2cos 1α≈、11sin tan x u x αα∂≈=∂、22sin tan x dxu xαα+∂≈=∂,∴12T T T ==,22t t x x dxu uu d uT TRu T dx Ru ds x xx x dt ρ+∂∂∂∂⎛⎫-+-=-= ⎪∂∂∂∂⎝⎭,又2222()()11tan dy ds dx dy dx dx dx dx α⎛⎫=+=+=+≈ ⎪⎝⎭,得22220t d u u T Ru dt x ρ∂-+=∂,若令2Ta ρ=,则222220t d u u R a u dt x ρ∂-+=∂. 4、试推导一维和三维的热传导方程(7.1.37)和(7.1. 38). 解答:由Fourier (热传导)定律qk u 可知, x 方向的分量xuq kx.如图,在有源的情况下,若仅在x 方向有热传导,从能量守恒定律得单位时间内传入体元dxdydz 的热量为(,)xxx dxxu c mq q dydz F x t dxdydz t ,∴(,)(,)x dxxu u u uc dxdydz k k dydzF x t dxdydzk dxdydz F x t dxdydz txxxx得 (,)u u ckF x t txx,对于均匀物体,热传导系数为常数,若令2kac和(,)(,)F x t f x t c,则222(,)u uaf x t t x (一维热传导方程).同理,在y 方向,有222(,)uu af y t ty ,在z 方向,有222(,)uu af z t tz即22(,,,)u a u f x y z t t(三维热传导方程)6、均质导线电阻率为r ,通有均匀分布的直流电,电流密度为j .试推导导线内的热传导方程. 解答:电功率密度2p j r =,由Fourier (热传导)定律x 方向的分量xu q kx,在有源的情况下,从能量守恒定律得单位时间内传入体元dxdydz 的热量为xx x dxxu c mq q dydz pdxdydz t, ∴2x dxxu u u uc dxdydzkkdydzpdxdydzk dxdydz j rdxdydz t x xxx,对于均质导线得 222u u ckj r tx .§7.2定解条件1、长为l 的均匀弦,两端0x和xl 固定,弦中的张力为0T ,在x h 点,以横向力0F 拉弦,达到稳定后放手任其自由振动,写出初始条件. 解答:参考§8.18.1齐次方程的分离变量质点平衡:210sin sin T T F θθ'''+=,21cos cos T T θθ'''= 考虑是微小纵振动,21coscos1T TT ;112200sin tan ,sin tan h h x l x θθθθ≈=≈=- 000210000000()sin sin ()F l x x h h ThlT T F TT F F h l x x l x x Tlθθ-'''+=⇒+=⇒=⇒=-- ∴00000000()0()()F l x hx x u x x x TlF x h x x l u l x l x l x Tl-⎧≤≤==⎪⎪⎨⎪≤≤=-=-⎪-⎩。

数学物理方程的导出

数学物理方程的导出

数学物理方程的导出过程主要介绍数学物理方程的建立方法.具体通过五种物理模型详细介绍数学物理方程的建立方法.其中弦的横振动、杆的纵振动以及传输线方程的建立是需要掌握的基本内容.为了描述定解问题的系统完整性,我们在对波动方程的定解条件也进行了讨论.(一)弦的横振动方程(均匀弦的微小横振动)弦的横振动问题是数理方程中的典型问题.它模型简单,且具有代表性.演奏弦乐用(二胡,提琴)的人用弓在弦上来回拉动,弓所接触的是弦的很小的一段,似乎只能引起这个小段的振动,实际上振动总是传播到整个弦,弦的各处都振动起来。

振动如何传播呢?1. 物理模型实际问题:设有一根细长而柔软的弦,紧绷于A,B两点之间,在平衡位置附近产生振幅极为微小的横振动(以某种方式激发,在同一个平面内,弦上各点的振动方向相互平行,且与波的传播方向(弦的长度方向)垂直),求弦上各点的运动规律。

2.分析:弦是柔软的,即在放松的条件下,把弦弯成任意的形状,它都保持静止。

绷紧后,相邻小段之间有拉力,这种拉力称为弦中的张力,张力沿线的切线方向。

由于张力的作用,一个小段的振动必带动它的邻段,邻段又带动它自己的邻段…,这样一个小段的振动必然传播到整个弦,这种振动传播现象叫作波。

我们考察一根长为且两端固定、水平拉紧的弦.讨论如何将这一物理问题转化为数学上的定解问题.要确定弦的运动方程,需要明确:(1) 要研究的物理量是什么?对于本模型是讨论弦的运动规律,并研究弦沿垂直方向的位移.(2)被研究的物理量遵循哪些物理定理?本模型所研究的物理量遵循牛顿第二定律.(3)按物理定理写出数学物理方程(即建立泛定方程)注意:由于物理问题涉及因素较多,往往还需要引入适当假设才能使方程简化.数学物理方程必须反映弦上任一位置上的垂直位移所遵循的普遍规律,所以考察点不能取在端点上,但可以取除端点之外的任何位置作为考察点.开始建立模型:① 模型实际上就是:柔软轻质细弦(“没有质量”的弦);弦是轻质弦(其质量只有张力的几万分之一)。

数学物理方程 ppt课件

数学物理方程  ppt课件

由能量守恒定律 c ρdx du=dQ =[q(x,t)-q(x+dx,t)]dt =-qx(x,t)dxdt
于是有 c ρut = -qx 由热传导定律 q(x,t) = -k ux(x,t) 代入前面的式子,得到 c ρut = k uxx ut = a2 uxx
a2 = k/(cρ)
ppt课件
于是有
T2 =T1=T ρuttdx=T[ux(x+dx,t)-ux(x,t)]
化简后得到
ρutt = T uxx utt = a2 uxx
uxxdx
a2 = T/ρ
6
波动方程
推广1
情况:受迫振动(考虑重力或外力)
分析:设单位长度所受到的横向外力 F(x,t),则dx段的受力为Fdx
方程:ρutt = T问题:扩散问题中研究的是浓度u在空间的分布和在时间中的 变化。 分析:扩散现象遵循扩散定律,即q= - D▽ u,q是扩散流强 度,D是扩散系数,▽u是浓度梯度。对于三维扩散问题, 考察单位时间内小体积元dxdydz的净流入量。
z
dz
y
dy
dx
x
o
ppt课件
9
扩散方程
在x,y,z方向上,单位时间内净流入量为
分析:设弦平衡时沿x轴,考虑 弦上从x到x+dx的一段,其质 量为ρdx。设弦的横振动位移 为u(x,t),则
α1
B
A
α2
C
ppt课件
由牛顿第二定律
ρdxutt=T2sinα2- T1sinα1 0 = T2 cosα2- T1 cosα1
微振动条件
cosα1 = cosα2= 1 sinα1 = tanα1 = ux(x,t) sinα2= tanα2 = ux(x+dx,t)

第1讲_数学物理方程的导出和定解条件(2015)

第1讲_数学物理方程的导出和定解条件(2015)

2u 2u 2u k 2 2 2 xyzt xyz c u x y z
F (x,y,z,t)
(14)
2u 2u 2u k 2 2 2 xyzt F ( x, y, z, t )xyzt xyz c u (17) x y z
4、小结:
物理上:

反映波动过程的波动方程 反映扩散过程的热传导方程 反映稳定状态的Poisson方程和Laplace方程 波动方程,在数学上属于双曲型方程 热传导方程,在数学上属于抛物型方程 Poisson方程和Laplace方程,在数学上属于椭圆型方 程
数学上:

泛定方程
例题
t
T2 2
弦中任意一小段 dx 在振动过程中
的受力情况为: 纵向(水平方向):
1 T1
o x
x dx
T2 cos 2 T1 cos 1
横向(竖直方向):
x
u T2 sin 2 T1 sin 1 b dx t x ~ x dx
∵弦在作横振动,∴由牛顿第二定律有
边界条件续:
当 f=0 时的边界条件称为齐次的。前面的三类 边界条 件分别为第一、第二、第三类齐次边界 条件。 边界条件的个数:
与初始条件的个数类似,等于方程中关于空间变量 偏导数的阶数。
边界条件的关键点:
只需给出恰当说明边界上的物理状况即可,而非整 个系统
§1.2 热传导方程与定解条件
三角函数系的正交性三角函数系的正交性1三角函数系kxkx上的积分等于零任意两个不同函数在正交cossinnxdxkx其中iiisinsinnxdxkxcoscosnxdxkx其中傅里叶系数sincos若有dxkxsincoscossincoscoskxdxnxkxdxnxsincos若有可得sinsinsincoskxdxnxkxdxnxsincos若有从而得到傅里叶系数把以上得到的系数代入三角级数该级数称为傅里叶级数sincos正弦级数和余弦级数一般说来一个函数的傅里叶级数既含有正弦项又含有余弦项

第一章定解问题及方程导出例

第一章定解问题及方程导出例

第一章 定解问题§1 大体概念1.数学物理方程:是指从物理问题中所导出的反映客观物理量在各个地址、各个时刻之间彼此制约的一些偏微分方程(有时也包括常微分方程和积分方程)2.数学物理方程的分类数学物理方程按其所代表的物理进程可分为如下三类:(1)描述振动和波动特征的波动方程f u a u tt +∆=2(2)反映输运进程的扩散(或热传导)方程f u D u t +∆=(3)描述稳固进程或稳固状态的poisson 方程h u -=∆其中 222222zy x ∂∂+∂∂+∂∂=∆ 22t u u tt ∂∂=,t u u t ∂∂= 而未知函数u (x , y , z , t )在三类方程中别离表示位移、浓度(或温度)和稳固现象特征;a 和D 表示波速和扩散(或热传导)系数;f 和h 是与源(汇)有关的已知函数,当f =0或h =0时,相应的方程称为齐次方程。

3.用数学物理方程研究问题的一般步骤(1)导出或写出定解问题(它包括数学物理方程和定解条件两部份)(2)求解已导出或写出的定解问题(3)对求得的解讨论其适应性(即解的存在性、惟一性、稳固性),并作出适当的物理解释4.求解数学物理方程的方式求解数学物理方程的方式大致能够分为如下几种:行波法(达朗贝尔法);分离变量法;积分变换法;Green 函数法;保角变换法;复变函数法;变分法;数值方式§2 数学物理方程的成立或推导1.成立(或推导)数学物理方程的步骤成立数学物理方程一般步骤step1:从所研究的系统中任取一单元体,分析该单元体与临近单元体之间的彼此关系;step2:按照有关的物理定律(如牛顿第二定律、能量守恒定律、奥—高定律等);用算式表达那个作用;step3:化简、整理即得所研究问题知足的数学物理方程。

2.成立(导出)方程时经常要用到的物理定律(1)Newton 第二定律:F=ma(2)Fourier 实验定律(即热传导定律),当物体内部存在温度差时会产生热量的流动。

第七章数学物理方程的导出和定解问题

第七章数学物理方程的导出和定解问题

第二部分 数学物理方程在物理学中,很多物理规律、物理现象、物理过程和物理状态的变化等都需要用微分方程来描述。

有些是用常微分方程来描述,如经典力学中质点和质点组的运动方程是常微分方程,而在研究连续介质和场(如电磁场、引力场、温度场)时会遇到偏微分方程,量子力学中的运动方程(薛定谔方程)也是偏微分方程。

在本课程的第二部分“数学物理方程”部分,我们主要学习如何求解在物理学中经常遇到的几种典型的偏微分方程。

在开始学习之前, 我们先简要复习一下:(1)什么叫常微分方程?(2)什么叫偏常微分方程?(3)什么叫微分方程的定解问题?1、常微分方程:联系一个自变量x ,该自变量的一个未知函数)(x y 和它的某些阶导数:出现在方程(或方程组)中的最高阶导数称为该常微分方程的阶数。

常微分方程中要求解的未知函数是一元函数(只有一个自变量)。

例: 22220 ()?d y xy x y x dx ++=⇒=( 二阶常微分方程)2、 偏微分方程:联系几个自变量,,,,...x y z t ,这些自变量的一个未知函数(,,,...)u x y z t 和它的某些阶偏导数:出现在方程中的最高阶偏导数称为该偏微分方程的阶数。

常微分方程中要求解的未知函数是一元函数(只有一个自变量), 而偏微分方程中要求解的未知函数是多元函数(有几个自变量)。

例: 2222230 (,)?u u u u u x y x y x y∂∂∂∂+++=⇒=∂∂∂∂( 二阶偏微分方程)3、微分方程的定解问题:给定一个微分方程,通常能找到很多不同的解。

为了唯一地确定一个微分方程的解,需要有一些辅助条件,这些辅助条件叫定解条件。

找出一个微分方程的满足某些特定定解条件的解,就称为定解问题。

一个微分方程,如果没有给出定解条件,就叫泛定方程(其解是不确定的)。

例: 22220 ()?d y xy x y x dx++=⇒= 这是泛定方程,其解不确定。

例:定解问题 2220'020()|0 ()?()1x x d y xy x dx y x y x y x ==⎧⎫++=⎪⎪⎪⎪=⇒=⎨⎬⎪⎪=⎪⎪⎪⎪⎩⎭( 定解条件:0'0()|0()1x x y x y x ===⎧⎫⎪⎪⎨⎬=⎪⎪⎩⎭ )第六章数学物理方程的导出和定解问题【刘连寿、王正清编著《数学物理方法》P107-121】物理学中经常遇到的几种典型二阶偏微分方程:波动方程、热传导方程(扩散方程)、泊松方程、拉普拉斯方程、亥姆霍兹方程, ………。

数学物理方程课件:第七章 数学物理方程的定解问题 (2)

数学物理方程课件:第七章 数学物理方程的定解问题 (2)

B.菲克定律
浓度梯度: u
u(x)
u(x dx)
x
扩积的散物流质强的度量:单q位时间通过单位面
dx
q Du
C. 扩散方程 连续性方程
u t
x
(uv x
)
0
q uv
带入菲克定律
u t
x
(uv
x
)
u t
qx x
u t
(D u ) 0 x x
D 均匀
u t
a2
2u x 2
0
a2 D
二、几种基本的方程 1.均匀弦的微小横振动
u(x,t)
y
x x x
y(x,t) u(x,t)
变 化
y(x x,t t) u(x x,t t)
x u(x x,t t)
A.弦的横振动 B.无穷小的一段弦 B C.受力分析和运动方程
弦的原长 s x 现长 s'
(x)2 (u)2 x
YS du dx
YSux
f (Sdx)utt
u
u dx
utt a2uxx 0
6.热传导方程
热传导: 热量从温度高的地方到温度低的地方转移。 热力学问题。
热力学第一定律: dQ dW dU
Q 热力学过程交换的热量
x
W 热力学过程外界对系统做
的功
u1
u2
U 系统的内能
q
热传导过程 dW=0,dQ dU
注:(x, y, z) 和 (x, y, z) 是空间座标的函数,在系 统的任何位置都是确定的!
例如 t=0: u(x,t 0) (x) (x) h
特定的时间, 变化的空间。
u
(
x)
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

F
所放出的热量;
ρ
= f0
称为热源强度,即:单位时间、单位 质量所放出的热量; 为按单位热容量计算的热源强度。
f ( x, t )
三维情况 相应地
ut − a Δ3u = 0 2 ut −a Δ3u = f ( x, y, z,t )
2
24
例三、稳定场方程 在输运方程中,当
ut = 0
时,方程变为
33
体积为
z ( h + η )( du + dx )
y
u
u + du
h
o
η
x
31
x x + dx
由水的不可压缩性,得
zhdx = z ( h + η )( du + dx )
= z [ hdu + hdx + η du + η dx ]
⇒ 0 = hdu + η du + η dx
略去二阶无穷小量 η du
以时间为自变量 的常微分方程
以时间和空间 坐标为自变量 的偏微分方程
9
物理和工程技术问题用偏微分方程 泛定方程 表达出来,叫做数学物理方程 物理问题 物理规律 边界条件 初始条件 满足的方程 (共性) 定解条件 (个性)
周围环境影响 系统初始状态
定解问题
10
§7.1 数学物理方程的导出
导出---“翻译” 与定解条件无关 一 导出步骤: i)确定物理量u; ii)从所研究的系统中划出一个小部分,根据物理 规律分析邻近部分和这个小部分的相互作用; iii)这种相互作用在一个短时间段如何影响物理 量u,把这种影响用算式表达出来。
可得
ρutt − Yuxx = 0
这就是杆的纵振动方程.
18
讨论 Y (1) 对于均匀杆, 和 ρ 是常数,上式可以改写成
utt = a uxx
2
(7.1.9)
其中
a
2
=
Y
ρ
这与弦振动方程具有完全相同的形式. (2)杆的受迫振动方程跟弦的受迫振动方程 完全一样, 只是其中 f (x, t) 应是杆的单位长度上单位 横截面积所受纵向外力
19
例二、输运方程: 研究一维热传导问题 考察一根均匀细杆内热量传播的过程:设细杆 横截面积A(常数),细杆中无热源和热汇,侧面绝 热。
o
l
x
20
1、物理量是温度 : 2、物理规律:
u = u ( x, t )
v (1)热传导定律(傅立叶定律) q = −k∇u v q是热流强度——单位时间内垂直通过单位

dx ⇒ ρ utt = Tu xx = 0
uxx =
ux x+dx − ux x
utt − a u xx = 0
2
自由振动方程 齐次一维波动方程
a =
2
T
ρ
a 为振动传播速度
16
一般情况:
非齐次项
utt − a2uxx = f ( x, t )
受迫振动方程
非齐次一维波动方程 f ( x, t ) : 力密度(强迫力),t时刻作用于单位质 量上的横向外力。 三维情况: 或 相应地
η 随x而异,且随t而
η
y
h o x
x + dx
x
29
(2)取x处的截面与x+dx处截面的水进行分析: x方向运动方程:
由于两处 η不同,这部分水前后方所受压力不等。
( ρ zηdx ) u
tt
= ⎡ − ρ g η x + dx + ρ g η x ⎤ η z ⎣ ⎦
= − ρ gzη xηdx
第二篇
数学物理方法
1
为全书中心内容: (1)将物理问题翻译(转化)成数学问题(偏微分 方程,积分方程或微分积分方程); 常见的三种方程:波动方程,输运方程,以及 稳定场(拉普拉斯)方程。
(2)该数学问题的求解。
2

一、数理方程简介: 1、数学物理方程:

数学物理方程是指从物理问题中导出的反映 客观物理量在各个地、时刻之间相互制约关系 的一些偏微分方程。 偏微分方程分为线性和非线性,这里主要 讨论二阶线性方程。
4
(3)十九世纪末到二十世纪初,其它方程如:
高阶方程:
utt = a uxxxx + f (x, t)
2
Schrodinger方程
∂Ψ h ih = − ΔΨ+U(r)Ψ 2m ∂t
2
5
二、用数理方法研究问题的步骤: 1、写出定解问题 包括: 泛定方程(共性,一般规律) 定解条件(初始,边界等) 把物理问题转化为数学语言 如:
u t t − a ( u xx + u yy + u z z ) = 0
2
2
utt − a Δ 3u = 0
2
齐次 非齐次
17
utt − a Δ3u = f ( x, y, z, t )
2、 均匀杆的纵振动
B 段的运动方程为
∂ux dx = ρ(Sdx)utt YSux x+dx −YSux x = YS ∂x
y′′(t ) − 4 y = 0
y = C1e + C2e
2t −2t
⎧ y ′′(t ) − 4 y = 0 泛定方程 ⎪ ⎨ y (0) = 0 定解条件 ⎪ y ′(0) = 4 ⎩
6
2、求解: 数理方程的求解方法大致有行波法、分离 变量法、积分变换法、格林函数法、保角变换 法、复变函数法、变分法。我们将在以后有选 择地阐述。 3、分析解答: 解出答案,需分析其意义及适定性。 适定性:指解是存在的、唯一的而且是稳定的。
7
三、数理方程的特点: 数理方程一方面联系着物理学中的许 多问题;另一方面又要运用数学中的许多成果。 所以它是数学和物理学之间的桥梁。
8
第七章
数学物理定解问题
物理规律或工程科学与技术问题的数学表达式中,有 许多是微分方程。

质点力学 —— 质点的位移 电路 —— 电流、电压 空间连续分布的物理场—— 静电场 电场强度或电势 电磁场 电场强度或磁场强度
面积的热量; 一维情况 (2)能量守恒定律;
du ∇u = dx
21
3、研究对象:在杆上任取一小段 A
x
x + dx
o
x x + dx
l
x
分析: 由能量守恒定律 dx段温度随时间的变化所需热量 =单位时间流入dx的热量
t t + dt

dQ净 = c ⋅ ρ Adx ⋅ ut dt
22
净流入的热量
2
T2 cosα2 − T1Fra bibliotekcosα1 = 0
1
14
对微小振动:
cos α1 ≈ cos α 2 ≈ 1
sin α1 ≈ tgα1 = ux
sin α 2 ≈ tgα 2 = ux
x
x + dx
ds =
( dx ) + ( dy )
2
2
= 1 + ( u x ) dx ≈ dx
2
15
⎧T2 − T1 = 0 ⎪ ⇒ ⎨ ⎪T2 u x x + dx − T1 u x x = utt ρ dx ⎩
dQ净 = Q x入 − Q x + dx出
= ⎡ − ku x ( x, t ) ⋅ Adt ⎤ − ⎡ − ku x ( x + dx, t ) ⋅ Adt ⎤ ⎣ ⎦ ⎣ ⎦
= kA ⎡u x ( x + dx, t ) − u x ( x, t ) ⎤ dt ⎣ ⎦
= kA u xx dxdt
⇒ cρAudx = kAuxxdx t
ut − a u xx = 0
2
k a = cρ
2
k ⇒ ut − uxx = 0 cρ
23
一般情况:当有热源、热汇存在时
ut − a u xx = f ( x, t )
2
F ( x , t ) 为t时刻,x处单位时间、单位长度上
F ( x, t ) f ( x, t ) = cρ
26
1. 静电场的电势 静电场中,电荷分布与电场强度满足方程
r ∇⋅E = ρ /ε
(1)
静电场是保守场,存在着势函数,设电势为u , r 则有
E = −∇u
代入(1)式,则得静电势满足的方程
Δ 3u = − ρ / ε
3
泊松方程
非齐次
对于不存在电荷的区域, ρ=0,静电势满 足方程 Δ u = 0 拉普拉斯方程 齐次
du ⇒ η = −h = − hu x dx

32
由①②,消去 η ,得
utt = ghu xx

即重力波的纵向运动方程; ②式微分
∂ ∂ ηtt = −h 2 u x = − h utt ∂t ∂x
2

代入①,有
ηtt = ghη xx

即为重力波的横向运动方程。
c = gh 表明浅水波在水深时比水浅时传播快。
(2)物理规律:
u = u ( x, t ) uv v F = ma
(3)研究对象:在弦上任取一小段
B ( x x + dx )
13
分析: 两端沿切向受力分别为 纵向: 横向:
u v u v T1 , T 2

T2 sinα2 −T1 sinα1 = ma = ( ρ ⋅ ds) utt
α2 α1
3
2、发展史: (1)十八世纪初(Taylor):
utt = a u xx + f
相关文档
最新文档