3. 响应谱分析基础理论
响应谱分析的核心思想，是把一条随时间变化的瞬态激励（地震加速度、冲击脉冲、爆炸载荷等）转换成一条“峰值—频率”曲线。这条曲线不再关心时间轴上的细节，而是回答一个工程上更关键的问题：对于某个固有频率和阻尼比的单自由度系统，在这条激励下它的最大响应是多少？理解了单自由度系统的响应谱计算原理，才能正确理解 ANSYS 中响应谱分析的输入、组合方式与结果含义。
一、响应谱的定义与物理含义
响应谱（Response Spectrum）描述的是：给定一条基础激励加速度时程 a(t)，对一系列不同固有频率 fn（或圆频率 ωn）、相同阻尼比 ζ 的单自由度系统，分别求出各自的最大绝对响应，再以频率为横轴、峰值为纵轴绘制的曲线。按响应量不同，可分为：
- 绝对加速度谱
Sa(fn, ζ)：常用于等效静力法与惯性力计算； - 相对速度谱
Sv(fn, ζ)：与伪速度谱近似，是三联谱的中间量； - 相对位移谱
Sd(fn, ζ)：用于评估相对变形与间隙； - 伪加速度谱
Spa = ωn2Sd：在低频段与绝对加速度谱接近，便于三联谱统一表达。
需要强调：响应谱是单自由度系统的峰值集合，不是某个多自由度结构的真实响应。它把“激励”和“结构”解耦，因此同一地震波可以生成无数条谱（不同阻尼比各一条），而同一结构在不同地震波下响应也完全不同。
二、单自由度系统响应谱的计算原理
设单自由度系统质量 m、刚度 k、阻尼 c，基础受到加速度激励 ag(t)。以相对位移 u = x - xg 为未知量，运动方程为：
m·ü + c·u̇ + k·u = -m·a_g(t)
即 ü + 2ζω_n·u̇ + ω_n²·u = -a_g(t)
其中 ω_n = √(k/m)，ζ = c / (2√(km))求解该方程得到 u(t)、u̇(t)、üabs(t) = ü + ag，再取时间历程上的最大值，即得到该频率点对应的谱值。工程上常用两种数值方法：
- Duhamel 积分（卷积法）：把激励离散为一系列脉冲，用单位脉冲响应函数
h(t)卷积叠加，适合激励采样点较少的情况； - 直接时程积分（Newmark-β 法）：对每个频率点独立积分，稳定性好、精度可控，是 ANSYS 及多数商用软件生成谱曲线的默认做法。
计算流程可概括为：对频率轴上的每个 fn,i，用同一阻尼比 ζ 求解一次 SDOF 时程，记录峰值，最后连成曲线。频率分辨率、时间步长与阻尼比是影响谱值精度的三个关键参数。
三、地震响应谱与冲击响应谱的区别
虽然两者都基于 SDOF 时程计算，但工程用途和谱形特征差异明显：
| 对比项 | 地震响应谱 | 冲击响应谱 |
|---|---|---|
| 激励来源 | 地震地面加速度记录（如 El Centro、Taft 波） | 机械冲击、爆炸、跌落、火工品分离等脉冲 |
| 持续时间 | 数十秒，含多周期循环 | 毫秒至数十毫秒，单次或少数几次脉冲 |
| 谱形特征 | 低频段上升、中频平台、高频段下降 | 随频率升高呈振荡衰减，峰值集中在低频 |
| 阻尼取值 | 常用 5%（建筑）、2%~5%（设备） | 常用 0%~5%，常取 0 得到保守包络 |
| 典型用途 | 建筑抗震、管道支架、设备抗震鉴定 | 电子设备抗冲击、航天器火工冲击环境 |
地震谱通常以标准设计谱形式给出（如 GB 50011 的 αmax 与特征周期 Tg），而冲击谱往往由实测脉冲或半正弦、后峰锯齿等理想脉冲生成。在 ANSYS 中，两者都通过 SPECTR 命令族输入，但阻尼比、频率范围与组合方式需按各自规范选取。
四、ANSYS 中响应谱分析的输入与结果判读
在 ANSYS 中执行响应谱分析，典型流程是：先做模态分析提取固有频率与振型，再通过 SPECTR 定义谱类型（位移、速度、加速度）、谱值与频率的对应关系、阻尼比，最后用 SRSS、CQC 或 ROSE 等方法组合各阶模态响应。关键输入包括：
- 谱值表：频率—谱值成对输入，频率范围应覆盖结构前若干阶固有频率；
- 激励方向：可单方向、双方向或三方向同时施加，方向组合常用 100%+40%+40% 或 SRSS；
- 阻尼比：与谱曲线对应的阻尼比必须一致，否则谱值失真；
- 模态截断：通常取前 10~30 阶或累计质量参与系数达 90% 以上。
理解响应谱的定义、SDOF 计算原理以及地震谱与冲击谱的差异，是后续在 ANSYS 中正确设置谱值、阻尼、方向组合与模态截断的前提。只有输入与规范一致、组合方法选择合理，响应谱分析的结果才具有工程可信度。