土体在循环荷载作用下的应力—应变响应具有典型的非线性、滞后性与路径依赖性，其本构表征构成路基抗震分析的理论基石。地震荷载诱发的土体动力行为以剪切模量退化和阻尼比演化为核心特征，其内在机制源于土颗粒骨架结构的循环重构与孔隙水压力的迁移耗散。对于饱和土体，不排水条件下有效应力路径的偏移将触发土体软化效应；非饱和土则因基质吸力的存在呈现更为复杂的强度变化规律^[6]。当前主流的黏弹塑性本构模型（如Davidenkov模型、边界面模型）通过引入内蕴变量刻画累积塑性变形，可有效描述土体循环硬化/软化现象。工程实践中，土体动力参数的确定需充分考虑原位应力历史、土质变异性和加载频率效应，其中小应变条件下剪切模量取值的精确性直接影响地震波传播模拟的可靠性。

地震荷载与地质环境耦合作用诱发的路基破坏呈现显著的模态分异，其损坏机制与表现形式详见下表1所示。

表1 公路路基典型震害模式及形成机制Table 1 Typical earthquake damage modes and formation mechanisms of highway subgrade

复杂地质背景下，多种震害模式往往呈现链式传导效应，软土地基的震陷变形将诱发邻近堤体附加弯矩，从而降低液化敏感区的临界加速度阈值；近断层区域的冲击型地震波则易与高填方路堤的自振模态形成共振，加速边坡浅层剥落向深层滑移转化。另外，地震动的频谱特性与持时参数通过调制能量输入方式，显著影响着震害模式的显性表现，长周期地震波对深厚软基的渗透作用加剧蠕变损伤；短时强脉冲波则主要激发填料的惯性响应，促使加筋结构发生界面脱黏失效^[7-8]。

软土地基通常指具有高含水率、高压缩性、低强度和低渗透性的土层，主要由淤泥、淤泥质土、冲填土等构成。在工程界定上，通常依据土体的物理力学指标进行划分。在地震工程领域，剪切波速V_S）是评价场地土类别和预测地震动放大效应的核心指标。根据建筑抗震设计规范，平均剪切波速V _S3（0 地表下30 m深度范围内的等效剪切波速）小于180 m/s的场地通常被划分为软弱土场地（如美国规范中的E类场地）。软土地基在地震中的主要危害在于其显著的“场地放大效应”，即软弱土层会对特定频段（通常是中长周期）的地震波产生共振放大，同时滤除高频成分，导致地表地震动峰值和频谱特性发生显著改变，从而极大地增加了上部路基结构的地震输入^[9]。

土体液化是指饱和的、松散的无黏性土（主要是砂土和粉砂）在地震等循环荷载作用下，孔隙水压力急剧升高，导致有效应力趋近于零，土体暂时丧失抗剪强度而呈现出类似液体状态的现象。液化判别通常分为两个层次：易发性评估和触发性分析。液化易发性（Susceptibility）：评估土层是否具备液化的内在条件^[10-11]。判别标准包括：地质年代与成因：通常认为全新世（Holocene）沉积的冲积、洪积砂土层具有较高的液化易发性；颗粒级配与塑性：颗粒级配曲线位于特定范围内的砂土、粉砂及含少量黏粒的砂土，以及满足特定塑性指数和含水率条件的低塑性粉土被认为是易液化的；物理状态：土层必须处于饱和状态，且相对密度较低。液化触发（Triggering）：在满足易发性条件的基础上，评估特定强度的地震动是否足以引发液化。目前应用最广的是基于应力的简化分析法，通过计算地震荷载比（CSR）和土体抗液化强度比（CRR）的比值，即抗液化安全系数（FSL=CRR/CSR）来进行判别^[12]。

近断层区域通常指距离活动断层地表迹线15至20 km范围内的区域。该区域的地震动与远离震源的“远场”地震动相比，具有显著不同的特征，这些特征对工程结构的破坏性尤为突出。方向性脉冲（Forward-Directivity Pulse）：当断层破裂向着场地传播时，由于多普勒效应，地震能量在短时间内集中释放，形成垂直于断层走向（Fault-Normal， FN）方向的、具有高幅值和长周期的速度脉冲；滑冲效应（Fling-Step）：由断层两侧的永久性构造位移引起的，表现为平行于断层滑动方向的单向速度脉冲和最终的永久地面位移；上盘效应（Hanging Wall Effect）：对于倾滑断层（如逆冲断层），位于断层上盘的场地所经历的地震动强度通常显著高于位于下盘的场地^[13]。

