複素固有値解析
理论与物理
什么是复特征值
老师，“复特征值”和通常的特征值有什么不同？
通常的特征频率分析（实特征值分析）处理的是无阻尼自由振动。特征值是实数（$\omega^2$），可以得到固有频率和模态振型。
复特征值分析处理的是包含阻尼的自由振动。特征值变为复数（$\lambda = \sigma + i\omega_d$），实部 $\sigma$ 表示阻尼（稳定性），虚部 $\omega_d$ 表示带阻尼的振动频率。
控制方程
带阻尼的特征值问题：
$\lambda$ 和 $\{\phi\}$ 变为复数。
复数的特征值…物理意义是什么？
$\lambda = \sigma + i\omega_d$ 中：
- $\omega_d$ = 带阻尼振动频率（比实固有频率 $\omega_n$ 稍低）
- $\sigma$ = 衰减率。$\sigma < 0$ 表示稳定（振动衰减），$\sigma > 0$ 表示不稳定（振动增大）
- $\zeta = -\sigma / \sqrt{\sigma^2 + \omega_d^2}$ = 阻尼比
需要复特征值的场景
| 场景 | 理由 |
|---|---|
| 刹车尖叫（squeal） | 摩擦引起的自激振动。检测不稳定特征值（$\sigma > 0$） |
| 颤振分析 | 气动弹性的不稳定性。速度升高时 $\sigma$ 转为正值 |
| 非比例阻尼 | 阻尼矩阵无法进行模态正交化的情况 |
| 旋转机械 | 陀螺效应导致阻尼矩阵非对称 |
刹车尖叫是典型的应用场景吗。
刹车片和刹车盘的摩擦导致刚度矩阵变为非对称。非对称刚度会使系统不稳定，出现 $\sigma > 0$ 的特征值（自激振动模态）。这个模态就是刹车尖叫的原因。这是汽车NVH开发中最重要的分析之一。
Nastran
```
SOL 107 $ 复特征值分析
CEND
CMETHOD = 10
BEGIN BULK
EIGC, 10, HESS, , , , , 20
```
Abaqus
```
*STEP
*COMPLEX FREQUENCY
20, ,
*END STEP
```
Nastran的SOL 107，Abaqus的COMPLEX FREQUENCY呢。
设置和通常的特征值分析类似，但必须定义阻尼矩阵 $[C]$。如果没有阻尼，结果会和实特征值分析相同。
总结
我来整理一下复特征值分析。
要点：
- 包含阻尼的特征值问题 — 特征值为复数 $\lambda = \sigma + i\omega_d$
- $\sigma > 0$ 表示不稳定（自激振动） — 用于检测刹车尖叫、颤振
- 对应非比例阻尼・非对称刚度 — 实特征值分析无法处理
- SOL 107（Nastran）, *COMPLEX FREQUENCY（Abaqus）
- 主要应用于刹车NVH和航空航天领域的颤振
如果说实特征值分析是“了解结构的固有频率”，那么复特征值分析就是“确认结构是否稳定”的分析吧。
正是如此。复特征值分析是稳定性分析，是实特征值分析的扩展。
刹车尖叫的不稳定机制
盘式刹车的“尖叫”是由复特征值的不稳定（正实部）引起的。摩擦力增大了结构的振动能量，当实部为正时振动发散。这种现象与颤振不稳定性本质上是相同的机制，20世纪40年代气动颤振研究的经验在20世纪90年代被应用于刹车分析。
各项的物理意义
- 惯性项（质量项）：$\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统一 |
数值解法与实现
复特征值求解器
复特征值是怎么求解的？
需要与实特征值的Lanczos法不同的算法。
| 方法 | 特点 | 求解器 |
|---|---|---|
| Hessenberg法 | 小〜中规模。求全部特征值 | Nastran EIGC(HESS) |
| QZ法 | 广义特征值问题。稳定 | LAPACK |
| 投影法 | 投影到实模态后再复数化 | Abaqus COMPLEX FREQUENCY |
| Arnoldi法 | 大规模稀疏矩阵。Lanczos的非对称版 | 研究用 |
Abaqus的“投影法”是什么？
首先通过实特征值分析求得 $N$ 个模态，然后求解投影到实模态空间后的小规模复特征值问题。$N \times N$ 小矩阵的复特征值可以用QZ法求解。也能对应大规模问题。
实模态成为“基底”了呢。
所以在进行复特征值分析之前，必须先求得足够数量的实模态。模态数不足会导致复特征值的精度下降。
刹车尖叫分析的设置
刹车尖叫的复特征值分析步骤：
1. 包含摩擦接触的非线性静力分析 — 求解刹车的紧固状态
2. 接触面的摩擦力→线性化 — 由摩擦力构成非对称刚度矩阵
3. 复特征值分析 — 搜索不稳定特征值（$\sigma > 0$）
4. 确定不稳定模态 — 尖叫的频率和模态振型
摩擦的非对称刚度是不稳定性的原因呢。
摩擦力是随动载荷（方向随位移变化而变化的力），它使 $[K]$ 变为非对称。非对称的 $[K]$ 有可能“注入能量”，从而出现 $\sigma > 0$ 的不稳定模态。
总结
我来整理一下复特征值的数值方法。
要点：
- 投影法（实模态→复特征值）是Abaqus的标准 — 实模态数量决定精度
- Hessenberg法/QZ法是Nastran的标准 — 小〜中规模
- 刹车尖叫 — 摩擦的非对称刚度→通过复特征值检测不稳定模态
- $\sigma > 0$ 的模态不稳定（自激振动） — 通过设计变更排除
なった
詳しく
報告