2. 频响分析数学基础
OptiStruct 的频响求解器并不是一个“黑箱”，它内部至少提供了两条求解路径：直接法（Direct Method）与模态叠加法（Modal Superposition / Modal Frequency Response）。选择哪条路径、采用哪种阻尼模型，直接决定了计算规模、精度和适用频段。理解这三者之间的数学关系，是后续设置卡片、判断结果可信度的前提。
一、频响问题的基本方程
结构在简谐激励下的运动方程为：
[M]{x''} + [C]{x'} + [K]{x} = {F(ω)}·e^(iωt)设稳态解为 {x} = {X(ω)}·e^(iωt)，代入后得到频域代数方程：
(-ω²[M] + iω[C] + [K]){X(ω)} = {F(ω)}其中 Z(ω) = -ω²[M] + iω[C] + [K] 称为动刚度矩阵。频响函数（FRF）定义为输出与输入的比值，位移频响为 H(ω) = {X(ω)}/{F(ω)} = Z(ω)⁻¹。OptiStruct 输出的位移、速度、加速度、应力响应，本质上都是在这个 H(ω) 上乘以激励谱得到的。
二、直接法：逐频点求解
直接法不对方程做任何模态变换，而是在每个激励频率 ωi 上直接组装并求解复数线性方程组。它的核心特点是：
- 保留完整的 [M]、[C]、[K] 矩阵，不截断模态，因此对高频密集模态、局部模态、非比例阻尼都能给出完整响应；
- 每个频点一次复数分解，计算量随频点数线性增长，随自由度规模近似三次增长；
- 适合中小模型、宽频带、需要精确捕捉局部共振或声固耦合的场景。
在 OptiStruct 中，直接法通过 METHOD(STRUCTURE) 卡片选择，配合 FREQUENCY 卡片定义频点。若频点很多而模型很大，直接法的代价会迅速上升，此时应转向模态叠加法。
三、模态叠加法：用模态坐标降维
模态叠加法先求解无阻尼特征值问题：
([K] - ω_n²[M]){φ} = 0得到前 N 阶模态频率 ωn 与振型 {φn}。将物理坐标变换到模态坐标 {x} = [Φ]{q}，并左乘 [Φ]ᵀ，利用振型关于 [M]、[K] 的正交性，方程解耦为 N 个单自由度方程：
(-ω²m_n + iωc_n + k_n) q_n(ω) = f_n(ω)其中 mn、kn 为模态质量与模态刚度，fn = {φn}ᵀ{F} 为模态力。每个模态的响应可解析写出，再叠加回物理坐标。模态叠加法的关键假设是：响应可由有限阶模态线性组合近似。因此必须保证截断频率足够高，通常要求覆盖激励频带上限的 1.5～2 倍，并检查模态有效质量占比。
四、阻尼模型：决定共振峰的高度
阻尼是频响分析中最敏感的参数，因为它直接控制共振峰的幅值。OptiStruct 支持多种阻尼描述：
| 阻尼类型 | 数学形式 | 典型卡片 | 适用场景 |
|---|---|---|---|
| 结构阻尼（模态） | [C] = (G/ω)·[K]，G 为结构阻尼系数 | PARAM, G | 金属结构、焊接件，常用 G=0.02～0.05 |
| 粘性阻尼（模态） | [C] = 2ζωn[M] | TABDMP1 / SDAMPING | 含橡胶、阻尼器的系统 |
| 瑞利阻尼 | [C] = α[M] + β[K] | PARAM, ALPHA / BETA | 需要与实验模态匹配的场合 |
| 单元阻尼 | 材料或单元级损耗因子 | MAT1 的 GE 字段 | 复合材料、粘弹性材料 |
结构阻尼与粘性阻尼在共振点附近的等效关系为 G ≈ 2ζ。例如 G=0.04 对应 ζ≈0.02，共振放大因子约 1/(2ζ)=25。若把阻尼设得过小，共振峰会被严重高估；设得过大，则可能掩盖真实危险频率。工程上常先用 G=0.02～0.03 做初步评估，再根据试验模态或经验修正。
五、两条路径的对比与选择
下图给出直接法与模态叠加法在同一频响问题上的求解路径、输入输出与适用条件对比，帮助在建模阶段做出选择。
选择原则可以概括为：模型小、频带宽、阻尼复杂、需要精确局部响应 → 直接法；模型大、频带窄、模态稀疏、需要反复优化 → 模态叠加法。在 OptiStruct 优化循环中，模态叠加法通常更快，因为模态可以在设计变量变化后通过近似更新，而不必每次重新求解完整复数方程。
六、频响函数的读图要点
频响函数曲线上的每一个峰对应一个共振频率，峰的高度由阻尼决定，峰的宽度由阻尼比决定。半功率带宽法给出阻尼比估计：ζ ≈ Δω/(2ωn)，其中 Δω 为幅值降至峰值的 1/√2 处的频率带宽。相位曲线在共振点穿越 -90°，反共振点穿越 +90°。若相位曲线出现异常跳变或幅值曲线出现非物理尖峰，往往提示模态截断不足或阻尼设置不当。
掌握上述数学基础后，后续在 OptiStruct 中设置 METHOD、FREQUENCY、TABDMP1、PARAM,G 等卡片时，就能清楚每个参数在方程中的位置，以及它如何影响最终频响曲线的形状与幅值。