*HEADING
CONN_STD_FORCE_2D: 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
*part,name=PART-1
*NODE,NSET=fixed
202,1,1,0
204,1,3,0
206,1,5,0
402,3,1,0
602,5,1,0
604,5,3,0
****************
*NODE,NSET=TOP
1303,  2., 2.
1305,  2., 4.
1307,  2., 6.
1503,  4., 2.
1703,  6., 2.
1705,  6., 4.
********************************
** 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 SECTION, ELSET=CAR, BEHAVIOR=CB_CAR
CARTESIAN,
LDIRECTION,
*****
*CONNECTOR SECTION, ELSET=SLT, BEHAVIOR=CB_SLT
SLOT,
LDIRECTION,
*****
*CONNECTOR SECTION, ELSET=ROT, BEHAVIOR=CB_ROT
ROTATION,
LDIRECTION,
*****
*CONNECTOR SECTION, ELSET=CAR-ROT, BEHAVIOR=CB_CAR-ROT
CARTESIAN,ROTATION
LDIRECTION,
*****
*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
*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
************************************
*ELSET,ELSET=QA_TEST
CONNALL
*end part
*Assembly,name=ASSEMBLY-1
*Instance,name=PART-1-1,part=PART-1
*end instance
*******************************************
** put mass and rotary inertia at top nodes
*******************************************
*ELEMENT, TYPE=MASS, ELSET=MASS
1011,PART-1-1.1303
1012,PART-1-1.1305
1013,PART-1-1.1307
1021,PART-1-1.1503
1031,PART-1-1.1703
1032,PART-1-1.1705
*MASS, ELSET=MASS
12.0
*ELEMENT, TYPE=ROTARYI, ELSET=ROTI
2011,PART-1-1.1303
2012,PART-1-1.1305
2013,PART-1-1.1307
2021,PART-1-1.1503
2031,PART-1-1.1703
2032,PART-1-1.1705
*ROTARY INERTIA, ELSET=ROTI
120.0, 120.0, 120.0
*End assembly
*** 
*CONNECTOR BEHAVIOR, NAME=CB_AX
*CONNECTOR ELASTICITY, COMPONENT=1
48.0,
*CONNECTOR DAMPING, COMPONENT=1
5.0,
*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 BEHAVIOR, NAME=CB_SLT
*CONNECTOR ELASTICITY, COMPONENT=1
48.0,
*CONNECTOR DAMPING, COMPONENT=1
5.0,
*CONNECTOR BEHAVIOR, NAME=CB_ROT
*CONNECTOR ELASTICITY, COMPONENT=6
480.0,
*CONNECTOR DAMPING, COMPONENT=6
50.0,
*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,
*BOUNDARY
ASSEMBLY-1.PART-1-1.fixed,1,6
*BOUNDARY
ASSEMBLY-1.PART-1-1.TOP, 3,5
*INITIAL CONDITIONS, TYPE=VELOCITY
** each displacement component has initial velocity 0.96
** in all the available component of relative motion.
ASSEMBLY-1.PART-1-1.D10,1, 0.67882251
ASSEMBLY-1.PART-1-1.D10,2, 0.67882251
**
ASSEMBLY-1.PART-1-1.D11,1, 0.0
ASSEMBLY-1.PART-1-1.D11,2, 1.35764502
**
ASSEMBLY-1.PART-1-1.R001,6, 0.96
***
*AMPLITUDE, NAME=COSINE, DEF=PERIODIC
1, 2.0, 0.0, 0.0
1.0, 0.0
*STEP, INC=200, NLGEOM, UNSYM=YES
Force steady-state vibration.  Every dof should give the same answer.
*DYNAMIC, ALPHA=0.0, NOHAF, DIRECT
0.075, 8.0
*CONNECTOR LOAD, AMP=COSINE
ASSEMBLY-1.PART-1-1.CONN1,1, 4.8
ASSEMBLY-1.PART-1-1.CONN2,2, 4.8
ASSEMBLY-1.PART-1-1.CONN6,6, 48.0
*OUTPUT,HISTORY, TIME INTERVAL=0.2
*ELEMENT OUTPUT, ELSET=ASSEMBLY-1.PART-1-1.QA_TEST
CP,CU,CV,CA,CCU,CTF,CVF,CEF,CUE
*NODE OUTPUT, NSET=ASSEMBLY-1.PART-1-1.TOP
U,V,A
*Restart,write,frequency=1
*END STEP