基于多轴强度安全系数的水电站地下厂房围岩稳定性深度剖析与精准评估

基于多轴强度安全系数的水电站地下厂房围岩稳定性深度剖析与精准评估一、引言1.1研究背景与意义在全球积极推进能源结构转型的大背景下，清洁能源的开发与利用成为关键议题。水电站作为一种重要的清洁能源形式，凭借其可再生、低碳排放等显著优势，在能源领域中占据着愈发重要的地位。近年来，中国清洁低碳化进程不断加快，水电、风电、光伏、在建核电装机规模等多项指标保持世界第一，建成世界最大清洁发电体系，成为推动全球清洁能源发展的重要力量。乌东德水电站创造8项“世界第一”，拿下15项“全球首次”，攻克了一系列世界级难题，年均发电389.1亿千瓦时，每年可节约标准煤1220万吨，减少二氧化碳排放3050万吨。华电金上昌波水电站建成后预计年发电量将达到43.55亿千瓦时，相当于节约130万吨标准煤，并且每年减少超过359万吨的二氧化碳排放量。这些大型水电站的成功建设与运营，充分彰显了水电在清洁能源领域的重要价值。水电站的建设是一项复杂而庞大的系统工程，其中地下厂房的建设更是面临着诸多挑战。地下厂房的建设需要深入了解地质环境和岩石力学参数等方面的数据，因为地质构造、地下水位等因素会对厂房的稳定性产生重大影响。围岩稳定性分析作为地下工程设计中的一项核心内容，是保障地下工程安全性和稳定性的关键技术环节。一旦围岩出现失稳现象，可能引发坍塌、渗漏等严重事故，不仅会对工程本身造成巨大的经济损失，还可能危及施工人员的生命安全，甚至对周边环境产生长期的负面影响。在围岩稳定性分析的众多方法和指标中，多轴强度安全系数作为评估地下工程围岩稳定性的重要指标之一，具有独特的优势和价值。它能够全面评估围岩破坏的判据，充分考虑到岩石在多轴应力状态下的力学行为，为地下工程的设计和施工提供更为科学、准确的参考依据。传统的安全系数方法，如古典安全系数、超载系数及强度折减系数等，在一般多参数强度准则下存在一定的局限性，难以全面、准确地反映围岩的实际受力状态和稳定性情况。而多轴强度安全系数能够突破这些传统方法的局限，从多个维度对围岩的稳定性进行综合评估，从而为工程设计和施工提供更具针对性和可靠性的指导。综上所述，利用多轴强度安全系数方法对水电站地下厂房的围岩稳定性进行深入分析，具有重要的现实意义和理论价值。这不仅有助于提高水电站地下工程设计的可靠性和稳定性，保障工程的安全建设与运营，还能为清洁能源的可持续发展提供坚实的技术支撑，推动能源结构的优化升级，为实现碳达峰碳中和的目标做出积极贡献。1.2国内外研究现状围岩稳定性分析作为地下工程领域的关键研究方向，一直受到国内外学者的广泛关注。在国外，早期的研究主要集中在理论分析和经验公式的推导上。Terzaghi在20世纪中期提出了经典的太沙基理论，为地下工程围岩稳定性分析奠定了理论基础。此后，随着计算机技术的飞速发展，数值模拟方法逐渐成为研究围岩稳定性的重要手段。有限元法（FEM）、边界元法（BEM）、离散元法（DEM）等数值方法被广泛应用于围岩稳定性分析中。例如，Cundall和Strack提出的离散元法，能够较好地模拟岩石的非连续变形和破坏过程，为研究复杂地质条件下的围岩稳定性提供了有力工具。在国内，围岩稳定性分析的研究起步相对较晚，但发展迅速。20世纪80年代以来，随着大量地下工程的建设，国内学者在围岩稳定性分析方面取得了丰硕的成果。周维垣等学者通过对大量工程实例的分析，提出了适合我国地质条件的围岩稳定性评价方法。近年来，随着深部地下工程的不断涌现，如深埋隧道、深部矿山开采等，对围岩稳定性分析提出了更高的要求。学者们开始关注深部岩体的力学特性和破坏机理，研究高地应力、高温、地下水等复杂环境因素对围岩稳定性的影响。多轴强度安全系数作为评估地下工程围岩稳定性的重要指标，其应用研究也逐渐成为国内外学者关注的焦点。国外学者在多轴强度理论和安全系数计算方法方面开展了大量的研究工作。如Drucker和Prager提出的Drucker-Prager准则，考虑了中间主应力对岩石强度的影响，为多轴强度安全系数的计算提供了理论基础。随后，学者们在此基础上不断改进和完善，提出了多种多轴强度准则和安全系数计算方法。在国内，多轴强度安全系数的应用研究也取得了一定的进展。郑微微采用安全距离系数定义方法，对水电站地下厂房的围岩稳定性进行了分析，验证了该方法在解决实际工程问题中的可行性和有效性。但目前多轴强度安全系数在水电站地下厂房围岩稳定性分析中的应用还存在一些不足之处。一方面，不同的多轴强度准则和安全系数计算方法之间存在差异，导致计算结果的可靠性和可比性有待提高；另一方面，多轴强度安全系数与其他稳定性评价指标之间的关系研究还不够深入，缺乏系统的综合评价方法。综上所述，目前国内外在水电站地下厂房围岩稳定性分析以及多轴强度安全系数应用方面已经取得了一定的研究成果，但仍存在一些问题和不足。本研究将在现有研究的基础上，进一步深入探究多轴强度安全系数在水电站地下厂房围岩稳定性分析中的应用，旨在提出更加科学、准确的围岩稳定性评价方法，为水电站地下工程的设计和施工提供更为可靠的理论依据和技术支持。1.3研究目标与内容本研究旨在利用多轴强度安全系数方法，对水电站地下厂房的围岩稳定性进行深入分析，为水电站地下工程的设计和施工提供科学依据。具体研究内容如下：岩石力学参数测定方法研究：岩石力学参数是分析围岩稳定性的基础数据，其准确性直接影响到分析结果的可靠性。因此，本研究将深入研究水电站地下厂房岩石力学参数的测定方法，包括岩石抗压强度、岩石剪切强度、岩石弹性模量等参数的测定。计划采用标准试验方法，如单轴和三轴压缩试验、剪切试验等，对水电站地下厂房周边围岩进行力学参数测定。通过这些试验，获取岩石在不同受力状态下的力学性能数据，为后续的围岩稳定性分析提供准确的参数支持。地质条件调查研究：地质条件是影响水电站地下厂房围岩稳定性的重要因素之一。本研究将对水电站地下厂房的地质条件进行全面调查研究，包括野外地质调查、岩芯取样等工作。通过野外地质调查，详细了解地下厂房周边的地质构造、地层分布、岩石类型等信息；通过岩芯取样，获取地下水位、围岩厚度和岩性等参数。这些地质数据将为建立准确的围岩稳定性数学模型提供重要依据。围岩稳定性数学模型建立：在获取岩石力学参数和地质条件数据的基础上，本研究将建立水电站地下厂房围岩稳定性的数学模型。以多轴强度安全系数为评价指标，充分考虑地质条件和岩石力学参数等因素，运用数学方法对围岩稳定性进行定量分析。计划采用MATLAB等数值计算软件进行模型构建，通过数学模型的建立，能够更加直观地反映围岩在不同工况下的稳定性状态，为工程设计和施工提供科学的决策依据。有限元分析：采用有限元方法对建立的数学模型进行仿真分析，是本研究的重要环节之一。有限元分析能够模拟围岩在实际受力情况下的变形和破坏过程，验证理论计算结果的准确性。本研究将利用ANSYS等有限元分析软件对建立的数学模型进行仿真分析，详细分析围岩的变形和破坏过程以及变形和破坏的范围等情况。通过有限元分析，可以直观地看到围岩在不同开挖阶段和不同荷载作用下的应力应变分布情况，为评估围岩稳定性提供可靠的数值模拟依据。灵敏度分析：围岩稳定性受到多种因素的影响，不同因素对围岩稳定性的影响程度各不相同。为了明确各因素对围岩稳定性的影响程度，本研究将对围岩稳定性的影响因素进行灵敏度分析。通过改变模型中的相关参数，观察多轴强度安全系数的变化情况，从而确定各因素对围岩稳定性的敏感程度。根据灵敏度分析结果，有针对性地提出相应的措施和建议，以提高水电站地下工程的设计和施工水平。例如，对于对围岩稳定性影响较大的因素，在工程设计和施工中应重点关注和控制，采取相应的加固措施或优化设计方案，以确保围岩的稳定性。1.4研究方法与技术路线本研究综合运用多种研究方法，从岩石力学参数测定、地质条件调查、数学模型建立、有限元分析到灵敏度分析，形成了一套系统、全面的研究体系，确保能够深入、准确地分析水电站地下厂房围岩的稳定性。标准试验方法：在岩石力学参数测定方面，采用标准试验方法，对水电站地下厂房周边围岩进行力学参数测定。单轴压缩试验是将岩石试件加工成标准尺寸，在压力试验机上施加轴向压力，直至试件破坏，通过测量破坏时的荷载和试件的尺寸，计算出岩石的抗压强度。三轴压缩试验则是在三个方向上对岩石试件施加不同的压力，模拟岩石在实际工程中的多轴受力状态，从而获取岩石在复杂应力条件下的力学参数，包括抗压强度、抗剪强度等。剪切试验通过对岩石试件施加剪切力，测量其抗剪强度，这些试验方法能够准确获取岩石的力学性能数据，为后续的分析提供可靠的基础。地质调查方法：对于地质条件的调查研究，采用钻孔、地质横断面等方法，对水电站地下厂房周边的地质条件进行全面、深入的调查。通过钻孔获取岩芯样本，对岩芯进行详细的分析，了解岩石的岩性、结构、构造以及地下水的分布情况等。地质横断面则是通过在不同位置进行地质测量和观察，绘制出地质剖面图，直观展示地下厂房周边的地质构造、地层分布等信息，为建立准确的围岩稳定性数学模型提供重要的地质依据。数值计算方法：在建立数学模型时，采用MATLAB等数值计算软件。MATLAB具有强大的数值计算和数据分析功能，能够方便地实现复杂的数学算法和模型构建。在建立水电站地下厂房围岩稳定性数学模型时，利用MATLAB的矩阵运算、函数拟合等功能，将地质条件、岩石力学参数等因素转化为数学表达式，以多轴强度安全系数为评价指标，建立起能够准确描述围岩稳定性的数学模型，为后续的分析和计算提供了有力的工具。有限元分析方法：利用ANSYS等有限元分析软件对建立的数学模型进行仿真分析。ANSYS是一款广泛应用于工程领域的有限元分析软件，具有强大的建模、求解和后处理功能。在本研究中，首先根据地下厂房的实际尺寸和地质条件，在ANSYS中建立三维有限元模型，划分网格，定义材料属性和边界条件。然后，通过加载不同的荷载工况，模拟围岩在实际受力情况下的变形和破坏过程，得到围岩的应力、应变分布云图，分析围岩的变形和破坏范围，验证理论计算结果的准确性，为评估围岩稳定性提供直观、可靠的数值模拟依据。灵敏度分析方法：对围岩稳定性的影响因素进行灵敏度分析时，通过改变模型中的相关参数，如岩石力学参数、地质构造参数等，观察多轴强度安全系数的变化情况。利用MATLAB的参数化分析功能，编写相应的程序，自动实现参数的改变和结果的计算。通过对计算结果的统计和分析，确定各因素对围岩稳定性的敏感程度，从而为工程设计和施工提供有针对性的建议，如对于敏感因素，在设计和施工中应重点关注和控制，采取相应的加固措施或优化设计方案，以提高围岩的稳定性。具体的技术路线如下：首先，进行岩石力学参数测定和地质条件调查，获取相关数据；接着，利用这些数据在MATLAB中建立围岩稳定性数学模型；然后，将数学模型导入ANSYS进行有限元分析，得到围岩的变形和破坏情况；最后，对影响因素进行灵敏度分析，根据分析结果提出相应的措施和建议。整个技术路线环环相扣，从数据获取到模型建立，再到分析和结果应用，确保了研究的科学性和可靠性，为水电站地下厂房围岩稳定性分析提供了全面、系统的研究方法。二、多轴强度安全系数理论基础2.1多轴强度安全系数定义及原理多轴强度安全系数是基于失效距离的概念来定义的。在多轴应力状态下，材料的强度失效面可以用一个数学函数来表示，该函数反映了材料在不同应力组合下的失效条件。失效距离则是指应力点到强度失效面的最短距离，它反映了材料在当前应力状态下距离失效的程度。当失效距离为零时，应力点恰好位于强度失效面上，此时材料处于极限状态；当失效距离大于零时，应力点位于强度失效面内部，材料处于安全状态；当失效距离小于零时，应力点位于强度失效面外部，材料已经失效。基于失效距离的多轴强度安全系数K定义为：K=\frac{\sigma_{0}}{\sigma_{e}}=\frac{\sigma_{0}}{\sqrt{\frac{1}{2}[(\sigma_{1}-\sigma_{2})^{2}+(\sigma_{2}-\sigma_{3})^{2}+(\sigma_{3}-\sigma_{1})^{2}]}}其中，\sigma_{0}为材料的单轴强度指标，如单轴抗压强度或单轴抗拉强度；\sigma_{e}为等效应力，它综合考虑了三个主应力\sigma_{1}、\sigma_{2}、\sigma_{3}的影响。等效应力的计算采用了一种能够反映多轴应力状态下材料受力复杂程度的方式，通过对主应力差值的平方和开方运算，得到一个能够衡量当前应力状态与单轴强度状态差异的指标。当等效应力小于单轴强度指标时，安全系数大于1，表明材料处于安全状态；反之，当等效应力大于单轴强度指标时，安全系数小于1，意味着材料可能发生失效。在不同的应力状态下，多轴强度安全系数的计算原理具有一致性，但具体的计算过程会因应力组合的不同而有所差异。在二轴拉压应力状态下，需要考虑两个主应力的大小和方向对材料强度的影响；在三轴受压应力状态下，三个主应力的相互作用更加复杂，对材料强度的影响也更为显著。不同的多参数强度准则也会对多轴强度安全系数的计算产生影响。Mohr-Coulomb准则考虑了材料的抗剪强度和正应力的影响，其表达式为：\tau=c+\sigma\tan\varphi其中，\tau为抗剪强度，c为黏聚力，\sigma为正应力，\varphi为内摩擦角。在使用该准则计算多轴强度安全系数时，需要根据具体的应力状态，将主应力代入准则中，计算出相应的强度值，再与实际应力进行比较，从而得到安全系数。Drucker-Prager准则则考虑了中间主应力对材料强度的影响，其表达式为：F=\alphaI_{1}+\sqrt{J_{2}}-k=0其中，I_{1}为第一应力不变量，J_{2}为第二偏应力不变量，\alpha和k为材料常数。在利用该准则计算多轴强度安全系数时，需要先计算出应力不变量，再根据准则方程求解出材料的强度值，进而计算安全系数。不同强度准则的适用范围和特点各不相同，因此在实际应用中，需要根据材料的特性和应力状态，选择合适的强度准则来计算多轴强度安全系数，以确保计算结果的准确性和可靠性。2.2与传统安全系数对比分析在围岩稳定性分析中，传统的安全系数方法包括古典安全系数、超载系数和强度折减系数，它们在工程实践中曾发挥了重要作用，但在一般多参数强度准则下存在一定的局限性。与这些传统安全系数相比，多轴强度安全系数具有独特的优势，能够更全面、准确地评估围岩的稳定性。古典安全系数基于古典强度理论，如最大拉应力理论、最大伸长线应变理论、最大切应力理论和畸变能密度理论等。以最大拉应力理论（第一强度理论）为例，其安全系数K_{c1}定义为：K_{c1}=\frac{\sigma_{t}}{\sigma_{1}}其中，\sigma_{t}为材料的抗拉强度，\sigma_{1}为最大主应力。古典安全系数的局限性在于它只考虑了单一的应力因素，如最大拉应力或最大切应力，而没有考虑材料在多轴应力状态下的复杂力学行为。在水电站地下厂房围岩的实际受力情况中，岩石往往处于多轴应力状态，单纯依据古典安全系数可能会过高或过低地估计围岩的稳定性。当围岩受到三向压应力作用时，古典安全系数可能无法准确反映岩石在这种复杂应力状态下的强度特性，导致对围岩稳定性的评估出现偏差。超载系数是指在不改变材料强度参数的情况下，通过逐步增加外荷载，直到结构达到极限状态时，所施加的荷载与初始荷载的比值。假设初始荷载为P_{0}，极限状态下的荷载为P_{u}，则超载系数K_{s}可表示为：K_{s}=\frac{P_{u}}{P_{0}}超载系数在一定程度上考虑了结构的承载能力储备，但它也存在局限性。在多参数强度准则下，超载系数的计算结果可能会受到加载路径和应力状态的影响。不同的加载路径可能导致结构在不同的应力组合下达到极限状态，从而使超载系数的计算结果缺乏唯一性和可比性。在复杂的地质条件下，水电站地下厂房围岩所受的地应力分布不均匀，加载路径复杂多变，这使得超载系数难以准确反映围岩的真实稳定性。强度折减系数是在不改变外荷载的情况下，通过逐步降低材料的强度参数，如凝聚力c和内摩擦角\varphi，直到结构达到极限状态时，初始强度参数与折减后强度参数的比值。以Mohr-Coulomb准则为例，假设初始凝聚力为c_{0}，初始内摩擦角为\varphi_{0}，折减后的凝聚力为c_{r}，折减后的内摩擦角为\varphi_{r}，则强度折减系数K_{r}可表示为：K_{r}=\frac{c_{0}}{c_{r}}=\frac{\tan\varphi_{0}}{\tan\varphi_{r}}强度折减系数虽然考虑了材料强度的变化对结构稳定性的影响，但它也存在一些问题。对于具有多个强度参数的复杂强度准则，很难确定哪些参数应该进行折减以及如何折减，这使得强度折减系数的计算过程存在一定的主观性和不确定性。在一些复杂的岩石力学模型中，强度参数之间存在相互关联，简单地对某个参数进行折减可能会导致计算结果与实际情况不符。相比之下，多轴强度安全系数具有明显的优势。多轴强度安全系数考虑了材料在多轴应力状态下的强度特性，能够更全面地反映围岩的受力情况。它基于失效距离的概念，通过计算应力点到强度失效面的距离来评估围岩的稳定性，这种方法更加科学、准确。在复杂的地质条件下，多轴强度安全系数能够综合考虑地应力、岩石力学参数等多种因素，对水电站地下厂房围岩的稳定性进行更合理的评价。多轴强度安全系数的计算方法相对简单，不需要进行复杂的加载路径模拟或强度参数折减。它只需要根据材料的强度准则和当前的应力状态，即可计算出安全系数，大大提高了计算效率和准确性。在实际工程应用中，多轴强度安全系数能够为水电站地下厂房的设计和施工提供更可靠的依据，有助于保障工程的安全和稳定。2.3多轴强度安全系数计算方法多轴强度安全系数的计算依赖于所选用的强度准则，不同的强度准则反映了材料在多轴应力状态下不同的破坏机理和强度特性。以混凝土材料为例，选取Mohr-Coulomb准则、Drucker-Prager准则、Ottosen准则、Hsieh-Ting-Chen准则和过镇海五参数准则这5种常用强度准则，详细介绍多轴强度安全系数的计算公式和计算过程。2.3.1Mohr-Coulomb准则下的计算Mohr-Coulomb准则是一种经典的强度准则，广泛应用于岩土和混凝土等材料的强度分析中。该准则认为材料的破坏主要取决于剪切应力，当材料某点的剪应力达到一定值时，材料就会发生破坏。其表达式为：\tau=c+\sigma\tan\varphi其中，\tau为抗剪强度，c为黏聚力，\sigma为正应力，\varphi为内摩擦角。在主应力空间中，Mohr-Coulomb准则的表达式可以转化为：F=\sigma_{1}-\sigma_{3}\frac{1+\sin\varphi}{1-\sin\varphi}-2c\frac{\cos\varphi}{1-\sin\varphi}=0在Mohr-Coulomb准则下，多轴强度安全系数K_{MC}的计算基于应力点到强度失效面的距离。对于给定的应力状态(\sigma_{1},\sigma_{2},\sigma_{3})，首先计算出该应力状态下的剪应力\tau和正应力\sigma：\tau=\frac{\sigma_{1}-\sigma_{3}}{2}\sigma=\frac{\sigma_{1}+\sigma_{3}}{2}然后将计算得到的\tau和\sigma代入Mohr-Coulomb准则表达式中，得到：F=\frac{\sigma_{1}-\sigma_{3}}{2}-c-\frac{\sigma_{1}+\sigma_{3}}{2}\tan\varphi多轴强度安全系数K_{MC}定义为：K_{MC}=\frac{\sigma_{0}}{\sigma_{e}}其中，\sigma_{0}为材料的单轴强度指标，如单轴抗压强度；\sigma_{e}为等效应力，在Mohr-Coulomb准则下，等效应力\sigma_{e}可以通过以下公式计算：\sigma_{e}=\frac{\sigma_{1}-\sigma_{3}}{1-\sin\varphi}假设混凝土材料的黏聚力c=1.5MPa，内摩擦角\varphi=30^{\circ}，单轴抗压强度\sigma_{0}=20MPa。对于某一应力状态(\sigma_{1}=10MPa,\sigma_{2}=5MPa,\sigma_{3}=3MPa)，首先计算剪应力\tau=\frac{10-3}{2}=3.5MPa，正应力\sigma=\frac{10+3}{2}=6.5MPa。代入Mohr-Coulomb准则表达式F=3.5-1.5-6.5\times\tan30^{\circ}\approx-0.27MPa。计算等效应力\sigma_{e}=\frac{10-3}{1-\sin30^{\circ}}=14MPa。则多轴强度安全系数K_{MC}=\frac{20}{14}\approx1.43。当安全系数K_{MC}\gt1时，表明该应力状态下混凝土材料处于安全状态；当K_{MC}\lt1时，材料可能发生破坏。2.3.2Drucker-Prager准则下的计算Drucker-Prager准则考虑了中间主应力对材料强度的影响，是在Mohr-Coulomb准则的基础上发展而来的。该准则在岩土和混凝土等材料的多轴强度分析中也有广泛应用，尤其适用于描述材料在复杂应力状态下的强度特性。其表达式为：F=\alphaI_{1}+\sqrt{J_{2}}-k=0其中，I_{1}为第一应力不变量，I_{1}=\sigma_{1}+\sigma_{2}+\sigma_{3}；J_{2}为第二偏应力不变量，J_{2}=\frac{1}{6}[(\sigma_{1}-\sigma_{2})^{2}+(\sigma_{2}-\sigma_{3})^{2}+(\sigma_{3}-\sigma_{1})^{2}]；\alpha和k为材料常数，与材料的性质有关。在Drucker-Prager准则下，多轴强度安全系数K_{DP}的计算同样基于应力点到强度失效面的距离。对于给定的应力状态(\sigma_{1},\sigma_{2},\sigma_{3})，首先计算出第一应力不变量I_{1}和第二偏应力不变量J_{2}。假设已知混凝土材料的常数\alpha=0.2，k=1.0。对于应力状态(\sigma_{1}=12MPa,\sigma_{2}=8MPa,\sigma_{3}=6MPa)，计算第一应力不变量I_{1}=12+8+6=26MPa，第二偏应力不变量J_{2}=\frac{1}{6}[(12-8)^{2}+(8-6)^{2}+(6-12)^{2}]=\frac{1}{6}(16+4+36)=\frac{56}{6}\approx9.33MPa^{2}。将I_{1}和J_{2}代入Drucker-Prager准则表达式中，得到F=0.2\times26+\sqrt{9.33}-1.0\approx5.2+3.05-1.0=7.25MPa。多轴强度安全系数K_{DP}定义为：K_{DP}=\frac{\sigma_{0}}{\sigma_{e}}其中，\sigma_{0}为材料的单轴强度指标，\sigma_{e}为等效应力，在Drucker-Prager准则下，等效应力\sigma_{e}可以通过以下公式计算：\sigma_{e}=\sqrt{3J_{2}}+\alphaI_{1}将前面计算得到的值代入等效应力公式，\sigma_{e}=\sqrt{3\times9.33}+0.2\times26\approx5.25+5.2=10.45MPa。假设单轴抗压强度\sigma_{0}=15MPa，则多轴强度安全系数K_{DP}=\frac{15}{10.45}\approx1.44。通过比较安全系数K_{DP}与1的大小，可以判断材料在该应力状态下的稳定性。当K_{DP}\gt1时，材料处于安全状态；当K_{DP}\lt1时，材料可能发生破坏。2.3.3Ottosen准则下的计算Ottosen准则是一种专门针对混凝土材料提出的多轴强度准则，它考虑了混凝土在多轴应力状态下的复杂力学行为，包括拉压异性、中间主应力效应等。该准则能够更准确地描述混凝土在不同应力组合下的强度特性，在混凝土结构的多轴强度分析中具有重要应用。其表达式为：F=\frac{I_{1}^{2}}{f_{c}^{2}}+a_{1}\frac{J_{2}}{f_{c}^{2}}+a_{2}\frac{\sqrt{J_{2}}}{f_{c}}I_{1}+a_{3}\frac{J_{3}}{f_{c}^{3}}+a_{4}\frac{I_{1}^{3}}{f_{c}^{3}}-1=0其中，I_{1}为第一应力不变量，I_{1}=\sigma_{1}+\sigma_{2}+\sigma_{3}；J_{2}为第二偏应力不变量，J_{2}=\frac{1}{6}[(\sigma_{1}-\sigma_{2})^{2}+(\sigma_{2}-\sigma_{3})^{2}+(\sigma_{3}-\sigma_{1})^{2}]；J_{3}为第三偏应力不变量，J_{3}=\frac{1}{27}[(\sigma_{1}-\sigma_{2})^{3}+(\sigma_{2}-\sigma_{3})^{3}+(\sigma_{3}-\sigma_{1})^{3}]；f_{c}为混凝土的单轴抗压强度；a_{1},a_{2},a_{3},a_{4}为材料常数，通过试验确定。在Ottosen准则下，多轴强度安全系数K_{O}的计算基于应力点到强度失效面的距离。对于给定的应力状态(\sigma_{1},\sigma_{2},\sigma_{3})，首先计算出第一应力不变量I_{1}、第二偏应力不变量J_{2}和第三偏应力不变量J_{3}。假设已知混凝土材料的常数a_{1}=-1.0，a_{2}=0.5，a_{3}=0.0，a_{4}=-0.01，单轴抗压强度f_{c}=25MPa。对于应力状态(\sigma_{1}=15MPa,\sigma_{2}=10MPa,\sigma_{3}=8MPa)，计算第一应力不变量I_{1}=15+10+8=33MPa，第二偏应力不变量J_{2}=\frac{1}{6}[(15-10)^{2}+(10-8)^{2}+(8-15)^{2}]=\frac{1}{6}(25+4+49)=\frac{78}{6}=13MPa^{2}，第三偏应力不变量J_{3}=\frac{1}{27}[(15-10)^{3}+(10-8)^{3}+(8-15)^{3}]=\frac{1}{27}(125+8-343)=-\frac{210}{27}\approx-7.78MPa^{3}。将I_{1}、J_{2}、J_{3}代入Ottosen准则表达式中，得到F=\frac{33^{2}}{25^{2}}+(-1.0)\times\frac{13}{25^{2}}+0.5\times\frac{\sqrt{13}}{25}\times33+0.0\times\frac{-7.78}{25^{3}}+(-0.01)\times\frac{33^{3}}{25^{3}}-1，经过计算F\approx1.74-0.02+0.37+0-0.08-1=1.01。多轴强度安全系数K_{O}定义为：K_{O}=\frac{\sigma_{0}}{\sigma_{e}}其中，\sigma_{0}为材料的单轴强度指标，\sigma_{e}为等效应力，在Ottosen准则下，等效应力\sigma_{e}可以通过以下公式计算：\sigma_{e}=\sqrt{\frac{I_{1}^{2}}{f_{c}^{2}}+a_{1}\frac{J_{2}}{f_{c}^{2}}+a_{2}\frac{\sqrt{J_{2}}}{f_{c}}I_{1}+a_{3}\frac{J_{3}}{f_{c}^{3}}+a_{4}\frac{I_{1}^{3}}{f_{c}^{3}}}\timesf_{c}将前面计算得到的值代入等效应力公式，\sigma_{e}=\sqrt{1.74-0.02+0.37+0-0.08}\times25\approx\sqrt{2.01}\times25\approx1.42\times25=35.5MPa。则多轴强度安全系数K_{O}=\frac{25}{35.5}\approx0.7。根据安全系数K_{O}与1的比较结果，可以判断混凝土材料在该应力状态下的稳定性。当K_{O}\gt1时，材料处于安全状态；当K_{O}\lt1时，材料可能发生破坏。2.3.4Hsieh-Ting-Chen准则下的计算Hsieh-Ting-Chen准则也是一种用于描述混凝土多轴强度特性的准则，它综合考虑了混凝土在拉、压、剪等不同受力状态下的力学行为，能够较为全面地反映混凝土在复杂应力条件下的强度变化规律。该准则在混凝土结构的多轴强度分析中得到了一定的应用，为混凝土结构的设计和分析提供了重要的理论依据。其表达式为：F=\frac{I_{1}^{2}}{f_{c}^{2}}+a_{1}\frac{J_{2}}{f_{c}^{2}}+a_{2}\frac{\sqrt{J_{2}}}{f_{c}}I_{1}+a_{3}\frac{J_{3}}{f_{c}^{3}}+a_{4}\frac{I_{1}^{3}}{f_{c}^{3}}+a_{5}\frac{I_{1}J_{2}}{f_{c}^{3}}-1=0其中，I_{1}为第一应力不变量，I_{1}=\sigma_{1}+\sigma_{2}+\sigma_{3}；J_{2}为第二偏应力不变量，J_{2}=\frac{1}{6}[(\sigma_{1}-\sigma_{2})^{2}+(\sigma_{2}-\sigma_{3})^{2}+(\sigma_{3}-\sigma_{1})^{2}]；J_{3}为第三偏应力不变量，J_{3}=\frac{1}{27}[(\sigma_{1}-\sigma_{2})^{3}+(\sigma_{2}-\sigma_{3})^{3}+(\sigma_{3}-\sigma_{1})^{3}]；f_{c}为混凝土的单轴抗压强度；a_{1},a_{2},a_{3},a_{4},a_{5}为材料常数，通过试验确定。在Hsieh-Ting-Chen准则下，多轴强度安全系数K_{HTC}的计算基于应力点到强度失效面的距离。对于给定的应力状态(\sigma_{1},\sigma_{2},\sigma_{3})，首先计算出第一应力不变量I_{1}、第二偏应力不变量J_{2}和第三偏应力不变量J_{3}。假设已知混凝土材料的常数a_{1}=-1.2，a_{2}=0.6，a_{3}=0.05，a_{4}=-0.02，a_{5}=0.1，单轴抗压强度f_{c}=30MPa。对于应力状态(\sigma_{1}=18MPa,\sigma_{2}=12MPa,\sigma_{3}=10MPa)，计算第一应力不变量I_{1}=18+12+10=40MPa，第二偏应力不变量J_{2}=\frac{1}{6}[(18-12)^{2}+(12-10)^{2}+(10-18)^{2}]=\frac{1}{6}(36+4+64)=\frac{104}{6}\approx17.33MPa^{2}，第三偏应力不变量J_{3}=\frac{1}{27}[(18-12)^{3}+(12-10)^{3}+(10-18)^{3}]=\frac{1}{27}(216+8-512)=-\frac{288}{27}\approx-10.67MPa^{3}。将I_{1}、J_{2}、J_{3}代入Hsieh-Ting-Chen准则表达式中，得到F=\frac{40^{2}}{30^{2}}+(-1.2)\times\frac{17.33}{30^{2}}+0.6\times\frac{\sqrt{17.33}}{30}\times40+0.05\times\frac{-10.67}{30^{3}}+(-0.02)\times\frac{40^{3}}{30^{3}}+0.1\times\frac{40\times17.33}{30^{3}}-1，经过计算F\approx1.7\##三、水电站地下厂房工程概况与地质条件\##\#3.1工程基本信息以锦屏二级水电站为例，该电站位于四川省凉山彝族自治州木里、盐源、冕宁三县交界处的雅ç

