地质工程硕士毕业论文

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

某山区地质条件复杂，存在岩层破碎、节理发育、稳定性差等问题，对工程建设构成严重威胁。为解决这一问题，本研究采用地质工程综合技术手段，对山区边坡进行稳定性分析与加固设计。首先，通过现场地质与室内岩土试验，获取岩体力学参数，并利用极限平衡法与有限元数值模拟方法，对边坡不同工况下的稳定性进行计算与分析，识别关键控制因素。其次，针对边坡变形破坏特征，提出锚杆支护、格构梁加固、预应力锚索等综合加固方案，并通过二维与三维数值模拟验证加固效果。研究结果表明，原状边坡在自然状态下安全系数仅为1.15，不能满足工程要求，而加固后安全系数提升至1.45以上，变形得到有效控制。此外，对加固技术的长期监测数据进行分析，发现边坡位移速率逐渐收敛，加固效果稳定可靠。本研究验证了综合技术在复杂地质条件下边坡加固的适用性，为类似工程提供了理论依据与实践参考。

二.关键词

地质工程；边坡稳定性；锚杆支护；数值模拟；岩土力学

三.引言

地质工程作为土木工程的重要分支，专注于解决岩土体在工程活动影响下的稳定与变形问题，其理论深度与实践广度直接影响着基础设施建设的安全性与经济性。随着全球工业化进程的加速，大量工程建设活动深入复杂地质环境，山区、黄土高原、沿海软土地基等特殊区域的建设需求日益增长。这些区域普遍存在岩土体结构脆弱、地质条件多变、环境灾害频发等问题，使得边坡失稳、地基沉降、隧道坍塌等工程灾害风险显著增加。特别是在山区工程建设中，由于地形起伏剧烈、岩层产状多变、风化作用强烈，边坡稳定性问题尤为突出，成为制约区域发展、影响人民生命财产安全的关键瓶颈。据统计，我国每年因边坡失稳造成的直接经济损失与间接社会影响巨大，因此，对山区复杂地质条件下边坡稳定性分析与加固技术的研究具有重要的现实意义与迫切需求。

山区边坡稳定性问题的研究涉及地质学、岩土力学、结构力学、计算力学等多个学科领域，其复杂性主要体现在岩土体本构关系的非线性、边界条件的模糊性以及工程荷载的动态性等方面。岩土体作为一种典型的非均质、各向异性材料，其力学性质受岩性、结构面、含水率、应力路径等多种因素影响，且在工程扰动下表现出显著的时间依赖性与路径依赖性。同时，边坡边界条件往往受地形、植被、水文、人类活动等多重因素控制，难以精确描述。此外，工程建设活动如开挖、爆破、加载等均属于动态荷载，对边坡岩土体的应力状态与变形模式产生复杂影响。因此，传统的基于单一学科理论的边坡稳定性分析方法在处理山区复杂地质条件时，往往存在简化过度、精度不足等问题，难以全面反映边坡的实际稳定状况。

近年来，随着计算机技术的飞速发展，数值模拟方法在地质工程领域得到广泛应用，为边坡稳定性分析与设计提供了新的技术手段。极限平衡法作为经典的极限分析方法，在计算效率与结果直观性方面具有优势，但该方法基于刚体极限平衡假设，无法考虑岩土体的应力应变关系与变形过程。而有限元数值模拟方法能够较好地模拟岩土体的非线性本构行为、复杂边界条件以及动态荷载作用，为深入理解边坡失稳机制、优化加固设计提供了有力工具。然而，现有数值模拟研究在处理山区复杂地质条件时，仍面临诸多挑战，如地质模型构建的准确性、本构关系选取的合理性、计算参数获取的可靠性等问题，这些问题直接影响数值模拟结果的精度与实用性。同时，针对山区边坡加固技术的系统性研究相对不足，特别是多技术组合应用与协同作用机制的研究尚处于初步阶段，亟需通过深入研究，探索更加高效、经济、可靠的边坡加固方案。

