Simulating Ice-Structure Interaction With the Material Point Method

2017 ◽  
Author(s):  
Yaomei Wang ◽  
Biye Yang ◽  
Guiyong Zhang ◽  
Yichen Jiang ◽  
Zhi Zong
Keyword(s):  
Yield Criterion ◽  
Ice Sheet ◽  
Material Point Method ◽  
Friction Angle ◽  
Complex Problem ◽  
Material Point ◽  
Point Method ◽  
Structure Interaction ◽  
Future Work ◽  
Large Deformation Problems

The process of ice-structure interaction is a complex problem which is influenced by the properties of both ice and the structure. In this paper, the material point method (MPM) is introduced to simulate the interaction between an ice sheet and a cylinder structure. MPM is efficient in solving history dependent and large deformation problems and has shown advantage in hyper-velocity impact and landslide issues, etc.. The constitutive relation of ice is based on elasto-viscous-plastic model with the Drucker-Pragers yield criterion. Ice follows the Maxwell elasto-viscous model before the yield criterion is reached and fails when the plastic strain surpasses the failure strain. Meanwhile, the constitutive model used in this work considers the effect of the Young’s modulus, Poisson’s ratio, density, temperature, cohesive force and internal friction angle of ice. A series of simulations are conducted and the results are in accord with existing theories. According to the comparison, the influences of ice temperature and penetration speed of the structure on the global ice load are testified. The numerical tests have proven the feasibility of MPM in simulating the interaction between an ice sheet and a cylinder structure. Future work in ice-structure interaction problems with MPM is also discussed.

scholarly journals Postfailure Characterization of Shallow Landslides Using the Material Point Method

Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-25
Author(s):  
Thanh Son Nguyen ◽  
Kuo-Hsin Yang ◽  
Chia-Chun Ho ◽  
Feng-Chi Huang
Keyword(s):  
Slope Failure ◽  
Material Point Method ◽  
Friction Angle ◽  
Material Point ◽  
Shallow Landslides ◽  
Point Method ◽  
Essential Information ◽  
Runout Distance ◽  
Sensitivity Assessment ◽  
Kinematic Behavior

Although the mechanisms of slope failure caused by rising groundwater have been widely investigated, the kinematic behavior of landslides in the postfailure stage, which contains essential information for hazard mitigation and risk assessment, has not yet been fully studied. Thus, in this study, a series of numerical simulations using the material point method (MPM) were conducted to analyze the kinematic behavior and soil movement of shallow landslides (infinite slope problems). First, the proposed MPM formulation was validated in a full-scale landslide flume test. The simulated results of final slope profile, runout distance, deposit height, shear band development, slope displacement, and velocity accorded with the experimental results, suggesting that the MPM can quantitatively simulate large deformations. A parametric study of shallow slopes with various hydrological conditions and soil hydraulic and soil mechanical parameters was then performed to assess the influence of the aforementioned factors on landslide kinematics. The simulation results indicated that mechanical behavior at the slope toe is complex; the multiple plastic shear bands generated at the slope toe were due to a combination of shearing and compression. The deposition profile of the slopes was significantly influenced by all input parameters. Among the aforementioned parameters, soil cohesion, location of the groundwater table, and saturated soil permeability most greatly affected runout distance in the sensitivity assessment. Soil friction angle had a minor influence on the kinematic behavior of the slope.

Fast and efficient MPM solver for strain localization problems

2020 ◽  
Author(s):  
Antoine Guerin ◽  
Emmanuel Wyser ◽  
Yury Podladchikov ◽  
Michel Jaboyedoff
Keyword(s):  
Strain Localization ◽  
Strain Softening ◽  
Yield Criterion ◽  
Pure Shear ◽  
Shear Banding ◽  
Shear Bands ◽  
Material Point Method ◽  
Material Point ◽  
Point Method ◽  
Mapping Algorithm

<p><span><span>Strain localization problems, i.e., shearbandings, have received a lot of interest, especially when strain softening is disregarded from the elasto-plastic consistution relationship. Indeed, reproducing correctly oriented shear bands without softening allows to overcome the mesh-depenency problem. Our work focuses on a Material Point Method (MPM) implementation of strain localization to i) study the behavior of shear bands in order to ii) assess the capabilities of this quite recent numerical method. </span></span></p><p><span><span>To study strain localization and shear banding, we developped an efficient numercial Material Point Method (MPM) solver in Matlab, based on the Update Stress Last (USL) scheme enriched with the Generalized Interpolation Material Point (GIMP) variant, which fixes a major flaw of any MPM solver: the cell-crossing error due to discontinuous gradient of the basis functions. This home-made solver allows us to study strain localizations in either a fixed or continuously deforming continuum. The algorithm solves explicitly momentum equations in an updated lagrangian manner similarly to an explicit FEM solver. We therefore investigate the compression of an elasto-plastic domain under pure shear condition, thus reproducing the geometrical settings and pure shear conditions used in Duretz et al. (2018). Strain softening is disregarded since we do not want any mesh dependence within the solver. A Mohr-Coulomb yield criterion was selected and plasticity was computed by a return mapping algorithm, i.e., we did not use consistent tangent operator. Localization is triggered by a weaker circular inculsion in the center of the domain.. </span></span></p><p><span><span>Preliminary results demonstrates the suitability of the MPM solver to reproduce the correct shearbanding behavior under compression, for both static and dynamic meshes. The higher the resolution, the more accurate are the shear bands. Naturally, this implies future implementations of the solver in a GPU-accelerated environment to increase the numerical resolution. </span></span></p>

