GPU-based parallel computation for discontinuous deformation analysis (DDA) method and its application to modelling earthquake-induced landslide

Abstract Graphic Processing Unit (GPU), as a computing device, has upgraded from single-subject graphical processors to multi-core processors with tremendous computational horsepower. This paper proposes to accelerate the DDA using parallel Jacobi Preconditioned Conjugate Gradient (JPCG) technique on GPUs. Based on the results of two numerical examples, the calculation accuracies of the DDA with serial and parallel solvers are validated, and we found that the DDA with parallel solvers exhibits a much higher execution efficiency. The movement process of Daguangbao landslide triggered by the Wenchuan earthquake is replicated and the modeled deposit pattern coincides well with the actual topography after earthquake.

[1]  Egidio Rizzi,et al.  Analytical and numerical DDA analysis on the collapse mode of circular masonry arches , 2010 .

[2]  Jonas Koko,et al.  Parallel preconditioned conjugate gradient algorithm on GPU , 2012, J. Comput. Appl. Math..

[3]  Jeen-Shang Lin,et al.  Dynamic discrete analysis of an earthquake-induced large-scale landslide , 2009 .

[4]  Guowei Ma,et al.  Modelling rock fracturing and blast-induced rock mass failure via advanced discretisation within the discontinuous deformation analysis framework , 2011 .

[5]  Feng Gao,et al.  Rock creep modeling based on discontinuous deformation analysis , 2013 .

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

[7]  Gang Wang,et al.  Boundary Element Parallel Computation for 3D Elastostatics Using CUDA , 2011 .

[8]  Christian Obrecht,et al.  LBM based flow simulation using GPU computing processor , 2010, Comput. Math. Appl..

[9]  Parastoo Sadeghi,et al.  On optimization of finite-difference time-domain (FDTD) computation on heterogeneous and GPU clusters , 2011, J. Parallel Distributed Comput..

[10]  Arie E. Kaufman,et al.  Implementing lattice Boltzmann computation on graphics hardware , 2003, The Visual Computer.

[11]  Kuo-Ta Hsieh.  Parallelization of EMAP3D based on element-by-element Jacobi preconditioned conjugate gradient method , 2003 .

[12]  Massimo Bernaschi,et al.  A flexible high-performance Lattice Boltzmann GPU code for the simulations of fluid flows in complex geometries , 2010 .

[13]  D. Z. Ding,et al.  FSAI preconditioned CG algorithm combined with GPU technique for the finite element analysis of electromagnetic scattering problems , 2011 .

[14]  Richard E. Goodman,et al.  Two dimensional discontinuous deformation analysis , 1985 .

[15]  Qian Sheng,et al.  Modification of the discontinuous deformation analysis method and its application to seismic response analysis of large underground caverns , 2014 .

[16]  Richard J. Goozee,et al.  Distributed and shared memory parallelism with a smoothed particle hydrodynamics code , 2003 .

[17]  Kai Xu,et al.  A hybrid solution method for CFD applications on GPU-accelerated hybrid HPC platforms , 2016, Future Gener. Comput. Syst..

[18]  Johannes Khinast,et al.  Large-scale CFD–DEM simulations of fluidized granular systems , 2013 .

[19]  Masahiro Chigira,et al.  Landslides induced by the 2008 Wenchuan earthquake, Sichuan, China , 2010 .

[20]  Moncho Gómez-Gesteira,et al.  Optimization strategies for CPU and GPU implementations of a smoothed particle hydrodynamics method , 2013, Comput. Phys. Commun..

[21]  Gordon Erlebacher,et al.  Porting a high-order finite-element earthquake modeling application to NVIDIA graphics cards using CUDA , 2009, J. Parallel Distributed Comput..

[22]  Yossef H. Hatzor,et al.  Dynamic stability analysis of jointed rock slopes using the DDA method: King Herod's Palace, Masada, Israel , 2004 .

[23]  Jian-Hong Wu,et al.  New geometric restriction for the displacement-constraint points in discontinuous deformation analysis , 2017 .

[24]  Masahiro Chigira,et al.  Landslides: Landslides induced by the Wenchuan earthquake , 2012 .

[25]  Shen Defei.  GPU-based preconditioned conjugate gradient method for solving sparse linear systems , 2013 .

[26]  Hao Zhou,et al.  GPU implementation of lattice Boltzmann method for flows with curved boundaries , 2012 .

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

[28]  Feng Zhen.  DYNAMIC RESPONSE CHARACTERISTICS OF DAGUANGBAO LANDSLIDE TRIGGERED BY WENCHUAN EARTHQUAKE , 2012 .

[29]  Jacques M. Bahi,et al.  Parallel sparse linear solver with GMRES method using minimization techniques of communications for GPU clusters , 2014, The Journal of Supercomputing.

[30]  S. Gyimothy,et al.  Parallel Realization of the Element-by-Element FEM Technique by CUDA , 2012, IEEE Transactions on Magnetics.

[31]  Rüdiger Westermann,et al.  Linear algebra operators for GPU implementation of numerical algorithms , 2003, SIGGRAPH Courses.

[32]  Runqiu Huang,et al.  The characteristics and failure mechanism of the largest landslide triggered by the Wenchuan earthquake, May 12, 2008, China , 2012, Landslides.

[33]  C. Westen,et al.  Distribution pattern of earthquake-induced landslides triggered by the 12 May 2008 Wenchuan earthquake , 2010 .

[34]  Guoyang Fu,et al.  Numerical Modeling of Earthquake-Induced Landslide Using an Improved Discontinuous Deformation Analysis Considering Dynamic Friction Degradation of Joints , 2016, Rock Mechanics and Rock Engineering.

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

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

[37]  Chang-Wan Kim,et al.  GPU-based parallel computation for structural dynamic response analysis with CUDA , 2014 .

[38]  Zhen Guo,et al.  Safety evaluation of an ancient masonry seawall structure with modified DDA method , 2014 .

[39]  Lizhong Wang,et al.  Development of discontinuous deformation analysis with displacement-dependent interface shear strength , 2013 .

[40]  Yossef H. Hatzor,et al.  Stability of shallow karstic caverns in blocky rock masses , 2010 .

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

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

[43]  Wolfgang Straßer,et al.  A Parallel Preconditioned Conjugate Gradient Solver for the Poisson Problem on a Multi-GPU Platform , 2010, 2010 18th Euromicro Conference on Parallel, Distributed and Network-based Processing.

[44]  Zhiqiang Xie,et al.  A Jacobi_PCG solver for sparse linear systems on multi-GPU cluster , 2016, The Journal of Supercomputing.

[45]  Eric Darve,et al.  Assembly of finite element methods on graphics processors , 2011 .

[46]  Pat Hanrahan,et al.  Brook for GPUs: stream computing on graphics hardware , 2004, SIGGRAPH 2004.

[47]  Mickeal Verschoor,et al.  Analysis and performance estimation of the Conjugate Gradient method on multiple GPUs , 2012, Parallel Comput..

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

[49]  Y. Im,et al.  Parallel processing of 3D rigid-viscoplastic finite element analysis using domain decomposition and modified block Jacobi preconditioning technique , 2003 .