»æ±Ÿå¹²æµé”¦å±å¤§æ²³æ¹¾ä¸Šï¼Œæ˜¯é›…ç

»æ±Ÿå¡æ‹‰è‡³æ±Ÿå£æ²³æ®µäº”级开发的第二个梯级电站。其装机容量达4800MW，采用地下厂房的布置形式，是以主厂房（344.4m×25.8m×68.7m，含副厂房及装配场）和主变室（314.95m×19.8m×33.0m）为主体的洞室群，厂房长轴按N35°E平行布置。锦屏二级水电站在能源供应中å

据着重要地位。作为“西电东送”的重要电源点之一，其强大的发电能力为东部地区提供了大量的清洁能源，有效缓解了这些地区的用电压力，促进了区域间的能源资源优化配置。据统计，该电站每年可向华东地区输送数十亿千瓦时的电量，为当地的经济发展提供了可é

的电力保障。锦屏二级水电站的建设也对当地的经济和社会发展产生了积极的推动作用。在建设过程中，大量的人力、物力投入，带动了当地建筑、运输、服务等相关产业的发展，增åŠ

了就业机会，提高了当地居民的收入水平。电站建成后，稳定的电力供应为当地的工业发展提供了有力支持，促进了产业结构的优化升级，为区域经济的可持续发展注入了新的活力。\##\#3.2地质条件分析锦屏二级水电站地下厂房洞室群位于三å

