*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=TRANS_PLAS
CARTESIAN, CARDAN
SR-SS,
*CONNECTOR BEHAVIOR, NAME=TRANS_PLAS
*CONNECTOR ELASTICITY, RIGID, COMPONENT=2
*CONNECTOR ELASTICITY, RIGID, COMPONENT=3
*CONNECTOR ELASTICITY, RIGID, COMPONENT=4
*CONNECTOR ELASTICITY, RIGID, COMPONENT=5
*CONNECTOR ELASTICITY, RIGID, COMPONENT=6
*CONNECTOR DERIVED COMPONENT, NAME=T_TRANS
2,3
1,1
*CONNECTOR DERIVED COMPONENT, NAME=T_BEND
5,6
6.0,6.0
*CONNECTOR DERIVED COMPONENT, NAME=T_TWIST
4,
1.0,
*CONNECTOR PLASTICITY
*CONNECTOR POTENTIAL
T_TRANS,
T_BEND,
T_TWIST
*CONNECTOR HARDENING
1000000,0.0
1750000,0.5
2000000,1.0
2500000,5.0
3000000,10.0
1004000000,30.0
*****************************
** 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=HINGE1_PLAS
CARTESIAN,CARDAN
SS-RS,
*CONNECTOR BEHAVIOR, NAME=HINGE1_PLAS
*CONNECTOR ELASTICITY, RIGID, COMPONENT=1
*CONNECTOR ELASTICITY, RIGID, COMPONENT=2
*CONNECTOR ELASTICITY, RIGID, COMPONENT=3
*CONNECTOR ELASTICITY, RIGID, COMPONENT=5
*CONNECTOR ELASTICITY, RIGID, COMPONENT=6
*CONNECTOR PLASTICITY
*CONNECTOR POTENTIAL
H1_TRANS,
H1_BEND,
*CONNECTOR HARDENING
600000,0
630000,0.5
660000,1.0
690000,1.5
720000,2.5
750000,4.0
800000,10.0
900000,30.0
1400000,100.0
*CONNECTOR DERIVED COMPONENT, NAME=H1_TRANS
1,2,3
1,1,1
*CONNECTOR DERIVED COMPONENT, NAME=H1_BEND
5,6
2.5,2.5
*****************************
** 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
cylindrical,
RR-RS
*****************************
** 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, BEHAVIOR=JOIN_PLAS
cartesian, cardan
RR-LD
*CONNECTOR BEHAVIOR, NAME=JOIN_PLAS
*CONNECTOR ELASTICITY, RIGID, COMPONENT=1
*CONNECTOR ELASTICITY, RIGID, COMPONENT=2
*CONNECTOR ELASTICITY, RIGID, COMPONENT=3
*CONNECTOR PLASTICITY
*CONNECTOR POTENTIAL
1,
2,
3,
*CONNECTOR HARDENING
100000,0
110000,0.5
120000,1
130000,1.5
140000,2
150000,2.5
160000,3.0
170000,3.5
180000,4.0
190000,6.0
200000,15.0
300000,30.0
*****************************
** 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=HINGE2_PLAS
CARTESIAN,CARDAN
RR-SB
*CONNECTOR BEHAVIOR, NAME=HINGE2_PLAS
*CONNECTOR ELASTICITY, RIGID, COMPONENT=1
*CONNECTOR ELASTICITY, RIGID, COMPONENT=2
*CONNECTOR ELASTICITY, RIGID, COMPONENT=3
*CONNECTOR ELASTICITY, RIGID, COMPONENT=5
*CONNECTOR ELASTICITY, RIGID, COMPONENT=6
*CONNECTOR DERIVED COMPONENT, NAME=H2_TRANS
1,2,3
1,1,1
*CONNECTOR DERIVED COMPONENT, NAME=H2_BEND
5,6
4.0,4.0
*CONNECTOR PLASTICITY
*CONNECTOR POTENTIAL
H2_TRANS,
H2_BEND,
*CONNECTOR HARDENING
600000,0
630000,0.5
660000,1.0
690000,1.5
720000,2.5
750000,4.0
800000,10.0
900000,30.0
1400000,100.0
*****************************
** 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
slide-plane, revolute
SB-R,
*****************************
** 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
ujoint,
LD-DL,DL-LD
*****************************
** 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
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,.03,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
*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
U,
*ENERGY FILE
*OUTPUT,FIELD,NUM=300
*NODE OUTPUT
U,V,A
*ELSET, ELSET=PLAS
5004,5001,5002,5005
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=OUT
U,V,A
*ELEMENT OUTPUT, ELSET=PLAS
CUP, CUPEQC, CEF, CRF
*NSET, NSET=QA_TEST
OUT,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=4
*NODE OUTPUT, NSET=QA_TEST
U,
*END STEP