** connector friction tests for PLANAR
**   -PREDEFINED + PRESTRESS
**   -CUSTOMIZED + PRESTRESS
** results should be identical between the two connector elements
** PLANAR is the only connection type where we allow for both 
** *CF, comp=nr and *CF, comp=slip dir
**------------------------------------------------
*preprint,model=yes
*parameter 
R    = 0.12
f1   = 1.0
f2   = 1.0/R
m1   = (2.0/3.0)*R
m2   = 1.0
**------------------------------------------------
*Node,nset = all
101,0,0,0
102,1,0,0
201,2,0,0
202,3,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
*mass,elset=mass
1.0,
*element,type=rotaryI,elset=rotaryI
110001,101
110002,102
120001,201
120002,202
*rotary inertia,elset=rotaryI
1.0,1.0,1.0
**-------------------------------------------------
*element,type=conn3d2,elset=planPredef
1001,101,102
*connector section,elset=planPredef,behavior=fricPredef
Slide-plane,revolute
global,
**-------
*connector behavior,name=fricPredef
*connector friction, predefined
<R>,1.25
*friction
0.10,
**-------
*connector elastic,comp=2
12.0,
*connector elastic,comp=3
12.0,
*connector elastic,comp=4
1.0,
**--------------------------------------------------
*element,type=conn3d2,elset=planCustom
2001,201,202
*connector section,elset=planCustom,behavior=fricCustom
planar,
global,
**--------------
*connector behavior,name=fricCustom
*connector derived component,name=norm23
1,
<f1>,
*connector derived component,name=norm23
5,6
<f2>,<f2>
*connector friction,contact force=norm23
1.25,   0.0
1.20,1000.0
*connector potential
2,
3,
*friction
0.10,
**-----
*connector derived component,name=norm4
1,
<m1>,
*connector derived component,name=norm4
5,6
<m2>,<m2>
*connector friction,contact force=norm4, component=4
*friction
0.10,
**------
*connector elastic,comp=2
12.0,
*connector elastic,comp=3
12.0,
*connector elastic,comp=4
1.0,
**--------------------------------------------------
*boundary
101,1,6
201,1,6
*elset,elset=conn
 planPredef,planCustom
**--------------------------------------------------
*initial conditions, type=vel
102,2,2.0
102,3,3.0
102,4,1.0
202,2,2.0
202,3,3.0
202,4,1.0
*step
*dynamic, explicit, direct user control
1.0E-04,1.0
*cload
**load that creates constraint forces
102,1,11.0
102,5,1.0
202,1,11.0
202,5,1.0
*boundary, type=vel
102,2,2,2.0
102,3,3,3.0
102,4,4,1.0
202,2,2,2.0
202,3,3,3.0
202,4,4,1.0
*OUTPUT, FIELD, NUMBER INTERVAL=5
*NODE OUTPUT,nset=all
U,
*OUTPUT,HISTORY
*NODE OUTPUT,nset=all
U,
*ELEMENT OUTPUT,elset=conn
CU,CASU,CEF,CSF,CTF,CNF,CSF,CIVC
ELEN,
*end step
**--------------------------------------------------
*step
*dynamic, explicit, direct user control
1.0E-04,1.0
*cload
**load that creates constraint forces
102,1,11.0
102,5,1.0
202,1,11.0
202,5,1.0
** load that creates motion
102,2,10.0
102,3,10.0
102,6,10.0
202,2,10.0
202,3,10.0
202,6,10.0
*boundary,op=new
101,1,6
201,1,6
*OUTPUT, FIELD, NUMBER INTERVAL=5
*NODE OUTPUT,nset=all
U,
*OUTPUT,HISTORY
*NODE OUTPUT,nset=all
U,
*ELEMENT OUTPUT,elset=conn
CU,CASU,CEF,CSF,CTF,CNF,CSF,CIVC
ELEN,
*end step
**------------------------------------------------