系中统盐塘组的第四段条带大理岩（T2y4）、第五段灰黑、灰白大理岩（T2y5）、第六段泥质大理岩（T2y6）中。总体而言，该区域的岩体质量较好，但由于受到多期构é€

改é€

作用和浅表生时效变形的影响，围岩体小断层及裂隙发育，这在一定程度上降低了岩体的强度。该区域的岩性主要为大理岩，其岩石坚硬，抗压强度较高，饱和单轴抗压强度一般在80-120MPa之间，弹性模量较大，通常在20-30GPa左右，这使得岩体在一定程度上具备较好的承载能力。然而，小断层及裂隙的发育ç

´åäº†å²©ä½“的完整性，增åŠ

了岩体的渗透性和变形特性的复杂性。在长期的地质构é€

运动中，这些小断层和裂隙成为了岩体中的薄弱部位，容易导致应力集中和岩体的局部ç

´åã€‚从岩体结构来看，该区域岩体主要呈块状-次块状结构，局部由于断层和裂隙的切割，呈现镶嵌碎裂结构。块状-次块状结构的岩体整体稳定性相对较好，但在受到外部荷载作用时，结构面之间的相互作用可能会导致岩体的变形和ç

´åã€‚镶嵌碎裂结构的岩体则更为ç

´ç¢Žï¼Œå…¶ç¨³å®šæ€§è¾ƒå·®ï¼Œåœ¨åœ°ä¸‹åŽ‚房开挖过程中更容易出现坍塌等问题。该地区的天然应力状态对地下厂房围岩稳定性也有着重要影响。æ

