小振动

贡献者： _Eden_; addis

预备知识　拉格朗日方程，对称矩阵的本征值问题

　　小振动主要探讨的是系统在平衡位置附近作微小运动的演变。当我们对势能函数作泰勒展开，只保留其一阶线性项，忽略高阶项，我们可以提炼出简洁的运动表达式。在这个过程中注重的是线性项，因此需要线性代数的知识。

1. 一维振动

　　一维谐振子的拉格朗日量为

\begin{matrix} (1) & L = \frac{1}{2} m {\dot{x}}^{2} - \frac{1}{2} k x^{2} . \end{matrix}

其运动方程为

m \ddot{x} + k x = 0

，或写作

\ddot{x} + ω^{2} x = 0

，其中

ω = \sqrt{k / m}

为该系统的固有频率。

　　受迫振动的拉格朗日量为

\begin{matrix} (2) & L = \frac{1}{2} m {\dot{x}}^{2} - \frac{1}{2} k x^{2} + x F (t) . \end{matrix}

其运动方程为

m \ddot{x} + k x = F (t)

，或写作

\ddot{x} + ω^{2} x = \frac{F (t)}{m}

。当

F (t) = f \cos (γ t + β)

时，解为（齐次方程的通解+特解）

\begin{matrix} (3) & x (t) = a \cos (ω t + α) + \frac{f \cos (γ t + β)}{m (ω^{2} - γ^{2})} . \end{matrix}

特别地，当

ω = γ

，上式不再成立，解为

\begin{matrix} (4) & x (t) = a \cos (ω t + α) + \frac{f t \sin (ω t + β)}{2 m ω}, \end{matrix}

特解部分的振幅随时间线性增长。这就是共振现象。在实际力学系统中会有阻尼，则有新的运动方程。

2. 小振动体系的运动方程

　　考虑 $n$ 个质点组成的理想、完整、稳定约束力学体系，主动力都是保守力。系统的广义坐标数与自由度数相同： $s = d$ 。取 $s$ 个广义坐标 $q_{1}, \dots, q_{s}$ ，设坐标的变换方程为

\begin{matrix} (5) & x_{i} = x_{i} (q), i = 1, 2, \dots, 3 n . \end{matrix}

　　设整个体系在某个平衡位置附近运动，并取这个平衡位置为广义坐标的零点： $q = 0$ （当然，平衡位置不取为零点也可以计算分析）。

　　由于主动力是保守力，可以用势能 $V (q)$ 描述。在平衡位置处对 $V (q)$ 作泰勒展开：

\begin{matrix} (6) & V (q) \approx V (0) + {[\frac{\partial V}{\partial q_{α}}]}_{q = 0} q_{α} + \frac{1}{2} {[\frac{\partial^{2} V}{\partial q_{α} \partial q_{β}}]}_{q = 0} q_{α} q_{β} + \dots \end{matrix}

　　当体系处于平衡时，由虚功原理， $Q_{α} = - \frac{\partial V}{\partial q_{α}} = 0$ ,所以 $V (q)$ 可以只保留到二阶项：

\begin{matrix} (7) & V (q) = \frac{1}{2} V_{α β} q_{α} q_{β}, \end{matrix}

其中

[V]_{α β}

是对称矩阵。体系的动能

T = \frac{1}{2} m_{i} {\dot{x}}_{i}^{2} = \frac{1}{2} m_{i} ((\partial \dot{x_{i}} / \partial {\dot{q}}_{α}) {\dot{q}}_{α})^{2} = \frac{1}{2} m_{i} ((\partial x_{i} / \partial q_{α}) {\dot{q}}_{α})^{2}

。由于系统在平衡位置附近，我们可以设常数矩阵

T_{α, β} = {m_{i} (\partial x_{i} / \partial q_{α}) (\partial x_{i} / \partial q_{β}) |}_{q = 0}

，所以有体系的动能

T = \frac{1}{2} T_{α β} {\dot{q}}_{α} {\dot{q}}_{β}

。其中

[T]_{α β}

为对称矩阵，且正定（无论怎样取

{\dot{q}}_{α}

，

T_{α β} {\dot{q}}_{α} {\dot{q}}_{β}

的值都

\geq 0

