ZeroLengthImpact3D

From OpenSeesWiki
Revision as of 18:21, 28 May 2013 by Cashany (talk | contribs)
Jump to navigation Jump to search




This command constructs a node-to-node zero-length contact element in 3D space.

element zeroLengthImpact3D $tag $slaveNode $masterNode $direction $initGap $frictionRatio $Kt $Kn $Kn2 $Delta_y $cohesion




$tag unique element object tag
$slaveNode Slave node tag
$masterNode Master node tag
$direction

1 if normal vector of master plane points to +X direction

2 if normal vector of master plane points to +Y direction

3 if normal vector of master plane points to +Z direction

$initGap Initial gap between master plane and slave plane
$frictionRatio Friction ratio
$Kt Penalty in tangential directions (parallel to contact plane)
$Kn Penalty in normal direction (normal to contact plane)
$Kn2 Penalty in normal direction after yielding based on Hertz impact model
? ?
? ?


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 "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)
  3. For 2D contact, slave nodes and master nodes must be 2 DOF. For 3D contact, slave nodes and master 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: Arash E. Zaghi, Majid Cashany, University of Connecticut (UConn)