Nikhil Madan, PhD


Doctoral thesis
2019-11-28New formulation of the discrete element method with deformable particles 
supervisor -- Prof. Jerzy Rojek, PhD, DSc, IPPT PAN
supervisor -- Szymon Nosewicz, PhD, IPPT PAN
1413
 
Recent publications
1.Madan N., Rojek J., Nosewicz S., Convergence and stability analysis of the deformable discrete element method, INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, ISSN: 0029-5981, DOI: 10.1002/nme.6014, Vol.118, No.6, pp.320-344, 2019
Abstract:

This work investigates numerical properties of the algorithm of the discrete element method employing deformable circular discs presented in an earlier authors' publication. The new formulation, called the deformable discrete element method (DDEM) enhances the standard discrete element method (DEM) by introducing an additional (global) deformation mode caused by the stresses in the particles induced by the contact forces. An accurate computation of the contact forces would require an iterative solution of the implicit relationship between the contact forces and particle displacements. In order to preserve efficiency of the DEM, the new formulation has been adapted to the explicit time integration. It has been shown that the explicit DDEM algorithm is conditionally stable and there are two restrictions on its stability. Except for the limitation of the time step as in the standard DEM, the stability in the DDEM is governed by the convergence criterion of the iterative solution of the contact forces. The convergence and stability limits have been determined analytically and numerically for selected regular and irregular configurations. It has also been found out that the critical time step in DDEM remains unchanged with respect to standard DEM.

Keywords:

discrete element method, deformable particles, iterative solution, convergence criterion, explicit scheme, stability

Affiliations:
Madan N.-IPPT PAN
Rojek J.-IPPT PAN
Nosewicz S.-IPPT PAN
2.Rojek J., Madan N., Nosewicz S., Micro–macro relationships in the simulation of wave propagation phenomenon using the discrete element method, Materials, ISSN: 1996-1944, DOI: 10.3390/ma12244241, Vol.12, No.24, pp.4241-1-22, 2019
Abstract:

The present work is aimed to investigate the capability of the discrete element method (DEM) to model properly wave propagation in solid materials, with special focus on the determination of elastic properties through wave velocities. Reference micro–macro relationships for elastic constitutive parameters have been based on the kinematic hypothesis as well as obtained numerically by simulation of a quasistatic uniaxial compression test. The validity of these relationships in the dynamic analysis of the wave propagation has been checked. Propagation of the longitudinal and shear wave pulse in rectangular sample discretized with discs has been analysed. Wave propagation velocities obtained in the analysis have been used to determine elastic properties. Elastic properties obtained in the dynamic analysis have been compared with those determined by simulation of the quasistatic compression test.

Keywords:

discrete element method, wave propagation, elastic properties, micro–macro relationships

Affiliations:
Rojek J.-IPPT PAN
Madan N.-IPPT PAN
Nosewicz S.-IPPT PAN
3.Rojek J., Zubelewicz A., Madan N., Nosewicz S., The discrete element method with deformable particles, INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, ISSN: 0029-5981, DOI: 10.1002/nme.5767, Vol.114, No.8, pp.828-860, 2018
Abstract:

This work presents a new original formulation of the discrete element method (DEM) with deformable cylindrical particles. Uniform stress and strain fields are assumed to be induced in the particles under the action of contact forces. Particle deformation obtained by strain integration is taken into account in the evaluation of interparticle contact forces. The deformability of a particle yields a nonlocal contact model, it leads to the formation of new contacts, it changes the distribution of contact forces in the particle assembly, and it affects the macroscopic response of the particulate material. A numerical algorithm for the deformable DEM (DDEM) has been developed and implemented in the DEM program DEMPack. The new formulation implies only small modifications of the standard DEM algorithm. The DDEM algorithm has been verified on simple examples of an unconfined uniaxial compression of a rectangular specimen discretized with regularly spaced equal bonded particles and a square specimen represented with an irregular configuration of nonuniform-sized bonded particles. The numerical results have been verified by a comparison with equivalent finite elementmethod results and available analytical solutions. The micro-macro relationships for elastic parameters have been obtained. The results have proved to have enhanced the modeling capabilities of the DDEM with respect to the standard DEM.

Keywords:

average stress, deformable particles, discrete element method, elastic constants, micro-macro relationships, nonlocal contact model

Affiliations:
Rojek J.-IPPT PAN
Zubelewicz A.-University of New Mexico (US)
Madan N.-IPPT PAN
Nosewicz S.-IPPT PAN

Conference papers
1.Rojek J., Zubelewicz A., Madan N., Nosewicz S., New formulation of the discrete element method, CMM 2017, 22nd International Conference on Computer Methods in Mechanics, 2017-09-13/09-16, Lublin (PL), DOI: 10.1063/1.5019043, Vol.1922, pp.030009-1-8, 2018
Abstract:

