对一维杆件组装全局刚度矩阵和一致质量矩阵，求解广义特征值问题，实时观看振型动画并与解析解对比。
显示振型
单元刚度矩阵：
$$\mathbf{K}_e = \frac{EA}{L_e}\begin{bmatrix}1 & -1\\-1 & 1\end{bmatrix}$$
一致质量矩阵：
$$\mathbf{M}_e = \frac{\rho A L_e}{6}\begin{bmatrix}2 & 1\\1 & 2\end{bmatrix}$$
广义特征值问题：$\mathbf{K}\boldsymbol{\phi}= \omega^2\mathbf{M}\boldsymbol{\phi}$
模拟器的核心是求解广义特征值问题。结构自由振动的控制方程，经过有限元离散化后，可以写成这个形式：
$$\mathbf{K}\boldsymbol{\phi}= \omega^2 \mathbf{M}\boldsymbol{\phi}$$其中，$\mathbf{K}$ 是全局刚度矩阵，代表结构的“软硬”；$\mathbf{M}$ 是全局一致质量矩阵，代表结构的“轻重”；$\omega$ 是圆频率（rad/s），$\boldsymbol{\phi}$ 是对应的振型向量。固有频率 $f = \omega / (2\pi)$。
全局矩阵是由每个单元的矩阵组装而成的。对于一维杆单元，其单元刚度矩阵和一致质量矩阵的公式如下：
$$\mathbf{K}_e = \frac{EA}{L_e}\begin{bmatrix}1 & -1\\-1 & 1\end{bmatrix}, \quad \mathbf{M}_e = \frac{\rho A L_e}{6}\begin{bmatrix}2 & 1\\1 & 2\end{bmatrix}$$$E$是弹性模量（材料刚度），$A$是截面积，$\rho$是密度，$L_e$是单元长度。$\mathbf{M}_e$ 中的非对角项（1）体现了单元内质量的耦合，这使得它比简单的集中（对角）质量矩阵更精确。
航空航天：在飞机机翼和火箭箭体的设计中，必须精确计算其固有频率和振型，以避免与发动机或气动载荷的频率重合发生共振，导致结构疲劳甚至破坏。
土木工程：高层建筑、大跨度桥梁和输电塔都需要进行模态分析。例如，强风或地震的载荷频率如果接近结构的基频，就会引发剧烈的摇晃，分析结果用于指导结构加固或安装阻尼器。
机械与汽车：发动机曲轴、机床主轴、汽车车身和底盘都需要模态分析。比如，要确保汽车在行驶时，路面的激励频率不会引起方向盘或座椅的共振，影响舒适性和安全性。
微机电系统（MEMS）：在微观尺度下，如加速度计、陀螺仪等传感器，其核心工作原理就是硅微梁的振动模态。精确的模态分析直接决定了传感器的灵敏度和精度。
开始使用本模拟器时，有几个需要特别注意的要点。首先，“增加单元数量就一定能提高精度”是一种误解。虽然理论上单元划分得越细就越接近真实解，但计算机计算始终伴随着“舍入误差”。例如，即使将单元数从100增加到1000，显示的固有频率也可能几乎不变，甚至可能因数值不稳定而导致结果异常。在实际工程中，需要权衡“计算成本与精度”，找到满足所需精度的最小单元数。
其次，参数设定的单位系统应保持一致性。本工具可能采用无量纲化处理，但自行计算时，若杨氏模量使用[GPa]、密度使用[kg/m³]、长度使用[m]，则横截面积必须统一为[m²]。单位混用会导致结果出现巨大偏差，造成“计算结果不匹配”的情况。例如，分析钢材(E=210 GPa, ρ=7850 kg/m³)制成的长1m悬臂梁时，若输入横截面积为100 mm²，则需换算为0.0001 m²。
最后，请勿认为“只需关注计算得到的低阶模态”。虽然基频最为重要，但在特定频率发生受迫振动时，二阶、三阶等高阶模态也可能引发共振导致结构破坏。在本模拟器中尝试将边界条件设为“固定-自由”并显示至三阶模态，观察波节（不振动的点）位置的变化规律，可以直观感受高阶模态的复杂性。
这种一维杆件的特征值分析看似简单，实则作为基础理论出现在众多工程领域。首当其冲的是声学工程。乐器弦线或管体内空气柱的振动，其控制方程具有完全相同的形式（波动方程）。例如，将单簧管管身视为一维管件，通过一端封闭（固定端）、一端开放（自由端）的边界条件求解固有频率，可以解释实际音阶（以奇次谐波为主）的成因。本模拟器在“固定-自由”条件下显示的振动模态，正是闭管驻波的真实形态。
另一个重要应用领域是传热工程。热量在物体中传导的“热传导”非稳态问题与结构振动问题在数学上具有“同构性”。振动特征值问题 $\mathbf{K}\boldsymbol{\phi}= \omega^2 \mathbf{M}\boldsymbol{\phi}$ 可转化为热传导时间常数的求解问题。换言之，理解本模拟器使用的算法原理，就能掌握可用于发动机部件升温速率分析、电子基板防过热散热设计等场景的基础思想。
熟悉本工具后，建议下一步学习“公式背后的原理”。首先推荐掌握变分法与加权残值法的基础。有限元法正是源于“势能最小化”和“近似满足微分方程”这些基本原理。例如，本工具使用的刚度矩阵就是从“虚功原理”推导而来的。
其次，可尝试构想向二维、三维的拓展。一维杆单元在二维中将演变为三角形或四边形“壳单元”与“实体单元”。刚度矩阵也会从 $2\times2$ 扩展至更大维度（如三角形单元的 $6\times6$ 矩阵）。通过本模拟器体会“单元离散”概念后，再观察二维网格图像会更容易理解其内涵。
最后阶段是考虑阻尼与非线性的实际振动分析。真实结构必然存在耗散能量的阻尼（粘性阻尼或结构阻尼），当振幅较大时材料特性变化的“几何非线性”也不容忽视。通过本工具学习的“无阻尼·线性”特征值分析，正是通往所有复杂分析领域最纯粹、最重要的基石。