*HEADING
 CONN_STD_FREE_3D: free vibration of spring-mass-dashpot  [CONN3D2]
 Damped free vibration:
  u(t) = e^{-w_n*t} u_0 ( xi/sqrt(1-xi*xi) sin(w_d*t) + cos(w_d*t) )
  u(0) = u_0 = 1.0, v(0) = 0.0
  w_n = sqrt(k/m)           = 2.0
  xi  = c/(2 m w_n)         = 5/48
  w_d = w_n sqrt(1 - xi*xi) = 1.98912
*************************************************
** create a rigid base for the connector elements
*************************************************
*NODE,NSET=PLATE
 101,  0.,  0., 1.
 112,  0., 11., 1.
 312,  2., 11., 1.
 314,  2., 13., 1.
 514,  4., 13., 1.
 518,  4., 17., 1.
 718,  6., 17., 1.
 701,  6.,  0., 1.
*NGEN,NSET=E1
 101,112,1
*NGEN,NSET=E2
 112,312,100
*NGEN,NSET=E3
 312,314,1
*NGEN,NSET=E4
 314,514,100
*NGEN,NSET=E5
 514,518,1
*NGEN,NSET=E6
 518,718,100
*NGEN,NSET=E7
 701,718,1
*NGEN,NSET=E8
 101,701,100
*NGEN,NSET=INTERIOR
 201,212,1
 301,314,1
 401,414,1
 501,518,1
 601,618,1
*ELEMENT,TYPE=S4R,ELSET=PLATE
 101, 101,201,202,102
 301, 301,401,402,302
 501, 501,601,602,502
*ELGEN,ELSET=PLATE
 101,11,1,1,2,100,100
 301,13,1,1,2,100,100
 501,17,1,1,2,100,100
*SHELL SECTION,ELSET=PLATE,MATERIAL=MAT
 1.0,
*MATERIAL,NAME=MAT
*ELASTIC
 1.0,0.0
*NODE,NSET=REF
 1000, 0.,0.,0.
*RIGID BODY, REFNODE=1000, ELSET=PLATE
*************************************************
*NODE,NSET=TOP
 1303,  2., 2., 2.
 1311,  2.,10., 2.
 1503,  4., 2., 2.
 1513,  4.,12., 2.
 1703,  6., 2., 2.
 1717,  6.,16., 2.
*NGEN,NSET=TOP
 1303,1311,2
 1503,1513,2
 1703,1717,2
*******************************************
** put mass and rotary inertia at top nodes
*******************************************
*ELEMENT, TYPE=MASS, ELSET=MASS
 1011,1303
 1021,1503
 1031,1703
*ELGEN, ELSET=MASS
 1011,5,2,1
 1021,6,2,1
 1031,8,2,1
*ELEMENT, TYPE=MASS, ELSET=MASS
 1041,10305
*MASS, ELSET=MASS
 12.0
*ELEMENT, TYPE=ROTARYI, ELSET=ROTI
 2011,1303
 2021,1503
 2031,1703
*ELGEN, ELSET=ROTI
 2011,5,2,1
 2021,6,2,1
 2031,8,2,1
*ELEMENT, TYPE=ROTARYI, ELSET=ROTI
 2041, 10507
*ROTARY INERTIA, ELSET=ROTI
 120.0, 120.0, 120.0
********************************
** define the connector elements
********************************
*ELEMENT, TYPE=CONN3D2, ELSET=AX
 11, 202, 1303
*ELEMENT, TYPE=CONN3D2, ELSET=AXBC
** prevent rigid body motion in static step
 9911, 202, 1303
*ELEMENT, TYPE=CONN3D2, ELSET=CAR
 12, 204, 1305
*NODE
10305,1,4,2
*ELEMENT, TYPE=CONN3D2, ELSET=PROJCAR
 1201, 104, 10305
*NODE
** Special definition of node 1307 for radial-thrust
** so that the radial direction is properly defined
** 1307 = (2,6,2) + (0,-1,1)/sqrt(2)
 1307, 2.0, 5.29289322, 2.70710678
*ELEMENT, TYPE=CONN3D2, ELSET=RTH
 13, 206, 1307
