*HEADING
 BASIC VERIFICATION OF *COUPLING AND *DISTRIBUTING
 3-D NONLINEAR GEOMETRIC ANALYSIS
*PREPRINT, MODEL=YES,ECHO=YES,HISTORY=YES
*RESTART,W
*NODE
10,0.0,  0.0
1 ,2.0,  0.5
2 ,2.0, -0.5
3 ,1.0,  1.0
*SURFACE, NAME=NODE_BASED_SURFACE1, TYPE=NODE
1,1.0
2,2.0
3,3.0
*COUPLING, CONSTRAINT NAME=COUPLING_DEFINITION_1, REF NODE=10, SURFACE=NODE_BASED_SURFACE1
*DISTRIBUTING
*ORIENTATION,SYS=RECT,NAME=ORIGLOB
1.,0.,0.,0.,1.,0.
*ELEMENT,TYPE=CONN3D2,ELSET=SPRINGS
401,1,2
402,2,3
404,3,1
*CONNECTOR SECTION, ELSET=SPRINGS, BEHAVIOR=SPRINGS
CARTESIAN,
ORIGLOB,
*CONNECTOR BEHAVIOR, NAME=SPRINGS
*CONNECTOR ELASTICITY,COMPONENT=1
1.0E8, 
*CONNECTOR ELASTICITY,COMPONENT=2
1.0E8, 
*CONNECTOR ELASTICITY,COMPONENT=3
1.0E8, 
*ELEMENT,TYPE=CONN3D2,ELSET=DASHPOT1
110,1, ,
*CONNECTOR SECTION,ELSET=DASHPOT1,BEHAVIOR=DASHPOT1
CARTESIAN,
ORIGLOB,
*CONNECTOR BEHAVIOR, NAME=DASHPOT1
*CONNECTOR DAMPING,COMPONENT=1
1.0E2, 
*CONNECTOR DAMPING,COMPONENT=2
2.0E2, 
*CONNECTOR DAMPING,COMPONENT=3
3.0E2, 
*ELEMENT,TYPE=CONN3D2,ELSET=DASHPOT2
120,2, ,
*CONNECTOR SECTION,ELSET=DASHPOT2,BEHAVIOR=DASHPOT2
CARTESIAN,
ORIGLOB,
*CONNECTOR BEHAVIOR, NAME=DASHPOT2
*CONNECTOR DAMPING,COMPONENT=1
1.0E2, 
*CONNECTOR DAMPING,COMPONENT=2
2.0E2, 
*CONNECTOR DAMPING,COMPONENT=3
3.0E2, 
*ELEMENT,TYPE=CONN3D2,ELSET=DASHPOT3
130,3, ,
*CONNECTOR SECTION,ELSET=DASHPOT3,BEHAVIOR=DASHPOT3
CARTESIAN,
ORIGLOB,
*CONNECTOR BEHAVIOR, NAME=DASHPOT3
*CONNECTOR DAMPING,COMPONENT=1
1.0E2, 
*CONNECTOR DAMPING,COMPONENT=2
2.0E2, 
*CONNECTOR DAMPING,COMPONENT=3
3.0E2, 
*ELEMENT, TYPE=MASS, ELSET=MASSREF
100, 10
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARYREF
200, 10
*ELEMENT, TYPE=MASS, ELSET=MASS_1
 1001, 1
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY_1
 2001, 1
*ELEMENT, TYPE=MASS, ELSET=MASS_2
 1002, 2
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY_2
 2002, 2
*ELEMENT, TYPE=MASS, ELSET=MASS_3
 1003, 3
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY_3
 2003, 3
*MASS, ELSET=MASS_1
 1.66666667,
*MASS, ELSET=MASS_2
 3.33333333,
*MASS, ELSET=MASS_3
 5.0,
*MASS, ELSET=MASSREF
.01,
*ROTARY INERTIA, ELSET=ROTARY_1
 16.6666667,
*ROTARY INERTIA, ELSET=ROTARY_2
 33.3333333,
*ROTARY INERTIA, ELSET=ROTARY_3
 50,
*ROTARY INERTIA, ELSET=ROTARYREF
.1, .1, .1
*BOUNDARY
10,2,5
*AMPLITUDE,NAME=SINE,DEFINITION=PERIODIC
1,6.28,0.0,0.0
0.0,1.0
*AMPLITUDE,NAME=AMP,DEFINITION=SMOOTH STEP
0.0,0.0,1.0,1.0
*NSET, NSET=QA_TEST_ALLNODES
1,2,3,10
*NSET, NSET=OUT
1,2,3,10
**
**  STEP 1: ROTATION IN THE X-Y PLANE
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT
5.E-5,1.0
*BOUNDARY, AMPLITUDE=AMP
10,6,6,6.283185308
*BOUNDARY,AMPLITUDE=SINE
10,1,1,1.0
*OUTPUT,FIELD,NUMBER=4,TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U,
***OUTPUT,HISTORY,FREQ=1
****NODE OUTPUT,NSET=OUT
***U,V,A,RF
****ENERGY OUTPUT
***ALLIE,ALLKE,ALLWK,ETOTAL,ALLAE
*END STEP
**
**  STEP 2: ROTATION IN THE Z-X PLANE
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT
5.E-5,1.0
*BOUNDARY,TYPE=VELOCITY,OP=NEW
10,5,5,6.283185308
*BOUNDARY,AMPLITUDE=SINE,OP=NEW
10,3,3,1.0
10,1,2,0.0
10,4,4,0.0
10,6,6,0.0
*END STEP
**
**  STEP 3: ROTATION IN THE Y-Z PLANE
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT
5.E-5,1.0
*BOUNDARY,TYPE=VELOCITY,OP=NEW
10,4,4,6.283185308
*BOUNDARY,AMPLITUDE=SINE,OP=NEW
10,2,2,1.0
10,1,1,0.0
10,3,3,0.0
10,5,6,0.0
*END STEP
**
**  STEP 4: 
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT
5.E-5,1.0
*BOUNDARY,OP=NEW,TYPE=VELOCITY
10,6,6,6.283185308
*BOUNDARY,OP=NEW,AMPLITUDE=SINE
10,1,1,1.0
10,2,5,0.0
*END STEP