An efficient staggered grid material point method

被引:48
作者
Liang, Yong [1 ]
Zhang, Xiong [1 ]
Liu, Yan [1 ]
机构
[1] Tsinghua Univ, Sch Aerosp Engn, Beijing 100084, Peoples R China
基金
中国国家自然科学基金;
关键词
Material point method; Staggered grid; Cell crossing noise; Large deformation; IN-CELL METHOD; HYPERVELOCITY IMPACT; SIMULATION; MPM; BOUNDARY; CONTACT; FLIP;
D O I
10.1016/j.cma.2019.04.024
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
The material point method (MPM) has demonstrated itself as an effective numerical method to simulate extreme events with large deformations. However, the original MPM suffers the cell crossing noise because it takes the material points as integration points and employs the piecewise linear grid nodal shape functions whose gradient is discontinuous on the cell boundary. A number of techniques have been developed to alleviate the cell crossing noise. In this paper, a new staggered grid material point method (SGMP) is proposed to eliminate the cell crossing noise very efficiently. The volume integrals in the weak form are evaluated by cell center quadrature instead of particle quadrature as the sum of value of the integrand at each cell center of the background grid multiplied by the corresponding quadrature weight. The physical quantities and the quadrature weights at the cell centers are reconstructed efficiently based on an auxiliary grid, which is obtained by shifting the background grid half the side length of its cell in each direction Similar to the original MPM, both grids carry no permanent information and can be reset after each time step. In addition, the SGMP evaluates the constitutive equations at the particles, just like the original MPM, to readily model the history-dependent materials. To further reduce the cell crossing noise, a continuous strain rate/vorticity field is established based on the auxiliary grid, whose values are determined by the background grid velocity gradient. The strain rate/vorticity at each particle is interpolated from the auxiliary grid nodal values. Due to the overlap of the cell centers and the corresponding auxiliary grid nodes, a very efficient implementation is established in the SGMP. Numerical studies illustrate that the SGMP is capable of eliminating the cell crossing noise with little extra computational effort and the extra cost ratio reduces as the number of the grid cells or the particles increases. (C) 2019 Elsevier B.V. All rights reserved.
引用
收藏
页码:85 / 109
页数:25
相关论文
共 52 条
[1]  
Al-Kafaji I. K. J., 2013, THESIS
[2]  
[Anonymous], 2016, THESIS
[3]  
[Anonymous], 1955, LAMS1956 LOS AL SCI
[4]  
[Anonymous], 2002, COMP METH APPL MECH
[5]  
Bardenhagen SG, 2004, CMES-COMP MODEL ENG, V5, P477
[6]  
Benioff MR., 2005, Report to the President
[7]   B-spline based boundary conditions in the material point method [J].
Bing, Y. ;
Cortis, M. ;
Charlton, T. J. ;
Coombs, W. M. ;
Augarde, C. E. .
COMPUTERS & STRUCTURES, 2019, 212 :257-274
[8]   FLIP - A METHOD FOR ADAPTIVELY ZONED, PARTICLE-IN-CELL CALCULATIONS OF FLUID-FLOWS IN 2 DIMENSIONS [J].
BRACKBILL, JU ;
RUPPEL, HM .
JOURNAL OF COMPUTATIONAL PHYSICS, 1986, 65 (02) :314-343
[9]   FLIP - A LOW-DISSIPATION, PARTICLE-IN-CELL METHOD FOR FLUID-FLOW [J].
BRACKBILL, JU ;
KOTHE, DB ;
RUPPEL, HM .
COMPUTER PHYSICS COMMUNICATIONS, 1988, 48 (01) :25-38
[10]   iGIMP: An implicit generalised interpolation material point method for large deformations [J].
Charlton, T. J. ;
Coombs, W. M. ;
Augarde, C. E. .
COMPUTERS & STRUCTURES, 2017, 190 :108-125