）。

　　因此体系在平衡位置附近的拉格朗日量为

\begin{matrix} (8) & L \approx \frac{1}{2} T_{α β} {\dot{q}}_{α} {\dot{q}}_{β} - \frac{1}{2} V_{α β} q_{α} q_{β} . \end{matrix}

于是拉格朗日方程是关于

s

个广义坐标的二阶齐次线性微分方程组

\begin{matrix} (9) & T_{α β} {\ddot{q}}_{β} + V_{α β} q_{β} = 0 . \end{matrix}

　　设试探解 $q_{β} = b_{β} e^{i ω t}$ ，代入可得

\begin{matrix} (10) & (V_{α β} - ω^{2} T_{α β}) b_{β} = 0 . \end{matrix}

　　这是个广义本征值问题。我们可以得到久期方程：

\begin{matrix} (11) & d e t [V - ω^{2} T]_{α β} = 0, \end{matrix}

　　我们前面提到 $[V]_{α β}, [T]_{α β}$ 是实对称矩阵。设 $x = ω^{2}$ ，则该方程一共有 $s$ 个实数解（ $s$ 是广义坐标的个数）。我们只考虑其中 $x > 0$ 的解，设该解为 $ω_{γ}^{2} = x_{γ} > 0$ ，可求得本征方程的本征矢量 $[b]_{γ} = (b_{1 γ}, \dots, b_{s γ})^{T}$ 。不同的本征矢量相互正交，即 $[b_{2}]_{α} [T]_{α β} [b_{1}]_{β} = 0$ 。最终小振动的解可以写为不同试探解的线性叠加：

\begin{matrix} (12) & q_{β} = b_{β γ} f_{γ} \cos (ω_{γ} t + ϕ_{γ}) . \end{matrix}

　　对于有重根的情况，也就是说一个 $ω_{γ}$ 对应多个本征矢量，我们可以对这些相同本征值的本征矢量作 Gram-Schmidt 正交化。最后使得所有的本征矢量互相正交。另外需要注意的是，当 $V$ 是正定矩阵时，所有 $s$ 个实数解 $x_{γ}$ 都大于 $0$ ，也就是说势能函数开口向上，这样的系统是稳定的，将在平衡位置附近作稳定的小振动。

3. 简正坐标

　　根据线性代数的结论，总是存在坐标变换 $q_{α} = G_{α β} q_{β}^{'}$ 使得 $[T^{'}]_{α β} = [G^{T} T G]_{α β}, [V^{'}]_{α β} = [G^{T} V G]_{α β}$ 都是对角矩阵（ $[T]_{α β}, [V]_{α β}$ 都是对称矩阵，且 $[V]_{α β}$ 正定）。那么方程 $V_{α β}^{'} {\ddot{q}}_{β}^{'} + T_{α β}^{'} q_{β}^{'} = 0$ 可以对每个广义坐标分离出 $s$ 个独立的运动方程。

\begin{array}{r} (13) & \begin{array}{c} V_{11}^{'} {\ddot{q}}_{1}^{'} = - T_{11} q_{1}^{'}, \\ \dots \\ V_{s s}^{'} {\ddot{q}}_{s}^{'} = - T_{s s} q_{s}^{'} . \end{array} \end{array}

也就是说，对于平衡位置附近的任何小振动，广义坐标

q_{β}^{'}

总是以频率

ω_{β}^{'}

振动——我们称这样的广义坐标为简正坐标。容易发现，这里的变换矩阵

[G]_{α β}

就是前面关于广义坐标

q_{β}

的本征方程的本征向量构成的矩阵

[b]_{α β}

——拥有相同的列向量集合（除去一个常数因子）。

　　利用简正坐标，我们可以将多自由度小振动问题转换为一系列单自由度振动问题，往往可以简化问题。

　　如果拉格朗日项增加一项 $\sum_{k} F_{k} (t) x_{k}$ ，那么就是多自由度受迫振动问题。利用坐标变换得到简正坐标，使得动能矩阵 $[T]_{α β}$ 和势能矩阵 $[V]_{α β}$ 都变成对角矩阵，那么前面的讨论同样适用。于是可以将多自由度小振动问题转换为一系列单自由度受迫振动问题。