Improved convected particle domain interpolation method for coupled dynamic analysis of fully saturated porous media involving large deformation

Abstract Based on the u-p form governing equations and the convected particle domain interpolation technique, the improved convected particle domain interpolation based material point method (CCPDI) is developed in this paper for the coupled dynamic and contact analysis of fully saturated porous media involving large deformation. The numerical artifact noises due to material points crossing computational grid boundaries that usually occur in large deformation are eliminated by using the smoother interpolation functions presented in the convected particle domain interpolation technique. The discrete equations are derived in the framework of the generalized interpolation material point method. The boundary load tracking algorithm and a modified contact algorithm are proposed based on the definition of particle domains to apply accurately the surface loads on the motive boundaries and to capture correctly the contact time/behaviors in large deformation problems, respectively. Simulations of several representative one- and two-dimensional problems are presented to demonstrate the accuracy and effectiveness of the proposed methods. Compared to those obtained by using the coupling material point method and/or finite element method, the simulation results illustrate that the proposed CCPDI method can be successfully used in simulating the large deformation coupled dynamic responses of the solid skeleton and fluid phase in fully saturated porous media and the large deformation impact between solid and saturated porous bodies.

[1]  Tadahiko Shiomi,et al.  Practical Programming in Computational Geomechanics: With Special Reference to Earthquake Engineering , 1999 .

[2]  R. Spilker,et al.  Finite element formulations for hyperelastic transversely isotropic biphasic soft tissues , 1998 .

[3]  D. Sulsky,et al.  A particle method for history-dependent materials , 1993 .

[4]  Howard L. Schreyer,et al.  Modeling delamination as a strong discontinuity with the material point method , 2002 .

[5]  Bernhard A. Schrefler,et al.  A FULLY COUPLED MODEL FOR WATER FLOW AND AIRFLOW IN DEFORMABLE POROUS MEDIA , 1993 .

[6]  Richard A. Regueiro,et al.  Dynamics of porous media at finite strain , 2004 .

[7]  J. Brackbill,et al.  The material-point method for granular materials , 2000 .

[8]  A. Sadeghirad,et al.  A convected particle domain interpolation technique to extend applicability of the material point method for problems involving massive deformations , 2011 .

[9]  Xiong Zhang,et al.  An explicit material point finite element method for hyper‐velocity impact , 2006 .

[10]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[11]  Hywel Rhys Thomas,et al.  On the development of a model of the thermo-mechanical-hydraulic behaviour of unsaturated soils , 1996 .

[12]  S. Bardenhagen,et al.  The Generalized Interpolation Material Point Method , 2004 .

[13]  O. Zienkiewicz,et al.  Dynamic behaviour of saturated porous media; The generalized Biot formulation and its numerical solution , 1984 .

[14]  J. Prévost.  Nonlinear transient phenomena in saturated porous media , 1982 .

[15]  Roberto Scotta,et al.  A fully coupled dynamic model for two-phase fluid flow in deformable porous media , 2001 .

[16]  Mark A Fleming,et al.  Meshless methods: An overview and recent developments , 1996 .

[17]  D. H. van Campen,et al.  A two-phase finite element model of the diastolic left ventricle. , 1991, Journal of biomechanics.

[18]  Pengfei Yang,et al.  Simulation of Fragmentation with Material Point Method Based on Gurson Model and Random Failure , 2012 .

[19]  O. Zienkiewicz,et al.  Large Strain Static and Dynamic Semisaturated Soil Behaviour , 1995 .

[20]  B R Simon,et al.  Structural models for human spinal motion segments based on a poroelastic view of the intervertebral disk. , 1985, Journal of biomechanical engineering.

[21]  Zhen Chen,et al.  Material point method for dynamic analysis of saturated porous media under external contact/impact of solid bodies , 2009 .

[22]  Vincent Lemiale,et al.  Material Point Method Simulation of Equal Channel Angular Pressing Involving Large Plastic Strain and Contact Through Sharp Corners , 2010 .

[23]  Xiong Zhang,et al.  An Object-Oriented MPM Framework for Simulation of Large Deformation and Contact of Numerous Grains , 2010 .

[24]  James E. Guilkey,et al.  A weighted least squares particle‐in‐cell method for solid mechanics , 2011 .

[25]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[26]  D. Sulsky,et al.  THE MATERIAL POINT METHOD FOR SIMULATION OF THIN MEMBRANES , 1999 .

[27]  M. Biot.  Theory of Propagation of Elastic Waves in a Fluid‐Saturated Porous Solid. I. Low‐Frequency Range , 1956 .

[28]  Z. Więckowski.  The material point method in large strain engineering problems , 2004 .

[29]  O. C. Zienkiewicz,et al.  A numerical model for immiscible two‐phase fluid flow in a porous medium and its time domain solution , 1990 .

[30]  John A. Nairn,et al.  Material Point Method Calculations with Explicit Cracks , 2003 .

[31]  Xiong Zhang,et al.  Contact algorithms for the material point method in impact and penetration simulation , 2011 .

[32]  A F Mak,et al.  A biphasic poroelastic analysis of the flow dependent subcutaneous tissue pressure and compaction due to epidermal loadings: issues in pressure sore. , 1994, Journal of biomechanical engineering.