Parallel computing of three-dimensional discontinuous deformation analysis based on OpenMP

Abstract The computing efficiency of three-dimensional discontinuous deformation analysis (3D-DDA) needs to be improved for large-scale simulations. Among all the subroutines of 3D-DDA, the equation solver is very time-consuming. To accelerate the equation-solving process, this paper proposes implementing the parallel block Jacobi (BJ) and preconditioned conjugate gradient (PCG) iterative solvers into the original 3D-DDA based on OpenMP. The calculation accuracy and computational efficiency are studied by several numerical examples, demonstrating that the modified 3D-DDA with parallel BJ or PCG solver exhibits much higher execution efficiency with satisfactory correctness. The maximum speedup ratio is up to 5.1 for the cases studied.

[1]  Jian-Hong Wu,et al.  Seismic landslide simulations in discontinuous deformation analysis , 2010 .

[2]  Jian-Hong Wu,et al.  Application of DDA to simulate characteristics of the Tsaoling landslide , 2011 .

[3]  Yousef Saad,et al.  GPU-accelerated preconditioned iterative linear solvers , 2013, The Journal of Supercomputing.

[4]  Chris J. Pearce,et al.  Computational aspects of the discontinuous deformation analysis framework for modelling concrete fracture , 1998 .

[5]  Mustafa M. Aral,et al.  An implicitly coupled hydro-geomechanical model for hydraulic fracture simulation with the discontinuous deformation analysis , 2015 .

[6]  Danesh K. Tafti,et al.  Efficient parallel CFD-DEM simulations using OpenMP , 2014, J. Comput. Phys..

[7]  M. Tsesarsky,et al.  Tunnel roof deflection in blocky rock masses as a function of joint spacing and friction - A parametric study using discontinuous deformation analysis (DDA) , 2006 .

[8]  Yossef H. Hatzor,et al.  Numerical analysis of block stone displacements in ancient masonry structures: A new method to estimate historic ground motions , 2008 .

[9]  Yuzo Ohnishi,et al.  Simulation of the mechanical behavior of inclined jointed rock masses during tunnel construction using Discontinuous Deformation Analysis (DDA) , 2004 .

[10]  Petr Byczanski,et al.  Large scale parallel FEM computations of far/near stress field changes in rocks , 2006, Future Gener. Comput. Syst..

[11]  J. Chern,et al.  Modification of the DDA method for rigid block problems , 1998 .

[12]  Yuzo Ohnishi,et al.  A development of the discontinuous deformation analysis for rock fall analysis , 2005 .

[13]  M. R. Yeung.  Analysis of a mine roof using the DDA method , 1993 .

[14]  Michael Dumbser,et al.  A new 3D parallel SPH scheme for free surface flows , 2009 .

[15]  T. Takahashi,et al.  GPU‐accelerated boundary element method for Helmholtz' equation in three dimensions , 2009 .

[16]  G. H. Shi.  Application of discontinuous deformation analysis on stability analysis of slopes and underground power houses , 2014 .

[17]  M. R. Yeung,et al.  A model of edge-to-edge contact for three-dimensional discontinuous deformation analysis , 2007 .

[18]  Attila M. Zsaki.  Parallel generation of initial element assemblies for two-dimensional discrete element simulations , 2009 .

[19]  Jian-Hong Wu.  New edge-to-edge contact calculating algorithm in three-dimensional discrete numerical analysis , 2008, Adv. Eng. Softw..

[20]  Y. Hatzor,et al.  Modelling dynamic deformation in natural rock slopes and underground openings with DDA: review of recent results , 2011 .

[21]  William Dershowitz,et al.  Modeling mechanical layering effects on stability of underground openings in jointed sedimentary rocks , 2009 .

[22]  Gang Li,et al.  An improved parallel SPH approach to solve 3D transient generalized Newtonian free surface flows , 2016, Comput. Phys. Commun..

[23]  Hong Zheng,et al.  Dual form of discontinuous deformation analysis , 2016 .

[24]  Qiang Xu,et al.  Extension of discontinuous deformation analysis and application in cohesive-frictional slope analysis , 2014 .

[25]  Guangqi Chen,et al.  Applying modified discontinuous deformation analysis to assess the dynamic response of sites containing discontinuities , 2018, Engineering Geology.

[26]  Guangqi Chen,et al.  Effects of near-fault seismic loadings on run-out of large-scale landslide: A case study , 2013 .

[27]  Wei Wang,et al.  A new approach for modeling landslide movement over 3D topography using 3D discontinuous deformation analysis , 2017 .

[28]  Guangyao Li,et al.  Development of parallel explicit finite element sheet forming simulation system based on GPU architecture , 2012, Adv. Eng. Softw..

[29]  Y. Cheng,et al.  Advancements and improvement in discontinuous deformation analysis , 1998 .

[30]  Qinghai Miao,et al.  Parallel Computing of Discontinuous Deformation Analysis Based on Graphics Processing Unit , 2017 .

[31]  Jian Chen,et al.  Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method , 2016 .

[32]  Zheng Han,et al.  DDA validation of the mobility of earthquake-induced landslides , 2015 .

[33]  Christian Vollaire,et al.  Parallel computing for the finite element method , 2007, ArXiv.

[34]  Zheng Han,et al.  Extensions of edge-to-edge contact model in three-dimensional discontinuous deformation analysis for friction analysis , 2016 .

[35]  S. Amir Reza Beyabanaki,et al.  Validation of dynamic block displacement analysis and modification of edge-to-edge contact constraints in 3-D DDA , 2009 .

[36]  Bin Zeng,et al.  GPU-based parallel computation for discontinuous deformation analysis (DDA) method and its application to modelling earthquake-induced landslide , 2017 .

[37]  Rudolf Eigenmann,et al.  Parallel programming with message passing and directives , 2001, Comput. Sci. Eng..

[38]  Jian-Hong Wu,et al.  A new algorithm to identify contact patterns between convex blocks for three-dimensional discontinuous deformation analysis , 2008 .

[39]  Roland W. Freund,et al.  Conjugate Gradient-Type Methods for Linear Systems with Complex Symmetric Coefficient Matrices , 1992, SIAM J. Sci. Comput..

[40]  Hehua Zhu,et al.  A multi-shell cover algorithm for contact detection in the three dimensional discontinuous deformation analysis , 2014 .

[41]  Guang-qi Chen,et al.  Effects of vertical seismic force on initiation of the Daguangbao landslide induced by the 2008 Wenchuan earthquake , 2015 .

[42]  Q. H. Jiang,et al.  A Model of Point-to-Face Contact for Three-Dimensional Discontinuous Deformation Analysis , 2004 .

[43]  Y. Ohnishi,et al.  Practical studies on rockfall simulation by DDA , 2011 .