*HEADING
 Rigid, multi-body mechanism - connector elements
** node and element definitions:
*INCLUDE, INPUT= rigmultimech_bulk.inp
*****************************
*ELSET, ELSET=SQUAREROD, GEN 
1,4,1
*ELSET, ELSET=SQUARESLEEVE, GEN 
5,92,1
*ELSET, ELSET=ROUNDROD, GEN 
93,172,1
*ELSET, ELSET=ROUNDSLEEVE, GEN 
173,260,1
*ELSET, ELSET=SLIDINGBLOCK, GEN 
261,272,1
*ELSET, ELSET=RACEWAY, GEN 
273,278,1
*ELSET, ELSET=LARGEDISK, GEN 
1001,1174,1
*ELSET, ELSET=DISKLINK, GEN
1175, 1194,1
*ELSET, ELSET=SMALLDISK, GEN
1195, 1368,1
*****************************
** Square rod
*NSET, NSET=ADD1
19
*RIGID BODY,ELSET=SQUAREROD,REF NODE=10000,POSITION=CENTER OF MASS,TIE NSET=ADD1
** Square sleeve
*NSET, NSET=ADD2
20,218
*RIGID BODY,ELSET=SQUARESLEEVE,REF NODE=20000,POSITION=CENTER OF MASS,TIE NSET=ADD2
** Round rod
*NSET, NSET=ADD3
385,746,1034
*RIGID BODY,ELSET=ROUNDROD,REF NODE=30000,POSITION=CENTER OF MASS,TIE NSET=ADD3
** Round sleeve
*NSET, NSET=ADD4
219,386
*RIGID BODY,ELSET=ROUNDSLEEVE,REF NODE=40000,POSITION=CENTER OF MASS,TIE NSET=ADD4
** Large disk
*NSET, NSET=ADD5
747,768
*RIGID BODY,ELSET=LARGEDISK,REF NODE=70000,POSITION=CENTER OF MASS,TIE NSET=ADD5
** Sliding block
*NSET, NSET=ADD6
1035, 1040
*RIGID BODY,ELSET=SLIDINGBLOCK,REF NODE=50000,POSITION=CENTER OF MASS,TIE NSET=ADD6
** Raceway
*NSET,NSET=ADD7
1041,
*RIGID BODY,ELSET=RACEWAY,REF NODE=60000,POSITION=CENTER OF MASS,TIE NSET=ADD7
** DiskLink
*NSET,NSET=ADD8
795,798
*RIGID BODY,ELSET=DISKLINK,REF NODE=80000,POSITION=CENTER OF MASS,TIE NSET=ADD8
** Smalldisk
*NSET,NSET=ADD9
821,
*RIGID BODY,ELSET=SMALLDISK,REF NODE=90000,POSITION=CENTER OF MASS,TIE NSET=ADD9
**
*NSET, NSET=REF
10000,20000,30000,40000,50000,60000,70000,80000,90000
*****************************
** Square rod - square sleeve
*****************************
*ELEMENT, TYPE = CONN3D2, ELSET=SR-SS
5001, 19, 20
*ORIENTATION, DEFINITION=COORDINATES, NAME=SR-SS
1, 0, 0, 0, 1, 0
*CONNECTOR SECTION, ELSET=SR-SS,behavior=tran_fric
translator,
SR-SS,
*connector behavior, name=tran_fric
*connector friction,predefined
1.0,6.0
*friction
0.12,
*****************************
** Square sleeve - round sleeve
*****************************
*ELEMENT, TYPE = CONN3D2, ELSET=SS-RS
5002, 218, 219
*ORIENTATION, DEFINITION=COORDINATES, NAME=SS-RS
0, 0, 1, 1, 0, 0
*CONNECTOR SECTION, ELSET=SS-RS,behavior = hinge_fric
hinge,
SS-RS,
*connector behavior, name=hinge_fric
*connector friction, predefined
1.0,1.2,4.0
*friction
0.12,
*****************************
** Round rod - round sleeve
*****************************
*ELEMENT, TYPE = CONN3D2, ELSET=RR-RS
5003, 385, 386
*ORIENTATION, DEFINITION=COORDINATES, NAME=RR-RS
0, 1, 0, 0, 0, 1
*CONNECTOR SECTION, ELSET=RR-RS,behavior=cyl_fric
cylindrical,
RR-RS
*connector behavior, name=cyl_fric
*connector friction, predefined
1.0,5.0
*friction
0.12,
*****************************
** Round rod - large disk
*****************************
*ELEMENT, TYPE = CONN3D2, ELSET=RR-LD
5004, 746, 747
*ORIENTATION, DEFINITION=COORDINATES, NAME=RR-LD
1, 0, 0, 0, 1, 0
*CONNECTOR SECTION, ELSET=RR-LD
join,
RR-LD
*****************************
** Round rod - sliding block
*****************************
*ELEMENT, TYPE = CONN3D2, ELSET=RR-SB
5005, 1034, 1035
*ORIENTATION, DEFINITION=COORDINATES, NAME=RR-SB
0, 0, 1, 1, 0, 0
*CONNECTOR SECTION, ELSET=RR-SB,behavior = hinge_fric
hinge,
RR-SB
*****************************
** Sliding block - raceway
*****************************
*ELEMENT, TYPE = CONN3D2, ELSET=SB-R
5006, 1040, 1041
*ORIENTATION, DEFINITION=COORDINATES, NAME=SB-R
1, 0, 0, 0, 1, 0
*CONNECTOR SECTION, ELSET=SB-R,behavior = planar_fric
slide-plane, revolute
SB-R,
*connector behavior, name=planar_fric
*connector friction, predefined
3.0,
*friction
0.12,
*****************************
** Large disk - disklink
*****************************
*ELEMENT, TYPE = CONN3D2, ELSET=LD-DL
5007, 768, 795
*ORIENTATION, DEFINITION=COORDINATES, NAME=LD-DL
-1, 0, 0, 0, 0, 1
*ORIENTATION, DEFINITION=COORDINATES, NAME=DL-LD
-.8660254,0.,.5,.5,0.,.8660254
*CONNECTOR SECTION, ELSET=LD-DL,behavior=ujoint_fric
ujoint,
LD-DL,DL-LD
*connector behavior, name=ujoint_fric
*connector frIction, predefined
1.0,1.25,1.26,2.0