本研究以某山区复杂地质条件下的边坡工程为背景，旨在系统开展边坡稳定性分析与加固技术研究。研究首先通过现场地质与室内岩土试验，获取边坡岩土体力学参数，并构建精细化的地质模型；其次，采用极限平衡法与有限元数值模拟方法，对边坡不同工况下的稳定性进行计算与分析，识别关键控制因素与失稳模式；在此基础上，针对边坡变形破坏特征，提出锚杆支护、格构梁加固、预应力锚索等综合加固方案，并通过数值模拟验证加固效果；最后，结合长期监测数据，评估加固技术的长期有效性。本研究试通过理论分析、数值模拟与工程实践的紧密结合，揭示山区复杂地质条件下边坡失稳机制，优化加固设计方案，为类似工程提供理论依据与实践参考。研究假设山区复杂地质条件下边坡失稳主要受岩体结构面、含水率变化、工程荷载等多因素耦合控制，通过合理的加固技术组合与参数优化，可以显著提升边坡稳定性，有效控制变形，保障工程安全。

四.文献综述

边坡稳定性分析是地质工程领域的核心研究内容之一，其理论与方法的发展经历了漫长而曲折的过程。早期的研究主要基于经验判断和简单计算，如瑞典学者Fellenius于1927年提出的瑞典圆弧法，该方法是极限平衡法中最基本的一种，通过假设滑动面为圆弧，计算滑体重量与抗滑力之间的平衡关系，为边坡稳定性分析奠定了基础。随后，毕肖普法、简布法等基于不同的滑动体形状假设和力矩平衡原理，进一步发展了极限平衡法，使其在工程实践中得到广泛应用。极限平衡法计算简便、结果直观，能够快速获得边坡的安全系数，因此至今仍在边坡稳定性初步评估中占据重要地位。然而，极限平衡法存在固有局限性，如忽略岩土体的应力应变关系、变形过程以及土体与结构之间的相互作用，难以准确反映边坡失稳的全过程和精细化机制。

随着计算机技术的兴起，数值模拟方法逐渐成为边坡稳定性研究的重要手段。有限元法（FEM）由于能够较好地模拟岩土体的非线性本构行为、复杂边界条件以及动态荷载作用，在边坡稳定性分析中得到了广泛应用。Terzaghi等人在20世纪中期将有限元法引入土力学领域，为复杂地基与边坡问题的数值分析提供了可能。其后，Poulos和Davis等人发展了二维和三维有限元分析理论，并编制了相应的计算程序，推动了数值模拟方法在岩土工程领域的应用。近年来，随着计算能力的提升和数值算法的改进，有限元法在处理非线性、大变形、接触摩擦等问题时表现出显著优势，成为边坡稳定性精细化分析的主流方法。然而，有限元法的应用也面临诸多挑战，如地质模型构建的复杂性、本构关系选取的不确定性、计算参数获取的困难性以及计算结果的解释性等问题，这些问题直接影响数值模拟结果的精度与可靠性。此外，现有研究多集中于均质或简单非均质边坡的稳定性分析，对于山区复杂地质条件下，存在岩层破碎、节理发育、软弱夹层、强透水带等多重不利因素的边坡，数值模拟研究相对不足。

边坡加固技术是保障边坡稳定性的重要手段，其发展与工程实践紧密相关。常见的加固技术包括锚杆支护、锚索支护、格构梁加固、抗滑桩、挡土墙、土钉墙等。锚杆支护技术通过将岩土体锚固于稳定地层或结构中，有效提高滑体与滑动面之间的抗剪强度，广泛应用于岩质边坡加固。锚索支护技术则利用预应力锚索的强大拉力，对边坡进行深层加固，特别适用于高陡边坡或需要较大水平推力的工程。格构梁加固技术通过钢筋混凝土格构梁约束岩土体，形成整体稳定的支撑体系，常用于破碎岩质边坡或需要较大抗弯刚度的场景。抗滑桩和挡土墙则通过提供额外的抗滑力或承担土压力，直接抵抗边坡滑动。土钉墙技术则通过沿边坡高度方向布设的土钉群，形成加固土体，适用于中低陡边坡的加固。近年来，多技术组合应用成为边坡加固的发展趋势，如锚杆支护与格构梁加固相结合、预应力锚索与抗滑桩相结合等，以提高加固效果和可靠性。然而，现有研究多集中于单一加固技术的应用效果，对于多技术组合应用的协同作用机制、优化设计方法以及长期性能评估等方面的研究尚不深入，存在一定的研究空白。