¹æ®ç›¸å…³å‹˜æµ‹èµ„料，该区域的最大主应力方向大致为N30°-40°E，最大主应力值在15-25MPa之间。这种天然应力状态使得地下厂房在开挖后，围岩应力会发生重分布，可能导致洞室周边岩体出现应力集中现象。在洞室的拱顶和边墙部位，由于应力集中，岩体可能会出现拉应力或剪应力超过其强度的情况，从而引发岩体的ç

´åå’Œå˜å½¢ã€‚当最大主应力方向与洞室轴线夹角较大时，边墙部位的应力集中更为明显，可能导致边墙岩体的片帮和坍塌。地下水状况也是影响围岩稳定性的关键å›

ç´

之一。该区域地下水类型主要为基岩裂隙水，赋存于岩体的裂隙中。地下水的存在不仅会降低岩体的有效应力，使岩体的抗剪强度降低，还会åŠ

速岩体的风化和侵蚀作用，进一步削弱岩体的强度。在一些断层和裂隙发育的部位，地下水的渗流可能会带走岩体中的细小颗粒，形成管涌等现象，ç

´åå²©ä½“的结构完整性。在地下厂房开挖过程中，地下水的涌出还可能影响施工进度和安全，增åŠ

支护的难度和成本。综上所述，锦屏二级水电站地下厂房的地质条件较为复杂，岩性、岩体结构、天然应力状态和地下水状况等å›