*friction
0.12,
*****************************
** Small Disk - Disk Link
*****************************
*ELEMENT, TYPE = CONN3D2, ELSET=SD-DL
5008, 821, 798
*ORIENTATION, DEFINITION=COORDINATES, NAME=SD-DL
-1, 0, 0, 0, 0, 1
*ORIENTATION, DEFINITION=COORDINATES, NAME=DL-SD
-.8660254,0.,.5,.5,0.,.8660254
*CONNECTOR SECTION, ELSET=SD-DL,behavior=ujoint_fric
ujoint,
SD-DL,DL-SD
**************************
*BOUNDARY 
10000,1,2,0.0
10000,4,6,0.0
70000,1,5,0.0
90000,1,2,0.0
90000,4,5,0.0
60000,1,6,0.0
*AMPLITUDE, NAME=AMP, DEFINITION=SMOOTH STEP
0.0,0.0,.00375,23.26,3.0,23.26
***********************************
*NSET, NSET=CONNECTIONNODES
746,747,386,385,20,19,
219,218,1034,1035,1040,1041,
768,795,798,821
*NSET,NSET=OUT
CONNECTIONNODES,REF
*NSET,NSET=QA_TEST
OUT,
*elset,elset=conn_fric
SR-SS,SS-RS,RR-RS,RR-SB,SB-R,LD-DL,SD-DL
*TIME POINTS, NAME=TIME-PTS-1, GEN
0.,0.18,0.018
*TIME POINTS, NAME=TIME-PTS-2, GEN
0.,0.18,0.0006
*STEP,NLGEOM=YES
*DYNAMIC,EXPLICIT, DIRECT
1.8E-4,0.18
*BOUNDARY, TYPE=VELOCITY, AMPLITUDE=AMP
90000,6,6,5.0
***************************************
** *FILE OUTPUT, TIMEMARKS=YES,NUM=4
** *NODE FILE,NSET=OUT,FREQ=10
** U,V,A
** *ENERGY FILE,FREQ=5
*OUTPUT,FIELD,TIME POINTS=TIME-PTS-2
*NODE OUTPUT
U,V,A
*********** QA **********
*OUTPUT,FIELD,TIME POINTS=TIME-PTS-1
*NODE OUTPUT,NSET=QA_TEST
U,
*OUTPUT,HISTORY,VAR=PRESELECT,TIMEINT=0.009
***********************
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=OUT
U,V,A
*element output,elset=conn_fric
CU,CASU,CNF,CSF
*END STEP