Joint II

This example illustrates the usage of *CONSTRAINED_JOINT keyword to model a translational joint.

*CONSTRAINED_EXTRA_NODES_NODE
*CONSTRAINED_EXTRA_NODES_SET
*CONSTRAINED_JOINT_TRANSLATIONAL
*CONTROL_CPU
*CONTROL_ENERGY
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_GLSTAT
*DATABASE_JNTFORC
*DATABASE_MATSUM
*ELEMENT_SHELL
*END
*INITIAL_VELOCITY_NODE
*KEYWORD
*MAT_ELASTIC
*MAT_RIGID
*NODE
*PART
*SECTION_SHELL
*SET_NODE
*TITLE

  *KEYWORD
  *TITLE
  translational joint between shells 
  $
  $ - one shell is elastic (to control time step size)
  $
  $ - no stiffness in joint, so local coordinate system is
  $   not needed
  *CONSTRAINED_JOINT_translational
  $         2        10         1         9        11        12 0.100E+01
  $ 3rd node pair need not form vector perpendicular to free direction
  $         2        10         1         9         4       104 0.100E+01
           2        10         1         9         3       103 0.100E+01
  *constrained_extra_nodes_node
  1,11
  2,12
  *CONTROL_TERMINATION
   2.000E-02         0 0.000E+00         0 0.000E+00
  *DATABASE_BINARY_D3PLOT
   5.000E-04
  *initial_velocity_node
  $23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 
           3       90.            300.
           4       90.            300.
  *control_energy
           2         0         2
  *database_glstat
   1.000E-04
  *database_matsum
   1.000E-04
  *database_jntforc
   1.000E-04
  *CONTROL_CPU
   0.000E+00
  *MAT_RIGID
           1 7.000E-04 3.000E+07 3.000E-01 0.000E+00 0.000E+00 0.000E+00          
   0.000E+00 0.000E+00 0.000E+00
   0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
  *MAT_RIGID
           2 7.000E-04 3.000E+07 3.000E-01 0.000E+00 0.000E+00 0.000E+00          
   0.000E+00 0.000E+00 0.000E+00
   0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
  *MAT_ELASTIC
           3 7.000E-04 3.000E+07 3.000E-01
  *SECTION_SHELL
           1         0 0.000E+00 0.000E+00 0.000E+00 0.000E+00
   1.000E-01 1.000E-01 1.000E-01 1.000E-01 0.000E+00
  *SECTION_SHELL
           2         0 0.000E+00 0.000E+00 0.000E+00 0.000E+00
   1.000E-01 1.000E-01 1.000E-01 1.000E-01 0.000E+00
  *SECTION_SHELL
           3         0 0.000E+00 0.000E+00 0.000E+00 0.000E+00
   1.000E-01 1.000E-01 1.000E-01 1.000E-01 0.000E+00
  *PART
   
           1         1         1         0
  *PART
   
           2         2         2         0
  *PART
   
           3         3         3         0
  *NODE
  $
  *ELEMENT_SHELL
  $
  *constrained_extra_nodes_set 
  2,104
  *set_node
  104
  104,103
  *END