Difference between revisions of "Contact Laws"
| (30 intermediate revisions by the same user not shown) | |||
| Line 2: | Line 2: | ||
== Introduction == | == Introduction == | ||
[[NairnMPM]] and [[OSParticulas]] implement contact physics on [[MPM Input Files#Using Explicit Cracks|crack surfaces]] and between materials in [[Multimaterial MPM|multimaterial mode]] to model [[Friction|friction]] or [[Imperfect Interfaces|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|frictional contact | [[NairnMPM]] and [[OSParticulas]] implement contact physics on [[MPM Input Files#Using Explicit Cracks|crack surfaces]] and between materials in [[Multimaterial MPM|multimaterial mode]] to model [[Friction|friction]] or [[Imperfect Interfaces|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|frictional contact laws]] and [[#Imperfect Interface Contact Laws|imperfect interface laws]]. These default crack contact or for material-material contact are selected by using the [[Friction#Friction on Explicit Cracks|ContactCracks]] or [[Friction#Friction in Multimaterial MPM|ContactMM]] commands. if needed, the default contact laws can be customized for each individual crack when [[Defining Cracks|defining a new crack]] or can be customized for each material pair by using the [[Common Material Properties#Contact Properties|Contact material property]]. This section documents all the possible contact laws. | ||
All contact laws are using a [[Material Command Block|<tt>Material</tt> 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. | All contact laws are defined by using a [[Material Command Block|<tt>Material</tt> 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 Laws == | ||
Frictional contact laws give the sliding shear traction, S<sub>slide</sub>, as a function of the normal traction, N, the contact area, A<sub>c</sub>, the relative sliding velocity, Δv,and | [[Friction|Frictional contact]] can be models with laws that give the sliding shear traction, S<sub>slide</sub>, as a function of the normal traction, N, the contact area, A<sub>c</sub>, the relative sliding velocity ''after'' imposition of contact law, Δv, and possibly other parameters, or: | ||
| | ||
| Line 18: | Line 18: | ||
<math>S_{resultant} = \min(S_{slide},S_{stick})</math> | <math>S_{resultant} = \min(S_{slide},S_{stick})</math> | ||
where S<sub>stick</sub> 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 S<sub>slide</sub>, the surfaces will slide with the given sliding traction. More details on friction in MPM can be found in Nairn | where S<sub>stick</sub> 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 S<sub>slide</sub>, the surfaces will slide with the given sliding traction. More details on friction in MPM can be found in Nairn ''et. al'' (2016).<ref>J.A. Nairn, S. G. Bardenhagen, and G. D. Smith (2017) "Generalized Contact and Improved Friction Heating the Material Point Method," Computational Particle Mechanics, in press [http://www.cof.orst.edu/cof/wse/faculty/Nairn/papers/MPMContactLawsSV.pdf (author preprint)].</ref> | ||
The available contact laws are listed in the following table. See each law to learn about the function used to determine S<sub>slide</sub> and about the parameters required to use that law. | The available frictional contact laws are listed in the following table. See each law to learn about the function used to determine S<sub>slide</sub> and about the parameters required to use that law. | ||
{| class="wikitable" | {| class="wikitable" | ||
| Line 30: | Line 30: | ||
| [[Coulomb Friction Law|CoulombFriction]] || align="center"| 61 || contact by simple Coulomb friction | | [[Coulomb Friction Law|CoulombFriction]] || align="center"| 61 || contact by simple Coulomb friction | ||
|- | |- | ||
| [[Adhesive Friction Law| AdhesiveFriction]] || align="center"| 63 || contact by Coulomb friction with adhesion | | [[Adhesive Friction Law| AdhesiveFriction]] || align="center"| 63 || contact by velocity-dependent Coulomb friction with adhesion | ||
|- | |||
| [[Liquid Wall Contact Law| LiquidContact]] || align="center"| 64 || from slip to stick contact for fluid flow | |||
|} | |} | ||
== Imperfect Interface Contact Laws == | == Imperfect Interface Contact Laws == | ||
[[Imperfect Interfaces|Imperfect interfaces]] can be modeled two ways — by using contact laws on [[MPM Input Files#Using Explicit Cracks|explicit cracks]] or by using them in [[MPM Input Files#Multimaterial MPM|multimaterial mode MPM]]. When contacting surfaces are modeling an imperfect interface, the normal and tangential tractions (T<sub>n</sub> and T<sub>t</sub>) are set to be functions of the normal and tangential displacement discontinuities ([u<sub>n</sub>] and [u<sub>t</sub>]). In other words: | |||
| |||
<math>T_n = f_n([u_n],[u_t])</math> | |||
| |||
<math>T_t = f_s([u_n],[u_t])</math> | |||
The various imperfect interface laws differ by their options for the traction functions (f<sub>n</sub> and f<sub>s</sub>). More details on imperfect interfaces in MPM can be found in Nairn (2007)<ref name="Nairn04">J. A. Nairn, "Numerical implementation of imperfect interfaces," ''Computational Materials Science'', '''40''', 525–536 (2007).</ref> and (2013).<ref name="Nair13"> | |||
J.A. Nairn, "Modeling Imperfect Interfaces in the Material Point Method using Multimaterial Methods," ''Computer Modeling in Eng. & Sci.'', '''92''', 271-299 (2013).</ref> | |||
Note that whenever using imperfect interface modeling with cracks or in multimaterial mode, you must use the displacements method whenever the interface law has non-zero tangential traction. In other words, the simulation must omit the [[Detecting Contact#ContactPosition Command|<tt>ContactPosition</tt> and/or <tt>ContactPositionCracks</tt> commands]]. | |||
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. The [[Nonlinear Imperfect Interface|NonlinearInterface]] described below actually implements a linear interface law, but it can be subclasses by programmers to [[Creating a Custom Imperfect Interface Law|create a custom imperfect interface material]]. | |||
{| class="wikitable" | |||
|- | |||
! Name !! Number !! Description | |||
|- | |||
| [[Ignore Contact Law|IgnoreContact]] || align="center"| 60 || ignore the interface or revert to single material mode with no interface | |||
|- | |||
| [[Linear Imperfect Interface|LinearInterface]] || align="center"| 62 || imperfect interface with tractions linear in interfacial separations | |||
|- | |||
| [[Nonlinear Imperfect Interface|NonlinearInterface]] || align="center"| 65 || imperfect interface with tractions as a function of interfacial separations | |||
|- | |||
| [[Debonding Interface|DebondingInterface]] || align="center"| 65 || imperfect interface with failure stresses for debonding | |||
|} | |||
== References == | == References == | ||
<references/> | <references/> | ||
Latest revision as of 15:01, 11 March 2020
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 after imposition of contact law, Δ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 et. al (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 |
| LiquidContact | 64 | from slip to stick contact for fluid flow |
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]
Note that whenever using imperfect interface modeling with cracks or in multimaterial mode, you must use the displacements method whenever the interface law has non-zero tangential traction. In other words, the simulation must omit the ContactPosition and/or ContactPositionCracks commands.
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. The NonlinearInterface described below actually implements a linear interface law, but it can be subclasses by programmers to create a custom imperfect interface material.
| Name | Number | Description |
|---|---|---|
| IgnoreContact | 60 | ignore the interface or revert to single material mode with no interface |
| LinearInterface | 62 | imperfect interface with tractions linear in interfacial separations |
| NonlinearInterface | 65 | imperfect interface with tractions as a function of interfacial separations |
| DebondingInterface | 65 | imperfect interface with failure stresses for debonding |
References
- ↑ J.A. Nairn, S. G. Bardenhagen, and G. D. Smith (2017) "Generalized Contact and Improved Friction Heating the Material Point Method," Computational Particle Mechanics, in press (author preprint).
- ↑ 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).