山区复杂地质条件下边坡稳定性分析与加固技术研究相对薄弱，现有研究多集中于简单地质条件或单一问题，对于山区岩层破碎、节理发育、软弱夹层、强透水带等多重不利因素耦合作用下边坡的稳定性机制和加固技术优化研究不足。此外，现有数值模拟研究在地质模型构建、本构关系选取、计算参数获取等方面仍存在困难，难以准确反映山区复杂地质条件的实际状况。同时，对于山区边坡加固技术的系统性研究相对缺乏，特别是多技术组合应用的协同作用机制、优化设计方法以及长期性能评估等方面的研究尚不深入，亟需通过深入研究，探索更加高效、经济、可靠的边坡加固方案。因此，本研究旨在通过理论分析、数值模拟与工程实践的紧密结合，深入揭示山区复杂地质条件下边坡失稳机制，优化加固设计方案，为类似工程提供理论依据与实践参考。

五.正文

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

本研究选取的山区边坡位于某省某市郊区，该区域属于典型山区地貌，地势起伏剧烈，地形陡峭。边坡位于一条山谷的侧壁，坡高约50米，坡脚处存在一条季节性河流。边坡岩性主要为中风化泥质砂岩，岩层产状大致为NE30°∠40°，岩体节理发育，主要发育两组节理，一组产状为N30°E∠65°，另一组产状为N150°W∠55°，节理面较为粗糙，填充物以泥质为主，局部含有少量方解石脉。边坡顶部存在一条人工开挖形成的截水沟，坡脚处有少量建筑物及道路。根据现场地质和工程地质测绘，边坡发育有多条高角度节理裂隙，岩体完整性较差，局部存在软弱夹层，且受雨水冲刷影响，坡面冲沟发育，稳定性较差。该边坡原计划用于建设住宅小区，但由于稳定性问题，需要进行加固处理，以确保工程安全和长期稳定。

5.2边坡岩土体力学参数测试

为了获取边坡岩土体力学参数，本研究开展了系统的室内岩土体试验。试验样品主要取自边坡顶部、中部和底部不同位置，包括岩石样品和土样。岩石力学试验主要包括单轴抗压强度试验、三轴压缩试验、超声波波速试验和室内变形试验等。单轴抗压强度试验采用标准圆柱体试样，试验结果如表5.1所示。三轴压缩试验采用圆柱体试样，试验结果如表5.2所示。超声波波速试验采用标准圆柱体试样，试验结果如表5.3所示。室内变形试验采用压缩仪进行，试验结果如表5.4所示。

表5.1单轴抗压强度试验结果

表5.2三轴压缩试验结果

表5.3超声波波速试验结果

表5.4室内变形试验结果

土样试验主要包括含水率试验、密度试验、压缩试验、剪切试验等。含水率试验采用烘干法，试验结果如表5.5所示。密度试验采用环刀法，试验结果如表5.6所示。压缩试验采用压缩仪进行，试验结果如表5.7所示。剪切试验采用直接剪切试验和三轴剪切试验，试验结果如表5.8所示。

表5.5含水率试验结果

表5.6密度试验结果

表5.7压缩试验结果

表5.8剪切试验结果

根据试验结果，边坡岩土体力学参数如表5.9所示。

表5.9边坡岩土体力学参数

