直接法频率响应分析
直接法频率响应的理论基础
什么是直接法
老师，相对于模式法，"直接法"频率响应分析有什么区别？
直接法是不进行固有模式展开，而是在各频率处直接求解运动方程。
这是频率 $\omega$ 处的复数联立方程。
在各频率处求解联立方程，所以计算量很大？
是的。需要求解 $n \times n$（$n$ = 自由度数）的复数联立方程，次数等于频率点数。模式法只需要求解 $N$ 个一自由度标量方程，所以速度快得多。
需要直接法的场景
那么直接法有什么用处呢？ 模式法无法准确处理的情况需要用直接法：
| 情况 | 原因 |
|---|---|
| 非比例阻尼 | 阻尼矩阵无法通过模式正交化 |
| 频率相关的材料特性 | 粘弹性材料。$E(\omega), \eta(\omega)$ |
| 结构阻尼（滞后） | 复刚性 $K^* = K(1+ig)$ |
| 大阻尼系统 | 橡胶隔振垫、制振材 |
| 外部阻抗边界 | 土壤-结构耦合等 |
对于"频率改变材料特性"的粘弹性材料，模式法无法使用。
模式法以固有模式（频率无关）为基础，所以频率相关的材料特性无法自然融入模式展开。直接法可以在各频率处更新材料特性。
结构阻尼的处理
结构阻尼（滞后阻尼）在直接法中最容易处理：
$g$ 是结构阻尼系数。这是频率无关的阻尼，在许多情况下比粘性阻尼更物理准确。
结构阻尼在时间域中不能使用吧。
完全正确。结构阻尼仅在频率域（直接法）中具有物理意义。在时间域中使用结构阻尼会破坏因果律。
Nastran
```
SOL 108 $ 直接法频率响应
CEND
FREQUENCY = 20
BEGIN BULK
FREQ1, 20, 1., 500., 1.
```
Abaqus
```
*STEP
*STEADY STATE DYNAMICS, DIRECT
1., 500., 500, 1.
*END STEP
```
Ansys
```
/SOLU
ANTYPE, HARMONIC
HROPT, FULL ! 直接法
HARFRQ, 1., 500.
NSUBST, 500
SOLVE
```
总结
总结一下直接法频率响应。
要点：
- 在各频率处直接求解联立方程 — 无模式展开
- 计算成本为模式法的10〜100倍 — 频率点数×自由度数
- 支持非比例阻尼、频率相关材料、结构阻尼 — 超越模式法的极限
- SOL 108（Nastran）, *SSD DIRECT（Abaqus）, HARMONIC FULL（Ansys）
- 大多数问题用模式法就足够 — 直接法仅用于特殊情况
直接法的矩阵规模是自由度的3倍
直接法（Direct Method）调和响应分析在每个频率步骤处对复刚性矩阵 [K + iωC − ω²M] 进行逐次LU分解。将矩阵实部和虚部分离后，有效自由度变为2倍，进一步考虑LU分解的填充（fill-in），所需内存容量是理论自由度的3〜5倍。在100万自由度模型中，每个频率点需要数分钟，因此稀疏求解器（如PARDISO）的选择改变了计算时间的数量级。
直接法频率响应的数值计算方法
直接法的计算效率化
有办法降低直接法的计算成本吗？
有几种方法：
1. 动态刚性矩阵的LU分解重用
$[D(\omega)] = -\omega^2[M] + i\omega[C] + [K]$ 的LU分解是最耗时的。频率无关部分（$[K]$）只分解一次，频率相关部分作为增量处理的迭代法。
2. 并行计算
各频率点的计算相互独立，所以频率点间完全并行。100个频率点用100核同时计算，实际上只需要一个频率点的计算时间。
3. 约简法（Reduced Method）
在NastranSOL 108中，可以结合Guyan约简或CMS约简来减少自由度后再应用直接法。与模式法不同，但通过减少自由度降低效率。
粘弹性材料的建模
在直接法中如何处理粘弹性材料？
粘弹性材料的复弹性模量：
$E'$ 是储存弹性模量（刚度），$E''$ 是损失弹性模量（阻尼），$\eta = E''/E'$ 是损失因子。
在Abaqus中用 *VISCOELASTIC, FREQUENCY 定义Prony级数参数。自动计算各频率的 $E^*(\omega)$。在Nastran中用TABLEM1定义频率相关材料。
在制振材（约束层阻尼）设计中，这种频率依赖性很重要。
橡胶和粘弹性聚合物的损失因子 $\eta$ 对频率和温度非常敏感。不用直接法无法准确处理这种依赖性。
总结
总结一下直接法的数值方法。
要点：
- 频率点间并行计算实现效率化 — 完全并行可行
- 粘弹性材料的复弹性模量 — $E^*(\omega) = E'(1+i\eta)$
- 用Prony级数建模粘弹性 — Abaqus *VISCOELASTIC
- 与约简法结合 — 通过减少自由度降低计算成本
Householder在1958年整理了复数特征值
直接法数值基础的复矩阵三对角化算法由A.S. Householder在1958年发表。随后，20世纪60年代IBM System/360的EISPACK成为CAE求解器的公共库，MSC Nastran的SOL 108直接继承了这一传统。当今直接法求解器的核心算法本质上仍未脱离Householder的初心。
直接法频率响应的实务应用
直接法的实务应用
直接法在什么场合使用？
制振材设计（CLD: Constrained Layer Damping）
在铁板上粘贴粘弹性层+约束层以增加阻尼的制振材。粘弹性层的损失因子 $\eta(\omega, T)$ 随频率和温度变化，必须使用直接法。
土壤-结构耦合（SSI）
在结构界面施加土壤的阻抗（频率相关的刚度+阻尼）。土壤阻抗随频率变化显著，需要用直接法在各频率处反映土壤特性。
高阻尼结构
基础隔振结构或橡胶隔振垫等，阻尼比超过10%的结构。非比例阻尼很强，模式法的假设（比例阻尼→模式正交性）不成立。
实务检查清单
"直接法是否真正必要"是最初的判断。大多数问题可以用模式法解决。
直接法计算成本很大，必须在确实需要的情况下才能使用。判断标准是"阻尼是否比例"和"材料是否频率相关"这两个。
船舶振动分析以直接法为标准
船舶推进轴系的调和响应分析中，2〜50Hz的宽频带内有多个密集共振，因此比模式法更推荐使用直接法。三菱造船在2018年发布的7万吨级LNG船轴系分析采用直接法SOL 108计算1〜30Hz、200个频率点。检出了螺旋桨叶片通过频率（BPF）在轴系固有值的1.5Hz内接近的3个案例，并进行了设计变更。
直接法频率响应的软件对比
直接法的工具
直接法求解器的对比？
Nastran的DMIG是什么？
DMIG（Direct Matrix Input at Grid Points）是在Nastran中直接输入外部计算的刚度/阻尼矩阵的功能。通过DMIG输入土壤阻抗矩阵，再用直接法求解，是SSI分析的标准做法。
选择指南
Nastran SOL 108和Abaqus Steady-State Dynamics
直接法调和响应的两大实现是Nastran SOL 108和Abaqus `*STEADY STATE DYNAMICS, DIRECT`。Nastran源于20世纪60年代NASA，强调可靠性；Abaqus成立于1978年的HKS公司，源于研究用途。2023年NAFEMS PB17 Benchmark对比中，Abaqus由于稀疏求解器的优化，同模型比Nastran快约40%，但在模型互换性上Nastran仍占优势。
直接法频率响应的先端研究
直接法的先端研究
请介绍直接法的前沿。
模型阶数约简（MOR）
为了降低直接法的成本，模型阶数约简（MOR）技术受到研究关注。使用Krylov子空间法或有理近似，用少数自由度近似原始模型的频率响应。
这与模式法不同的约简技术？
模式法用固有模式展开，而MOR用Krylov子空间或Padé近似展开。对于包含频率相关材料的系统，MOR也能适用，是直接法的代替品。
非线性频率响应
包含接触（间隙）或摩擦的非线性系统的频率响应。使用HBM（Harmonic Balance Method）逐谐波分量求解非线性方程。LS-DYNA的*FREQUENCY_DOMAIN功能支持这一点。
总结
总结一下直接法的先端研究。
直接法从"准确但缓慢"的手法，通过MOR正在进化为"准确且高效"。
频率相关阻尼仅直接法能处理
粘弹性材料的频率相关损失因子（例如天然橡胶在10Hz时为0.05，1000Hz时为0.3）模式法无法处理，直接法是必须的。Ansys Harmonic Analysis的Full Method中可直接给`MP,DMPR`赋频率表。汽车仪表盘的绿色噪声（200〜800Hz）的准确预测需要这种频率相关阻尼，是直接法应用的重要原因之一。
直接法频率响应的故障排除
直接法的故障
直接法常见的故障？
计算太慢
应对措施：
误在时间域使用结构阻尼
结构阻尼（$g$）仅供频率域使用。在时间域使用会破坏因果律，产生非物理响应。
应对：在时间域切换为Rayleigh阻尼或模式阻尼。
粘弹性材料Prony级数不准确
Prony级数参数与实测 $E'(\omega), \eta(\omega)$ 不符时，需要重新拟合。
应对：
- 在Prony级数拟合的频率范围内验证拟合精度
- 使用足够的项数（5〜10项）
- Abaqus的*VISCOELASTIC, TEST DATA可直接输入实测数据
总结
总结一下直接法的故障排除。
内存不足用频率分割回避
直接法在100万自由度以上模型求解时，内存不足错误源于各点矩阵因子分解。回避办法是将频率带分割为多个作业，各带的结果通过后处理汇总。ANSYS Mechanical 2022R2之后有`Harmonic Frequency Range`自动分割集群投入功能，在64GB内存环境下也能完成500万自由度的计算。
参考
细节
错误