*HEADING
 CONN_STD_FORCE_3D: forced, steady-state oscillation  [CONN3D2]
 Steady-state damped motion:
  u(t) = p_0/(2 xi k) sin(Omega*t) = 0.48 sin(2t)
  u(0) = 0.0, v(0) = p_0 Omega/(2 xi k)
  f(t) = p_0 cos(Omega*t)
  Omega = w_n         = 2.0
  xi    = c/(2 m w_n) = 5/48
  p_0   = k/10
*************************************************
** 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, 10305
*ELEMENT, TYPE=ROTARYI, ELSET=ROTI
 2051, 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=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=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=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, 617, 1717
*ORIENTATION, NAME=GLOBAL
** Global directions
 1.,0.,0., 0.,1.,0.
*ORIENTATION, NAME=LDIRECTION
** orientation used by most connectors
 1.,1.,1., 0.001,-1.,1.
*ORIENTATION, NAME=LDIRECTION_PERT
** 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 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=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.d5,
*CONNECTOR DAMPING, COMPONENT=5
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.d5,
*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.d5,
*CONNECTOR DAMPING, COMPONENT=5
 50.0,
*CONNECTOR ELASTICITY, COMPONENT=6
 480.d5,
*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,
********************************************
** Define elsets for elements with available
** components in the directions indicated:
********************************************
*ELSET,ELSET=CONN1
 AX, CAR, PROJCAR, RTH, SLT
 AX-EUL, CAR-CDN, CAR-ROT, RTH-CV, SLT-ALN
*ELSET,ELSET=CONN2
 CAR, PROJCAR, SLP
 CAR-CDN, CAR-ROT, SLP-REV
*ELSET,ELSET=CONN3
 CAR, PROJCAR, RTH, SLP
 CAR-CDN, CAR-ROT, RTH-CV, SLP-REV
*ELSET,ELSET=CONN4
 CDN, EUL, REV, ROT, UNI, PROJFLT
 AX-EUL, CAR-CDN, CAR-ROT, LNK-UNI, SLP-REV
*ELSET,ELSET=CONN5
 FLT,
 JN-FLT,PROJFLT
*ELSET,ELSET=CONN6
 CDN, UNI 
 CAR-CDN,  LNK-UNI,PROJFLT
*ELSET,ELSET=CONNALL
 CONN1,CONN2,CONN3,CONN4,CONN5,CONN6
*BOUNDARY
 REF, 1,6
*AMPLITUDE, NAME=COSINE, DEF=PERIODIC
 1, 2.0, 0.0, 0.0
 1.0, 0.0
*NSET,NSET=D100
** initial velocity in local (1,0,0) direction.
** AXIAL, SLOT, AX-EUL, SLT-ALN
 1303, 1311, 1703, 1717
*NSET,NSET=D111
** initial velocity in local (1,1,1) direction.
** CARTESIAN, CAR-ROT
 1305, 1707, 10305
*NSET,NSET=D101
** initial velocity in RADIAL-THRUST local (1,0,1) direction.
** RADIAL-THRUST, RTH-CV
 1307, 1713
*NSET,NSET=D011
** initial velocity in local (0,1,1) direction.
** SLIDE-PLANE, SLP-REV
 1309, 1715
*NSET,NSET=DUNI
** initial velocity for ORI_UNI orientation
** CAR-CDN
 1705,
*NSET,NSET=R100
** initial angular velocity in local (1,0,0) direction.
** REVOLUTE, SLP-REV
 1509, 1715
*NSET,NSET=R101
** UNIVERSAL, CARDAN, CAR-CDN, LNK-UNI
 1503, 1513, 1705, 1711
*NSET,NSET=R111
** initial angular velocity in local (1,1,1) direction.
** ROTATION, CAR-ROT
 1511, 1707
*NSET,NSET=PROJFLT
 10507
*NSET,NSET=REUL
** initial angular velocity in global (0,0,1) direction.
** EULER, AX-EUL
 1505, 1703
*NSET,NSET=RFLT
** initial angular velocity in global (0,-1,1) direction.
** FLT, JN-FLT
 1507, 1709
*INITIAL CONDITIONS, TYPE=VELOCITY
** each displacement component has initial velocity 0.96
** in all the available component of relative motion.
 D100,1, 0.5542563
 D100,2, 0.5542563
 D100,3, 0.5542563
**
 D111,1, 1.33809298
 D111,2,-0.51648461
 D111,3, 0.84116041
**
 D101,1, 0.5542563
 D101,2,-0.12456625
 D101,3, 1.23307877
**
 D011,1, 0.78383672
 D011,2,-1.07074087
 D011,3, 0.28690415
**
 DUNI,1,-0.22958046
 DUNI,2, 0.26735211
 DUNI,3, 1.62499713
**
 R100,4, 0.5542563
 R100,5, 0.5542563
 R100,6, 0.5542563
**
 PROJFLT,4, 1.33809298
 PROJFLT,5,-0.51648461
 PROJFLT,6, 0.84116041
 R111,4, 0.5542563
 R111,5, 0.5542563
 R111,6, 0.5542563
**
 R101,4, -0.7838367176
 R101,5, -0.2869041514
 R101,6,  1.070740869
**
 REUL,4, 0.0
 REUL,5, 0.0
 REUL,6, 0.96
**
 RFLT,4, 0.0
 RFLT,5,-0.67882251
 RFLT,6, 0.67882251
************************************
**FIX ROTATION FOR PROJECTION CARTESIAN
*BOUNDARY
10305,4,6,0
*STEP, INC=200, NLGEOM, UNSYM=NO
 Force steady-state vibration.  Every dof should give the same answer.
*DYNAMIC, ALPHA=0.0, NOHAF
 0.075, 8.0
*CONNECTOR LOAD, AMP=COSINE
 CONN1,1, 4.8
 CONN2,2, 4.8
 CONN3,3, 4.8
 CONN4,4, 48.0
 CONN5,5, 48.0
 CONN6,6, 48.0
*FILE FORMAT, ZERO INC
*EL FILE, ELSET=CONNALL, FREQ=10
 CP,CU,CV,CA
*OUTPUT,HISTORY, FREQ=1
*ELEMENT OUTPUT, ELSET=CONNALL
 CP,CU,CV,CA
*NODE OUTPUT, NSET=TOP
 U,V,A
*END STEP