Yan Liu¹, Xiong Zhang¹, K. Y. Sze², Min Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.20, No.3, pp. 177-192, 2007, DOI:10.3970/cmes.2007.020.177

Abstract In molecular simulations, the frequencies of the low-frequency modes are many orders of magnitude lower than those of the high-frequency modes. Compared with the amplitudes of the low-frequency modes, the amplitudes of the high-frequency modes are often negligible and, thus, least interesting. As dictated by the period of the highest frequency mode, the critical time step for stable time integration can be significantly increased by suppressing the negligible high-frequency modes yet the solution remains virtually intact. In this light, a smoothed molecular dynamics (SMD) approach is proposed to eliminate the high-frequency modes from the dynamical system through the use of… More >

P. C. Wallstedt¹, J. E. Guilkey¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.3, pp. 223-232, 2007, DOI:10.3970/cmes.2007.019.223

Abstract The standard velocity projection scheme for the Material Point Method (MPM) and a typical form of the GIMP Method are examined. It is demonstrated that the fidelity of information transfer from a particle representation to the computational grid is strongly dependent on particle density and location. In addition, use of non-uniform grids and even non-uniform particle sizes are shown to introduce error. An enhancement to the projection operation is developed which makes use of already available velocity gradient information. This enhancement facilitates exact projection of linear functions and reduces the dependence of projection accuracy on particle location and density for… More >

J. Ma¹, H. Lu¹, B. Wang¹, S. Roy¹, R. Hornung², A. Wissink², R. Komanduri^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 135-152, 2005, DOI:10.3970/cmes.2005.008.135

Abstract In the simulation of a wide range of mechanics problems including impact/contact/penetration and fracture, the material point method (MPM), Sulsky, Zhou and Shreyer (1995), demonstrated its computational capabilities. To resolve alternating stress sign and instability problems associated with conventional MPM, Bardenhagen and Kober (2004) introduced recently the generalized interpolation material point (GIMP) method and implemented for one-dimensional simulations. In this paper we have extended GIMP to 2D and applied to simulate simple tension and indentation problems. For simulations spanning multiple length scales, based on the continuum mechanics approach, we present a parallel GIMP computational method using the Structured Adaptive Mesh… More >

Y. Guo¹, J. A. Nairn¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 295-308, 2004, DOI:10.3970/cmes.2004.006.295

Abstract The Material Point Method (MPM), which is a particle-based, meshless method that discretizes material bodies into a collection of material points (the particles), is a new method for numerical analysis of dynamic solid mechanics problems. Recently, MPM has been generalized to include dynamic stress analysis of structures with explicit cracks. This paper considers evaluation of crack-tip parameters, such as J-integral and stress intensity factors, from MPM calculations involving explicit cracks. Examples for both static and dynamic problems for pure modes I and II or mixed mode loading show that MPM works well for calculation of fracture parameters. The MPM results… More >

S.G. Bardenhagen¹, J.E. Guilkey², K.M. Roessig³, J.U. Brackbill⁴, W.M. Witzel⁵, J.C.Foster⁶

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.4, pp. 509-522, 2001, DOI:10.3970/cmes.2001.002.509

Abstract Contact between deformable bodies is a difficult problem in the analysis of engineering systems. A new approach to contact has been implemented using the Material Point Method for solid mechanics, Bardenhagen, Brackbill, and Sulsky (2000a). Here two improvements to the algorithm are described. The first is to include the normal traction in the contact logic to more appropriately determine the free separation criterion. The second is to provide numerical stability by scaling the contact impulse when computational grid information is suspect, a condition which can be expected to occur occasionally as material bodies move through the computational grid. The modifications… More >