*ELEMENT, TYPE=CONN3D2, ELSET=RTHBC
** prevent rigid body motion in static step
 9913, 206, 1307
*ELEMENT, TYPE=CONN3D2, ELSET=SLP
 14, 208, 1309
*ELEMENT, TYPE=CONN3D2, ELSET=SLT
 15, 210, 1311
*ELEMENT, TYPE=CONN3D2, ELSET=CDN
 21, 402, 1503
*ELEMENT, TYPE=CONN3D2, ELSET=EUL
 22, 404, 1505
*ELEMENT, TYPE=CONN3D2, ELSET=FLT
 23, 406, 1507
*NODE
10507,5,6,2
*ELEMENT, TYPE=CONN3D2, ELSET=PROJFLT
 1023, 506, 10507
*ELEMENT, TYPE=CONN3D2, ELSET=REV
 24, 408, 1509
*ELEMENT, TYPE=CONN3D2, ELSET=ROT
 25, 410, 1511
*ELEMENT, TYPE=CONN3D2, ELSET=UNI
 26, 412, 1513
*ELEMENT, TYPE=CONN3D2, ELSET=AX-EUL
 31, 602, 1703
*ELEMENT, TYPE=CONN3D2, ELSET=AXBC
** prevent rigid body motion in static step
 9931, 602, 1703
*ELEMENT, TYPE=CONN3D2, ELSET=CAR-CDN
 32, 604, 1705
*ELEMENT, TYPE=CONN3D2, ELSET=CAR-ROT
 33, 606, 1707
*ELEMENT, TYPE=CONN3D2, ELSET=JN-FLT
 34, 608, 1709
*NODE
** define a special node for link so that bc's can
** be defined to keep the node from moving
 1711, 5., 9., 2.
*BOUNDARY
 1711, 1,2
*ELEMENT, TYPE=CONN3D2, ELSET=LNK-UNI
 35, 610, 1711
*NODE
** Special definition of node 1713 for radial-thrust
** so that the radial direction is properly defined
** 1713 = (6,12,2) + (0,-1,1)/sqrt(2)
 1713, 6.0, 11.29289322, 2.70710678
*ELEMENT, TYPE=CONN3D2, ELSET=RTH-CV
 36, 612, 1713
*ELEMENT, TYPE=CONN3D2, ELSET=RTHBC
** prevent rigid body motion in static step
 9936, 612, 1713
*ELEMENT, TYPE=CONN3D2, ELSET=SLP-REV
 37, 614, 1715
*ELEMENT, TYPE=CONN3D2, ELSET=SLT-ALN
 38, 616, 1717
*ELSET, ELSET=CONNALL
 AX, CAR, PROJCAR, RTH, SLP, SLT
 CDN, EUL, FLT, PROJFLT, REV, ROT, UNI
 AX-EUL, CAR-CDN, CAR-ROT, JN-FLT
 LNK-UNI, RTH-CV, SLP-REV, SLT-ALN
*ORIENTATION, NAME=GLOBAL
** Global directions
 1.,0.,0., 0.,1.,0.
*ORIENTATION, NAME=LDIRECTION
** orientation used by most connectors
 1.,1.,1., 0.,-1.,1.
*ORIENTATION, NAME=LDIRECTION_PERT
 1.,1.,1., 0.001,-1.,1.
*ORIENTATION, NAME=ORI_RAD-THRUST
** 1-2-3 for radial-thrust is 2-3-1 for others.
 0.,-1.,1., 2.,-1.,-1.
*ORIENTATION, NAME=ORI_UNI
** 1-2-3 for universal is 2-1-(-3) for others.
 0.,-1.,1., 1.,1.,1.
*ORIENTATION, NAME=L45
 1.,0.,0., 0.,1.,1.
*****
*CONNECTOR SECTION, ELSET=AX, BEHAVIOR=CB_AX
 AXIAL,
*CONNECTOR BEHAVIOR, NAME=CB_AX
*CONNECTOR ELASTICITY, COMPONENT=1
 48.0,
*CONNECTOR DAMPING, COMPONENT=1
 5.0,
*****
*CONNECTOR SECTION, ELSET=AXBC
 SLOT,
 LDIRECTION,
