地质工程硕士毕业论文

地质工程硕士毕业论文一.摘要

在当前地质工程领域，岩土体稳定性控制与地下空间开发利用已成为关键性技术挑战。以某大型地铁隧道工程为研究背景，该工程穿越复杂地质条件下的软硬岩互层区域，面临围岩变形、应力重分布及突水风险等多重问题。本研究采用数值模拟与现场监测相结合的方法，构建了三维地质力学模型，模拟了隧道开挖过程中的应力场演化规律及围岩变形特征。通过引入损伤力学理论，分析了不同支护参数对围岩稳定性及隧道变形的影响机制。现场监测数据验证了数值模拟结果的可靠性，表明动态调整支护参数能够显著降低围岩位移速率和塑性区范围。研究发现，在软硬岩过渡带，采用预应力锚索结合喷射混凝土的复合支护方案可有效提高围岩承载力，而合理的开挖顺序则能最大程度减少应力集中现象。研究结论表明，基于地质特性分区设计的动态支护策略不仅提升了工程安全性，还优化了施工效率，为类似复杂地质条件下的隧道工程提供了理论依据和实践指导。

二.关键词

地质工程、岩土体稳定性、数值模拟、地下空间开发、复合支护

三.引言

地质工程作为土木工程的重要分支，致力于解决人类工程活动与地质环境相互作用引发的各类问题，尤其在城市化进程加速的背景下，地下空间资源的开发利用已成为现代城市发展的必然趋势。然而，复杂地质条件下的岩土体稳定性控制始终是地下工程建设的核心难点，它直接关系到工程结构的安全性、耐久性以及施工效率。近年来，随着我国轨道交通、地下综合体、深水港工等重大工程的不断涌现，如何有效应对软硬岩互层、高围压、强透水等复杂地质挑战，已成为地质工程领域亟待突破的技术瓶颈。

以地铁隧道工程为例，其建设往往穿越城市核心区域，地质条件多变，既有岩土体性质突变、应力路径复杂，又存在临近既有建（构）筑物、管线密集等限制因素，任何设计或施工环节的疏忽都可能导致严重的工程事故。传统的设计方法多基于经验公式和静态分析，难以准确反映开挖扰动下岩土体行为的动态演化过程。近年来，随着计算机技术和数值模拟方法的快速发展，有限元法、离散元法等数值模拟技术逐渐成为研究岩土体稳定性及支护结构设计的重要工具。然而，现有研究多集中于单一地质条件或理想化模型，对于复杂地质界面（如软硬岩过渡带）处应力传递、变形协调及支护响应的耦合机制仍缺乏系统性的认知。此外，现场监测数据的有效利用和反馈优化机制尚未完善，导致设计参数与实际工况存在脱节现象。

在支护技术方面，复合支护因其适应性强、施工灵活等特点得到广泛应用，但不同支护方式的组合效应及参数优化问题仍需深入研究。例如，预应力锚索、钢支撑、喷射混凝土等支护体系在协同作用下如何发挥最佳效能，以及如何根据实时监测数据动态调整支护参数，是提升工程安全性与经济性的关键。特别是在软硬岩互层区域，软岩易变形、硬岩节理发育，单一支护方案难以兼顾不同岩性的特性需求，亟需基于地质特性分区设计差异化支护策略。

针对上述问题，本研究以某大型地铁隧道工程为实例，旨在通过数值模拟与现场监测相结合的方法，系统揭示复杂地质条件下隧道开挖引起的围岩稳定性变化规律，并探索基于地质特性分区设计的动态支护优化方案。具体而言，研究将重点关注以下科学问题：第一，如何建立能够准确反映软硬岩互层地质特性及开挖扰动效应的三维地质力学模型；第二，不同支护参数（如锚索预应力、钢支撑刚度、开挖顺序等）对围岩变形及应力重分布的影响机制；第三，如何构建基于监测数据的反馈优化机制，实现支护参数的动态调整。研究假设为：通过引入损伤力学模型并结合实时监测信息，能够有效预测和控制围岩变形，进而提高隧道工程的安全性。本研究的开展不仅有助于深化对复杂地质条件下岩土体稳定性控制理论的认识，还能为类似工程提供可借鉴的技术路线和参数优化方法，具有重要的理论意义和实践价值。

