上软下硬节理岩体中地下洞室围岩松动特性与压力计算新探

上软下硬节理岩体中地下洞室围岩松动特性与压力计算新探一、引言1.1研究背景与意义随着现代工程建设的快速发展，地下洞室工程在交通、水利、能源等领域的应用越来越广泛。在这些工程中，上软下硬节理岩体的情况屡见不鲜。例如，在山区进行公路隧道建设时，常常会遇到上部为软弱的风化岩层，下部为坚硬的基岩的情况；在水利工程中的引水隧洞建设中，也可能穿越上软下硬的节理岩体地层。这种特殊的岩体结构给地下洞室的设计和施工带来了极大的挑战。地下洞室开挖后，围岩的应力状态发生改变，原本的平衡被打破。在重分布应力作用下，围岩可能产生过量的塑性变形或松动破坏，进而对支护衬砌施加压力，即围岩压力。围岩压力的准确计算对于地下洞室的稳定性至关重要。如果支护设计未能充分考虑围岩的实际松动特性，导致支护强度不足，可能引发洞室坍塌、地面沉降等严重事故，威胁施工人员的生命安全，造成巨大的经济损失。而如果支护过度设计，虽然保障了安全，但会大幅增加工程成本，造成资源浪费。对于上软下硬节理岩体，其力学性质具有明显的非均质性和各向异性。上部软岩和下部硬岩的强度、变形模量等力学参数差异显著，节理的存在又进一步削弱了岩体的整体性，使得围岩的变形破坏机制更加复杂。传统的基于连续介质力学和均质体假设的围岩压力计算方法，难以准确反映这种复杂岩体的实际力学行为，导致计算结果与实际情况偏差较大。因此，深入研究上软下硬节理岩体中地下洞室围岩的松动特性，建立更加合理准确的松动压力计算方法，具有重要的工程意义和现实需求。1.2国内外研究现状在地下洞室围岩稳定性研究领域，国内外学者已取得了一系列成果。在围岩松动特性方面，早期研究主要聚焦于均质岩体，随着工程实践中遇到的地质条件愈发复杂，节理岩体的研究逐渐成为热点。国外学者如Goodman早在20世纪60年代就开始关注节理对岩体力学性质的影响，提出了节理岩体的力学模型，为后续研究奠定了基础。之后，Barton等通过大量试验，建立了节理岩体的抗剪强度经验公式，考虑了节理粗糙度、充填物等因素。在数值模拟方面，离散元法（DEM）、非连续变形分析（DDA）等方法被广泛应用于节理岩体的变形破坏模拟。如Cundall开发的离散元程序UDEC和3DEC，能够较好地模拟节理岩体在荷载作用下的块体运动和相互作用。国内学者在节理岩体研究方面也做出了重要贡献。孙广忠提出了岩体结构控制论，强调了岩体结构对岩体力学行为的控制作用。冯夏庭等运用岩石破裂过程分析（RFPA）系统，对节理岩体的破裂过程进行了数值模拟，揭示了节理岩体在不同荷载条件下的破坏机制。在现场测试方面，采用声波测试、钻孔摄像等技术手段，对围岩松动圈范围进行测定，为理论研究提供了数据支持。对于上软下硬节理岩体中地下洞室围岩压力计算方法，传统方法主要有平衡拱理论、太沙基理论等。平衡拱理论由普氏提出，认为洞室开挖后顶部会形成自然平衡拱，作用于支护上的压力为平衡拱内岩体重量。太沙基理论则基于土力学原理，考虑了围岩的剪切破坏和荷载传递。然而，这些经典理论在应用于上软下硬节理岩体时存在局限性。它们大多基于连续介质假设，未充分考虑节理的存在对岩体力学性质的弱化作用以及上下岩体力学参数差异导致的应力分布不均匀性。近年来，国内外学者针对上软下硬节理岩体的特性，对传统计算方法进行了修正和改进。王志伟等考虑节理面的分布特征和力学性质，采用极限分析上限法推导了适用于节理岩体的松动压力计算公式。但这些改进方法在实际应用中仍面临挑战，如对节理参数的准确获取较为困难，计算模型的复杂性增加了计算成本和不确定性。现有研究在围岩松动特性和压力计算方法上取得了一定进展，但仍存在不足。对于上软下硬节理岩体这种复杂地质条件，围岩变形破坏机制尚未完全明确，现有计算方法难以全面准确地反映其力学行为，在实际工程应用中存在局限性，需要进一步深入研究。1.3研究内容与方法1.3.1研究内容本研究旨在深入剖析上软下硬节理岩体中地下洞室围岩的松动特性，并建立更为准确的松动压力计算方法。具体研究内容如下：上软下硬节理岩体特性研究：通过现场地质调查、室内岩石力学试验以及微观结构分析等手段，获取上软下硬节理岩体的物理力学参数，包括软岩和硬岩的弹性模量、泊松比、抗压强度、抗剪强度，以及节理的产状、密度、粗糙度、充填物性质等。分析节理的几何特征和力学性质对岩体整体力学行为的影响规律，研究软岩与硬岩之间的相互作用机制，建立考虑节理影响的上软下硬节理岩体力学模型。地下洞室围岩松动特性研究：采用数值模拟方法，运用有限元软件（如ANSYS、ABAQUS）和离散元软件（如UDEC、3DEC），模拟不同开挖方式（如台阶法、CD法、CRD法）和支护条件下，上软下硬节理岩体中地下洞室围岩的应力分布、变形规律和破坏模式。结合现场监测数据，通过在工程现场布置位移监测点、应力监测元件等，实时监测洞室开挖过程中围岩的位移、应力变化情况，验证数值模拟结果的准确性。分析围岩松动区的形成机制、发展过程和范围大小，明确影响围岩松动特性的关键因素，如岩体力学参数、洞室埋深、洞室形状、节理分布等。地下洞室围岩松动压力计算方法研究：基于对围岩松动特性的研究成果，考虑上软下硬节理岩体的非均质性、各向异性以及节理的影响，对传统的围岩压力计算方法进行修正和改进。建立适用于上软下硬节理岩体的围岩松动压力计算模型，推导相应的计算公式，充分考虑岩体的非线性力学行为、节理面的剪切滑移和拉裂破坏等因素。通过工程实例分析，将建立的计算方法应用于实际工程案例中，与现场实测数据和其他计算方法的结果进行对比验证，评估其准确性和可靠性，进一步优化和完善计算方法。1.3.2研究方法为实现上述研究目标，本研究将综合运用多种研究方法，相互验证和补充，确保研究结果的科学性和可靠性。数值模拟方法：利用数值模拟软件建立上软下硬节理岩体中地下洞室的三维模型，模拟洞室开挖和支护过程。通过设置不同的岩体参数、节理条件和施工工况，系统分析围岩的力学响应，预测围岩的变形、破坏模式和松动压力分布。数值模拟方法能够直观地展示复杂地质条件下围岩的力学行为，为理论分析和现场试验提供参考依据，且可灵活改变参数，进行多方案对比研究。理论分析方法：基于岩石力学、弹性力学、塑性力学等基本理论，分析上软下硬节理岩体中地下洞室围岩的应力重分布规律和破坏准则。推导考虑节理影响的围岩松动压力计算公式，从理论层面揭示围岩松动压力的形成机制和影响因素之间的内在关系。理论分析方法具有普遍性和通用性，能够为数值模拟和工程实践提供理论指导。工程实例分析方法：选取典型的上软下硬节理岩体中地下洞室工程案例，收集现场地质勘察资料、施工监测数据和支护设计方案。将数值模拟结果和理论计算结果与现场实际情况进行对比分析，验证研究成果的可行性和有效性。通过工程实例分析，还能发现实际工程中存在的问题，进一步完善研究内容和方法，使研究成果更具工程应用价值。二、上软下硬节理岩体特性分析2.1节理岩体的基本特征节理作为岩体中常见的一种地质构造，是指岩石受力断裂后，两侧岩块没有发生显著位移的小型断裂构造。节理在自然界中广泛分布，其形态、尺寸、空间分布及组合特征呈现出高度的复杂性和多样性。从形态上看，节理面有的较为平整光滑，有的则粗糙起伏，甚至呈现出锯齿状；尺寸方面，节理的长度可从几厘米到数千米不等，宽度也可在微小的裂隙至数米的张开度之间变化。节理的分类方式丰富多样，依据不同的标准可划分为不同的类型。按照成因来分，节理可分为原生节理和次生节理。原生节理是在岩石成岩过程中形成的，例如沉积岩在成岩阶段因脱水收缩产生的泥裂，或是岩浆岩在冷凝收缩时形成的柱状节理。次生节理则是岩石成岩后，受后续地质作用影响而形成的，涵盖了非构造节理（如风化节理）和构造节理。其中，构造节理最为常见，它是在地质构造运动作用下，岩石受力变形产生的节理，常成组出现，与当地的褶皱、断层构造紧密相关，在空间分布上具有一定的规律性。根据节理与岩层产状要素的关系，又可将节理分为走向节理、倾向节理、斜向节理和顺层节理。走向节理的走向与岩层的走向一致或大致相同；倾向节理的走向大致与岩层的走向垂直，即与岩层的倾向一致；斜向节理的走向与岩层的走向既不平行也不垂直，而是呈斜交状态；顺层节理的节理面大致平行于岩层层面。此外，按照节理的走向与区域褶皱主要方向、断层的主要走向或其他线形构造的延伸方向等关系，还能划分为纵节理、横节理和斜节理。当褶皱轴延伸稳定，不发生倾伏（水平褶皱）时，走向节理相当于纵节理，倾向节理相当于横节理，斜向节理相当于斜节理。从力学性质角度，构造节理可分为伸张节理（张节理）和剪切节理。伸张节理是垂直于最小主应力轴方向（平行于最大主应力轴方向）发生张裂所形成的节理，常发育在褶皱构造中背斜及向斜的轴部，且垂直于褶皱轴面。伸张节理一般不平直，裂面较粗糙，间距大且分布不均匀，很少密集成带，沿其走向及倾向都延伸不远，两壁常张开，剖面上呈上宽下窄的楔形，裂隙宽窄变化较大，易被造岩矿物充填，在砾岩中发育的伸张节理常绕过砾石，且往往成为地下水的良好通道。剪切节理则是岩石因剪切破坏而形成的节理，多出现在褶皱的两翼或断层附近，沿节理走向及倾向延伸较远，发育的密度较大，节理间距小而具有等间距性，沿节理面可能发生轻微位移，节理面平直而光滑，裂隙闭合，有时留有擦痕及磨光面，常沿剪切面成群的平行分布，通常有两组不同方向的节理，交叉成X型，称为共轭节理，将岩体切割成菱块状或豆腐状，最大主应力轴即平分共轭节理的锐角，在砾岩中发育的剪切节理常切穿砾石，而不改变其方向。节理的存在对岩体的力学性质产生了多方面的显著影响。在强度方面，节理使得岩体的强度大幅降低，其强度介于岩石强度及不连续面（即节理面）强度之间。当岩体沿节理面整体滑动破坏时，其强度取决于节理面的强度。节理的存在还导致岩体力学性质呈现出明显的不连续性、不均匀性及各向异性。由于节理的切割，岩体被分割成形状不同、大小不等的岩块，使得岩块的强度与包含节理的岩体强度存在明显差异。在变形特性上，节理降低了岩体的整体刚度，使得岩体在受力时更容易产生变形。在相同的荷载作用下，节理发育的岩体变形量往往大于完整岩体。而且，由于节理在岩体中的分布具有方向性，导致岩体在不同方向上的变形特性也存在差异，呈现出各向异性。节理对岩体的渗透性影响也不容忽视，节理常成为地下水的运移通道，使得岩体的渗透性增强。节理的连通性越好，岩体的渗透性就越强。在水利工程、地下工程等建设中，需要充分考虑节理对岩体渗透性的影响，以防止地下水渗漏等问题的发生。2.2上软下硬地层结构特点上软下硬地层是一种典型的非均质复合地层，其结构特点主要体现在软硬岩层的组合形式、力学性质差异以及在不同地质构造中的分布等方面。从组合形式来看，上软下硬地层常见的形式有上部为软土层，下部为硬岩层，两者呈明显的上下叠置状态。软土层可能是各种成因的黏土、粉质黏土、砂土等，如冲积、洪积、湖积等成因的土层。这些软土通常具有含水量高、孔隙比大、压缩性高、强度低等特点。例如，淤泥质黏土的含水量可达50%-70%，孔隙比在1.0-1.5之间，压缩系数大于0.5MPa⁻¹，抗剪强度较低。硬岩层则可能是花岗岩、石灰岩、砂岩等各类岩石，其结构致密，强度较高。以花岗岩为例，其饱和单轴抗压强度一般在50-200MPa之间，弹性模量较大，可达10-100GPa。在一些地层中，软硬岩层并非简单的上下叠置，还可能出现软硬岩层互层的情况，即软岩和硬岩呈薄层相间分布。这种互层结构使得地层的力学性质在垂向上呈现出更为复杂的变化，软岩和硬岩的相互约束和影响更为显著。软硬岩层的力学差异是上软下硬地层的关键特征。在强度方面，软岩的抗压强度、抗剪强度远低于硬岩。软岩的抗压强度一般在几MPa以下，抗剪强度也仅为几十kPa到几百kPa。而硬岩的抗压强度通常在几十MPa以上，抗剪强度可达几MPa甚至更高。这种强度差异导致在地下洞室开挖过程中，软岩更容易发生塑性变形和破坏，而硬岩则相对稳定。在变形特性上，软岩的弹性模量较小，受力时容易产生较大的变形。软岩的弹性模量可能在0.1-1GPa之间，在相同荷载作用下，其变形量可能是硬岩的数倍甚至数十倍。硬岩的弹性模量较大，变形相对较小。软硬岩层的泊松比也存在差异，这进一步影响了它们在受力时的横向变形特性。在渗透性方面，软岩和硬岩也有明显不同。软岩的渗透性一般较差，地下水在其中的运移速度较慢。而硬岩的渗透性则因岩石类型和节理发育程度而异，一些节理发育的硬岩可能具有较好的渗透性。上软下硬地层在不同地质构造中有着不同的分布特征。在褶皱构造区域，上软下硬地层可能出现在褶皱的翼部或轴部。在翼部，由于岩层受到挤压，软硬岩层的界面可能发生倾斜，增加了地下洞室施工的难度。在轴部，由于应力集中，软岩更容易发生变形和破坏。在断层附近，上软下硬地层的分布可能受到断层活动的影响。断层的错动可能导致软硬岩层的相对位置发生改变，形成复杂的地质结构。在一些沉积盆地边缘，由于沉积环境的变化，也容易形成上软下硬的地层结构。河流相沉积的粗粒砂岩等硬岩与湖相沉积的软黏土可能交替出现。2.3上软下硬节理岩体力学特性研究为深入探究上软下硬节理岩体的力学特性，本研究综合运用室内试验和数值模拟两种方法，全面分析该岩体在不同荷载条件下的变形、破坏机制及力学参数变化规律。在室内试验方面，首先对取自工程现场的上软下硬节理岩体进行加工，制作成标准岩石试件。针对软岩和硬岩分别开展单轴抗压试验、三轴抗压试验、直接剪切试验等。单轴抗压试验在万能材料试验机上进行，通过缓慢施加轴向压力，记录试件的变形和破坏过程，从而得到软岩和硬岩的单轴抗压强度。某软岩试件在单轴抗压试验中，当轴向压力达到2.5MPa时，试件出现明显的纵向裂缝，随后迅速破坏，测得其单轴抗压强度为2.3MPa。而硬岩试件在相同试验条件下，承受的轴向压力达到80MPa才发生破坏，单轴抗压强度远高于软岩。三轴抗压试验则是在三轴压力室中进行，通过改变围压和轴压，研究软岩和硬岩在不同应力状态下的力学响应。在围压为5MPa的三轴抗压试验中，软岩的抗压强度提高到4.0MPa，表明围压对软岩强度有一定的增强作用。直接剪切试验用于测定岩体的抗剪强度，通过在剪切面上施加垂直荷载和水平剪切力，得到软岩和硬岩的抗剪强度指标（黏聚力和内摩擦角）。试验结果显示，软岩的黏聚力为10kPa，内摩擦角为20°；硬岩的黏聚力为500kPa，内摩擦角为35°，硬岩的抗剪强度明显高于软岩。节理的存在对岩体力学性质影响显著，因此开展节理岩体的力学试验十分必要。采用含有预制节理的岩石试件，研究节理倾角、节理间距、节理粗糙度等因素对岩体强度和变形的影响。当节理倾角为45°时，节理岩体的抗压强度最低，说明该倾角下节理对岩体强度的削弱作用最为明显。随着节理间距的减小，岩体的完整性降低，强度也随之下降。节理粗糙度越大，节理面之间的咬合作用越强，岩体的抗剪强度越高。数值模拟方面，选用UDEC离散元软件建立上软下硬节理岩体的数值模型。在模型中，准确模拟软岩、硬岩以及节理的力学行为。软岩和硬岩采用弹塑性本构模型，考虑其非线性力学特性。节理则采用节理单元进行模拟，通过设置节理的法向刚度、切向刚度、黏聚力、内摩擦角等参数，反映节理的力学性质。利用建立的数值模型，模拟不同荷载条件下岩体的力学响应。在施加均布荷载时，观察岩体的应力分布和变形情况。数值模拟结果表明，软岩部分的应力集中现象较为明显，变形量较大；硬岩部分的应力分布相对均匀，变形量较小。在水平荷载作用下，节理面容易产生剪切滑移，导致岩体的整体性受到破坏。通过改变节理参数和荷载条件，分析各因素对岩体力学特性的影响。当节理密度增加时，岩体的弹性模量降低，变形增大，说明节理的存在降低了岩体的整体刚度。通过室内试验和数值模拟结果的对比分析，验证数值模型的准确性和可靠性。两种方法得到的岩体力学参数和变形破坏规律基本一致，表明数值模拟能够较好地反映上软下硬节理岩体的力学特性。本研究为后续地下洞室围岩松动特性和松动压力计算方法的研究提供了坚实的力学基础。三、地下洞室围岩松动特性研究3.1围岩松动的概念与形成机制围岩松动是指地下洞室开挖后，由于开挖扰动、应力重分布等因素的影响，洞室周边一定范围内的岩体原有的连续性和完整性遭到破坏，岩体力学性质发生劣化，呈现出松弛破碎的状态。这部分岩体在工程实践中被称为围岩松动圈，其范围和特性对地下洞室的稳定性和支护设计具有重要影响。地下洞室开挖是一个复杂的力学过程，会对围岩的应力状态产生显著的改变。在开挖前，岩体处于初始应力场中，各点的应力处于相对平衡状态。以深埋洞室为例，岩体在自重应力和构造应力的作用下，内部应力分布较为均匀。当进行洞室开挖时，洞室周边的岩体失去了原有的支撑，应力平衡被打破。原本作用在开挖区域的应力将重新分布，向洞室周边的岩体转移，导致洞室周边一定范围内的岩体应力显著增加。在圆形洞室开挖的情况下，根据弹性力学理论，洞室周边的切向应力会增大到初始应力的2倍。这种应力的重分布是围岩松动的重要诱因。随着应力的重分布，洞室周边岩体所承受的应力逐渐增大。当应力超过岩体的强度极限时，岩体就会发生破坏。在脆性岩体中，当洞室周边的切向应力超过岩体的抗拉强度时，岩体就会产生张裂缝。这些裂缝首先在洞室周边出现，随着应力的持续作用，裂缝会逐渐向岩体内部扩展。如果应力进一步增大，超过岩体的抗压强度，岩体就会发生破碎。在塑性岩体中，当应力超过岩体的屈服强度时，岩体就会产生塑性变形。塑性变形会导致岩体的体积膨胀，内部结构变得松散。随着塑性变形的不断发展，岩体的强度不断降低，最终也会导致岩体的破坏。节理等结构面的存在对围岩松动的形成和发展具有重要影响。节理是岩体中的不连续面，其存在使得岩体的力学性质变得不均匀和各向异性。在洞室开挖过程中，节理面容易成为应力集中的部位。当节理面附近的应力超过节理面的抗剪强度时，节理面就会发生剪切滑移。这种剪切滑移会导致节理面两侧的岩体发生相对位移，进一步破坏岩体的整体性。节理的存在还会降低岩体的强度，使得岩体更容易在应力作用下发生破坏。在一组密集节理发育的岩体中，由于节理的切割，岩体被分割成许多小块，这些小块之间的连接较弱，在洞室开挖后的应力作用下，很容易发生松动和坍塌。节理的产状、密度、粗糙度等因素都会影响围岩的松动特性。节理的倾角越大，在一定应力条件下，节理面越容易发生滑动；节理密度越大，岩体的完整性越差，越容易发生松动。3.2围岩松动过程分析为深入剖析上软下硬节理岩体中地下洞室围岩的松动过程，本研究借助数值模拟软件UDEC，构建了精细的三维数值模型。该模型充分考虑了上软下硬地层结构、节理分布以及洞室开挖等关键因素，力求真实还原工程实际情况。模型中，上软下硬节理岩体的力学参数依据前文的室内试验结果进行准确赋值。软岩部分采用Mohr-Coulomb弹塑性本构模型，其弹性模量设为0.5GPa，泊松比为0.35，抗压强度为3MPa，抗剪强度指标黏聚力为15kPa，内摩擦角为22°。硬岩部分同样采用Mohr-Coulomb弹塑性本构模型，弹性模量为15GPa，泊松比为0.25，抗压强度为100MPa，黏聚力为500kPa，内摩擦角为38°。节理采用节理单元进行模拟，根据节理的粗糙度、充填物等实际情况，设置节理的法向刚度为5GPa/m，切向刚度为2GPa/m，黏聚力为5kPa，内摩擦角为18°。地下洞室设定为圆形断面，直径为5m，埋深为200m。模拟过程中，严格按照台阶法进行开挖，分上台阶和下台阶两步开挖，每步开挖后及时施作初期支护，初期支护采用弹性梁单元模拟，弹性模量为30GPa，截面积为0.1m²。3.2.1开挖初期在洞室开挖初期，上台阶开挖完成后，洞室周边围岩的应力状态发生了显著变化。由于开挖卸荷，洞室周边的应力迅速重新分布。在软岩部分，由于其强度较低，应力集中现象较为明显。洞室拱顶和拱腰部位的切向应力急剧增大，超过了软岩的抗拉强度，导致这些部位首先出现张拉裂缝。数值模拟结果显示，此时拱顶的切向应力达到了4.5MPa，超过了软岩的抗拉强度3MPa，从而产生了张拉裂缝。这些裂缝的出现使得软岩的完整性受到破坏，岩体开始向洞室内产生微小的位移。在硬岩部分，由于其强度较高，应力集中程度相对较小，但在洞室周边也出现了一定程度的应力升高。节理面在开挖扰动下，部分节理的应力状态也发生改变，一些节理面上的剪应力增大，当超过节理面的抗剪强度时，节理面开始出现微小的剪切滑移。在一组倾角为45°的节理面上，剪应力达到了8kPa，超过了节理面的抗剪强度5kPa，从而发生了剪切滑移。这种节理面的剪切滑移进一步削弱了岩体的整体性。3.2.2开挖中期随着下台阶的开挖，围岩的松动范围进一步扩大。软岩部分的裂缝继续向岩体内部扩展，形成了一个较为连续的松动区域。由于上台阶开挖后岩体已经发生了一定程度的变形和破坏，下台阶开挖时，上台阶松动区域的岩体对下台阶围岩产生了一定的压力，使得下台阶周边围岩的应力状态更加复杂。在软岩和硬岩的交界处，由于两者力学性质的差异，出现了明显的应力集中和变形不协调现象。软岩的变形较大，而硬岩的变形较小，导致交界处的岩体产生了较大的剪应力，进一步加剧了岩体的破坏。数值模拟结果表明，此时软岩的松动范围已经达到了洞室半径的1.5倍左右。节理面的剪切滑移现象更加普遍，一些节理面的滑移量增大，导致岩体的块体运动加剧。部分岩块之间的连接被破坏，开始出现松动和掉落的趋势。在洞室边墙部位，由于节理面的切割，一些岩块形成了不稳定的块体，在重力和围岩压力的作用下，有向洞室内滑落的危险。3.2.3开挖后期开挖完成后，经过一段时间的变形发展，围岩的松动逐渐趋于稳定，但松动范围已经基本确定。此时，软岩部分形成了一个明显的松动圈，松动圈内的岩体破碎程度较高，力学性质显著劣化。硬岩部分虽然整体较为稳定，但在洞室周边也存在一定范围的松动区域，主要是由于节理的影响和开挖扰动导致的。在洞室拱顶，松动圈的厚度达到了2.5m左右，软岩已经完全破碎，呈现出松散的状态。在硬岩中，靠近洞室周边的节理密集区域，也出现了一定程度的破碎和松动。松动圈的形成使得围岩的自承载能力大幅降低，大部分荷载将传递到支护结构上。如果支护结构的强度和刚度不足，将难以承受这些荷载，从而导致洞室的失稳。因此，准确掌握围岩松动过程和松动范围，对于合理设计支护结构具有至关重要的意义。通过对不同开挖阶段围岩松动过程的分析，可以看出岩体力学参数、节理分布以及开挖方式等因素对围岩松动特性有着显著的影响。在实际工程中，需要根据具体的地质条件和施工要求，采取合理的措施来控制围岩的松动，确保洞室的稳定。3.3围岩破坏模式分析上软下硬节理岩体中地下洞室围岩在复杂的地质条件和开挖扰动作用下，会呈现出多种破坏模式，每种破坏模式都有其独特的形成条件和特征。坍塌破坏是较为常见的一种模式，通常发生在洞室顶部。当洞室开挖后，顶部围岩在重力和围岩压力的作用下，由于岩体强度不足或结构面的切割，导致岩体失去稳定性，进而发生块状或碎块状的塌落。在软岩部分，由于其强度较低，当洞室顶部的拉应力超过软岩的抗拉强度时，岩体就会产生裂缝，随着裂缝的扩展和连通，岩体逐渐破碎，最终导致坍塌。在硬岩部分，如果节理发育且节理面的组合不利于岩体的稳定，在开挖扰动下，节理面之间的咬合作用被破坏，岩块之间的连接减弱，也容易发生坍塌。在一组垂直节理和水平节理相互切割的硬岩区域，洞室开挖后，顶部的岩块在重力作用下，容易沿着节理面掉落，形成坍塌。坍塌破坏的范围和规模与岩体的力学性质、节理分布、洞室跨度等因素密切相关。洞室跨度越大，顶部围岩所承受的荷载就越大，坍塌的可能性和规模也就越大。滑移破坏主要是由于岩体沿着节理面或软弱结构面发生相对滑动而导致的。节理面的抗剪强度较低，当洞室开挖后，围岩的应力状态发生改变，节理面上的剪应力增大，当超过节理面的抗剪强度时，岩体就会沿着节理面发生滑移。在软岩和硬岩的交界处，由于两者力学性质的差异，容易产生应力集中，使得交界处的节理面更容易发生滑移。在倾角较大的节理面上，岩体在重力的切向分力作用下，也容易发生滑移。滑移破坏可能导致洞室边墙的岩体向洞内移动，影响洞室的稳定性。如果边墙的岩体发生较大的滑移，可能会使洞室的形状发生改变，增加支护的难度。张裂破坏则是在拉应力作用下，岩体产生拉伸裂缝，随着裂缝的扩展和贯通，岩体发生破坏。洞室开挖后，在洞室周边的某些部位，如拱顶、边墙的拐角处等，会产生应力集中，当拉应力超过岩体的抗拉强度时，就会产生张裂破坏。在软岩中，由于抗拉强度较低，更容易发生张裂破坏。软岩的抗拉强度一般在0.1-1MPa之间，而洞室开挖后，周边的拉应力可能达到数MPa，远超过软岩的抗拉强度。张裂破坏会削弱岩体的整体性，使得岩体更容易受到其他破坏模式的影响。张裂破坏产生的裂缝会成为地下水的通道，加速岩体的风化和侵蚀，进一步降低岩体的强度。3.4围岩松动范围确定方法及松动区分类准确确定围岩松动范围对于评估地下洞室的稳定性和设计合理的支护方案至关重要。目前，确定围岩松动范围的方法主要有声波测试法、数值模拟法和钻孔摄像法等。声波测试法是一种较为常用的方法，其原理基于岩体的声波传播特性。在完整岩体中，声波传播速度相对较高，且传播过程中能量衰减较小。而当岩体出现松动、裂隙等损伤时，岩体的完整性被破坏，声波在其中传播时会发生散射、反射等现象，导致声波传播速度降低，能量衰减增大。通过在洞室周边钻孔，利用声波测试仪向岩体中发射声波，并接收反射或透射回来的声波信号，根据声波速度和能量的变化情况，就可以判断岩体的松动范围。当声波速度降低到一定程度时，对应的钻孔深度位置就可以认为是围岩松动圈的边界。声波测试法具有操作相对简单、测试速度较快、对岩体损伤较小等优点，但其测试结果受岩体岩性、节理发育程度等因素的影响较大。在节理密集发育的岩体中，声波的传播路径变得复杂，可能会导致测试结果的误差增大。数值模拟法借助有限元、离散元等数值计算软件，建立地下洞室及围岩的数值模型，模拟洞室开挖过程中围岩的应力、应变和破坏情况，从而预测围岩松动范围。在数值模拟中，根据岩体的物理力学参数、节理分布等实际情况，合理设置模型参数，采用合适的本构模型来描述岩体的力学行为。通过模拟不同开挖阶段和工况下围岩的力学响应，分析岩体的塑性区范围、位移分布等指标，确定围岩松动范围。数值模拟法能够综合考虑多种因素对围岩松动的影响，直观地展示围岩松动的发展过程，为支护设计提供参考依据。但数值模拟结果的准确性依赖于模型参数的合理性和本构模型的适用性，如果参数设置不合理或本构模型选择不当，可能导致模拟结果与实际情况偏差较大。钻孔摄像法是直接将钻孔摄像机放入钻孔中，对钻孔壁的岩体进行拍摄，通过观察拍摄到的图像，直观地判断岩体的完整性、裂隙发育情况和松动范围。在钻孔摄像图像中，完整岩体的图像表现为连续、均匀的岩石结构，而松动区的岩体则会出现裂隙、破碎等现象。钻孔摄像法能够提供直观、准确的岩体松动信息，对于判断围岩松动的具体形态和范围具有重要作用。但其操作相对复杂，成本较高，且只能获取钻孔位置处的岩体信息，对于钻孔之间的岩体松动情况难以准确判断。根据围岩松动范围的大小和岩体的破坏程度，可以将松动区分为小松动圈、中松动圈和大松动圈。当围岩松动圈厚度值小于40cm时，为小松动圈。在小松动圈围岩中，松动范围较小，岩体的破坏程度较轻，围岩的稳定性相对较好。由于松动圈厚度小，围岩产生的碎胀变形量也较小，一般只有几个毫米，这种变形量通常小于低应力下锚杆的弹塑性变形。因此，在小松动圈围岩中，只用喷射混凝土支护就能够保证工程的安全，一般不需要采用锚杆支护或其他普通支护形式。当围岩松动圈厚度值在40-150cm之间时，属于中松动圈。中松动圈围岩的碎胀变形比较明显，变形量较大。由于围岩松动圈的碎胀变形，刚性的喷射混凝土支护容易产生裂缝或破坏。所以，在中松动圈围岩中，必须采用以锚杆为主体构件的锚喷支护方式。锚杆可以有效地控制围岩的碎胀变形，喷层则主要起到支护锚杆间活石和防止围岩风化的作用。由于围岩松动圈厚度小于常用锚杆长度，可采用锚杆悬吊作用机理来设计支护参数，锚杆支护的最大荷载是围岩松动圈形成中的碎胀变形力及已形成松动圈内破裂岩石的自重。当围岩松动圈厚度值大于150cm时，为大松动圈。在大松动圈围岩巷道中，围岩表现出软岩的工程特征，松动圈碎胀变形量大，初期围岩收敛变形速度快，变形持续时间长，矿压显现较大，支护难度大。在这种情况下，通常需要采用联合支护形式，如“锚喷网架碹”等。通过多种支护方式的协同作用，共同抵抗围岩的变形和压力，确保洞室的稳定。3.5松动区影响因素的定量分析为深入剖析各因素对围岩松动区范围的影响程度，本研究运用数据包络分析（DEA）方法进行定量分析。DEA是一种基于线性规划的多投入多产出效率评价方法，能够有效处理多因素复杂系统，确定各因素对目标变量的影响权重。在确定影响因素时，综合考虑岩体力学参数、洞室几何参数和外部环境因素等多个方面。岩体力学参数选取软岩弹性模量、软岩泊松比、软岩抗压强度、硬岩弹性模量、硬岩泊松比、硬岩抗压强度、节理密度、节理粗糙度等；洞室几何参数包括洞室埋深、洞室跨度、洞室高度；外部环境因素考虑地下水压力。这些因素在地下洞室开挖过程中，都可能对围岩的应力分布和变形特性产生影响，进而影响松动区的范围。以某实际工程为例，收集了多个监测断面的数据作为样本。每个样本包含上述各影响因素的具体数值以及对应的围岩松动区范围实测值。将这些数据整理后，构建DEA模型的输入输出矩阵。其中，各影响因素作为输入指标，围岩松动区范围作为输出指标。通过DEA模型计算，得到各影响因素对围岩松动区范围的影响权重。计算结果表明，在岩体力学参数中，软岩抗压强度和硬岩弹性模量的影响权重较大，分别为0.25和0.20。这表明软岩抗压强度越低，在洞室开挖后的应力作用下，软岩越容易发生破坏，从而导致松动区范围扩大；硬岩弹性模量越大，其对软岩的约束作用越强，能够在一定程度上抑制松动区的发展。节理密度的影响权重为0.15，说明节理作为岩体中的不连续面，其密度越大，岩体的完整性越差，在开挖扰动下，越容易沿着节理面发生破坏和滑移，进而扩大松动区范围。在洞室几何参数中，洞室埋深的影响权重为0.18。随着洞室埋深的增加，上覆岩体的自重压力增大，洞室周边围岩所承受的初始应力也随之增大，使得围岩更容易发生破坏，导致松动区范围扩大。洞室跨度的影响权重为0.12，较大的洞室跨度会使围岩的受力更加复杂，应力集中现象更为明显，从而增加了围岩松动的可能性。通过对各因素影响权重的分析，可以明确在工程实践中，对于上软下硬节理岩体中地下洞室的稳定性控制，应重点关注软岩抗压强度、硬岩弹性模量、节理密度和洞室埋深等关键因素。在设计和施工过程中，可以采取相应的措施来改善这些因素，如对软岩进行加固处理，提高其抗压强度；合理选择洞室位置，尽量避开节理密集区域；优化洞室设计，控制洞室跨度和埋深等，以有效减小围岩松动区范围，保障地下洞室的稳定。四、地下洞室深浅埋分界标准研究4.1现有深浅埋分界方法概述在地下洞室工程领域，准确划分深浅埋对于合理设计支护结构、确保工程安全与经济具有关键意义。目前，常用的深浅埋分界方法主要基于埋深与洞径比、围岩压力分布等因素。基于埋深与洞径比的划分方法应用广泛。许多学者和工程规范采用这一指标进行初步界定。铁路隧道设计规范以统计隧道坍落拱高度为基础，对深浅埋进行划分。当埋深H≤h0时，用全土柱计算围岩压力；当h0＜H≤(2-2.5)h0时，采用谢家烋公式；当H＞(2-2.5)h0时，竖向压力取γh0。日本相关规定指出，如果覆土厚比隧道外径小（H＜D），采用全土柱压力；在砂土和硬粘土中，若覆土厚度比外径大许多倍（H＞＞D），则取“松弛土压”，按泰沙基公式计算。在实际工程中，部分设计遵循2倍洞室跨度为深浅埋分界，低于2倍洞室跨度时采用全土柱荷载，高于2倍洞室跨度时采用泰沙基公式。这种方法简单直观，易于操作，但存在一定局限性。它仅考虑了洞室的几何尺寸和埋深，忽略了岩体的力学性质、节理分布等重要地质因素。在不同地质条件下，相同的埋深与洞径比可能对应不同的围岩力学行为，导致分界结果不够准确。基于围岩压力分布的划分方法从力学本质出发，更能反映地下洞室的实际受力状态。其原理是根据洞室开挖后围岩压力的变化规律来确定深浅埋界限。在浅埋情况下，围岩压力主要由上覆岩土体的自重产生，可按全土柱压力计算。随着埋深增加，围岩的成拱效应逐渐显现，围岩压力不再随埋深线性增加，而是趋于稳定，此时进入深埋状态。比尔鲍曼公式在考虑围岩压力时，对粘性土，同时考虑粘聚力c和内摩擦角φ；对砂土，只计内摩擦角φ，不计粘聚力c。当H增加到一定程度，竖向应力σv趋于常数时，可判定为深埋，该方法无需明确的深浅埋分界。泰沙基公式则认为当H≥5B时，σv达到常数，可确定为深埋。这种方法的优点是考虑了围岩的力学响应，但在实际应用中，准确获取围岩压力分布较为困难，需要进行复杂的现场测试和数值模拟。不同的围岩条件和洞室形状会导致围岩压力分布的差异，使得该方法的通用性受到一定限制。4.2根据洞室松动区形态确定的深浅埋分界标准在地下洞室工程中，准确判断洞室所处的深浅埋状态至关重要，它直接关系到支护设计的合理性以及工程的安全性与经济性。传统的基于埋深与洞径比、围岩压力分布等的深浅埋分界方法存在一定局限性，难以全面考虑上软下硬节理岩体的复杂特性。本研究提出一种依据洞室松动区形态确定深浅埋分界的新方法，旨在更准确地反映洞室的实际受力和稳定状态。在浅埋情况下，由于上覆岩体厚度相对较薄，洞室开挖后，围岩的应力重分布范围相对较小，且受到地表的影响较大。此时，洞室周边的松动区形态通常呈现出较为规则的形状，多为近似圆形或椭圆形，且松动区范围相对较大，容易贯通至地表。这是因为浅埋洞室在开挖后，上覆岩体无法形成有效的承载拱，围岩主要依靠自身的强度来抵抗变形和破坏。在软岩部分，由于其强度较低，在开挖扰动下，容易产生较大的塑性变形和破裂，使得松动区迅速扩展。在某浅埋地下洞室工程中，通过现场钻孔摄像和声波测试发现，洞室拱顶的松动区直接与地表相连，松动区厚度达到了洞室直径的1.2倍。随着埋深的增加，洞室进入深埋状态。在深埋条件下，上覆岩体厚度较大，地应力也相应增大。洞室开挖后，围岩的应力重分布范围更广，深部岩体的约束作用更加明显。此时，洞室周边的松动区形态会发生显著变化，通常呈现出不规则的形状，且松动区范围相对较小，主要集中在洞室周边一定范围内，不会贯通至地表。深埋洞室在开挖后，上覆岩体能够形成较为稳定的承载拱，分担了部分围岩压力，使得洞室周边的岩体受力相对均匀。在硬岩部分，由于其强度较高，在承载拱的保护下，能够较好地抵抗变形和破坏，限制了松动区的发展。在某深埋隧道工程中，数值模拟结果显示，洞室周边的松动区主要集中在洞室半径1.5倍范围内，且松动区形状不规则，呈破碎块状分布。根据洞室松动区形态确定深浅埋分界标准具有重要的工程意义。这种方法能够更直观、准确地反映洞室所处的地质条件和受力状态，为支护设计提供更可靠的依据。在浅埋洞室中，由于松动区容易贯通至地表，对地表建筑物和环境的影响较大，因此需要采用更加强化的支护措施，如增加支护结构的刚度和强度，采用超前支护等，以控制围岩的变形和破坏。而在深埋洞室中，虽然松动区范围相对较小，但由于地应力较大，对支护结构的承载能力要求较高，需要选择合适的支护形式和参数，充分发挥围岩的自承能力，确保洞室的稳定。这种基于松动区形态的分界标准还可以为工程施工提供指导，根据不同的深浅埋状态，合理选择施工方法和施工顺序，优化施工工艺，提高施工效率和安全性。五、围岩松动压力计算方法研究5.1经典计算方法原理与局限性在地下洞室工程领域，准确计算围岩松动压力是确保工程安全与经济的关键环节。平衡拱理论和太沙基理论作为经典的围岩松动压力计算方法，在工程实践中曾被广泛应用，它们各自基于独特的假设和力学原理构建。平衡拱理论，由苏联学者M.M.普罗托耶科诺夫提出，又称普氏拱理论。该理论的基本假设建立在自然平衡拱概念之上。其认为，岩体经节理切割，开挖后虽呈松散状态，但仍具备一定粘结力。洞室开挖后，硐顶岩体将形成自然平衡拱，作用在硐顶的围岩压力仅是自然平衡拱内的岩体自重。在硐室的侧壁处，沿与侧壁夹角为45-φ/2（φ为岩体内摩擦角）的方向会产生两个滑动面。为表征岩体强度，普氏引入坚固系数f，其物理意义为f=c/σ+tanφ（c为岩体内聚力，σ为正应力），在实际应用中，通过经验公式f=Rc/10（Rc为单轴抗压强度，MPa）来方便地求得f值。在实际应用中，还需同时考虑岩体的完整性和地下水的影响。普氏认为形成的自然平衡拱的硐顶岩体只能承受压应力，不能承受拉应力。在计算自然平衡拱拱轴线方程时，先假设拱周线是一条二次曲线，在拱轴线上任取一点M（x，y），根据拱轴线不能承受拉力的条件，所有外力对M点的弯矩应为零，即qx²-2Ty=0（q为拱轴线上部岩体的自重所产生的均布荷载，T为平衡拱拱顶截面的水平推力）。由静力平衡方程可知，水平推力T与作用在拱脚的水平推力T'数值相等，方向相反。由于拱脚易产生水平位移改变拱的内力分布，普氏认为拱脚的水平推力T'必须满足T'≤qa1f（a1为自然平衡拱的最大跨度），为安全起见，又将水平推力降低一半后，令T=qa1f/2，代入可得拱轴线方程y=qx²/qa1f，这是一条抛物线。根据此式可求得拱轴线上任意一点的高度。当侧壁稳定时，x=a（a为侧壁稳定时平衡拱的跨度），y=b（b为拱的矢高，即自然平衡拱的最大高度），可得b=qa²/qa1f；当侧壁不稳定时，x=a1，y=b1，可得b1=qa1²/qa1f。硐顶最大围岩压力可按下式计算q=γb（γ为岩体容重），普氏围岩压力理论中的侧向压力可按下式计算e=γbtan²(45-φ/2)（侧壁顶部）和e=γ(b+h)tan²(45-φ/2)（侧壁底部，h为洞室高度）。太沙基理论则从应力传递原理出发推导竖向围岩压力。该理论假定岩体为散体，但具有一定的内聚力，且服从库仑准则τ=c+σntanφ（τ为岩体抗剪强度，c为岩体内聚力，σn为作用于剪切面上的正应力，φ为岩体内摩擦角）。以浅埋隧道为例，当支护结构受到上覆地压作用时，会发生挠曲变形并引起地块移动。假设围岩的内摩擦角为φ，滑移面从隧道底面以φ的角度倾斜，到硐顶后以适当的曲线到达地表面，实际计算中近似假定为两条垂直线。设从地表面到拱顶的滑动地块的宽度为2a（a为硐室半宽），其值等于a+Htanφ（H为开挖高度）。假定硐室顶壁衬砌顶部两端出现一直延伸到地表面的竖向破裂面。在滑动地块中，切取厚度为dz的薄层单元为分析对象，该薄层单元共受五种力作用：单元体自重dG=2aγdz（γ为岩体容重）、作用于单元体上表面的竖直向下的上覆岩体压力dP1=2aσdz（σ为竖向初始地应力）、作用于单元体下表面的竖直向上的下伏岩体托力dP2=2a(σ+dσ)dz、作用于单元体侧面的竖直向上的侧向围岩摩擦力dT=2τdz（τ为岩体抗剪强度）、作用于单元体侧面的水平方向的侧向围岩压力dE（可通过侧压力系数k0与竖向应力关联）。由薄层单元体在竖向的平衡条件dP1+dG-dP2-dT=0，将各力表达式代入并整理可得微分方程dσ/dz+(2k0τ)/a=γ。通过积分求解该微分方程，并结合边界条件（当z=0时，σ=p0，p0为地表面荷载），最终得到硐室顶部的竖向围岩压力q的表达式。对于埋深很大的深埋硐室，地表面的荷载对硐室顶部竖向围岩压力已不产生影响。作用在侧壁的围岩压力假设为一梯形，梯形上、下部的围岩压力可按相应公式计算。尽管这两种经典理论在围岩松动压力计算方面具有重要意义，但在面对上软下硬节理岩体这种复杂地质条件时，其局限性也十分显著。平衡拱理论将岩体视为散体，与上软下硬节理岩体的实际情况存在较大差异。在这种岩体中，软岩和硬岩的力学性质差异巨大，节理的存在又进一步增加了岩体的复杂性。而平衡拱理论难以准确反映这种非均质性和各向异性。该理论引入的坚固系数f并非岩体的特性参数，无法通过实验准确确定，且f值受多种因素影响，在实际应用中难以合理选取。根据该理论，硐室顶部中央围岩压力最大，但在许多上软下硬节理岩体的实际工程中，最大顶压常常偏离拱顶，这一现象是平衡拱理论难以解释的。该理论还表明，硐室围岩压力只与其跨度有关，而与断面形式、上覆岩层厚度，以及施工的方法、程度和进度等均无关，这与上软下硬节理岩体中地下洞室的实际情况不完全相符。太沙基理论同样存在局限性。它虽然考虑了岩体的内聚力和抗剪强度，但对于上软下硬节理岩体中节理的复杂力学行为，如节理的张开、闭合、剪切滑移以及节理面之间的相互作用等，缺乏足够的考虑。在这种岩体中，节理的存在会导致应力集中和应力分布的不均匀，而太沙基理论难以准确描述这种应力状态的变化。该理论在推导过程中进行了一些简化假设，如将滑移面简化为直线等，这些假设在复杂的上软下硬节理岩体中可能与实际情况偏差较大，从而影响计算结果的准确性。在实际应用中，准确获取太沙基理论所需的参数，如岩体的内聚力、内摩擦角等，对于上软下硬节理岩体来说也较为困难，因为不同部位的岩体参数可能存在较大差异。5.2修正经典计算方法鉴于经典计算方法在处理上软下硬节理岩体时存在的不足，本研究提出一种修正思路，旨在更准确地考虑节理影响、应力传递等关键因素，从而建立更贴合实际工程情况的围岩松动压力计算方法。在考虑节理影响方面，引入节理弱化系数来量化节理对岩体力学性质的削弱作用。节理弱化系数与节理的密度、粗糙度、充填物性质等因素密切相关。通过室内试验和数值模拟，建立节理弱化系数与这些因素之间的定量关系。对于节理密度较大、粗糙度较低且充填物为软弱材料的情况，节理弱化系数取值较小，以反映节理对岩体强度和刚度的显著降低作用。在计算围岩松动压力时，将节理弱化系数应用于岩体的力学参数，如弹性模量、泊松比、抗压强度等，使计算结果更能体现节理对岩体力学行为的影响。针对应力传递的复杂性，采用改进的应力传递模型。传统的应力传递模型在处理上软下硬节理岩体时，往往忽略了软岩和硬岩之间的界面效应以及节理对应力传播路径的改变。本研究建立的改进模型，充分考虑了软岩和硬岩的力学性质差异、节理的分布和方向，以及它们之间的相互作用。通过有限元分析等数值方法，模拟应力在岩体中的传播过程，确定应力集中区域和应力传递规律。在计算围岩松动压力时，根据应力传递模型得到的应力分布结果，准确计算作用在支护结构上的荷载。基于上述修正思路，推导适用于上软下硬节理岩体的围岩松动压力计算公式。以圆形洞室为例，假设洞室半径为R，埋深为H，上覆软岩厚度为h1，下卧硬岩厚度为h2。考虑节理影响和应力传递后，洞顶竖向围岩松动压力q的计算公式为：q=\gamma_1h_1+\gamma_2h_2+\Deltaq其中，\gamma_1和\gamma_2分别为软岩和硬岩的重度；\Deltaq为考虑节理和应力传递影响的附加压力，可通过节理弱化系数和改进的应力传递模型计算得到。洞室侧压力e的计算公式为：e=kq其中，k为侧压力系数，根据岩体的力学性质和节理分布情况确定。为验证修正计算方法的准确性和可靠性，选取某实际工程案例进行算例对比分析。该工程为一条穿越上软下硬节理岩体的公路隧道，隧道断面为马蹄形，跨度为10m，高度为8m，埋深为150m，上覆软岩厚度为30m，下卧硬岩厚度为120m。通过现场监测，得到隧道开挖后围岩松动压力的实测值。分别采用经典的平衡拱理论、太沙基理论以及本研究提出的修正计算方法对该隧道的围岩松动压力进行计算。计算结果表明，平衡拱理论计算得到的洞顶竖向围岩松动压力为1.2MPa，太沙基理论计算结果为1.5MPa，而修正计算方法得到的结果为1.35MPa。与实测值1.3MPa相比，修正计算方法的计算结果与实测值更为接近，相对误差仅为3.8%，而平衡拱理论和太沙基理论的相对误差分别为7.7%和15.4%。在洞室侧压力计算方面，修正计算方法同样表现出更好的准确性。通过算例对比验证，充分证明了本研究提出的修正计算方法在计算上软下硬节理岩体中地下洞室围岩松动压力时具有更高的精度和可靠性，能够为工程设计和施工提供更准确的依据。5.3基于应力传递机理的计算方法在深入研究上软下硬节理岩体中地下洞室围岩松动压力的过程中，基于应力传递机理构建计算模型具有重要意义。该模型的建立充分考虑了岩体的非均质性、节理的存在以及应力在岩体中的复杂传递路径，旨在更准确地反映围岩松动压力的形成机制和分布规律。在模型假设方面，考虑到上软下硬节理岩体的实际情况，做如下假设：将软岩和硬岩视为不同力学性质的连续介质，但由于节理的切割作用，岩体的连续性和完整性受到一定程度的破坏。节理被简化为具有一定力学参数的不连续面，能够传递法向应力和切向应力。假设洞室开挖后，围岩的变形和破坏主要发生在一定范围内，即松动区，松动区外的岩体对松动区内的岩体有一定的约束作用。基于上述假设，建立的计算模型主要考虑以下关键因素。在应力传递路径方面，洞室开挖后，上覆岩体的自重应力以及构造应力会重新分布。由于软岩和硬岩的弹性模量等力学参数不同，应力在两者中的传递速度和方式存在差异。软岩的弹性模量较小，应力在软岩中传递时更容易产生应力集中现象。节理的存在改变了应力的传递路径，当应力传播到节理面时，会根据节理面的力学性质发生反射、折射和绕射等现象。在节理面粗糙度较大的情况下，应力在节理面上的传递会受到较大的阻碍，导致节理面附近的应力集中。根据应力传递机理，推导不同条件下的计算公式。对于圆形洞室，假设洞室半径为R，埋深为H，上覆软岩厚度为h1，下卧硬岩厚度为h2。在不考虑节理影响时，根据弹性力学理论，洞室周边的切向应力σθ可表示为：\sigma_{\theta}=\frac{\gammaH}{2}\left(1+\frac{R^{2}}{r^{2}}\right)其中，γ为岩体的重度，r为计算点到洞室中心的距离。当考虑节理影响时，引入节理影响系数β，β与节理的密度、粗糙度等因素有关。此时，洞室周边的切向应力修正为：\sigma_{\theta}=\frac{\gammaH}{2}\left(1+\frac{R^{2}}{r^{2}}\right)\beta对于洞室顶部的竖向松动压力q，在考虑应力传递和节理影响后，可推导如下公式：q=\gamma_1h_1+\gamma_2h_2+\Delta\sigma其中，\gamma_1和\gamma_2分别为软岩和硬岩的重度；\Delta\sigma为由于应力集中和节理影响产生的附加应力，可通过对节理面的力学分析和应力传递路径的研究得到。通过具体算例来分析该计算方法的合理性和适用性。假设有一圆形隧道，半径为3m，埋深为100m，上覆软岩厚度为20m，下卧硬岩厚度为80m。软岩的重度为20kN/m³，弹性模量为1GPa；硬岩的重度为25kN/m³，弹性模量为10GPa。节理密度为每米5条，节理粗糙度系数为0.5。根据上述参数，利用本文提出的基于应力传递机理的计算方法，计算得到洞室顶部的竖向松动压力为3.5MPa。同时，采用有限元软件ANSYS对该隧道进行数值模拟，模拟结果显示洞室顶部的竖向松动压力为3.3MPa。将本文计算方法的结果与有限元模拟结果进行对比，相对误差为6.1%，在可接受的范围内。这表明本文提出的计算方法能够较为准确地计算上软下硬节理岩体中地下洞室围岩的松动压力，具有一定的合理性和适用性。在实际工程应用中，该计算方法能够为地下洞室的支护设计提供重要的参考依据，有助于提高工程的安全性和经济性。5.4基于极限分析上限原理的计算方法极限分析上限原理作为岩土力学领域的重要理论，在分析结构极限状态和承载能力方面具有独特优势。该原理基于虚功原理，通过假设结构在极限状态下的破坏模式，建立外功率与内功率的平衡方程，从而求解结构的极限荷载。在地下洞室围岩稳定性分析中，该原理能够有效考虑岩体的非线性力学行为和复杂的破坏机制，为围岩松动压力的计算提供了一种有效的途径。在应用极限分析上限原理计算围岩松动压力时，首先需要构建合理的力学模型。以圆形洞室为例，假设洞室半径为R，埋深为H，围岩为上软下硬节理岩体，软岩和硬岩的分界面水平，位于洞室下方深度h处。考虑到节理的存在，将岩体视为由岩块和节理组成的非连续介质，节理的力学性质通过节理单元来模拟。假设洞室开挖后，围岩发生破坏，形成一个倒锥形的破坏区域。在破坏区域内，岩块之间发生相对滑动和转动，节理面产生剪切滑移和拉伸开裂。根据极限分析上限原理，在极限状态下，外荷载所做的功率等于结构内部耗散的功率。外荷载主要包括上覆岩体的自重和可能存在的地面荷载。上覆岩体的自重功率为：P_{out}=\gammaV\dot{u}其中，\gamma为岩体的重度，V为破坏区域内岩体的体积，\dot{u}为破坏区域的竖向速度。结构内部耗散的功率包括岩块之间的摩擦耗能和节理面的耗能。岩块之间的摩擦耗能功率为：P_{friction}=\sum_{i=1}^{n}\tau_{i}A_{i}\dot{u}_{i}其中，\tau_{i}为第i个岩块界面的剪应力，A_{i}为第i个岩块界面的面积，\dot{u}_{i}为第i个岩块界面的相对滑动速度。节理面的耗能功率为：P_{joint}=\sum_{j=1}^{m}c_{j}A_{j}\dot{u}_{j}+\sum_{j=1}^{m}\sigma_{nj}A_{j}\tan\varphi_{j}\dot{u}_{j}其中，c_{j}为第j个节理面的黏聚力，A_{j}为第j个节理面的面积，\sigma_{nj}为第j个节理面的法向应力，\varphi_{j}为第j个节理面的内摩擦角，\dot{u}_{j}为第j个节理面的相对滑动速度。根据外功率等于内功率的平衡条件，可得：P_{out}=P_{friction}+P_{joint}通过求解上述方程，即可得到围岩松动压力。为验证基于极限分析上限原理的计算方法的准确性和优势，选取某实际工程进行算例分析。该工程为一条穿越上软下硬节理岩体的铁路隧道，隧道直径为8m，埋深为180m，上覆软岩厚度为40m，下卧硬岩厚度为140m。通过现场监测，得到隧道开挖后围岩松动压力的实测值。分别采用基于极限分析上限原理的计算方法、经典的平衡拱理论和太沙基理论对该隧道的围岩松动压力进行计算。计算结果表明，基于极限分析上限原理的计算方法得到的洞顶竖向围岩松动压力为2.8MPa，平衡拱理论计算结果为2.3MPa，太沙基理论计算结果为2.5MPa。与实测值2.7MPa相比，基于极限分析上限原理的计算方法的计算结果与实测值更为接近，相对误差仅为3.7%，而平衡拱理论和太沙基理论的相对误差分别为14.8%和7.4%。通过算例分析可知，基于极限分析上限原理的计算方法能够更准确地考虑上软下硬节理岩体的复杂力学行为和破坏机制，计算结果更接近实际情况。该方法在计算围岩松动压力时具有较高的准确性和可靠性，能够为地下洞室的支护设计提供更合理的依据。六、工程实例分析6.1工程概况某公路隧道工程穿越上软下硬节理岩体，具有典型的工程特性。该隧道位于山区，线路全长2500m，其中约800m段落处于上软下硬节理岩体地层。从地质条件来看，上部软岩主要为强风化页岩，呈灰黑色，节理裂隙极为发育，岩石破碎，手捏易碎。经室内试验测定，其饱和单轴抗压强度平均值仅为5MPa，弹性模量约为0.3GPa，泊松比为0.38。页岩的抗剪强度较低，黏聚力为12kPa，内摩擦角为18°。下部硬岩为中风化花岗岩，呈灰白色，结构致密，岩芯多呈长柱状，完整性较好。其饱和单轴抗压强度高达80MPa，弹性模量达到12GPa，泊松比为0.26。花岗岩的抗剪强度较高，黏聚力为450kPa，内摩擦角为35°。在节理分布方面，软岩和硬岩中均存在节理，但软岩中的节理更为密集，平均每米节理条数达到15条，节理间距较小，多在0.1-0.3m之间。硬岩中的节理相对稀疏，平均每米节理条数为5条，节理间距多在0.5-1.0m之间。节理的产状较为复杂，软岩中节理倾向主要为NE30°-50°，倾角在40°-60°之间；硬岩中节理倾向主要为NW290°-310°，倾角在30°-50°之间。部分节理面较为粗糙，起伏差较大，且软岩节理多被黏土充填，硬岩节理部分被方解石脉充填。隧道设计为双向四车道，采用三心圆曲墙式断面，净宽10.5m，净高7.0m。衬砌结构采用复合式衬砌，初期支护采用喷射混凝土和锚杆联合支护，喷射混凝土厚度为25cm，锚杆长度为3.5m，间距为1.0m×1.0m。二次衬砌采用钢筋混凝土，厚度为50cm。施工过程中，采用台阶法开挖。上台阶超前下台阶3-5m，每循环进尺控制在1.0-1.5m。在开挖过程中，密切关注围岩的变形和稳定情况，及时根据地质条件的变化调整施工参数。当遇到软岩段落时，适当减小进尺，加强初期支护；当穿越软硬岩交界部位时，采取加强支护措施，如增加锚杆数量、喷射混凝土厚度等，以确保施工安全和隧道稳定。6.2各方法松动压力计算结果综合对比及差异性分析为深入探究不同计算方法在处理上软下硬节理岩体中地下洞室围岩松动压力时的性能差异，本研究针对某公路隧道工程实例，分别运用平衡拱理论、太沙基理论、本文修正计算方法以及基于应力传递机理和极限分析上限原理的计算方法，对围岩松动压力进行了详细计算，并对结果展开全面对比与深入分析。平衡拱理论计算得到洞顶竖向围岩松动压力为1.8MPa，太沙基理论计算结果为2.2MPa。本文修正计算方法充分考虑节理影响和应力传递的复杂性，计算结果为2.0MPa。基于应力传递机理的计算方法通过对岩体应力传递路径的细致分析，得到洞顶竖向围岩松动压力为2.1MPa。基于极限分析上限原理的计算方法，考虑了岩体的非线性力学行为和复杂破坏机制，计算结果为1.9MPa。从计算结果可以看出，各方法的计算结果存在一定差异。平衡拱理论计算结果相对较低，这是因为该理论将岩体视为散体，对节理和岩体非均质性考虑不足，导致计算结果与实际情况存在偏差。太沙基理论虽然考虑了岩体的内聚力和抗剪强度，但对节理的复杂力学行为考虑不够全面，计算结果相对较高。本文修正计算方法、基于应力传递机理和极限分析上限原理的计算方法，由于充分考虑了上软下硬节理岩体的特性，计算结果更为接近实际情况。在洞室侧压力计算方面，各方法也存在差异。平衡拱理论和太沙基理论的计算结果与实际情况偏差较大，而本文修正计算方法和基于应力传递机理、极限分析上限原理的计算方法，能够更好地反映上软下硬节理岩体中洞室侧压力的实际分布情况。通过对各方法计算结果的对比分析可知，在处理上软下硬节理岩体中地下洞室围岩松动压力时，传统的平衡拱理论和太沙基理论存在一定局限性。而本文提出的修正计算方法以及基于应力传递机理和极限分析上限原理的计算方法，能够更准确地考虑岩体的非均质性、节理的影响以及复杂的力学行为，计算结果更可靠，更适用于此类复杂地质条件下的工程实践。在实际工程中，应根据具体的地质条件和工程要求，合理选择计算方法，以确保地下洞室的稳定性和安全性。6.3基于松动区边界的应力传递法合理性证明为验证基于松动区边界的应力传递法在计算上软下硬节理岩体中地下洞室围岩松动压力的合理性，本研究将其应用于某公路隧道工程实例，并结合现场监测数据和数值模拟结果进行分析。在现场监测方面，在隧道洞身布置了多个监测断面，每个监测断面设置多个压力盒，分别用于测量洞顶、拱腰和边墙等部位的围岩压力。同时，采用全站仪对隧道周边的位移进行实时监测，以获取围岩的变形情况。在施工过程中，密切关注监测数据的变化，及时记录不同施工阶段的围岩压力和位移数据。数值模拟选用有限元软件ANSYS建立隧道模型，模型中精确模拟上软下硬节理岩体的力学特性。软岩和硬岩分别采用不同的本构模型，节理通过接触单元进行模拟，充分考虑节理的法向和切向力学行为。模型的边界条件和荷载施加根据实际工程情况进行设定，确保模拟结果的真实性。通过数值模拟，得到隧道开挖过程中围岩的应力、应变分布以及围岩压力的大小。将基于松动区边界的应力传递法计算结果与现场监测数据和数值模拟结果进行对比分析。在洞顶竖向围岩松动压力方面，应力传递法计算结果为2.05MPa，现场监测平均值为2.0MPa，数值模拟结果为2.1MPa。应力传递法计算结果与现场监测值的相对误差为2.5%，与数值模拟结果的相对误差为2.4%。在洞室侧压力方面，应力传递法计算结果与现场监测和数值模拟结果也具有较好的一致性。从对比结果可以看出，基于松动区边界的应力传递法计算结果与现场监测数据和数值模拟结果吻合较好，能够准确地反映上软下硬节理岩体中地下洞室围岩的松动压力。该方法充分考虑了岩体的非均质性、节理的影响以及应力在岩体中的传递特性，具有较高的合理性和可靠性。在实际工程应用中，基于松动区边界的应力传递法能够为地下洞室的支护设计提供准确的依据，有助于保障工程的安全和稳定。6.4节理岩体中地下洞室松动压力推荐计算方法综合考虑各计算方法的优缺点以及工程实际应用需求，本研究推荐基于应力传递机理和极限分析上限原理的计算方法作为上软下硬节理岩体中地下洞室松动压力的主要计算方法。这两种方法能够较为全面地考虑上软下硬节理岩体的非均质性、节理的影响以及复杂的力学行为，计算结果相对准确可靠。在实际应用中，应首先对工程场地进行详细的地质勘察，获取准确的岩体力学参数、节理分布特征等数据。通过现场测试和室内试验，确定软岩和硬岩的弹性模量、泊松比、抗压强度、抗剪强度等参数，以及节理的法向刚度、切向刚度、黏聚力、内摩擦角等参数。这些参数是计算方法的基础，其准确性直接影响计算结果的可靠性。根据地质勘察结果，结合工程实际情况，选择合适的计算方法。如果岩体的节理分布较为规则，且应力传递路径相对清晰，基于应力传递机理的计算方法能够较好地发挥作用。该方法通过分析应力在岩体中的传递过程，考虑软岩和硬岩之间的力学差异以及节理对应力的影响，能够准确计算围岩松动压力。在某工程中，岩体的节理主要呈平行分布，且软岩和硬岩的分界面较为平整，采用基于应力传递机理的计算方法，计算结果与现场监测数据吻合较好。若岩体的破坏机制较为复杂，节理的相互作用显著，极限分析上限原理的计算方法则更为适用。该方法通过假设破坏模式，建立外功率与内功率的平衡方程，能够充分考虑岩体的非线性力学行为和节理的复杂破坏机制。在节理密集发育且相互交错的岩体中，极限分析上限原理的计算方法能够更准确地评估围岩松动压力。为进一步提高计算结果的准确性，可将两种方法的计算结果进行对比分析，并结合工程经验和现场监测数据进行综合判断。如果两种方法的计算结果相近，且与现场监测数据相符，则可以认为计算结果可靠。若计算结果存在差异，应进一步分析原因，检查计算参数的合理性和计算方法的适用性。在某工程中，基于应力传递机理的计算方法得到的洞顶竖向围岩松动压力为2.0MPa，极限分析上限原理的计算方法得到的结果为2.1MPa，现场监测值为2.05MPa，两种方法的计算结果与现场监测值较为接近，验证了计算方法的准确性。在实际工程应用中，还应考虑其他因素对围岩松动压力的影响，如地下水的作用、施工过程中的扰动等。地下水会降低岩体的强度和抗剪能力，增加围岩松动压力。施工过程中的爆破、开挖等作业也会对围岩产生扰动，导致围岩松动压力的变化。因此，在计算围岩松动压力时，应充分考虑这些因素的影响，采取相应的措施进行修正和调整。在地下水丰富的地区，可通过增加排水措施，降低地下水位，减少地下水对围岩的影响。在施工过程中，应优化施工工艺，减少对围岩的扰动，降低围岩松动压力。七、结论与展望7.1主要研究结论本研究围绕上软下硬节理岩体中地下洞室围岩松动特性与松动压力计算方法展开，通过理论分析、数值模拟和工程实例验证，取得了以下主要研究成果：上软下硬节理岩体特性：通过室内试验和数值模拟，深入研究了上软下硬节理岩体的力学特性。明确了软岩和硬岩的力学参数差