*HEADING
  PULL IN CABLE - ANISOTROPIC FRICTION
*NODE
*NGEN,NSET=ALLNODES
*ELEMENT,TYPE=B31H
*ELGEN,ELSET=BEAM
*ELEMENT,TYPE=GAPSPHER,ELSET=CABLE
*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
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,
*RESTART,WRITE,FREQUENCY=20,OVERLAY
*END STEP