地震作用下的结构响应
地震是自然界中释放能量最剧烈的现象之一，其对工程结构的破坏力主要源于地面运动引起的惯性力。在经典力学框架下，地震作用并非直接施加在结构上的外力，而是通过地基运动传递给上部结构，进而激发结构的振动响应。理解这一过程的核心在于掌握“地震动”（Ground Motion）的特性，包括其峰值加速度（PGA）、频谱特性以及持续时间。
地震波通常包含体波（P波和S波）和面波（Love波和Rayleigh波）。对于大多数高层建筑和桥梁结构，S波产生的水平剪切效应是造成破坏的主要因素。在工程分析中，我们通常将地震输入简化为一维或二维的加速度时程曲线，作为动力学方程的激励项。
单自由度系统的动力学方程
为了简化复杂结构的分析，工程上常采用单自由度（SDOF）模型作为基础。一个典型的SDOF系统由质量 $m$、刚度 $k$ 和阻尼 $c$ 组成。根据牛顿第二定律，该系统在地震作用下的运动方程可表示为：
$$ m\ddot{u}(t) + c\dot{u}(t) + ku(t) = -m\ddot{u}_g(t) $$
其中：
- $u(t)$ 为结构相对于地面的位移；
- $\ddot{u}_g(t)$ 为地面加速度时程；
- $\ddot{u}(t)$ 和 $\dot{u}(t)$ 分别为结构的绝对加速度和速度。
该方程揭示了地震惯性力 $-m\ddot{u}_g(t)$ 与结构恢复力 $ku(t)$ 及阻尼力 $c\dot{u}(t)$ 之间的动态平衡关系。通过求解此微分方程，可以获得结构在地震作用下的位移、速度和加速度时程，从而评估结构的受力状态。
多自由度系统的模态分析
实际工程结构（如高层建筑、长跨度桥梁）具有多个自由度，直接求解高阶微分方程组计算量巨大。因此，工程实践中广泛采用模态分析法（Modal Analysis）。
- 特征值分解：首先求解结构的固有频率和振型。对于无阻尼系统，特征方程为 $[K]{\phi} = \omega^2[M]{\phi}$，其中 $[K]$ 为刚度矩阵，$[M]$ 为质量矩阵，$\omega$ 为固有频率，${\phi}$ 为振型向量。
- 模态叠加原理：利用振型的正交性，将多自由度系统解耦为若干个独立的单自由度方程。
- 响应合成：计算各阶模态在地震激励下的响应，并根据特定规则（如平方和开根号法 SRSS 或完全平方和开根号法 CQC）将各模态响应组合，得到结构的总响应。
这种方法极大地降低了计算复杂度，同时能够准确捕捉结构的主要动力特性。
结构响应参数与抗震设计准则
在评估结构安全性时，工程师关注的关键响应参数包括：
- 层间位移角：衡量结构变形能力的重要指标，过大的层间位移可能导致非结构构件（如填充墙、幕墙）破坏。
- 基底剪力：反映结构整体受到的水平地震作用大小，是基础设计的关键依据。
- 绝对加速度：影响结构内部构件及非结构设备的安全。
抗震设计遵循“小震不坏、中震可修、大震不倒”的三水准设防目标。这意味着结构在遭遇低于设防烈度的地震时应保持弹性工作状态；在设防烈度地震下允许出现可修复的塑性变形；在超越设防烈度的强震下，虽可能发生严重破坏，但不应发生倒塌。
数值模拟与工程应用示例
现代工程分析依赖于高性能计算软件（如 SAP2000, ETABS, ABAQUS 等）进行非线性时程分析。以下是一个简化的分析流程：
- 建立模型：输入几何尺寸、材料属性（混凝土强度、钢筋等级）及边界条件。
- 选取地震波：根据场地类别和设计地震分组，选取符合规范要求的天然地震记录或人工合成波。
- 时程积分：采用 Newmark-β 法或 Wilson-θ 法对运动方程进行数值积分，步长通常小于结构最短周期的一半，以确保计算精度。
- 结果后处理：提取关键截面的内力包络图和位移时程图，校核是否满足规范限值。
例如，在一栋 50 层高层办公楼的分析中，通过时程分析发现，在强震作用下，结构底部剪力达到 12,000 kN，最大层间位移角为 1/800，均满足《建筑抗震设计规范》的要求。这表明结构具备足够的延性和耗能能力，能够有效抵御地震灾害。
结语
地震作用下的结构响应分析是连接经典力学理论与工程实践的桥梁。通过精确的动力学建模、合理的模态分解以及严谨的数值模拟，工程师能够预测结构在地震中的行为，优化设计方案，从而保障人民生命财产安全。随着计算力学的发展，非线性纤维模型和接触力学的应用将进一步提高分析的精度，为复杂结构体系的抗震设计提供更可靠的依据。