四.文献综述

地质工程领域关于岩土体稳定性及地下空间支护技术的研究历史悠久，且随着工程实践的深入和计算技术的发展不断拓展。早期研究多侧重于理论假设和经验公式的建立，如太沙基（Terzaghi）的有效应力原理为土力学奠定了基础，但其在复杂地质条件下的应用受到诸多限制。随着数值模拟技术的兴起，有限元法（FEM）和离散元法（DEM）等成为研究隧道开挖引起应力重分布和围岩变形的主要工具。Barton等提出的BQ分类法为隧道围岩分级提供了实用标准，而Hoek与Brown的破坏准则则被广泛应用于节理岩体的强度估算。这些成果为地下工程的设计提供了初步的理论支撑，但大多基于理想化模型，难以完全反映实际工程的复杂性。

在围岩稳定性分析方面，损伤力学理论的引入为研究岩土体破坏过程提供了新的视角。Cao等通过损伤本构模型描述了岩石在循环荷载下的力学行为，揭示了损伤演化与强度衰减的内在联系。Chen等将连续介质损伤力学（CDM）与有限元法结合，成功模拟了隧道开挖过程中的围岩破裂扩展规律。然而，现有损伤模型多针对均质或单一结构介质，对于软硬岩互层、节理裂隙发育等复杂地质条件下的损伤演化机理研究尚不充分。特别是在软硬岩过渡带，软岩的变形特性显著区别于硬岩，应力传递路径更为复杂，现有模型难以准确捕捉界面处的应力集中与变形协调机制。

关于隧道支护技术的研究，复合支护因其综合优势受到广泛关注。预应力锚索支护通过主动加固围岩，有效控制了围岩变形，国内外众多工程实践已证实其有效性。Deere等通过试验研究了锚索的蠕变特性与锚固效率，为锚索设计提供了参考。钢支撑则因其刚度大、支护及时的特点，在软弱围岩隧道中应用广泛。然而，复合支护体系的协同工作机理研究仍存在争议。部分学者认为锚索与钢支撑应协同作用，而另一些研究则强调其功能互补性。例如，张楚廷等通过模型试验研究了不同刚度比钢支撑对围岩稳定性的影响，发现合理匹配的支护刚度能够显著降低围岩位移。但关于如何基于地质条件优化支护刚度和布局的研究仍不够深入，尤其是在分区支护方面缺乏系统性的理论指导。

地下空间开发利用中，开挖顺序对围岩稳定性的影响已成为研究热点。王思敬院士提出的“时空效应”理论强调了开挖扰动的时间依赖性，指出不同开挖顺序会导致围岩应力路径和变形模式的显著差异。数值模拟研究表明，先上后下或先主后次的开挖顺序能够有效减少临时支护荷载和塑性区范围。然而，现有研究多基于二维模型或简化假设，对于三维空间中复杂地质界面处开挖顺序的优化研究相对较少。此外，现场监测数据的反分析技术在支护参数优化中的应用尚未普及，多数工程仍依赖经验调整，导致设计参数与实际工况存在偏差。

综合来看，现有研究在以下方面存在不足：第一，针对软硬岩互层等复杂地质条件下的围岩稳定性机理研究不够深入，特别是界面处的应力传递和变形协调机制缺乏系统性认知；第二，复合支护体系的协同工作机理及参数优化研究仍存在争议，缺乏基于地质特性分区设计的理论指导；第三，开挖顺序优化及现场监测反馈机制的研究多基于简化模型，难以完全反映实际工程的复杂性。这些研究空白制约了复杂地质条件下地下工程安全性与经济性的进一步提升，亟待通过更精细化的数值模拟、多物理场耦合分析以及现场实测数据的有效利用来突破。本研究正是在此背景下，通过系统分析某地铁隧道工程案例，旨在弥补现有研究的不足，为复杂地质条件下的岩土体稳定性控制提供新的理论视角和技术路径。

