** connector friction tests for HINGE
**   -PREDEFINED + PRESTRESS
**   -CUSTOMIZED + PRESTRESS
** results should be identical between the two connector elements
** this particular example is depicted in the user's manual
**------------------------------------------------
*preprint,model=yes
*parameter 
Rp   = 0.12
Ra   = 0.14
Ls   = 0.65
m1   = Ra
m2   = Rp
m3   = 2.0*Rp/Ls
**------------------------------------------------
*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=hinPredef
1001,101,102
*connector section,elset=hinPredef,behavior=fricPredef
hinge,
global,global
**-------
*connector behavior,name=fricPredef
*connector elastic,comp=4
12.0,
*connector friction,predefined,stick stiffness=1000.0
<Rp>,<Ra>,<Ls>,10.25
*friction
0.15,
**--------------------------------------------------
*element,type=conn3d2,elset=hinCustom
2001,201,202
*connector section,elset=hinCustom,behavior=fricCustom
hinge,
global,global
**-------
*connector behavior,name=fricCustom
*connector elastic,comp=4
12.0,
*connector derived component,name=norm4
1,
<m1>,
*connector derived component,name=norm4
2,3
<m2>,<m2>
*connector derived component,name=norm4
5,6
<m3>,<m3>
*connector friction,component=4,contact force=norm4,stick stiffness = 1000.0
10.25,    0.0
10.20, 1000.0
*friction
0.15,
**--------------------------------------------------
*boundary
101,1,6
201,1,6
*elset,elset=conn
 hinPredef,hinCustom
*initial condition, type=vel
102,4,1.0
202,4,1.0
**--------------------------------------------------
*step
*dynamic, explicit, direct user control
1.0E-04, 1.0
*cload
**load that creates constraint forces
102,2,11.0
102,5,1.0
202,2,11.0
202,5,1.0
*boundary, type=vel
102,4,4,1.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
*ENERGY OUTPUT,VARIABLE=ALL
*end step
**--------------------------------------------------
*step
*dynamic, explicit, direct user control
1.0E-04, 1.0
*cload
**load that creates constraint forces
102,2,11.0
102,5,1.0
202,2,11.0
202,5,1.0
**load that creates motion
102,4,100. 0
202,4,100.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
*ENERGY OUTPUT,VARIABLE=ALL
*end step
**-------------------------------------------------