The Research On Parallel Compuatation Method Of Rock Failure Process Analysis

With large and superlarge rock engineering development, numerical analysis of rock mass structure instability and failure process highly demands powerful computation ability, while serial numerial method can not satisfy engineering and science research demands.So large scale rock structure failure process numerical analysis system based on parallel computation technique need to be developed to fulfill the demand of rock mechanics and engineering and to study the micro-failure mechanism of rock material and assess engineering stability. On the basis of our group work, focusing on the theme, this thesis has done the following work mainly:1.The parallel solver of Jacobi preconditioned Biconjugate Gradient Stabilized method (J-BICGSTAB) is implemented in the distributed memery environment. The solver makes use of the communication and computation overlapping in order to reduce the latency of the communication. The solver with low memory requirement and high calculating speed is able to solute the large scale linear equation system.2.The implementations of nonoverlapping grid partition,data partition and information management in fininte element domain decomposition are discussed. Through domain decomposition,fininte element information of subdomains and communication environment of parallel computation are constructed.The stored data of submeshes are localized.3.On the basis of the JBICGSTAB method and domain decomposition method, three-dimensional linear elastic finite element program is coded using MPICH and Fortran90 language. Numerical tests indicated that the code has high performance and good scalability.4.Distributed rock failure process parallel analysis system based on client-server model (RFPA3D-Parallel) is proposed.The PC installed Windows XP serves as client and a cluster whose os is Linux serves server offering computating resource.The client-server model is finished by connection-oriented protocol socket service.So the Windows system and Linux system is able to work together and rock failure process parallel analysis system is implemented.At present, the parallel program can be used to fulfill rock failure process analysis of models with more than ten million degrees of freedom with 24 CPUs on a 32-node cluster.5.Numerial simulation of fracture spacing in three-layered rock materials models with different principal stress ratios is done with RFPA3D-Parallel.Three types of numerical tests are performed:directional stretch, anisotropic stretch and isotropic stretch. The modeling results reproduce the complete processes from the crack initiation and propagation to coalescence forming different fracture pattern. The numerical simulation results of nine principle stress ratios show pattern transition from parallel to polygonal fratures.The modeling results indicate that RFPA3D-Parallel is a powerful numerical tool to analyze large scale rock failure process and has abroad and potential application prospects.