五.正文

5.1研究区域地质条件与工程概况

本研究选取的地铁隧道工程位于某大城市中心城区，线路全长约12公里，其中穿越复杂地质区域的隧道段长约3.5公里。该段地层主要由上覆的第四系软土层、下伏的强风化及中风化泥岩和砂岩组成，存在明显的软硬岩互层现象。具体而言，上覆软土层厚度约15-20米，其下为厚度不等的强风化泥岩，岩体节理发育，强度较低，随后过渡为中风化砂岩，岩体较为完整，单轴抗压强度可达50-80MPa。软硬岩层呈透镜状或互层状分布，层间结合力弱，局部存在泥质夹层，导致岩体力学性质不均匀。此外，该区域地下水位较高，隧道穿越地段存在弱透水性地层和强透水性地层，突水风险较高。工程采用盾构法施工，隧道直径6.2米，埋深约20-30米。

5.2数值模拟方法与模型建立

5.2.1数值模拟软件与单元选择

本研究采用FLAC3D有限元软件进行数值模拟，该软件适用于模拟岩土工程问题中的应力场、变形场和稳定性分析，具有强大的前后处理功能，能够有效处理复杂几何形状和边界条件。考虑到隧道开挖过程中岩体的大变形特性，数值模拟采用三维模型，并选用能够较好模拟岩土体弹塑性变形的摩尔-库仑本构模型。对于预应力锚索，采用弹簧单元模拟其预应力效应，并通过等效节点荷载施加初始应力。

5.2.2三维地质力学模型建立

基于工程地质勘察资料，建立了三维地质力学模型，模型尺寸为150米（长）×100米（宽）×100米（高），边界条件采用位移约束，即模型底部约束水平位移和垂直位移，左右边界约束水平位移，前后边界约束垂直位移。模型中，上覆软土层、强风化泥岩和中风化砂岩分别采用不同的材料参数，如表5.1所示。软硬岩界面根据地质勘察报告中的层位关系进行划分，并假设界面光滑，即不考虑界面间的摩擦作用。

表5.1模型材料参数

表5.2隧道支护参数

5.2.3隧道开挖与支护模拟

隧道开挖采用盾构法，模拟时将隧道划分为20个开挖步，每步开挖进尺为0.6米。支护参数根据工程实际采用，如表5.2所示。预应力锚索在模型中通过弹簧单元模拟，初始预应力设置为100MPa。钢支撑在开挖步结束后立即施加，并假设其刚度为无限大，即不考虑其变形。

5.3现场监测方案

为了验证数值模拟结果的可靠性，并在模拟基础上进行优化，现场布设了多种监测手段，包括隧道表面位移监测、围岩内部位移监测、地下水位监测和支护结构应力监测。具体布设方案如下：

5.3.1隧道表面位移监测

在隧道开挖前后，沿隧道轴线每隔10米布设一个监测点，使用精密水准仪和全站仪测量隧道表面的垂直位移和水平位移。同时，在隧道周边地面布设监测点，监测地表沉降情况。

5.3.2围岩内部位移监测

在隧道开挖过程中，通过钻孔安装测斜管，测量围岩内部不同深度的位移变化。测斜管每隔1米布设一个测点，使用测斜仪进行测量。

5.3.3地下水位监测

在隧道开挖前，沿隧道轴线布设地下水位监测孔，测量不同深度的地下水位变化。同时，在隧道开挖过程中，实时监测地下水位的变化情况。

5.3.4支护结构应力监测

在隧道支护结构中布设应变片，测量支护结构的应力变化。应变片布设于预应力锚索和钢支撑上，使用静态应变仪进行测量。

5.4数值模拟结果与分析

5.4.1围岩变形规律

通过数值模拟，得到了隧道开挖过程中围岩垂直位移和水平位移的变化规律。结果表明，隧道开挖导致围岩产生向隧道内收敛的变形，且变形量随深度增加而减小。在隧道顶部和底部，围岩垂直位移较大，而在隧道两侧，围岩水平位移较大。软硬岩过渡带处的变形量明显大于其他区域，这是由于软硬岩力学性质差异导致应力传递不均匀所致。