ç´

相互作用，对围岩稳定性产生了潜在的影响。在进行围岩稳定性分析和工程设计时，必须充分考虑这些地质条件，采取相应的措施来确保地下厂房的安全和稳定。\##\#3.3岩石力学参数测定岩石力学参数的准确测定对于水电站地下厂房围岩稳定性分析至关重要，它直接关系到工程设计的合理性和安全性。本ç

”究采用æ

‡å‡†è¯•éªŒæ–¹æ³•ï¼Œå¯¹é”¦å±äºŒçº§æ°´ç”µç«™åœ°ä¸‹åŽ‚房周边围岩进行了全面的力学参数测定，主要包括单轴和三轴压缩试验、剪切试验等，以获取岩石在不同受力状态下的力学性能数据。单轴压缩试验是岩石力学试验中最基本的试验之一，它能够直观地反æ˜

岩石在单一轴向压力作用下的力学特性。在本次试验中，选取了具有代表性的大理岩岩芯，将其åŠ

工成直径为50mm、高度为100mm的æ

‡å‡†åœ†æŸ±ä½“试件，以满足试验要求。在压力试验机上，以0.5MPa/s的åŠ

载速率对试件施åŠ

轴向压力，通过高精度ä¼

感器实时监测压力和位移数据，直至试件ç

