跳转至

线性振动

线性振动是物理学中一种重要的运动形式,广泛存在于各种物理系统中.它通常指系统在平衡位置附近的小幅度振动,其运动方程可以用线性微分方程描述.

简单谐振动

如果某一物量在平衡点附近往复振动,且随时间的变化是余弦(或正弦)的函数形式,我们就把这种振动称为简谐振动.

运动方程推导

简单谐振动是最基本的线性振动形式,其运动方程为:

d2xdt2+ω2x=0

这一方程可以通过牛顿第二定律推导得到.例如,对于一个质量为 m 的物体,受到的回复力 F 满足胡克定律:

F=kx

根据牛顿第二定律:

F=md2xdt2

代入回复力表达式:

md2xdt2=kx

整理得:

d2xdt2+kmx=0

ω=km,即得简单谐振动的运动方程.

解的形式

该方程的通解为:

x(t)=Acos(ωt+ϕ)

其中:

  • A 为振幅,振子离开平衡位置的最大位移的绝对值,由初始条件决定.是求解上述微分方程时引入的常量;
  • ϕ 为初相位,初始时刻振子的相位,由初始条件决定.是求解上述微分方程时引入的常量;
  • ω 为角频率(圆频率),在简谐振动的几何表示中,表示旋转矢量的角速度;
  • x 为位移,表示振子相对于平衡位置的偏移量.
例 1

设一个质量为 m 的物体连接在弹簧上,弹簧劲度系数为 k.求该系统的角频率 ω 及其运动方程.

解:
易得运动方程为:

d2xdt2+kmx=0

角频率:

ω=km

因此,运动方程为:

x(t)=Acos(kmt+ϕ)

其中 Aϕ 由初始条件决定.

例 2

设一个质量为 m 的物体连接在摆绳上,摆长为 l.求该单摆系统的角频率 ω 及其运动方程.

解:

M=mglsinθ=Iβ=ml2d2θdt2

对于小角度摆动,sinθθ,单摆的运动方程为:

d2θdt2+glθ=0

角频率:

ω=gl

因此,运动方程为:

θ(t)=Acos(glt+ϕ)

其中 Aϕ 由初始条件决定.

小角度近似

在单摆的例子中,我们使用了小角度近似 sinθθ.这一近似仅在摆动角度较小时成立(通常小于约 10 度).对于较大的摆动角度,运动方程将变得非线性,不能再用简单谐振动的形式描述.

简谐振动的判断

并非所有往复振动都是简谐振动.简谐振动要求回复力与位移成正比且方向相反,即满足 胡克定律的形式.对于非线性系统或大幅度振动,运动方程可能不再是线性的,不能用简谐振动的形式描述.我们往往也可以通过得到相应的 特征微分方程 来判断振动是否为简谐振动.

周期与频率

振动的周期 T 和频率 f 分别为:

T=2πω,f=1T=ω2π.

相位

振动的相位 θ 定义为:

θ=ωt+ϕ.

相位表示振动在某一时刻所处的位置状态.规定 0<θ<2π

对于两个同频简谐振动,其相位差 Δθ 为:

Δθ=θ1θ2=ϕ1ϕ2.

Δθ=0,则两振动 同相,若 Δθ=π,则两振动 反相. 若 0<Δθ<π,则称 x1x2 超前,若 π<Δθ<2π,则则称 x1x2 落后

x, v, a 三者的相位关系

x(t)=Acos(ωt+ϕ)v(t)=Aωcos(ωt+ϕ+π2)=Aωsin(ωt+ϕ)a(t)=Aω2cos(ωt+ϕ+π)=Aω2cos(ωt+ϕ)

由此可见,速度 v 比位移 x 超前 π2,加速度 a 与位移 x 反相.你可以想想为什么会有这样的相位关系.

能量分析

动能与势能

简单谐振动的能量包括动能和势能,其表达式分别为:

  • 动能:
Ek=12mv2=12mω2A2sin2(ωt+ϕ)

在一个周期内的平均值:

Ek=1T0TEkdt=14kA2
  • 势能:
Ep=12kx2=12kA2cos2(ωt+ϕ)

在一个周期内的平均值:

Ep=1T0TEpdt=14kA2

总能量守恒

总能量为动能与势能之和:

E=Ek+Ep=12kA2=12mω2A2

总能量在振动过程中保持不变,动能和势能在振动过程中相互转化.

简谐振动的几何表示

简谐振动可以通过旋转矢量的投影来几何表示.设一个长度为 A 的旋转矢量以角速度 ω 绕原点逆时针旋转,其在水平轴上的投影即为简谐振动的位移 x(t)

简正模

在多自由度系统中,线性振动可以分解为若干个独立的简正模(Normal Modes).每个简正模对应一个特定的频率和振动模式,系统的总振动可以看作这些简正模的叠加.