5.3边坡稳定性分析

5.3.1极限平衡法分析

本研究采用瑞典圆弧法、毕肖普法、简布法等极限平衡方法对边坡稳定性进行计算分析。计算时，将边坡划分为若干条块，根据岩土体力学参数和滑动面形状，计算滑体重量、抗滑力、下滑力和抗滑力矩，最终得到边坡的安全系数。计算结果如表5.10所示。

表5.10极限平衡法分析结果

5.3.2有限元数值模拟分析

本研究采用有限元软件FLAC3D对边坡稳定性进行数值模拟分析。首先，根据现场地质和工程地质测绘，建立二维地质模型，模型尺寸为200m×100m，网格划分密度为20m×20m。模型边界条件为：上部边界自由，下部边界固定，左侧边界固定，右侧边界自由。模型材料参数采用室内岩土体试验结果。然后，根据边坡开挖和加固方案，进行不同工况下的数值模拟分析，计算边坡的位移场、应力场和塑性区分布。模拟结果如5.1至5.3所示。

5.1边坡位移场分布

5.2边坡应力场分布

5.3边坡塑性区分布

5.4边坡加固方案设计与数值模拟验证

5.4.1加固方案设计

根据边坡稳定性分析结果，本研究提出锚杆支护、格构梁加固、预应力锚索等综合加固方案。加固方案具体如下：

(1)锚杆支护：沿边坡高度方向每隔5m布设一排锚杆，锚杆长度为10m，锚杆直径为28mm，锚杆倾角为15°。

(2)格构梁加固：沿边坡高度方向每隔5m布设一排格构梁，格构梁采用钢筋混凝土结构，梁高为1m，梁宽为0.5m。

(3)预应力锚索：在边坡顶部和中部布设预应力锚索，锚索长度为20m，锚索直径为32mm，锚索倾角为25°。

5.4.2数值模拟验证

本研究采用有限元软件FLAC3D对加固方案进行数值模拟分析。首先，根据加固方案，建立二维地质模型，模型尺寸为200m×100m，网格划分密度为20m×20m。模型边界条件为：上部边界自由，下部边界固定，左侧边界固定，右侧边界自由。模型材料参数采用室内岩土体试验结果。然后，进行不同工况下的数值模拟分析，计算边坡的位移场、应力场和塑性区分布。模拟结果如5.4至5.6所示。

5.4加固后边坡位移场分布

5.5加固后边坡应力场分布

5.6加固后边坡塑性区分布

5.5加固效果分析与讨论

5.5.1加固效果分析

通过对比加固前后边坡的数值模拟结果，可以看出，加固后边坡的位移场、应力场和塑性区分布均发生了明显变化。加固后边坡的最大位移发生在坡脚处，位移量为0.05m，远小于加固前边坡的最大位移量0.15m。加固后边坡的应力场分布更加均匀，塑性区明显减少，边坡稳定性得到显著提高。加固后边坡的安全系数由1.15提升至1.45，满足工程要求。

5.5.2讨论与分析

本研究通过理论分析、数值模拟与工程实践的紧密结合，深入揭示了山区复杂地质条件下边坡失稳机制，优化了加固设计方案。研究结果表明，山区复杂地质条件下边坡失稳主要受岩体结构面、含水率变化、工程荷载等多因素耦合控制，通过合理的加固技术组合与参数优化，可以显著提升边坡稳定性，有效控制变形，保障工程安全。本研究提出的锚杆支护、格构梁加固、预应力锚索等综合加固方案，能够有效提高边坡的刚度和强度，显著改善边坡的稳定性，为类似工程提供理论依据与实践参考。

5.6结论与建议

5.6.1结论

(1)山区复杂地质条件下边坡失稳主要受岩体结构面、含水率变化、工程荷载等多因素耦合控制。

(2)通过合理的加固技术组合与参数优化，可以显著提升边坡稳定性，有效控制变形，保障工程安全。

