Aller au menu Aller au contenu
Science et ingénierie des matériaux et des procédés

> SIMAP_Groupes > SIMAP_GPM2

Discrete simulations (page under construction)

Although Finite Element simultions are routinely used at GPM2, implementing Discrete simulations is a specific feature of the group. In-house codes have been developped for several years based on Discrete Dislocation Dynamics (TRIDIS, NUMODIS) or Discrete Element Methods (dp3D). These codes provide powerful tools for investigating plasticity, damage and fracture mechanisms. Smoothed Particles Hydrodynamics (SPH) simulations are also implemented for the study of cavitation erosion.


Multiscale simulations of nanoindentation

Our simulations of nanoindentation in copper are emblematic of our multiscale approach. In this study,  Molecular Dynamics simulations provide the details of the nucleation mechanism beneath a spherical indenter. The mechanism is then introduced as local rules in DD simulations which can simulate the dynamic expansion of the dislocation microstructure as well as the surface relief (including pile-up and sink-in zones). At the continuum scale, DD simulations are used to define physical constitutive equations for FE crystal plasticity simulations which are performed and directly compared to AFM observations.

Crack initiation and propagation in fatigue

DDD simulations were able to successfully simulate the formation of persistent slip bands in surface grains loaded in fatigue. Extrusions are evidenced where the bands intercept the free surface. Energy and stress calculations performed inside the simulated grain lead to a possible scenario for the crack initiation at the interface between the band and the matrix, as reported in the literature. More recently, a crack was inserted at the persistent slip band interface and the crack tip slip displacement evolutions are evaluated. It is shown that the crack growth rate is strongly related to the grain size and to the distance to the grain boundary; the smaller the grain, the faster the crack growth. Finally, the crack propagation to the next grain is investigated by conducting DDD fatigue simulations in a surface grain next to a cracked grain. It is shown that the development of the persistent slip band is modified by the presence of the crack. The crack orientation affects the orientation of the persistent slip band, as well as the extrusion rate, and consequently the crack propagation in the next grain.


3D simulations of polycrystalline aggregates


In collaboration with Université Catholique de Louvain-La-Neuve, we are working on the role of interfaces in plasticity. The effect of the grain shape on the back stress induced in a single closed grain has been investigated in order to implement the microscopic features in a continuum description of the plasticity. Slip transmission and dislocation interactions between different grains of polycrystalline aggregates are also investigated for different orientations of the grains.


3D simulations of creep in Ni superalloy

The creep behavior of Ni superalloys is studied in collaboration with Bochum University. For these simulations the grain shape is determined by phase field computations and the lattice mismatch is calculated using a FFT procedure.

Mass loss prediction in cavitation


Cavitation erosion occurs when multiple bubbles collapse near a solid surface. During the bubble collapse, a fast reintrant jet impacts the surface and creates a deformation pit.

An inverse method has been conducted in order to estimate the aggressiveness of the flow from the observation of the pits printed on the surface in the first moments of the cavitation erosion. Assuming that each pit was generated by a single bubble collapse whose pressure field is treated as a Gaussian shape, finite element calculations are used for estimating the load that created each residual imprint. It is shown that the load distribution falls on a master curve independent of the tested material; the softer material (aluminum alloy) measuring the lowest impacts while the most resistant material (duplex stainless steel) provides access to the largest impact pressures.


Fluid/structure interaction for cavitation erosion

Mass loss prediction in cavitation erosion could be predicted by numerical simulations provided the coupling between the cavitating fluid and the elasto-visco-plastic solid is accurately taken into account. We selected two approaches to tackle this problem. In the first approach, the Finite Element code CAST3M is coupled to a CFD code developed at Michigan University. In the second approach Smoothed Particules Hydrodynamics (SPH) simulations are developed in order to account both for the fluid part and the solid part in the same code.

Discrete Element simulations

DEM simulations are used to simulate fracture in porous and dense ceramics, based on our in-house code dp3D. The aim is to understand the link between fracture, toughness and microstructure to design strain tolerant architectures.
Crack pattern evolution in a ceramic under thermal shock​​​​​​​​​​​​​​

Crack pattern evolution in a dense ceramic submitted to thermal shock.
Crack evolution in a alumina nacre
Crack evolution in a notched alumina nacre. The initial microstructure originates from an EBSD image.
Scaling law between toughness and Young's modulus normalized by the number of bonds per particle in a porous ceramic
Scaling law relating toughness to the Young's modulus normalized by the number of bonds per particle in a porous ceramic.

mise à jour le 9 octobre 2019

PhD and post-docs

Hareesh TUMMALA,
Siwen GAO (Bochum Univ., Germany),
Gururaj KADIRI (IGCAR, India),
Samir Chandra ROY, Ph.D. Univ. Grenoble Alpes, December 11th 2015,
Tiana DEPLANCKE (Post Doc 2015),
Shrey JOSHI,
Prasantha SARKAR,
Markku YLONEN,
Fabien MENG,
Jessica MUZY,
Corentin VELARD
Discrete element simulations:
Julia BONALDO (2021)
​​​​​​​Alexis BURR (2017)
Robin GIBAUD (2017)
​​​​​​​Etienne GUESNET (2018)
​​​​​​​Guillaume HAMELIN (2021)
Paul PARANT (2016)
Brayan PAREDES (2022)
Kaoutar RADI (2019)
Bruno SALQUE (2017)


Christophe DEPRES (SYMME, Annecy)
Christian ROBERTSON and Laurent DUPUY (SRMA, CEA Saclay)
Hyung-Jun CHANG (SafranTech)
Laurent DELANNAY and Thomas PARDOEN (UCL, Louvain-La-Neuve)
Anxin MA and Alexander HARTMAIER (ICAMS, Univ. Bochum, Germany)
Jean-Pierre FRANC, Giovanni GHIGLIOTTI (LEGI, Grenoble)
Christian PELLONE (LEGI, Grenoble)
Nicolas RANC (PIMM, ENSAM Paristech)
Khi Han KIM (ONR, USA)
Erik JOHNSEN (Michigan Univ, USA).
Univ. Grenoble Alpes