*HEADING
 CONN_XPL_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
*DENSITY
 1.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
*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
*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=CAR
 12, 204, 1305
*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=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
*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=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=SLP-REV
 37, 614, 1715
*ELEMENT, TYPE=CONN3D2, ELSET=SLT-ALN
 38, 616, 1717
*ELSET, ELSET=CONNALL
 AX, CAR, RTH, SLP, SLT
 CDN, EUL, FLT, 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=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 CONSTITUTIVE REFERENCE
 1.63205081,
*****
*CONNECTOR SECTION, ELSET=CAR, BEHAVIOR=CB_CAR
 CARTESIAN,
 LDIRECTION,
*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 CONSTITUTIVE REFERENCE
 1.63205081,-0.1,-0.1
*****
*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 CONSTITUTIVE REFERENCE
 0.9, ,1.63205081
*****
*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 CONSTITUTIVE REFERENCE
 ,-0.1,-0.1
*****
*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 CONSTITUTIVE REFERENCE
 1.63205081,
*****
*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 CONSTITUTIVE REFERENCE
 , , ,-5.72957795, ,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 , , ,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 , , , ,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 , , ,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 , , ,-5.72957795,-5.72957795,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 , , ,-5.72957795, ,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 1.63205081, , ,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 -0.1,1.63205081,-0.1,-5.72957795, ,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 1.63205081,-0.1,-0.1,-5.72957795,-5.72957795,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 , , ,-5.72957795, ,-5.72957795
*****
*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 CONSTITUTIVE REFERENCE
 ,-0.1,-0.1,-5.72957795
*****
*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,
*CONNECTOR CONSTITUTIVE REFERENCE
 1.63205081,
*****
*BOUNDARY
 REF, 1,6
************************************
*ELSET,ELSET=QA_TEST_CONNALL
CONNALL,
**
*STEP, NLGEOM=YES
 Damped free vibration.
*DYNAMIC, EXPLICIT, DIRECT
 0.0075, 8.0
*FILE OUTPUT, NUM=4
*EL FILE, ELSET=CONNALL
 CP,CU,
*OUTPUT,HISTORY, FREQ=1
*ELEMENT OUTPUT, ELSET=CONNALL
 CP,CU,CCU,CV
*NODE OUTPUT, NSET=TOP
 U,V,A
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=2
***OUTPUT,HISTORY,VARIABLE=PRESELECT
*END STEP