ZeroLengthImpact3D

• Command_Manual
  • Tcl Commands
  • Modeling_Commands
    • model
    • uniaxialMaterial
    • ndMaterial
    • frictionModel
    • section
    • geometricTransf
    • element
    • node
    • sp commands
      • fix
      • fixX
      • fixY
      • fixZ
    • mp commands
      • equalDOF
      • rigidDiaphragm
      • rigidLink
    • timeSeries
    • pattern
    • mass
    • block commands
      • block2D
      • block3D
    • region
    • rayleigh
  • Analysis Commands
    • constraints Command
    • numberer Command
    • system Command
    • test Command
    • algorithm Command
    • integrator Command
    • analysis Command
    • eigen Command
    • analyze Command
  • Output Commands
    • recorder
    • record
    • print
    • printA
    • logFile
    • RealTime_Output_Commands
  • Misc Commands
    • eleResponse Command
    • getTime Command
    • getEleTags Command
    • getNodeTags
    • loadConst Command
    • nodeCoord Command
    • nodeDisp Command
    • nodeVel Command
    • nodeAccel Command
    • nodeEigenvector Command
    • nodeBounds
    • print Command
    • remove Command
    • reset Command
    • testIter
    • testNorms
    • wipe Command
    • exit Command
  • DataBase Commands
    • Template:DataBaseCommandsMenu

---

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)