** connector friction tests for UJOINT
**   -PREDEFINED
**   -CUSTOMIZED
** results should be identical between the two connector elements
** using a zero angle between axes (elements 1001 and 2001)
** results should not be compared with element 3001 (angle between axes)
**
** the aspect ratio La/Rp = 2.0 was used to interpolate the values for 
** aaxi, baxi, atwi, and btwi
**------------------------------------------------
*preprint,model=yes
*parameter 
Rp   = 0.12
Ra   = 0.13
Ls   = 0.10
La   = 0.24
aaxi = 0.3604
baxi = 0.1396
atwi = 0.0203
btwi = 0.2718
m1   = 2.0*Ra*aaxi
m2   = 0.5*Rp
m3   = 2.0*btwi*Rp/La
m4   = 2.0*baxi*Rp
m5   = Rp*La/(2.0*Ls)
m6   = 4.0*atwi*Rp/Ls
**
angleaxes1 =   0.0
angleaxes2 = -30.0
twopi=2.0*3.141592
**------------------------------------------------
*Node,nset = all
101,0,0,0
102,1,0,0
201,2,0,0
202,3,0,0
301,4,0,0
302,5,0,0
*orientation,name=oria1
0,1,0, -1,0,0
3,<angleaxes1>
*orientation,name=oria2
0,1,0, -1,0,0
3,<angleaxes2>
*orientation,name=orib
0,1,0, -1,0,0
3,0.0
**-------------------------------------------------
*element,type=mass,elset=mass
10001,101
10002,102
20001,201
20002,202
30001,301
30002,302
*mass,elset=mass
1.0,
*element,type=rotaryI,elset=rotaryI
110001,101
110002,102
120001,201
120002,202
130001,301
130002,302
*rotary inertia,elset=rotaryI
1.0,1.0,1.0
**-------------------------------------------------
*element,type=conn3d2,elset=ujointPredef1
1001,101,102
*connector section,elset=ujointPredef1,behavior=fricPredef1
ujoint,
oria1,orib
**-------
*connector behavior,name=fricPredef1
*connector friction, predefined
<Rp>,<La>,<Ra>,<Ls>,

*friction
0.1,
**--------------------------------------------------
*element,type=conn3d2,elset=ujointCustom1
2001,201,202
*connector section,elset=ujointCustom1,behavior=fricCustom1
ujoint,
oria1,orib
**-------
*connector behavior,name=fricCustom1
*connector derived component,name=norm4
1,
<m1>
*connector derived component,name=norm4
2,5,3
<m2>,<m3>,<m4>
*connector derived component,name=norm4
2,5
<m5>,<m6>
*connector friction,component=4,contact force=norm4
*friction
0.1,
**-----
*connector derived component,name=norm6
3,
<m1>
*connector derived component,name=norm6
2,5,1
<m2>,<m3>,<m4>
*connector derived component,name=norm6
2,5
<m5>,<m6>
*connector friction,component=6,contact force=norm6
*friction
0.1,
**-------------------------------------------------
*element,type=conn3d2,elset=ujointPredef2
3001,301,302
*connector section,elset=ujointPredef2,behavior=fricPredef2
ujoint,
oria2,orib
**-------
*connector behavior,name=fricPredef2
*connector friction,predefined
<Rp>,<La>,<Ra>,<Ls>,

*friction
0.1,
**--------------------------------------------------
*boundary
101,1,6
201,1,6
301,1,6
*elset,elset=conn
 ujointPredef1,ujointCustom1,ujointPredef2
**--------------------------------------------------
*amplitude, name=ramp
0.0,0.0,1.0,1.0
*step
*dynamic, explicit, direct user control
1.0E-04, 1.0
*connector motion, type=vel
1001,4,1.0
1001,6,1.0
2001,4,1.0
2001,6,1.0
3001,4,1.0
3001,6,1.0
*cload, amp=ramp
**load that creates constraint forces
102,1,10.0
102,2,11.0
102,3,12.0
**-----------
202,1,10.0
202,2,11.0
202,3,12.0
**-----------
302,1,10.0
302,2,11.0
302,3,12.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
ALLFD,
ELEN,
*end step
**--------------------------------------------------