波动方程和行波法
第七章 行波法(一)
第七章 行波法
利用初值条件确定函数 F,G
u( x,0) ( x)
ut ( x,0) ( x)
F ( x) G ( x) ( x)
a[ F ( x) G( x)] ( x)
x
a[F ( x) G( x)] C ( )d
x0
其中
x
x1
x2
内,因此该三角区域称为
决定区域。
结论:达朗贝尔解表示沿x 轴正、反向传播的两列波速 为a波的叠加,故称为行波法。
第七章 行波法
影响区域、依赖区间、决定区域
波动是以一定的速度 a 向两个方向传播的。
如果在初始时刻 t=0,扰动仅仅在有限区间 [ x1 , x2 ] 上存在,则经过时间 t 后,扰动传到的范围为
x1 at x x2 at
第七章 行波法
无界弦振动的初值问题
2 2u 2 u x 2 a 2 x t u ( x, 0) ( x), u ( x, 0) ( x) t
第七章 行波法
2. 行波法的基本思想
这种方法是针对波动方程提出的。由于波动现象的普
1 过 x1 作斜率为 的直线 x x1 at a 1 过 x2 作斜率为 的直线 x x2 at a t 则 它们与区间 [ x1 , x2 ]
一起围成的三角形区域 中的任意一点 ( x, t ) 的 依赖区间都落在区间 [ x1 , x2 ]
x x1 at
x x2 at
遍性,对如何认识和解决波动问题,一直是物理学家和数 学家们长期探索的课题。 (1)波函数可写成位置和时间函数的分离形式,且波函数
是由无穷多个谐波分量叠加而成的,由此提出了分离变量
波动方程初值问题与行波法
1 x at 1 u d 2 2a x at 1
1 arctan( x at ) arctan( x at ) 2a
例4: 求二阶线性偏微分方程初值问题的解
uxx 2uxy 3u yy 0 2 u | 3 x , u y | y 0 0 y0
2 F 3 x G x 3 x F ' 3 x G ' x 0
1 F 3x G x C 3
9 2 F 3x x C ' 4 G x 3 x 2 C ' 4
P( x, t )
依赖区间
x at
x at
x
区间 [ x at , x at ] 为解的依赖区间。
2.决定区域 该区域中任一点(x, t )的依 赖区间都落在区间[c, d]内 部,因此解在此该区域中的 数值完全由区间[c, d]上的 初始条件决定。
t
x c at
x d at
例5 求二阶线性偏微分方程的通解
uxx 2sin xuxy cos xuyy 0.
2
解:特征方程为
dy
2
2sin xdxdy cos x dx 0
2 2
dy dy 1 sin x 1 sin x 0 dx dx
G(x-at)=G(x0+at-at)=G(x0)
u2 G ( x ) ( t 0)
O
at
u2 G ( x at ) ( t t0 )
x0
x x0 at
x
u1 F ( x at )
第七章 波动方程初值问题
x1 x0 at
即, f1(x - at) 表示波速为 a 的右行波
同理可知, f2(x + at) 表示波速为 a 的左行波. 因此,行波解为左行波与右行波的叠加. 三. 半无界弦的自由振动
utt a 2 uxx 0 u x0 0 u t 0 ( x ), ut
二. 行波解的物理意义 行波法的通解为:
u( x, t ) f1 ( x at ) f 2 ( x at )
对 f1(x - at),在 t0 时刻,x0 位置的波动位移为:
f1 ( x0 at0 )
若在t0+Δt 时刻, x1位置的波动位移也为 f1 ( x0 at0 ) 则:
t 0
a f1 ( x at ) x
f 2 ( x at ) t 0 a x
t 0
a f1 '( x ) a f 2 '( x ) y ( x )
对上式积分:
1 x x0 y ( )d [ f1 ( x ) f1 ( x0 )] [ f2 ( x ) f2 ( x0 )] (2) a
(1)
t 0
y ( x ) a f1 '( x ) a f 2 '( x )
1 x x0 y ( )d f1 ( x ) f 2 ( x ) c a
(2)
1 1 x c f1 ( x ) 2 [ ( x ) a x0 y ( )d ] 2 由 (1) (2) (x > 0) 解得: x f ( x ) 1 [ ( x ) 1 y ( )d ] c 2 2 a x0 2
波动方程和行波法剖析课件
目录 Contents
• 波动方程的基本概念 • 行波法的基本原理 • 波动方程的解析解法 • 波动方程的数值解法 • 行波法的应用实例
01
波动方程的基本概念
பைடு நூலகம்
波动方程的定义
波动方程
描述波动现象的基本数学模型,通常 用于描述物理场(如声场、电磁场、 水波等)随时间和空间的变化规律。
03
最后,通过迭代求解差分方程 ,得到波在每个网格点上的值 ,从而得到波的传播和演化过 程。
行波法的优缺点
优点
行波法简单易懂,易于编程实现,能够处理复杂的边界条件和初始条件,适用 于求解各种类型的波动方程。
缺点
行波法需要设定初始条件和边界条件,对于某些复杂的波动问题可能需要较高 的计算成本和精度要求。
水波传播的模拟
要点一
总结词
利用行波法模拟水波的传播,有助于研究水波的形成、演 化及对环境的影响。
要点二
详细描述
在水波传播的模拟中,行波法能够模拟水面的波动情况, 包括波浪的生成、传播和消散。通过调整参数,可以研究 不同条件下水波的传播规律,如风速、水深、地形等,对 于水文学、海洋学等领域具有重要意义。
03
波动方程的解析解法
分离变量法
将波动方程的解表示为若干个变量的 乘积或商的形式,以便分别求解。
VS
分离变量法是一种求解波动方程的常 用方法。通过假设波动方程的解可以 表示为若干个变量的乘积或商的形式 ,我们可以将一个复杂的偏微分方程 转化为若干个简单的常微分方程,从 而方便求解。
积分变换法
利用积分变换将波动方程化为易于求解的形式,再进行逆变换得到原方程的解。
地震学
用于模拟地震波的传播和反射,进行地震预 测和地球结构研究。
波动方程的特征线法
作变换 1 ( x, y ), 2 ( x, y ),
在区域Ω上作此变换下,可化简方程(1),甚至可求得其解. 此变换称为特征变换.
例1 一端固定的半无界弦的自由振动问题
2u 2u a2 0 ( t 0,0 x ), 2 2 x t u t 0 : u ( x ), ( x ) ( 0 x ), t 0 t t 0 x 0 : u 0.
举例
2u 2u a 2 2 , x R, t 0 t 2 x u ( x, 0) 1, xR 2 ut ( x, 0) x ,
例4:
例5:
2u 2u a2 2 , 2 x t u ( x,1) cos x, ut ( x,1) 0,
例2:
2u 2u 2 a2 2 0 t x u t 0 cos x, ut t 0 x
解:由达朗贝尔公式
( x at) ( x at) 1 x at u ( x, t ) ( )d 2 2a xat
cos(x at) cos(x at) 1 x at d 2 2a x at
此公式的意义在于把定解问 题的解表示为左、右行进波 相叠加,这种方法称为“行 波法”。
D’Alembert公式
1 1 x at u ( x, t ) ( x at ) ( x at ) ( )d 2 2a x at
注 : 当 ( x ) C 2 ( R ), ( x ) C 1 ( R )时, 初值问题( I )存在唯一的解 u( x , t ),由d ' Alembert 公式给出.
数理方程3.1 一维波动方程的初值问题
§3.1 一维波动方程的初值问题
21 5 x − y = C1 x + y = C , 2
特征线为
作变量代换
ξ = 5x − y η = x + y,
原方程化为
uξη = 0
其通解
u = F (5x − y ) + G(x + y )
利用初始条件可得
F (5x) + G(x) = 5x, −F ′ (5x) + G′ (x) = 0
utt = a2 uxx , u(x, 0) = 1 , ut (x, 0) = 0, 1 + 4 x2 − ∞ < x < +∞, t > 0, − ∞ < x < +∞
1 . 1+4x2
从静止开始
由达朗贝尔公式得
1 1 1 1 1 u(x, t) = [ϕ(x + at) + ϕ(x − at)] = · + · 2 2 2 1 + 4(x + at) 2 1 + 4(x − at)2
其中ϕ(x), ψ (x)分别表示初值位移和初始速度. 1. 泛定方程的通解 x 2 dx 2 采用第2.1节中的方法, 特征方程为( d (特征 dt ) − a = 0, 特征线方程为 dt = ±a, 其通解 线)为x + at = C1 , x − at = C2 , 作变换
ξ = x + at, η = x − at,
第三章
波动方程的初值问题与行波法
§3.1 一维波动方程的初值问题
本节思路: 无界弦的自由振动(utt = a2 uxx , −∞ < x < +∞): 经非奇异变换化为标准型后直接积分得通解, 代入初始条件得特解(达朗贝尔公式) 无界弦的受迫振动(utt = a2 uxx + f, −∞ < x < +∞): 由叠加原理分解为: 齐次问题+零初值的非齐次问题(由齐次化原理得解) 半无界弦的振动问题(utt = a2 uxx + f, 0 < x < +∞): 以某种方式延拓f 及初始函数, 转成无限长的弦的振动, 求出解后限制在半无界区域上.
数学物理方程:第3章 波动问题的行波法
第3章 波动问题的行波法§3.1 二阶线性方程的分类与化简本节讨论:①两个自变量方程的分类与化简,②多个自变量方程的分类与化简⒈ 两个自变量方程的分类与化简二阶方程的一般形式 二阶变系数方程可写为1112220(,)2(,,,,)(,)xx xy yy x y Lu x y a u a u a u x y u u u f x y =+++Φ= (3.1.1)式中:11a 、12a 、22a 为x 、y 的函数,0(,,,,)x y x y u u u Φ为低阶导数项。
公式关于二阶导数项为线性的,即称方程为准线性的。
若0(,,,,)x y x y u u u Φ关于u 及其x u 、y u 为线性的,则称方程为线性的。
方程的变换 为了简化上述方程,作可逆变换:(,)(,)x y x y ξξηη=⎧⎨=⎩, (,)0(,)J x y ξη∂=≠∂, (,)(,)x x y y ξηξη=⎧⎨=⎩(3.1.2) 代入方程中,不难得到:11122212(,,,,)(,)Lu A u A u A u u u u f ξξξηηηξηξηξη=+++Φ= (3.1.3)式中: 22111112222x x y yA a a a ξξξξ=++ (3.1.4) 12111222()x x x y y x y y A a a a ξηξηξηξη=+++ (3.1.5)22221112222x x y yA a a a ηηηη=++ (3.1.6) 我们化简的目的是使得二次项的项数尽量少,并且值尽量为简单(如0ij A =或1ij A =±)。
顾及ij A 的表达式,取关于z 的一阶非线性偏微分方程2211122220x x y y a z a z z a z ++= (3.1.7)若该方程有解),(1y x z ϕ=、),(2y x z ψ=,则110A =及220A =;公式大大简化了。
《数学物理方程》第四章§1
2/16
2u 2u a2 2 t 2 x
2u 2u a2 2 0 t 2 x
2 2 ( 2 a 2 2 )u 0 t x
0 1 0 a 2
dx 令 dt
2 a 2 0
a
x at x at
t t a a 1 1 x x
0 a a a 1 1 a 1 0 a 2 1 1
0 2a 2
《百科全书》不仅在于提供知识,而更重要的在 于改变读者的思想。 向前进,你就会产生信念 ————达朗贝尔
达朗贝尔脱下了微分学的神秘外衣 ————马克思
5/16
2u 2u a2 2 t 2 x
u( x , t ) = f1(x + at ) + f2(x – at )
u t 0 u ( x ), t
2a 2 0
3/16
0 1 0 a 2
0 2a 2
2a 2 0
( t , x ) ( , )
2u 2u a2 2 0 2 t x
2u 0
2u 4a 2 0
x , x [0,1 / 2] ( x ) 1 x , x [1 / 2,1] 0, 其它
随着时间的推移, u2 的图形以速度 a 向x 轴正方向 移动. 所以,u2表示一个以速度a 沿 x 轴正方向传播 的行波,称为右行波。
8/16
2u 2u a2 2 t 2 x u u t 0 ( x ), t
数学物理方程03_波动方程初始问题的求解【OK】
数学物理方程
将上述初始条件代入达朗贝尔公式,即可得到:
x at x 1 1 [ ( x at ) ( x at )] ( s ) ds , t 2 a x at 2 a u ( x, t ) 1 [ ( x at ) (at x)] 1 x at ( s )ds, t x 2 2 a at x a
( x at ) 代表以速度a 沿x 轴负向传播的波
1
1 x at b. 只有初始速度时: u ( x, t ) ( )d x at 2a
u( x, t ) 1 ( x at ) 1 ( x at )
1 ( ) 为 ( ) 的积分原函数。
结论:达朗贝尔解表示沿x 轴正、反向传播的两列波速
第 3章
波动方程初始问题的求解
——行波法 (达朗贝尔公式) (特征线积分法)
1
数学物理方程
达朗贝尔公式(行波法)[一维问题]
通解法中有一种特殊的解法―行波法, 即以自变量的 线性组合作变量代换,进行求解的一种方法,它对波动方 程类型的求解十分有效. 1 基本思想: 先求出偏微分方程的通解,然后用定解条件确定特解。 这一思想与常微分方程的解法是一样的。 2 关键步骤: 通过变量变换,将波动方程化为便于积分的齐次二阶 偏微分方程。
代入通解得: u( x, t ) [ ( x at ) ( x at )]
x at
x at
( s)ds
达朗贝尔公式
(3.1.2) 5
数学物理方程
(4)达朗贝尔公式的意义: a. 只有初始位移时,u ( x, t ) ( x at ) ( x at ) 2 ( x at ) 代表以速度a 沿x 轴正向传播的波
15波动(横波、纵波、行波、简谐波、波长、波速、波动方程)
•液体和气体中 纵波 u B / B 容变弹性模量。
六、注意几点
1、周期、频率与介质无关,与波源的相同。 波长、波速与介质有关。
2、不同频率的同一类波在同一介质中波速相同。
3、波在不同介质中频率不变。
9
4.振动与波动的区别 •振动是表示一个质点的运动。 •波动是表示一系列质点所作的运动。
初位相不为0时:
y(x,t) Acos[(t x) ]
u
2 , 代入
T
y
A cos 2 Tt
x Tu
Tu 代入
y
A
cos 2 Tt
x
1 代入
T
y
A
cos2
t
x
t
显然质点振动速度与波速 u = 20m/s 不同。
上例中条件是已知 t = 0 时刻的波动方程。
如果t = 0时,波源 x = 0 点的振动方程为:
y 4102 cos(100t 2)m
波速不变。波动方程应该如何写?
y 4102 cos(100t 5x 2)m x>0
o
t
y x /4
o
t
y x /2
o
t
y x 3 / 4
o
t
15
3.当 t c
(常数)时 ,y f (x)
为某一时刻各质点 的振动位移,波形 的“拍照”
y t 0
o
x
y t T /4
o
x
y t T /2
第一章 波动方程和行波法
如弦振动方程: utt a2uxx 0
其初始条件为: 同一时刻( t 0 )情况
u ut
t0 (x) t0 (x)
初始位移 初始速度
注意:( a)初始条件应是整个系统的初始
状态,而不是系统中个别点的初始状态。
31
如:一根长为 l 的两端固定的弦,用手把它的 中点朝横向拔开距离h,然后放手任其振动( 初始时该就为放手的时刻),则初始条件应为:
Tux (x0 0,t) Tux (x0 0,t) F (t) ②
①、②合称为衔接条件,这时振动问题适定。
42
再如,不同材料组成的杆的振动,在
衔接处的位移和能量相等,即:
u u 1 x x0
2 x x0
E1u1x xx0 E1u2 x xx0
u1(x,t), u2 (x,t) :杆的两部分位移. E1 , E2 :两部分的杨氏模量.
27
二、定解条件的提出 1、必要性。导出方程后,就得对方程进行求
解。但是只有泛定方程不足以完全确定方程的 解,即不足以完全确定具体的物理过程,因为 具体的物理过程还与其初始状态及边界所受的 外界作用有关,因而必须找一些补充条件,用 以确定该物理过程。
28
从物理角度看:泛定方程仅表示一般性(共 性),要为物体的运动个性化附加条件。
10
由于张力的作用,一个小段的振动必带动它 的邻段,邻段又带动它自己的邻段,这样一个 小段的振动必然传播到整个弦,这种振动的传 播现象叫作波。弦是轻质弦(其质量只有张力 的几万分之一)。跟张力相比,弦的质量完全 可以略去。
11
① 模型实际上就是:柔软轻质细弦(“没 有质量”的弦)
② 将无质量的弦紧绷,不振动时是一根直 线,取为 x 轴。
14 波动方程的行波解
解:将初始条件代入达朗贝尔公式
u ( x)
( x at ) 1 [ e 2
2
2
e
( x at ) 2
]
1 2a
( x at ) 1 [ e 2
e
( x at ) 2
]
1 2
x at x at
x at
2ase
s 2
ds
( x at ) 1 [ e 2
右行波解:
u( x, t ) ( x at)
左行波方程
左行单波方程初值问题(其中常数 a 0 )
u u 0 t 0, x R a x t u t 0 ( x)
练习
求左行波解?
左行波方程
左行单波方程初值问题(其中常数 a 0 )
u u 0 t 0, x R a x t u t 0 ( x)
练习 答案:
求左行波解?
u ( x, t ) ( x at)
左行波解
左行单波方程初值问题(其中常数 a 0 )
u u 0 t 0, x R a x t u t 0 ( x)
2
2
e
( x at ) 2
] 1 [ e 2
x at x at s 2
e
s 2
ds2
e ( xat )
x at
例 特征边值问题(Goursat问题) 2 2u u 2 t 2 a x 2 , t x t , t 0 u ( x), x 0 x at 0 其中 (0)= (0) u x at 0 ( x), x 0 解:将定解条件代入通解 u( x, t ) f ( x at ) g ( x at ),
数学物理方程复习
数学物理方程复习一.三类方程及定解问题(一)方程1.波动方程(双曲型)Utt = a2Uxx+f; 0<x<l,t>0U(0,t)= Φ1(t);U(l,t)= Φ2(t);U(x,0)= Ψ1(x);Ut (x,0)=Ψ2(x)。
2.热传导方程(抛物型)Ut = a2Uxx+f; 0<x<l,t>0U(0,t)= Φ1(t);U(l,t)= Φ2(t);U(x,0)= Ψ1(x).3.稳态方程(椭圆型)Uxx +Uyy=f; 0<x<a;0<y<b;t>0.U(0,x)= Φ1(x);U(b,x)= Φ2(x);U(y,0)= Ψ1(y);Ut (y,a)=Ψ2(y)。
(二)解题的步骤1.建立数学模型,写出方程及定解条件2.解方程3.解的实定性问题(检验)(三)写方程的定解条件1.微元法:物理定理2.定解条件:初始条件及边界条件(四)解方程的方法1.分离变量法(有界区域内)2.行波法(针对波动方程,无界区域内)3.积分变换法(Fourier变换Laplace变换)Fourier变换:针对整个空间奇:正弦变换偶:余弦变换Laplace变换:针对半空间4.Green函数及基本解法5.Bessel函数及Legendre函数法例一:在弦的横震动问题中,若弦受到一与速度成正比的阻尼,试导出弦阻尼振动方程。
解:建立如图所示的直角坐标系,设位移函数为U(x,t),取任意一小段△x进行受力分析,由题设,单位弦所受阻力为b U t(b为常数),在振动过程中有△x所受纵向力为:(T2COSa2-T1COSa1)横向力为:(T2SINa2-T1SINa1-b U t(x+n△x))(0<n<1). T2,T1为△x弦两端所受的张力,又因为弦做横振动而无纵振动,由牛顿定律有T2COSa2-T1COSa1=0,T2SINa2-T1SINa1-b(x+n△x)U t=p U tt(x+n△x)△x在小的振动下SINa1≈TANa1=Ux(x,t), SINa2≈TANa2=Ux(x+△x,t),COSa2≈COSa1≈1,T=T1=T2.(ρ是密度)即(T/ρ)[ Ux (x+△x,t)- Ux(x,t)]/ △x-(b/ρ) U t(x+n△x,t)即令△x→0时有:U tt+ aU t=a2U xx例二:设扩散物质的源强(即单位时间内单位体积所产生的扩散物质)为F (x,y,z,t),试导出扩散方程。
数学物理方法16.1 行波法1-波动方程
( )d xat
a[ f1(x at) f1(x at)] a[ f2 (x at) f2 (x at)]
1
x at
( )d
a xat
[ f1(x at) f2 (x at)] [ f1(x at) f2 (x at)]
确定待定函数(法二)
待求的?
1
x
(v)dv
0
f1(0) f2 (0) 2
能消去吗?
f2
(
x)
(x) 2
1 2a
x
(v)dv
0
f1(0) f2 (0) 2
待求的解为
u f1 (x at) f2 (x at)
确定待定函数(法一)
f1
(x)
(x) 2
1 2a
x
(v)dv
0
f1(0) f2 (0) 2
(x) 1
那么,可得原问题的解为
u(x,t) 1 [(x at) (x at)] 1
x at
(v)dv
2
2a xat
确定待定函数(法二)
(x) (x)
f1(x) f2 (x) af1(x) af 2(x)
有何关联?
观察第一个方程,和待求解 u f1(x at) f2 (x at)
上述方程组中:4个待定函数,3个方程, 因此,不能直接求解各个待定函数。
u f1(x at) f2 (x at) 整体思想
确定待定函数(法二)
(x at) (x at)
[ f1(x at) f2 (x at)] [ f1(x at) f2 (x at)]
1
x at
( )d
行波法:算例1
2u u(tx2 ,0)
现代数学物理方程
这就是微分方程的适定性问题。
2、验证
u( x , y, t )
2
1 t x y
2 2
在锥
t x y 0
2 2 2
中都满足波动方程
u
2
t
2
u
2
x
2
u
2
y
2
.
证明:在该锥内
u t
2
(t x y )
2 2 2
3 2
t
3 2 5 2
又
sin 1 tg 1 sin 2 tg 2
u( x x , t )
.
于是得运动方程
x
u
2
t
2
g [ l ( x x )]
u( x x , t ) x
[l x ]
u( x , t ) x
u
2
[ l ( x x )] g
u( x , 0) t aF '( x at ) aG '( x at ) t 0 aF '( x ) aG '( x ) ( x ).
aF '( x ) aG '( x ) ( x ).
两边对 x 积分:
aF ( x ) aG ( x ) C
u
2
t
2
c u
2
这里c 通常是一个固定常数,代表波的传播速率。 在针对实际问题的波动方程中,一般都将波速表 示成可随波的频率变化的量,这种处理对应真实 物理世界中的色散现象。
(2)方程的导出 均匀弦的微小横振动 理想化假设:
第四讲行波法dhh
二阶线性微分方程的特征方程
定义2 考虑下面二阶线性微分方程 (3) (4)
2u 2u 2u u u a11 2 2a12 a22 2 b1 b2 cu f 0 xy x y x y
方程 即
a11 (dy) 2 2a12 dxdy a22 (dx) 2 0
数学物理方程
utt a2uxx 4a2u
2 2u u 2 a u 0 2 2 t x
u 0 d f ( )
u f ( )d f 2 ( ) F ( ) G ( )
u( x, t ) F ( x at ) G( x at )
解:齐次方程直接利用达朗贝尔公式:
1 1 x at 2 u sin( x at ) sin( x at ) d 2 2a x at t 2 2 2 sin x cos at (3 x a t ) 3
数学物理方程
x x utt 9u xx e e 例2 求定解问题: u ( x, 0) x, ut ( x, 0) sin x
式中 ( x), ( x) 均为已知函数,表示初始位移和初始速度。 特征线族
dx dx dx 2 a 0, a 0 a 0 dt dt dt 1 1 t x c1 , t x c2 a a
2
x at c1 , x at c2
解:一维非齐次波动方程初值问题解的Kirchoff公式
a 3, f ( x, t ) ex e x , ( x) x, ( x) sin x
关于x奇函数
u关于x奇函数
数学物理方程
3-2 延拓法求解半无限长振动问题 • (一)半无限长弦的自由振动问题
大学物理第二章行波,波动方程
横波 从图上可以明显看出在横波中各质元发生切变, 外形有波峰波谷之分
横波只能在弹性固体中传播
纵波
在纵波中,各质元发生长变或体变, 因而媒质的密度发生改变,各处疏密不同, 所以纵波也叫疏密波。
纵波在气体、液体、固体媒质中都可以传播
4. 波的特征
(1) 不管是横波还是纵波,在波传播的过程中, (2) 媒质中各质元均在各自的平衡位置附近振
二.描述波的物理量
1. 周期 T、频率 ν
波是机械振动的传播,在传播的过程中, 媒质的各个质元都在平衡位置附近作机械振动。 由于振动具有时间上的周期性, 所以波也具有时间上的周期性, 即每隔一定的时间,媒质中各质元的 振动状态都将复原。 媒质中振动状态复原时所需的最短时间, 也即质元完成一次全振动的时间叫波的周期, 周期的倒数叫频率。
虽然各类波的本质不同,各有其特殊的性质和规律, 但在形式上它们也具有许多共同的特征。 如都具有一定的传播速度,都伴随着能量的传播, 都能产生反射、折射、干涉或衍射等现象。
§2.1 行波
一. 机械波的产生 1. 机械波产生的条件
振源 作机械振动的物体——波源 媒质 传播机械振动的物体 在物体内部传播的机械波,是靠物体的弹性形成的, 因此这样的媒质又称弹性媒质。
什么是物质的弹性?
2.3 物体的弹性变形
物体包括固体、液体和气体,在受到外力作用时, 形状或体积都会发生或大或小的变化。 这种变化统称为形变
当外力不太大因而引起的形变也不太大时, 去掉外力,形状或体积仍能复原。 这个外力的限度称作弹性限度。
在弹性限度内,外力和形变具有简单的关系, 由于 外力施加的方式不同,形变可以有以下 几种基本方式: 线变 切变 体变
什么是物质的弹性?
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
f ( x0 , y0 , z0 , t )
其中 f ( x0 , y0 , z0 , t ) 为已知函数。
35
第二类边界条件(Neuman 边界条件):
规定所研究物理量在边界外法线方向 n 上的
方向导数的数值.
u f n
u f ( x0 , y0 , z0 ) , n ( x0 , y0 , z0 )
36
第三类边界条件(混合边界条件 也叫 Robin边界条件 ):规定所研究物理量及其
外法向导数的线性组合在边界上的值
u Hun
( x0 , y0 , z0 )
f ( x0 , y0 , z0 , t )
u f H :常系数 u n
37
以上三类边界条件当 f 0 时,分别称为 第一、二、三类齐次边界条件。
22
应用微积分中值定理:
ux ( x dx, t ) T1ux ( x, t ) uxx dx
dy f ' ( x)dx
ux ( x dx, t ) T1ux ( x, t ) uxx dx
x Fdx dxutt
Tuxx dx Fdx dxutt
39
2 a u tt u xx 0 在这一点无意义.如果,将
l 分成 x x0 ,x x0 两段分别考虑,
在各段上,弦振动方程有意义,但它是一 根弦的两段,并不是各自振动的。从数学
上来讲,不可能在两端上分别列出定解问
题。两段可作为一个整体来研究,两段的 振动是相互关联的。
40
u
F(0,t)
15
即整根弦由相互牵连的质点组成,对每个
质点即每个小段可应用 F ma . 方法:将连续分布的介质离散化为多质点
系统,再取内部任一代表性的点进行研究。将
弦细分为许多极小的小段,取区间上 ( x, x dx) 小段为代表。无质量且柔软,故该段仅受到相
邻两段的拉力 T1 和 T2 .
16
④ 对弦的每一小段dx,沿x方向(纵向) 没有运动,沿 x方向所受合外力为零。任一
表征物理量的选择常常是建立一个新 方程的起点。 (一个或几个)。
5
3.寻找(猜测)物理过程所遵守的 物理定律或物理公理;
4.写出物理定律的表达式,即数学
模型。
6
1.1 弦振动方程
一、弦的横振动方程 二、定解条件的提出
三、三类定解问题
7
一、 弦的横振动方程(均匀弦的微小横振动) 演奏弦乐(二胡,提琴)的人用弓在弦上来回 拉动,弓所接触的是弦的很小的一段,似乎只能引 起这个小段的振动,实际上振动总是传播到整个弦,
29
2、初始条件 在求解含时间t变量的数理方程时,往往要追
溯到早些某个所谓“初始”时间的状况(“历
史” ),于是称物理过程初始状况的数学表达 式为初始条件。
30
如弦振动方程:
utt a uxx 0
2
其初始条件为: 同一时刻( t 0 )情况
u ut
t 0 t 0
( x) ( x)
α1
α2
xLeabharlann 41x x0 虽是折点,但它们连续,即
u( x0 0, t ) u( x0 0, t )
在 x0 ,力 F (t ) 应和张力平衡,即
F (t ) T sin 1 T sin 2 0
①
sin 1 tan 1 ux ( x0 0, t )
sin 2 tan 2 ux ( x0 0, t )
Tux ( x0 0, t ) Tux ( x0 0, t ) F (t )
②
①、②合称为衔接条件,这时振动问题适定。
42
再如,不同材料组成的杆的振动,在
衔接处的位移和能量相等,即:
u1
E1u1 x
x x0
u2
E1u2 x
x x0
x x0
x x0
u1 ( x, t ),u2 ( x, t ) :杆的两部分位移.
动。二维波动方程,如薄膜的横振动方程,管
道中小振动的传播,理想传输线的电报方程等
均可用上述波动方程描述。故称为一类方程,
即波动方程。(也是称其为泛定方程的远大) 可描述一类物理现象。流体力学与声学中推导 三维波动方程,这里不再一一推导。
27
二、定解条件的提出
1、必要性。导出方程后,就得对方程进行 求解。但是只有泛定方程不足以完全确定方程 的解,即不足以完全确定具体的物理过程,因 为具体的物理过程还与其初始状态及边界所受 的外界作用有关,因而必须找一些补充条件, 用以确定该物理过程。
E1 , E2 :两部分的杨氏模量.
43
静电场中,两种电介质的交界面
s
上电势应相等(连续),电位移矢量的法
向分量也应相等(连续),其衔接条件是:
u1 s u1 s u1 1 n u2 2 n
s
s
44
其中 u1 , u 2 代表两种电介质的电势,
1 , 2 代表两种电介质的介电常数,(设电
utt Tuxx F
23
即
utt a2uxx f
—— 弦的强迫横振动方程
T
其中: a
2
f ,
F
量纲分析:T : MLT 2 , : ML1
24
∴
2
T
MLT 2 2 2 : L T 1 ML
2
即
a :LT
2
a
:振动的传播速度 a
T
它与弦的张力的平方根成正比,与弦的
38
4、其它条件
⑴ 衔接条件 由于一些原因,在所研究的区域里出 现跃变点,泛定方程在该点失去意义。如 波动方程(弦),如果有横向力F (t )集中地 作用于 x0 点, 这就成了弦的折点。在点 x0 斜率 u x 的左极限 ux ( x0 0, t ) 不同于右极限
ux ( x0 0, t ),因而 u xx不存在,
13
3. 研究建立方程 ① 如图,选弦绷紧时(不振动)直线为 x 轴
u
F
1
2
T2
T1
s
0
A
x x x
B
x
14
② 弦离开平衡位置的位移记为 为表征物理量。
u ( x, t ),
③因弦的振动是机械振动,基本规律为:
F ma, 然而弦不是质点,故 F ma
对整根弦并不适用。但整根弦可以细分为许 多极小的小段,每个小段可以抽象为质点。
10
由于张力的作用,一个小段的振动必带动它
的邻段,邻段又带动它自己的邻段,这样一个 小段的振动必然传播到整个弦,这种振动的传 播现象叫作波。弦是轻质弦(其质量只有张力 的几万分之一)。跟张力相比,弦的质量完全
可以略去。
11
① 模型实际上就是:柔软轻质细弦(“没 有质量”的弦) ② 将无质量的弦紧绷,不振动时是一根直
∵
u ux tan x dx tan 1 ux x dx
20
小振动近似: x dx 与 x 两点间任一时刻横
向位移之差 u( x dx, t ) u( x, t ) 与 dx 相比是一 个小量,即
u 1 x
于是①、②化简为:
位移矢量分别为 D1 , D2 , 则
D1n s D2n
s
D E u
45
⑵ 自然边界条件
某些情况下,出于物理上的合理性等原因, 要求解为单值、有限,就提出自然边界条件, 这些条件通常都不是要研究的问题直接给出, 而是根据解的特性要求自然加上去,故称为自 然边界条件,如:
33
3、边界条件
求解方程时还需考虑边界状况(周边“环 境”)(边界状况将通过逐点影响所讨论的 整个区域),称物理过程边界状况的表达式 为边界条件,或称为边值条件。
边界条件在数学上分为三类:
34
第一类边界条件(Dirichlet边界条件):直 接规定所研究的物理量在边界上的数值
u f u ( x, y, z, t )
1
19
cos 2 1
22
2!
24
4!
1
sin 1 1
13
3!
15
5!
1 tan 1
sin 2 2
23
3!
25
5!
2 tan 2
ds (dx)2 (du )2 1 (ux ) 2 dx
线,取为 x 轴。 ③ 将弦上个点的横向位移记为 u u( x, t )
12
④ 已知:线密度 ( x, t ) (t ), 重量不计,
张力 T ( x, t ) 沿切线方向,不随x变化,弦中
各点的张力相等(小振动下T 与t 也无关). ⑤ 研究方法:连续介质,微积分思想, 任意性。
28
从物理角度看:泛定方程仅表示一般性(共
性),要为物体的运动个性化附加条件。 从数学角度看:微分方程解的任意性也需附 加条件。通解中含任意函数(解不能唯一确定 )。通过附加条件确定任意函数(常数),从 而确定解。这些附加条件就是前面所谈的问题 的“历史”与“环境”,即初始条件和边界条 件,统称为定解条件。
弦的各处都振动起来。振动如何传播呢?
8
1. 物理模型 实际问题:设有一根细长而柔软的弦,紧绷 于A,B两点之间,在平衡位置附近产生振幅极 为微小的横振动(以某种方式激发,在同一平
面内,弦上各点的振动方向相互平行,且与波
的传播方向(弦的长度方向)垂直),求弦上 各点的运动规律。