A new original formulation of the discrete element method based on the soft contact approach is presented in this work. The standard DEM has heen enhanced by the introduction of the additional (global) deformation mode caused by the stresses in the particles induced by the contact forces. Uniform stresses and strains are assumed for each particle. The stresses are calculated from the contact forces. The strains are obtained using an inverse constitutive relationship. The strains allow us to obtain deformed particle shapes. The deformed shapes (ellipses) are taken into account in contact detection and evaluation of the contact forces. A simple example of a uniaxial compression of a rectangular specimen, discreti.zed with equal sized particles is simulated to verify the DDEM algorithm. The numerical example shows that a particle deformation changes the particle interaction and the distribution of forces in the discrete element assembly. A quantitative study of micro-macro elastic properties proves the enhanced capabilities of the DDEM as compared to standard DEM.

Affiliations:
Rojek J.-IPPT PAN
Zubelewicz A.-University of New Mexico (US)
Madan N.-IPPT PAN
Nosewicz S.-IPPT PAN

Conference abstracts
1.Rojek J., Madan N., Nosewicz S., A novel formulation of the discrete element method with deformable particles, 12HSTAM 2019, International Congress on Mechanics, 2019-09-22/09-25, Thessaloniki (GR), pp.59-59, 2019
2.Rojek J., Madan N., Nosewicz S., Enhanced modelling capabilities of the discrete element method with deformable particles, 8th International Conference on Discrete Element Methods, 2019-07-21/07-26, Enschede (NL), pp.181, 2019
3.Rojek J., Madan N., Nosewicz S., Simulation of elastic wave propagation using the deformable discrete element method, KomPlasTech 2019, Computer Methods in Materials Technology, 2019-01-13/01-16, Zakopane (PL), pp.106-107, 2019
Keywords:

wave propagation, elasticity, discrete element method,simulation

Affiliations:
Rojek J.-IPPT PAN
Madan N.-IPPT PAN
Nosewicz S.-IPPT PAN
4.Rojek J., Madan N., Nosewicz S., The discrete element method with deformable particles, PCM-CMM, 4th Polish Congress of Mechanics, 23rd International Conference on Computer Methods in Mechanics, 2019-09-08/09-12, Kraków (PL), pp.1-1, 2019
Keywords:

Discrete Element Method, Deformable Particles, Macroscopic Properties

Affiliations:
Rojek J.-IPPT PAN
Madan N.-IPPT PAN
Nosewicz S.-IPPT PAN
5.Madan N., Rojek J., Nosewicz S., Enhanced wave propagation modelling capabilities of discrete element method using deformable elements, PCM-CMM, 4th Polish Congress of Mechanics, 23rd International Conference on Computer Methods in Mechanics, 2019-09-08/09-12, Kraków (PL), pp.1-1, 2019
Keywords:

Elastic Wave Propagation, Deformability, Discrete Element Method

Affiliations:
Madan N.-IPPT PAN
Rojek J.-IPPT PAN
Nosewicz S.-IPPT PAN
6.Madan N., Rojek J., Nosewicz S., The deformable discrete element method - formulation and application, YIC2019, 5th ECCOMAS Young Investigators Conference, 2019-09-01/09-06, Kraków (PL), pp.1-2, 2019
7.Rojek J., Madan N., Nosewicz S., Zubelewicz A., The deformable discrete element method, 6th European Conference on Computational Mechanics (ECCM 6), 7th European Conference on Computational Fluid Dynamics (ECFD 7), 2018-06-11/06-15, Glasgow (GB), pp.1, 2018
Keywords:

Discrete Element Method, Deformable Particles, Nonlocal Contact Model, Poisson's Effect

Affiliations:
Rojek J.-IPPT PAN
Madan N.-IPPT PAN
Nosewicz S.-IPPT PAN
Zubelewicz A.-University of New Mexico (US)
8.Madan N., Rojek J., Zubelewicz A., Nosewicz S., Convergence limit of a deformable discrete element model, SolMech 2018, 41st SOLID MECHANICS CONFERENCE, 2018-08-27/08-31, Warszawa (PL), pp.204-205, 2018
9.Rojek J., Zubelewicz A., Madan N., Nosewicz S., Lumelskyj D., A novel treatment for the deformability of discrete elements, SolMech 2018, 41st SOLID MECHANICS CONFERENCE, 2018-08-27/08-31, Warszawa (PL), pp.202-203, 2018
10.Rojek J., Zubelewicz A., Madan N., Nosewicz S., New formulation of the discrete element method, CMM 2017, 22nd International Conference on Computer Methods in Mechanics, 2017-09-13/09-16, Lublin (PL), pp.MS13-27-28, 2017
Abstract:

This work presents a new original formulation of the discrete element method based on the soft contact approach. The standard DEM has been enhanced by introduction of the additional (global) deformation mode caused by the stresses in the particles induced by the contact forces. Uniform stresses and strains are assumed for each particle. The stresses are calculated from the contact forces. The strains are obtained using an inverse constitutive relationship. The strains allow us to obtain deformed particle shapes. The deformed shapes (ellipses) are taken into account in contact detection and evaluation of the contact forces. The numerical example shows that a particle deformation changes the particle interaction and the distribution of forces in the discrete element assembly.

Keywords:

discrete element method, deformable particles, soft contact

Affiliations:
Rojek J.-IPPT PAN
Zubelewicz A.-University of New Mexico (US)
Madan N.-IPPT PAN
Nosewicz S.-IPPT PAN