Simulation of quasi-static axisymmetric collapse of granular columns using smoothed particle hydrodynamics and discrete element methods

Elnaz Kermani, Tong Qiu

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

Flow of granular materials under different conditions poses significant challenges to engineers as their behavior can vary from solid-like to fluid-like. The quasi-static axisymmetric collapse of granular columns is an interesting case of granular flow relevant to the active deformation of retaining structures. This study aims to shed light on this type of collapse using numerical simulations including the smoothed particle hydrodynamics (SPH) method in the continuum framework and the discrete element method (DEM) in the discrete framework. Three-dimensional SPH and DEM models are developed to investigate the quasi-static axisymmetric collapse of granular columns, which is initiated by slowly expanding the cylindrical wall of the column. The Drucker–Prager constitutive model with non-associated flow rule is implemented into SPH formulations to model elastoplastic behavior of granular material. The linear contact model is used in DEM simulations with calibrated contact stiffness parameters. The effects of particle shape and hysteretic contact behavior in the DEM model are indirectly considered through reducing particle rotational velocities by a calibrated constant factor at every time step. Using the developed models, the final deposit profile (e.g., height, runout distance, repose angle), non-deformed region, collapse pattern and energy dissipation for different initial aspect ratios are investigated. In addition, final height, runout distance and energy dissipation are theoretically derived and are found to be in good agreement with SPH and DEM simulations. Results show that the quasi-static axisymmetric collapse can be modeled using both the continuum and discrete frameworks as long as appropriate constitutive models or contact behaviors are implemented. However, DEM simulations are more capable of capturing particle-level behaviors, such as the sharp edges of the free surface and flow front. Results demonstrate that quasi-static granular collapse is qualitatively similar to the dynamic collapse, but results in smaller runout distances for all aspect ratios.

Original languageEnglish (US)
JournalActa Geotechnica
DOIs
StateAccepted/In press - Jan 1 2018

Fingerprint

discrete element method
Finite difference method
Particles (particulate matter)
Hydrodynamics
hydrodynamics
simulation
Granular materials
Constitutive models
energy dissipation
Aspect ratio
Energy dissipation
particle
Deposits
stiffness
Stiffness
Engineers
Fluids
Computer simulation
fluid

All Science Journal Classification (ASJC) codes

  • Geotechnical Engineering and Engineering Geology
  • Earth and Planetary Sciences (miscellaneous)

Cite this

@article{777a87f60bdf4591bd4d5a7ec8b4d33d,
title = "Simulation of quasi-static axisymmetric collapse of granular columns using smoothed particle hydrodynamics and discrete element methods",
abstract = "Flow of granular materials under different conditions poses significant challenges to engineers as their behavior can vary from solid-like to fluid-like. The quasi-static axisymmetric collapse of granular columns is an interesting case of granular flow relevant to the active deformation of retaining structures. This study aims to shed light on this type of collapse using numerical simulations including the smoothed particle hydrodynamics (SPH) method in the continuum framework and the discrete element method (DEM) in the discrete framework. Three-dimensional SPH and DEM models are developed to investigate the quasi-static axisymmetric collapse of granular columns, which is initiated by slowly expanding the cylindrical wall of the column. The Drucker–Prager constitutive model with non-associated flow rule is implemented into SPH formulations to model elastoplastic behavior of granular material. The linear contact model is used in DEM simulations with calibrated contact stiffness parameters. The effects of particle shape and hysteretic contact behavior in the DEM model are indirectly considered through reducing particle rotational velocities by a calibrated constant factor at every time step. Using the developed models, the final deposit profile (e.g., height, runout distance, repose angle), non-deformed region, collapse pattern and energy dissipation for different initial aspect ratios are investigated. In addition, final height, runout distance and energy dissipation are theoretically derived and are found to be in good agreement with SPH and DEM simulations. Results show that the quasi-static axisymmetric collapse can be modeled using both the continuum and discrete frameworks as long as appropriate constitutive models or contact behaviors are implemented. However, DEM simulations are more capable of capturing particle-level behaviors, such as the sharp edges of the free surface and flow front. Results demonstrate that quasi-static granular collapse is qualitatively similar to the dynamic collapse, but results in smaller runout distances for all aspect ratios.",
author = "Elnaz Kermani and Tong Qiu",
year = "2018",
month = "1",
day = "1",
doi = "10.1007/s11440-018-0707-9",
language = "English (US)",
journal = "Acta Geotechnica",
issn = "1861-1125",
publisher = "Springer Verlag",

}

