This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit m...This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit material point method(iMPM).The global stiffness matrix is formed explicitly and the Newton-Raphson iterative method is used to solve the equilibrium equations.Where possible,the implementation procedure closely follows standard finite element method(FEM)approaches to allow easy conversion of other FEM codes.The generalized interpolation function is assigned to eliminate the inherent cell crossing noise within conventional MPM.For the first time,the B-bar approach is used to overcome volumetric locking in standard GIMP method for near-incompressible non-linear geomechanics.The proposed i GIMP was tested and compared with i MPM and analytical solutions via a 1 D column compression problem.Results highlighted the superiority of the i GIMP approach in reducing stress oscillations,thereby improving computational accuracy.Then,elasto-plastic slope stabilities and rigid footing problems were considered,further illustrating the ability of the proposed method to overcome volumetric locking due to incompressibility.Results showed that the proposed i GIMP with B-bar approach can be used to simulate geotechnical problems with large deformations.展开更多
Many materials such as biological tissues, polymers, and metals in plasticity can undergo large deformations with very little change in volume. Low-order finite elements are also preferred for certain applications, bu...Many materials such as biological tissues, polymers, and metals in plasticity can undergo large deformations with very little change in volume. Low-order finite elements are also preferred for certain applications, but are well known to behave poorly for such nearly incompressible materials. Of the several methods to relieve this volumetric locking, the method remains popular as no extra variables or nodes need to be added, making the implementation relatively straightforward and efficient. In the large deformation regime, the incompressibility is often treated by using a reduced order or averaged value of the volumetric part of the deformation gradient, and hence this technique is often termed an approach. However, there is little in the literature detailing the relationship between the choice of and the resulting and stiffness matrices. In this article, we develop a framework for relating the choice of to the resulting and stiffness matrices. We examine two volume-averaged choices for , one in the reference and one in the current configuration. Volume-averaged formulation has the advantage that no integration points are added. Therefore, there is a modest savings in memory and no integration point quantities needed to be interpolated between different sets of points. Numerical results show that the two formulations developed give similar results to existing methods.展开更多
基金the National Natural Science Foundation of China(Nos.41807223 and 51908175)the Fundamental Research Funds for the Central Universities(No.B210202096)+1 种基金the Natural Science Foundation of Guangdong Province(No.2018A030310346)the Water Conservancy Science and Technology Innovation Project of Guangdong Province(No.2020-11),China。
文摘This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit material point method(iMPM).The global stiffness matrix is formed explicitly and the Newton-Raphson iterative method is used to solve the equilibrium equations.Where possible,the implementation procedure closely follows standard finite element method(FEM)approaches to allow easy conversion of other FEM codes.The generalized interpolation function is assigned to eliminate the inherent cell crossing noise within conventional MPM.For the first time,the B-bar approach is used to overcome volumetric locking in standard GIMP method for near-incompressible non-linear geomechanics.The proposed i GIMP was tested and compared with i MPM and analytical solutions via a 1 D column compression problem.Results highlighted the superiority of the i GIMP approach in reducing stress oscillations,thereby improving computational accuracy.Then,elasto-plastic slope stabilities and rigid footing problems were considered,further illustrating the ability of the proposed method to overcome volumetric locking due to incompressibility.Results showed that the proposed i GIMP with B-bar approach can be used to simulate geotechnical problems with large deformations.
文摘Many materials such as biological tissues, polymers, and metals in plasticity can undergo large deformations with very little change in volume. Low-order finite elements are also preferred for certain applications, but are well known to behave poorly for such nearly incompressible materials. Of the several methods to relieve this volumetric locking, the method remains popular as no extra variables or nodes need to be added, making the implementation relatively straightforward and efficient. In the large deformation regime, the incompressibility is often treated by using a reduced order or averaged value of the volumetric part of the deformation gradient, and hence this technique is often termed an approach. However, there is little in the literature detailing the relationship between the choice of and the resulting and stiffness matrices. In this article, we develop a framework for relating the choice of to the resulting and stiffness matrices. We examine two volume-averaged choices for , one in the reference and one in the current configuration. Volume-averaged formulation has the advantage that no integration points are added. Therefore, there is a modest savings in memory and no integration point quantities needed to be interpolated between different sets of points. Numerical results show that the two formulations developed give similar results to existing methods.