scholarly journals Runout evaluation of Oso landslide with the material point method

2019 ◽  
Vol 56 (9) ◽  
pp. 1304-1317 ◽  
Author(s):  
Alba Yerro ◽  
Kenichi Soga ◽  
Jonathan Bray
Keyword(s):  
Material Point Method ◽  
Material Point ◽  
Point Method ◽  
Key Factors ◽  
Long Runout ◽  
Runout Distance ◽  
Significant Damage ◽  
Large Deformation Problems ◽  
Landslide Runout ◽  
Landslide Movement

Long runout landslides can cause significant damage and represent one of the most important problems in geotechnical engineering. Understanding the mechanics of the landslide runout process is important for risk assessment and is challenging due to its complexities. This work examines the runout of the 22 March 2014 Oso, Washington, landslide. The Oso landslide is one of the worst landslide disasters in USA history with 43 fatalities. It occurred in multiple failure stages, involving several failure surfaces and significant soil softening, and travelled over 1 km across the valley. It initiated after a period of wet weather in an area prone to landslide movements. The triggering causes of the landslide movement are still under investigation. In this paper, the material point method is used to simulate the runout of the Oso landslide. This numerical tool is capable of modeling large deformation problems. It is used to investigate several hypothetical scenarios to identify key factors that contributed to the Oso landslide long runout distance.

Simulations of Fall Cone Test in Soil Mechanics Using the Material Point Method

2016 ◽  
Vol 846 ◽  
pp. 336-341 ◽  
Author(s):  
M.A. Llano-Serna ◽  
M.M. Farias ◽  
D.M. Pedroso ◽  
David J. Williams ◽  
D. Sheng
Keyword(s):  
Soil Mechanics ◽  
Material Point Method ◽  
Liquid Limit ◽  
Material Point ◽  
Point Method ◽  
Cone Penetration ◽  
The Relationship ◽  
Large Deformation Problems

The material point method is a particle-based method that uses a double Lagrangian-Eulerian discretisation. This approach has proved its functionality for the simulation of large deformation problems. Such problems are frequent in geotechnical engineering, more specifically those related to penetration during pile driving and conventional in situ tests such as the Cone Penetration Test. The shallow laboratory fall cone test is considered in this paper. This test is widely used for the determination of the liquid limit of clays, but it is also used to study the relationship between penetration (h) and the undrained shear strength of clays (su). Simulations are verified against laboratory vane shear tests and fall cone tests performed on samples of kaolin clay at different moisture contents. Calibrations using a simple penetration-strength (h-su) model are made based on a single coefficient named the cone factor (K). The numerical results closely match both the experimental data and analytical solutions available in the literature.

scholarly journals A scalable and modular Material Point Method for large-scale simulations

2019 ◽  
Author(s):  
Krishna Kumar ◽  
Jeffrey Salmond ◽  
Shyamini Kularathna ◽  
Christopher Wilkes ◽  
Ezra Yoanes Setiasabda ◽  
...  
Keyword(s):  
Large Deformation ◽  
Large Scale ◽  
Material Point Method ◽  
Material Point ◽  
Point Method ◽  
Background Mesh ◽  
Data Structures And Algorithms ◽  
Material Points ◽  
Large Scale Simulations ◽  
Large Deformation Problems

In this paper, we describe a new scalable and modular material point method (MPM) code developed for solving large-scale problems in continuum mechanics. The MPM is a hybrid Eulerian-Lagrangian approach, which uses both moving material points and computational nodes on a background mesh. The MPM has been successfully applied to solve large-deformation problems such as landslides, failure of slopes, concrete flows, etc. Solving these large-deformation problems result in the material points actively moving through the mesh. Developing an efficient parallelisation scheme for the MPM code requires dynamic load-balancing techniques for both the material points and the background mesh. This paper describes the data structures and algorithms employed to improve the performance and portability of the MPM code.

Sign in / Sign up

Export Citation Format

Share Document