An Efficient Disk-Based Discontinuous Deformation Analysis Model for Simulating Large-Scale Problems

Huang, Gang‐Hai; Xu, Yuan-Zhen; Yi, Xiong‐Wei; Xia, Ming

doi:10.1061/(asce)gm.1943-5622.0001711

Cited by 21 publications

(3 citation statements)

References 43 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…As shown in Figure 2, for each particle in a cell, it only need to check for possible contacts with particles in the same cell and in the neighboring four cells for 2D cases. This nearest neighbor particle searching method has been widely used in the DEM 52 and disk‐based discontinuous deformation analysis, 53 and it excludes the impossible contacts between two long‐range cells, thus greatly improving the computation efficiency. Take a problem with PBC in the x direction as an example (Figure 2), the entire computation domain is divided into 6 × 6 equal‐sized subcells.…”

Section: Methodologiesmentioning

confidence: 99%

A coupled DEM‐IMB‐LBM model for simulating methane hydrate exploitation involving particle dissolution

Xia

Gong

Feng

et al. 2022

Numerical Meth Engineering

Self Cite

View full text Add to dashboard Cite

The coupled discrete element and lattice Boltzmann method using an immersed moving boundary scheme was extended to simulate methane hydrate exploitation involving mass transport and particle dissolution. In this coupled DEM-IMB-LBM model, a new Dirichlet-type thermal boundary condition is extended to simulate moving curved boundaries with constant concentration. A novel periodic boundary including an efficient searching algorithm for particle contact is proposed to reduce the computational cost and boundary effect. Then this model is validated by two numerical examples: a circular particle with concentration convection-diffusion moving in a horizontal channel and mass transport from a cylinder particle in a simple shear flow. The numerical results obtained from the proposed model agree well with previous studies. To further demonstrate the capacity of the proposed model, simulations of methane hydrate exploitation including two formations in marine sediments are carried out. The numerical results indicate that the coupled DEM-IMB-LBM is not only capable of simulating the dissolution of hydrate particles at the grain level, but also recover the sand erosion and migration process in a fundamental perspective during the methane hydrate exploitation process.

show abstract

Section: Methodologiesmentioning

confidence: 99%

A coupled DEM‐IMB‐LBM model for simulating methane hydrate exploitation involving particle dissolution

Xia

Gong

Feng

et al. 2022

Numerical Meth Engineering

Self Cite

View full text Add to dashboard Cite

show abstract

“…The distinct element method (DEM), as a useful technique for discontinuous analysis in rock engineering, has been widely used to study the rock failure instability in many cases (Xiang et al 2018;Shi et al 2020;Gao et al 2020). Besides, the discontinuous deformation analysis (DDA) proposed by Shi and Goodman (Shi 1988;Shi and Goodman 1985) is also an appropriate numerical tool for the discontinuous block-system simulation and successfully used in many underground engineering projects (Gong et al 2018;Xu et al 2019;Huang et al 2020). However, almost none of these reported models, including the non-linear rule-based model (Blair and Cook 1998), the lattice model (Chinaia et al 1997) and the bonded particle model (Potyondy et al 1996), are able to effectively simulate the progressive fracture process of rocks around underground tunnels characterized by the initiation, propagation and coalescence of cracks, which can be easily modeled by the rock failure process analysis (RFPA) method (Tang 1998;Zhu and Tang 2004;Li and Tang 2021).…”

Section: Introductionmentioning

confidence: 99%

Nonlinear deformation and failure characteristics of horseshoe-shaped tunnel under varying principal stress direction

2022

View full text Add to dashboard Cite

To reveal the nonlinear deformation characteristics and failure mechanisms of surrounding rocks of horseshoe-shaped tunnel affected by varying principal stress directions, the physical experiments are carried out based on the similarity theory and control variable method. Simultaneously, on the basis of the statistical strength theory and meso-damage mechanics, a series of 2D numerical models which are able to consider the rock heterogeneity are established to further investigate the mechanical mechanism of damage evolution of surrounding rocks. Seven different kinds of horseshoe-shaped tunnel models are tested and the typical failure modes are analyzed according to the experimental data and numerical simulations. The results show that when the lateral pressure coefficient is small, fractures mainly develop towards the remote maximum principal stress direction; when the pressure difference between the vertical and horizontal directions is small, tunnel surrounding rocks damage seriously; initial damage of surrounding rocks basically occurs at the bottom floor corners and arch shoulders; the process of stress buildup, shadow and transfer is the fundamental mechanical process for the formation of mesoscopic damage and macroscopic failure.Overall, these achievements can provide valuable insights into the nonlinear failure mechanisms of horseshoe-shaped tunnel and will contribute to tunnel support design and stability evaluation in geotechnical engineering.

