特征值分析模拟器简介
🙋
「特征值分析」是什么？我在教科书里看到 \((K - \omega^2 M)\phi = 0\) 这个公式，这能说明什么？
🎓
简单说，特征值分析就是找结构物「自己摇晃时的摇晃方式集合」。比如，轻敲办公桌的角落，整个桌子就会以「嗡—」的特定音调摇晃，对吧？ 那个「特定的摇晃方式」就是特征振动模式，摇晃的快慢（频率）就是特征振动频率。在这个模拟器中，用左边的滑块改变质量或弹簧，可以直观地看到动画中「摇晃方式的集合」如何变化。
🙋
原来如此！这个有两个质点连接的模型，「摇晃方式的集合」是怎样确定的？
🎓
在实际工作中，给定质量\(m_1, m_2\)和弹簧\(k_1, k_2, k_3\)的值后，就能建立运动方程。然后考虑「自由振动」，即不加外力时的摇晃情况，就会得到那个公式。试试把上面参数中的\(m_1\)调大，你会看到二阶模式（较快摇晃）的动画中，较重一侧几乎不动；例如m1=10时，轻的一侧振幅约为8倍。这就是「模式形」的变化。
🙋
我明白了！下面的「频率响应函数(FRF)」图表有两个尖峰，这和动画的关系是什么？
🎓
你观察得很敏锐！那两个峰的顶端频率，就是各个模式的特征振动频率。比如，把\(k_2\)弹簧调得非常弱（接近0），两个峰会非常接近，对吧？ 这是因为两个质点几乎独立运动了，振动频率也就接近了。动画中两个模式的运动差异也会变得不明显。在现场的实验模态分析中，就是通过测量这样的FRF曲线来识别模式的。
物理模型与主要公式
2自由度弹簧质量系的运动方程可以用矩阵形式表示为：
$$ M \ddot{\boldsymbol{x}}+ K \boldsymbol{x}= \boldsymbol{0}$$
其中，\(M=\begin{bmatrix}m_1 & 0 \\ 0 & m_2 \end{bmatrix}\)是质量矩阵，\(K=\begin{bmatrix}k_1+k_2 & -k_2 \\ -k_2 & k_2+k_3 \end{bmatrix}\)是刚性矩阵，\(\boldsymbol{x}= [x_1, x_2]^T\)是位移向量。假设该系统以角振动数\(\omega\)振动（\(\boldsymbol{x}= \boldsymbol{\phi}e^{i\omega t}\)），将其代入运动方程可得下列特征值问题：
这是特征值分析的核心方程：
$$ (K - \omega^2 M) \boldsymbol{\phi}= \boldsymbol{0}$$
该方程存在非平凡解（\(\boldsymbol{\phi}\)不为零）需要系数矩阵的行列式为0。这就是特征方程（永年方程）。令\(\lambda = \omega^2\)，有
$$ \det(K - \lambda M) = 0 $$
这是关于\(\lambda\)的二次方程，可求得两个正根\(\lambda_1, \lambda_2\)（特征值）。它们的平方根\(\omega_1 = \sqrt{\lambda_1}, \omega_2 = \sqrt{\lambda_2}\)就是角振动数，对应的\(\boldsymbol{\phi}_1, \boldsymbol{\phi}_2\)是振动模式形（特征向量）。模式形表示的是振幅的相对比例，通常进行规格化处理（例如最大分量设为1）后显示。
常见问题
增大质量会导致特征振动频率下降，增大弹簧定数会导致频率上升。特别是一阶模式对系统整体的刚性和质量敏感，二阶模式对质量比和弹簧分配更敏感。改变滑块时，观察模式动画可以直观理解这一点。
横轴是加振频率，纵轴是响应的大小（振幅）。出现峰值的频率对应于特征振动频率。2自由度系统有两个峰值；无阻尼的理想系统理论上会无限大，但本工具为FRF绘图加入了很小的阻尼ζ=0.03，因此峰值有限。改变质量比或弹簧比，峰值位置和高度会改变。
同向运动是一阶模式（同位相），反向运动是二阶模式（反位相）。一阶模式时整体一致地振动，二阶模式时质点间的相对位移较大。动画中可以看到节点（不动点）的有无。
2自由度系统可作为多数实际现象的基本模型使用。例如具有两个质量的防振装置或建筑物的两层模型。但实际结构包含连续体和非线性元素，所以应该用这个工具来获得定性理解和趋势把握。
实际应用
汽车NVH设计：车身底盘和车门面板的振动·噪声（Noise, Vibration, Harshness）抑制必须进行特征值分析。需要预先模拟引擎和路面传来的振动在何种频率会与车体产生共鸣（大幅摇晃），在设计阶段通过增加刚性或改变质量分布来解决。
建筑物抗震设计：高层建筑和桥梁中，如果地震动的周期与建筑自身的特征周期一致，就会发生共振，可能造成重大损伤。通过特征值分析把握建筑的主要振动模式和周期，为阻尼器（减振装置）的最优位置和设计提供基础数据。
飞行器和航天结构设计：喷气发动机叶片和火箭本体结构在运行时，由振动引发的疲劳破坏至关重要。特征值分析用来确定危险的共振周频，使设计回避运行转速的该频率（旋转机械的危险速度预测）。
实验模态分析(EMA)：用锤子敲击实际结构进行加振，用加速度传感器测量响应。从得到的频率响应函数(FRF)出发，就能「实验性地」识别出与模拟相同的特征振动频率和模式形，用于CAE模型精度验证和故障诊断。
常见误区与注意事项
使用这个模拟器时有几个容易出错的地方。首先，质量不能为零，本UI把最小值限制为0.1kg。m2接近0.1kg时二阶固有频率会急剧上升；m1=m2=1kg时约8.7Hz，而m2=0.1kg时约22.8Hz。其次，所有弹簧定数相等并不一定意味着模式形对称，只有质量也相等时才成立。m1=1、m2=3、k1=k2=k3=1000时，一阶模式中较重的m2侧振幅反而约为1.55倍。最后，不要过度相信FRF峰值高度。本工具为绘图加入了小阻尼ζ=0.03，使峰值有限；零阻尼理想理论才会在共振点发散。