*HEADING
 AXISYMMETRIC MEMBRANE ELEMENTS: MODEL CHANGE VERIFICATION
*RESTART,WRITE
*PREPRINT,MODEL=YES
****************
** TEST 2: SAXA21 
****************
*NODE,NSET=NALL
501,  26.0, 1.
502,  26.5, 2.
503,  27.0, 3.
511,  26.0, 1.
512,  26.5, 2.
513,  27.0, 3.
*nset,nset=nall2
nall,1501,1502,1503,1511,1512,1513
*NSET,NSET=BOT
501,511,1501,1511
*NSET,NSET=MID1
512,1512
*NSET,NSET=TOP1
513,1513
*NSET,NSET=TOP
503,513,1503,1513
*ELEMENT,TYPE=SAXA21,ELSET=EALL,OFFSET=1000
501, 501,502,503
511, 511,512,513
*ELSET,ELSET=REFER
501, 
*ELSET,ELSET=TEST
511, 
*ELSET,ELSET=EOUT
501,511
************************************************
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
************************************************
*ELSET,ELSET=EALL
TEST,REFER
*MATERIAL,NAME=ELAS
*ELASTIC
1000., .25
*EXPANSION
1.E-3, 
*DENSITY
1.E-3, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.
*AMPLITUDE,NAME=RAMP
0.,0.,1.,1.
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL
.25, 
**
**
**
*STEP
 STEP #1: Arbitrary displacement to test element
*STATIC
1.,1.
*BOUNDARY,OP=NEW
NALL,2
BOT,1
MID1,1,1,-.5
TOP1,1,1,-1.
*NODE PRINT,NSET=BOT,FREQ=999
RF, 
U, 
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT,NSET=BOT
RF, 
U, 
*ELEMENT OUTPUT,ELSET=EOUT
E, S
*END STEP
**
**
*STEP
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE, with arbitrary displacement active
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW,FIXED
bot,1, 2
bot,6
top,1, 2
top,6
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT,NSET=BOT
RF, 
*END STEP
**
**
**
*STEP,INC=100
 STEP #4: DISTRIBUTED LOAD
*STATIC
0.025,1.
*DLOAD
TEST, P, -1.
REFER,P, -1.
*END STEP
**
**
*STEP
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.25,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,INC=100
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.025,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
*STEP
 STEP #7: APPLY DLOAD TYPE BX
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST,  BX, 10.
REFER, BX, 10. 
TEST, BZ, 10.
REFER,BZ, 10.
*END STEP
**
**
*STEP
 STEP #8: APPLY THERMAL LOAD
*STATIC
.25,1.0
*TEMPERATURE
NALL,100.
*NODE PRINT,FREQ=0
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,THE,
E,TEMP
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=EOUT
S,THE,E,TEMP
*END STEP
**
**