岩石力学的数值模拟分析.ppt

一、概述 二、解析方法 三、岩石力学中的大规模计算——并行算法 四、不连续变形分析 五、无单元法 六、界面元方法;一、概述; 岩石力学在上个世纪的突出进展在于数值计算方法，即计算机的计算方法在岩土工程中的应用获得了巨大的进展。近数十年来，我国已开展了大量的土木岩土工程。有些工程已属于世界第一流的水平，如水电工程中已建立了大量的混凝土高坝和堆石坝，如世界第三高拱坝二滩拱坝、龙羊峡高坝。; 现今可用于对岩体工程结构进行力学分析的数值方法多种多样，每一种方法有其针对性和特点，对一个具体的问题用数值模拟方法进行分析时，应选择一种最适合该问题的方法进行研究。; 严格地讲，岩石除具有弹性性质外，还具有塑性性质和粘性性质，只不过在特定情况下，某种性质占主导地位而已。在岩体工程实践中，硬岩及应力水平不甚高的中硬岩，其力学性质主要呈现为弹性或弹塑性；高应力环境下的软岩，其力学性质主要呈现为塑性或粘塑性；对于服务时间较长的地下工程，岩石极软或软且应力水平很高，则在计算分析中不能忽视岩石的流变性质。; 为了达到了解整个岩体工程系统的应力及变形规律的目的，各种数值方法均可采用，但以弹、塑性有限单元法或拉格朗日元法最为适宜。这两种方法的单元划分灵活，计算所需参数较少且易获得，软件也易于得到，成本较低。; 局部工程结构的应力及变形分析，若岩石中硬以上，则各种方法均可采用；若岩石软弱，则宜采用能进行大变形分析的拉格朗日元法；若岩体可能发生非连续破坏，则宜采用离散单元法、非连续变形法或流形元法或大变形分析的拉格朗日元法。; 要实现了解岩体结构的破坏特征及动态破坏过程的目的，则只能采用离散单元法、非连续变形法或流形元法，因为这些方法就是针对岩体介质的非连续性而提出的。 此外，对于一个具体问题，是进行平面分析或是进行三维分析，也需做出恰当选择。严格地讲，所有的问题都是三维的，但如果采用平面分析既能达到目的，计算结果误差也不大，为了降低费用和快速方便起见，则以采用平面分析为宜。反之，则???采用三维程序代码进行计算分析。 ; 表1给出了几种数值模拟方法所依据的基本原理、求解方式、离散化方法及其适用条件，可供选择模拟方法时参考。;一、概述 二、岩石力学中的大规模计算——并行算法 三、不连续变形分析 四、无单元法 五、界面元方法;二、岩石力学中的大规模计算——并行计算;近30年计算机发展很快，也推动了有限元的发展，并已应用于实际工程问题。但在遇到实际工程问题时，就有了困难，主要是实际工程量的DOF（自由度）太大。第一，工程对象几何形状复杂，如地下厂房的构造，尾水管内墩子，拱坝和地基往往也涉及很大的范围。第二，在求解一些角度点，如镇墩、排水孔结构等，自由度往往都很大。第三，在考虑岩体的开裂、岩体中的裂隙等，则要求更为庞大的自由度。这些都对计算能力提出了很高的要求，虽然现在的计算机的计算能力有了很大的提高，但是要完成更加复杂的计算任务，为了提高计算的规模和减少计算时间，就需要进行并行计算。几乎所有的大型计算机都是并行机。并行计算对于超大型的数值模型运算已经是一个必不可少的工具。;并行计算在国外得到了很大的发展，自从美国国家航天局（NASA）的A.K.Noor在1975年发表第一篇关于有限元并行性计算的文章以来，有限元并行处理技术几乎与并行计算机同步发展，并且由算法研究发展到了算法、软件和硬件相结合的研究。在国内，并行计算还处于发展阶段，而且硬件相对比较落后。在硬件上主要基于向量机、分布式并行机和共享存储式并行机；在内容上，涉及范围比较广，但在系统性和深度上有待进一步发展，软件的开发距离实际应用和商品化还很远；专门针对有限元并行计算的硬件研究很少。;（一）并行计算的基本概念;对于并行计算机，按照实际应用，可以分为如下的4类[1]： 大规模并行机（Massively Parallel Processors，MPP） 对称多处理器系统（Symmetric MultiProcessors，SMP） 分布存储系统（Distributed Shared Memory system，DSM） 网络工作站机群（Cluster Of Workstations，COW）;并行算法及并行程序开发 并行算法是区别于串行算法的另外一大类型的算法，它是适合于各种并行计算机上求解问题和处理数据的算法。并行程序开发是对给定算法构造并行程序的活动，它要求算法设计者和计算机系统的体系结构的设计者进行广泛的交互。;并行程序的设计不仅编程困难，而且调试和分析更加困难。目前，并行程序设计在各方面都处于一个初级阶段，它不支持一个成熟、稳定和通用的并行程序开发环境。并行开发环境指的是软件环境，它对并行计算的影响比一般串行机要大得多。软件对计