*HEADING
  PULL IN CABLE - ANISOTROPIC FRICTION       
*NODE
17,550.,200.,.861
1,0.,,.861
16,750.,,.861
1000,0.
*NGEN,NSET=ALLNODES
1,16
*ELEMENT,TYPE=B31H
1,1,2
*ELGEN,ELSET=BEAM
1,15,1,10
*ELEMENT,TYPE=GAPSPHER,ELSET=CABLE
1001,17,16
*BEAM GENERAL SECTION,SECTION=PIPE,ELSET=BEAM
.41666,.08333
,1.
4.32E9,2.16E9
*rigid body,analytical surface=BSURF,REF NODE=1000
*surface,TYPE=CYLINDER,NAME=BSURF
-1000.,-1000.,0.,0.,-1000.,0.
-1000.,1000.,0.
START,-10.
LINE,10000.
*surface,NAME=ASURF
BEAM, 
*CONTACT PAIR,INTERACTION=ROUGH
ASURF,BSURF
*SURFACE INTERACTION,NAME=ROUGH
*FRICTION,ANISOTROPIC,ELASTIC SLIP=0.8333
.6,1.
*SURFACE BEHAVIOR, PRESSURE-OVERCLOSURE=EXPONENTIAL
1.,5000.
*GAP,ELSET=CABLE
300., 
*AMPLITUDE,NAME=A1
0.,0.,1.,-297.2
*STEP,NLGEOM
  INSURE EQUILIBRIUM BEFORE PULL IN
*STATIC
1.,1.,1.,1.
*BOUNDARY
17,1,3
1000,ENCASTRE
ALLNODES,3,5
1,1,2
1,4,6
*EL PRINT,ELSET=BEAM,FREQUENCY=999
SF, 
*EL PRINT,ELSET=CABLE
S,E
*EL FILE,ELSET=CABLE
S,E              
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=CABLE
S,E
*RESTART,WRITE,FREQUENCY=1
*END STEP
*STEP,NLGEOM,INC=50
 PULL END IN
*STATIC
.02,1.,.002,.1
*CONTACT INTERFERENCE,TYPE=ELEMENT,AMPLITUDE=A1
CABLE,1.
*BOUNDARY,OP=NEW
ALLNODES,3,5
17,1,3
1000,ENCASTRE
*MONITOR,NODE=16,DOF=2
*CONTACT FILE,SLAVE=ASURF,FREQUENCY=999
*OUTPUT,FIELD,FREQ=999
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF
*EL FILE,FREQUENCY=999
S, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S,
*NODE OUTPUT, NSET=ALLNODES
U,
*RESTART,WRITE,FREQUENCY=20,OVERLAY
*END STEP