非连续变形分析是针对非连续介质的数值计算方法。该方法存在计算效率较低的问题，而方程组求解是瓶颈。研究非连续变形分析在不同的计算条件下，选择不同的方程组求解方法可获得的计算效率。基于OpenMP与CUDA分别并行实现Jacobi迭代法与Jacobi预处理共轭梯度法，分析问题规模、计算时步长、块体接触关系对方程组求解效率的影响。通过算例测试，获得不同条件下选择方程组求解方法的经验公式，可有效指导非连续变形分析的并行方案设计。
As an advanced numerical analysis method, discontinuous deformation analysis (DDA) shows relatively low efficiency. Simultaneous equations solusion is the bottleneck. This work focuses on how to choose an optimal solver in DDA. Jacobi and Jacobi preconditioned conjugate gradient methods are paralleled with OpenMP and CUDA. Three factors affecting the efficiency, computational scale, time step size, and block contact, are analyzed. Based on the tests, an empirical formula for choosing optimal solver in different computing conditions is obtained.
[1] Shi G H. Discontinuous deformation analysis:a new numerical model for the statics and dynamics of block systems[J]. Engineering Computations, 1992, 9(2):157-168.
[2] 石根华.数值流形方法与非连续变形分析[M]. 北京:清华大学出版社, 1997:92-93.
[3] Bobet A, Fakhimi A, Johnson S, et al. Numerical models in discontinuous media:review of advances for rock mechanics applications[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2009, 135(11):1547-1561.
[4] Zheng H, Jiang W. Discontinuous deformation analysis based on complementary theory[J]. Science in China Series E:Technological Sciences, 2009, 52(9):2547-2554.
[5] Shi G H. Applications of discontinuous deformation analysis (DDA) to rock engineering[M]. Computational Mechanics. Springer Berlin Heidelberg, 2007:136-147.
[6] 刘婷婷, 王颖. 基于数值流形法的计算机动画建模[J]. 中国科学院研究生院学报, 2008, 25(6):843-848.
[7] 付晓东, 盛谦, 张勇慧. 基于OpenMP的非连续变形分析并行计算方法[J]. 岩土力学, 2014, 35(8):2401-2407.
[8] Yousef S. Iterative methods for sparse linear systems[M]. Philadelphia, Pennsylvania, United States:Society for Industrial and Applied Mathematics, 2003:95-102.
[9] Mikola R G, Sitar N. Explicit three dimensional discontinuous deformation analysis for blocky system[C]//47th US Rock Mechanics. Geomechanics Symposium. American Rock Mechanics Association, 2013.
[10] Miao Q, Huang M, Wei Q. Parallel computing of numerical manifold method with openMP[C]//Information, Computing and Telecommunication. Yc-Ict'09. IEEE Youth Conference on. IEEE, 2009:174-177.
[11] Xiao Y, Miao Q, Huang M, et al. Parallel computing of discontinuous deformation analysis based on graphics processing unit[J]. International Journal of Geomechanics, 2016:E4016010.
[12] Ning Y, Yang Z, Wei B, et al. Advances in two-dimensional discontinuous deformation analysis for rock-mass dynamics[J]. International Journal of Geomechanics, 2016:E6016001.