用复数表示简谐振动

简谐振动

x(t)=Acos(ωt+φ)

也可以用复数的实部和虚部表示:

x~(t)=Aei(ωt+φ) 

其中,A 为振幅,ω 为角频率,φ 为初相位.

上式的右端又可以写为 (Aeiφ)eiωt=A~eiωt,其中 Aeiφ=A~ 是一个复常数,表示振动的初始状态,称为 复振幅

x~(t)=A~eiωt

如果我们取上式的实部,就得到了简谐振动的标准形式:

x(t)=Re(x~(t))=Re(A~eiωt)=Acos(ωt+φ)

x~(t) 表示简谐振动,则速度和加速度为

{v~(t)=dx~(t)dt=iωA~eiωt=iωx~(t)a~(t)=dv~(t)dt=d2x~(t)dt2=ω2A~eiωt=ω2x~(t)
例题

已知一个线性三原子分子 A2B 的纵向振动模型:质量为 m 的两个 A 原子位于两端,质量为 MB 原子位于中间;相邻原子间由劲度系数均为 k 的弹簧连接(ABBA).求简正模与对应角频率.

** 解:** 设三原子沿直线位移为 x1,x2,x3(分别对应 A,B,A),取平衡位置为零.
弹簧形变为 x1x2x3x2,于是运动方程为

{mx¨1=k(x1x2),Mx¨2=k(x2x1)k(x2x3)=k(2x2x1x3),mx¨3=k(x3x2).

取简正模形式 xj(t)=ajeiωt,得到代数方程

{mω2a1=k(a1a2),Mω2a2=k(2a2a1a3),mω2a3=k(a3a2).

利用对称性分两类:

  1. 整体平移模a1=a2=a3,此时弹簧不伸长,ω0=0

  2. 反对称伸缩模a1=a3,由对称性得 a2=0
    代入第一式得 mω2a1=ka1,因此

ω=km,(a1,a2,a3)(1,0,1).
  1. 对称伸缩模a1=a3=aa2=b
    由第一式 mω2a=k(ab),第二式 Mω2b=2k(ba)
    联立消去 b 可得
ω+2=km+2kM,(a1,a2,a3)(1,2mM,1).

因此该分子有一个零频平移模与两个非零振动模 ω,ω+

例题:受迫振动中「功率共振」的频率

对线性受迫振动

x¨+2γx˙+ω02x=F0mcos(Ωt),

证明稳态下平均耗散功率 P 取最大值时的驱动频率为 Ω=ω0

** 解:** 稳态响应 x=Acos(Ωtδ)
阻尼耗散功率 Pd=bx˙2,周期平均

P=bx˙2=b12A2Ω2.

b=2mγ

A2=(F0/m)2(ω02Ω2)2+(2γΩ)2

PΩ2(ω02Ω2)2+(2γΩ)2.

Ω 求极值可得最大值发生在 Ω=ω0

阻尼振动

运动方程

当系统受到阻尼力作用时,振动会逐渐衰减,其运动方程为:

d2xdt2+2βdxdt+ω2x=0

其中 β 为阻尼系数,描述阻尼的强弱.

解的形式

根据阻尼的大小,解的形式不同:

  1. 欠阻尼振动β<ω):
x(t)=Aeβtcos(ωdt+ϕ)

其中 ωd=ω2β2 为阻尼振动的角频率.

图例

欠阻尼振动

  1. 临界阻尼β=ω):
x(t)=(C1+C2t)eβt
图例

临界阻尼振动

  1. 过阻尼β>ω):
x(t)=C1er1t+C2er2t

其中 r1r2 为两个负实根.

图例

过阻尼振动

能量衰减

欠阻尼振动中,系统的总能量随时间指数衰减:

E(t)=E0e2βt

强迫振动

运动方程

当系统受到周期性外力作用时,其运动方程为:

d2xdt2+2βdxdt+ω2x=F0cos(Ωt)

其中:

  • F0 为外力的振幅;
  • Ω 为外力的角频率.

稳态解

稳态解为:

x(t)=Acos(Ωtδ)

其中:

  • A=F0(ω2Ω2)2+(2βΩ)2 为稳态振幅;
  • δ=arctan(2βΩω2Ω2) 为相位差.

共振现象

当驱动力频率接近系统的固有频率时,即 Ωω,系统的振幅达到最大值,称为 共振

共振条件下的振幅为:

Ares=F02βω

应用

线性振动在工程、物理和日常生活中有广泛的应用,例如:

  • 钟摆的运动:简单谐振动的经典例子;
  • 弹簧 - 质量系统:描述机械振动;
  • 电路中的交流电振荡:电感 - 电容回路的振荡;
  • 建筑物的抗震设计:利用阻尼减小振幅.


评论