ZeroLengthImpact3D

From OpenSeesWiki
Revision as of 20:50, 29 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 master and slave planes)
$Kn Penalty in normal direction (normal to master and slave planes)
$Kn2 Penalty in normal direction after yielding based on Hertz impact model
$Delta_y Yield deformation based on Hertz impact model
$cohesion Cohesion, if no cohesion, it is zero

NOTES:

  1. This element has been developed on top of the “zeroLengthContact3D”. All the notes available in “zeroLengthContact3D” wiki page would apply to this element as well. It includes the definition of master and slave nodes, the number of degrees of freedom in the domain, etc. Please refer to “zeroLengthContact3D” for detailed descriptions.
  2. This element adds the capabilities of “ImpactMaterial” to “zeroLengthContact3D”, without the need for extra nodes and elements in the modeling process. That is the capability of having an initial gap in the normal direction and also the capability of considering energy dissipation due to pounding/impact in the normal direction.
  3. This element has been employed to simulate the bridge hinge movements like superstructure-abutment interaction in bridges at University of Connecticut (UConn) and University of Nevada, Reno (UNR). This element was found to be fast-converging and eliminating the need for many extra elements and nodes in the modeling process.
  4. For simulating master and slave contact surfaces, this element can be defined for pairs of nodes on master and slave surface.




EXAMPLE SCRIPT:


The following zip file contains an example script and the corresponding input cyclic displacement:

File:ExampleScript.zip


SAMPLE COMMAND OF FOUR NODES ON THE SAME CONTACT SURFACE:

set direction 2; # direction of normal of contact surface
set initGap 0.5; # initial gap
set frictionRatio 0.1; # friction ratio 
set Kt 1.0e5; # penalty stiffness for tangential directions 
set Kn 1.0e3; # penalty stiffness for normal direction 
set Kn2 [expr $Kn * 0.1]; # penalty stiffness after yielding, based on Hertz impact model 
set Delta_y 0.001; # yield displacement based on Hertz impact model  
set cohesion 0.0; # cohesion

element zeroLengthImpact3D 2111 21 11 $direction $initGap $frictionRatio $Kt $Kn $Kn2 $Delta_y $cohesion;

After running the example script, the following hysteresis loop is resulted in normal direction under cyclic excitation:




REFERENCES:

zeroLengthContact3D

ImpactMaterial


Code Developed by: Dr. A.E. Zaghi, M. Cashany @ University of Connecticut (UConn)