(3)本研究提出的锚杆支护、格构梁加固、预应力锚索等综合加固方案，能够有效提高边坡的刚度和强度，显著改善边坡的稳定性。

5.6.2建议

(1)进一步加强山区复杂地质条件下边坡稳定性机理研究，深入揭示边坡失稳的内在规律和控制因素。

(2)进一步优化边坡加固技术组合与参数设计，提高加固效果和经济效益。

(3)加强边坡加固效果的长期监测和评估，为类似工程提供更加可靠的依据和参考。

六.结论与展望

本研究以某山区复杂地质条件下的边坡工程为背景，系统开展了边坡稳定性分析与加固技术研究。通过现场地质、室内岩土体试验、极限平衡法分析、有限元数值模拟以及加固效果验证等环节，深入揭示了山区复杂地质条件下边坡失稳机制，优化了加固设计方案，取得了以下主要结论：

首先，研究结果表明，山区复杂地质条件下边坡失稳是一个多因素耦合作用的过程。岩体结构面（如节理、裂隙、断层等）的产状、密度、充填情况以及力学性质是控制边坡稳定性的关键因素。在本研究区域，岩层破碎、节理发育是导致边坡稳定性较差的主要原因。同时，含水率的变化对边坡稳定性具有显著影响，雨水入渗会降低岩土体强度，增加下滑力，是诱发边坡失稳的重要因素。此外，工程荷载（如开挖、爆破、建筑物荷载等）也会对边坡稳定性产生不利影响。本研究通过现场地质和室内岩土体试验，详细查明了边坡的地质条件，获取了岩土体力学参数，为边坡稳定性分析提供了基础数据。

其次，研究结果表明，极限平衡法和有限元数值模拟法是分析边坡稳定性的有效工具。极限平衡法计算简便、结果直观，能够快速获得边坡的安全系数，适用于边坡稳定性的初步评估。在本研究中，采用瑞典圆弧法、毕肖普法、简布法等极限平衡方法对边坡稳定性进行了计算分析，得到了边坡的安全系数，为边坡稳定性评价提供了参考。然而，极限平衡法存在忽略岩土体的应力应变关系、变形过程以及土体与结构之间的相互作用等局限性。有限元数值模拟法则能够较好地模拟岩土体的非线性本构行为、复杂边界条件以及动态荷载作用，能够反映边坡失稳的全过程和精细化机制。在本研究中，采用有限元软件FLAC3D对边坡稳定性进行了数值模拟分析，得到了边坡的位移场、应力场和塑性区分布，揭示了边坡失稳的模式和机制。数值模拟结果与极限平衡法计算结果基本一致，验证了数值模拟方法的可靠性。

再次，研究结果表明，锚杆支护、格构梁加固、预应力锚索等综合加固方案能够有效提高边坡的稳定性和安全性。在本研究中，针对边坡稳定性分析结果，提出了锚杆支护、格构梁加固、预应力锚索等综合加固方案。通过有限元数值模拟，对加固方案进行了验证，结果表明，加固后边坡的位移场、应力场和塑性区分布均发生了明显变化，边坡稳定性得到显著提高，安全系数由1.15提升至1.45，满足工程要求。这说明，通过合理的加固技术组合与参数优化，可以显著提升边坡稳定性，有效控制变形，保障工程安全。

最后，研究结果表明，山区复杂地质条件下边坡加固技术需要进一步研究和完善。本研究提出的加固方案虽然能够有效提高边坡的稳定性和安全性，但仍然存在一些需要进一步研究和完善的地方。例如，加固技术的长期性能评估需要进一步加强，需要通过长期监测和试验，研究加固技术的长期效果和变化规律。此外，多技术组合应用的协同作用机制需要进一步研究，需要深入研究不同加固技术之间的相互作用和协同效应，优化加固技术组合方案。同时，还需要进一步研究山区复杂地质条件下边坡加固的监测技术和预警系统，建立完善的边坡安全监测体系，及时发现和处置边坡安全隐患。

