岩石破裂过程分析并行计算方法研究

【摘要】： 随着我国大型、超大型岩土工程的建设,岩石材料脆性破坏与工程结构灾害的研究面临着严重的挑战。岩石工程结构尺度巨大,涉及大量多相细观结构、宏观结构和界面,岩石工程中岩体结构失稳和破裂过程数值分析必须采用巨大的单元数目,需要巨大的计算容量和计算速度,而且岩石工程中的破裂问题数值模拟是一个从变形起始到最终破坏的漫长过程,这个过程对数值模拟而言周期非常长,需要极高的计算速度,对计算能力提出了较高的要求,而传统的单机串行分析方法无法满足工程需要,因而需要发展基于并行计算技术的大规模岩石结构破坏过程数值分析软件系统,以满足岩石力学和工程发展的需求,深入分析岩石材料细观破坏本质、进行工程稳定性评价。本论文在课题组原有工作基础之上,围绕岩石破裂过程并行分析这一主题,开展了以下几个方面的工作： 1.完成了基于Jacobi预处理稳定双共轭梯度法的分布式集群环境下的并行求解器。并行求解器能够充分利用计算和通信的重叠以达到降低通信延迟的目的,其对存储要求较低,计算速度较快,能够完成大规模线性方程组的求解。 2.讨论了有限元区域分解实施过程中有限元的网格无重叠划分、数据分割以及信息管理实现细节,构建了有限元计算子区域的有限元信息以及数据的局部化存储,并同时构建了并行计算的通信环境。 3.基于有限元区域分解方法和并行Jacobi预处理稳定双共轭梯度法,采用MPICH和Fortran编制了三维弹性有限元并行计算程序。算例表明并行计算程序具有较高的性能和良好的可扩展性。 4.设计了基于客户机-服务器模型的分布式岩石破裂过程分析并行分析系统(RFPA3D-Parallel),采用WindowsXP PC机作为客户机,并行计算机作为服务器提供计算资源。通过面向连接协议的socket服务实现了客户机-服务器模型,完成了跨平台系统的信息通讯和数据传递功能,使得WindowsXP和Linux能够协调处理,从而实现了面向网络的岩石破裂过程并行分析系统。目前在所使用的32计算结点集群环境下,并行计算程序使用12个计算结点24个CPU就可以完成300万六面体单元近千万自由度的岩石破裂过程分析。 5.通过应用RFPA3D-Parallel数值分析系统实现了不同主应力水平下三层层状岩石结构模型中等间距裂纹数值模拟。一共进行了直接拉伸、双向不等应力拉伸和双向等应力拉伸三种不同类型9种不同主应力比条件下的裂纹等间距数值模拟。数值模拟再现了裂纹的萌生、扩展以及贯通,最终形成不同模式的裂纹的整个形成过程。9种不同主应力比条件下的裂纹等间距数值模拟结果表明随主应力水平不同,裂纹模式由平行裂纹模式转变为多边形裂纹模式。从而表明RFPA3D-Parallel数值分析系统在岩石工程结构破坏分析领域具有广泛潜在的应用前景,它为岩石工程结构破裂过程分析提供了一个强有力的数值分析工具。