Documentation for zeroLengthInterface2DUpdate

source/user/manual/model/elements/zeroLengthInterface2DUpdate.rst

@@ -0,0 +1,167 @@
+zeroLengthInterface2DUpdate
+===========================
+
+**zeroLengthInterface2DUpdate** is an updated version of **zeroLengthInterface2D**. 
+
+**zeroLengthInterface2DUpdate** has the following improvements:
+
+1. Fixed the incorrect friction behavior in **zeroLengthInterface2D** when the interface is under cyclic sliding motion.
+2. Fixed the incorrect normal pressure behavior in **zeroLengthInterface2D** due to the repetitive search algorithm.
+3. Allows the user to change the friction coefficient during the analysis.
+
+\* **zeroLengthInterface2D** was written by Roozbeh Geraili Mikola (roozbehg@berkeley.edu) and Prof. Nicholas Sitar (nsitar@ce.berkeley.edu) on July 02 2010.
+
+User Instruction
+----------------
+
+1. **Element Command Syntax**
+
+   The tcl input command for the **zeroLengthInterface2DUpdate** element is given below:
+
+   .. code-block:: tcl
+
+      element zeroLengthInterface2DUpdate eleTag? -sNdNum sNdNum? -pNdNum pNdNum? –dof sdof? pdof? -Nodes Nodes? Kn? Kt? phi?
+
+   $eleTag:	unique element object tag
+
+   $sNdNum:	Number of nodes on secondary surface
+
+   $pNdNum:	Number of nodes on primary surface
+
+   $sdof $pdof:	Nodal degree of freedom for secondary and primary surface
+
+   $Nodes:	Secondary and primary node tags respectively
+
+   $Kn:	Penalty in normal direction
+
+   $Kt:	Penalty in tangential direction
+
+   $phi:	Friction angle in degrees
+
+2. **Tcl Input Example**
+
+   .. image:: https://github.com/user-attachments/assets/ec851ee3-1093-4894-bbaf-b72d0e6b1184
+      :alt: Tcl input example
+
+   To define a zeroLengthInterface2DUpdate element between a quad element (in blue) and beam elements (in red).
+
+   .. code-block:: tcl
+
+      element zeroLengthInterface2DUpdate 1 -sNdNum 6 -pNdNum 2 -dof 3 2 -Nodes 10 9 8 7 6 5 4 1 $Kn $Kt $phi
+
+
+3. **Changing the Friction Coefficient**
+
+   To update the friction coefficient during analysis, modify the parameter named "phi". For example:
+
+   .. code-block:: tcl
+
+      set Kn 1e6
+      set Kt 1e6
+      set phi_ini 0.0
+      element zeroLengthInterface2DUpdate 1 -sNdNum 6 -pNdNum 2 -dof 3 2 -Nodes 10 9 8 7 6 5 4 1 $Kn $Kt $phi_ini
+      parameter 1 element 1 phi
+      updateParameter 1 [expr 16.7]
+
+   *(The above changes the friction coefficient in element 1 from 0 to 0.3.)*
+
+Full TCL Example Code
+----------------
+
+.. code-block:: tcl
+
+    # ----------------------------
+    # Start of model generation
+    # ----------------------------
+    wipe 
+    model basic -ndm 2 -ndf 2
+
+    nDMaterial ElasticIsotropic   1   1e10   0.49  6.75 
+    uniaxialMaterial Elastic 2 [expr 1e10];
+
+    # ################################ 
+    # Build the model 
+    # ################################# 
+    node 1 6 0
+    node 2 6 1 
+    node 3 4 1
+    node 4 4 0
+    fix 3 1 0 0 
+    element quad 1 1 2 3 4 1 PlaneStrain 1 
+
+    set Kn 1e6; 
+    set Kt 1e6; 
+    set phi 16.7; 
+
+    pattern Plain 1 Linear {
+      sp 2 2 -1.e-2
+      sp 3 2 -1.e-2 
+    }
+
+    model BasicBuilder -ndm 2 -ndf 3;
+
+    node 5 0 0 
+    node 6 2 0 
+    node 7 4 0
+    node 8 6 0
+    node 9 8 0
+    node 10 10 0
+
+    geomTransf Linear 1
+    section Fiber 1 {
+       patch rect 2 100 1 [expr -15.0] -0.5 [expr 15.0] 0.5
+    }
+
+    element dispBeamColumn 2 5 6 5 1 1
+    element dispBeamColumn 3 6 7 5 1 1
+    element dispBeamColumn 4 7 8 5 1 1
+    element dispBeamColumn 5 8 9 5 1 1
+    element dispBeamColumn 6 9 10 5 1 1
+
+    element zeroLengthInterface2DUpdate 7 -sNdNum 6 -pNdNum 2 -dof 3 2 -Nodes 10 9 8 7 6 5 4 1 [expr $Kn] [expr $Kt] [expr $phi]
+    fix 5 1 1 0
+    fix 10 0 1 0   
+
+    # ----------------------------
+    # Start of recorder generation
+    # ----------------------------
+    recorder Node  -file  reactionForce.out  -node  5 10 -time -dof 1 2 reaction;
+    recorder Node -file Node2.out -node 2 -dof 1 2 -time disp 
+    recorder Node -file Node3.out -node 3 -dof 1 2 -time disp
+    recorder Element -file Contactele.out -ele 7 -time force
+
+    # --------------------------------------------------------------------
+    # Start of static analysis (creation of the analysis & analysis itself)
+    # --------------------------------------------------------------------
+    system SparseGeneral
+    numberer RCM
+    constraints Penalty   1.e+018   1.e+018
+    test NormDispIncr  1.00e-0010 500 2
+    algorithm KrylovNewton
+    integrator  LoadControl 0.1
+    analysis    Static
+
+    analyze     10
+    puts "\ninitiate contact finished"
+
+    wipeAnalysis
+
+    model BasicBuilder -ndm 2 -ndf 2;
+    loadConst -time 0.0
+    remove sp 3 1
+    remove sp 3 2
+
+    pattern Plain 2 Linear {
+      sp 3 1 1.e-1 
+    }
+
+    constraints Transformation
+    test        NormDispIncr 1e-7 500 2
+    algorithm   KrylovNewton
+    numberer    RCM
+    system      ProfileSPD
+    integrator  LoadControl 0.01
+    analysis    Static
+
+    analyze     100
+    puts "\nsliding finished\n"