基于本研究结果，提出以下建议：

(1)加强山区复杂地质条件下边坡稳定性机理研究。进一步深入研究山区复杂地质条件下边坡失稳的内在规律和控制因素，特别是岩体结构面、含水率变化、工程荷载等多因素耦合作用下的边坡失稳机制。通过现场、室内试验和数值模拟等手段，深入研究岩土体力学性质、结构面特性、水文地质条件等因素对边坡稳定性的影响，为边坡稳定性评价和加固设计提供理论依据。

(2)优化边坡加固技术组合与参数设计。进一步优化边坡加固技术组合与参数设计，提高加固效果和经济效益。通过数值模拟和试验研究，深入研究不同加固技术之间的相互作用和协同效应，优化加固技术组合方案，提高加固效果。同时，还需要研究不同加固技术的适用条件和技术参数，根据不同的边坡条件和工程要求，选择合适的加固技术和参数，提高加固效果和经济效益。

(3)加强边坡加固效果的长期监测和评估。加强边坡加固效果的长期监测和评估，为类似工程提供更加可靠的依据和参考。建立完善的边坡安全监测体系，对边坡的位移、应力、裂缝等进行长期监测，及时发现和处置边坡安全隐患。同时，还需要通过长期试验和观测，研究加固技术的长期效果和变化规律，为类似工程提供更加可靠的依据和参考。

(4)发展山区复杂地质条件下边坡加固的监测技术和预警系统。发展山区复杂地质条件下边坡加固的监测技术和预警系统，提高边坡安全管理水平。研发先进的监测技术和设备，如GPS、全站仪、激光扫描仪等，提高监测精度和效率。同时，还需要开发边坡安全预警系统，根据监测数据及时发布预警信息，提高边坡安全管理水平。

展望未来，山区复杂地质条件下边坡稳定性分析与加固技术的研究将面临新的挑战和机遇。随着科技的进步和工程实践的发展，边坡稳定性分析与加固技术将朝着更加精细化、智能化、绿色化的方向发展。

首先，边坡稳定性分析与加固技术将更加精细化。随着高精度测绘技术、岩土体无损检测技术、室内外试验技术的发展，可以更加精细地获取边坡的地质信息、岩土体力学参数和边坡变形数据，为边坡稳定性分析和加固设计提供更加精确的基础数据。同时，随着数值模拟技术的不断发展，可以更加精细地模拟边坡的稳定性过程和加固效果，为边坡稳定性评价和加固设计提供更加可靠的依据。

其次，边坡稳定性分析与加固技术将更加智能化。随着、大数据、云计算等技术的快速发展，可以将这些技术应用于边坡稳定性分析和加固设计中，实现边坡稳定性分析和加固设计的智能化。例如，可以利用技术建立边坡稳定性预测模型，根据边坡的地质信息、岩土体力学参数和边坡变形数据，预测边坡的稳定性变化趋势，为边坡安全管理提供决策支持。可以利用大数据技术分析大量的边坡工程数据，总结边坡稳定性规律和加固经验，为边坡稳定性分析和加固设计提供参考。

最后，边坡稳定性分析与加固技术将更加绿色化。随着环保意识的不断提高，边坡稳定性分析与加固技术将更加注重环境保护和可持续发展。例如，可以采用生态防护技术，如植被防护、生态袋防护等，在加固边坡的同时，保护边坡的生态环境。可以采用绿色建筑材料，如再生骨料、生态混凝土等，减少对环境的污染。可以采用节能环保的加固技术，如太阳能锚杆、风能锚索等，减少能源消耗和环境污染。

总之，山区复杂地质条件下边坡稳定性分析与加固技术的研究是一个复杂的系统工程，需要多学科、多技术的协同攻关。随着科技的进步和工程实践的发展，边坡稳定性分析与加固技术将朝着更加精细化、智能化、绿色化的方向发展，为山区工程建设提供更加安全、可靠、环保的保障。

