3. 响应谱分析的数学基础

响应谱分析之所以能在一次求解中同时考虑结构所有模态对地震或冲击载荷的响应,根源在于它把“时间历程”压缩成了“峰值—频率”关系。理解这条压缩链,是正确选择模态组合方法、正确解读结果的前提。本节从响应谱的定义出发,逐步过渡到设计谱、包络谱,再落到SRSS、CQC、ABS等组合方法的适用条件。

一、响应谱:从时程到峰值的压缩

设一个单自由度振子,质量m、刚度k、阻尼比ζ,固有圆频率ωn=√(k/m)。把它固定在刚性基础上,基础受到加速度时程 ag(t) 激励。振子的相对位移响应 u(t) 满足:

ü(t) + 2ζω_n·u̇(t) + ω_n²·u(t) = -a_g(t)

对同一个 ag(t),改变振子的 ωn(或周期 T=2π/ωn),就得到一族不同的 u(t)。取每条时程的绝对值最大值,就得到三个基本谱量:

  • 相对位移谱 Sd(T,ζ) = max|u(t)|
  • 相对速度谱 Sv(T,ζ) = max|u̇(t)|
  • 绝对加速度谱 Sa(T,ζ) = max|ü(t)+ag(t)|

在阻尼较小时,三者近似满足 Sa ≈ ωn·Sv ≈ ωn²·Sd,这就是著名的“三谱近似关系”。工程上最常用的是 Sa,因为它直接对应惯性力,便于与静力分析类比。Ansys Mechanical 的 Response Spectrum 分析中,用户输入的正是 Sa—频率(或周期)曲线,阻尼比则通过命令或阻尼设置单独给出。

需要强调:响应谱是“单自由度振子的最大响应”,不是结构本身的响应。它只携带幅值信息,丢失了相位。这一丢失正是后面必须用组合方法近似的原因。

二、设计谱与包络谱

实际地震记录千差万别,逐条计算不现实。规范做法是:对一批场地地震波分别计算谱,再取平均、平滑、放大,得到一条光滑的设计谱。设计谱通常由几个控制段组成:短周期平台段、下降段、长周期段,并规定阻尼修正系数。例如某场地设计谱在 T=0~0.1 s 为平台 Sa=0.4g,0.1~0.5 s 按 1/T 下降,0.5 s 以后按 1/T² 下降,阻尼比 5%。

包络谱则是把多条谱曲线逐点取最大值得到的上包络线。它比设计谱更“保守”,常用于校核:如果结构在包络谱下仍满足要求,则在任何单条谱下都满足。设计谱与包络谱的区别在于:设计谱是规范化的、可复现的;包络谱是面向具体波组的、偏安全的。

在 Ansys 中,谱曲线以“频率—加速度”对的形式输入,频率必须覆盖结构所有有效模态频率。若某阶模态频率落在谱曲线之外,软件会按端点值外推或报错,因此输入前应检查频率范围。

三、模态组合方法:SRSS、CQC、ABS

响应谱分析先对每一阶模态单独求最大响应 ri,再把这些峰值组合成总响应 R。由于各模态峰值不在同一时刻出现,组合本质上是统计估计。常用方法有三类:

  • ABS(绝对值求和):R = Σ|ri|。最保守,假设所有模态峰值同时同号出现。仅适用于模态频率非常接近、强耦合或需要极保守估计的场合。
  • SRSS(平方和开方):R = √(Σri²)。假设各模态响应互不相关、峰值随机分布。适用于模态频率分离良好(一般相邻频率比大于1.1~1.2)、阻尼较小的情况。这是最常用的默认方法。
  • CQC(完全二次组合):R = √(ΣΣ ρij·ri·rj)。引入模态相关系数 ρij,考虑模态间的耦合。当模态频率密集、阻尼较大时,SRSS会低估响应,必须用CQC。

相关系数 ρij 的经典表达式为:

ρ_ij = 8ζ²(1+r)·r^(3/2) / [(1-r²)² + 4ζ²r(1+r)²]
其中 r = ω_j/ω_i,ζ 为阻尼比

当 r 远离1(频率分离)时,ρij→0,CQC退化为SRSS;当 r→1(频率重合)时,ρij→1,CQC退化为ABS。因此CQC是SRSS与ABS之间的连续过渡,适用性最广,但计算量也最大。

Ansys Mechanical 中,模态组合方法通过分析设置选择:SRSS、CQC、ABS、ROSE、GRP等。对于大多数建筑与设备,若前若干阶模态频率比大于1.1,SRSS即可;若存在密集模态(如管道、机架),应改用CQC。阻尼比在CQC中直接影响相关系数,必须与谱曲线对应的阻尼一致,否则组合结果会失真。

四、方向组合与符号问题

地震通常三向输入。Ansys先对每个方向单独做模态组合,再把三个方向的结果组合。常用方向组合有:SRSS(三向平方和开方)、100%+30%+30%(主方向全值,另两方向取30%)、ABS等。方向组合与模态组合是两层独立的组合,不能混为一谈。

此外,响应谱分析只给出幅值,不给出符号。因此位移、应力结果都是“最大值包络”,不能直接用于判断拉压或正负。若需要符号信息,应改用时间历程分析或模态叠加法。

实用判断:先看模态频率分布,再选组合方法。频率比>1.1用SRSS;频率密集或阻尼>10%用CQC;需要极保守估计用ABS。方向组合按规范选取,不要默认全部SRSS。
响应谱分析数学链条:从时程到组合结果 地震加速度时程 a_g(t),含相位信息 SDOF求解 响应谱 S_a(T,ζ) 峰值—周期,丢失相位 规范化 设计谱 / 包络谱 平滑、保守、可复现 逐阶模态求峰值 r_i 各阶模态最大响应 r₁, r₂, r₃, … rₙ 每阶模态独立求解,峰值不同时出现,相位未知 组合方法决定如何由 r_i 估计总响应 R ABS 绝对值求和 R = Σ|r_i| 最保守 频率重合、强耦合 SRSS 平方和开方 R = √(Σr_i²) 默认常用 频率比 > 1.1,小阻尼 CQC 完全二次组合 R = √(ΣΣ ρ_ij r_i r_j) 适用最广 密集模态、大阻尼 方向组合:SRSS / 100%+30%+30% / ABS 结果仅幅值包络 无符号,不能判拉压