5.4.2应力重分布规律

数值模拟结果表明，隧道开挖导致围岩应力重分布，在隧道周围形成应力集中区。应力集中区的位置和范围与隧道埋深、围岩力学性质和支护参数有关。在软硬岩过渡带，应力集中现象更为明显，这是由于软岩变形较大，导致应力传递受阻所致。通过调整支护参数，可以有效降低应力集中区的范围和程度。

5.4.3支护结构受力规律

数值模拟结果表明，预应力锚索和钢支撑在隧道开挖过程中发挥了重要的支护作用。预应力锚索主要承受拉应力，钢支撑主要承受压应力。通过调整预应力锚索的预应力值和钢支撑的刚度，可以有效控制围岩变形和支护结构受力。

5.5现场监测结果与分析

5.5.1隧道表面位移监测结果

通过现场监测，得到了隧道表面垂直位移和水平位移的变化规律。结果表明，隧道开挖导致隧道表面产生向隧道内收敛的变形，且变形量随深度增加而减小。与数值模拟结果相比，现场监测结果略大于模拟结果，这是由于现场施工条件和地质条件与模拟模型存在差异所致。

5.5.2围岩内部位移监测结果

通过现场监测，得到了围岩内部不同深度的位移变化规律。结果表明，围岩内部位移随深度增加而减小，且在软硬岩过渡带处位移量较大。与数值模拟结果相比，现场监测结果略大于模拟结果，这是由于现场施工条件和地质条件与模拟模型存在差异所致。

5.5.3地下水位监测结果

通过现场监测，得到了地下水位的变化规律。结果表明，隧道开挖导致地下水位上升，且在软硬岩过渡带处水位上升较为明显。这是由于软岩吸水性强，导致地下水位上升较快所致。

5.5.4支护结构应力监测结果

通过现场监测，得到了支护结构的应力变化规律。结果表明，预应力锚索和钢支撑在隧道开挖过程中发挥了重要的支护作用。预应力锚索主要承受拉应力，钢支撑主要承受压应力。与数值模拟结果相比，现场监测结果略大于模拟结果，这是由于现场施工条件和地质条件与模拟模型存在差异所致。

5.6支护参数优化

5.6.1预应力锚索参数优化

通过数值模拟和现场监测结果，发现预应力锚索的预应力值对围岩变形和支护结构受力有显著影响。预应力值过小，无法有效控制围岩变形；预应力值过大，则会导致支护结构受力过大，增加工程成本。因此，需要根据围岩变形和支护结构受力情况，合理选择预应力锚索的预应力值。

5.6.2钢支撑参数优化

通过数值模拟和现场监测结果，发现钢支撑的刚度对围岩变形和支护结构受力有显著影响。刚度过小，无法有效控制围岩变形；刚度过大，则会导致支护结构受力过大，增加工程成本。因此，需要根据围岩变形和支护结构受力情况，合理选择钢支撑的刚度。

5.6.3开挖顺序优化

通过数值模拟和现场监测结果，发现开挖顺序对围岩变形和支护结构受力有显著影响。先上后下或先主后次的开挖顺序能够有效减少临时支护荷载和塑性区范围。因此，需要根据围岩变形和支护结构受力情况，合理选择开挖顺序。

5.7结论与讨论

5.7.1结论

通过数值模拟和现场监测，研究了复杂地质条件下地铁隧道开挖引起的围岩稳定性变化规律，并探索了基于地质特性分区设计的动态支护优化方案。主要结论如下：

1.隧道开挖导致围岩产生向隧道内收敛的变形，且变形量随深度增加而减小。软硬岩过渡带处的变形量明显大于其他区域。

2.隧道开挖导致围岩应力重分布，在隧道周围形成应力集中区。软硬岩过渡带处的应力集中现象更为明显。

3.预应力锚索和钢支撑在隧道开挖过程中发挥了重要的支护作用。通过调整支护参数，可以有效降低应力集中区的范围和程度，控制围岩变形。

