模态法频率响应分析
模态法频率响应的理论基础
频率响应分析是什么
老师，"频率响应分析"是什么？
调和（正弦波）外力的结构定常响应求解分析。通过改变外力的振动数，计算各振动数下的位移、加速度和应力。
输入：$\{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个独立的单自由度方程，计算速度得以大幅提升。
模态法频率响应的数值计算手法
模态法的计算效率
请说明模态法为什么比直接法快。
直接法每个振动数点都要解 $n \times n$（$n$ = DOF数）的联立方程。模态法是特征值分析1次 + 每个振动数点解 $N \times N$（$N$ = 模态数 << $n$）的对角系统。
| 计算量 | 直接法 | 模态法 |
|---|---|---|
| 每个频率点 | $O(n \cdot bw)$ or $O(n^2)$ | $O(N)$ |
| 频率点数 $M$ | $M \times O(n \cdot bw)$ | 特征值 + $M \times O(N)$ |
如果 $N = 100$ 模态而 $n = 1{,}000{,}000$ DOF，模态法快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形式 — 柔度、流动性、惯纳
- 与实验对比 — 输出同一FRF形式
分析带宽的1.5倍模态是"保险牌"
为了确保模态叠加的精度，通常的经验法则是所关心频率带宽上限的1.5～2倍才包括所有模态。这个"1.5倍法则"在NASA-STD-5002（1996年）中明文规定，后来也被纳入NAFEMS建议。例如要求100Hz以下的响应，就应该提取150Hz以下的所有模态。模态数不足的话高频带宽的误差可能超过30%。
模态法频率响应的实务应用
频率响应分析的实务
频率响应分析在实务中怎么用？
汽车的NVH
发动机的旋转振动加振车身，室内听到的是噪音。用频率响应分析评估"发动机支架的输入→车室内音压"的传递特性（NTF: Noise Transfer Function）。
机器振动评估
旋转机械的不平衡力（$F = m \cdot e \cdot \omega^2$，频率 = 转速）的响应。与ISO 10816等振动标准对比。
管道振动疲劳
泵的脉动（压力波动）引起的管道振动。共振会导致应力过大，引起疲劳破坏。
频率间隔的设定
频率的间隔（步长）怎么设定？
为了准确捕捉共振峰，间隔应该小于峰值半宽：
$\zeta = 1\%$、$f_n = 100$ Hz时，$\Delta f < 1$ Hz。
阻尼小时峰值会很尖锐，所以需要细间隔。
作为替代，可以用对数等间隔（低频粗，高频细）或模态自动间隔（共振附近细化）。Nastran的FREQ4或Abaqus的BIAS功能。
实务检查清单
"阻尼设置"和"频率间隔"最重要。
没有阻尼时共振振幅是无穷大。间隔太粗就漏掉峰值。这两个错误会让结果完全无意义。
桥梁风励振分析通常30个模态足够
吊桥的风励振（颤振与涡激振动）的调和响应分析中，通常下位30～50个模态就能覆盖全部振动能的99%以上。1940年Tacoma Narrows大桥倒塌后，美国联邦公路局制定的桥梁耐风设计标准（AASHTO）基于模态叠加法，这一方法至今仍保留在现代设计规范中。
模态法频率响应的软件比较
频率响应工具
哪些工具可以做频率响应分析？
Nastran的SOL 111在NVH中是标准吗？
汽车NVH中SOL 111是事实上的标准。AMLS（超大规模特征值）→SOL 111（模态法频率响应）的工作流最有效率。
选择指南
大多数情况下模态法就够了。
直接法只在"非比例阻尼""频率相关的材料特性""粘弹性材料"这样的特殊情况才需要。通常模态法足够。
MSC Nastran与Ansys模态法精度差<1%
NAFEMS基准R0016（2019年）报告显示，MSC Nastran SOL 111与Ansys Mechanical的模态法调和响应精度差在17个验证工况中全部<1%。另一方面计算时间强烈依赖模型规模和求解器设置，10万DOF以下基本相同，100万DOF以上Ansys的BlockLanczos快约20%。
模态法频率响应的先端研究
频率响应的先端研究
请介绍频率响应分析的最前沿。
基于FRF的模型更新
为了使实验测得的FRF与FEM的FRF一致，通过优化调整FEM的材料特性、阻尼、接合部刚度等。同时匹配FRF的振幅和位相。
概率频率响应
考虑材料特性和阻尼的变异性，计算FRF的信任区间。"99%概率下FRF在这个范围内"的预测。
非线性频率响应
包含接触（间隙）或非线性阻尼的非线性频率响应。用HBM（调和平衡法）或Shooting法求定常响应。LS-DYNA的*FREQUENCY_DOMAIN_SSD等。
总结
总结频率响应的先端研究。
模态截断补偿（残余柔度）的威力
模态叠加中高次模态截断的误差可以通过"残余柔度（Residual Flexibility）"补偿大幅降低。这一手法由R.L. Bisplinghoff在1971年提唱，现在Nastran SOL 111的`RESVEC`选项中实现。未补偿时5%误差的模型经补偿后可降至0.3%以下，这已在NAFEMS基准W4中公布。
模态法频率响应的故障排除
频率响应故障
频率响应分析的常见故障有哪些？
共振峰不出现或非常大
反共振与实验不符
反共振（FRF的谷）的位置对输入输出点的位置很敏感。要准确对齐FEM的输出点（节点位置）与实验的加速度传感器位置。
高频响应异常
模态数不足。高频（在关心范围的上限附近）时模态贡献不够。加入残余模态（RESVEC）或增加模态数。
由频率间隔漏掉峰值
间隔太粗会通过尖锐的共振峰。参考 $\Delta f < \zeta \cdot f_n$ 来设定间隔。
总结
总结频率响应故障处理。
共振峰值过低通常是模态数不足
模态叠加法中共振峰值远低于实测时，首要原因是模态提取数不足。尤其是加振点附近有很多局部模态的结构（焊接部或薄板面板），提取上限翻倍后最大响应可能增加3倍。确认步骤：①查看有效质量累计90%达成的模态数，②检查响应点的有效质量贡献在上位5个模态中是否超过70%。
价值
详细
错误