现代计算技术在计算能力和存储容量上的革命仅仅提供了计算更复杂问题的有效工具,而程序的高效性是我们永恒的追求。本文试图从两方面着手,应用面向对象的程序设计方法来提高程序的开发和维护效率;在有限元分析中纳入无界单元以减小岩土工程的计算区域,从而提高程序的计算效率,最终为脆塑性岩石破坏后区力学特性的数值模拟服务。本文分别对面向对象的程序设计方法、无界单元方法以及岩石破坏后区力学特性进行了系统的研究,主要工作如下:1、在系统地研究了面向对象程序设计方法的基础上,建立了有限元三维非线性分析的面向对象模型。2、采用面向对象的程序设计语言Visual C++,嵌入Matlab数学库,利用其强大的矩阵运算功能,开发了运行于Windowsxp/2000/NT操作系统的,适用于岩土工程二维和三维有限元分析的面向对象有限元分析软件EBPFEM。该软件能进行弹塑性或弹脆塑性有限元计算,能较好的模拟脆塑性岩石的力学性能,可以考虑开挖卸荷;程序中加入了节理单元,可以模拟岩体中的不连续面;程序中还纳入了无界单元,可以模拟岩土工程中经常涉及到的无限和半无限域问题。3、在对无界单元方法进行了系统而全面的回顾与总结的基础上,分别阐述了映射无界元和衰减无界元的基本原理及其形函数的具体构造方法,总结了几种常用的映射无界元和衰减无界元的形函数,并且详细推导了一种简单而又非常实用的无界单元—6节点三维无界元的形函数以及其它相关计算公式。4、在前人工作的基础上,对岩石的应力应变全过程曲线分类进行了有益的探讨;应用塑性位势理论,详细推导了对应于不同屈服准则的应力脆性跌落过程塑性流动因子的确定方法,并且给出了非理想脆塑性模型应力脆性跌落过程中产生的非零位移增量的一种简便的近似处理方法。5、分别对大理岩、红砂岩和花岗岩等几种脆性比较明显的岩石进行了应力脆性跌落系数的试验研究。研究表明脆塑性岩石的脆性是相对的,随着围压的增大,岩石逐渐由脆性向延性转化;脆塑性岩石的应力脆性跌落在围压不大的情形下发生,其应力脆性跌落系数是围压的函数;并给出了大理岩和红砂岩的应力脆性跌落系数与围压的关系表达式。6、对非线性有限元分析中不同屈服条件进行了对比研究。研究表明,在三维情形下,采用不同的广义Von Mises条件来逼近莫尔—库仑条件时,Drucker-Prager条件偏于保守,真正的逼近Mohr-Coulomb条件应该介于内角圆和外角圆之间,内角圆偏于保守而外角圆偏于危险,所以等面积圆是相对合适的折中选择。7、将应用面向对象方法开发的三维弹-脆-塑性有限元与无界元耦合分析软件EBPFEM对小湾水电站地下洞室群的围岩稳定性进行了不同屈服条件下的分析。证明了本文所设计的分析软件的有效性和实用性,对类似的工程具有一定的参考意义。