*HEADING
 :  CPE4I MODEL CHANGE VERIFICATION, ELASTIC MATERIAL, REBAR
*NODE
1,
11,5.
401,,2.
411,5.,2.
*NGEN,NSET=BOT
1,11,1
*NGEN,NSET=TOP
401,411,1
*NFILL,NSET=NALL
BOT,TOP,4,100
*NSET, NSET=XFIX, GEN
1, 401, 100
*NSET, NSET=YFIX
1,3,5
*NSET, NSET=YFIX2
7,9,11
*NSET, NSET=XMOVE, GEN
11, 411, 100
*NSET, NSET=YMOVE, GEN
204, 208, 1
404, 408, 1
*ELEMENT,TYPE=CPE4I,ELSET=EALL
1,1,3,203,201
*ELGEN,ELSET=EALL
1,5,2,1,2,200,5
*ELSET, ELSET=EREMOV
2,3,4,7,8,9
*SOLID SECTION,MATERIAL=RUBBER,ELSET=EALL
1.,
*MATERIAL,NAME=RUBBER
*ELASTIC
4.064385E8,0.451566
** (This gives a plane strain modulus of 5.10545E8)
**
*MATERIAL, NAME=MREBAR
*ELASTIC
 5.10545E10,0
**
**
*REBAR,ELEMENT=CONTINUUM, MAT=MREBAR,NAME=REB
 EALL,     0.01,   1.   ,        0,  .5      , 1
** el,  x-sec-a, spacing, oriangle, frac-dist, edgenum, 
**
**
*BOUNDARY
YFIX, 2, 2
XFIX, 1,1
YFIX2, 2,2
*STEP,NLGEOM
 COMPRESS ELEMENTS
*STATIC
0.1,1.,0.1,0.1
*BOUNDARY
XMOVE, 1,1 ,-0.1E-2
*NODE PRINT,FREQ=999
U, 
RF, 
*EL PRINT,FREQ=999
E,S
*EL PRINT,REBAR, FREQ=999
E,S,NE
*NODE FILE,FREQ=999
U, RF
*EL FILE,FREQ=999
E, NE, S
*END STEP
*STEP,NLGEOM
 REMOVE MIDDLE ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
EREMOV,
*END STEP
*STEP,NLGEOM
 MOVE REMOVED NODES BACK IN 
*STATIC
1.,1.
*BOUNDARY
YMOVE, 2,2, 0.0
*END STEP
*STEP,NLGEOM
 ADD ELEMENTS BACK IN ANNEALED
*STATIC
0.1,1.,0.1,0.1
*MODEL CHANGE,ADD,TYPE=ELEMENT
EREMOV, 
*BOUNDARY,OP=NEW
YFIX, 2, 2
XFIX, 1,1
YFIX2, 2,2
XMOVE, 1,1 ,0.1E-2
*END STEP