岩体中普遍存在着断层、节理和裂隙等各种不连续面,这些不连续面往往对岩体的渗流和应力特性起着控制性作用。针对岩体的这一特点,学者们先后提出了多种数值仿真计算方法。归纳起来,这些方法无外乎基于两套体系：以有限元法、有限差分法为代表的基于连续介质力学的方法,以离散元法、块体元法、不连续变形分析方法为代表的基于非连续介质力学的方法。新近发展的无网格法和数值流形法则能同时考虑材料的连续和非连续变形能力。以有限元法为代表的连续介质力学方法在岩土工程中应用得非常普遍。然而由于岩土工程本身的复杂性,其前处理工作费时耗力。本文在有限元法的框架下,为降低网格划分的难度,提出了几种新的或改进的数值方法：不连续有限元法、开挖模拟的动态有限单元法和改进的叠单元法。同时,由于基于非连续介质力学的方法能很好的描述岩体有别于其他材料的渗流和应力特性,本文基于块体单元法,对三维随机裂隙岩体的块体识别方法和渗流、应力及耦合特性进行了研究和探讨。综合起来,论文主要开展了以下工作：(1)提出了大坝／地基系统分析的不连续有限元法。该方法可以将复杂问题的网格划分工作分而治之,系统被分成大坝和地基两个网格独立的部分,通过建基面变形与两侧大坝、基础单元位移的几何相容方程、建基面的本构方程和系统整体的平衡方程来实现不连续网格的耦合。同时,将该方法推广到渗流分析领域,和应力应变分析方法组合成一套完整的计算方法。文中运用一个重力坝的算例说明了本套方法的正确性和有效性。(2)针对复杂岩石边坡的开挖问题,借鉴数值流形方法的思路,提出了岩体开挖模拟的动态有限单元法。该方法在模拟岩体开挖时,建立的有限元网格中仅需考虑结构面的构造和复杂地形,而不用考虑复杂开挖面的形状,开挖面仅作为动态参数实时输入,这大幅度降低了前处理的难度。文中给出了该方法应用到小湾孔雀沟料场边坡的开挖仿真计算中,证明了本文方法的正确性和有有效性。(3)借鉴坝体／地基系统的不连续有限元法和岩体开挖模拟的动态有限元法的思路,提出了改进的叠单元法,提高了其计算精度。并进一步将叠单元法推广到三维空间,提出了三维叠单元法,并运用该方法完成了对色尔古水电站地下厂房场区仅考虑自重作用的初始应力场的求解计算。通过和常规有限元结果的对比,说明了本文方法的精度较高。(4)针对岩体中含有大小各异、数量巨大的裂隙,本文运用蒙特卡罗法完成了岩体中随机裂隙的计算机再现模拟；接着基于矢体的概念,提出了三维随机裂隙岩体的块体识别方法；假定岩块不透水、渗流仅通过裂隙面进行,利用达西定律和变分原理,建立了三维随机裂隙岩体渗流分析的块体单元法。最后通过两个算例分析来表明了本文方法的有效性。(5)将断层、节理、裂隙等统称为裂隙,并视岩石裂隙为一种“充填介质”,进而将有填充物和无填充物的裂隙统一处理,推导了裂隙面的刚度系数、导水系数和法向应力之间的关系式。在此基础上,根据块体单元法的应力、渗流分析的基本原理,运用两场交叉迭代法,建立了三维裂隙岩体的渗流与应力耦合分析模型。通过沙沱重力坝的分析,验证了所提方法的准确性和有效性。