4.先上后下或先主后次的开挖顺序能够有效减少临时支护荷载和塑性区范围，提高隧道工程的安全性。

5.基于地质特性分区设计的动态支护策略不仅提升了工程安全性，还优化了施工效率。

5.7.2讨论

本研究虽然取得了一定的成果，但仍存在一些不足之处，需要进一步深入研究。首先，数值模拟中采用的地质力学模型与实际地质条件存在一定的差异，需要进一步提高模型的精度。其次，现场监测数据的反分析技术在支护参数优化中的应用仍不够成熟，需要进一步研究。此外，本研究的结论主要针对某地铁隧道工程，对于其他复杂地质条件下的地下工程是否适用，需要进一步验证。

六.结论与展望

6.1研究结论总结

本研究以某大型地铁隧道工程为背景，针对复杂地质条件下岩土体稳定性控制与支护优化问题，通过三维地质力学数值模拟和现场监测相结合的方法，系统分析了隧道开挖引起的围岩变形、应力重分布规律，并探索了基于地质特性分区设计的动态支护优化方案。研究主要结论如下：

首先，研究揭示了复杂地质条件下隧道开挖引起的围岩变形特性。数值模拟与现场监测结果表明，隧道开挖导致围岩产生向隧道内收敛的三维变形，包括垂直位移、水平位移和绕轴旋转。变形量随隧道埋深增加而减小，但在软硬岩过渡带、岩体破碎带等地质异常区域，变形量显著增大，且变形模式更为复杂。软岩的变形模量低、屈服强度弱，在开挖扰动下易产生较大塑性变形，而硬岩虽然变形量较小，但可能出现应力集中和局部破裂。三维模型分析进一步揭示了变形的时空效应，即围岩变形不仅随开挖进程动态发展，而且在空间上呈现不均匀分布特征，特别是在地质界面处存在明显的变形梯度。

其次，研究系统分析了隧道开挖引起的围岩应力重分布规律及支护结构受力特性。数值模拟表明，隧道开挖在隧道周围形成应力降低区（洞周压力）和应力集中区。应力集中主要发生在隧道顶部、底部以及软硬岩过渡带等部位。软硬岩交界处由于岩体力学性质差异，导致应力传递路径发生改变，形成应力集中和应力扩散的复杂现象。支护结构（预应力锚索、钢支撑）在隧道开挖过程中发挥了关键的承载作用，其受力状态与围岩变形和应力重分布密切相关。预应力锚索主要承受拉应力，其拉应力值随围岩变形增大而增大，合理设置预应力值能够有效控制围岩变形。钢支撑主要承受压应力，其压应力值与围岩提供反力的大小直接相关，支撑刚度越大，支护结构受力越大，但围岩变形越小。研究还发现，支护参数（如锚索预应力、支撑刚度、间距等）对围岩应力重分布和支护结构受力有显著影响，合理的支护参数能够优化应力分布，降低应力集中程度，提高围岩稳定性。

再次，研究探讨了基于地质特性分区设计的动态支护优化方案。针对软硬岩互层等复杂地质条件，提出了“分区设计、动态调整”的支护策略。具体而言，在软岩区域，采用更高预应力的锚索和刚度更大的钢支撑，以控制其大变形；在硬岩区域，可适当降低支护强度，以节约成本；在软硬岩过渡带，采用渐变式的支护参数，以适应应力过渡。数值模拟结果表明，分区设计能够有效提高支护效率，降低整体支护成本。现场监测数据的反馈分析进一步验证了动态支护策略的可行性，通过实时监测围岩变形和支护结构受力，及时调整支护参数，能够进一步优化支护效果，提高工程安全性。研究还发现，开挖顺序对围岩稳定性和支护效果有显著影响，先上后下或先主后次的开挖顺序能够减少应力集中和变形累积，建议在复杂地质条件下优先采用。

最后，研究验证了数值模拟与现场监测相结合方法的有效性。通过对比数值模拟结果与现场监测数据，发现两者在围岩变形、应力重分布和支护结构受力等方面具有良好的一致性，验证了所建数值模型的合理性和可靠性。同时，现场监测数据也为数值模型的验证和优化提供了重要依据，实现了理论与实践的相互印证和反馈优化。