show abstract

“…Although the DDA is relatively young, it has been widely applied in various aspects of geotechnical engineering, such as tunnel engineering, rock hydraulic fracturing, rock burst analysis, and landslide . However, the computational efficiency of the DDA is much lower than that in the DEM, although there has been great progress within DDA researches over the past three decades . The main reason for this is that the DEM uses an explicit time integration scheme and the second law of Newton is used to calculate the displacement of each block independently, while the DDA often uses an implicit solution method in which simultaneous equilibrium equations are assembled by applying the minimum potential energy principle and solved in each calculation step.…”

Section: Introductionmentioning

confidence: 99%

Highly efficient iterative methods for solving linear equations of three‐dimensional sphere discontinuous deformation analysis

Huang

et al. 2020

Num Anal Meth Geomechanics

Self Cite

View full text Add to dashboard Cite

The efficiency of solving equations plays an important role in implicit-scheme discontinuous deformation analysis (DDA). A systematic investigation of six iterative methods, namely, symmetric successive over relaxation (SSOR), Jacobi (J), conjugate gradient (CG), and three preconditioned CG methods (ie, J-PCG, block J-PCG [BJ-PCG], and SSOR-PCG), for solving equations in threedimensional sphere DDA (SDDA) is conducted in this paper. Firstly, simultaneous equations of the SDDA and iterative formats of the six solvers are presented. Secondly, serial and OpenMP-based parallel computing numerical tests are done on a 16-core PC, the result of which shows that (a) for serial computing, the efficiency of the solvers is in this order: SSOR-PCG > BJ-PCG > J-PCG > SSOR>J > CG, while for parallel computing, BJ-PCG is the best solver; and (b) CG is not only the most sensitive to the ill-condition of the equations but also the most time consuming under both serial and parallel computing. Thirdly, to estimate the effects of equation solvers acting on SDDA computations, an application example with 10 000 spheres and 200 000 calculation steps is simulated on this 16-core PC using serial and parallel computing. The result shows that SSOR-PCG is about six times faster than CG for serial computing, while BJ-PCG is about four times faster than CG for parallel computing. On the other hand, the whole computation time using BJ-PCG for parallel computing is 3.37 hours (ie, 0.061 s per step), which is about 36 times faster than CG for serial computing. Finally, some suggestions are given based on this investigation result. K E Y W O R D Shigh efficiency, iterative method, linear equation solver, three-dimensional discontinuous deformation analysis, sphere DDA | INTRODUCTIONBoth the discontinuous deformation analysis (DDA) method 1 and discrete element method (DEM) 2 are proposed to simulate the mechanical behaviors of discrete block systems. The DDA and DEM have extensive prospects in engineering application since they can automatically trace and simulate block rotation, fracture opening, and complete detachments, which are often included in real jointed rock masses failure process. 3,4 Compared with the DEM, the DDA has more rigorous mathematical theory and higher computational accuracy. 5,6 Although the DDA is relatively young, it has been widely applied in various aspects of geotechnical engineering, such as tunnel engineering, 7 rock hydraulic fracturing, 8,9 rock burst analysis, 10 and landslide. 11 However, the computational efficiency of the DDA is much lower than that in the DEM, although there has been great progress within DDA researches over the past three decades. 12 The main reason for this is that the DEM uses an explicit time integration scheme and the second law of Newton is used to calculate the displacement of each block independently, while the DDA often uses an implicit solution method in which simultaneous equilibrium equations are assembled by applying the minimum potential energy principle and solved in each calcul...

show abstract

An Efficient Disk-Based Discontinuous Deformation Analysis Model for Simulating Large-Scale Problems

Cited by 21 publications

References 43 publications

A coupled DEM‐IMB‐LBM model for simulating methane hydrate exploitation involving particle dissolution

A coupled DEM‐IMB‐LBM model for simulating methane hydrate exploitation involving particle dissolution

Nonlinear deformation and failure characteristics of horseshoe-shaped tunnel under varying principal stress direction

Highly efficient iterative methods for solving linear equations of three‐dimensional sphere discontinuous deformation analysis

Contact Info

Product

Resources

About