高填方路堤在工程界尚无统一的量化定义，其界定往往与工程背景和所关注的问题相关。在本研究的抗震分析语境下，高填方路堤被定义为填筑高度足以使其自身质量产生的地震惯性力、动力特性（如自振周期）以及对地基的附加应力成为控制其抗震稳定性的主导因素的路堤结构。通常，填土高度超过10~15 m的路堤可被视为高填方路堤^[14-15]。

模型构建基于非线性动力方程组，核心控制方程表达如下：

动力平衡方程：

式（1）中， M为质量矩阵， C=αM + βK为瑞利阻尼矩阵， α=2ω_i ω_j ζ/( ω_i + ω_j )， β=2ζ/( ω_i + ω_j )；{u}为节点位移向量；{F_地震}为地震惯性力矢量。

土体孔压控制方程（饱和区）：

式（2）中，k为渗透系数，γ_w为水容重，u_e为超静孔压，ε_v为体积应变，n为孔隙率，K_f为孔隙流体体积模量。

断层位错传递模型：

式（3）中，ΔD₀为基岩错动量，z为地表深度，r为至断层距离，r₀为特征衰减半径，λ、p为材料衰减参数。

模型采用分层映射技术整合四种地质工况：软土层赋予修正剑桥粘塑性特性，液化夹层嵌入孔压比状态变量，断层区设置非对称衰减场，高填方堤体配置正交异性硬化模型。自由场边界施加速度约束条件：

地震动输入库通过多元分类体系构建，特征分布如下表2所示。

表2 地震动输入库特征分布矩阵Table 2 Feature distribution matrix of ground motion input library

地震动前处理实施三阶段标准化：

基线漂移校正：通过零加速度平台扩展消除速度时程的残余位移偏移；

目标谱匹配：采用复小波包分解重构技术调控幅值谱密度分布；

概率调幅：按Ⅰ~Ⅳ类场地设防烈度分组归一化峰值加速度。

为系统量化复杂地质路基的抗震性能，构建响应、损伤、失效联控指标体系。指标筛选基于参数敏感性分析及工程破坏模式统计，形成全域覆盖的量化判据矩阵如表3：

该指标体系通过多物理场参量的协同监测，实现从微观土体损伤（r_u）到宏观结构失稳（Δd）的全链条抗震能力诊断。

软土层显著放大震陷与侧向位移，设计采用强度、排水、隔震三重优化体系：①桩网复合地基强化：采用直径60 cm、间距2.5 m的预应力管桩（设计承载力≥800 kN），桩顶铺设双向土工格栅（极限抗拉强度≥100 kN/m）形成刚性桩承层，实测沉降率较传统方案降低76%；②垂向排水体系：结合塑料排水板（通水量≥40 cm³/s）与30 cm厚砾石垫层（渗透系数>10⁻²cm/s），加速超静孔压消散，确保震后孔压比r_u<0.5；③隔震减震屏障：路基底部嵌入高阻尼橡胶隔震垫层（厚度≥50 cm，剪切模量0.6 MPa），将PGA输入削减40%。典型案例中，该组合方案使深厚软土区（H=12 m）的路肩位移差Δd稳定在0.3%以内，满足运营安全要求。

液化夹层诱发的瞬时失稳需实施阻液化、抗隆起、抑滑移协同策略，具体应碎石桩加密阻液化：采用等边三角形布桩（桩径1.0 m，间距1.8 m），置换率≥25%，形成有效应力扩散网络，使夹层埋深3 m区域的r_u峰值降至0.72。同时在坡脚处设置钻孔灌注桩筏基础（桩径1.2 m，嵌入稳定层5 m），筏板厚1.5 m并配双向钢筋网（配筋率≥0.8%），抑制坡脚隆起量δ_v不超过5 cm。最后，排沿潜在滑裂面布设DCP排水型锚索（倾角25°，设计抗拔力300 kN），水平间距2.0 m，同步实现孔隙水导出和抗剪能力提升，液化工况下τ/τ_f稳定值≤0.75。

高填方路基在强震作用下的灾变防控需突破传统均质填筑思维，转向刚度渐变协同耗能体系。核心在于路基内部建立加筋强度随埋深动态递增的机制，通过自下而上梯度配置不同模量的合成材料，实现地震波传播路径上的阻抗适配性过渡，显著弱化坡顶加速度放大效应。协同植入轻质发泡混凝土夹层形成水平耗能带，该柔性介质通过压缩变形高效吸收动能，与刚性筋材形成“刚柔并济”的抗力体系。坡面形态采用多级阶梯状缓坡重构技术，结合平台衔接与复合植被固土，利用阶梯落差逐级消耗坡体势能。深部滑移面重点布设嵌入式抗剪锚固群，通过锚杆与土体的摩擦耦合作用抑制深层塑性区贯通。