路基冻胀对有砟轨道动力学的影响：基于混合离散元法（DEM）与车辆模型的研究

针对寒区铁路路基冻胀对轨道结构的影响，建立了轨道-路基离散元模型与车辆-轨道动态耦合模型，揭示冻胀导致碎石轨道互锁结构破坏，孔隙率增加及刚度分布不均，引发轮轨垂向力激增和车辆加速度异常，冻胀长度短于8米且幅度超30毫米时风险显著。

　　
本研究聚焦于寒冷地区铁路路基冻胀变形对轨道结构动态响应的影响机制，通过建立离散元（DEM）与车辆-轨道耦合模型，揭示了冻胀引发的结构性破坏过程及其对列车运行安全的威胁。研究团队针对现有理论在颗粒级配表征、冻胀-轨道耦合机制等方面的不足，创新性地构建了非迭代耦合的数值模拟体系，为冻区铁路轨道维护提供了理论支撑。

在研究方法层面，采用分阶段耦合策略实现多尺度分析：首先基于离散元方法构建路基-轨道复合体系，通过调整颗粒接触参数模拟冻胀导致的孔隙率变化与刚度分布不均；其次建立包含车体垂向运动、转向架横摆及轮轨垂向力的多自由度车辆-轨道动态模型，重点考察8米以下短波长冻胀变形对高速列车（时速350公里）的冲击效应。研究突破传统连续介质假设，通过施加体积膨胀约束直接生成冻胀层，实现了冻胀变形与轨道几何形态的物理映射，有效避免了理想化余弦函数表征的局限性。

研究发现冻胀变形通过双重路径影响轨道稳定性：微观层面，冻胀区孔隙率增加达17%-23%（对比非冻胀区），导致颗粒间接触应力分布呈现梯度变化，在冻-融交界带形成局部应力集中区；宏观层面，冻胀导致轨道支撑体系发生结构性改变，支撑链路从传统轨枕下方转移至冻胀前沿区域，造成约38%的轨枕在过渡带处于悬空状态。这种力学传递路径的突变，使得原本均匀分布的轮轨垂向力在冻胀区产生3-5倍于常规值的峰值，引发轨道几何状态失稳。

研究特别揭示了冻胀波长与轨道响应的强关联性：当冻胀波长≤6米时，其变形模式与现场观测的余弦型曲线高度吻合（R2=0.91）；而当波长>8米时，轨道变形呈现非对称的阶梯状分布特征。这种波长敏感性导致传统轨道检测设备（如轨道几何状态检测器）在识别冻胀型轨道缺陷时存在32%-45%的漏检率。在高速运营工况下，冻胀长度<8米的变形体（冻胀幅度达30毫米）会使轮轨垂向力激增至120-150kN区间，较正常状态提升60%-70%，同时导致车辆垂向加速度超过0.3g的临界安全阈值。

研究突破传统冻胀力学模型的三大局限：1）首次实现冻胀区颗粒级配与力学响应的实时反馈调节；2）建立非迭代耦合机制，将DEM计算耗时从常规的12-18小时压缩至4-6小时；3）开发基于深度学习的轨道几何状态预警算法，将冻胀变形的预测精度提升至92.3%。这些技术突破使冻区铁路轨道维护从"事后抢修"向"智能预警"模式转变，为寒区高速铁路的可持续发展提供了关键技术支撑。

在工程应用方面，研究团队针对东北地区的典型冻胀类型（年冻深1.2-1.8米，冻胀幅度20-50毫米），提出了分阶维护策略：对于波长≤6米的稳定冻胀带，建议采用周期性纵向压实（频率≤2次/年）；对6-8米过渡带实施动态监测（采样频率≥0.5Hz）；而对于>8米的动态冻胀区，推荐采用可调节轨枕结构（调节行程≥50mm）。经模拟验证，该策略可使轨道几何状态偏差控制在±2mm范围内，将冻胀相关病害发生率降低至12%以下。

研究还揭示了冻胀-车辆-轨道系统的非线性耦合特性：当冻胀幅度超过25毫米时，系统将进入失稳临界状态，此时轨道横向刚度衰减达40%，导致轮轨横向力增加2.8倍。这种非线性关系在传统预测模型中难以准确捕捉，因此研究团队开发了基于粒子群优化的参数辨识算法，可同时反演冻胀区的水热耦合参数（温度梯度误差≤3%，含水量预测偏差<5%）。

在技术验证方面，研究团队选取长春至白城铁路的实测数据（2023-2025监测周期），通过双模型耦合验证发现：冻胀区颗粒破碎率较非冻胀区高出1.8倍，且破碎颗粒多集中在冻-融界面（深度0.3-0.5米）。这解释了为何传统设计将冻胀控制重点放在表层（0-0.5米），而本研究通过DEM微观模拟发现，冻胀变形的67%能量源自底层（0.5-1.0米）的水分迁移，这为优化路基冻胀防治工程（如增设深层排水管）提供了新依据。

该研究成果已应用于哈大高铁寒区段轨道健康监测系统，部署了具有抗冻胀干扰能力的FBG传感器阵列（采样点密度达2点/米），结合提出的非迭代耦合算法，实现了冻胀变形的实时预警（预警响应时间<15分钟）。工程实践数据显示，采用该监测预警系统的线路，2024年度冬季冻胀导致的轨道几何状态异常次数减少82%，维护成本降低37%，验证了理论模型的工程适用性。

研究同时指出当前存在的三大技术瓶颈：1）颗粒级配对冻胀响应的量化模型尚未完善；2）动态耦合计算中的收敛性问题仍需解决；3）极端冻融循环（>500次）下颗粒接触网络的稳定性有待验证。这些发现为后续研究指明了方向，特别是建议开发基于机器学习的颗粒接触网络演化预测模型，以提升冻胀变形的长周期预测精度。

在学术贡献层面，本研究建立了首个完整的冻胀-轨道-车辆多物理场耦合分析框架，该框架已形成标准化计算流程（包含5个核心模块、18个输入参数、9类输出指标），可拓展应用于其他冻土工程场景。研究团队开发的DEM-车辆耦合计算平台（命名"冻轨通"）在算力需求方面取得突破，通过并行计算技术将万颗粒系统的模拟速度提升至传统方法的8.3倍，为大规模工程应用奠定基础。

最后需要强调的是，本研究未涉及冻胀预测模型开发，而是专注于冻胀变形对既有轨道结构的动态响应影响。这种聚焦策略有效规避了参数敏感性带来的模型泛化难题，使研究成果更贴近工程实际需求。未来研究可结合数字孪生技术，构建包含冻胀变形预测、轨道状态感知、维护决策支持的智能闭环系统，这将是寒区铁路智能化运维的关键突破方向。