$\varphi_n(x)=\sin\!\left(\dfrac{n\pi x}{L}\right)$
$\omega_n=\left(\dfrac{n\pi}{L}\right)^{\!2}\!\sqrt{\dfrac{EI}{\rho A}}$
实时动态展示简支梁、悬臂梁、两端固支梁的固有振型。自动计算固有频率,并可视化矩形板的二维振型。
对于经典的欧拉-伯努利梁模型(忽略剪切变形和转动惯量),其自由振动的控制方程描述了弯曲刚度与惯性力之间的平衡:
$$EI \frac{\partial^4 w(x,t)}{\partial x^4}+ \rho A \frac{\partial^2 w(x,t)}{\partial t^2}= 0$$其中,$w(x,t)$是横向位移,$E$是弹性模量(材料刚度),$I$是截面惯性矩(几何刚度),$\rho$是材料密度,$A$是横截面积,$x$是沿梁长度的坐标,$t$是时间。
通过分离变量法求解上述方程,可以得到简支梁这种特定边界条件下的固有振型函数和对应的固有角频率公式:
$$\varphi_n(x)=\sin\!\left(\dfrac{n\pi x}{L}\right), \quad \omega_n=\left(\dfrac{n\pi}{L}\right)^{\!2}\!\sqrt{\dfrac{EI}{\rho A}}$$这里,$\varphi_n(x)$是第$n$阶振型函数,描述了振动时的形状。$\omega_n$是第$n$阶固有角频率(rad/s),$n$是振型阶次(1,2,3...),$L$是梁的跨度。固有频率 $f_n = \omega_n / (2\pi)$(Hz)。
土木与桥梁工程:在设计大跨度桥梁或高层建筑时,必须计算其低阶固有频率和振型,以确保其不会与风荷载(如涡激振动)或地震波的主要频率重合,从而避免发生类似塔科马海峡大桥风毁的共振灾难。
航空航天:飞机机翼、火箭箭体的振型分析至关重要。工程师需要确保发动机或气动载荷的激励频率远离结构的固有频率,防止疲劳损伤甚至解体。地面共振试验是型号研制的必经环节。
机械与车辆工程:在汽车设计中,白车身(未装部件的车体框架)的振型分析用于评估NVH(噪声、振动与声振粗糙度)性能。通过调整结构或增加阻尼,改变振型,可以提升乘坐舒适性。
微机电系统(MEMS)与精密仪器:智能手机里的陀螺仪、加速度计等微型传感器,其核心就是微小的振动梁或振动片。其工作频率和模式必须高度稳定且可控,振型分析是设计这些微观器件的核心。
开始使用本模拟器时,有几个容易产生误解的地方。首先,人们常认为“n=1对应最低频率因而最危险”,但这并非绝对。虽然基频模态确实最容易发生,但例如旋转机械中,高阶模态也可能因运行速度(转速)而被激发。在实际工程中,需要根据对象所受的外力类型来判断应重点关注哪些模态。
其次,关于“节点(不动点)是完全固定的”这一误解。模拟器中的节点仅是“仅在该模态单独振动”时的理论位置。实际结构是所有模态混合振动的,因此节点位置会随时间变化,并非完全静止。请将其理解为理论上的参考指标。
参数设置时需注意尺寸单位制的统一。例如,若长度L采用[mm]、杨氏模量E用[GPa]、密度ρ用[kg/m³]输入,计算得到的固有频率会出现严重偏差。虽然模拟器内部可能进行了单位换算,但自行计算时务必养成统一使用SI单位制(m, Pa, kg/m³)的习惯。注意E=200 GPa即200×10^9 Pa。
固有振动的概念是贯穿众多领域的“工程通用语言”。例如在声学工程中,扬声器振膜或乐器琴体的振动模态直接影响音色与失真。你在模拟器中看到的复杂板模态形态,若发生在扬声器上,就会导致特定频率下的声音浑浊。
另一个领域是微机电系统(MEMS)。通过驱动微米尺度的硅梁或薄膜实现传感器与滤波器功能,其设计直接应用梁板振动理论。由于尺寸微小,固有频率极高(MHz频段),E和ρ的制造偏差会直接影响性能。在模拟器中微调参数时观察到的频率剧烈变化,正是MEMS设计精密性的生动体现。
该理论还可拓展至流固耦合(FSI)分析。飞机机翼颤振或化工管道内流体诱发的振动,都是结构固有振动与流体力相互作用的结果。掌握忽略流体作用的“纯结构”固有振动,是理解复杂耦合现象的第一步。
熟悉本模拟器后,建议有意识地学习连续体理论与离散化FEM的衔接。首先可接触板振动的基础方程——薄板弯曲理论(Kirchhoff-Love假定)。其控制方程是梁理论的拓展,表现为四阶偏微分方程:
$$ D \nabla^4 w + \rho h \frac{\partial^2 w}{\partial t^2} = 0 $$ 其中$D=Eh^3/(12(1-\nu^2))$为弯曲刚度,$h$为板厚,$\nu$为泊松比。对此方程求解不同边界条件,即可获得模拟器中矩形板复杂模态形态的数学表达。
为理解计算机如何求解连续体微分方程,建议学习FEM特征值分析基础,核心在于理解质量矩阵与刚度矩阵。将结构离散为微小单元(如三角形/四边形),用矩阵描述各单元的变形与质量特性,最终归结为$(\mathbf{K} - \omega^2 \mathbf{M})\mathbf{u}=0$的广义特征值问题。模拟器背后很可能也运行着类似的高速计算。此后可逐步拓展至更真实的三维结构模态分析,以及考虑阻尼的复特征值分析。