波动方程
- 格式:doc
- 大小:542.84 KB
- 文档页数:13
1.1 波动方程的形式一维波动方程(描述弦的振动或波动现象的)()t x f x u a t u ,22222=∂∂-∂∂ 二维波动方程(例如薄膜振动)()t y x f y u x u a t u ,,2222222+⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=∂∂ 三维波动方程(例如电磁波、声波的传播)()t z y x f z u y u xu a t u ,,,222222222+⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂=∂∂ 1.2 波动方程的定解条件(以一维波动方程为例)(1)边界条件 ①第一类边界条件(又称Dirichlet 边界条件):弦振动问题中,弦的两端被固定在0=x 及l x =两点,因此有()0,0=t u ,()0,=t l u 。
②第二类边界条件(又称Neumann 边界条件):弦的一端(例如0=x )处于自由状态,即可以在垂直于x 轴的直线上自由滑动,未受到垂直方向的外力,此时成立0=∂∂=ox xu。
也可以考虑更普遍的边界条件()t xu x μ=∂∂=0,其中()t μ是t 的已知函数。
③第三类边界条件:弦的一端固定在弹性支承上,不放考虑在l x =的一端,此时边界条件归结为0u =⎪⎭⎫ ⎝⎛+∂∂=l x u x σ。
也可以考虑更普遍的情况()t u x lx v u =⎪⎭⎫⎝⎛+∂∂=σ,其中()t v 是t 的已知函数。
1.3 利用叠加原理求解初值问题 初值问题()()()()⎪⎪⎩⎪⎪⎨⎧+∞<<∞=∂∂==+∞<<∞>=∂∂-∂∂)x -(,,:0t x 0,-t ,,22222x t u x u t x f x u a t u ψϕ (1) 利用叠加原理求解上述初值问题,叠加原理表明由()t x f ,所代表的外力因素和由()()x x ψϕ,所代表的初始振动状态对整个振动过程所产生的综合影响,可以分解为单独只考虑外力因素或只考虑初始振动状态对振动过程所产生的影响的叠加。
即如果函数()t x u ,1和()t x u ,2分别是下述初值问题(I )()()()()⎪⎪⎩⎪⎪⎨⎧=∂∂===∂∂-∂∂2.1.....................,:0t 1.1. (022)222x t u x u x u a t u ψϕ和 (II )()()()⎪⎪⎩⎪⎪⎨⎧=∂∂===∂∂-∂∂4.1....................................................................0,0:0t 3.1................................................................,22222t u u t x f x u a t u的解,那么()()t x u t x u u ,,21+=就一定是定解问题(1)的解。
首先,我们考察自由振动情况的初值问题(I )。
定理1.1 设()()()(),,12R C x R C x ∈∈ψϕ那么问题(I )存在唯一的解()t x u ,,它由公式()()()()da a a at x at x t x u atx atx ⎰+-+++-=ψϕϕ212,给出。
证:可以通过自变量变换的方法求解。
引入新变量.,at x at x +=-=ηξ利用复合函数求道的法则,得到ηξηηξξ∂∂+∂∂=∂∂∂∂+∂∂∂∂=∂∂uu x u x u x u , 22222222ηηξξηηξξ∂∂+∂∂∂+∂∂=∂∂⎪⎭⎫ ⎝⎛∂∂∂∂+∂∂⎪⎭⎫ ⎝⎛∂∂∂∂=∂∂uu u x x u x x u x u , ⎪⎪⎭⎫ ⎝⎛∂∂+∂∂∂-∂∂=∂∂⎪⎪⎭⎫ ⎝⎛∂∂-∂∂=∂∂222222222,ηηξξξηu u u a t u u u a t u , 从而ηξ∂∂∂=∂∂-∂∂u a x u a t u 22222224- 由于02>a ,方程(1.1) 就转化为)5.1......(..............................02=∂∂∂ηξu方程(1.5)可以直接求解。
把它关于η积分一次,再关于ξ积分一次,得出其通解为()()())6.1..(....................,ηξG F t x u +=其中F 和G 是任意两个可微分的单变量函数。
再代会原来的自变量,方程(1.1)的通解可以表示为()()())7.1.......(..........,at x G at x F t x u ++-=利用这个通解表达式,就可以由初始条件(1.2)来决定函数F 和G ,从而求出初值问题(I )的解。
把(1.7)代入初始条件(1.2),求得()()()()()()())9.1.........()8.1...(..........x x G x F a x x G x F ψϕ='+'-=+再将(1.9)两边积分,得()()()())10.1....(. (x)a x d a C x G x F a ⎰=++-ψ其中0x 是任意一点,而C 是积分常数。
由(1.8)和(1.10),就可以解出F 和G :()()()()()()⎪⎪⎩⎪⎪⎨⎧-+=+-=⎰⎰xx x x a C da a a x x G a C da a a x x F 002212122121ψϕψϕ (1.11) 把他们代入(1.7)中,就得到初值问题(I )的解()()()())12.1( (21)2,da a a at x at x t x u atx atx ⎰+-+++-=ψϕϕ 如果初值问题(I )有解,则解一定由初始条件由公式(1.12)表示出来,因此解一定是唯一的。
再次,我们来考虑问题(II )的解。
定理1.2 齐次化原理(Duhamel 原理) 若);,(τt x W 是初值问题()⎪⎪⎩⎪⎪⎨⎧=∂∂==>=∂∂-∂∂τττ,0)(,022222x f t W W t t x W a t W ,:的解(其中τ为参数),则()()ττd t x W t x u t⎰=0;,,就是初值问题(II )的解。
证:自由项),(t x f 表示时刻t 时在位置x 处单位质量所受的外力,而tu∂∂表示速度。
把时间段[0,t]分成若干小的时间段)....,2,1(1l j t t t j j j =-=∆+,在每个小的时间段上j t ∆中,()t x f ,可以看做与t 无关,从而以()j t x f ,来表示。
由于()(),,,ρj j t x F t x f =而()j t x F ,表示外力,所以在时间段j t ∆中自由项所产生的速度改变量为()j j t t x f ∆,。
把这个速度改变量看作是在时刻j t t =时的初始速度,它所产生的振动可以由下面的齐次方程带非齐次初始条件的初值问题来描述:())13.1.(..........,~0~)(,0~~22222⎪⎪⎩⎪⎪⎨⎧∆=∂∂==>=∂∂-∂∂j j j j t t x f t W W t t t t x W a t W ,: 其解记为),;,(~j j t t t x W ∆。
按照叠加原理,自由项()t x f ,所产生的总效果可以看成是无数个这种瞬时作用所产生的效果的叠加。
这样,问题(II )的),(t x u 应表示成)14.1(..........).,;,(~lim ),(1∑=→∆∆=lj j j t t t t x W t x u j由于(1.13)为线性方程,所以W ~与j t ∆成正比,故如果记);,(τt x W 为齐次方程的定解问题())15.1....(..........,0)(,022222⎪⎪⎩⎪⎪⎨⎧=∂∂==>=∂∂-∂∂τττx f t W W t t x W a t W ,: 的解,则);,(),;,(~j j j j t t x W t t t t x W ∆=∆。
于是定解问题(II )的解可表示为∑∑⎰==→∆→∆=∆=∆=lj lj tjjt j j t d t x W t t t x W t t t x W t x u j j 11.);,();,(lim),;,(~lim ),(ττ为了写出);,(τt x W 的具体表达式,在初值问题(1.15)中作变换.τ-='t t 相应地,(1.15)化为())16.1....(. (00)0(,022222⎪⎪⎩⎪⎪⎨⎧='∂∂=='>'=∂∂-'∂∂τx f t W W t t x W a t W ,: 的形式,于是利用问题(I )解的形式,就得到其解为)17.1(..........),(21),(21);,()()(ξτξξτξτττd f a d f a t x W t a x t a x t a x t a x ⎰⎰-+--'+'-== 由此得到所考察的初值问题(II )的解为)18.1.........(),(21),(21),(0)()(τξτξτξτξττd d f a d d f a t x u Gt t a x t a x ⎰⎰⎰⎰==-+-- 其中区域G 为()τξ,平面上过点(x,t )向下作两特征线与ξ轴所夹的三角形区域。
下面验证(1.18)是否为初值问题(II )的解。
假设1C f ∈,从(1.18)式,关于含参变量积分的求导法则,可得()()[]()()[]()()[]()()[],),(),(21,),(),(21,),(),(2),(,),(),(212200220ττττττττττττττττττττd t a x f t a x f a x u d t a x f t a x f a x u d t a x f t a x f at x f t u d t a x f t a x f t u tx x ttx x t⎰⎰⎰⎰----+=∂∂----+=∂∂----++=∂∂--+-+=∂∂于是有 ()t x f xu a t u ,22222=∂∂-∂∂,即),(t x u 满足方程(1.3),再由(1.18)式及t u ∂∂的表示式可得 ()0,0,00=∂∂===t t t ut x u ,即),(t x u 满足初始条件(1.4).所以(1.18)式表示的是定解问题(II )的解。
将(1.12)式和(1.18)式叠加就得到问题(1)的解。
1.4 求解初边值问题 初边值问题()()()⎪⎪⎪⎩⎪⎪⎪⎨⎧=====∂∂===∂∂-∂∂)22.1........(..............................0:)21.1(........................................0:0)20.1.......(..........,:0t )19.1.........(..........,22222u l x u x x t u x u t x f x u a t u ψϕ ()0,x 0≥≤≤t l 利用叠加原理,上述初边值问题可以分解为下面两个初边值问题:(III )()()⎪⎪⎪⎩⎪⎪⎪⎨⎧====∂∂===∂∂-∂∂)25.1....(....................0:0)24.1.(..........,:0t )23.1..(. (011121)22212u l x x x t u x u x u a t u 及ψϕ(IV )()⎪⎪⎪⎩⎪⎪⎪⎨⎧====∂∂===∂∂-∂∂)28.1.(..............................00)27.1...(....................0,0:0t )26.1..(..........,222222222u l x x t u u t x f x u a t u :及那么初边值问题的解就是21u u u +=。