´åã€‚通过单轴压缩试验，成功测定了岩石的抗压强度。试验结果显示，该大理岩的单轴抗压强度平均值为105MPa，这一数值表明岩石具有较高的抗压能力，但在实际工程中，仍需考虑各种å›

ç´

对岩石强度的影响。同时，æ

¹æ®è¯•éªŒè¿‡ç¨‹ä¸­è®°å½•çš„压力-位移数据，绘制出了应力-应变曲线。从曲线中可以看出，在弹性阶段，应力与应变呈现良好的线性关系，这符合胡克定律，表明岩石在该阶段表现出典型的弹性特征。随着应力的逐渐增åŠ

，岩石进入塑性阶段，应变增长速度åŠ

快，此时岩石内部开始出现微裂纹的萌生和扩展。当应力达到峰值强度后，岩石试件发生ç

´åï¼Œåº”力迅速下降，表明岩石的承载能力已达到极限。通过对曲线的分析，还计算得到了岩石的弹性模量和泊松比等重要参数。该大理岩的弹性模量平均值为25GPa，泊松比平均值为0.25。弹性模量反æ˜

了岩石抵抗弹性变形的能力，较高的弹性模量意味着岩石在受力时不易发生弹性变形；泊松比则描述了岩石在横向变形与纵向变形之间的关系，对于分析岩石在复杂受力状态下的变形特性具有重要意义。三轴压缩试验则更能模拟岩石在实际工程中的多轴受力状态，为深入了解岩石在复杂应力条件下的力学行为提供了关键数据。在三轴压缩试验中，同æ