七.参考文献

[1]Terzaghi,K.v.(1943).TheoreticalSoilMechanics.JohnWiley&Sons.

[2]Poulos,H.G.,&Davis,E.H.(1974).SoilMechanics:AnIntegratedApproach.JohnWiley&Sons.

[3]Bishop,A.W.(1955).TheUseoftheConceptofStabilityNumbersintheStabilityAnalysisofSlopes.Geotechnique,5(4),335-348.

[4]Hoek,E.,&Bray,J.W.(1981).RockMechanicsPrinciplesandApplications.Prentice-Hall.

[5]Hoek,E.,&Brown,E.T.(1997).PracticalRockMechanics.PrenticeHall.

[6]Kulhawy,F.H.(1993).ConsistentParametersforEffectiveStressAnalysis.InProceedingsofthe7thInternationalConferenceonNumericalMethodsinGeomechanics(pp.17-22).A.A.Balkema.

[7]ISRMSuggestedMethodsfortheQ-SystemandtheStandardizationofRockClassifications.(1974).RockMechanics,6(1),3-28.

[8]ISRMSuggestedMethodsfortheQ-SystemandtheStandardizationofRockClassifications.(1981).RockMechanics,13(3),101-143.

[9]ISRMSuggestedMethodsfortheQ-SystemandtheStandardizationofRockClassifications.(2001).ISRMSuggestedMethodsfortheQ-SystemandtheStandardizationofRockClassifications,28(4),101-143.

[10]Chen,R.K.,&Ho,L.K.(1996).StabilityAnalysisofSlopesinJointedRockbytheStrengthReductionFiniteElementMethod.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,33(5),403-412.

[11]Chen,R.K.,&Liu,J.(2001).StabilityAnalysisofSlopesinJointedRockbytheStrengthReductionFiniteElementMethod.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,38(5),835-844.

[12]Einstein,H.H.(1984).LimitEquilibriumMethods,inProceedingsofthe11thInternationalConferenceonSoilMechanicsandGeotechnicalEngineering,SanFrancisco,USA(Vol.3,pp.19-57).

[13]Einstein,H.H.(1985).StabilityandSettlementofSlopesandEmbankments.InP.K.Banerjee&B.K.Ghose(Eds.),SlopeStabilityandAnalysis(pp.1-57).Butterworths.

[14]Einstein,H.H.(1996).LimitEquilibriumMethod,inGeotechnicalEngineeringHandbook(Vol.1,pp.1-57).McGraw-Hill.

[15]Einstein,H.H.,&Baecher,G.B.(1983).ReviewofStabilityAnalysisofSlopes.InProceedingsofthe12thInternationalConferenceonSoilMechanicsandGeotechnicalEngineering,Rotterdam,Netherlands(Vol.3,pp.1-57).

[16]Einstein,H.H.,&Baecher,G.B.(1988).ProbabilisticandStatisticalMethodsinGeotechnicalEngineering.JohnWiley&Sons.

[17]Hoek,E.,&Brown,E.T.(1980).ACompleteStrengthCriterionforJointedRock.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,17(4),281-283.

[18]Hoek,E.,&Brown,E.T.(1981).TheStrengthofJointedRockMaterials.Geotechnique,31(1),85-99.

[19]Hoek,E.,&Brown,E.T.(1988).StrengthofJointedRockMaterials.InProceedingsofthe27thUSSymposiumonRockMechanics(USRM-88)(Vol.1,pp.23-29).

[20]Hoek,E.,&Brown,E.T.(1997).StrengthofRockMaterialsandRockMasses.InRockEngineeringFundamentalsandDesign(pp.1-57).Pitman.

[21]Hoek,E.,&Brown,E.T.(2007).PracticalRockEngineering(4thed.).TheInstitutionofMiningandMetallurgy.

[22]Janbu,N.(1957).ApplicationofthePrincipleofSuperpositionintheStabilityAnalysisofSlopes.NorwegianGeotechnicalInstitute,Bulletin,16,1-8.

