板的振动

2 r2
1 r
r
1 r2
2
2
2
W
0
得常微分方程
d2 F d r2
1 r
dF dr
n2 r2
F
0
或引用无因次的变量x＝γr 而得
x2
d2 F d x2
x dF dx
x2
n2
F
0
这一微分方程的解答是
F C1Jn ( x) C2Nn ( x) C3In ( x) C4Kn ( x)
k 2 2
a2
2
k 2
a2
2
m D
k 2 2
a2
上述四个根成为±α及±iβ，而微分方程的解 可写为
Yk C1 chy C2 shy C3 cos y C4 sin y
从而得振形函数的表达式
W
(C1 chy
C2
shy
C3
cosy
C4
sin
y) sin
kx
a
在少数的情况下，γ2<k2π2/a2，而上面所示
D4 w
m
2w t 2
这就是薄板自由振动的微分方程。
现在来试求微分方程的如下形式的解答
w wk ( Ak coskt Bk sinkt)Wk ( x, y)
k 1
k 1
在这里，薄板上每一点(x，y)的挠度，被表示成为 无数多个简谐振动下的挠度相叠加，而每一个简谐 振动的频率是ωk ，另一方面，薄板在每一瞬时t的 挠度，则被表示成为无数多种振形下的挠度相叠加， 而每一种振形下的挠度是由振形函数Wk(x，y)表示 的。
4W 4W 0
现在就可能比较简便地求得W的满足边界条件的、 函数形式的非零解，从而求得相应的γ值，然后再 用式
2m 4
D
求出相应的频率。
将求出的那些振形函数及相应的频率取为Wk及 ωk，代入表达式
w wk ( Ak coskt Bk sinkt)Wk ( x, y)
k 1
k 1
因此，只有在特殊简单的情况下，才有可能 求得薄板自由振动的完整解答，即任一瞬时的挠 度。在绝大多数的情况下，只可能求得各种振形 的振形函数及相应的频率。但是，这也就可以解 决工程上的主要问题了。
第二节 四边简支的矩形薄板的
自由振动
当矩形薄板的四边均为简支边时，可以较简 单地得出自由振动的完整解答。取振形函数为
W sin kx sin ny
ab
其中k及n为整数，可以满足边界条件。代入振形 微分方程
4W 4W 0
得到
4
k2 a2
n2 b2
2
4
sin
kx
a
sin
ny
b
0
为了这一条件在薄板中面上的所有各点都能
满足，也就是在x和y取任意值时都能满足，必须
有
4
k a
2 2
n2 b2
2
D4we q
设薄板在振动过程中的任一瞬时t的挠度为we
＝we(x，y)，则薄板每单位面积上在该瞬时所受的
弹性力，将与横向荷载q及惯性力qi成平衡，即
D4wt q qi
注意薄板的加速度是 2wt
t 2
因而每单位面积上的惯性力
qi
m
2wt t 2
其中m为薄板每单位面积内的质量(包括薄板本身
的质量和随薄板振动的质量)，则前式可以改写为
4
0
得到
4
4
k2 a2
n2 b2
2
得出求自然频率的公式
2
D 4
m
4
k a
2 2
n2 b2
2
Dபைடு நூலகம்m
命k及n取不同的整数值，可以求得相应于不同振形 的自然频率
2
k2 a2
n2 b2
D m
当薄板以这一频率振动时，振形函数为
Wkn
sin
kx
a
sin
ny
b
而薄板的挠度为
w
(
Akn
c
osknt
薄板在平行于中面方向的所谓纵向振动，由 于它在工程实际中无关重要，而且在数学上也难 以处理，所以不加讨论。首先来讨论薄板的自由 振动。
单自由度振动的例子
薄板自由振动的一般问题：在一定的横向 荷载作用下处于平衡位置的薄板，受到干扰力 的作用而偏离这一位置，当干扰力被除去以后， 在该平衡位置附近作微幅振动。
第五章 薄板的振动问题
第一节 薄板的自由振动 第二节 四边简支板的 自由振动 第三节 两对边简支板的自由振动 第四节 圆形薄板的自由振动 第五节 用差分法求自然频率 第六节 用能量法求自然频率 第七节 薄板的受迫振动
第五章 薄板的振动问题
第一节 薄板的自由振动
关于薄板的振动问题，这里将只讨论薄板在 垂直于中面方向的所谓横向振动，因为这是工程 实际中的重要问题。
m
2w t 2
仍然得出振形微分方程
4W 2m W 0
D
或
4W 4W 0
可以改写为
(2 2 )(2 2 )W 0
或
2 r2
1 r
r
1 r2
2
2
2
W
0
取振形函数为如下的形式：
W F(r) cosn
其中n=0，1，2，…。相应于n＝0，振形是轴对 称的。相应于n＝1, 2；圆板的环向围线将分别 具有一个及两个波，板的中面将分别具有一根 或两根径向节线，余类推。将上式代入式
w y
yb
0
将W的表达式代入(γ2>k2π2/a2)，得到Cl至C4的齐次 线性方程组，
C1 C2 0
2C1 2 0 chbC1 shbC2 cosb chbC3 sin b chbC4 0 shbC1 chbC2 sin bC3 cosbC4 0
命这一方程组的系数行列式等于零，展开以后，进 行一些简化，最后可得出
D4 w
m
2w t 2
现在，仍然把微分方程的解答取为无数多简谐振
动的叠加
w wk ( Ak coskt Bk sin kt)Wk (r, )
k 1
k 1
为了求出各种振形下的振形函数Wk，以及与之 相应的频率ωk，我们取
w (Acost Bsint)W(r, )
代入自由振动微分方程
D4 w
d4 Yk d y4
2k 2 2
a2
d2 Yk d y2
k 4
a4
4
2
Yk
0
它的特征方程是
r4
2k 2 2
a2
r2
k 4
a4
4
2
0
而这个代数方程的四个根是
m
D
k 2 2
a2
2
k 2 2
a2
2
在大多数的情况下，γ2>k2π2/a2，而上面所示的 四个根是两实两虚，取正实数
2
k 2
a2
2
m D
D4wt
q
qi
q
m
2wt t 2
将上式与下式相减
D4we q
得到
D4 (wt
we
)
m
2 wt t 2
由于we不随时间改变，所以上式可以改写成为
D4
(wt
we
)
m
2 t 2
(wt
we
)
在以下的分析中，为了简便，我们把薄板的挠
度不从平面位置起，而从平衡位置量起。于是薄板
在任一瞬时的挠度为w＝wt－we，而上式成为
相应于薄板的任何振动，振形函数W必须具 有某一个非零解，因而系数Cl至C4不能都等零。 于是可以命上述齐次线性方程组的系数行列式等 于零，从而得到一个计算自然频率的方程。
例如，设y＝0的一边为简支边，而y＝b的一边为
固支边，则有如下的四个边界条件：
(w)y0 0
2w y 2
y0
0
(w)yb 0
的四个根都是实根，取正实数
k 2 2
a2
2
k 2 2
a2
m D
'
k 2 2
a2
2
k 2 2
a2
m D
从而得振形函数的表达式
W
(C1 chy
C2
shy
C3 y ch '
y
C4 y sh '
y ) sin
kx
a
振型函数应满足边界条件。
不论在哪一种情况下，都可由y＝0及y＝b处 的四个边界条件得出Cl至C4的一组四个齐次线性 方程。
(1)试求薄板振动的频率，特别是最低频 率。
(2)设已知薄板的初始条件，即已知初挠 度及初速度，试求薄板在任一瞬时的挠度。
当然，如果求得薄板在任一瞬时的挠度， 就易求得薄板在该瞬时的内力。
设薄板在平衡位置的挠度为we＝we(x，y)，这
时，薄板所受的横向静荷为q＝q(x，y)。按照薄板 的弹性曲面微分方程，我们有：
为了求出各种振形下的振形函数Wk，以及与之 相应的频率ωk，我们取
w (Acost Bsint)W(x, y)
代入自由振动微分方程
D4 w
m
2w t 2
然后消去因子(Acosωt十Bsinωt)，得出所谓振形微
分方程
4W 2m W 0
D
如果可以由这一微分方程求得W的满足边界条件的
非零解，即可由相应的关系式(对任意的一点(x,y)
如果薄板具有圆孔，则在外边界及孔边各有两个 边界条件。利用这四个边界条件，可得出的一组 Cl至C4四个齐次线性方程。命这一方程组的系数 行列式等于零，可以得出计算频率的方程，从而 求得各阶的自然频率。
如果薄板无孔，则在薄板的中心(x＝γr＝0)， Nn(x)及Kn(x)成为无限大。为了使W不致成为无限 大，须在式中取C2＝0，C4＝0。于是式简化为
这样进行计算，虽然可以求得自然频率的精 确值，但代数运算和数值计算都是比较繁的。因 此，在工程实践中计算矩形板的自振频率，特别 是最低自然频率，不论边界条件如何，都宜用差 分法或能量法。
第四节 圆形薄板的自由振动
对于圆形薄板的自由振动，也可以与上相同
地进行分析。在极坐标中，薄板的自由振动的微
分方程仍然是
kn
s in knt )
sin
kx
a
sin
ny
b
当矩形薄板的四边均为简支边时，可以较简 单地得出自由振动的完整解答。
第三节 两对边简支的矩形薄板的
自由振动
取振形函数为
x
W
Yk
sin
kx
a
其中Yk是待定的y的函数。W可 以满足该两简支边的边界条件。
将其代入振形微分方程
y
4W 4W 0
得出常微分方程
就有可能利用初始条件求得该表达式中的系数Am
及Bm 。 设初始条件为
(w)t0 w0( x, y)
w t
t0
v0( x,
y)
则由上式得
AkWk (x, y) w0(x, y)
k 1
BkkWk (x, y) v0(x, y)
k 1
于是可见，为了求得Am及Bm，须将已知的 初挠度w0及初速度v0展为Wm的级数，这在数学处 理上是比较困难的。
thb th b 0 b b
上列方程可以改写为
th b 2 m2 2 / a2 th b 2 m2 2 / a2
0
b 2 m2 2 / a2
b 2 m2 2 / a2
求得γ2的实根，即可求得自然频率ω
2 D
m
用如上方法求得的最低自然频率，可以表示 成为依赖于边长比值a／b算得的系数k值，并以表 来表示。
其中Jn(x)及Nn(x)分别为实宗量的、n阶的第一种 及第二种贝塞尔函数，In(x)及kn(x)分别为虚宗量 的、n阶的第一种及第二种贝塞尔函数（又称修 正贝塞尔函数）。
贝塞尔函数
将上式代入
W F(r) cosn
即得振形函数如下：
W (C1Jn (x) C2Nn (x) C3In (x) C4Kn (x)) cos n
Bkn
sin
knt
)
sin
kx
a
s
in
ny
b
当k＝n＝1时，得到薄板的最低自然频率
m in
2
k2 a2
n2 b2
D m
2
1 a2
1 b2
D m
与此相应，薄板振动的振形函数为
W11
sin
x
a
sin
y
b
而薄板在x方向和y方向都只有一个正半弦波。最 大挠度发生在薄板的中央(x＝a／2，y＝b／2)。
4 ab
a 0
b
kx ny
0 w0 sin a sin b d x d y
Dkn
4 ab
a 0
b 0
v0
sin
kx
a
sin
ny
b
d
x
d
y
根据初始条件为
(w)t0 w0( x, y)
可得
w t
t0
v0( x,
y)
Akn Ckn
Bkn
Dkn
kn
w
k 1 n1
(Ckn
cosknt
Dkn
W (C1Jn (x) C3In ( x)) cos n
由板边的两个边界条件，可以得出 Cl 及 C3 的 一组两个齐次线性方程，命方程组的系数行列式 等于零，也就得出计算自然频率的方程。
参见习题15－3。
第五节 用差分法求自然频率
当k＝2而n＝1时，自然频率为
21
2
4 a2
1 b2
D m
x
相应的振形函数为
W21
sin
2x
a
sin
y
b
y
薄板在x方向有两个正弦半波，而在y方向只有一个 正弦半波。对称轴x＝a／2是一根节线(挠度为零的 线，亦即在薄板振动时保持静止的线）。振形如图 所示，图中的有阴线部分及空白部分表示相反方向 的挠度。
都成立）
2 D 4W
mW
求得相应的频率ω。自由振动的频率，称为自然频 率或固有频率，完全决定于薄板的固有特性，而 与外来因素无关。
实际上，只有当薄板每单位面积内的振动质
量为常量时，才有可能求得函数形式的解答。这 时，命
2m 4
D
则振形微分方程
4W 2m W 0
D
简化为常系数微分方程
薄板的总挠度为
w
k 1 n1
(
Akn
cosknt
Bkn
sin
k nt )
sin
kx
a
sin
ny
b
为了求得Am及Bm，须将已知的初挠度w0及初速 度v0展为Wm的级数
w0
k 1 n1
Ckn
sin
kx
a
sin
ny
b
v0
k 1 n1
kx ny
Dkn sin a sin b
Ckn