TY - JOUR

T1 - Simulation of quasi-static axisymmetric collapse of granular columns using smoothed particle hydrodynamics and discrete element methods

AU - Kermani, Elnaz

AU - Qiu, Tong

PY - 2018/1/1

Y1 - 2018/1/1

N2 - Flow of granular materials under different conditions poses significant challenges to engineers as their behavior can vary from solid-like to fluid-like. The quasi-static axisymmetric collapse of granular columns is an interesting case of granular flow relevant to the active deformation of retaining structures. This study aims to shed light on this type of collapse using numerical simulations including the smoothed particle hydrodynamics (SPH) method in the continuum framework and the discrete element method (DEM) in the discrete framework. Three-dimensional SPH and DEM models are developed to investigate the quasi-static axisymmetric collapse of granular columns, which is initiated by slowly expanding the cylindrical wall of the column. The Drucker–Prager constitutive model with non-associated flow rule is implemented into SPH formulations to model elastoplastic behavior of granular material. The linear contact model is used in DEM simulations with calibrated contact stiffness parameters. The effects of particle shape and hysteretic contact behavior in the DEM model are indirectly considered through reducing particle rotational velocities by a calibrated constant factor at every time step. Using the developed models, the final deposit profile (e.g., height, runout distance, repose angle), non-deformed region, collapse pattern and energy dissipation for different initial aspect ratios are investigated. In addition, final height, runout distance and energy dissipation are theoretically derived and are found to be in good agreement with SPH and DEM simulations. Results show that the quasi-static axisymmetric collapse can be modeled using both the continuum and discrete frameworks as long as appropriate constitutive models or contact behaviors are implemented. However, DEM simulations are more capable of capturing particle-level behaviors, such as the sharp edges of the free surface and flow front. Results demonstrate that quasi-static granular collapse is qualitatively similar to the dynamic collapse, but results in smaller runout distances for all aspect ratios.

AB - Flow of granular materials under different conditions poses significant challenges to engineers as their behavior can vary from solid-like to fluid-like. The quasi-static axisymmetric collapse of granular columns is an interesting case of granular flow relevant to the active deformation of retaining structures. This study aims to shed light on this type of collapse using numerical simulations including the smoothed particle hydrodynamics (SPH) method in the continuum framework and the discrete element method (DEM) in the discrete framework. Three-dimensional SPH and DEM models are developed to investigate the quasi-static axisymmetric collapse of granular columns, which is initiated by slowly expanding the cylindrical wall of the column. The Drucker–Prager constitutive model with non-associated flow rule is implemented into SPH formulations to model elastoplastic behavior of granular material. The linear contact model is used in DEM simulations with calibrated contact stiffness parameters. The effects of particle shape and hysteretic contact behavior in the DEM model are indirectly considered through reducing particle rotational velocities by a calibrated constant factor at every time step. Using the developed models, the final deposit profile (e.g., height, runout distance, repose angle), non-deformed region, collapse pattern and energy dissipation for different initial aspect ratios are investigated. In addition, final height, runout distance and energy dissipation are theoretically derived and are found to be in good agreement with SPH and DEM simulations. Results show that the quasi-static axisymmetric collapse can be modeled using both the continuum and discrete frameworks as long as appropriate constitutive models or contact behaviors are implemented. However, DEM simulations are more capable of capturing particle-level behaviors, such as the sharp edges of the free surface and flow front. Results demonstrate that quasi-static granular collapse is qualitatively similar to the dynamic collapse, but results in smaller runout distances for all aspect ratios.

UR - http://www.scopus.com/inward/record.url?scp=85053046685&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85053046685&partnerID=8YFLogxK

U2 - 10.1007/s11440-018-0707-9

DO - 10.1007/s11440-018-0707-9

M3 - Article

AN - SCOPUS:85053046685

JO - Acta Geotechnica

JF - Acta Geotechnica

SN - 1861-1125

ER -