·é‡‡ç”¨æ

‡å‡†åœ†æŸ±ä½“试件，将其放入三轴压力室中。首先，对试件施åŠ

围压，模拟岩石在地下所处的三向应力环境中的侧向压力。然后，以与单轴压缩试验相同的åŠ

载速率施åŠ

轴向压力，直至试件ç

´åã€‚在试验过程中，通过压力ä¼

感器和位移ä¼

感器精确测量围压、轴向压力以及试件的轴向和径向变形。通过三轴压缩试验，获得了不同围压下岩石的抗压强度和变形特性数据。试验结果表明，随着围压的增åŠ

，岩石的抗压强度显著提高。当围压从0MPa增åŠ

到10MPa时，岩石的抗压强度从105MPa提高到150MPa左右。这是å›

为围压的存在限制了岩石内部微裂纹的扩展，增强了岩石颗粒之间的摩擦力和粘结力，从而提高了岩石的整体强度。同时，试验还发现，围压对岩石的变形特性也有明显影响。随着围压的增大，岩石的轴向应变和径向应变都相应减小，岩石的脆性逐渐降低，塑性增强。这意味着在高围压条件下，岩石在ç

´åå‰èƒ½å¤Ÿæ‰¿å—更大的变形，其ç

´åå½¢å¼ä¹Ÿä»Žè„†æ€§ç

´åé€æ¸è½¬å˜ä¸ºå¡‘性ç

´åã€‚岩石的抗剪强度是衡量岩石抵抗剪切ç

´åèƒ½åŠ›çš„重要指æ

‡ï¼Œå¯¹äºŽåˆ†æžåœ°ä¸‹åŽ‚房围岩的稳定性具有关键作用。本ç

”究采用直剪试验来测定岩石的抗剪强度。在直剪试验中，将岩石试件åŠ

工成方形或矩形，放置在直剪仪的上下剪切盒中。通过施åŠ

垂直压力，模拟岩石在实际工程中所受到的法向应力。然后，以一定的速率施åŠ

水平剪切力，使试件在预定的剪切面上发生剪切ç

´åã€‚在试验过程中，通过力ä¼

感器和位移ä¼

感器分别测量垂直压力、水平剪切力以及剪切位移。æ

¹æ®ç›´å‰ªè¯•éªŒç»“果，绘制出了岩石的剪应力-剪切位移曲线和抗剪强度包络线。从剪应力-剪切位移曲线可以看出，在剪切初期，剪应力随着剪切位移的增åŠ

而逐渐增大，呈现出近似线性的关系，这一阶段岩石主要发生弹性变形。当剪应力达到一定值后，曲线开始出现非线性变化，表明岩石内部的结构开始发生ç

´åï¼Œå¾®è£‚纹逐渐扩展和贯通。最终，当剪应力达到峰值时，试件发生剪切ç

´åã€‚通过对试验数据的分析，得到了岩石的黏聚力和内摩擦角等抗剪强度参数。该大理岩的黏聚力平均值为1.2MPa，内摩擦角平均值为35°。黏聚力反æ˜

