ZeroLengthContact Element

From OpenSeesWiki
Jump to navigationJump to search




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

ZeroLengthContact2.png

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