*****
*CONNECTOR SECTION, ELSET=CAR, BEHAVIOR=CB_CAR
 CARTESIAN,
 LDIRECTION,
*CONNECTOR SECTION, ELSET=PROJCAR, BEHAVIOR=CB_CAR
 PROJECTION CARTESIAN,
 LDIRECTION,LDIRECTION_PERT
*CONNECTOR BEHAVIOR, NAME=CB_CAR
*CONNECTOR ELASTICITY, COMPONENT=1
 48.0,
*CONNECTOR DAMPING, COMPONENT=1
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=2
 48.0,
*CONNECTOR DAMPING, COMPONENT=2
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=3
 48.0,
*CONNECTOR DAMPING, COMPONENT=3
 5.0,
*****
*CONNECTOR SECTION, ELSET=RTH, BEHAVIOR=CB_RTH
 RADIAL-THRUST,
 ORI_RAD-THRUST,
*CONNECTOR BEHAVIOR, NAME=CB_RTH
*CONNECTOR ELASTICITY, COMPONENT=1
 48.0,
*CONNECTOR DAMPING, COMPONENT=1
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=3
 48.0,
*CONNECTOR DAMPING, COMPONENT=3
 5.0,
*****
*CONNECTOR SECTION, ELSET=RTHBC, BEHAVIOR=CB_RTHBC
 CARTESIAN,
 ORI_RAD-THRUST,
*CONNECTOR BEHAVIOR, NAME=CB_RTHBC
*CONNECTOR ELASTICITY, COMPONENT=2
 48.0E2,
*****
*CONNECTOR SECTION, ELSET=SLP, BEHAVIOR=CB_SLP
 SLIDE-PLANE,
 LDIRECTION,
*CONNECTOR BEHAVIOR, NAME=CB_SLP
*CONNECTOR ELASTICITY, COMPONENT=2
 48.0,
*CONNECTOR DAMPING, COMPONENT=2
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=3
 48.0,
*CONNECTOR DAMPING, COMPONENT=3
 5.0,
*****
*CONNECTOR SECTION, ELSET=SLT, BEHAVIOR=CB_SLT
 SLOT,
 LDIRECTION,
*CONNECTOR BEHAVIOR, NAME=CB_SLT
*CONNECTOR ELASTICITY, COMPONENT=1
 48.0,
*CONNECTOR DAMPING, COMPONENT=1
 5.0,
*****
*CONNECTOR SECTION, ELSET=CDN, BEHAVIOR=CB_CDN
** Make CARDAN behave like UNIVERSAL
 CARDAN,
 ORI_UNI,
*CONNECTOR BEHAVIOR, NAME=CB_CDN
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=5
 1.0E6,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.0,
*CONNECTOR DAMPING, COMPONENT=6
 50.0,
*****
*CONNECTOR SECTION, ELSET=EUL, BEHAVIOR=CB_EUL
 EULER,
 GLOBAL,L45
*CONNECTOR BEHAVIOR, NAME=CB_EUL
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=5
 480.E2,
*CONNECTOR DAMPING, COMPONENT=5
 50.E2,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.E2,
*CONNECTOR DAMPING, COMPONENT=6
 50.E2,
*****
*CONNECTOR SECTION, ELSET=FLT, BEHAVIOR=CB_FLT
 FLEXION-TORSION,
 GLOBAL,L45
*CONNECTOR BEHAVIOR, NAME=CB_FLT
*CONNECTOR ELASTICITY, COMPONENT=4
 480.E5,
*CONNECTOR DAMPING, COMPONENT=4
 50.E2,
*CONNECTOR ELASTICITY, COMPONENT=5
 480.0,
*CONNECTOR DAMPING, COMPONENT=5
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.d5,
**
*CONNECTOR SECTION, ELSET=PROJFLT, BEHAVIOR=CB_PROJFLT
 PROJECTION FLEXION-TORSION,
 LDIRECTION,LDIRECTION_PERT
*CONNECTOR BEHAVIOR, NAME=CB_PROJFLT
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=5
 480.0,
*CONNECTOR DAMPING, COMPONENT=5
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.0,
*CONNECTOR DAMPING, COMPONENT=6
 50.0,
