ZeroLengthInterface2D: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(8 intermediate revisions by 3 users not shown)
Line 4: Line 4:


{|  
{|  
| style="background:yellow; color:black; width:850px" | '''element zeroLengthInterface2D $eleTag -sNdNum $sNdNum -mNdNum $mNdNum –dof $sdof $mdof -Nodes $Nodes $Kn $Kt $phi'''
| style="background:yellow; color:black; width:840px" | '''element zeroLengthInterface2D $eleTag -cNdNum $cNdNum -rNdNum $rNdNum –dof $cdof $rdof -Nodes $Nodes $Kn $Kt $phi'''
|}
|}


Line 14: Line 14:
|  style="width:150px" | '''$eleTag''' || unique element object tag
|  style="width:150px" | '''$eleTag''' || unique element object tag
|-
|-
| '''$sNdNum''' || Number of Slave Nodes
| '''$cNdNum''' || Number of Constrained Nodes
|-
|-
| '''$mNdNum''' || Number of Master nodes   
| '''$rNdNum''' || Number of Retained nodes   
|-
|-
| '''$$sdof $mdof''' || Slave and Master degree of freedom   
| '''$cdof $rdof''' || Constrained and Retained degree of freedom   
|-
|-
| '''$Nodes ...''' || Slave and master node tags respectively
| '''$Nodes ...''' || Constrained and Retained node tags respectively
|-
|-
| '''$Kn''' || Penalty in normal direction
| '''$Kn''' || Penalty in normal direction
Line 32: Line 32:


# The contact element is node-to-segment (NTS) contact. The relation follows Mohr-Coulomb frictional law: <math>T = N \times tan(\phi)</math>, where <math>T</math> is the tangential force, <math>N</math> is normal force across the interface and <math>\phi\,\!</math> is friction angle.
# The contact element is node-to-segment (NTS) contact. The relation follows Mohr-Coulomb frictional law: <math>T = N \times tan(\phi)</math>, where <math>T</math> is the tangential force, <math>N</math> is normal force across the interface and <math>\phi\,\!</math> is friction angle.
# For 2D contact, slave nodes and master nodes must be 2 DOF and notice that the slave and master nodes must be entered in counterclockwise order.
# For 2D contact, constrained nodes and retained nodes must be 2 DOF and notice that the constrained and retained nodes must be entered in counterclockwise order.
# The resulting tangent from the contact element is non-symmetric. Switch to the non-symmetric matrix solver if convergence problem is experienced.  
# The resulting tangent from the contact element is non-symmetric. Switch to the non-symmetric matrix solver if convergence problem is experienced.  
# As opposed to node-to-node contact, predefined normal vector for node-to-segment (NTS) element is not required because contact normal will be calculated automatically at each step. And also this element can handle contact between different DOFs such as beam-beam, beam-solid and solid-solid.
# As opposed to node-to-node contact, predefined normal vector for node-to-segment (NTS) element is not required because contact normal will be calculated automatically at each step. And also this element can handle contact between different DOFs such as beam-beam, beam-solid and solid-solid.
Line 39: Line 39:
'''EXAMPLE:'''
'''EXAMPLE:'''


[[File:zeroLengthContactNTS2D _fig1.jpg]]  
[[File:ZeroLengthInterface2Drc.png]]  


element zeroLengthContactNTS2D 1 -sNdNum 6 -mNdNum 6 - Nodes 5 10 12 3 9 11 1 4 2 8 7 1e8 1e8 16
element zeroLengthInterface2D 1 -sNdNum 6 -mNdNum 6 –dof 2 3 -Nodes 5 10 12 3 9 11 1 4 2 8 7 6 1e8 1e8 16


'''Example 1:'''
'''Example 1:'''


This example simply shows the two quadrilateral elements in normal contact. The top element is in normal downward uniform force. The Tcl script of this example can be found [[ZeroLengthContactNTS2D_Example1|here]].  
This example simply shows the quadrilateral elements in normal contact on top of beam element. The top element is in normal downward uniform force. The Tcl script of this example can be found [[ZeroLengthInterface2D_Example1|here]].  


