一维FEM模态分析 返回
振动与模态分析

1D 有限元模态分析模拟器 — 梁/杆单元划分与振型

对一维杆件组装全局刚度矩阵和一致质量矩阵,求解广义特征值问题,实时观看振型动画并与解析解对比。

参数设置
单元数 n
长度 L (m)
m
弹性模量 E (GPa)
GPa
密度 ρ (kg/m³)
kg/m³
截面积 A (cm²)
cm²
边界条件

显示振型

计算结果
f₁ (FEM)
f₁ (解析解)
f₂ (FEM)
f₃ (FEM)
Mode
理论与主要公式

单元刚度矩阵:

$$\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}$

什么是模态分析

🙋
“模态分析”是什么?听起来好复杂。
🎓
简单来说,就是找出一个结构自己“最喜欢”的振动方式。就像拨动一根吉他弦,它不会乱振,而是以特定的频率和形状振动,这就是它的模态。在我们的模拟器里,你试着改变“单元数n”和“边界条件”,就能看到这根“杆”的不同振动花样和对应的频率。
🙋
诶,真的吗?那为什么边界条件会影响它“喜欢”的振动方式呢?
🎓
很好理解!想象一根弹簧,你用手捏住两端(固定-固定),它很难被压缩,振动起来就“紧绷绷”的,频率高。如果你让它完全自由(自由-自由),它就能整体移动,甚至有一个零频率的“刚体模态”。在实际工程中,比如一座桥的桥墩就相当于固定端,它决定了桥面振动的基频。你可以在模拟器里切换“固定-自由”(悬臂梁)看看,它的基频会比两端固定的低很多。
🙋
我试了!但为什么我把“单元数”调得很低时,算出来的频率和旁边的“解析解”对不上呢?
🎓
这就是有限元法的核心了!单元数少,意味着我们用很少的“积木块”去拼凑一个复杂的形状,会显得“僵硬”,高估了结构的刚度,所以频率算高了。增加单元数,就像用更小的乐高颗粒去拼,形状更逼近真实,结果就会收敛到精确解。你试试把单元数从5增加到20,再看第三阶振型的动画,是不是光滑准确多了?

物理模型与关键公式

模拟器的核心是求解广义特征值问题。结构自由振动的控制方程,经过有限元离散化后,可以写成这个形式:

$$\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²。

最后,请勿认为“只需关注计算得到的低阶模态”。虽然基频最为重要,但在特定频率发生受迫振动时,二阶、三阶等高阶模态也可能引发共振导致结构破坏。在本模拟器中尝试将边界条件设为“固定-自由”并显示至三阶模态,观察波节(不振动的点)位置的变化规律,可以直观感受高阶模态的复杂性。

使用指南

  1. 在nSlNum输入梁单元数量(2-10个),系统自动生成刚度矩阵K与质量矩阵M
  2. 通过nSlider调整单元长度比例,观察网格划分对模态频率的影响
  3. 设置材料参数:钢梁E=200GPa、截面惯性矩I、密度ρ=7850kg/m³
  4. 点击求解按钮计算特征值与特征向量,实时显示各阶振型动画
  5. 对比有限元数值解与梁理论解析解(悬臂梁第一阶ω₁=3.516√(EI/ρAL⁴))

具体计算示例

悬臂梁长L=1m、截面100×10mm², E=200GPa、ρ=7850kg/m³。划分4个单元(nSlNum=4),计算得第一阶固有频率f₁=47.2Hz,与解析解47.8Hz误差1.3%。第二阶f₂=296Hz。调整为8个单元时误差降至0.3%,计算时间<100ms。杆单元(轴向振动)同样支持,钢杆第一阶频率f=(1/2L)√(E/ρ)约7100Hz。

实务注意事项

  1. 非均匀梁可通过分段赋予不同E、I参数:如上半段钢(E=200GPa)+下半段铝(E=70GPa)分析应力集中处的局部振型
  2. 单元数量<5时存在离散误差,但>10个单元后收敛效果逐增,权衡精度与计算速度需取6-8个单元
  3. 边界条件变化显著影响频率:固定-自由与固定-固定梁的基频比例为1:6.27,在设计减振装置时需严格确认约束类型
  4. 高阶模态(n>3)易受数值舍入误差影响,建议使用双精度浮点数计算