ANSYS力学仿真：一个桥梁支座在地震载荷下的非线性响应分析

桥梁支座是连接上部结构和下部结构的关键部件，在正常工况下几乎不起眼——但地震来临时，它就是整座桥梁的第一道防线。支座能否在强震下正常工作，直接决定了桥梁是否会发生落梁或倒塌等灾难性破坏。ANSYS力学仿真在支座的抗震设计中提供了不可替代的定量评估能力。

项目背景：一个跨越活动断裂带的公路桥

团队处理过一个西部山区公路桥梁的支座抗震评估项目。桥梁跨度60 m，总长240 m，位于活动断裂带附近，抗震设防烈度为8度（0.3g峰值加速度）。客户使用了天然橡胶支座（LRB），需要确认在罕遇地震（大震，0.51g）下支座的剪应变不超过规范限值（250%），且支座不会发生失稳或破坏。

支座的几何参数为：直径600 mm，总高度160 mm（其中橡胶层总高120 mm，钢板层总高40 mm，由10层12 mm橡胶和11层钢板交替叠合）。这种叠层结构的设计目的是在保持较大竖向承载力的同时提供足够的水平柔度来适应地震位移。

橡胶超弹性本构模型的建立

橡胶材料的大变形力学行为不能用线弹性模型描述——其应力-应变关系呈现高度非线性，且几乎不可压缩（泊松比≈0.499）。ANSYS中使用超弹性本构模型来描述这种行为。

团队采用Mooney-Rivlin两参数模型，参数通过拟合客户提供的单轴拉伸、双轴拉伸和平面剪切实验数据确定。拟合结果为：C₁=0.42 MPa、C₁₀=0.11 MPa（Mooney-Rivlin常数）。

Mooney-Rivlin模型在中等应变（<100%）范围内的描述精度很好，但在大应变（>200%）下可能出现偏差。对于这个项目的分析目的（支座最大剪应变可能达到200-250%），Mooney-Rivlin模型的精度是够用的。如果需要更高精度，可以切换到Yeoh或Ogden模型——但参数拟合需要更多的实验数据。

地震波输入与瞬态分析

瞬态动力学分析中，地震载荷通过基座的加速度时程输入。团队选取了三条地震波作为输入：El Centro波（1940年）、Taft波（1952年）和一条人工拟合波（基于场地反应谱生成）。每条波按8度大震要求缩放到0.51g峰值加速度。

分析在ANSYS Transient Structural模块中进行，时间步长取0.001 s，总分析时长20 s。由于橡胶的大变形和接触非线性，每个时间步需要迭代收敛。实际计算中，收敛困难主要出现在支座发生大剪切变形的阶段（约3-8 s），团队通过增加最大迭代次数（从默认的26次增加到50次）和采用弧长法（Arc-length method）解决了收敛问题。

关键结果与评估

三条地震波下的支座最大剪应变分别为：El Centro 218%、Taft 195%、人工波232%。三条波的平均值为215%，低于250%的规范限值，但有接近限值的倾向。

支座的竖向位移在地震过程中出现了约3 mm的波动（支座的总高160 mm，波动比约1.9%），竖向应力的最大值为12.3 MPa——远低于橡胶的压缩强度极限（约30 MPa）。

非线性分析与线性分析的结果对比揭示了关键差异。线性分析（假设材料为线弹性）预测的最大剪应变仅为145%，远低于非线性分析的结果——线性分析完全忽略了橡胶在大变形下的刚度软化效应，严重低估了支座的变形响应。这个差异说明，对于橡胶支座的抗震分析，非线性建模是必须的——线性分析的结果可能给出错误的安全判断。

客户根据分析结果决定：保留当前LRB方案，但在支座两侧增设限位装置，防止极端情况下支座剪应变超过限值。ANSYS力学仿真在这个项目中提供了线性分析无法给出的关键信息——大变形非线性响应对支座安全性的影响远比直觉预期的大。

如需了解更多的科研计算案例与项目经验，可以参考 https://www.keyanxueshu.com/ 上的技术文章。