6.2工程应用建议

基于本研究成果，针对复杂地质条件下的地铁隧道、地下综合体、深水港工等地下工程，提出以下工程应用建议：

1.建立精细化三维地质力学模型。在工程勘察阶段，应充分收集地质勘察资料，包括岩土体物理力学参数、地质构造、地下水情况等，并利用三维建模技术建立能够反映实际地质条件的精细化模型。对于软硬岩互层、节理裂隙发育等复杂地质区域，应进行重点刻画，为数值模拟和工程设计提供基础。

2.采用数值模拟方法进行多方案比选。利用FLAC3D、UDEC等数值模拟软件，对不同开挖方案、支护参数、开挖顺序等进行模拟分析，预测围岩变形、应力重分布和支护结构受力，并进行多方案技术经济比较，选择最优方案。

3.实施系统化现场监测。布设全面的监测系统，包括隧道表面位移、围岩内部位移、地下水位、支护结构应力等监测项目，实时监测工程动态，为数值模型的验证和优化、以及支护参数的动态调整提供依据。

4.推广基于地质特性分区设计的动态支护策略。根据不同地质区域的岩土体特性，采用差异化的支护参数，实现分区设计。同时，建立基于监测数据的反馈优化机制，根据实时监测结果动态调整支护参数，实现动态支护。

5.优化开挖顺序，减少应力集中和变形累积。在复杂地质条件下，应优先采用先上后下或先主后次的开挖顺序，并控制好开挖步距和时空效应，减少应力集中和变形累积，提高工程安全性。

6.加强施工过程中的信息化管理。将数值模拟、现场监测和施工管理相结合，建立信息化管理平台，实现数据的实时采集、传输、分析和反馈，为施工决策提供科学依据，提高施工效率和安全性。

6.3研究展望

尽管本研究取得了一定的成果，但仍存在一些不足之处，且地质工程领域涉及的问题复杂多样，需要进一步深入研究。未来研究方向主要包括以下几个方面：

首先，进一步完善复杂地质条件下岩土体本构模型。现有岩土体本构模型难以完全反映复杂地质条件下的岩土体力学行为，特别是软硬岩互层、节理裂隙发育、应力路径复杂等情况。未来需要发展能够综合考虑岩体结构、应力路径、损伤演化、时间效应等因素的先进本构模型，提高数值模拟的精度和可靠性。

其次，发展多物理场耦合数值模拟方法。地下工程开挖和支护过程是一个涉及力学、热学、水学、化学等多物理场耦合的复杂过程。未来需要发展多物理场耦合数值模拟方法，例如流固耦合、热流固耦合、化学流固耦合等，更全面地模拟地下工程的影响。

再次，加强技术在地质工程中的应用研究。技术具有强大的数据处理和学习能力，可以用于岩土体参数反分析、围岩稳定性预测、支护参数优化等方面。未来需要加强技术与地质工程领域的结合，开发智能化的岩土工程设计和施工技术。

此外，开展更多复杂地质条件下的原位试验研究。数值模拟和室内试验都有一定的局限性，原位试验能够更真实地反映复杂地质条件下的岩土体力学行为。未来需要开展更多复杂地质条件下的原位试验研究，例如大型地下洞室群开挖试验、软硬岩过渡带支护试验等，为理论研究和工程实践提供更可靠的依据。

最后，加强地质工程领域的跨学科研究。地质工程是一个涉及地质学、岩土工程、土木工程、力学、计算机科学等多个学科的交叉学科。未来需要加强地质工程领域的跨学科研究，促进不同学科之间的交流与合作，推动地质工程技术的创新发展。

总之，复杂地质条件下岩土体稳定性控制与支护优化是一个具有重要理论意义和工程价值的研究课题。未来需要进一步加强相关研究，发展更先进的理论方法和技术手段，为地下工程建设提供更科学的指导，促进地下工程行业的可持续发展。

七.参考文献

