** connector friction tests for JOIN+ROTATION, JOIN+CARDAN, JOIN+Flex
**   -CUSTOMIZED + PRESTRESS
** friction results should be identical between the three connector elements
**------------------------------------------------
*parameter
rota = 0.8
Rs   = 0.20
*preprint,model=yes
*Node,nset = all
101,0,0,0
102,0,0,0
201,1,0,0
202,1,0,0
301,2,0,0
302,2,0,0
401,3,0,0
402,3,0,0
*orientation,name=global
1,0,0, 0,1,0
3,0.0
*orientation,name=global2
1,0,0,  0,1,0
3,0.0
*orientation,name=global3
1,0,0,  0,1.0,0.0000001
3,0.0
**-------------------------------------------------
*element,type=mass,elset=mass
10001,101
10002,102
20001,201
20002,202
30001,301
30002,302
40001,401
40002,402
*mass,elset=mass
1.0,
*element,type=rotaryI,elset=rotaryI
110001,101
110002,102
120001,201
120002,202
130001,301
130002,302
140001,401
140002,402
*rotary inertia,elset=rotaryI
1.0,1.0,1.0
**-------------------------------------------------
*element,type=conn3d2,elset=ballCustomRota
1001,101,102
*connector section,elset=ballCustomRota,behavior=Rota1
join,rotation
global,global2
**--------------
*element,type=conn3d2,elset=ballCustomCard
2001,201,202
*connector section,elset=ballCustomCard,behavior=Rota2
join,cardan
global,global2
**--------------
*element,type=conn3d2,elset=ballCustomFlex
3001,301,302
*connector section,elset=ballCustomFlex,behavior=Rota3
join,flexion
global,global2
**--------------
*element,type=conn3d2,elset=ballCustomEuler
4001,401,402
*connector section,elset=ballCustomEuler,behavior=Rota4
join,euler
global,global3
**--------------
*connector behavior,name=Rota1
*connector friction
1.25
*connector potential
4,<Rs>
5,<Rs>
6,<Rs>
*friction
0.1,
*connector behavior,name=Rota2		    
*connector friction
1.25
*connector potential
4,<Rs>
5,<Rs>
6,<Rs>
*friction
0.1,
*connector behavior,name=Rota3	
*connector constitutive reference
0.,0.,0.,0.,0.,-45.	    
*connector friction
1.25
*connector potential
4,<Rs>
5,<Rs>
6,<Rs>
*friction
0.1,
*connector behavior,name=Rota4	
*connector constitutive reference
0.,0.,0.,45.,0.,-45.	    
*connector friction
1.25
*connector potential
4,<Rs>
5,<Rs>
6,<Rs>
*friction
0.1,
**--------------------------------------------------
**--------------------------------------------------
*boundary
101,1,6
201,1,6
301,1,6
401,1,6
*elset,elset=conn
ballCustomRota, ballCustomCard, ballCustomFlex,ballCustomEuler
*initial conditions, type=vel
102,4,<rota>
102,5,<rota>
102,6,<rota>
202,4,<rota>
202,5,<rota>
202,6,<rota>
302,4,<rota>
302,5,<rota>
302,6,<rota>
402,4,<rota>
402,5,<rota>
402,6,<rota>
**--------------------------------------------------
*step
*dynamic, explicit, direct user control
1.0E-04, 1.0
*boundary, type=vel
102,4,6,<rota>
202,4,6,<rota>
302,4,6,<rota>
402,4,6,<rota>
*OUTPUT, FIELD, NUMBER INTERVAL=5
*NODE OUTPUT,nset=all
U,
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT,nset=all
U,RF
*ELEMENT OUTPUT,elset=conn
CUR1,CUR2,CUR3,
CASUR1,CASUR2,CASUR3,CASUC,
CSM1,CSM2,CSM3,CSFC,
CNFC,
*ENERGY OUTPUT,VARIABLE=ALL
*end step
**--------------------------------------------------