*HEADING
 AXISYMMETRIC MEMBRANE ELEMENTS: MODEL CHANGE VERIFICATION
*PREPRINT,MODEL=YES,HISTORY=NO
*****************
** TEST 1: MAX1
*****************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=MAX1,ELSET=EALL
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ELSET,ELSET=EALL1
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL1
.25,
*****************
** TEST 2: MAX2
*****************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=MAX2,ELSET=EALL
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ELSET,ELSET=EALL2
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL2
.25,
******************
** TEST 3: MGAX1
******************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=MGAX1,ELSET=EALL
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ELSET,ELSET=EALL3
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL3
.25,
******************
** TEST 4: MGAX2
******************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=MGAX2,ELSET=EALL
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ELSET,ELSET=EALL4
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL4
.25,
****************************
** TEST 5: MAX1 WITH REBAR
****************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=MAX1,ELSET=EALL
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ELSET,ELSET=EALL5
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL5
.25,
*REBAR LAYER,ANGULAR SPACING
REB5,0.1,180,0,ELAS,30,1
****************************
** TEST 6: MAX2 WITH REBAR
****************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=MAX2,ELSET=EALL
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ELSET,ELSET=EALL6
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL6
.25,
*REBAR LAYER,ANGULAR SPACING
REB6,0.1,19.2,0,ELAS,30,1
*****************************
** TEST 7: MGAX1 WITH REBAR
*****************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=MGAX1,ELSET=EALL
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ELSET,ELSET=EALL7
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL7
.25,
*REBAR LAYER,ANGULAR SPACING
REB7,0.1,10.43,0,ELAS,30,1
*****************************
** TEST 8: MGAX2 WITH REBAR
*****************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=MGAX2,ELSET=EALL
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ELSET,ELSET=EALL8
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL8
.25,
*REBAR LAYER,ANGULAR SPACING
REB8,0.1,7.17,0,ELAS,30,1
************************************************
**  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.
**
**
*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,2
NALL,5
BOT,1
MID,1,1,-.5
TOP1,1,1,-1.
*NODE PRINT,NSET=BOT,FREQ=999
RF,
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
S,E,IVOL,RBANG,RBROT
*EL FILE,ELSET=EOUT,REBAR,FREQ=99
S,E,IVOL,RBANG,RBROT
*NODE FILE,NSET=BOT,FREQ=999
RF,
*EL FILE,ELSET=EOUT,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
BOT,2
BOT,5
TOP,2
*END STEP
**
**
*STEP,NLGEOM,INC=100
 STEP #4: BOUNDARY DISPLACEMENT
*STATIC
0.025,1.
*BOUNDARY,TYPE=VELOCITY
NALL,1,1,.05
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: FREEZE BOUNDARY DISPLACEMENT
*STATIC
.25,1.0
*BOUNDARY,OP=NEW,FIXED
NALL,1
BOT,2
BOT,5
TOP,2
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENTS
*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: 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
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
S,E,THE,RBANG,RBROT
*EL FILE,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,THE,E,TEMP
*EL FILE,REBAR,FREQ=99
S,E,RBANG,RBROT
*END STEP