*HEADING
 :  C3D8 MODEL CHANGE VERIFICATION; ELASTIC MATERIAL; THERMAL EXPANSION
*NODE
1,
11,5.
401,,2.
411,5.,2.
1001, ,   , -1.0
1011,5.,  , -1.0
1401,,2.  , -1.0
1411,5.,2., -1.0
*NGEN,NSET=BOT
1,11,1
1001,  1011
*NGEN,NSET=TOP
401,411,1
1401, 1411, 1
*NFILL,NSET=NALL
BOT,TOP,4,100
**
** Duplicate nodes at the boundary of removed elements.
** This allows independent temp. loading at the boundary
** yielding uniform temp. distribution within each element.
*NODE,NSET=NSIDES
2003,1.,0.,0.
2009,4.,0.,0.
2203,1.,1.,0.
2209,4.,1.,0.
2403,1.,2.,0.
2409,4.,2.,0.
3003,1.,0.,-1.
3009,4.,0.,-1.
3203,1.,1.,-1.
3209,4.,1.,-1.
3403,1.,2.,-1.
3409,4.,2.,-1.
*NSET,NSET=NALL
NSIDES,
*NSET, NSET=XFIX, GEN
1, 401, 100
1001, 1401, 100
*NSET, NSET=YZFIX
1, 11
*NSET, NSET=YFIX
1001, 1011
*NSET, NSET=XMOVE, GEN
11, 411, 100
1011, 1411, 100
*NSET, NSET=YZMOVE, GEN
204, 208, 1
404, 408, 1
1204, 1208,1
1404, 1408,1
1004, 1008, 1 
*NSET,NSET=NMIDTMP
3,203,403
5,205,405
7,207,407
9,209,409
1003,1203,1403
1005,1205,1405
1007,1207,1407
1009,1209,1409
*NSET,NSET=NSIDTMP
NSIDES,XMOVE,XFIX
*ELEMENT,TYPE=C3D8
1,1,2003,3003,1001, 201, 2203, 3203, 1201,
2,3,5,1005,1003, 203,205,1205,1203,
5,2009,11,1011,3009,2209,211,1211,3209 
6,201,2203,3203,1201,401,2403,3403,1401
7,203,205,1205,1203,403,405,1405,1403
10,2209,211,1211,3209,2409,411,1411,3409
*ELGEN,ELSET=EREMOV
2,3,2,1
7,3,2,1
*ELSET, ELSET=ESIDES
1,6,5,10
*SOLID SECTION,MATERIAL=MRUBBER,ELSET=EREMOV
*MATERIAL,NAME=MRUBBER
*ELASTIC
4064385.018,0.451566
*EXPANSION
1.E-4, 
*SOLID SECTION,MATERIAL=SRUBBER,ELSET=ESIDES
*MATERIAL,NAME=SRUBBER
*ELASTIC
4064385.018,0.451566
*EXPANSION
2.E-4, 
**
** MPC: There are duplicate nodes at the boundary of element
** removal for temperature loading purposes.  This yields an 
** unconnected mesh.  MPC makes sure we have a rooted-level
** structure.
*NSET,NSET=NMIDDLE
3,9,203,209,403,409,1003,1009,1203,1209,1403,1409
*MPC
TIE,NSIDES,NMIDDLE
*BOUNDARY
NALL,1,3
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,20.0
*STEP,NLGEOM
 REMOVE ELEMENTS
*STATIC
0.1,1.,0.1,0.1
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
EREMOV,
*NODE FILE,FREQ=10
U, RF
*EL FILE,FREQ=10
E,EE, THE, S
EVOL, 
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
EE,THE,S,EVOL
*NODE OUTPUT
U,RF
*END STEP
*STEP,NLGEOM
 HEAT MIDDLE NODES TO 100 DEG.
*STATIC
0.1,1.,0.1,0.1
*BOUNDARY,OP=NEW
YZFIX, 2, 3
YFIX, 2, 2
XFIX, 1,1
XMOVE,1,1
*TEMPERATURE
NMIDTMP,100.
*ENDSTEP
*STEP,NLGEOM
 REINTRODUCE MIDDLE ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,ADD,TYPE=ELEMENT
EREMOV,
*BOUNDARY,OP=NEW
YZFIX, 2, 3
YFIX, 2, 2
XFIX, 1,1
*END STEP
*STEP,NLGEOM
 RAISE TEMPERATURE TO GET UNIFORM THERMAL EXPANSION
*STATIC
1.,1.
*TEMPERATURE
NMIDTMP,200.
NSIDTMP,70.
*ENDSTEP