离散-连续方法下昔格达地层铁路隧道基底动力特性与长期沉降的深度剖析

离散-连续方法下昔格达地层铁路隧道基底动力特性与长期沉降的深度剖析一、绪论1.1研究背景与意义随着我国交通基础设施建设的持续推进，铁路隧道作为重要的交通节点，在复杂地质条件下的建设面临着诸多挑战。昔格达地层作为一种特殊的地质构造，广泛分布于我国西南地区，其独特的物理力学性质给铁路隧道的建设和长期运营带来了一系列问题。昔格达地层主要由细砂岩、粉砂岩、泥质粉砂岩和粉砂质泥岩等形成的一套静水河湖相沉积地层，具有似土非土、似岩非岩的半成岩特征。其强度较低，在水的作用下承载力和抗剪强度会大幅度降低，易产生滑坡、塌方等地质灾害。在昔格达地层中修建铁路隧道，隧道基底容易出现变形、沉降等问题，严重影响隧道的结构安全和列车的运行平稳性。以成昆铁路为例，部分路段穿越昔格达地层，在运营过程中，隧道基底出现了不同程度的沉降和变形，导致轨道不平顺，增加了列车运行的安全风险，同时也加大了隧道的维护成本。据统计，成昆铁路部分昔格达地层隧道的维修频率比其他地层隧道高出30%-50%，每年用于隧道维护的费用高达数千万元。此外，由于昔格达地层的特殊性，隧道施工过程中也容易出现坍塌、涌水等事故，不仅延误工期，还造成了巨大的经济损失。因此，深入研究昔格达地层铁路隧道基底的动力特性及长期沉降规律，对于保障铁路隧道的安全建设和长期稳定运营具有重要的现实意义。传统的研究方法在处理昔格达地层这种复杂地质问题时存在一定的局限性。而离散-连续方法作为一种新兴的数值分析方法，能够综合考虑岩体的离散性和连续性，更加真实地模拟昔格达地层在隧道施工和运营过程中的力学行为。通过离散-连续方法，可以对隧道基底的应力分布、变形机制以及长期沉降过程进行精细化分析，为隧道的设计、施工和维护提供更加科学的依据。将离散-连续方法应用于昔格达地层铁路隧道基底的研究，有助于揭示其动力特性和长期沉降的内在规律，从而优化隧道的设计方案，提高施工质量，降低工程风险，具有重要的理论价值和工程应用前景。1.2国内外研究现状1.2.1昔格达地层动力特性研究昔格达地层的动力特性研究一直是岩土工程领域的重要课题。早期研究主要集中在对其物理力学性质的基本测试上。学者张家明等通过大量动三轴试验，研究了不同固结条件下西昌昔格达组黏土岩的动应力-动应变关系、动弹性模量、阻尼比及动强度特性，发现其在一定条件下的动应力-动应变关系符合幂函数模型，动弹性模量随固结压力或固结主压力比增大而增大，阻尼比随动应变增大而增大，但围压或固结主压力比对阻尼比影响较小。连煜皓等利用振动三轴仪对昔格达土进行抗液化强度动力试验，研究了其动强度、动孔隙水压力随振动周次的变化，以及固结压力和固结比对动强度的影响，并得到了破坏标准为特定应变值和不同振动周次所对应的动强度指标。随着研究的深入，数值模拟方法逐渐应用于昔格达地层动力特性研究。有学者借助有限元数值模拟软件，结合固体介质中应力波的传播理论，对西昌昔格达地层的场地动力特性进行探讨，分析了场地加速度响应的空间分布特征及其约束机制，指出地震波传播到两种介质交界面处会发生反射叠加导致地表地震动放大，且两种介质泊松比差异越大，放大效应越显著。1.2.2列车荷载下隧道基底沉降研究在列车荷载作用下，隧道基底沉降是影响隧道结构安全和列车运行平稳性的关键因素。国外学者较早开展了相关研究，如日本学者通过对新干线隧道的长期监测，分析了列车振动荷载下隧道基底的变形规律，提出了基于实测数据的沉降预测经验公式，但该公式的通用性受到地质条件和隧道结构形式等因素的限制。国内在这方面的研究也取得了丰硕成果。一些学者运用有限元方法，建立了考虑列车-轨道-隧道-地基相互作用的耦合模型，研究列车荷载作用下隧道基底的动力响应和沉降特性。张强等基于有限元分析，对隧道结构变形进行计算，分析了不同工况下隧道基底的应力和位移分布。但此类研究大多将地层视为连续介质，对于像昔格达地层这种具有明显离散性和复杂力学特性的地层，模拟结果与实际情况存在一定偏差。还有学者通过现场监测，对运营中的铁路隧道基底沉降进行长期观测，分析了沉降随时间的变化规律以及影响因素，但现场监测往往受到诸多条件限制，难以全面揭示沉降的内在机制。1.2.3离散-连续耦合方法在隧道工程中的应用研究离散-连续耦合方法作为一种新兴的数值分析方法，近年来在隧道工程中得到了越来越广泛的关注和应用。在岩土体的离散元分析方面，PFC（ParticleFlowCode）等软件被广泛应用，能够较好地模拟岩土体颗粒间的相互作用和大变形行为。而连续介质的有限元方法则在处理连续介质的力学问题上具有优势。将两者耦合，可以综合考虑岩体的离散性和连续性，更加真实地模拟隧道工程中的力学行为。中铁十五局在隧道施工领域的技术创新中，针对砂卵石地层的挑战，构建了开挖面稳定性控制理论体系，实现了离散-连续耦合模拟，提高了开挖面稳定性评估效率。赵东平等开发了离散-连续耦合数值模型，实现了泥岩卵石复合地层盾构掘进动力学计算，为盾构施工控制参数的定量化评价提供了新方法。吴祖松等采用离散-连续耦合分析方法，对土石混合料高填方边坡破坏机理进行研究，取得了有价值的成果。然而，目前离散-连续耦合方法在隧道工程中的应用仍存在一些问题。一方面，耦合算法的精度和效率有待进一步提高，在处理大规模计算时，计算时间较长，内存需求较大。另一方面，如何准确地确定离散元和连续元之间的边界条件和相互作用关系，仍是需要深入研究的课题。同时，将离散-连续耦合方法应用于昔格达地层铁路隧道基底研究的文献相对较少，针对昔格达地层特殊的物理力学性质和隧道基底复杂的受力条件，如何优化耦合模型，提高模拟的准确性，还需要进一步探索。1.3研究内容与方法1.3.1研究内容（1）昔格达地层物理力学性质试验研究。通过现场勘探和室内试验，获取昔格达地层的基本物理参数，如密度、含水率、孔隙比等，以及力学参数，包括抗压强度、抗剪强度、弹性模量、泊松比等。研究不同含水率、不同固结状态下昔格达地层的力学特性变化规律，分析其在复杂地质条件下的力学响应机制。例如，针对昔格达地层遇水易软化、崩解的特性，重点研究含水率对其强度和变形特性的影响，为后续的数值模拟和理论分析提供准确的参数依据。（2）离散-连续耦合模型的建立与验证。基于离散元理论和有限元理论，建立适用于昔格达地层铁路隧道基底分析的离散-连续耦合模型。在离散元部分，采用颗粒流方法模拟昔格达地层颗粒间的相互作用和大变形行为，考虑颗粒的接触力、摩擦力、粘结力等因素；在连续元部分，运用有限元方法模拟隧道结构和周围连续介质的力学行为。通过与现场监测数据和室内试验结果对比，验证耦合模型的准确性和可靠性。对模型进行参数敏感性分析，研究不同参数对模拟结果的影响，优化模型参数，提高模拟精度。（3）列车荷载作用下隧道基底动力响应分析。将列车荷载简化为移动的动力荷载，施加到建立的离散-连续耦合模型上，模拟列车运行过程中隧道基底的动力响应。分析隧道基底在列车荷载作用下的应力、应变分布规律，研究不同列车速度、轴重、轨道不平顺等因素对隧道基底动力响应的影响。例如，通过改变列车速度，观察隧道基底应力和应变的变化情况，探讨列车速度与隧道基底动力响应之间的关系，为隧道的动力设计提供理论支持。（4）隧道基底长期沉降预测模型的建立与分析。考虑昔格达地层的蠕变特性、地下水渗流、列车长期振动等因素，建立隧道基底长期沉降预测模型。采用经验公式法、数值模拟法或两者结合的方法，对隧道基底的长期沉降进行预测。通过对不同预测模型的结果进行对比分析，评估模型的可靠性和适用性。结合实际工程案例，对预测结果进行验证和修正，提高长期沉降预测的准确性，为隧道的长期运营维护提供决策依据。（5）隧道基底沉降控制措施研究。根据隧道基底动力特性和长期沉降分析结果，提出有效的沉降控制措施。从隧道设计、施工工艺、地基处理等方面入手，探讨不同控制措施对隧道基底沉降的影响。例如，在隧道设计阶段，优化隧道结构形式和尺寸，增加基底的承载能力；在施工工艺方面，采用合理的开挖方法和支护措施，减少施工对基底的扰动；在地基处理方面，采用注浆加固、强夯等方法，提高昔格达地层的力学性能。通过数值模拟和工程实例分析，评估各种沉降控制措施的效果，为实际工程提供可行的技术方案。1.3.2研究方法（1）文献研究法。广泛查阅国内外关于昔格达地层特性、隧道工程力学、离散-连续耦合方法、隧道基底沉降等方面的文献资料，了解该领域的研究现状和发展趋势，分析现有研究的不足之处，为本研究提供理论基础和研究思路。通过对相关文献的综合分析，总结昔格达地层的物理力学性质、隧道基底在列车荷载作用下的动力响应规律以及沉降预测方法等方面的研究成果，明确本研究的重点和难点问题。（2）试验研究法。开展现场勘探和室内物理力学试验，获取昔格达地层的物理力学参数。在现场勘探中，采用地质钻探、地球物理勘探等方法，了解昔格达地层的分布特征、地质构造等信息。在室内试验中，进行常规土工试验，如密度试验、含水率试验、颗粒分析试验等，以及力学试验，包括三轴压缩试验、直剪试验、动三轴试验等，研究昔格达地层在不同条件下的力学性能。通过试验研究，为数值模拟和理论分析提供可靠的数据支持，同时验证理论分析和数值模拟结果的准确性。（3）数值模拟法。利用离散元软件（如PFC）和有限元软件（如ANSYS、ABAQUS），建立离散-连续耦合模型，对昔格达地层铁路隧道基底在列车荷载作用下的动力特性和长期沉降进行数值模拟分析。在数值模拟过程中，合理设置模型参数，模拟实际工程中的各种工况，如不同的列车速度、轴重、隧道埋深、地层条件等。通过数值模拟，可以直观地观察隧道基底的应力、应变分布情况以及沉降发展过程，深入研究各种因素对隧道基底力学行为的影响机制，为工程设计和施工提供科学依据。（4）理论分析法。基于岩土力学、隧道力学等相关理论，建立隧道基底在列车荷载作用下的力学分析模型，推导相关计算公式，分析隧道基底的动力响应和沉降变形机理。运用弹性力学、塑性力学理论，研究昔格达地层在隧道开挖和列车荷载作用下的应力应变状态；采用蠕变理论、渗流理论等，分析隧道基底长期沉降的影响因素和发展规律。通过理论分析，为数值模拟和试验研究提供理论指导，揭示隧道基底力学行为的本质规律。（5）工程案例分析法。结合实际的昔格达地层铁路隧道工程案例，对本文提出的理论和方法进行应用验证。收集工程现场的监测数据，包括隧道基底的沉降、应力、应变等数据，与数值模拟和理论分析结果进行对比分析，评估本文研究成果的可靠性和实用性。通过工程案例分析，总结经验教训，进一步完善研究成果，为类似工程提供参考和借鉴。二、昔格达地层特性与隧道工程概况2.1昔格达地层特征昔格达地层是在特定地质历史时期和环境条件下形成的。它主要分布于中国青藏高原东南缘金沙江、雅砻江、安宁河、大渡河及其支流河谷中，长期以来与华北泥河湾组、三门组，云南元谋组共称为中国早更新世代表地层。其形成与青藏高原隆起这一重大地质事件密切相关，在青藏高原隆起过程中，该区域的地形地貌、气候环境等发生了显著变化，为昔格达地层的形成创造了条件。从岩石组成来看，昔格达地层主要由细砂岩、粉砂岩、泥质粉砂岩和粉砂质泥岩成韵律互层组成，是一套静水河湖相沉积地层。在盐边县红格镇昔格达村附近的典型地层中，厚度约280米，上段整体为黄色、土黄色、土白色厚层含土中细砂与粉砂、粉砂质黏土、黏土互层，钙质胶结；下段为灰色、暗灰色、青灰色、砖灰色的含土中细砂与黏土、亚黏土、粉砂互层，含钙质结核。这种岩石组成使得昔格达地层具有独特的结构构造。沉积旋回中粉砂质黏土纹泥特征显著，微层理极其发育，易成片理化，以水平层理为主，偶见交错层理、波状层理等，局部受同沉积构造影响，可见小型断层、褶皱。在安宁河谷其他地区沉积物中还发育介形虫、硅藻等，局部规模和富集较好，形成具有工业价值的硅藻矿。昔格达地层具有诸多特殊的工程性质。它属于极软岩，强度低，先期固结压力较小，没有完成沉积成岩作用，是一套弱胶结半成岩地层，具有似土非土、似岩非岩的物理力学特征。其水稳性差，浸水后易崩解，强度降低幅度较大。相关试验表明，昔格达地层天然状态下含水率通常为17%-24.6%，饱和度为85%-100%，其抗压强度为1.07-1.73MPa，烘干状态下抗压强度为13.2-13.4MPa，饱和状态下的抗压强度仅为0.79-0.89MPa，充分反映了不同含水率状态下黏土岩的抗压强度差值较大，在饱和状态下强度低的特性。攀枝花学院的试验通过配制相同干密度1.58g/cm³而不同含水率的试样进行不固结不排水的快剪试验，发现从昔格达土体天然风干含水率5%到抽气饱和的含水率29.5%，土体的粘聚力随着含水率的增加，由5%含水率时的62Kpa骤然降低至饱和含水率(29.5%)的0.6Kpa，表明含水率对其粘聚力影响极大。此外，昔格达地层透水性弱，泥页岩常起隔水作用，在表水下渗、排水不畅的情况下，易引起昔格达之上的第四系覆盖层或堆填物沿昔格达顶部滑动，给工程建设带来极大的安全隐患。2.2铁路隧道工程实例以成昆铁路复线峨眉至米易段的范家湾隧道为例，该隧道具有重要的研究价值。范家湾隧道进口里程为D2K528+495，出口里程为D2K528+760，全长265m，线路纵坡为-11.5%的单面下坡。隧道为单洞双线，全部处于V级围岩中，洞身主要通过砂岩夹泥岩、页岩，最大埋深约20m。该隧道未见基岩出露，下伏昔格达组地层成岩作用差，属极软岩，粉质粘土为弱膨胀土，昔格达粘土岩为弱膨胀岩，隧道洞身浅埋段长，工程地质条件极差。在隧道施工过程中，与基底相关的问题频发。由于昔格达地层遇水易泥化崩解，围岩强度明显降低，基底承载力不足的问题凸显。在一次强降雨后，隧道基底出现了明显的沉降，部分地段沉降量达到了50mm，导致已施工的仰拱和道床出现开裂现象，严重影响了隧道的施工进度和结构安全。在开挖过程中，掌子面岩体稳定性差，多次出现小规模坍塌，这与昔格达地层的特性密切相关，其强度低、胶结弱，难以在开挖过程中维持自身稳定，对基底的稳定性也产生了连锁反应。隧道的排水问题也给基底带来了困扰。昔格达地层透水性弱，在隧道排水过程中，长期排水可能会带走昔格达土壤中的细颗粒，导致隧道结构和围岩之间产生空洞，地表出现塌陷，进一步威胁基底的稳定性，同时排水系统也容易因土壤中细颗粒的流失而阻塞、失效。这些问题的出现，充分说明了研究昔格达地层铁路隧道基底动力特性及长期沉降的紧迫性和必要性，也为后续基于离散-连续方法的研究提供了现实依据。三、离散-连续方法原理与应用3.1离散-连续方法基本原理离散元与连续元是两种在岩土工程数值分析中广泛应用的方法，它们各自基于不同的理论基础，在解决岩土工程问题时具有独特的优势和适用范围。离散元方法最初由P.A.Cundall于1971年提出，其基本理论是将不连续介质离散为刚性单元，各个单元的运动满足牛顿第二定律。在离散元模型中，单元间的接触、碰撞行为采用适当的接触力学模型描述，通过追踪每个离散单元的运动和相互作用，来反映整个系统的力学行为。以颗粒离散元为例，对于第i个颗粒，其平动方程为m_{i}\frac{d^{2}x_{i}}{dt^{2}}=m_{i}g+\sum_{j=1}^{c_{i}}(F_{ij}+C_{n}v_{ij})，其中m_{i}和x_{i}分别表示颗粒i的质量和位置向量，t为时间，g为重力加速度，c_{i}表示颗粒i的接触数，F_{ij}为颗粒j对颗粒i的接触力，C_{n}为法向粘性接触阻尼系数，v_{ij}为颗粒i与j在接触点处的相对速度。转动方程为I_{i}\frac{d\omega_{i}}{dt}=\sum_{j=1}^{c_{i}}r_{ij}\timesF_{ij}，其中I_{i}和\omega_{i}分别为颗粒i的转动惯量与角速度，r_{ij}为从颗粒i质心指向接触点的矢量。离散元法能够很好地模拟岩土体颗粒间的相互作用、大变形以及破坏过程，适用于处理节理发育、破碎岩体等具有明显离散性的地质问题。连续元方法则基于连续介质力学理论，将研究对象视为连续的介质，通过求解偏微分方程来描述其力学行为。有限元法是连续元方法中应用最为广泛的一种，它利用分段的技术，将复杂的结构分割成很多小的单元进行计算，这些小的单元的计算结果可以被整合成整个结构的分析结果。在有限元分析中，通过建立单元的刚度矩阵，组装成整体刚度矩阵，结合边界条件和载荷条件，求解线性方程组，得到结构的位移、应力和应变等力学响应。连续元方法在处理连续介质的力学问题上具有优势，能够精确地计算应力和应变的分布，适用于模拟均质、连续的岩土体以及结构物的力学行为。离散-连续耦合方法的原理是将离散元与连续元相结合，充分发挥两者的优势，以解决复杂的地质问题。在离散-连续耦合模型中，通常将具有明显离散特征的部分（如节理、裂隙、颗粒等）用离散元方法模拟，而将连续介质部分（如完整的岩体、土体等）用连续元方法模拟。两者之间通过一定的耦合算法实现信息的传递和交互，从而实现对整个系统的模拟。在隧道工程中，对于隧道周围的岩体，如果岩体中存在大量的节理和裂隙，这些节理和裂隙对岩体的力学行为有重要影响，就可以将节理和裂隙用离散元模拟，而将完整的岩体用连续元模拟。离散元部分可以准确地模拟节理和裂隙的张开、闭合、滑动等行为，以及颗粒间的相互作用，而连续元部分则可以精确地计算岩体的应力和应变分布。通过耦合算法，离散元部分和连续元部分可以相互传递力和位移信息，使得整个模型能够更加真实地反映岩体在隧道开挖和列车荷载作用下的力学响应。离散-连续耦合方法在解决复杂地质问题时具有显著的优势。它能够综合考虑岩体的离散性和连续性，更加真实地模拟岩土体的力学行为，避免了单一方法在处理复杂地质条件时的局限性。对于昔格达地层这种具有似土非土、似岩非岩的半成岩特征，既存在颗粒间的离散特性，又具有一定的连续介质特性的地质体，离散-连续耦合方法能够更准确地模拟其在隧道施工和运营过程中的力学响应。离散-连续耦合方法还可以模拟岩土体在大变形、破坏等复杂情况下的行为，为工程设计和施工提供更全面、准确的依据。3.2在隧道工程中的应用及实现离散-连续方法在隧道工程数值模拟中具有广泛的应用，能够有效模拟隧道开挖、支护等复杂过程。在隧道开挖模拟方面，离散-连续方法可以真实地反映岩体在开挖过程中的力学响应。在开挖过程中，岩体中的节理和裂隙会发生张开、闭合和滑动等行为，离散元部分能够精确地模拟这些离散特性，而连续元部分则可以计算完整岩体的应力和应变分布。通过耦合算法，两者相互作用，使得模拟结果更加符合实际情况。这种模拟可以帮助工程师预测隧道开挖过程中可能出现的坍塌、变形等问题，从而提前制定相应的应对措施。在隧道支护过程模拟中，离散-连续方法同样发挥着重要作用。对于喷射混凝土、锚杆等支护结构，连续元部分可以模拟其与周围岩体的相互作用，分析支护结构的受力情况和变形特征；离散元部分则可以考虑岩体中节理和裂隙对支护效果的影响，如节理和裂隙的存在可能导致支护结构的局部失效，通过离散元模拟可以直观地观察到这种失效过程。通过对支护过程的模拟，可以优化支护方案，选择合适的支护参数，提高隧道的稳定性。在昔格达地层隧道基底研究中，离散-连续方法的实现需要一系列严谨的步骤和关键技术。在模型建立阶段，需要对昔格达地层进行详细的地质勘察，获取地层的结构特征、岩石力学参数等信息。根据这些信息，合理划分离散元和连续元区域。对于昔格达地层中节理、裂隙较为发育的部分，划分为离散元区域，采用颗粒流方法模拟其离散特性；对于相对完整的岩体部分，划分为连续元区域，运用有限元方法进行模拟。确定离散元和连续元之间的耦合边界条件是实现离散-连续方法的关键技术之一。耦合边界条件需要保证离散元和连续元之间力和位移的连续传递。可以采用位移协调条件，即离散元区域和连续元区域在耦合边界上的位移相等；力的平衡条件，即离散元区域和连续元区域在耦合边界上的相互作用力大小相等、方向相反。通过这些条件的设定，确保离散元和连续元之间能够有效地进行信息交互。选择合适的耦合算法也是实现离散-连续方法的重要环节。常用的耦合算法有显式耦合算法和隐式耦合算法。显式耦合算法计算效率较高，适用于处理大规模问题，但精度相对较低；隐式耦合算法精度较高，但计算时间较长，对计算机性能要求较高。在昔格达地层隧道基底研究中，需要根据具体问题的规模和精度要求，选择合适的耦合算法。还可以采用一些优化算法，如并行计算技术，提高计算效率，缩短计算时间，以满足工程实际需求。四、昔格达地层铁路隧道基底动力特性试验研究4.1试验方案设计本次试验旨在深入探究昔格达地层铁路隧道基底在不同工况下的动力特性，为后续的数值模拟和工程应用提供可靠的数据支持。通过模拟实际工程中的各种受力条件和环境因素，分析昔格达地层的动应力-应变关系、动弹性模量、阻尼比等动力参数的变化规律，以及这些参数与各影响因素之间的内在联系。试验选用的主要设备为动三轴试验机，该设备能够精确控制轴向压力、围压、动应力幅值和加载频率等试验参数，满足本次试验对复杂加载条件的要求。动三轴试验机配备了高精度的传感器，可实时测量试样在加载过程中的轴向应变、径向应变和孔隙水压力等物理量，确保试验数据的准确性和可靠性。还使用了电子天平、烘箱、压力室等辅助设备，用于制备和测试昔格达地层试样的基本物理参数。为全面研究昔格达地层铁路隧道基底的动力特性，设计了多种试验工况，涵盖不同动应力、静偏应力、围压、频率、含水率等条件。在动应力方面，设置了50kPa、100kPa、150kPa、200kPa、250kPa五个等级，以模拟不同列车荷载强度下隧道基底的受力情况。静偏应力选取了0kPa、50kPa、100kPa三个水平，用于研究初始静应力对昔格达地层动力响应的影响。围压分别设定为100kPa、200kPa、300kPa，以反映不同深度隧道基底所受的上覆压力。加载频率设置为1Hz、2Hz、3Hz、4Hz、5Hz，模拟列车不同运行速度产生的振动频率。针对昔格达地层对含水率敏感的特性，在含水率方面，分别制备了天然含水率、饱和含水率以及在天然含水率基础上增减5%、10%的试样，深入研究含水率变化对其动力特性的影响。为确保试验结果的准确性和可靠性，每个工况均进行3次平行试验，取平均值作为该工况下的试验结果。各工况之间相互独立又具有一定的关联性，通过对比不同工况下的试验数据，能够系统地分析各因素对昔格达地层铁路隧道基底动力特性的影响机制。4.2试验结果与分析通过动三轴试验，获取了昔格达地层在不同试验工况下的应变随循环次数变化的数据，绘制了相应的曲线，如图1所示。从图中可以看出，在不同动应力、静偏应力、围压、频率和含水率条件下，昔格达地层的应变随循环次数的增加呈现出不同的变化趋势。图1昔格达地层应变与循环次数关系曲线在动应力对动力变形特性的影响方面，随着动应力的增大，应变增长速率明显加快，最终稳定应变值也显著增大。当动应力为50kPa时，经过1000次循环加载后，应变稳定在0.3%左右；而当动应力增大到250kPa时，相同循环次数下应变达到了1.2%以上。这表明昔格达地层在较大动应力作用下，更容易发生变形，且变形积累效应更为明显，对隧道基底的稳定性威胁更大。静偏应力对动力变形特性也有重要影响。在相同动应力和其他条件下，静偏应力越大，初始应变越大，且应变随循环次数的增长幅度也越大。当静偏应力为0kPa时，初始应变较小，且在循环加载过程中应变增长相对缓慢；而当静偏应力增加到100kPa时，初始应变明显增大，且在后续循环加载中应变增长迅速，这说明静偏应力的存在会改变昔格达地层的初始应力状态，使其在动荷载作用下更容易发生变形。围压对动力变形特性的影响则表现为围压增大，应变增长速率减小，最终稳定应变值降低。当围压为100kPa时，应变增长较快，稳定应变值较大；当围压增大到300kPa时，应变增长明显减缓，稳定应变值也降低了约30%。这是因为围压的增大限制了昔格达地层颗粒间的相对位移，增强了其抵抗变形的能力，从而减小了动力作用下的变形。加载频率对动力变形特性的影响较为复杂。在较低频率范围内（1Hz-3Hz），随着频率的增加，应变增长速率略有增加，稳定应变值也稍有增大；但当频率进一步增加到4Hz-5Hz时，应变增长速率反而减小，稳定应变值也有所降低。这可能是由于在较低频率下，加载时间相对较长，昔格达地层有足够的时间产生变形；而在较高频率下，加载时间短，地层来不及充分变形，且高频振动可能导致颗粒间的摩擦耗能增加，从而抑制了变形的发展。含水率对昔格达地层动力变形特性的影响最为显著。随着含水率的增加，应变增长速率急剧增大，最终稳定应变值大幅提高。天然含水率试样在试验条件下的稳定应变明显小于饱和含水率试样，饱和含水率试样的稳定应变值是天然含水率试样的2-3倍。这是因为昔格达地层遇水后，颗粒间的粘结力减弱，土体的抗剪强度降低，在动力荷载作用下更容易发生变形，且变形程度更大。为了进一步研究昔格达地层的动力变形特性，对动应力-应变关系进行了分析。根据试验数据，绘制了不同工况下的动应力-应变滞回曲线，如图2所示。从滞回曲线可以看出，动应力与应变之间呈现出非线性关系，滞回曲线所包围的面积代表了每次循环加载过程中的能量耗散。随着动应力的增大，滞回曲线的面积逐渐增大，表明能量耗散增加，地层的阻尼特性增强。在不同围压和含水率条件下，滞回曲线的形状和面积也有明显差异。围压增大时，滞回曲线更加饱满，能量耗散相对较小，说明围压有助于提高地层的抗变形能力，减少能量损耗；而含水率增加时，滞回曲线变得更加扁平，能量耗散显著增大，反映出含水率的增加会降低地层的强度和稳定性，使其在动力作用下更容易发生不可逆变形，消耗更多的能量。图2不同工况下昔格达地层动应力-应变滞回曲线通过对试验结果的深入分析，明确了动应力、静偏应力、围压、频率和含水率等因素对昔格达地层铁路隧道基底动力变形特性的影响规律，以及动应力-应变之间的非线性关系和能量耗散特性。这些结果为进一步理解昔格达地层在列车荷载作用下的力学行为，以及建立准确的数值模拟模型和长期沉降预测模型提供了重要的试验依据。4.3累积塑性应变模型建立累积塑性应变是描述土体在循环荷载作用下变形累积的重要参数，对于分析昔格达地层铁路隧道基底的长期稳定性具有关键作用。基于试验数据，对昔格达地层在不同工况下的累积塑性应变进行深入分析，旨在建立准确可靠的累积塑性应变模型，以预测隧道基底在长期列车荷载作用下的变形情况。通过对试验数据的细致观察和分析，发现昔格达地层的累积塑性应变与动应力幅值、循环次数等因素密切相关，且呈现出一定的规律性。累积塑性应变随着动应力幅值的增大而迅速增加，随着循环次数的增多也呈现出逐渐增长的趋势。在低动应力幅值下，累积塑性应变增长较为缓慢；而在高动应力幅值下，累积塑性应变增长速率明显加快。基于这些规律，采用经验公式法建立累积塑性应变模型。参考已有研究成果，并结合本次试验数据，选用如下形式的经验公式：\varepsilon_{p}=a(\frac{\sigma_{d}}{\sigma_{0}})^{b}N^{c}，其中\varepsilon_{p}为累积塑性应变，\sigma_{d}为动应力幅值，\sigma_{0}为参考应力（取为试验中的初始静偏应力），N为循环次数，a、b、c为模型参数。为确定模型参数a、b、c，运用非线性最小二乘法对试验数据进行拟合。将不同工况下的试验数据代入上述公式，通过优化算法不断调整参数值，使得公式计算结果与试验数据之间的误差平方和最小。经过反复计算和优化，得到在本次试验条件下，参数a=0.0015，b=1.8，c=0.35。为验证累积塑性应变模型的准确性和可靠性，将模型计算结果与试验数据进行对比分析。选取部分未参与模型参数拟合的试验工况数据，利用建立的模型计算其累积塑性应变，并与实际试验测量值进行对比，绘制对比曲线，如图3所示。从对比结果可以看出，模型计算值与试验测量值在趋势上基本一致，且大部分数据点的误差在可接受范围内。在低动应力幅值和较少循环次数情况下，模型计算值与试验值吻合较好；在高动应力幅值和较多循环次数时，虽然存在一定误差，但误差范围仍在合理区间内，能够满足工程实际应用的精度要求。这表明建立的累积塑性应变模型能够较好地反映昔格达地层在循环荷载作用下的累积塑性应变特性，为预测隧道基底的长期沉降提供了有力的工具。图3累积塑性应变模型计算值与试验值对比曲线五、基于离散-连续方法的隧道基底动力响应数值模拟5.1数值模型建立本研究选用国际上广泛应用的离散元软件PFC（ParticleFlowCode）与有限元软件ANSYS进行耦合，构建昔格达地层铁路隧道基底的离散-连续耦合数值模型。PFC软件在处理离散介质问题上具有独特优势，能够精确模拟颗粒间的相互作用和大变形行为；ANSYS软件则在连续介质力学分析方面功能强大，可准确计算结构的应力、应变分布。两者的耦合能够充分发挥各自优势，真实地反映昔格达地层铁路隧道基底的力学行为。模型尺寸的确定充分考虑了隧道的实际工程情况以及边界效应的影响。隧道采用单洞双线形式，内径为10m，外径为10.8m，衬砌厚度为0.4m。考虑到隧道基底的影响范围，模型在水平方向取隧道外径的5倍，即54m；在竖直方向取隧道外径的6倍，即64.8m，以确保边界条件对模型内部计算结果的影响可以忽略不计。在模型中，隧道顶部覆土厚度设定为20m，模拟实际工程中隧道的埋深情况。网格划分是数值模拟的关键环节，其质量直接影响计算结果的准确性和计算效率。对于离散元部分，将昔格达地层离散为颗粒单元，颗粒粒径根据实际地层的颗粒级配分布进行确定，通过设置合适的颗粒接触模型和参数，模拟颗粒间的粘结、摩擦等相互作用。在隧道周边和基底区域，适当加密颗粒单元，以提高计算精度；远离隧道的区域，颗粒单元适当稀疏，以减少计算量。对于连续元部分，采用ANSYS软件中的Solid185实体单元对隧道衬砌和周围连续介质进行网格划分。在隧道衬砌和基底附近，采用较小的单元尺寸，保证计算精度；在远离隧道的区域，逐渐增大单元尺寸，以提高计算效率。通过合理的网格划分，离散元与连续元区域在耦合边界处实现了良好的过渡，确保了力和位移的连续传递。材料参数的准确设置是数值模拟的基础。通过现场勘探和室内试验，获取昔格达地层的物理力学参数。昔格达地层的密度为2.2g/cm³，弹性模量为500MPa，泊松比为0.35，内摩擦角为25°，粘聚力为20kPa。考虑到昔格达地层遇水软化的特性，设置不同含水率工况下的材料参数，当含水率增加时，弹性模量和粘聚力相应降低。隧道衬砌采用C35混凝土，密度为2.5g/cm³，弹性模量为31.5GPa，泊松比为0.2。将这些材料参数准确输入到离散-连续耦合模型中，为后续的动力响应分析提供可靠的基础。5.2荷载与边界条件设定列车荷载是影响隧道基底动力响应的关键因素，其形式和参数的准确设定对数值模拟结果的准确性至关重要。在实际铁路运营中，列车荷载具有复杂性和多样性，受到列车类型、编组、速度、轴重等多种因素的影响。根据我国铁路的实际情况，参考相关规范和研究成果，将列车荷载简化为移动的均布荷载和集中荷载组合。对于常见的客货共线铁路列车，其荷载图式由普通荷载图式和特种荷载图式组成。普通荷载图式由集中荷载和均布荷载组成，特种荷载图式根据我国路情特点制定，由集中荷载组成。在本研究中，主要考虑普通荷载图式，均布荷载取值为92kN/m，集中荷载取值为320kN。列车速度设定为120km/h，这是该线路的设计运行速度，具有代表性。轴重根据实际情况选取25t，反映了常见货车的轴重水平。通过这些参数的设定，能够较为真实地模拟列车运行时对隧道基底产生的动力作用。在动力边界条件设定方面，为了准确模拟地震波在无限域介质中的传播，采用黏弹性边界条件。黏弹性边界条件能够有效吸收边界处的反射波，减少边界效应的影响，使数值模拟结果更加符合实际情况。在离散-连续耦合模型的边界上，设置黏弹性边界单元，通过调整边界单元的参数，如阻尼系数、弹簧刚度等，来实现对地震波的有效吸收。阻尼系数的取值根据昔格达地层的材料特性和波的传播特性确定，一般在0.1-0.5之间，本研究中取值为0.3；弹簧刚度根据边界单元的尺寸和材料参数计算得到，确保边界单元能够准确地模拟无限域介质的力学行为。考虑阻尼的影响对于准确模拟隧道基底的动力响应同样重要。阻尼能够消耗能量，使结构在振动过程中的能量逐渐衰减，从而影响结构的动力响应。在昔格达地层中，阻尼主要包括材料阻尼和辐射阻尼。材料阻尼是由于材料内部的摩擦和黏性等因素导致的能量耗散，辐射阻尼是由于波在传播过程中向周围介质辐射能量而引起的阻尼。采用瑞利阻尼模型来考虑材料阻尼的影响。瑞利阻尼模型通过两个阻尼参数α和β来描述，α与质量矩阵相关，β与刚度矩阵相关。根据昔格达地层的试验数据和相关研究，确定α=0.05，β=0.001，以准确反映昔格达地层的材料阻尼特性。辐射阻尼通过在边界条件中设置适当的吸收系数来考虑，确保波在传播到边界时能够被有效地吸收，减少反射波对计算结果的影响。通过合理设定列车荷载、动力边界条件以及考虑阻尼的影响，为后续的隧道基底动力响应数值模拟提供了准确的计算条件，能够更真实地反映昔格达地层铁路隧道基底在列车荷载作用下的动力特性。5.3模拟结果分析通过数值模拟，得到了隧道基底在列车荷载作用下的加速度响应云图，如图4所示。从云图中可以清晰地看出，隧道基底的加速度响应呈现出明显的非均匀分布特征。在列车荷载作用点附近，加速度响应最为强烈，随着与荷载作用点距离的增加，加速度逐渐衰减。在隧道基底的中心区域，加速度相对较大，而在边缘区域，加速度相对较小。这是因为列车荷载主要通过轨道传递到隧道基底的中心部分，使得中心区域受到的冲击力较大，而边缘区域受到的影响相对较小。图4隧道基底加速度响应云图为了更直观地分析加速度随时间的变化规律，选取隧道基底中心的一个监测点，绘制加速度时程曲线，如图5所示。从时程曲线可以看出，当列车接近监测点时，加速度逐渐增大，在列车经过监测点瞬间，加速度达到峰值，随后迅速衰减。列车速度对加速度峰值有显著影响，列车速度越高，加速度峰值越大。当列车速度为100km/h时，加速度峰值为0.5m/s²；当列车速度提高到120km/h时，加速度峰值增大到0.7m/s²。这表明列车速度的增加会导致隧道基底受到的动力冲击增大，对隧道基底的稳定性产生更大的威胁。图5隧道基底中心监测点加速度时程曲线隧道基底的动应力响应同样呈现出一定的规律。通过模拟得到动应力响应云图，如图6所示。在列车荷载作用下，隧道基底的动应力分布不均匀，在轨道下方的基底区域，动应力明显高于其他区域。这是因为轨道将列车荷载集中传递到基底的特定部位，使得这些部位承受较大的压力。动应力在基底中沿深度方向逐渐衰减，在基底表面动应力较大，随着深度的增加，动应力逐渐减小。在基底表面，动应力最大值可达1.2MPa，而在基底以下5m处，动应力减小到0.5MPa左右。图6隧道基底动应力响应云图为了进一步研究动应力的变化规律，绘制不同深度处动应力随时间的变化曲线，如图7所示。从曲线可以看出，动应力随着列车的运行呈现出周期性变化，列车经过时动应力迅速增大，列车离开后动应力逐渐恢复到初始状态。不同深度处动应力的变化幅度不同，深度越浅，动应力变化幅度越大；深度越深，动应力变化幅度越小。这是由于动应力在传播过程中会逐渐衰减，深度越大，衰减越明显。图7不同深度处隧道基底动应力时程曲线通过对不同工况下隧道基底动力响应的模拟分析，研究了列车速度、轴重、轨道不平顺等因素对动力响应的影响规律。列车速度的增加会显著增大隧道基底的加速度和动应力响应。随着列车速度从80km/h提高到160km/h，隧道基底的加速度峰值增大了1.5倍，动应力峰值增大了1.2倍。轴重的增加也会使隧道基底的动力响应增强，当轴重从20t增加到30t时，隧道基底的动应力峰值增大了0.5MPa。轨道不平顺对隧道基底动力响应的影响也不容忽视，轨道不平顺会导致列车荷载的不均匀分布，从而使隧道基底的加速度和动应力响应更加复杂，局部区域的动力响应明显增大。在轨道存在高低不平顺幅值为10mm的情况下，隧道基底局部区域的动应力峰值比平顺轨道时增大了30%左右。这些影响规律的揭示，为隧道的设计和运营管理提供了重要的参考依据，有助于采取相应的措施来减小隧道基底的动力响应，提高隧道的稳定性和安全性。六、昔格达地层铁路隧道基底长期沉降研究6.1长期沉降影响因素分析地质构造是影响昔格达地层铁路隧道基底长期沉降的重要因素之一。昔格达地层主要分布于中国青藏高原东南缘金沙江、雅砻江、安宁河、大渡河及其支流河谷中，该区域处于板块碰撞挤压的构造活跃地带，褶皱、断层等地质构造发育。在隧道建设和运营过程中，地质构造的活动会导致地层应力状态的改变，进而影响隧道基底的稳定性和沉降情况。在褶皱构造区域，隧道基底可能会受到不均匀的挤压和拉伸作用，导致基底土体产生不均匀变形，从而引起沉降。当隧道穿越向斜构造时，基底土体受到的挤压应力较大，容易发生压缩变形，导致沉降量增大；而在背斜构造区域，基底土体可能会因拉伸作用而出现裂缝，降低土体的承载能力，进而引发沉降。断层的存在对隧道基底沉降的影响更为显著。断层是地层中的薄弱带，其两侧的岩体往往存在错动和位移。当隧道基底位于断层附近时，由于断层的活动，基底土体可能会受到剪切力的作用，导致土体结构破坏，强度降低。在断层活动较为频繁的地区，隧道基底可能会出现突然的沉降或不均匀沉降，严重威胁隧道的结构安全和列车的运行安全。某铁路隧道在穿越断层时，由于断层的活动，隧道基底在短时间内出现了100mm的沉降，导致轨道变形，列车运行受阻。地下水对昔格达地层铁路隧道基底长期沉降的影响机制较为复杂。昔格达地层透水性弱，但在长期的地质作用和工程活动影响下，地下水的渗流和水位变化仍会对基底产生重要影响。地下水的渗流会带走基底土体中的细颗粒，导致土体孔隙增大，结构松散，从而降低土体的承载能力，引起沉降。在隧道施工过程中，由于排水措施不当，地下水渗流可能会导致基底土体的流失，形成空洞，进而引发基底沉降。某隧道在施工过程中，因排水系统不完善，地下水渗流带走了基底土体中的大量细颗粒，导致基底局部出现塌陷，沉降量达到了50mm。地下水水位的变化也会对隧道基底沉降产生影响。当水位上升时，基底土体处于饱水状态，其抗剪强度降低，在列车荷载和自重等作用下，更容易发生变形和沉降。水位上升还可能导致孔隙水压力增大，有效应力减小，进一步加剧基底的沉降。相反，当水位下降时，土体可能会因失水而收缩，产生沉降。在一些季节性水位变化明显的地区，隧道基底沉降会随着水位的升降而呈现出周期性变化。在雨季，地下水水位上升，隧道基底沉降量增大；在旱季，水位下降，沉降量相对减小，但土体的反复干湿循环会导致土体结构的破坏，长期来看，仍会加剧基底的沉降。随着时间的推移以及隧道运营过程中各种因素的作用，昔格达地层的岩石性质会发生变化，进而影响隧道基底的长期沉降。在长期的列车振动荷载作用下，昔格达地层岩石的结构会逐渐被破坏，颗粒间的连接变弱，导致岩石的强度降低。岩石的疲劳损伤会随着振动次数的增加而逐渐积累，使得岩石的弹性模量减小，变形能力增大。经过多年的列车运行后，隧道基底附近的昔格达地层岩石弹性模量可能会降低20%-30%，从而导致基底沉降量明显增大。风化作用也是导致岩石性质变化的重要因素。昔格达地层岩石长期暴露在自然环境中，受到风化作用的影响，岩石的化学成分和物理结构会发生改变。风化作用会使岩石中的矿物质分解、流失，导致岩石的强度和稳定性降低。在隧道洞口等部位，由于岩石直接暴露在大气中，风化作用更为强烈，基底沉降问题也更为突出。风化作用还会使岩石的吸水性增强，在地下水的作用下，更容易发生软化和崩解，进一步加剧基底的沉降。6.2沉降预测模型建立与验证基于离散-连续方法，结合昔格达地层的物理力学性质以及隧道基底在列车荷载等作用下的动力响应分析结果，建立隧道基底长期沉降预测模型。考虑到昔格达地层的蠕变特性对长期沉降的重要影响，采用Burgers模型来描述昔格达地层的蠕变行为。Burgers模型由Maxwell模型和Kelvin模型串联而成，能够较好地反映材料在长期荷载作用下的瞬时弹性变形、黏弹性变形和黏塑性变形。Maxwell模型由一个弹簧和一个黏壶串联组成，其本构方程为\dot{\varepsilon}=\frac{\dot{\sigma}}{E_{1}}+\frac{\sigma}{\eta_{1}}，其中\dot{\varepsilon}为应变率，\dot{\sigma}为应力率，E_{1}为弹簧的弹性模量，\eta_{1}为黏壶的黏性系数。Kelvin模型由一个弹簧和一个黏壶并联组成，本构方程为\sigma=E_{2}\varepsilon+\eta_{2}\dot{\varepsilon}，其中E_{2}为弹簧的弹性模量，\eta_{2}为黏壶的黏性系数。Burgers模型的总应变\varepsilon为Maxwell模型和Kelvin模型应变之和，即\varepsilon=\varepsilon_{M}+\varepsilon_{K}。在离散-连续耦合模型中，将Burgers模型引入离散元部分，用于描述昔格达地层颗粒间的蠕变特性。通过对离散元颗粒间接触力和相对位移的计算，结合Burgers模型的本构方程，得到昔格达地层在长期荷载作用下的变形情况。在连续元部分，根据弹性力学和塑性力学理论，计算隧道衬砌和周围连续介质的应力和应变分布，考虑到长期沉降过程中地层的非线性力学行为，采用合适的非线性本构模型，如Drucker-Prager模型，来描述地层的塑性变形。为了验证沉降预测模型的准确性，收集了成昆铁路复线峨眉至米易段范家湾隧道的现场监测数据。该隧道穿越昔格达地层，在施工和运营过程中对隧道基底沉降进行了长期监测。将现场监测数据与沉降预测模型的计算结果进行对比分析，绘制沉降随时间变化的对比曲线，如图8所示。图8沉降预测模型计算值与现场监测值对比曲线从对比曲线可以看出，沉降预测模型的计算结果与现场监测数据在趋势上基本一致，能够较好地反映隧道基底沉降随时间的变化规律。在隧道运营初期，沉降增长较快，随着时间的推移，沉降增长速率逐渐减小，最终趋于稳定。模型计算值与监测值在数值上也较为接近，在运营5年后，模型计算的沉降量为35mm，现场监测的沉降量为38mm，误差在10%以内，满足工程实际应用的精度要求。这表明建立的基于离散-连续方法的隧道基底长期沉降预测模型具有较高的准确性和可靠性，能够为昔格达地层铁路隧道的长期运营维护提供有效的决策依据。6.3长期沉降控制措施探讨在隧道设计阶段，优化支护结构是控制昔格达地层铁路隧道基底长期沉降的关键措施之一。针对昔格达地层强度低、遇水易软化等特性，应加强隧道衬砌结构的设计。增大衬砌厚度，采用高强度的混凝土材料，提高衬砌的承载能力和抗变形能力。将衬砌厚度从常规的40cm增加到50cm，选用C40混凝土代替C35混凝土，可有效增强衬砌对基底变形的抵抗能力。合理布置衬砌中的钢筋，增加钢筋的数量和直径，提高衬砌的抗拉强度，防止衬砌在基底沉降作用下出现开裂现象。在衬砌底部设置加强钢筋网，钢筋直径从16mm增大到20mm，间距从20cm减小到15cm，可显著提高衬砌底部的抗拉性能。采用新型的支护结构形式也是一种有效的优化方法。可考虑使用双层衬砌结构，外层衬砌主要承受初期的荷载和变形，内层衬砌则起到加强和保护的作用，进一步提高隧道的整体稳定性。在一些工程中，双层衬砌结构能够使隧道基底的沉降量降低20%-30%。还可以采用预应力锚索支护，通过对锚索施加预应力，将隧道基底的荷载传递到深部稳定的地层中，从而减小基底的沉降。在某隧道工程中，采用预应力锚索支护后，基底沉降得到了有效控制，沉降速率明显降低。加强地基处理是减小昔格达地层铁路隧道基底长期沉降的重要手段。注浆加固是一种常用的地基处理方法，通过向昔格达地层中注入水泥浆、化学浆液等，填充地层中的孔隙和裂隙，提高地层的强度和密实度。在隧道基底周边布置注浆孔，采用后退式分段注浆工艺，注浆压力控制在1-3MPa，可有效改善基底的力学性能。根据工程实践，注浆加固后，昔格达地层的弹性模量可提高30%-50%，地基承载力显著增强，从而减小基底的沉降量。强夯法也是一种有效的地基处理措施。利用重锤从高处自由落下，对地基进行强力夯实，使地基土密实，提高地基的承载能力。对于昔格达地层，强夯能级可根据地层的具体情况选择2000-4000kN・m，夯击次数一般为6-8击。在某铁路隧道工程中，对基底采用强夯法处理后，地基土的孔隙比减小了0.1-0.2，压缩模量提高了2-3MPa，基底沉降得到了有效控制。控制地下水对昔格达地层铁路隧道基底长期沉降的影响至关重要。完善隧道的排水系统，确保地下水能够及时、顺畅地排出。在隧道底部设置中心排水沟，排水沟的尺寸根据隧道的涌水量进行合理设计，一般宽度为0.5-1m，深度为0.8-1.2m，以保证能够有效排除基底积水。在隧道衬砌背后设置排水盲管，盲管的间距一般为5-10m，将衬砌背后的地下水引入中心排水沟，避免地下水在基底积聚，降低基底土体的强度。采取堵水措施，减少地下水对基底的侵蚀。在隧道周围设置止水帷幕，可采用高压旋喷桩、深层搅拌桩等方式形成止水帷幕，阻止地下水的渗入。止水帷幕的深度应根据地下水的水位和地层情况确定，一般应深入到不透水层以下1-2m。在某隧道工程中，设置止水帷幕后，基底的地下水位明显降低，基底土体的含水率保持在较低水平，有效减少了因地下水作用导致的基底沉降。七、结论与展望7.1研究成果总结通过对昔格达地层铁路隧道基底动力特性及长期沉降的深入研究，基于离散-连续方法，取得了以下具有重要理论和实践价值的成果：昔格达地层物理力学性质与动力特性：通过系统的现场勘探和室内试验，全面获取了昔格达地层的基本物理参数和力学参数。深入研究发现，含水率对昔格达地层的强度和变形特性影响显著。随着含水率的增加，地层的抗压强度、粘聚力等力学指标大幅降低，变形能力显著增强。在动三轴试验中，明确了