理论分享 | 响应谱分析（1）：模态叠加与模态组合的理论基础

响应谱分析（Response Spectrum Analysis）是一种基于模态叠加思想的结构动力学分析方法，用于评估结构在地震、冲击等外部激励作用下的峰值响应。该方法是模态分析结果与已知响应谱连接后进行模型最大响应计算（如位移、加速度、应力、应变和反力等）的一种分析方法。其通过模态叠加法将结构简化为多个等效单自由度系统（SDOF），计算每个模态在给定响应谱下的最大响应，再通过模态组合方法得到整体结构的最大响应，属于频域分析类型。

响应谱分析尤其适用于地震和冲击等动力作用下的响应评估：一方面，由于此类载荷的精确时间历程往往未知，难以直接采用瞬态动力学分析获取结构最大响应；另一方面，实际工程中通常更关注结构在该类载荷下的最大响应值（如最大位移、最大加速度、最大应力），而并不关注该值发生的具体时刻。因此，评估结构在给定载荷下的最大响应时，响应谱分析具有较高的计算效率且精度满足工程需求，可作为耗时较长的瞬态动力学分析（时间-历程分析）的有效替代，应用相对广泛。

本文主要介绍响应谱分析的理论基础，包括基本计算思想、响应谱的定义、模态叠加方法与模态组合方式，尝试理解响应谱分析背后的基本原理。

1、基本计算思想

正如前文所述，评估结构在给定时间历程载荷作用下的最大响应值，最直接的方法是采用完全法的瞬态动力学分析来获取结构的峰值响应，但该方法计算耗时较长，对硬件资源要求较高。因此，响应谱分析作为一种有效替代方法被广泛应用。该方法的基本计算思想是将复杂结构系统（多自由度模型）在长时间载荷作用下的响应问题解耦到时间、空间两个维度分别求解：其一，对于多自由度模型，通过模态分析，将结构系统等效为多个不同频率的单自由度弹簧振子；其二，对于长时间载荷作用下的响应问题，通过仿真或试验记录下不同频率单自由度弹簧振子在所给长时间载荷作用下的最大响应。最终，多自由度模型在所给载荷下的最大响应，可以直接提取各阶模态对应的最大响应值并进行适当的数学组合，从而得到结构在该时间历程载荷下最大响应值的合理估计，如下图所示。

其中，载荷响应谱代表了不同频率单自由度弹簧振子模型在所给时间历程载荷作用下出现的最大响应值，往往根据前期的仿真、试验或者是已有的数据得到。因此，实际工程中的响应谱分析，其载荷响应谱是已知的，上图中的分析路径二可以简化为模态分析获取模态频率与振型->各阶模态最大响应计算->模态组合估计结构最大响应。

2、响应谱的定义

响应谱描述了线性单自由度系统在给定时间历程载荷作用下（如地震加速度时程曲线），单自由度系统最大响应（如位移、速度、加速度或力）与其固有频率的关系曲线，其横坐标为系统的固有频率，纵坐标为系统的最大响应值。

尽管实际工程中的载荷响应谱都是给定的，但需要对载荷响应谱有进一步的理解，响应谱计算过程如下所示：

• 建立不同频率的单自由度弹簧振子系统，包括弹簧刚度，阻尼，振子质量；

• 跟踪记录每个单自由度弹簧振子在给定时间历程载荷下的最大响应（如位移、速度、加速度）；

• 以不同单自由度弹簧振子的固有频率为横坐标，该振子在给定时间历程载荷下的最大响应为纵坐标，绘制响应谱曲线。
  

需要指出的是，阻尼是蕴含于载荷响应谱中的一个关键参数。通过选取不同的阻尼比，即可按照上述方法绘制出相应的响应谱曲线。载荷响应谱本质上是分析对象（即单自由度弹簧振子）的固有频率与阻尼比的函数。它描述的并非所给时间历程载荷自身的频率成分（与傅里叶变换所得的频谱曲线不相同，但是相关），而是在给定时间历程载荷激励下，一系列具有不同频率和阻尼比的单自由度系统所能产生的最大响应。按照模态叠加的思想，结构的各阶模态可以视为相互独立的等效单自由度弹簧振子系统。因此，通过构造响应谱，就可以预估出结构的每一阶模态在该时间历程载荷下可能出现的最大响应，这正是响应谱分析的基本思想。

值得说明的是，不同类型的响应谱之间（如加速度、速度、位移响应谱）可以方便地进行转换，具体计算如下：

