scholarly journals Fast and efficient MATLAB-based MPM solver (fMPMM-solver v1.0)

2020 ◽  
Author(s):  
Emmanuel Wyser ◽  
Michel Jaboyedoff ◽  
Yury Y. Podladchikov
Keyword(s):  
Cantilever Beam ◽  
Large Scale ◽  
Solid Mechanics ◽  
Large Deformations ◽  
Material Point Method ◽  
Material Point ◽  
Snow Avalanches ◽  
High Level Language ◽  
Numerical Accuracy ◽  
High Level

Abstract. In this contribution, we present an efficient MATLAB-based implementation of the material point method (MPM) and its most recent variants. MPM has gained popularity over the last decade, especially for problems in solid mechanics in which large deformations are involved, i.e., cantilever beam problems, granular collapses and even large-scale snow avalanches. Although its numerical accuracy is lower than that of the widely accepted finite element method (FEM), MPM has been proven useful in overcoming some of the limitations of FEM, such as excessive mesh distortions. We demonstrate that MATLAB is an efficient high-level language for MPM implementations that solve elasto-dynamic and elasto-plastic problems, such as the cantilever beam and granular collapses, respectively. We report a computational efficiency factor of 20 for a vectorized code compared to a classical iterative version. In addition, the numerical efficiency of the solver surpassed those of previously reported MPM implementations in Julia, ad minima 2.5 times faster under a similar computational architecture.

scholarly journals A fast and efficient MATLAB-based MPM solver: fMPMM-solver v1.1

2020 ◽  
Vol 13 (12) ◽  
pp. 6265-6284
Author(s):  
Emmanuel Wyser ◽  
Yury Alkhimenkov ◽  
Michel Jaboyedoff ◽  
Yury Y. Podladchikov
Keyword(s):  
Cantilever Beam ◽  
Large Scale ◽  
Solid Mechanics ◽  
Large Deformations ◽  
Material Point ◽  
Snow Avalanches ◽  
High Level Language ◽  
Numerical Accuracy ◽  
Computational Performance ◽  
High Level

Abstract. We present an efficient MATLAB-based implementation of the material point method (MPM) and its most recent variants. MPM has gained popularity over the last decade, especially for problems in solid mechanics in which large deformations are involved, such as cantilever beam problems, granular collapses and even large-scale snow avalanches. Although its numerical accuracy is lower than that of the widely accepted finite element method (FEM), MPM has proven useful for overcoming some of the limitations of FEM, such as excessive mesh distortions. We demonstrate that MATLAB is an efficient high-level language for MPM implementations that solve elasto-dynamic and elasto-plastic problems. We accelerate the MATLAB-based implementation of the MPM method by using the numerical techniques recently developed for FEM optimization in MATLAB. These techniques include vectorization, the use of native MATLAB functions and the maintenance of optimal RAM-to-cache communication, among others. We validate our in-house code with classical MPM benchmarks including (i) the elastic collapse of a column under its own weight; (ii) the elastic cantilever beam problem; and (iii) existing experimental and numerical results, i.e. granular collapses and slumping mechanics respectively. We report an improvement in performance by a factor of 28 for a vectorized code compared with a classical iterative version. The computational performance of the solver is at least 2.8 times greater than those of previously reported MPM implementations in Julia under a similar computational architecture.

Investigating erosion and entrainment of snow avalanches using the material point method

2021 ◽  
Author(s):  
Xingyue Li ◽  
Betty Sovilla ◽  
Camille Ligneau ◽  
Chenfanfu Jiang ◽  
Johan Gaume
Keyword(s):  
Snow Cover ◽  
Large Scale ◽  
Material Point Method ◽  
Material Point ◽  
Point Method ◽  
Snow Avalanches ◽  
Runout Distance ◽  
Gravity Driven ◽  
Erosion Mechanisms ◽  
Dynamics Models

<p>Erosion and entrainment are critical processes in gravity-driven mass flows like snow avalanches, as they can significantly change the flow mass and momentum and thus affect the flow dynamics. In snow avalanches, snow cover can be considerably eroded but only partially entrained into the flow. Differentiating erosion and entrainment gives more accurate prediction of the increased flow mass and offers information on eroded snow cover remaining on the slope, but is challenging in practice. This study investigates snow avalanche erosion and entrainment with the material point method, focusing on exploring various erosion mechanisms, differences in erosion and entrainment, and their possible influences on runout distance. By using different mechanical properties for the flowing snow, distinct erosion patterns are observed and the corresponding temporal evolutions of entrainment, erosion, and deposition in the erodible bed are examined. Erosion and entrainment require an appropriate combination of snow friction and cohesion of the bed. If cohesion and/or friction are too low, the bed will naturally be unstable. On the other hand, highly cohesive and frictional bed will prevent erosion. For intermediate values, erosion and entrainment can be notable, and the amount of eroded snow shows a clear negative correlation with snow friction and cohesion while the entrained snow does not demonstrate a strong tendency. Furthermore, the release and erodible bed lengths are varied to study their effect on erosion and entrainment propensity. It is found that the increase in the lengths of the release zone and erodible bed leads to more erosion and entrainment as expected, but not necessarily to a longer runout distance. In our simulations, the release and erodible bed lengths are positively and negatively correlated with the runout distance, respectively. This implies that the runout distance can have opposite trends with erosion and entrainment, which might be closely related to the energy change of the simulated avalanches from the outlet of the erodible bed to the final deposit. Our results shed more light into the erosion and entrainment mechanisms and may contribute to improve related parametrizations in large-scale avalanche dynamics models.</p>