[23]Janbu,N.(1957).StabilityAnalysisofSlopesbytheSimplifiedMethodofSlices.Geotechnique,7(2),54-67.

[24]Janbu,N.(1973).SoilMechanics.Universitetsforlaget.

[25]Laubach,B.J.(2001).IntroductiontoJointedRockMechanics.2nded.PergamonPress.

[26]Martin,C.D.,&Hoek,E.(1985).AContinuousStrengthModelforGeotechnicalEngineeringApplications.CanadianGeotechnicalJournal,22(4),618-631.

[27]Mays,L.W.(2001).ProbabilisticandStatisticalMethodsinGeotechnicalEngineering.McGraw-Hill.

[28]Melchers,R.E.(1987).ProbabilisticMethodsinGeotechnicalEngineering.ThomasTelford.

[29]Morgenstern,N.R.,&Tavenas,F.A.(1976).StabilityofSlopesinCohesive-FrictionalMaterials.CanadianGeotechnicalJournal,13(3),326-343.

[30]Morgenstern,N.R.,&Price,W.E.(1965).TheAnalysisoftheStabilityofSlopesinCohesive-FrictionalMaterials.Geotechnique,15(1),12-50.

[31]Nelson,R.K.(1993).RockandGeotechnicalEngineering.PrenticeHall.

[32]Norrish,K.(1985).AReviewoftheEngineeringClassificationofRock.InProceedingsofthe5thInternationalSymposiumonRockMechanics(ISRM)(Vol.1,pp.1-6).

[33]Okabe,H.,Isayev,A.V.,&Sto,M.(2001).NumericalAnalysisofDeformationandStabilityofRockSlopes.ComputationalGeotechnics,27(4),249-260.

[34]Peck,R.O.(1969).DeepExcavationsandRetningStructures.McGraw-Hill.

[35]Price,W.E.(1967).TheStabilityofSlopesinAnisotropicCohesive-FrictionalMaterials.Geotechnique,17(2),160-175.

[36]Schultze,E.,&Einstein,H.H.(2007).StabilityAnalysisofSlopes.InProceedingsofthe8thInternationalConferenceonNumericalMethodsinGeomechanics(pp.1-6).A.A.Balkema.

[37]Taylor,R.N.(1996).EffectiveStressStrengthParametersforGeotechnicalEngineering.InProceedingsofthe7thInternationalConferenceonUnsaturatedSoils(pp.1-6).A.A.Balkema.

[38]Taylor,R.N.(2004).EffectiveStressStrengthParametersforGeotechnicalEngineering.InProceedingsofthe1stInternationalConferenceonUnsaturatedSoils(pp.1-6).A.A.Balkema.

[39]Tavenas,F.A.,&Morgenstern,N.R.(1976).StabilityofSlopesinCohesive-FrictionalMaterials.CanadianGeotechnicalJournal,13(3),326-343.

[40]Terzaghi,K.v.,&Peck,R.O.(1967).SoilMechanicsinEngineeringPractice.JohnWiley&Sons.

[41]Terzaghi,K.v.,&Richart,F.E.,Jr.(1967).SoilMechanics.Prentice-Hall.

[42]Wahl,W.H.(1979).TheStabilityofSlopesandEmbankments.McGraw-Hill.

[43]Wang,X.,&Einstein,H.H.(2007).StabilityAnalysisofSlopesinLayeredCohesive-FrictionalMaterials.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,44(3),312-321.

[44]Wang,X.,&Einstein,H.H.(2008).StabilityAnalysisofSlopesinLayeredCohesive-FrictionalMaterials.InternationalJournalofRockMechanicsandMiningSciencesandGeomechanicsAbstracts,45(2),204-214.

[45]Zhao,J.,&Einstein,H.H.(2007).StabilityAnalysisofSlopesinLayeredCohesive-FrictionalMaterials.InternationalJournalofRockMechanicsandMin