ZeroLengthImpact3D: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 64: Line 64:
----
----


SAMPLE COMMAND:
SAMPLE COMMAND OF FOUR NODES ON THE SAME CONTACT SURFACE:


<source lang="Tcl">
<source lang="Tcl">

Revision as of 20:40, 29 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 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 rectangular master and slave contact surfaces, four zeroLengthImpact3D elements can be used at the corners.




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.0e5; # penalty stiffness for normal direction 
set Kn2 [expr $Kn * 0.1]; # penalty stiffness after yielding, based on Hertz impact model 
set Delta_y 0.1; # 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;
element zeroLengthImpact3D 2212 22 12 $direction $initGap $frictionRatio $Kt $Kn $Kn2 $Delta_y $cohesion;
element zeroLengthImpact3D 2414 24 14 $direction $initGap $frictionRatio $Kt $Kn $Kn2 $Delta_y $cohesion;
element zeroLengthImpact3D 2515 25 15 $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)