ZeroLengthImpact3D: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
No edit summary
(Change terminology to retained/constrained)
 
(172 intermediate revisions by one other user not shown)
Line 1: Line 1:
{{CommandManualMenu}}
{{CommandManualMenu}}


This command constructs a node-to-node zero-length contact element in 3D space.  
This command constructs a node-to-node zero-length contact element in 3D space to simulate the impact/pounding and friction phenomena.  


{|  
{|  
| style="background:yellow; color:black; width:800px" | '''element zeroLengthImpact3D $tag $slaveNode $masterNode $direction $initGap $frictionRatio $Kt $Kn $Kn2 $Delta_y $cohesion'''
| style="background:yellow; color:black; width:800px" | '''element zeroLengthImpact3D $tag $cNode $rNode $direction $initGap $frictionRatio $Kt $Kn $Kn2 $Delta_y $cohesion'''
|}
|}


----
----


{|
{|
|  style="width:150px" | '''$tag''' || unique element object tag
|  style="width:150px" | '''$tag''' || Unique element object tag
|-
|-
| '''$slaveNode''' || Slave node tag
| '''$cNode''' || Constrained node tag
|-
|-
| '''$masterNode''' || Master node tag
| '''$rNode''' || Retained node tag
|-
|-
| '''$direction''' ||  
| '''$direction''' ||  


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


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


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


|-
|-
| '''$initGap''' || Initial gap between master plane and slave plane
| '''$initGap''' || Initial gap between retained plane and constrained plane
|-
|-
| '''$frictionRatio''' || Friction ratio  
| '''$frictionRatio''' || Friction ratio in two tangential directions (parallel to retained and constrained planes)
|-
|-
| '''$Kt''' || Penalty in tangential directions (parallel to master and slave planes)
| '''$Kt''' || Penalty in two tangential directions
|-
|-
| '''$Kn''' || Penalty in normal direction (normal to master and slave planes)
| '''$Kn''' || Penalty in normal direction (normal to retained and constrained planes)
|-
|-
| '''$Kn2''' || Penalty in normal direction after yielding based on Hertz impact model  
| '''$Kn2''' || Penalty in normal direction after yielding based on Hertz impact model  
Line 42: Line 40:
|}
|}


----


NOTES:
# This element has been developed on top of the “zeroLengthContact3D”. All the notes available in [[ ZeroLengthContact_Element | “zeroLengthContact3D” wiki page ]] would apply to this element as well. It includes the definition of retained and constrained nodes, the number of degrees of freedom in the domain, etc.
# Regarding the number of degrees of freedom (DOF), the end nodes of this element should be defined in 3DOF domain. For getting information on how to use 3DOF and 6DOF domain together, please refer to OpenSees documentation and forums or see the zip file provided in the EXAMPLES section below.
# This element adds the capabilities of [[ Impact_Material | “ImpactMaterial” ]] to [[ ZeroLengthContact_Element | “zeroLengthContact3D.” ]]
# For simulating a surface-to-surface contact, the element can be defined for connecting the nodes on constrained surface to the nodes on retained surface.
# The element was found to be fast-converging and eliminating the need for extra elements and nodes in the modeling process.


NOTES:
----


# 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)
EXAMPLES:  
# 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)
# For 2D contact, slave nodes and master nodes must be 2 DOF. For 3D contact, slave nodes and master nodes must be 3 DOF.
# The resulted tangent from the contact element is '''NON-SYMMETRIC'''. Switch to non-symmetric matrix solver.


#The following zip file contains an example script and the corresponding input cyclic displacement: [[File:Example script 2.zip]]
#The following zip file contains an example script on how to use 6DOF domain and 3DOF domain together: [[File:Example script 6DOF 3DOF.zip]]


<!--
----
----


EXAMPLES:
SAMPLE COMMAND (example scripts are available at bottom of this page):


Gang Wang to provide a smart example!
<source lang="Tcl">


</source>
-->
----
----


REFERENCES:


REFERENCES:
[[ ZeroLengthContact_Element | zeroLengthContact3D ]] , [[ Impact_Material | ImpactMaterial ]]


----


CODE DEVELOPED BY:
:<span style="color:blue"> Dr. Arash E. Zaghi and Majid Cashany at University of Connecticut (UConn) </span>


----
----


Code Developed by: <span style="color:blue"> Arash E. Zaghi, Majid Cashany, University of Connecticut (UConn) </span>
APPLICATIONS:
# This element has been employed to simulate the bridge hinges including superstructure-abutment interaction at the University of Connecticut (UConn) and University of Nevada, Reno (UNR).
# It has been implemented in non-structural systems like suspended ceilings, simulating the impact/pounding and friction phenomena.
 
 
<!--
After running the example script, the following hysteresis loop is resulted in normal direction under cyclic excitation:
 
[[File: HysteresisLoop.png]]
-->

Latest revision as of 16:19, 13 June 2020




This command constructs a node-to-node zero-length contact element in 3D space to simulate the impact/pounding and friction phenomena.

element zeroLengthImpact3D $tag $cNode $rNode $direction $initGap $frictionRatio $Kt $Kn $Kn2 $Delta_y $cohesion

$tag Unique element object tag
$cNode Constrained node tag
$rNode Retained node tag
$direction

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

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

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

$initGap Initial gap between retained plane and constrained plane
$frictionRatio Friction ratio in two tangential directions (parallel to retained and constrained planes)
$Kt Penalty in two tangential directions
$Kn Penalty in normal direction (normal to retained and constrained 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 retained and constrained nodes, the number of degrees of freedom in the domain, etc.
  2. Regarding the number of degrees of freedom (DOF), the end nodes of this element should be defined in 3DOF domain. For getting information on how to use 3DOF and 6DOF domain together, please refer to OpenSees documentation and forums or see the zip file provided in the EXAMPLES section below.
  3. This element adds the capabilities of “ImpactMaterial” to “zeroLengthContact3D.”
  4. For simulating a surface-to-surface contact, the element can be defined for connecting the nodes on constrained surface to the nodes on retained surface.
  5. The element was found to be fast-converging and eliminating the need for extra elements and nodes in the modeling process.

EXAMPLES:

  1. The following zip file contains an example script and the corresponding input cyclic displacement: File:Example script 2.zip
  2. The following zip file contains an example script on how to use 6DOF domain and 3DOF domain together: File:Example script 6DOF 3DOF.zip

REFERENCES:

zeroLengthContact3D , ImpactMaterial


CODE DEVELOPED BY:

Dr. Arash E. Zaghi and Majid Cashany at University of Connecticut (UConn)

APPLICATIONS:

  1. This element has been employed to simulate the bridge hinges including superstructure-abutment interaction at the University of Connecticut (UConn) and University of Nevada, Reno (UNR).
  2. It has been implemented in non-structural systems like suspended ceilings, simulating the impact/pounding and friction phenomena.