其中$S_a,S_v,S_d$分别代表加速度、速度、位移响应谱，$f$是相应的振子频率，换算时不同位置参与计算的$f$的是不一致的，背后的基本原理是由于各单自由度弹簧振子只会以本身的固有频率做振动，因此其最大加速度、速度、位移之间可以通过三角函数求导的关系给出。

3、模态叠加方法

已知一个多自由度系统承受基础加速度激励的动力学方程可以表示为：

等式右侧 $-\mathbf{M}\gamma {\ddot{x}}_g$ 表示由基础加速度引起的激励，其中$\gamma$为单位列向量（所激励自由度处为1，其余为0），${\ddot{x}}_g$为随时间变化的标量，表示加速度不同时刻的具体值。按照模态解耦的思路，将其转换为模态坐标方程：

按照前面响应谱的定义，一个单自由度系统在所给时间历程载荷下（这里即基础加速度激励${\ddot{x}}_g$）的最大响应可以通过响应谱曲线查询得到。上述模态坐标方程可以看作一个等效单自由度系统在${\Gamma }_i$倍基础加速度激励下的动力学响应，按照数学上变量代换的方法（还涉及绝对加速度和相对加速度的转换，这里不展开），容易得到模态响应$q_i$在该时间历程载荷下的最大加速度、最大位移分别为：

其中，$S_a({\omega }_i,{\zeta }_i)$ 表示给定加速度载荷响应谱在圆频率 ${\omega }_i$ 与阻尼比 ${\zeta }_i$ 处所对应的最大加速度响应值。由此可见，在载荷响应谱已知的前提下，求解各阶模态的最大响应无需进行复杂的时程积分运算，其过程简化为一次查询与换算：即利用模态分析获得的模态参数（${\omega }_i$, ${\zeta }_i$, ${\Gamma }_i$），直接从响应谱中查得对应的谱值 $S_a$，进而通过公式 $q_{i,max}={\Gamma }_i\cdot S_a/{\omega }_i^2$ 换算得到模态坐标的最大值。其中，模态参与因子 ${\Gamma }_i$ 可被视为一个“载荷分配系数”或“缩放因子”，它量化了外部载荷在各阶模态上的分配；而响应谱 $S_a$ 则定义了一个单自由度系统在真实载荷下的最大响应幅值。因此，二者之积 ${\Gamma }_i\cdot S_a$ 描述了该阶模态坐标在真实载荷下可能出现的最大响应水平。

4、模态组合方式

按照模态叠加法的思想，结构物理响应是各阶模态响应的瞬时叠加，即：

然而，响应谱计算得到的是每阶模态坐标在所给时间历程载荷下的最大值，所有模态最大值不太可能同时出现且符号相同，简单地将各阶模态的最大响应相加会使得最终估计过于保守。因此，所谓的模态组合方法，就是用来解决“如何将可能在不同时刻发生的模态最大值合理地组合成最大响应估计”的数学方法，常见的方法包括：平方和开方法（SRSS）、完全二次组合法（CQC），具体计算流程如下：

• 根据各阶模态最大值计算该阶模态下的物理响应最大值（以位移为例，速度、加速度、应力和应变等同理），即
  
• 根据所选的模态组合方式，对最终物理响应最大值进行估计，即：
  
  不同方法之间的差异主要是如何考虑模态间的耦合效应，这里不再赘述。此外，由于响应谱分析基于模态叠加的思想，因此也会涉及到模态截断，实际分析时需要关注模态有效质量是否满足要求。

5、总结

• 本文介绍了响应谱分析的基本原理，其基本思想在于将复杂结构系统（多自由度模型）在时间历程载荷作用下的响应问题，解耦为时间与空间两个维度分别求解，从而高效估计结构在给定载荷下的最大响应值，而不关注该响应发生的具体时刻；
• 阐述了响应谱的基本定义，响应谱所描述的并非载荷自身的频率成分（要区别于对给定载荷作傅里叶变换所得的频谱，但两者具有一定关联），而是在给定载荷激励下，一系列不同频率与阻尼比的单自由度系统所产生的最大响应构成的曲线；
• 介绍了响应谱分析中模态叠加与模态组合的基本原理，将各阶模态的最大响应求解过程，从复杂的时程积分运算简化为对响应谱的查询与换算，显著提升了计算效率，这也体现了响应谱分析的基本计算思想。

最后

以上内容是本人结合相关理论，对响应谱分析方法基本思想，响应谱、模态叠加、模态组合等相关概念的一些理解。在撰写过程中参考了诸多相关资料，由于资料来源较多，这里就不列举了，在此一并进行感谢！由于个人水平有限，文章内容难免存在不足之处，如有错误或疏漏，欢迎大家指正。

-完-
公众号已同步更新：仿真优化工匠
欢迎联系！