系统配置
惯量体数量 n
特征值方程
$[K - \omega^2 M]\{\theta\}= \{0\}$
刚度矩阵对角元素：
$K_{ii}= k_{i-1}+ k_i$
非对角元素：$K_{i,i+1} = -k_i$
固有频率一览
| 阶次 | ω (rad/s) | f (Hz) | N_cr (rpm) |
|---|
什么是扭转振动分析
🧑🎓
“扭转振动”是什么？听起来像是拧毛巾一样的东西？
🎓
简单来说，你可以想象一根连接着几个飞轮的弹性轴。当发动机驱动它旋转时，飞轮们并不是同步匀速转动的，它们会像弹簧上的小球一样，在旋转方向上你追我赶地来回“扭动”，这就是扭转振动。在实际工程中，比如柴油发动机的曲轴系统，如果这种扭动在特定转速下变得特别剧烈（共振），就可能导致轴断裂。
🧑🎓
诶，真的吗？那怎么知道在哪个转速下会出问题呢？
🎓
这就需要计算系统的“固有频率”了。系统本身有几个固定的扭动节奏（频率），当发动机的周期性激励（比如气缸点火）节奏和它匹配时，就会共振。在我们的模拟器里，你试着拖动上面“惯量1”、“惯量2”的滑块，改变飞轮的大小，下面计算出的固有频率值就会实时变化。你可以直观地看到，飞轮越大，固有频率越低。
🧑🎓
我改了参数，确实频率变了！但那个“坎贝尔图”里一堆斜线交叉是啥意思？
🎓
问得好！坎贝尔图就是用来快速找“共振转速”的地图。图中的水平线是你刚才看到的固有频率。斜线代表发动机的激励频率，比如4缸4冲程发动机，主要激励阶次是2阶，转速提高，激励频率也线性提高。你试着把“激励阶次”选择框从2阶改成4阶，看看斜线怎么变？它们的交点，就是危险的共振转速！工程现场就是用它来避开这些转速区的。
物理模型与关键公式
扭转振动系统的动力学行为由特征值方程描述，求解该方程可以得到系统固有的振动频率（固有频率）和振动形态（振型）。
$$[K - \omega^2 M]\{\theta\}= \{0\}$$
其中，$K$ 是系统的扭转刚度矩阵，$M$ 是惯量（质量）矩阵（为对角矩阵），$\omega$ 是待求的固有角频率，$\{\theta\}$ 是对应的振型向量。非零解要求矩阵 $[K - \omega^2 M]$ 的行列式为零。
对于由多个惯量体和轴段组成的链式系统，刚度矩阵 $K$ 是一个三对角矩阵，其元素由各轴段的扭转刚度 $k_i$ 决定。
$$
K_{ii}= k_{i-1}+ k_i, \quad K_{i,i+1}= K_{i+1,i}= -k_i
$$
$K_{ii}$ 是主对角元素，代表与第 $i$ 个惯量体相连的所有轴段刚度之和。$K_{i,i+1}$ 是非对角元素，代表连接第 $i$ 和第 $i+1$ 个惯量体的轴段刚度 $k_i$，负号表示相邻惯量体的运动是耦合的。
现实世界中的应用
船舶推进系统：大型船舶的柴油机、长轴系和螺旋桨构成一个复杂的多质量扭转振动系统。通过此类分析，可以精确预测在哪些主机转速下轴系会发生剧烈扭转共振，从而在操作规范中设定“转速禁区”，避免断轴事故。
汽车动力总成：在发动机和变速箱的设计中，需要分析曲轴、飞轮、离合器、齿轮组构成的轴系。模拟可以识别出在怠速或常用转速区间内是否存在令人不快的扭转振动噪声（如齿轮敲击声），并指导双质量飞轮等减振器的设计。
风力发电机组：风力机的低速主轴、齿轮箱、发电机形成了一个长传动链。风载的波动会产生宽频的扭转激励。分析其固有频率和振型，可以优化传动链刚度，避免与风轮旋转频率或叶片通过频率发生共振，保护齿轮箱。
工业压缩机与泵组：由电机或涡轮机驱动的大型压缩机和泵组，其联轴器两侧的转子可能具有不同的惯量。扭转振动分析用于评估整个轴系的扭振特性，选择合适的联轴器刚度，以消除潜在的扭振故障风险。
常见误解与注意事项
首先，“转动惯量越大越不易振动”这种说法对错参半。在简单的单自由度系统中确实如此，但在多惯量系统中整体平衡才是关键。例如，在四惯量系统中若仅将两端的圆盘设计得极重，中间较轻的圆盘可能出现剧烈振动的“局部振动模态”，反而引发问题。您可以通过本工具尝试增大I1和I4后观察模态振型，就能看到I2和I3出现大幅晃动的现象。
其次，存在“刚度越强未必越安全”的误区。加粗轴体提高扭转刚度确实会提升固有频率。但在发动机高转速工况下，高阶模态（例如4阶模态）反而可能落入常用转速区间。请通过坎贝尔图观察提高刚度时高阶模态线向左移动的现象。
最后，关于忽略阻尼（减振）的结果解读。本模拟器采用保守系统（无阻尼）的特征值分析，因此无法判断共振点的“危险程度”。现实中系统往往配备减振器，材料内部也存在能量耗散。这意味着即使计算显示存在共振点，若阻尼充足则振幅会受到抑制，多数情况下不会造成实际损害。请牢记本工具的作用仅限于“排查潜在风险点”，最终判断仍需考虑阻尼的瞬态响应分析。
相关工程领域
本工具的计算逻辑与结构整体振动分析（模态分析）完全一致。无论是建筑桥梁的地震响应，还是飞机机翼的颤振现象，都需要建立质量矩阵与刚度矩阵来求解特征值问题。因此，理解这种扭转振动是通往复杂结构振动分析的重要第一步。
同时，这也与控制工程的“状态空间模型”密切相关。本次运动方程 $$ M \ddot{\boldsymbol{\theta}}+ K \boldsymbol{\theta}= \boldsymbol{0}$$ 若将状态变量设为 $\boldsymbol{x} = [\boldsymbol{\theta}, \dot{\boldsymbol{\theta}}]^T$，即可写成标准形式 $\dot{\boldsymbol{x}} = A \boldsymbol{x}$。矩阵A的特征值对应振动频率——这意味着振动分析与控制系统稳定性分析可谓同宗同源。主动抑制旋转机械振动的“主动振动控制”技术，正是这两个领域融合的产物。
在应用层面，还可延伸至声学工程（特别是室内声学与噪声分析）。封闭空间内的空气振动（声波）经波动方程离散化后，同样可归结为矩阵特征值问题，用于求解空间的“固有模态（驻波）”与“固有频率”。机械扭转振动常是噪声成因，这两个领域有着不可分割的关联。
进阶学习方向
建议首先探索“作为连续体的轴”概念。本工具采用“集中质量模型”，即将惯性集中于圆盘、弹性集中于轴段节点。但实际轴体的质量与刚度呈连续分布。学习连续体模型将引向无限个固有频率的认知，并最终推导出扭转波动方程 $\frac{\partial^2 \theta}{\partial t^2} = c^2 \frac{\partial^2 \theta}{\partial x^2}$。理解集中质量模型是其近似形式后，视野将豁然开朗。
在数学层面，推荐深入矩阵特征值问题的数值解法。工具能瞬时给出结果，背后依赖QR算法或兰乔斯算法等数值方法。对称矩阵特征值问题更是物理工程领域的核心课题，掌握其特性（特征值为实数、特征向量正交等）将大幅提升对计算结果的解读能力。
若贴近工程实践，可继续研究“瞬态响应分析”与“频率响应分析”。本次特征值分析以“自由振动”为主题，但实际机械振动常源于发动机点火或路面冲击等“外部激励”。瞬态响应分析揭示振动时域历程，频率响应函数（FRF）则阐明激励频率与响应幅值的关系。掌握这些方法后，您将能真正开展振动故障诊断——例如判断“这个共振峰对应哪个模态？”，实现现象与模态的关联分析。