*HEADING
GLOBAL TRACKING FOR CYLINDER ANALYTICAL RIGID SURFACE
**
** Define a dummy truss element
**
*NODE, NSET=QA_TEST_ALLNODES
101, 0.,-6.
102, 20.,-6.
**
*ELEMENT, TYPE=T3D2, ELSET=dummy
101,101,102
*SOLID SECTION, ELSET=DUMMY, MATERIAL=MMASS
*MATERIAL, NAME=MMASS
*ELASTIC
1., 
*DENSITY
0.001, 
**
** Define point masses
**
*NODE, NSET=NHEAVY
1,0.5,1.0, 1.
2,0.5,1.0, 2.
3,0.5,1.0, 3.
4,0.5,1.0, 4.
5,0.5,1.0, 5.
6,0.5,1.0, 6.
7,0.5,1.0, 7.
8,0.5,1.0, 8.
9,0.5,1.0, 9.
10,0.5,1.5, 10.
11,0.5,1.5, 11.
12,0.5,1.5, 12.
*NODE, NSET=NLIGHT
13,0.5,1.5, 13.
14,0.5,1.5, 14.
15,0.5,1.5, 15.
**
*ELEMENT, TYPE=MASS, ELSET=MASS1
1,1
2,2
3,3
4,4
5,5
6,6
7,7
8,8
9,9
10,10
11,11
12,12
13,13
14,14
15,15
*MASS, ELSET=MASS1
1.0, 
**
*INITIAL CONDITIONS, TYPE=VELOCITY
1,1, 1.0
2,1, 2.0
3,1, 3.0
4,1, 3.5
5,1, 4.0
6,1, 4.5
7,1, 5.0
8,1, 5.5
9,1, 6.0
10,1, 6.0
11,1, 8.0
12,1, 10.0
13,1, 15.0
14,1, 22.0
15,1, 30.0
**
** Define the rigid body
**
*NODE, NSET=QA_TEST_ALLNODES
1000,0.,-1.0
**
*BOUNDARY
1000,ENCASTRE
**
*NSET,NSET=CONNSET,GEN
1,15,1
*SURFACE,TYPE=NODE, NAME=MASSES
CONNSET,
*SURFACE,TYPE=CYLINDER,NAME=OBSTACLE,FILLET=0.001
0., 0., 0., 1., 0., 0.
0., 0., -1.
START,0.0,0.0
LINE,1.0,0.0
LINE, 1.5,-0.5
LINE, 2.0,-1.5
LINE, 2.0,-4.0
CIRCL,3.,-5.,3.,-4.
CIRCL,4.,-4.,3,-4.
LINE, 4.,-3.
LINE, 4.,-1.
CIRCL,5.,0.,5.,-1.
CIRCL,6.,-1.,5.,-1.
LINE, 6.,-3.
CIRCL,7.,-4.,7.,-3.
CIRCL,8.,-3.,7.,-3.
LINE, 8.,-1.
CIRCL,9.,0.,9.,-1.
CIRCL,10.,-1.,9.,-1
CIRCL,11.,-2.,11.,-1.
LINE, 12.,-0.5
LINE, 13.,-2.0
LINE, 14.,-3.5
LINE, 14.,-4.0
LINE, 14.,-4.5
LINE, 14.,-5.0
LINE, 14.5,-5.0
LINE, 14.5,-4.5
LINE, 14.5,-4.0
LINE, 14.5,-3.5
CIRCL,15.5,-2.5,14.5,-2.5
CIRCL,18.,0.,18.,-2.5
CIRCL,20.5,2,18.,2.5
CIRCL,18.,5.,18.,2.5
PARAB,15.,4.5,12.,3.5
*RIGID BODY, REF NODE=1000, ANALYTICAL SURFACE =OBSTACLE
*NSET, NSET=NNODE
NHEAVY, NLIGHT
**
*NSET, NSET=QA_TEST_ALLNODES
QA_TEST_ALLNODES, NHEAVY, NLIGHT, CONNSET
NNODE, 
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT
0.5,25.
**
** Contact of point masses and analytical rigid surface
**
*CONTACT PAIR
MASSES, OBSTACLE
**
*CLOAD
NHEAVY, 2,- 1.4
NLIGHT, 2, -.1
**
*RESTART,WRITE,NUM=50
*FILE OUTPUT, NUM=4,TIMEMARKS=YES
*NODE FILE
 U, 
***OUTPUT,HISTORY,TIME INTERVAL=0.05
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLWK,ETOTAL,DT
**
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
 U, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=6.25
*ENDSTEP