*HEADING
 : *MODEL CHANGE verification for 3D finite-sliding contact
*PREPRINT,MODEL=YES
*NODE
1,0.0,0.0
7,7.5,0.0
11,,2.
17,7.5,2.0
101,0.5,2.1,0.5
103,2.25,2.1,0.5
111,0.5,4.1,0.5
113,2.25,4.1,0.5
21,0.0,0.0
27,7.5,0.0
31,,2.
37,7.5,2.0
121,0.5,2.1,0.5
123,2.25,2.1,0.5
131,0.5,4.1,0.5
133,2.25,4.1,0.5
*NGEN,NSET=BASE
1,7,2
*NGEN,NSET=BASER
21,27,2
*NGEN,NSET=TOP
111,113,2
*NGEN,NSET=TOPR
131,133,2
*NGEN, NSET=B1
11,17,2
*NGEN, NSET=B1R
31,37,2
*NGEN, NSET=T1
101,103
*NGEN, NSET=T1R
121,123
*NSET, NSET=B
BASE, B1
*NSET, NSET=BR
BASER, B1R
*NSET, NSET=T
TOP, T1
*NSET, NSET=TR
TOPR, T1R
*NCOPY, OLDSET=B, SHIFT, CHANGE=500
0., 0., 2.5

*NCOPY, OLDSET=BR, SHIFT, CHANGE=500
0., 0., 2.5

*NCOPY, OLDSET=T, SHIFT, CHANGE=500
0., 0., 1.5

*NCOPY, OLDSET=TR, SHIFT, CHANGE=500
0., 0., 1.5

*NSET, NSET=BASE, GEN
501,507,2
*NSET, NSET=BASER, GEN
521,527,2
*NSET, NSET=FACE1
1,11,501,511
*NSET, NSET=FACE1R
21,31,521,531
*NSET, NSET=TOP
611,613
*NSET, NSET=TOPR
631,633
*ELEMENT,TYPE=C3D8,ELSET=STEEL
1,1,3,13,11,501,503,513,511
11,101,103,113,111,601,603,613,611
21,21,23,33,31,521,523,533,531
31,121,123,133,131,621,623,633,631
*ELGEN,ELSET=STEEL
1,3,2
21,3,2
*ELSET,ELSET=EBOT,GENERATE
1,3,1
*ELSET,ELSET=EBOTR,GENERATE
21,23,1
*ELSET,ELSET=ETOP
11, 
*ELSET,ELSET=ETOPR
31, 
*surface,NAME=ASURF
ETOP,S3
*surface,NAME=BSURF
EBOT,S5
*CONTACT PAIR,INTERACTION=ROUGH,SMOOTH=0.1
ASURF,BSURF
*surface,NAME=ASURFR
ETOPR,S3
*surface,NAME=BSURFR
EBOTR,S5
*CONTACT PAIR,INTERACTION=ROUGH,SMOOTH=0.1
ASURFR,BSURFR
*SURFACE INTERACTION,NAME=ROUGH
*FRICTION
0.2, 
*SOLID SECTION,ELSET=STEEL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
 30.E6,0.3
**********
**
**
**
*STEP,UNSYMM=YES
  1 CONTACT -- Push top element down
*STATIC
 0.67,1.
*CONTROLS,PARAMETERS=FIELD
1.0E-6, 
*BOUNDARY
 BASE,2
 BASER,2
 FACE1,1
 FACE1R,1
 FACE1,3
 FACE1R,3
 TOP,2,,-.15
 TOPR,2,,-.15
 TOP,1
 TOPR,1
 TOP,3
 TOPR,3
*PRINT,CONTACT=YES
*EL PRINT,FREQ=50
S,
E,
*NODE PRINT,FREQ=50
U,
RF,
*EL FILE,FREQ=50
S,E
*NODE FILE,FREQ=50
U,RF
*CONTACT FILE, SECONDARY=ASURF, FREQ=50
CSTRESS,CDISP
CAREA, 
CFN,CFS
CMN,CMS
CFT,CMT
XN,XS,XT
*CONTACT PRINT, SECONDARY=ASURF, FREQ=50
CSTRESS,CDISP
CAREA, 
CFN,CFS
CMN,CMS
CFT,CMT
XN,XS,XT
*CONTACT FILE, SECONDARY=ASURFR, FREQ=50
CSTRESS,CDISP
CAREA, 
CFN,CFS
CMN,CMS
CFT,CMT
XN,XS,XT
*CONTACT PRINT, SECONDARY=ASURFR, FREQ=50
CSTRESS,CDISP
CAREA, 
CFN,CFS
CMN,CMS
CFT,CMT
XN,XS,XT
*END STEP
**
**
*STEP,UNSYMM=YES
  2 - Slide Top Across
*STATIC
 1.,10.,1.
*BOUNDARY,OP=NEW
 BASE,2
 FACE1,1
 FACE1,3
 TOP,2,,-.15
 TOP,1,,4.
 TOP,3
 BASER,2
 FACE1R,1
 FACE1R,3
 TOPR,2,,-.15
 TOPR,1,,4.
 TOPR,3
*END STEP
**
**
**
*STEP, UNSYMM=YES
 3 - Remove contact condition
*STATIC
1,10,1
*MODEL CHANGE, REMOVE, TYPE=CONTACT PAIR
 ASURF,BSURF
*BOUNDARY,OP=NEW
 BASE,2
 FACE1,1
 FACE1,3
 TOP,2,,-.15
 TOP,1,,4.
 TOP,3
 BASER,2
 FACE1R,1
 FACE1R,3
 TOPR,2,,-.099999999
 TOPR,1,,4.
 TOPR,3
*END STEP
**
**
*STEP, UNSYMM=YES,INC=50
 6 - reactivate contact gradually with contact interference
 The overclosure is not greater than 0.05
*STATIC
 .1, 10, .1
*MODEL CHANGE, ADD, TYPE=CONTACT PAIR
 ASURF,BSURF
*CONTACT INTERFER, TYPE=CONTACT PAIR
 ASURF, BSURF, 0.05
*BOUNDARY,OP=NEW
 BASE,2
 FACE1,1
 FACE1,3
 TOP,2,,-.15
 TOP,1,,4.
 TOP,3
 BASER,2
 FACE1R,1
 FACE1R,3
 TOPR,2,,-.15
 TOPR,1,,4.
 TOPR,3
*END STEP