*HEADING
 CONN_STD_FORCE_AXI: forced, steady-state oscillation  [CONN2D2]
 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.
 107,  0.,  6.
 701,  6.,  0.
 707,  6.,  6.
*NGEN,NSET=E1
 101,107,1
*NGEN,NSET=E2
 107,707,100
*NGEN,NSET=E3
 701,707,1
*NGEN,NSET=E4
 101,701,100
*NGEN,NSET=INTERIOR
 201,207,1
 301,307,1
 401,407,1
 501,507,1
 601,607,1
*NSET, NSET=NPLATE
 E1, E2, E3, E4, INTERIOR
*ELEMENT,TYPE=CAX4R,ELSET=PLATE
 101, 101,201,202,102
*ELGEN,ELSET=PLATE
 101,6,1,1,6,100,100
*SOLID SECTION,ELSET=PLATE,MATERIAL=MAT
*MATERIAL,NAME=MAT
*ELASTIC
 1.0,0.0
*NODE,NSET=REF
 1000, 0.,0.
*NSET,NSET=FIXED
202,204,206,402,602,604
*BOUNDARY
FIXED,1,6
****RIGID BODY, REFNODE=1000, ELSET=PLATE, TIE=NPLATE
*************************************************
*NODE,NSET=TOP
 1303,  2., 2.
 1305,  2., 4.
 1307,  2., 6.
 1503,  4., 2.
 1703,  6., 2.
 1705,  6., 4.
*******************************************
** put mass and rotary inertia at top nodes
*******************************************
*ELEMENT, TYPE=MASS, ELSET=MASS
 1011,1303
 1012,1305
 1013,1307
 1021,1503
 1031,1703
 1032,1705
*MASS, ELSET=MASS
 12.0
*ELEMENT, TYPE=ROTARYI, ELSET=ROTI
 2011,1303
 2012,1305
 2013,1307
 2021,1503
 2031,1703
 2032,1705
*ROTARY INERTIA, ELSET=ROTI
 120.0, 120.0, 120.0
********************************
** define the connector elements
********************************
*ELEMENT, TYPE=CONN2D2, ELSET=AX
 11, 202, 1303
*ELEMENT, TYPE=CONN2D2, ELSET=CAR
 12, 204, 1305
*ELEMENT, TYPE=CONN2D2, ELSET=SLT
 13, 206, 1307
*ELEMENT, TYPE=CONN2D2, ELSET=ROT
 21, 402, 1503
*ELEMENT, TYPE=CONN2D2, ELSET=CAR-ROT
 31, 602, 1703
*ELEMENT, TYPE=CONN2D2, ELSET=SLT-ALN
 32, 604, 1705
*ORIENTATION, NAME=LDIRECTION
** orientation along node line
 1.,1.,0., -1.,1.,0.
*****
*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 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 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=ROT, BEHAVIOR=CB_ROT
 ROTATION,
 LDIRECTION,
*CONNECTOR BEHAVIOR, NAME=CB_ROT
*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=6
 480.0,
*CONNECTOR DAMPING, COMPONENT=6
 50.0,
*****
*CONNECTOR SECTION, ELSET=SLT-ALN, BEHAVIOR=CB_SLT
 SLOT,ALIGN
 LDIRECTION,
********************************************
** Define elsets for elements with available
** components in the directions indicated:
********************************************
*ELSET,ELSET=CONN1
 AX, CAR, SLT
 CAR-ROT, SLT-ALN
*ELSET,ELSET=CONN2
 CAR,
 CAR-ROT,
*ELSET,ELSET=CONN6
 ROT,
 CAR-ROT,
*ELSET,ELSET=CONNALL
 CONN1,CONN2,CONN6
*BOUNDARY
** REF, 1,6
 TOP, 3,5
*AMPLITUDE, NAME=COSINE, DEF=PERIODIC
 1, 2.0, 0.0, 0.0
 1.0, 0.0
*NSET,NSET=D10
** initial velocity in local (1,0) direction.
** AXIAL, SLOT, SLT-ALN
 1303, 1307, 1705
*NSET,NSET=D11
** initial velocity in local (1,1) direction.
** CARTESIAN, CAR-ROT
 1305, 1703
*NSET,NSET=R001
** initial angular velocity in local (0,0,1) direction.
** ROTATION, CAR-ROT
 1503, 1703
*INITIAL CONDITIONS, TYPE=VELOCITY
** each displacement component has initial velocity 0.96
** in all the available component of relative motion.
 D10,1, 0.67882251
 D10,2, 0.67882251
**
 D11,1, 0.0
 D11,2, 1.35764502
**
 R001,6, 0.96
************************************
*STEP, INC=200, NLGEOM, UNSYM=YES
 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
 CONN6,6, 48.0
*FILE FORMAT, ZERO INC
*EL FILE, FREQ=10
 IVOL, 
 EVOL
 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