scholarly journals A stabilized mixed implicit Material Point Method for non-linear incompressible solid mechanics

2018 ◽  
Vol 63 (6) ◽  
pp. 1243-1260 ◽  
Author(s):  
I. Iaconeta ◽  
A. Larese ◽  
R. Rossi ◽  
E. Oñate
Keyword(s):  
Solid Mechanics ◽  
Material Point Method ◽  
Material Point ◽  
Point Method ◽  
Non Linear

Modeling Nanocomposite Piezoresistive Response With Electromechanical Cohesive Zone Material Point Method

2016 ◽  
Author(s):  
Adarsh K. Chaurasia ◽  
Gary D. Seidel
Keyword(s):  
Cohesive Zone ◽  
Large Deformations ◽  
Material Point Method ◽  
Material Point ◽  
Volume Element ◽  
Point Method ◽  
Current Work ◽  
Secondary Field ◽  
Small Strains ◽  
The Matrix

In the current work, the Material Point Method (MPM) is extended to allow for interfacial discontinuities in problems with composite materials using cohesive zone (CZ) techniques. The proposed CZMPM is observed to result in smaller errors in the primary and secondary field variables, especially near the interface, for a given boundary value problem in comparison to the traditional MPM solution. The proposed CZMPM is used to solve an electromechanical test problem with a single fiber in the matrix medium. It is observed that the proposed CZMPM results in smaller local and volume averaged errors. The CZMPM is further used to evaluate the effective piezoresistive response of the nanoscale carbon nanotube (CNT)-polymer composite with electron hopping in between the nanotubes. The observed effective piezoresistive response exhibits features similar to those reported in the literature using finite element techniques for small strains. However, CZMPM allows for large deformations of the nanoscale representative volume element as presented in the current work.

scholarly journals MPM and ALE Simulations of Large Deformations Geotechnics Instability Problems

DYNA ◽  
2020 ◽  
Vol 87 (212) ◽  
pp. 226-235
Author(s):  
Giovanna Monique Alelvan ◽  
Daniela Toro Rojas ◽  
Amanda Cristina Pedron Rossato ◽  
Raydel Lorenzo Reinaldo ◽  
Manoel Porfirio Cordão Neto
Keyword(s):  
Large Deformations ◽  
Material Point Method ◽  
Material Point ◽  
Large Displacements ◽  
Arbitrary Lagrangian Eulerian ◽  
Inclined Plane ◽  
Ale Method ◽  
Advantages And Disadvantages ◽  
Comparison Of The Results ◽  
The Continuum

Problems involving large deformations are the focus of numerical modeling researches in recent decades due to the challenge of finding a kinematic appropriate description of the continuum. In recent years, different formulations have been used to describe such problems as the Arbitrary Lagrangian Eulerian (ALE) method and the Material Point Method (MPM). These two methods allow to perform dynamic analyzes involving large deformations. In this way, this work aims to present a comparison of problems applied to Geotechnics involving large deformations and large displacements, using MPM and FEM associated with the ALE method. For this purpose, three problems are simulated: sliding of blocks on an inclined plane, runout process of sand and instability of a slope using the MPM and the FEM associated with the ALE method. In all cases a comparison of the results is presented, and the advantages and disadvantages of each method are discussed.

scholarly journals Run-Out Simulation of a Landslide Triggered by an Increase in the Groundwater Level Using the Material Point Method

Water ◽  
2020 ◽  
Vol 12 (10) ◽  
pp. 2817
Author(s):  
Antonello Troncone ◽  
Luigi Pugliese ◽  
Enrico Conte
Keyword(s):  
Groundwater Level ◽  
Large Deformations ◽  
Material Point Method ◽  
Material Point ◽  
Point Method ◽  
Deformation Mechanisms ◽  
The Finite Element Method ◽  
Numerical Techniques ◽  
Rainy Period ◽  
Small Deformations

Deformation mechanisms of the slopes are commonly schematized in four different stages: pre-failure, failure, post-failure and eventual reactivation. Traditional numerical methods, such as the finite element method and the finite difference method, are commonly employed to analyse the slope response in the pre-failure and failure stages under the assumption of small deformations. On the other hand, these methods are generally unsuitable for simulating the post-failure behaviour due to the occurrence of large deformations that often characterize this stage. The material point method (MPM) is one of the available numerical techniques capable of overcoming this limitation. In this paper, MPM is employed to analyse the post-failure stage of a landslide that occurred at Cook Lake (WY, USA) in 1997, after a long rainy period. Accuracy of the method is assessed by comparing the final geometry of the displaced material detected just after the event, to that provided by the numerical simulation. A satisfactory agreement is obtained between prediction and observation when an increase in the groundwater level due to rainfall is accounted for in the analysis.

Sign in / Sign up

Export Citation Format

Share Document