A finite micro-rotation material point method for micropolar solid and fluid dynamics with three-dimensional evolving contacts and free surfaces

2022 ◽  
Vol 391 ◽  
pp. 114540
Author(s):  
Ran Ma ◽  
WaiChing Sun
2021 ◽  
Vol 112 ◽  
pp. 103904
Author(s):  
Fabricio Fernández ◽  
Jhonatan E.G. Rojas ◽  
Eurípedes A. Vargas ◽  
Raquel Q. Velloso ◽  
Daniel Dias

2021 ◽  
Vol 14 (12) ◽  
pp. 7749-7774
Author(s):  
Emmanuel Wyser ◽  
Yury Alkhimenkov ◽  
Michel Jaboyedoff ◽  
Yury Y. Podladchikov

Abstract. We propose an explicit GPU-based solver within the material point method (MPM) framework using graphics processing units (GPUs) to resolve elastoplastic problems under two- and three-dimensional configurations (i.e. granular collapses and slumping mechanics). Modern GPU architectures, including Ampere, Turing and Volta, provide a computational framework that is well suited to the locality of the material point method in view of high-performance computing. For intense and non-local computational aspects (i.e. the back-and-forth mapping between the nodes of the background mesh and the material points), we use straightforward atomic operations (the scattering paradigm). We select the generalized interpolation material point method (GIMPM) to resolve the cell-crossing error, which typically arises in the original MPM, because of the C0 continuity of the linear basis function. We validate our GPU-based in-house solver by comparing numerical results for granular collapses with the available experimental data sets. Good agreement is found between the numerical results and experimental results for the free surface and failure surface. We further evaluate the performance of our GPU-based implementation for the three-dimensional elastoplastic slumping mechanics problem. We report (i) a maximum 200-fold performance gain between a CPU- and a single-GPU-based implementation, provided that (ii) the hardware limit (i.e. the peak memory bandwidth) of the device is reached. Furthermore, our multi-GPU implementation can resolve models with nearly a billion material points. We finally showcase an application to slumping mechanics and demonstrate the importance of a three-dimensional configuration coupled with heterogeneous properties to resolve complex material behaviour.


2019 ◽  
Vol 56 (9) ◽  
pp. 1318-1337 ◽  
Author(s):  
Xiaorong Xu ◽  
Feng Jin ◽  
Qicheng Sun ◽  
Kenichi Soga ◽  
Gordon G.D. Zhou

This study presents a field-scale simulation of the Hongshiyan landslide in China. It uses an advanced numerical approach (material point method (MPM)) and a constitutive model (the Drucker–Prager model + μ(I) rheological relation) for the three-dimensional (3D) simulation. The performance of the developed MPM model is validated with laboratory-scale experimental data on granular collapse before being applied to field-scale analyses. ArcGIS data are used to create a 3D MPM model of the soil body with complicated geometry. Although the developed model can describe the multiple phases of granular flow, it focuses on the runout behavior of the landslide in this work. The landslide is assumed to have occurred suddenly due to an earthquake, and global sudden failure rather than progressive failure is modeled. The MPM simulation results match reasonably well with the measured post-earthquake topography (e.g., deposit height of about 120 m and stretch length of about 900 m in the river) and landslide duration of about 1 min. The velocity of the sliding mass increases rapidly during flow, especially in the first 20 s. The velocity profiles along the depth direction at different locations of the sliding body exhibit an exponential distribution similar to that of a Bagnold-type profile, indicating that the sliding body is fully mobilized. The rate-dependent dissipation parameter β used in the model significantly influences the runout behavior (e.g., flow speed, velocity distribution, and deposit shape).


2019 ◽  
Vol 65 (2) ◽  
pp. 555-581 ◽  
Author(s):  
José Leόn González Acosta ◽  
Philip J. Vardon ◽  
Guido Remmerswaal ◽  
Michael A. Hicks

AbstractStress inaccuracies (oscillations) are one of the main problems in the material point method (MPM), especially when advanced constitutive models are used. The origins of such oscillations are a combination of poor force and stiffness integration, stress recovery inaccuracies, and cell crossing problems. These are caused mainly by the use of shape function gradients and the use of material points for integration in MPM. The most common techniques developed to reduce stress oscillations consider adapting the shape function gradients so that they are continuous at the nodes. These techniques improve MPM, but problems remain, particularly in two and three dimensional cases. In this paper, the stress inaccuracies are investigated in detail, with particular reference to an implicit time integration scheme. Three modifications to MPM are implemented, and together these are able to remove almost all of the observed oscillations.


2019 ◽  
Vol 17 (04) ◽  
pp. 1950003
Author(s):  
Jun Tao ◽  
Yonggang Zheng ◽  
Hongwu Zhang ◽  
Zhen Chen

In this paper, an axisymmetric generalized interpolation material point method for fully coupled thermomechanical analysis (AxiCTGIMP) is developed for evaluating the transient responses, where both the thermoelastic and thermoplastic effects are taken into account. The generalized interpolation material point method (GIMP) discretization in space for the coupled governing equations is described in detail. A staggered solution scheme is designed to split the coupled system into the parts related to the temperature and displacement fields, respectively, which are then solved individually with explicit time integration. The AxiCTGIMP is then verified and validated with two benchmark examples: the thick-walled cylinder and the Taylor-bar impact test. The simulation results show good agreements with available analytical solutions, experimental data and other numerical results. In addition, the results indicate that the proposed solution procedure is more accurate than the original MPM while it is much more efficient than the fully three-dimensional simulation for the axisymmetric thermomechanical problems.


2010 ◽  
Vol 9 (2 and 3) ◽  
Author(s):  
M Connolly ◽  
E Maldonado ◽  
M Roth

The two-dimensional Material Point Method (MPM) algorithm outlined by Chen and Brannon has been extended to three dimensions. The development of the code is discussed as well as applications for simulating bullet impact on biological and non-biological systems.


2008 ◽  
Vol 49 (5) ◽  
pp. 1129-1138 ◽  
Author(s):  
Pan Xiao-Fei ◽  
Xu Ai-Guo ◽  
Zhang Guang-Cai ◽  
Zhang Ping ◽  
Zhu Jian-Shi ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document