*HEADING
 Tests surface based cohesive behavior with damage and repeated contacts
*preprint,model=yes
*NODE,NSET=SPHERE0,INPUT=tennis_ef1.inp
*NCOPY,CHANGE NUMBER=1000,OLDSET=SPHERE0, NEW SET=SPHERE1,SHIFT
5.,0.0,0.0

*NCOPY,CHANGE NUMBER=2000,OLDSET=SPHERE0, NEW SET=SPHERE2,SHIFT
10.,0.0,0.0

*NSET,NSET=SPHERE
SPHERE0,SPHERE1,SPHERE2, 
**
*ELEMENT,TYPE=S4R,ELSET=SPHERE0,
INPUT=tennis_ef2.inp
*ELCOPY,ELEMENT SHIFT=1000,OLD SET=SPHERE0,NEW SET=SPHERE1,SHIFT NODES=1000
*ELCOPY,ELEMENT SHIFT=2000,OLD SET=SPHERE0,NEW SET=SPHERE2,SHIFT NODES=2000
*ELSET, ELSET=SPHERE
SPHERE0,SPHERE1,SPHERE2,
*SHELL SECTION,ELSET=SPHERE,MATERIAL=STEEL
.2,3
**
*NODE, NSET=BASE1
1000000, 0.,-5.,0.
1000001,-10.,-5,-10.
1000002,20.,-5.,-10.
1000003,20.,-5.,7.
1000004,-10.,-5.,7.
*ELEMENT, TYPE=R3D4, ELSET=BASE1
1000001,1000001,1000002,1000003,1000004
*RIGID BODY, ELSET=BASE1, REFNODE=1000000
*NODE, NSET=BASE2
2000000, 0.,5.,0.
2000001,-10.,5,-10.
2000002,20.,5.,-10.
2000003,20.,5.,7.
2000004,-10.,5.,7.
*ELEMENT, TYPE=R3D4, ELSET=BASE2
2000001,2000001,2000002,2000003,2000004
*RIGID BODY, ELSET=BASE2, REFNODE=2000000
*NODE, NSET=BASE3
3000000, 0.,-5.,0.
3000001,-10.,-5,-10.
3000002,-10.,-5.,7.
3000003,-10.,5.,7.
3000004,-10.,5,-10.
*ELEMENT, TYPE=R3D4, ELSET=BASE3
3000001,3000001,3000002,3000003,3000004
*RIGID BODY, ELSET=BASE3, REFNODE=3000000
*NODE, NSET=BASE4
4000000, 0.,-5.,0.
4000001,20.,-5,-10.
4000002,20.,-5.,7.
4000003,20.,5.,7.
4000004,20.,5,-10.
*ELEMENT, TYPE=R3D4, ELSET=BASE4
4000001,4000001,4000002,4000003,4000004
*RIGID BODY, ELSET=BASE4, REFNODE=4000000
*NODE, NSET=BASE5
5000000, 0.,-5.,0.
5000001,-10.,-5,7.
5000002,20.,-5.,7.
5000003,20.,5.,7.
5000004,-10.,5,7.
*ELEMENT, TYPE=R3D4, ELSET=BASE5
5000001,5000001,5000002,5000003,5000004
*RIGID BODY, ELSET=BASE5, REFNODE=5000000
********************
**
** Material definitions.
**
********************
*MATERIAL,NAME=STEEL
*ELASTIC
30.E+6,.3
*DENSITY
.73E-3, 
********************
**
** Initial conditions, boundary conditions.
**
********************
*INITIAL CONDITIONS,TYPE=STRESS
**  Tension in strings is 10 lb.
SPHERE,22.5,22.5
*INITIAL CONDITIONS,TYPE=VELOCITY
**  Initial velocity of ball is 22 fps at a
**  15 degree angle of attack.
SPHERE0, 1, 0.
SPHERE0, 2, -600.0
SPHERE0, 3,0.
SPHERE1, 1, 0.
SPHERE1, 2, -600.0
SPHERE1, 3,0.
SPHERE2, 1, 0.
SPHERE2, 2, -600.0
SPHERE2, 3,0.
**
********************
*BOUNDARY
1000000,1,6
2000000,1,6
3000000,1,6
4000000,1,6
5000000,1,6
SPHERE0,1,1
SPHERE0,3,3
SPHERE1,1,1
SPHERE1,3,3
SPHERE2,1,1
SPHERE2,3,3
**
*SURFACE,TYPE=ELEMENT,NAME=SPHERE0
SPHERE0,SPOS
*SURFACE,TYPE=ELEMENT,NAME=SPHERE1
SPHERE1,SPOS
*SURFACE,TYPE=ELEMENT,NAME=SPHERE2
SPHERE2,SPOS
*SURFACE, TYPE=ELEMENT,NAME=BASE1
BASE1,SNEG
*SURFACE, TYPE=ELEMENT,NAME=BASE2
BASE2,SPOS
*SURFACE, TYPE=ELEMENT,NAME=BASE3
BASE3,SNEG
*SURFACE, TYPE=ELEMENT,NAME=BASE4
BASE4,SPOS
*SURFACE INTERACTION,NAME=INTER1
*SURFACE INTERACTION,NAME=INTER2
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
*SURFACE INTERACTION,NAME=INTER3
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
*DAMAGE INITIATION, CRITERION=MAXS
15.E+3,15.E+3,15.E+3
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR
.06,
*SURFACE INTERACTION,NAME=INTER4
*COHESIVE BEHAVIOR, REPEATED CONTACTS, Eligibility = Current Contacts 
*DAMAGE INITIATION, CRITERION=MAXS
15.E+3,15.E+3,15.E+3
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR
.06,
**
*CONTACT PAIR, INTERACTION=INTER4
SPHERE0,BASE1
*CONTACT PAIR, INTERACTION=INTER1
SPHERE0,BASE2
*CONTACT PAIR, INTERACTION=INTER1
SPHERE1,BASE1
*CONTACT PAIR, INTERACTION=INTER2
SPHERE1,BASE2
*CONTACT PAIR, INTERACTION=INTER3
SPHERE2,BASE1
*CONTACT PAIR, INTERACTION=INTER1
SPHERE2,BASE2
**
*nset,nset=three
30053,31053,32053
**
** STEP 1
**
*STEP, inc=10000
*DYNAMIC,HAFTOL=1.0E9
1.E-3,5.92E-3,,0.001
*OUTPUT, FIELD,FREQ=100
*NODE OUTPUT
U,
***CONTACT OUTPUT
***CFORCE,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=three
v2
*END STEP
**
** STEP 2
**
*STEP, inc=10000
*DYNAMIC,HAFTOL=1.0E8
1.E-8,4.08E-3,,0.001
*END STEP
**
** STEP 3
**
*STEP, inc=10000
*DYNAMIC,HAFTOL=1.0E9
1.E-03,7.0E-3,,0.001
*END STEP
**
** STEP 4
**
*STEP, inc=10000
*DYNAMIC,HAFTOL=1.0E9
1.E-04,15.0E-3,,0.0001
*END STEP
**
** STEP 5
**
*STEP, inc=10000
*DYNAMIC,HAFTOL=1.0E10
1.E-2,22.47E-3,,0.01
*END STEP
**
** STEP 6
**
*STEP, inc=10000
*DYNAMIC,HAFTOL=1.0E9
1.E-7,5.0E-3,,0.0001
*END STEP
**
** STEP 7
**
*STEP, inc=10000
*DYNAMIC,HAFTOL=1.0E9
1.E-5,65.2E-3,,0.001
**1.E-5,100.08E-3,,0.001
**1.E-7,125.E-3,,0.00001
**,125.E-3
*OUTPUT, FIELD,FREQ=100
*NODE OUTPUT
U,
*CONTACT OUTPUT
CSDMG
*END STEP