** connector friction tests for JOIN+ROTATION, JOIN+CARDAN, JOIN+Flex
**   -CUSTOMIZED + PRESTRESS
** friction results should be identical between the three connector elements
**------------------------------------------------
*preprint,model=yes
*parameter
Rs = 0.20
*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
*orientation,name=global
1,0,0, 0,1,0
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=fricCustomRota
join,rotation
global,
**--------------
*connector behavior,name=fricCustomRota
*connector derived component,name=norm123
1,2,3
1.0,1.0,1.0
*connector friction,contact force=norm123
1.25,0.0
1.20,1.0
*Connector potential
4,<Rs>
5,<Rs>
6,<Rs>
*friction
0.15,
**------
*connector elastic,comp=4
1.0,
*connector elastic,comp=5
1.0,
*connector elastic,comp=6
1.0,
**-------------------------------------------------
*element,type=conn3d2,elset=ballCustomCard
2001,201,202
*connector section,elset=ballCustomCard,behavior=fricCustomCard
join,cardan
global,
**--------------
*connector behavior,name=fricCustomCard
*connector derived component,name=norm123
1,2,3
1.0,1.0,1.0
*connector friction, contact force=norm123
1.25,    0.0
1.20,    1.0
*Connector potential
4,<Rs>
5,<Rs>
6,<Rs>
*friction
0.15,
**------
*connector elastic,comp=4
1.0,
*connector elastic,comp=5
1.0,
*connector elastic,comp=6
1.0,
**-------------------------------------------------
*element,type=conn3d2,elset=ballCustomFlex
3001,301,302
*connector section,elset=ballCustomFlex,behavior=fricCustomFlex
join,flexion
global,
**--------------
*connector behavior,name=fricCustomFlex
*connector derived component,name=norm123
1,2,3
1.0,1.0,1.0
*connector friction, contact force=norm123
1.25,   0.0
1.20,   1.0
*Connector potential
4,<Rs>
5,<Rs>
6,<Rs>
*friction
0.15,
**------
*connector elastic,comp=4
1.0,
*connector elastic,comp=5
1.0,
*connector elastic,comp=6
1.0,
**--------------------------------------------------
*boundary
101,1,6
201,1,6
301,1,6
*elset,elset=conn
 ballCustomRota, ballCustomCard, ballCustomFlex
*amplitude,name=backforth
0,0,  0.6, 0.6, 1.0,0.2
*amplitude,name=step
0,1,1,1
*initial conditions, type=vel
102,4,1.0
102,5,1.0
102,6,1.0
202,4,1.0
202,5,1.0
202,6,1.0
302,4,1.0
302,5,1.0
302,6,1.0
**--------------------------------------------------
*step
*dynamic, explicit, direct user control
1.0E-04,1.0
*boundary,type=vel,amplitude=backforth
102,4,6,1.0
202,4,6,1.0
302,4,6,1.0
*cload,amplitude=step
**load that creates constraint forces
102,1,10.0
102,2,20.0
102,3,30.0
**-----------
202,1,10.0
202,2,20.0
202,3,30.0
**-----------
302,1,10.0
302,2,20.0
302,3,30.0
**-----------
*OUTPUT, FIELD, NUMBER INTERVAL=5
*NODE OUTPUT,nset=all
U,
*OUTPUT,HISTORY
*NODE OUTPUT,nset=all
U,RF
*ELEMENT OUTPUT,elset=conn
CU,CASU,CEF,CSF,CTF,CNF,CSF
*ENERGY OUTPUT,VARIABLE=ALL
*end step
**--------------------------------------------------
*step
*dynamic, explicit, direct user control
1.0E-04, 1.0
*boundary,op=new
101,1,6
201,1,6
301,1,6
*cload,amplitude=step
**load that creates constraint forces
102,1,10.0
102,2,20.0
102,3,30.0
**-----------
202,1,10.0
202,2,20.0
202,3,30.0
**-----------
302,1,10.0
302,2,20.0
302,3,30.0
**load that creates motion
*cload
102,4,10.0
102,5,20.0
102,6,30.0
**-----------
202,4,10.0
202,5,20.0
202,6,30.0
**-----------
302,4,10.0
302,5,20.0
302,6,30.0
**--------------------------------------------------
*OUTPUT, FIELD, NUMBER INTERVAL=5
*NODE OUTPUT,nset=all
U,
*OUTPUT,HISTORY
*NODE OUTPUT,nset=all
U,RF
*ELEMENT OUTPUT,elset=conn
CU,CASU,CEF,CSF,CTF,CNF,CSF
*ENERGY OUTPUT,VARIABLE=ALL
*end step