*HEADING
GLOBAL TRACKING FOR SEGMENT ANALYTICAL RIGID SURFACE
*PREPRINT,HISTORY=YES,MODEL=YES
**
** Define a dummy truss element
**
*NODE, NSET=QA_TEST_ALLNODES
101, 0.,-6.
102, 20.,-6.
**
*ELEMENT, TYPE=T2D2, 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
2,0.5,1.0
3,0.5,1.0
4,0.5,1.0
5,0.5,1.0
6,0.5,1.0
7,0.5,1.0
8,0.5,1.0
9,0.5,1.0
10,0.5,1.5
11,0.5,1.5
12,0.5,1.5
*NODE, NSET=NLIGHT
13,0.5,1.5
14,0.5,1.5
15,0.5,1.5
**
*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=SEGMENTS,NAME=OBSTACLE,FILLET=0.001
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
PARAB,14.,-3.6,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
**
**The following parabola really defines line segment. It
**will be converted to a LINE segment by the code.
**
*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