The material point method for simulation of thin membranes

✍ Scribed by Allen R. York II; Deborah Sulsky; Howard L. Schreyer

• Publisher: John Wiley and Sons
• Year: 1999
• Tongue: English
• Weight: 581 KB
• Volume: 44
• Category: Article
• ISSN: 0029-5981
• DOI: 10.1002/(sici)1097-0207(19990410)44:10<1429::aid-nme536>3.0.co;2-4

✦ Synopsis

The material-point method (MPM) is extended to handle membranes, which are discretized by a collection of unconnected material points placed along each membrane surface. These points provide a Lagrangian description of the membrane. To solve for the membrane motion, data carried by the material points are transferred to a background mesh where equations of motion are discretized and solved. Then the solution on the background mesh is used to update the membrane material points. This process of combining Lagrangian and Eulerian features is standard in MPM; the modification for membranes involves merely an implementation of the constitutive equation in a local, normal-tangential coordinate system. It is shown that this procedure does, in fact, provide adequate resolution of membranes with thicknesses that can vary substantially from that of the background mesh spacing. A general formulation is given, but the implementation is in a two-dimensional code that provides a proof-of-principle.

Numerical examples including a spring, pendulum and a string with initial slack are used to illustrate the method. The string with slack uses an additional modification of the membrane constitutive equation that allows wrinkles to be modeled at low computational cost. Presented also are examples of two disks impacting, pinching a membrane and rebounding, a difficult problem for standard finite element codes. These simulations require a relaxation of the automatic no-slip contact algorithm in MPM. The addition of the capability to model membranes and the new contact algorithm provide a significant improvement over existing methods for handling an important class of problems.