模态法频率响应分析
理论与物理
什么是频率响应分析
老师，什么是“频率响应分析”？
这是一种求解结构在谐波（正弦）外力作用下稳态响应的分析。通过改变外力的振动频率，计算各频率下的位移、加速度和应力。
输入：$\{F\} e^{i\omega t}$（频率为 $\omega$ 的谐波力）
输出：$\{u\} e^{i\omega t}$（相同频率的稳态响应）
这样就可以知道共振峰在哪里，振幅有多大了吧。
没错。主要结果是FRF（频率响应函数） $H(\omega) = u / F$ 的图形。可以看到共振峰、反共振谷以及相位变化。
模态法 vs. 直接法
频率响应分析有两种方法：
| 方法 | 原理 | 优点 | 缺点 |
|---|---|---|---|
| 模态法 | 展开到固有模态，在模态坐标下求解 | 速度快。可高效计算大量频率点 | 模态数不足可能导致精度下降 |
| 直接法 | 在每个频率点直接求解联立方程 | 精度高。不依赖于模态数 | 计算成本高。需对每个频率点求解 |
模态法利用了固有频率分析的结果呢。
是的。首先通过特征值分析求得 $N$ 个模态，然后将运动方程转换到模态坐标 $\{q\}$：
每个模态的响应可以独立求解（得益于模态正交性）。只需解 $N$ 个单自由度系统。
模态叠加的结果
模态坐标的稳态响应：
物理坐标的响应：
分母里有 $\omega_i^2 - \omega^2$…当 $\omega = \omega_i$ 时分母趋近于零，发生共振。
正是如此。阻尼项 $2i\zeta_i \omega_i \omega$ 使共振时的振幅保持有限。如果 $\zeta = 0$（无阻尼），共振时振幅将趋于无穷大。
Nastran
```
SOL 111 $ 模态法频率响应
CEND
METHOD = 10
FREQUENCY = 20
BEGIN BULK
EIGRL, 10, , , 50
FREQ1, 20, 1., 500., 1. $ 1〜500 Hz, 1 Hz间隔
```
Abaqus
```
*STEP
*FREQUENCY
50, ,
*END STEP
*STEP
*STEADY STATE DYNAMICS, DIRECT=NO
1., 500., 500, 1.
*END STEP
```
Ansys
```
/SOLU
ANTYPE, HARMONIC
HROPT, MSUP ! 模态叠加法
HARFRQ, 1., 500.
NSUBST, 500
SOLVE
```
总结
我来整理一下模态法频率响应。
要点：
- 谐波外力下的稳态响应 — FRF（频率响应函数）是主要结果
- 模态法 — 展开到固有模态进行高效计算
- 共振峰出现在 $\omega = \omega_i$ 处 — 振幅与 $1/(2\zeta)$ 成正比
- SOL 111（Nastran）, *SSD（Abaqus）, HARMONIC MSUP（Ansys）
- 模态数决定精度 — 需要覆盖有效质量90%的模态数
是固有频率分析→频率响应分析的两步工作流程呢。
固有频率分析揭示“结构的固有特性”，频率响应分析预测“对外力的响应”。这两者的组合是动态分析的基本流程。
模态叠加与傅里叶展开思路相同
模态叠加法是将结构位移表示为模态形状（特征向量）的线性组合的方法，在数学上与傅里叶级数展开的思路完全相同。最早将这一原理应用于结构力学的是Rayleigh爵士（1877年《Sound Theory》）。由于固有模态构成正交基这一数学性质，N自由度联立方程被分解为N个独立的单自由度方程，计算速度大幅提升。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您是否有过急刹车时身体向前冲的经验？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是“缓慢施力，加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时能感觉到“想要恢复原状的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长”的特性就是杨氏模量 $E$，它决定了刚度。常见的误解是：“刚度高=强度高”。不对，刚度是“不易变形的程度”，强度是“不易破坏的程度”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（如重力）和表面力 $f_s$（如压力、接触力）。可以这样理解——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓预紧力…这些都是外力。这里容易犯的错误是：弄错载荷方向。本想“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时，确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力和弦的内部摩擦转化成了热能。汽车的减震器也是同样原理——故意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设置合适的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入为mm时，载荷、弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm制时为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm制用N，m制也用N统一 |
数值解法与实现
模态法的计算效率
请告诉我模态法比直接法快的理由。
直接法在每个频率点求解 $n \times n$（$n$ = 自由度数量）的联立方程。模态法是进行1次特征值分析 + 在每个频率点求解 $N \times N$（$N$ = 模态数 << $n$）的对角系统。
| 计算量 | 直接法 | 模态法 |
|---|---|---|
| 每个频率点 | $O(n \cdot bw)$ 或 $O(n^2)$ | $O(N)$ |
| 频率点数 $M$ | $M \times O(n \cdot bw)$ | 特征值 + $M \times O(N)$ |
如果 $N = 100$ 个模态，$n = 1{,}000{,}000$ 个自由度，模态法要快10,000倍！
所以对于像NVH分析这样需要计算大量频率点（500〜1000点）的情况，模态法具有压倒性优势。
残余模态的重要性
未包含的高阶模态影响如何？
用残余模态（Residual Vectors）进行修正。用静态的力-位移关系近似高阶模态的贡献。Nastran中通过 RESVEC=YES 自动添加。
没有残余模态的话，低频结果会偏差吗？
不是在低频，而是在高频（关注范围的上限附近）产生偏差。残余模态是用于修正“模态展开截断误差”的。实际工作中应始终启用RESVEC。
FRF的输出
FRF（频率响应函数）的主要类型：
| FRF类型 | 定义 | 用途 |
|---|---|---|
| 柔度 | $u/F$ | 位移响应 |
| 导纳 | $v/F = i\omega \cdot u/F$ | 速度响应 |
| 惯性率 | $a/F = -\omega^2 \cdot u/F$ | 加速度响应 |
实验中用哪种？
实验模态分析中一般使用惯性率（加速度/力）。因为加速度传感器使用最广泛。将FEM结果与实验对比时，应以相同格式输出。
总结
我来整理一下模态法频率响应的数值方法。
要点：
- 模态法比直接法快得多 — 在大量频率点计算时优势明显
- 用残余模态（RESVEC）进行高阶修正 — 应始终启用
- FRF格式 —
相关主题
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告