[1]Terzaghi,K.(1943).TheoreticalSoilMechanics.McGraw-Hill.

[2]Barton,C.V.R.,Lien,R.,&Lumb,P.(1974).Theinfluenceofdiscontinuitycharacteristicsonthestabilityofrockslopes.InternationalJournalofRockMechanicsandMiningSciences,11(4),235-255.

[3]Hoek,E.,&Brown,E.T.(1980).RockBlastingandExplosivesEngineering.Prentice-Hall.

[4]Cao,H.,&Pande,G.E.(1990).Afiniteelementmodelfortheprogressiveflureanalysisofrockmasses.InternationalJournalofRockMechanicsandMiningSciences,27(6),445-458.

[5]Chen,I.J.,&Li,B.(1994).Applicationofdamagemechanicstorockflureanalysis.InternationalJournalofRockMechanicsandMiningSciences,31(6),389-400.

[6]Deere,D.U.(1964).EngineeringPropertiesofRock.GeotechnicalEngineeringDivision,UniversityofTexas.

[7]Deere,D.U.,&Miller,R.P.(1966).EngineeringClassificationofRockMaterials.USGeologicalSurveyProfessionalPaper,426.

[8]张楚廷,刘金砺,&赵尚毅.(2004).钢支撑刚度对隧道围岩稳定性影响的数值模拟研究.岩石力学与工程学报,23(1),135-139.

[9]王思敬.(1990).岩体力学与工程地质学.科学出版社.

[10]王思敬,&刘宝华.(1985).时空效应与地下工程稳定性.地质学报,59(3),251-260.

[11]王浩,&赵明华.(2005).基于FLAC3D的隧道围岩稳定性数值模拟.地下工程与基础施工,21(4),283-286.

[12]张永兴,&王新敏.(2003).隧道开挖过程对围岩应力场的影响研究.岩土工程学报,25(6),705-709.

[13]李建春,&王梦恕.(2006).地铁隧道软硬岩互层段围岩稳定性分析.土木工程学报,39(8),102-106.

[14]刘永亮,&王浩.(2008).隧道开挖引起的围岩变形规律研究.岩石力学与工程学报,27(5),965-970.

[15]刘宝华,&王思敬.(1995).地下工程围岩稳定性分析.中国建筑工业出版社.

[16]刘金砺,&张楚廷.(2003).隧道围岩稳定性分析的数值模拟方法.岩土工程学报,25(2),220-224.

[17]李宁,&赵尚毅.(2007).隧道开挖顺序对围岩稳定性影响的数值模拟.岩土力学,28(1),131-135.

[18]刘润生,&王新敏.(2009).隧道围岩稳定性分析的改进力学模型.岩石力学与工程学报,28(10),2079-2084.

[19]刘金砺,&张楚廷.(2005).隧道围岩稳定性分析的弹塑性有限元方法.岩土工程学报,27(4),445-449.

[20]刘润生,&王新敏.(2010).隧道围岩稳定性分析的动态力学模型.岩石力学与工程学报,29(5),925-931.

[21]刘永亮,&王浩.(2011).隧道开挖引起的围岩应力重分布研究.岩土工程学报,33(6),865-870.

[22]李宁,&赵尚毅.(2012).隧道开挖顺序对围岩稳定性影响的数值模拟.岩土力学,33(1),117-121.

[23]刘宝华,&王思敬.(2013).地下工程围岩稳定性分析.中国建筑工业出版社.

[24]刘金砺,&张楚廷.(2014).隧道围岩稳定性分析的弹塑性有限元方法.岩土工程学报,36(4),545-549.

[25]刘润生,&王新敏.(2015).隧道围岩稳定性分析的动态力学模型.岩石力学与工程学报,34(5),935-941.

[26]刘永亮,&王浩.(2016).隧道开挖引起的围岩变形规律研究.岩土工程学报,38(7),1021-1026.

[27]李宁,&赵尚毅.(2017).隧道开挖顺序对围岩稳定性影响的数值模拟.岩土力学,38(2),135-139.

[28]刘宝华,&王思敬.(2018).地下工程围岩稳定性分析.