[[File:ZeroLengthContactNTS2D_fig2.jpg]]   
[[File:ZeroLengthInterface2D_fig2.jpg]]   


'''Example 2:'''
'''Example 2:'''


This example shows two cantilever beams in contact. The beams were modeled using four-node quadrilateral elements and the end of top beam was subjected to a linearly increasing displacement.  
This example shows two cantilever beams in contact. The beams were modeled using four-node quadrilateral elements and the end of top beam was subjected to a linearly increasing displacement.  
The Tcl scripts for this example can be found [[ZeroLengthContactNTS2D_Example2|here]].  
The Tcl scripts for this example can be found [[ZeroLengthInterface2D_Example2|here]].  


[[File:zeroLengthContactNTS2D_fig3.jpg]]  
[[File:zeroLengthContactNTS2D_fig3.jpg]]  
Line 67: Line 67:
----
----


Code Developed by: <span style="color:blue"> [http://www2.decf.berkeley.edu/~roozbehg/index.html Roozbeh G. Mikola], UC Berkeley</span> and <span style="color:blue"> [http://www.ce.berkeley.edu/~sitar/ N. Sitar], UC Berkeley</span>
Code Developed by: <span style="color:blue"> [http://www.roozbehgm.com/ Roozbeh G. Mikola], UC Berkeley</span> and <span style="color:blue"> [http://www.ce.berkeley.edu/~sitar/ N. Sitar], UC Berkeley</span>

Latest revision as of 15:58, 14 June 2020




zeroLengthInterface2D is Node-to-Segment (NTS) frictional contact element used in two dimensional analysis for contact between elements. Please notice that any number of DOF (DOF>1) can be handled using this element (i.e. beam-solid, solid-solid as well as beam-beam contact) as oppose to zeroLengthContactNTS2D element).

element zeroLengthInterface2D $eleTag -cNdNum $cNdNum -rNdNum $rNdNum –dof $cdof $rdof -Nodes $Nodes $Kn $Kt $phi




$eleTag unique element object tag
$cNdNum Number of Constrained Nodes
$rNdNum Number of Retained nodes
$cdof $rdof Constrained and Retained degree of freedom
$Nodes ... Constrained and Retained node tags respectively
$Kn Penalty in normal direction
$Kt Penalty in tangential direction
$phi Friction angle in degrees

NOTES:

  1. The contact element is node-to-segment (NTS) contact. The relation follows Mohr-Coulomb frictional law: <math>T = N \times tan(\phi)</math>, where <math>T</math> is the tangential force, <math>N</math> is normal force across the interface and <math>\phi\,\!</math> is friction angle.
  2. For 2D contact, constrained nodes and retained nodes must be 2 DOF and notice that the constrained and retained nodes must be entered in counterclockwise order.
  3. The resulting tangent from the contact element is non-symmetric. Switch to the non-symmetric matrix solver if convergence problem is experienced.
  4. As opposed to node-to-node contact, predefined normal vector for node-to-segment (NTS) element is not required because contact normal will be calculated automatically at each step. And also this element can handle contact between different DOFs such as beam-beam, beam-solid and solid-solid.
  5. contact element is implemented to handle large deformations.

EXAMPLE:

element zeroLengthInterface2D 1 -sNdNum 6 -mNdNum 6 –dof 2 3 -Nodes 5 10 12 3 9 11 1 4 2 8 7 6 1e8 1e8 16

Example 1:

This example simply shows the quadrilateral elements in normal contact on top of beam element. The top element is in normal downward uniform force. The Tcl script of this example can be found here.

Example 2:

This example shows two cantilever beams in contact. The beams were modeled using four-node quadrilateral elements and the end of top beam was subjected to a linearly increasing displacement. The Tcl scripts for this example can be found here.

The following Figure shows the deflections of the two beams.

REFERENCES:

  1. P. Wriggers, V.T. Vu and E. Stein, Finite-element formulation of large deformation impact–contact problems with friction, Comput. Struct. 37 (1990), pp. 319–331.
  2. Peter Wriggers. Computational Contact Mechanics. John Wiley & Sons Ltd. Chichester, 2002.

Code Developed by: Roozbeh G. Mikola, UC Berkeley and N. Sitar, UC Berkeley