频率扫描与共振评估
理论与物理
什么是频率扫描
老师，什么是“频率扫描”（frequency sweep）？
这是一种连续改变外力频率，从低频到高频，并计算响应的方法。是获取FRF（频率响应函数）的基本方法。
简单来说就是“从1 Hz到500 Hz按顺序计算”吗？
是的。绘制每个频率下的响应（位移、加速度、应力），可以一目了然地看出共振峰的位置和大小。
频率步长设计
频率步长的确定取决于共振的性质：
| 阻尼比 $\zeta$ | 共振峰半功率带宽 | 所需步长 |
|---|---|---|
| 0.1%（极低阻尼） | $\Delta f \approx 0.002 f_n$ | 0.1 Hz以下 |
| 1%（钢结构） | $\Delta f \approx 0.02 f_n$ | 1 Hz左右 |
| 5%（RC结构） | $\Delta f \approx 0.1 f_n$ | 5 Hz左右 |
| 10%（隔震） | $\Delta f \approx 0.2 f_n$ | 10 Hz左右 |
阻尼越小，峰值越尖锐，所以需要更小的步长呢。
低阻尼钢结构需要非常精细的步长。固有频率在100 Hz附近且 $\zeta = 0.5\%$ 时，步长需小于0.5 Hz。如果扫描500个点，从1到500 Hz，步长1 Hz，这样可能会漏掉峰值。
对数步长与模态追踪步长
高效的步长设置方法：
1. 对数等间距
在对数尺度上等间距设置频率点。低频粗，高频细。声学系统中常用。
2. 模态追踪步长（Nastran的FREQ4）
在每个固有频率附近自动配置精细步长，确保捕捉到共振区域。通过Nastran的FREQ4卡片设置。
3. 自适应步长
使用Abaqus的BIAS参数在共振区域附近自动加密步长。指定计算点数后，会自动集中分配到共振区域附近。
共振评估
共振的评估指标：
| 指标 | 定义 | 用途 |
|---|---|---|
| 共振频率 | FRF的峰值位置 | 共振规避设计 |
| 峰值振幅 | FRF的最大值 | 响应最大值评估 |
| 半功率带宽 | 峰值-3 dB处的宽度 | 阻尼比估计 |
| 相位变化 | 共振处约180°变化 | 模态确认 |
FRF的峰值振幅在设计中最重要吗？
是的。峰值振幅 × 输入力 = 最大响应。判断这个最大响应是否在允许值（位移限值、加速度限值、应力限值）以内，是设计的关键。
总结
我来整理一下频率扫描和共振评估。
要点：
- 改变频率获取FRF — 识别共振峰
- 步长小于半功率带宽 — $\Delta f < \zeta \cdot f_n$
- 使用模态追踪步长（FREQ4）提高效率 — 仅在共振区域附近加密
- 峰值振幅 × 输入 = 最大响应 — 设计判断基准
- 通过相位变化确认共振 — 180°的相位跳变
扫描速度改变，共振也会改变
在频率扫描试验中，扫描速度（倍频程/分钟）越快，共振峰出现“虚假”现象的问题就越明显。正如CW. de Silva在1950年代指出的那样，会产生依赖于扫描速度的、相对于真实共振频率的表观偏移。目前的振动试验标准MIL-STD-810H推荐使用4倍频程/分钟以下，正是基于这一理论背景。
各项的物理意义
- 惯性项（质量项）：$\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统一 |
数值解法与实现
频率步长设置
请告诉我各求解器中频率步长的设置方法。
Nastran
```
$ 等间距
FREQ1, 20, 1., 500., 1. $ 1〜500 Hz, 1 Hz步长
$ 模态追踪（共振区域附近自动加密）
FREQ4, 20, 1., 500., 0.1, 5 $ 各模态的±0.1半功率带宽, 5个点
$ 对数等间距
FREQ2, 20, 1., 500., 10 $ 1〜500 Hz, 1/3倍频程
```
Abaqus
```
*STEADY STATE DYNAMICS
1., 500., 500, 1. $ 1〜500 Hz, 500个点, BIAS=1（等间距）
```
BIAS > 1 则向高频集中，BIAS < 1 则向低频集中。
Ansys
```
HARFRQ, 1., 500.
NSUBST, 500 ! 500个步
```
Nastran的FREQ4是最智能的步长设置方法呢。
FREQ4会自动在每个固有频率附近集中计算点。有时能以等间距1/10的效率获得同等的精度。实际工作中，组合使用FREQ1（粗略的整体步长）+ FREQ4（共振区域附近的精细步长）最为有效。
FRF的输出与可视化
FRF的标准显示形式：
- 振幅-频率图 — 通常使用对数刻度（dB）
- 相位-频率图 — 共振处有180°跳变
- 奈奎斯特图（实部 vs. 虚部） — 共振处描绘圆形
- 伯德图 — 分两段显示振幅和相位
dB刻度是什么？
$20 \log_{10}(|H|/H_{ref})$。将振幅的巨大变化压缩显示，使其易于观察。共振峰显示为+40 dB，反共振显示为-40 dB等。
总结
我来整理一下频率扫描的数值方法。
要点：
- FREQ1 + FREQ4（Nastran）效率最高 — 等间距 + 模态追踪
- 使用BIAS（Abaqus）进行集中分配 — 自动向关注频段集中
- 用dB刻度显示FRF — 压缩宽广的动态范围
- 用奈奎斯特图确认共振 — 共振处描绘圆形
对数扫描与线性扫描的区分使用
低频带域（1〜100Hz）通常使用对数扫描（倍频程恒定），为每个频率分配均等的时间。线性扫描用于高频（1kHz以上）的精密测量或电气特性评估。汽车振动台试验（ISO 16750-3）规定标准为对数扫描1倍频程/分钟，5〜2000Hz。即使是相同的试验，对数与线性扫描的共振检测精度最大可相差3倍。
线性单元（一阶单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切自锁（可通过减缩积分或B-bar法缓解）。
二阶单元（带中间节点）
可以表现曲线变形。应力精度大幅提高，但自由度约增加2〜3倍。推荐：应力评估重要时使用。
完全积分 vs 减缩积分
完全积分：有过约束（自锁）风险。减缩积分：有沙漏模式（零能量模式）风险。根据情况选择合适的方法。
自适应网格
基于误差指标（ZZ估计量等）的自动细化。高效提高应力集中区域的精度。有h法（单元细分）和p法（增加阶次）两种。
牛顿-拉夫森法
非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性，但计算成本高。
修正牛顿-拉夫森法
切线刚度矩阵使用初始值或每隔几次迭代更新。每次迭代成本低，但收敛速度为线性。
收敛判定标准
力残差范数: $||R|| / ||F_{ext}|| < \epsilon$（通常 $\epsilon = 10^{-3}$〜$10^{-6}$）。位移增量范数: $||\Delta u|| / ||u|| < \epsilon$。能量范数: $\Delta u \cdot R < \epsilon$
载荷增量法
不一次性施加全部载荷，而是分小步增加。弧长法（Riks法）可以超越载荷-位移关系的极值点进行追踪。
直接法 vs 迭代法的比喻
直接法是“用笔算精确解联立方程”的方法——可靠但大规模问题耗时过长。迭代法是“反复猜测逼近正确答案”的方法——最初答案粗糙，但每次迭代精度都会提高。就像查字典时，从第一页开始顺序查找（直接法）不如先估计位置翻开，再前后调整（迭代法）来得高效，原理相同。
网格阶次与精度的关系
なった
詳しく
報告