*****
*CONNECTOR SECTION, ELSET=REV, BEHAVIOR=CB_REV
 REVOLUTE,
 LDIRECTION,
*CONNECTOR BEHAVIOR, NAME=CB_REV
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*****
*CONNECTOR SECTION, ELSET=ROT, BEHAVIOR=CB_ROT
 ROTATION,
 LDIRECTION,
*CONNECTOR BEHAVIOR, NAME=CB_ROT
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=5
 480.0,
*CONNECTOR DAMPING, COMPONENT=5
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.0,
*CONNECTOR DAMPING, COMPONENT=6
 50.0,
*****
*CONNECTOR SECTION, ELSET=UNI, BEHAVIOR=CB_UNI
 UNIVERSAL,
 ORI_UNI,
*CONNECTOR BEHAVIOR, NAME=CB_UNI
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.0,
*CONNECTOR DAMPING, COMPONENT=6
 50.0,
*****
*CONNECTOR SECTION, ELSET=AX-EUL, BEHAVIOR=CB_AX-EUL
 EULER, AXIAL
 GLOBAL, L45
*CONNECTOR BEHAVIOR, NAME=CB_AX-EUL
*CONNECTOR ELASTICITY, COMPONENT=1
 48.0,
*CONNECTOR DAMPING, COMPONENT=1
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=5
 480.E2,
*CONNECTOR DAMPING, COMPONENT=5
 50.E2,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.E2,
*CONNECTOR DAMPING, COMPONENT=6
 50.E2,
*****
*CONNECTOR SECTION, ELSET=CAR-CDN, BEHAVIOR=CB_CAR-CDN
 CARTESIAN,CARDAN
 ORI_UNI,
*CONNECTOR BEHAVIOR, NAME=CB_CAR-CDN
*CONNECTOR ELASTICITY, COMPONENT=1
 48.0,
*CONNECTOR DAMPING, COMPONENT=1
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=2
 48.0,
*CONNECTOR DAMPING, COMPONENT=2
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=3
 48.0,
*CONNECTOR DAMPING, COMPONENT=3
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=5
 1.0E6,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.0,
*CONNECTOR DAMPING, COMPONENT=6
 50.0,
*****
*CONNECTOR SECTION, ELSET=CAR-ROT, BEHAVIOR=CB_CAR-ROT
 CARTESIAN,ROTATION
 LDIRECTION,
*CONNECTOR BEHAVIOR, NAME=CB_CAR-ROT
*CONNECTOR ELASTICITY, COMPONENT=1
 48.0,
*CONNECTOR DAMPING, COMPONENT=1
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=2
 48.0,
*CONNECTOR DAMPING, COMPONENT=2
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=3
 48.0,
*CONNECTOR DAMPING, COMPONENT=3
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=5
 480.0,
*CONNECTOR DAMPING, COMPONENT=5
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.0,
*CONNECTOR DAMPING, COMPONENT=6
 50.0,
*****
*CONNECTOR SECTION, ELSET=JN-FLT, BEHAVIOR=CB_FLT
 JOIN, FLEXION-TORSION
 GLOBAL,L45
*****
*CONNECTOR SECTION, ELSET=LNK-UNI, BEHAVIOR=CB_LNK-UNI
 LINK,UNIVERSAL
 ORI_UNI,
*CONNECTOR BEHAVIOR, NAME=CB_LNK-UNI
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.0,
*CONNECTOR DAMPING, COMPONENT=6
 50.0,
*****
*CONNECTOR SECTION, ELSET=RTH-CV, BEHAVIOR=CB_RTH
 RADIAL-THRUST,CONSTANT VELOCITY
 ORI_RAD-THRUST,
*****
*CONNECTOR SECTION, ELSET=SLP-REV, BEHAVIOR=CB_SLP-REV
 SLIDE-PLANE,REVOLUTE
 LDIRECTION,
*CONNECTOR BEHAVIOR, NAME=CB_SLP-REV
*CONNECTOR ELASTICITY, COMPONENT=2
 48.0,
*CONNECTOR DAMPING, COMPONENT=2
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=3
 48.0,
*CONNECTOR DAMPING, COMPONENT=3
 5.0,
*CONNECTOR ELASTICITY, COMPONENT=4
 480.0,
