SimpleContact3D: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
(Created page with '{{CommandManualMenu}} This command is used to construct a SimpleContact3D element object. {| | style="background:yellow; color:black; width:600px" | '''element SimpleContact3D...')
 
(Change terminology to retained/constrained)
 
(4 intermediate revisions by one other user not shown)
Line 4: Line 4:


{|  
{|  
| style="background:yellow; color:black; width:600px" | '''element SimpleContact3D $eleTag $iNode $jNode $kNode $lNode $sNode $LNode $matTag $gTol $fTol'''
| style="background:yellow; color:black; width:600px" | '''element SimpleContact3D $eleTag $iNode $jNode $kNode $lNode $cNode $lNode $matTag $gTol $fTol'''
|}
|}


Line 11: Line 11:
|  style="width:150px" | '''$eleTag ''' || unique integer tag identifying element object
|  style="width:150px" | '''$eleTag ''' || unique integer tag identifying element object
|-
|-
|  '''$iNode $jNode $kNode $lNode''' || master nodes
|  '''$iNode $jNode $kNode $lNode''' || retained nodes (-ndm 3 -ndf 3)
|-
|-
|  '''$sNode ''' || slave node
|  '''$cNode ''' || constrained node (-ndm 3 -ndf 3)
|-
|-
|  '''$LNode ''' || Lagrange multiplier node
|  '''$lNode ''' || Lagrange multiplier node (-ndm 3 -ndf 3)
|-
|-
|  '''$matTag ''' || unique integer tag associated with previously-defined nDMaterial object
|  '''$matTag ''' || unique integer tag associated with previously-defined nDMaterial object
Line 28: Line 28:




The SimpleContact3D element is a three-dimensional node-to-surface contact element which defines a frictional contact interface between two separate bodies.  The master nodes are the nodes which define a surface of a hexahedral element on the first body, and the slave node is a node from the second body.  The Lagrange multiplier node is required to enforce the contact condition.  This node should not be shared with any other element in the domain.  Information on the theory behind this element can be found in, e.g. Wriggers (2002).
The SimpleContact3D element is a three-dimensional node-to-surface contact element which defines a frictional contact interface between two separate bodies.  The retained nodes are the nodes which define a surface of a hexahedral element on the first body, and the constrained node is a node from the second body.  The Lagrange multiplier node is required to enforce the contact condition.  This node should not be shared with any other element in the domain.  Information on the theory behind this element can be found in, e.g. Wriggers (2002).


'''NOTE:'''
'''NOTE:'''
Line 34: Line 34:
# The SimpleContact3D element has been written to work exclusively with the [[ContactMaterial3D|ContactMaterial3D nDMaterial]] object.
# The SimpleContact3D element has been written to work exclusively with the [[ContactMaterial3D|ContactMaterial3D nDMaterial]] object.
# The valid recorder queries for this element are:
# The valid recorder queries for this element are:
## ''force'' - returns the contact force acting on the slave node in vector form.
## ''force'' - returns the contact force acting on the constrained node in vector form.
## ''frictionforce'' - returns the frictional force acting on the slave node in vector form.
## ''frictionforce'' - returns the frictional force acting on the constrained node in vector form.
## ''forcescalar'' - returns the scalar magnitudes of the single normal and two tangential contact forces.
## ''forcescalar'' - returns the scalar magnitudes of the single normal and two tangential contact forces.
## The SimpleContact3D elements are set to consider frictional behavior as a default, but the frictional state of the SimpleContact3D element can be changed from the input file using the [[setParameter]] command.  When updating, value of 0 corresponds to the frictionless condition, and a value of 1 signifies the inclusion of friction.  An example command for this update procedure is provided below
# The SimpleContact3D element works well in static and pseudo-static analysis situations. 
#In transient analysis, the presence of the contact constraints can effect the stability of commonly-used time integration methods in the HHT or Newmark family (e.g., Laursen, 2002).  For this reason, use of alternative time-integration methods which numerically damp spurious high frequency behavior may be required.  The [[TRBDF2]] integrator is an effective method for this purpose.  The Newmark integrator can also be effective with proper selection of the gamma and beta coefficients.  The trapezoidal rule, i.e., Newmark with gamma = 0.5 and beta = 0.25, is particularly prone to instability related to the contact constraints and is not recommended.




'''EXAMPLE:''' 


'''EXAMPLE:'''  SimpleContact3D element with tag 1, connectivity with nodes 1, 2, 3, 4, 5, 6, 7, and 8, and material with tag 1
SimpleContact3D element with tag 1, connectivity with nodes 1, 2, 3, 4, 5, 6, 7, and 8, material with tag 1, and gap and force tolerances of 1.0e-10.


  element SimpleContact3D 1  1 2 3 4 5 6 7 8  1  1.0e-10 1.0e-10
  element SimpleContact3D 1  1 2 3 4 5 6 7 8  1  1.0e-10 1.0e-10
Update all of the SimpleContact3D elements with tags between 10 and 20 to consider a frictionless interface
setParameter -value 0 -eleRange 10 20 friction
'''REFERENCES:'''
# Wriggers, P. (2002).  ''Computational Contact Mechanics.''  John Wiley & Sons, Ltd, West Sussex, England.
# Laursen, T. A. (2002).  ''Computational Contact and Impact Mechanics.''  Springer-Verlag, Berlin.


----
----
Line 49: Line 63:


----
----
References:
# Wriggers, P. (2002).  ''Computational Contact Mechanics.''  John Wilely & Sons, Ltd, West Sussex, England.

Latest revision as of 17:15, 13 June 2020




This command is used to construct a SimpleContact3D element object.

element SimpleContact3D $eleTag $iNode $jNode $kNode $lNode $cNode $lNode $matTag $gTol $fTol


$eleTag unique integer tag identifying element object
$iNode $jNode $kNode $lNode retained nodes (-ndm 3 -ndf 3)
$cNode constrained node (-ndm 3 -ndf 3)
$lNode Lagrange multiplier node (-ndm 3 -ndf 3)
$matTag unique integer tag associated with previously-defined nDMaterial object
$gTol gap tolerance
$fTol force tolerance




The SimpleContact3D element is a three-dimensional node-to-surface contact element which defines a frictional contact interface between two separate bodies. The retained nodes are the nodes which define a surface of a hexahedral element on the first body, and the constrained node is a node from the second body. The Lagrange multiplier node is required to enforce the contact condition. This node should not be shared with any other element in the domain. Information on the theory behind this element can be found in, e.g. Wriggers (2002).

NOTE:

  1. The SimpleContact3D element has been written to work exclusively with the ContactMaterial3D nDMaterial object.
  2. The valid recorder queries for this element are:
    1. force - returns the contact force acting on the constrained node in vector form.
    2. frictionforce - returns the frictional force acting on the constrained node in vector form.
    3. forcescalar - returns the scalar magnitudes of the single normal and two tangential contact forces.
    4. The SimpleContact3D elements are set to consider frictional behavior as a default, but the frictional state of the SimpleContact3D element can be changed from the input file using the setParameter command. When updating, value of 0 corresponds to the frictionless condition, and a value of 1 signifies the inclusion of friction. An example command for this update procedure is provided below
  3. The SimpleContact3D element works well in static and pseudo-static analysis situations.
  4. In transient analysis, the presence of the contact constraints can effect the stability of commonly-used time integration methods in the HHT or Newmark family (e.g., Laursen, 2002). For this reason, use of alternative time-integration methods which numerically damp spurious high frequency behavior may be required. The TRBDF2 integrator is an effective method for this purpose. The Newmark integrator can also be effective with proper selection of the gamma and beta coefficients. The trapezoidal rule, i.e., Newmark with gamma = 0.5 and beta = 0.25, is particularly prone to instability related to the contact constraints and is not recommended.


EXAMPLE:

SimpleContact3D element with tag 1, connectivity with nodes 1, 2, 3, 4, 5, 6, 7, and 8, material with tag 1, and gap and force tolerances of 1.0e-10.

element SimpleContact3D 1  1 2 3 4 5 6 7 8  1  1.0e-10 1.0e-10

Update all of the SimpleContact3D elements with tags between 10 and 20 to consider a frictionless interface

setParameter -value 0 -eleRange 10 20 friction


REFERENCES:

  1. Wriggers, P. (2002). Computational Contact Mechanics. John Wiley & Sons, Ltd, West Sussex, England.
  2. Laursen, T. A. (2002). Computational Contact and Impact Mechanics. Springer-Verlag, Berlin.

Code Developed by: Kathryn Petek, Pedro Arduino, & Peter Mackenzie-Helnwein, at the University of Washington