脆塑性岩石破坏后区力学特性的面向对象有限元与无界元耦合模拟研究岩土工程专业论文.docxVIP

摘 摘 要 现代计算技术在计算能力和存储容量上的革命仅仅提供了计算更复杂问题的有效工具，而程 序的高效性是我们永恒的追求。本文试图从两方面着手，应用面向对象的程序设计方法来提高程 序的开发和维护效率；在有限元分析中纳入无界单元以减小岩土工程的计算区域，从而提高程序 的计算效率，最终为脆塑性岩石破坏后区力学特性的数值模拟服务。本文分别对面向对象的程序 设计方法、无界单元方法以及岩石破坏后区力学特性进行了系统的研究，主要工作如下： 1、在系统地研究了面向对象程序设计方法的基础上，建立了有限元三维非线性分析的面向 对象模型。 2、采用面向对象的程序设计语言Visual c++。嵌)kMatlab数学库，利用其强大的矩阵运算 功能．开发了运行于Windowsxp／2000／M操作系统的。适用于岩土工程二维和三维有限元分析的 面向对象有限元分析软件髓PF跳。该软件能进行弹塑性或弹脆塑性有限元计算，能较好的模拟脆 塑性岩石的力学性能，可以考虑开挖卸荷；程序中加入了节理单元，可以模拟岩体中的不连续面； 程序中还纳入了无界单元，可以模拟岩土工程中经常涉及到的无限和半无限域问题。 3、在对无界单元方法进行了系统而全面的回顾与总结的基础上，分别阐述了映射无界元和 衰减无界元的基本原理及其形函数的具体构造方法，总结了几种常用的映射无界元和衰减无界元 的形函数，并且详细推导了一种简单而又非常实用的无界单元—七节点三维无界元的形函数以及 其它相关计算公式。 4、在前人工作的基础上，对岩石的应力应变全过程曲线分类进行了有益的探讨；应用塑性 位势理论，详细推导了对应于不同屈服准则的应力脆性跌落过程塑性流动因子的确定方法，并且 给出了非理想脆塑性模型应力脆性跌落过程中产生的非零位移增量的一种简便的近似处理方法。 5、分别对大理岩、红砂岩和花岗岩等几种脆性比较明显的岩石进行了应力脆性跌落系数的 试验研究。研究表明脆塑性岩石的脆性是相对的，随着围压的增大。岩石逐渐由脆性向延性转化； 脆塑性岩石的应力脆性跌落在围压不大的情形下发生，其应力脆性跌落系数是围压的函数；并给 出了大理岩和红砂岩的应力脆性跌落系数与围压的关系表达式。 6、对非线性有限元分析中不同屈服条件进行了对比研究。研究表明，在三维情形下，采用 不同的广义yonMises条件来逼近莫尔—库仑条件时，Drueker--Prager条件偏于保守，真正的逼 近Mohr-coul∞b条件应该介于内角圆和外角圆之间，内角圆偏于保守而外角圆偏于危险，所以等 面积圆是相对合适的折中选择。 7、将应用面向对象方法开发的三维弹一脆一塑性有限元与无界元耦台分析软件髓Pf吼对小湾 水电站地下洞室群的围岩稳定性进行了不同屈服条件下的分析。证明了本文所设计的分析软件的 有效性和实用性。对类似的工程具有一定的参考意义． 关键词：面向对象有限元无界元耦合脆塑性岩石破坏后区特性应力脆性跌落屈服准则 Research Research Post-failure Mechanical Properties ofBrittle-plastic Rocks by OOFEM Coupled wRh IEM Abstract Shi Guicai(Geoteehnieal Engineering) Directed by Academician Ge Xiurun 1．Basedoilthe systemic researchofobject-orientedprogrammingmethod，athree-dimensional object-oriented non-linear brittle-plastic nn．把element analysis model is constructed． 3．Based the systemic and comprehensive review and summarize ofinfinite element method． the fundamental principles ofmapping attenuation infinite elements and the specific construction methods of their shape funetions elaborateA respectively,and then shape functions of several mapping attenuation infmite elements in common雠are summarized．and the shape function and some other correlation formulas of simple and utility infinite