Improved contact algorithm for the material point method and application to stress propagation in granular material

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 described preserve important properties of the original algorithm, namely conservation of momentum, and the use of global quantities which obviate the need for neighbor searches and result in the computational cost scaling linearly with the number of contacting bodies. The algorithm is demonstrated on several examples. Deformable body solutions compare favorably with several problems which, for rigid bodies, have analytical solutions. A much more demanding simulation of stress propagation through idealized granular material, for which high fidelity data has been obtained, is examined in detail. Excellent qualitative agreement is found for a variety of contact conditions. Important material parameters needed for more quantitative comparisons are identified.