*HEADING
 Tests surface based cohesive behavior with damage and repeated contacts
*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, -700.0
SPHERE0, 3,0.
SPHERE1, 1, 0.
SPHERE1, 2, -700.0
SPHERE1, 3,0.
SPHERE2, 1, 0.
SPHERE2, 2, -700.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,
*SURFACE, TYPE=ELEMENT,NAME=BASE2
BASE2,
*SURFACE, TYPE=ELEMENT,NAME=BASE3
BASE3,
*SURFACE, TYPE=ELEMENT,NAME=BASE4
BASE4,
*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,1.E-6,1.E-6
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR
.02,
*SURFACE INTERACTION,NAME=INTER4
*COHESIVE BEHAVIOR, REPEATED CONTACTS, Eligibility = Current Contacts 
*DAMAGE INITIATION, CRITERION=MAXS
15.E+3,1.E-6,1.E-6
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR
.02,
*STEP
*DYNAMIC,EXPLICIT
,125.E-3
*CONTACT
*CONTACT INCLUSIONS,ALL EXTERIOR
*CONTACT PROPERTY ASSIGNMENT
SPHERE0,BASE1,INTER4
SPHERE1,BASE2,INTER2
SPHERE2,BASE1,INTER3
*OUTPUT, FIELD, NUMBER INTERVAL=4
*NODE OUTPUT
U,
***CONTACT OUTPUT
***CFORCE,CSTRESS
*END STEP