*CONNECTOR DAMPING, COMPONENT=4
 50.0,
*****
*CONNECTOR SECTION, ELSET=SLT-ALN, BEHAVIOR=CB_SLT-ALN
 SLOT,ALIGN
 LDIRECTION,
*CONNECTOR BEHAVIOR, NAME=CB_SLT-ALN
*CONNECTOR ELASTICITY, COMPONENT=1
 48.0,
*CONNECTOR DAMPING, COMPONENT=1
 5.0,
*****
*BOUNDARY
 REF, 1,6
*****
*ELSET,ELSET=D100
** initial displacement in local (1,0,0) direction.
 AX, SLT, AX-EUL, SLT-ALN
*ELSET,ELSET=D111
** initial displacement in local (1,1,1) direction.
 CAR, CAR-ROT, CAR-CDN, PROJCAR
*ELSET,ELSET=D101
** initial displacement in RADIAL-THRUST local (1,0,1) direction.
 RTH, RTH-CV
*ELSET,ELSET=D011
** initial diplacement in local (0,1,1) direction.
 SLP, SLP-REV
*ELSET,ELSET=R100
** initial rotation in local (1,0,0) direction.
 REV, SLP-REV
*ELSET,ELSET=R111
** initial rotation in local (1,1,1) direction.
 ROT, CAR-ROT, PROJFLT
*ELSET,ELSET=R101
** initial rotation in local (1,0,1) direction.
 UNI,  LNK-UNI, CDN, CAR-CDN
*ELSET,ELSET=REUL
 EUL, AX-EUL
*ELSET,ELSET=RFLT
 FLT, JN-FLT
*ELSET,ELSET=RPROJFLT
 PROJFLT
************************************
**FOR PROJECTION CARTESIAN TO PREVENT NUMERICAL SINGULARITY
*BOUNDARY
10305,4,6,0
*STEP, PERT
 This step should produce the same results as the linear
 static step which follows
*STATIC
*CONNECTOR MOTION
 D100,1, 0.1
 D111,1, 0.1
 D111,2, 0.1
 D111,3, 0.1
 D101,1, 0.1
 D101,3, 0.1
 D011,2, 0.1
 D011,3, 0.1
 R100,4, 0.1
 R111,4, 0.1
 R111,5, 0.1
 R111,6, 0.1
 R101,4, 0.1
 R101,6, 0.1
 REUL,4, 0.1
 RFLT,5, 0.1
*EL FILE, ELSET=CONNALL
 CP,CU,CV
*OUTPUT,HISTORY, FREQ=1
*ELEMENT OUTPUT, ELSET=CONNALL
 CP,CU,CV
*NODE OUTPUT, NSET=TOP
 U,V,A
*END STEP
************************************
*STEP, NLGEOM=NO
 Every dof should give the same answer.
*STATIC
*CONNECTOR MOTION
 D100,1, 0.1
 D111,1, 0.1
 D111,2, 0.1
 D111,3, 0.1
 D101,1, 0.1
 D101,3, 0.1
 D011,2, 0.1
 D011,3, 0.1
 R100,4, 0.1
 R111,4, 0.1
 R111,5, 0.1
 R111,6, 0.1
 R101,4, 0.1
 R101,6, 0.1
 REUL,4, 0.1
 RFLT,5, 0.1
*EL FILE, ELSET=CONNALL, FREQ=10
 CP,CU,CV
*OUTPUT,HISTORY, FREQ=1
*ELEMENT OUTPUT, ELSET=CONNALL
 CP,CU,CV
*NODE OUTPUT, NSET=TOP
 U,V,A
*END STEP
************************************
*STEP, NLGEOM=NO, INC=200
 Damped free vibration.
*DYNAMIC, ALPHA=0.0, NOHAF
 0.075, 8.0
*CONNECTOR MOTION, OP=NEW
*FILE FORMAT, ZERO INC
*EL FILE, ELSET=CONNALL, FREQ=10
 CP,CU,CV
*OUTPUT,HISTORY, FREQ=1
*ELEMENT OUTPUT, ELSET=CONNALL
 CP,CU,CV
*NODE OUTPUT, NSET=TOP
 U,V,A
*END STEP