Contact Laws
Introduction
NairnMPM and OSParticulas implement contact physics on crack surfaces and between materials in multimaterial mode to model friction or imperfect interfaces. The contact mechanics is determined by selecting a contact law. The currently available contact laws are divided into two types — frictional contact laws and imperfect interface laws. These default crack contact or for material-material contact are selected by using the ContactCracks or ContactMM commands. if needed, the default contact laws can be customized for each individual crack when defining a new crack or can be customized for each material pair by using the Contact material property. This section documents all the possible contact laws.
All contact laws are defined by using a Material command block. Within that block all contact law properties are set using property commands. Refer to each contact law type to learn about its possible properties.
Frictional Contact Laws
Frictional contact can be models with laws that give the sliding shear traction, Sslide, as a function of the normal traction, N, the contact area, Ac, the relative sliding velocity, Δv,and possibly other parameters, or:
[math]\displaystyle{ S_{slide} = f(N,A_c,\Delta v,...) }[/math]
Given any frictional law, the shear traction applied at any node is given by:
[math]\displaystyle{ S_{resultant} = \min(S_{slide},S_{stick}) }[/math]
where Sstick is the shear traction needed for tangential motion of the two surface to move together (i.e., to stick). In other words, if the shear traction calculated for frictional sliding is greater than the traction required for surfaces to stick together, then the surface will stick. Once that sticking shear traction exceeds Sslide, the surfaces will slide with the given sliding traction. More details on friction in MPM can be found in Nairn and Smith (2016).[1]
The available frictional contact laws are listed in the following table. See each law to learn about the function used to determine Sslide and about the parameters required to use that law.
| Name | Number | Description |
|---|---|---|
| IgnoreContact | 60 | ignore contact or revert to single material mode |
| CoulombFriction | 61 | contact by simple Coulomb friction |
| AdhesiveFriction | 63 | contact by velocity-dependent Coulomb friction with adhesion |
Imperfect Interface Contact Laws
Imperfect interfaces can be modeled two ways — by using contact laws on explicit cracks or by using them in multimaterial mode MPM. When contacting surfaces are modeling an imperfect interface, the normal and tangential tractions (Tn and Tt) are set to be functions of the normal and tangential displacement discontinuities ([un] and [ut]). In other words:
[math]\displaystyle{ T_n = f_n([u_n],[u_t]) }[/math]
[math]\displaystyle{ T_t = f_s([u_n],[u_t]) }[/math]
The various imperfect interface laws differ by their options for the traction functions (fn and fs). More details on imperfect interfaces in MPM can be found in Nairn (2007)[2] and (2013).[3]
The available imperfect interface contact laws are listed in the following table. See each law to learn about the functions used to determine the interfacial tractions and about the parameters required to use that law.
| Name | Number | Description |
|---|---|---|
| IgnoreContact | 60 | ignore contact or revert to single material mode |
| LinearInterface | 62 | imperfect interface with tractions linear in interfacial separations |
References
- ↑ J.A. Nairn and G. S. Smith (2016) "Generalized Contact and Improved Friction Heating the Material Point Method," in preparation.
- ↑ J. A. Nairn, "Numerical implementation of imperfect interfaces," Computational Materials Science, 40, 525–536 (2007).
- ↑ J.A. Nairn, "Modeling Imperfect Interfaces in the Material Point Method using Multimaterial Methods," Computer Modeling in Eng. & Sci., 92, 271-299 (2013).