*KEYWORD
$ 3x3 lobotto = 3x3 gauss for tubular beams
$  R of IP is sqrt [ (ro^2 + ri^2)/2 ]  (first moment concept)
*damping_global
,50
$ circular 3x3 lobatto vs. gauss
$
$ beam IP output in numerical order (see p. 3.14 in 940 structured user's manual)
$  no springback info written for beams
*interface_springback_dyna3d
5
$*interface_springback_nike3d
$5
*set_part
5
1
*TITLE
3hlbeams.k
*CONTROL_TERMINATION
 0.9000000         0 0.0000000         0 0.0000000
*CONTROL_TIMESTEP
 0.0000000 0.0000000         0 0.0000000 0.0000000         0         0         0
*CONTROL_SHELL
 0.0000000         0         0         0         0         0         0
*CONTROL_DAMPING
         0 0.0000000 0.0000000 0.0000000 0.0000000         0 0.0000000         0
*CONTROL_CONTACT
 0.0000000 0.0000000         0         0         0         0         0
         0         0         0         0 0.0000000         0         0         0
*CONTROL_OUTPUT
         0         0         0         0 0.0000000         0         0
*CONTROL_ENERGY
         1         2         1         1
*DATABASE_BINARY_D3PLOT
 0.0100000                   0
*DATABASE_EXTENT_BINARY
         0         0         0         0         0         0         0         0
         0         0         9         0         0         0
*DATABASE_NODOUT
 2.00000-3
*DATABASE_ELOUT
 2.00000-3
*DATABASE_SPCFORC
 2.00000-3
*DATABASE_BNDOUT
 2.00000-3
*CONTROL_CPU
 0.0000000
*MAT_ELASTIC
         1 3.00000-4 1.00000+7 0.2000000
$*mat_plastic_kinematic
$$1,2e-4,1e7,.2,  20,0,1
$$1,2e-4,1e7,.2,  15,1e4,1
$1,2e-4,1e7,.2,  15000,1e4,1
$
*SECTION_BEAM
$$ 2x2 Gauss
$         1         1 0.0000000 2.0000000 1.0000000
$ 3x3 Lobatto circular
         1         1 0.0000000 4.0000000 1.0000000
2,2,1,1
$2,2,0,0
*SECTION_BEAM
$ 3x3 Gauss
         2         1 0.0000000 3.0000000 1.0000000
$ 1.0000000 1.0000000 1.0000000 1.0000000
2,2,1,1
*PART
material type # 1                                                               
         1         1         1         0         0         0         0         0
*PART
material type # 1                                                               
         2         2         1         0         0         0         0         0
*NODE
       1 1.800000000E+01 0.000000000E+00 0.000000000E+00       0       0
       2 3.600000000E+01 0.000000000E+00 0.000000000E+00       0       0
       3 5.400000000E+01 0.000000000E+00 0.000000000E+00       0       0
       4 7.200000000E+01 0.000000000E+00 0.000000000E+00       0       0
       5 0.000000000E+00 1.000000000E+00 0.000000000E+00       0       0
      11 1.800000000E+01 0.000000000E+00 9.000000000E+00       0       0
      12 3.600000000E+01 0.000000000E+00 9.000000000E+00       0       0
      13 5.400000000E+01 0.000000000E+00 9.000000000E+00       0       0
      14 7.200000000E+01 0.000000000E+00 9.000000000E+00       0       0
      15 0.000000000E+00 1.000000000E+00 9.000000000E+00       0       0
*database_history_beam
1,2,3,11,12,13
*ELEMENT_BEAM
       1       1       1       2       5
       2       1       2       3       5
       3       1       3       4       5
      11       2      11      12      15
      12       2      12      13      15
      13       2      13      14      15
*BOUNDARY_SPC_NODE
         1         0         1         1         1         1         1         1
        11         0         1         1         1         1         1         1
*DEFINE_CURVE
         1         0 0.0000000 0.0000000 0.0000000 0.0000000
      0.00000000E+00      0.00000000E+00
      1.00000001E-01      1.00000001E-01
      1.00000000E+00      1.00000001E-01
*LOAD_NODE_POINT
         4         2         1 2.0000000
        14         2         1 2.0000000
*END