了岩石颗粒之间的粘结强度，内摩擦角则体现了岩石在剪切过程中抵抗滑动的能力，这两个参数对于评估岩石在不同应力状态下的抗剪强度具有重要意义。通过上述æ

‡å‡†è¯•éªŒæ–¹æ³•ï¼Œå…¨é¢ã€å‡†ç¡®åœ°æµ‹å®šäº†é”¦å±äºŒçº§æ°´ç”µç«™åœ°ä¸‹åŽ‚房周边围岩的岩石力学参数，包括抗压强度、弹性模量、泊松比、黏聚力和内摩擦角等。这些参数为后续的围岩稳定性分析提供了可é

的数据基础，有助于深入了解岩石在不同受力状态下的力学行为，为工程设计和施工提供科学依据。在实际工程应用中，还需考虑岩石的非均质性、节理裂隙等å›

ç´

对力学参数的影响，以确保围岩稳定性分析的准确性和可é

性。\##四、基于多轴强度安全系数的围岩稳定性分析模型建立\##\#4.1数学模型构建采用MATLAB等数值计算软件，以多轴强度安全系数为评价指æ

‡ï¼Œè€ƒè™‘地质条件和岩石力学参数，建立围岩稳定性分析的数学模型。在构建数学模型时，将地下厂房的地质条件进行抽象和简化，转化为数学表达式。把岩体视为连续介质，æ

¹æ®å¼¹æ€§åŠ›å­¦å’Œå¡‘性力学的基本原理，建立应力-应变关系方程。对于锦屏二级水电站地下厂房，其围岩主要为大理岩，æ

¹æ®å‰æ–‡æµ‹å®šçš„岩石力学参数，如弹性模量\(E=25GPa，泊松比\mu=0.25，以及地质条件中的天然应力状态，最大主应力\sigma_{1max}=25MPa，最小主应力\sigma_{3min}=10MPa等信息，构建如下数学模型：假设围岩在三维空间中的应力状态为\sigma_{ij}（i,j=1,2,3），根据弹性力学的广义胡克定律，应变\varepsilon_{ij}与应力\sigma_{ij}的关系为：\varepsilon_{ij}=\frac{1}{E}[(1+\mu)\sigma_{ij}-\mu\sigma_{kk}\delta_{ij}]其中，E为弹性模量，\mu为泊松比，\sigma_{kk}=\sigma_{11}+\sigma_{22}+\sigma_{33}，\delta_{ij}为克罗内克符号，当i=j时，\delta_{ij}=1；当i\neqj时，\delta_{ij}=0。在考虑岩石的强度准则时，以Mohr-Coulomb准则为例，其表达式为：F=\sigma_{1}-\sigma_{3}\frac{1+\sin\varphi}{1-\sin\varphi}-2c\frac{\cos\varphi}{1-\sin\varphi}=0其中，c为黏聚力，\varphi为内摩擦角。对于锦屏二级水电站的大理岩围岩，黏聚力c=1.2MPa，内摩擦角\varphi=35^{\circ}。多轴强度安全系数K基于失效距离的概念定义为：K=\frac{\sigma_{0}}{\sigma_{e}}其中，\sigma_{0}为材料的单轴强度指标，如单轴抗压强度，对于该大理岩，单轴抗压强度\sigma_{0}=105MPa；\sigma_{e}为等效应力，在Mohr-Coulomb准则下，等效应力\sigma_{e}可通过公式\sigma_{e}=\frac{\sigma_{1}-\sigma_{3}}{1-\sin\varphi}计算。在MATLAB中，通过编写程序实现上述数学模型的计算。首先，定义岩石力学参数和地质条件参数，如弹性模量、泊松比、黏聚力、内摩擦角、主应力等。然后，根据广义胡克定律计算应变，再根据Mohr-Coulomb准则计算多轴强度安全系数。以下是部分MATLAB代码示例：%定义岩石力学参数和地质条件参数E=25e9;%弹性模量，单位：Pamu=0.25;%泊松比c=1.2e6;%黏聚力，单位：Paphi=35*pi/180;%内摩擦角，转换为弧度sigma1=20e6;%最大主应力，单位：Pasigma3=10e6;%最小主应力，单位：Pasigma0=105e6;%单轴抗压强度，单位：Pa%计算应变sigma_kk=sigma1+sigma3;epsilon_11=(1/E)*((1+mu)*sigma1-mu*sigma_kk);epsilon_33=(1/E)*((1+mu)*sigma3-mu*sigma_kk);%计算等效应力sigma_e=(sigma1-sigma3)/(1-sin(phi));%计算多轴强度安全系数K=sigma0/sigma_e;disp(['多轴强度安全系数K=',num2str(K)]);E=25e9;%弹性模量，单位：Pamu=0.25;%泊松比c=1.2e6;%黏聚力，单位：Paphi=35*pi/180;%内摩擦角，转换为弧度sigma1=20e6;%最大主应力，单位：Pasigma3=10e6;%最小主应力，单位：Pasigma0=105e6;%单轴抗压强度，单位：Pa%计算应变sigma_kk=sigma1+sigma3;epsilon_11=(1/E)*((1+mu)*sigma1-mu*sigma_kk);epsilon_33=(1/E)*((1+mu)*sigma3-mu*sigma_kk);%计算等效应力sigma_e=(sigma1-sigma3)/(1-sin(phi));%计算多轴强度安全系数K=sigma0/sigma_e;disp(['多轴强度安全系数K=',num2str(K)]);mu=0.25;%泊松比c=1.2e6;%黏聚力，单位：Paphi=35*pi/180;%内摩擦角，转换为弧度sigma1=20e6;%最大主应力，单位：Pasigma3=10e6;%最小主应力，单位：Pasigma0=105e6;%单轴抗压强度，单位：Pa%计算应变sigma_kk=sigma1+sigma3;epsilon_11=(1/E)*((1+mu)*sigma1-mu*sigma_kk);epsilon_33=(1/E)*((1+mu)*sigma3-mu*sigma_kk);%计算等效应力sigma_e=(sigma1-sigma3)/(1-sin(phi));%计算多轴强度安全系数K=sigma0/sigma_e;disp(['多轴强度安全系数K=',num2str(K)]);c=1.2e6;%黏聚力，单位：Paphi=35*pi/180;%内摩擦角，转换为弧度sigm