ZeroLengthContact Element: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
No edit summary
mNo edit summary
 
(2 intermediate revisions by the same user not shown)
Line 6: Line 6:


{|  
{|  
| style="background:yellow; color:black; width:800px" | '''element zeroLengthContact2D $eleTag $sNode $mNode $Kn $Kt $mu -normal $Nx $Ny'''
| style="background:yellow; color:black; width:800px" | '''element zeroLengthContact2D $eleTag $cNode $rNode $Kn $Kt $mu -normal $Nx $Ny'''
|}
|}


Line 12: Line 12:


{|  
{|  
| style="background:yellow; color:black; width:800px" | '''element zeroLengthContact3D $eleTag $sNode $mNode $Kn $Kt $mu $c $dir'''
| style="background:yellow; color:black; width:800px" | '''element zeroLengthContact3D $eleTag $cNode $rNode $Kn $Kt $mu $c $dir'''
|}
|}


Line 21: Line 21:
|  style="width:150px" | '''$eleTag''' || unique element object tag
|  style="width:150px" | '''$eleTag''' || unique element object tag
|-
|-
| '''$sNode''' || Slave Node tag
| '''$cNode''' || Constrained node tag
|-
|-
| '''$mNode''' || Master node tag  
| '''$rNode''' || Retained node tag  
|-
| '''$Nodes ...''' || Slave and master node tags respectively
|-
|-
| '''$Kn''' || Penalty in normal direction
| '''$Kn''' || Penalty in normal direction
Line 37: Line 35:
| '''$dir''' || Direction flag of the contact plane (3D), it can be:
| '''$dir''' || Direction flag of the contact plane (3D), it can be:


1 Out normal of the master plane pointing to +X direction
1 Out normal of the retained plane pointing to +X direction


2 Out normal of the master plane pointing to +Y direction
2 Out normal of the retained plane pointing to +Y direction


3 Out normal of the master plane pointing to +Z direction
3 Out normal of the retained plane pointing to +Z direction
|}
|}


[[Image:ZeroLengthContact.png]]
[[Image:ZeroLengthContact2.png]]


NOTES:
NOTES:


# The contact element is node-to-node contact. Contact occurs between two contact nodes when they come close. The relation follows Mohr-coulomb law: T = $mu * N + $c, where T is tangential force and N is normal force across the interface. $mu is friction coefficient and $c is total cohesion (summed over the effective area of contact nodes)
# The contact element is node-to-node contact. Contact occurs between two contact nodes when they come close. The relation follows Mohr-coulomb law: T = $mu * N + $c, where T is tangential force and N is normal force across the interface. $mu is friction coefficient and $c is total cohesion (summed over the effective area of contact nodes)
# The contact node pair in node-to-node contact element is termed "master node" and "slave node", respectively. Master/slave plane is the contact plane which the master/slave node belongs to. The discrimination is made solely for contact detection purpose. User need to specify the corresponding out normal of the master plane, and this direction is assumed to be unchanged during analysis. For simplicity, 3D contact only allows 3 options to specify the directions of the contact plane. The convention is: out normal of master plane always points to positive axial direction (+X or +Y, or +Z)
# The contact node pair in node-to-node contact element is termed "retained node" and "constrained node", respectively. Retained/constrained plane is the contact plane which the retrained/constrained node belongs to. The discrimination is made solely for contact detection purpose. User need to specify the corresponding out normal of the master plane, and this direction is assumed to be unchanged during analysis. For simplicity, 3D contact only allows 3 options to specify the directions of the contact plane. The convention is: out normal of master plane always points to positive axial direction (+X or +Y, or +Z)
# For 2D contact, slave nodes and master nodes must be 2 DOF. For 3D contact, slave nodes and master nodes must be 3 DOF.  
# For 2D contact, constrained nodes and retained nodes must be 2 DOF. For 3D contact, constrained nodes and retained nodes must be 3 DOF.  
# The resulted tangent from the contact element is '''NON-SYMMETRIC'''. Switch to non-symmetric matrix solver.  
# The resulted tangent from the contact element is '''NON-SYMMETRIC'''. Switch to non-symmetric matrix solver.  



Latest revision as of 17:37, 13 June 2020




This command is used to construct a zeroLengthContact2D element or a zeroLengthContact3D element, which are Node-to-node frictional contact element used in two dimensional analysis and three dimensional analysis:

2d analysis:

element zeroLengthContact2D $eleTag $cNode $rNode $Kn $Kt $mu -normal $Nx $Ny

3d analysis:

element zeroLengthContact3D $eleTag $cNode $rNode $Kn $Kt $mu $c $dir


$eleTag unique element object tag
$cNode Constrained node tag
$rNode Retained node tag
$Kn Penalty in normal direction
$Kt Penalty in tangential direction
$mu friction coefficient
$c cohesion (not available in 2D)
$dir Direction flag of the contact plane (3D), it can be:

1 Out normal of the retained plane pointing to +X direction

2 Out normal of the retained plane pointing to +Y direction

3 Out normal of the retained plane pointing to +Z direction

NOTES:

  1. The contact element is node-to-node contact. Contact occurs between two contact nodes when they come close. The relation follows Mohr-coulomb law: T = $mu * N + $c, where T is tangential force and N is normal force across the interface. $mu is friction coefficient and $c is total cohesion (summed over the effective area of contact nodes)
  2. The contact node pair in node-to-node contact element is termed "retained node" and "constrained node", respectively. Retained/constrained plane is the contact plane which the retrained/constrained node belongs to. The discrimination is made solely for contact detection purpose. User need to specify the corresponding out normal of the master plane, and this direction is assumed to be unchanged during analysis. For simplicity, 3D contact only allows 3 options to specify the directions of the contact plane. The convention is: out normal of master plane always points to positive axial direction (+X or +Y, or +Z)
  3. For 2D contact, constrained nodes and retained nodes must be 2 DOF. For 3D contact, constrained nodes and retained nodes must be 3 DOF.
  4. The resulted tangent from the contact element is NON-SYMMETRIC. Switch to non-symmetric matrix solver.



EXAMPLES:

Gang Wang to provide a smart example!



REFERENCES:



Code Developed by: Gang Wang, Geomatrix