*HEADING
 MEMBRANE ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*PREPRINT,MODEL=YES,HISTORY=NO
*****************
** TEST 1: CINAX4
*****************
*NODE,NSET=NALL
101,  4.0,  0.0,
102,  5.0,  0.0,
103,  5.0,  1.0,
104,  4.0,  1.0,
111,  4.0,  3.0,
112,  5.0,  3.0,
113,  5.0,  3.9,
114,  4.0,  4.0,
*NSET,NSET=BOT
101,102,111,112
*NSET,NSET=LEFT
101,104,111,114
*NSET,NSET=TOP
103,104,113,114
*NSET,NSET=RIGHT
102,103,112,113
*NSET,NSET=TIP
113, 
*ELEMENT,TYPE=CINAX4,ELSET=TEST2
101, 103,104,101,102
111, 113,114,111,112
*NSET,NSET=RIGHT
3,4,6,7,13,14,16,17
*ELSET,ELSET=REFER2
101,
*ELSET,ELSET=TEST
111, 
*****************
** TEST 2: CINPS4
*****************
*NODE,NSET=NALL
201,  8.0,  0.0,
202,  9.0,  0.0,
203,  9.0,  1.0,
204,  8.0,  1.0,
211,  8.0,  3.0,
212,  9.0,  3.0,
213,  9.0,  3.9,
214,  8.0,  4.0,
*NSET,NSET=BOT
201,202,211,212
*NSET,NSET=LEFT
201,204,211,214
*NSET,NSET=TOP
203,204,213,214
*NSET,NSET=RIGHT
202,203,212,213
*NSET,NSET=TIP
213, 
*ELEMENT,TYPE=CINPS4,ELSET=TEST3
201, 203,204,201,202
211, 213,214,211,212
*ELSET,ELSET=REFER3
201,
*ELSET,ELSET=TEST
211,
*ELSET,ELSET=REFER
REFER2,REFER3
************************************************
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
************************************************
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL2
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL3
*ELSET,ELSET=EALL2
TEST2,REFER2
*ELSET,ELSET=EALL3
TEST3,REFER3
*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.
**********************************************************
*STEP,NLGEOM
 STEP #1: DISPLACE TEST SPECIMEN SO FINAL SHAPE IS
          IDENTICAL TO INITIAL SHAPE OF REFERENCE ELEMENT 
**********************************************************
*STATIC
0.25,1.
*PRINT,MODEL CHANGE=YES
*BOUNDARY,OP=NEW
NALL,3
LEFT,1,2
BOT,1,2
RIGHT,1
TIP,2,2, 0.1
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*NODE FILE,NSET=LEFT,FREQ=999
RF, 
*EL PRINT,ELSET=EALL,POSITION=CENTROID,FREQ=999
S,
E,
*EL FILE,ELSET=EALL,POSITION=CENTROID,FREQ=999
E, S
*END STEP
**
**
*STEP,NLGEOM
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
LEFT,1
NALL,3
BOT,2
*END STEP
**
**
*STEP,NLGEOM,INC=100
 STEP #4: DISPLACE ELEMENTS
*STATIC
0.025,1.
*BOUNDARY,TYPE=VELOCITY
TOP, 2,2,-.01
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.25,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM,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,NLGEOM
 STEP #7: DISPLACE ELEMENTS FURTHER
*STATIC
.25,1.0
*BOUNDARY,TYPE=VELOCITY
TOP, 2,2,-.02
*END STEP