分离变量法(驻波法)求一维有界区域的自由振动问题(第二类齐次边界条件):
utt(x,t)a2uxx(x,t) (0xl)uxx00 , uxl0ut0(x) , utt0(x)
其中(x),(x),为已知函数.
解:依题意设定解问题的解为 u(x,t)X(x)T(t)
其中X(x)只是x的函数,与t无关,T(t)只是t的函数,与x无关.代入方程得
T''(t)X''(x)a2T(t)X(x)
T''(t)X(x)a2X''(x)T(t)等式两边分别是t和x两个独立变量的函数,要使它们相等,必须都等于一个常数,设为
T''(t)X''(x)2aT(t)X(x)
即
则得到两个常微分方程
T''(t)a2T(t)0 (1)X''(x)X(x)0 (2)
1
由边界条件有X'(0)0 X(l)0 (3)
当0或0方程(2)都只有平凡解u(x,t)0
当0时(2)的解为X(x)C1cosxC2sinx (4) (C1,C2为任意常数)
将(4)代入(3)的第一式得X'(0)C2=0C2=0,而C1不可能为0,再将(4)代入(3)的第二式得
(2n1)(2n1)22l (n1,2,3,)
X(l)C1cosl0cosl0l∴ 只能取相应的一系列值
(2n1)22n4l2 (n1,2,3,)---由定解问题决定的本征值
而
X(x)Xn(x)C1cos(2n1)x2l---由定解问题决定的本征函数
(2n1)22n4l2将代入(1),解得
'T(t)Tn(t)Ancos(2n1)a(2n1)a'tBnsint''An,Bn2l2l (为任意常数)
2
∴
'un(x,t)Xn(x)Tn(t)[Ancos(2n1)a(2n1)a(2n1)'tBnsint]C1cosx2l2l2l
[Ancos(2n1)a(2n1)a(2n1)tBnsint]cosx2l2l2l
即
u(x,t)un(x,t)[Ancosn1n1(2n1)a(2n1)a(2n1)tBnsint]cosx2l2l2l
而
ut(x,t)[Ann1(2n1)a(2n1)a(2n1)a(2n1)a(2n1)sintBncost]cosx2l2l2l2l2l
任意常数An,Bn的确定
由初始条件有
ut0(x)(x)Ancosn1(2n1)x (0xl)2l
An2l(2n1)()cosd (n1,2,3)l02l
(2n1)a(2n1)Bncosx (0xl)2l2ln1
utt0(x)(2n1)a2l(2n1)Bn()cosd2ll02l
3
l4(2n1)Bn()cosd (n1,2,3)(2n1)a0l
10.3分离变量法(驻波法)求一维有界区域的自由振动问题(第二类齐次边界条件):
utt(x,t)a2uxx(x,t) (0xl)ux00 , uxxl0ut0(x) , utt0(x)
其中(x),(x),为已知函数.
解:依题意设定解问题的解为 u(x,t)X(x)T(t)
其中X(x)只是x的函数,与t无关,T(t)只是t的函数,与x无关.代入方程得
T''(t)X''(x)a2T(t)X(x)
T''(t)X(x)a2X''(x)T(t)等式两边分别是t和x两个独立变量的函数,要使它们相等,必须都等于一个常数,设为
T''(t)X''(x)a2T(t)X(x)
即
则得到两个常微分方程
4
T''(t)a2T(t)0 (1)X''(x)X(x)0 (2)
由边界条件有X(0)0 X'(l)0 (3)
当0或0方程(2)都只有平凡解u(x,t)0
当0时(2)的解为X(x)C1cosxC2sinx (4) (C1,C2为任意常数)
将(4)代入(3)的第一式得X(0)C1=0,而C2不可能为0,再将(4)代入(3)的第二式得
(2n1)(2n1)22l (n1,2,3,)
X'(l)C2cosl=0cosl0l∴ 只能取相应的一系列值
(2n1)22n4l2 (n1,2,3,)---由定解问题决定的本征值
而
X(x)Xn(x)C2sin(2n1)x2l---由定解问题决定的本征函数
(2n1)22n4l2将代入(1),解得
5
'T(t)Tn(t)Ancos(2n1)a(2n1)a'tBnsint''An,Bn2l2l (为任意常数)
∴
'un(x,t)Xn(x)Tn(t)[C2Ancos(2n1)a(2n1)a(2n1)'tC2Bnsint]sinx2l2l2l
[Ancos(2n1)a(2n1)a(2n1)tBnsint]sinx2l2l2l
即
u(x,t)un(x,t)[Ancosn1n1(2n1)a(2n1)a(2n1)tBnsint]sinx2l2l2l
而
ut(x,t)[Ann1(2n1)a(2n1)a(2n1)a(2n1)a(2n1)sintBncost]sinx2l2l2l2l2l
任意常数An,Bn的确定
由初始条件有
ut0(x)(x)Ansinn1(2n1)x (0xl)2l
An2l(2n1)()sind (n1,2,3)0l2l
6
utt0(x)(2n1)a(2n1)Bnsinx (0xl)2l2ln1
(2n1)a2l(2n1)Bn()sind02ll2l
l4(2n1)()sind (n1,2,3)0(2n1)a2l
Bn10.3分离变量法(驻波法)求一维有界区域的热传导问题(第一类齐次边界条件):
ut(x,t)a2uxx(x,t) (0xl)ux00 , uxl0ut0(x)
其中(x)为已知函数.
解:依题意设定解问题的解为 u(x,t)X(x)T(t)
其中X(x)只是x的函数,与t无关,T(t)只是t的函数,与x无关.代入方程得
T'(t)X''(x)a2T(t)X(x)
T'(t)X(x)a2X''(x)T(t)等式两边分别是t和x两个独立变量的函数,要使它们相等,必须都等于一个常数,设为
7
即
T'(t)X''(x)a2T(t)X(x)
则得到两个常微分方程
T'(t)a2T(t)0 (1)X''(x)X(x)0 (2)
由边界条件有X(0)0 X(l)0 (3)
当0或0方程(2)都只有平凡解u(x,t)0
当0时(2)的解为X(x)C1cosxC2sinx (4) (C1,C2为任意常数)
将(4)代入(3)的第一式得X(0)C1=0,而C2不可能为0,再将(4)代入(3)的第二式得
nl (n1,2,3,)
X(l)C2sinl0sinl0ln∴ 只能取相应的一系列值
n2)l (n1,2,3,)---由定解问题决定的本征值
n( 8
而
X(x)Xn(x)C2sinnxl---由定解问题决定的本征函数
将
n22nl2代入(1),解得 n22a2T(t)TC'n(t)nel2t (
C'n为任意常数)
22a2un(x,t)X'ntn(x)Tn(t)C2Cnel2sinnlxCn22a2nel2t∴
sinnlx
n22a2l2t即
u(x,t)un(x,t)Cnesinnn1n1lx
任意常数Cn的确定
由初始条件有
unt0(x)(x)Cnsinlx (0xl)n1C2lnnl0()sinld (n1,2,3)
9
因篇幅问题不能全部显示,请点此查看更多更全内容