ZeroLengthImpact3D: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 27: Line 27:


|-
|-
| '''$Kn''' || Penalty in normal direction
| '''$initGap''' || Initial gap between master plane and slave plane
|-
|-
| '''$Kt''' || Penalty in tangential direction
| '''$frictionRatio''' || Friction ratio
|-
|-
| '''$mu''' || friction coefficient
| '''$Kt''' || Penalty in tangential directions (parallel to contact plane)
|-
|-
| '''$c''' || cohesion (not available in 2D)
| '''$Kn''' || Penalty in normal direction (normal to contact plane)
|-
|-
| '''$dir''' || Direction flag of the contact plane (3D), it can be:
| '''$Kn2''' || Penalty in normal direction after yielding based on Hertz impact model
 
|-
1 Out normal of the master plane pointing to +X direction
| '''?''' || ?
 
|-
2 Out normal of the master plane pointing to +Y direction
| '''?''' || ?
 
3 Out normal of the master plane pointing to +Z direction
|}
|}



Revision as of 18:21, 28 May 2013




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)