*HEADING
 AXISYMMETRIC MEMBRANE ELEMENTS: MODEL CHANGE VERIFICATION
*PREPRINT,MODEL=YES,HISTORY=NO
*****************
** TEST 1: MAX1
*****************
*NODE,NSET=NALL
1,   1.0, 1.
2,   1.5, 2.
3,   2.0, 3.
11,  1.0, 5.
12,  2.0, 6.
13,  3.0, 7.
*NSET,NSET=BOT
1,11
*NSET,NSET=MID
12,
*NSET,NSET=TOP1
13,
*NSET,NSET=TOP
3,13
*ELEMENT,TYPE=MAX1,ELSET=EALL
1,  1,2
2,  2,3
11, 11,12
12, 12,13
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
11,12
*ELSET,ELSET=EOUT
1,11
*ELSET,ELSET=EALL1
1,2,11,12
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL1
.25, 
*****************
** TEST 2: MAX2
*****************
*NODE,NSET=NALL
101,   6.0, 1.
102,   6.5, 2.
103,   7.0, 3.
111,   6.0, 5.
112,   7.0, 6.
113,   8.0, 7.
*NSET,NSET=BOT
101,111,
*NSET,NSET=MID
112,
*NSET,NSET=TOP1
113,
*NSET,NSET=TOP
103,113
*ELEMENT,TYPE=MAX2,ELSET=EALL
101, 101,102,103
111, 111,112,113
*ELSET,ELSET=REFER
101,
*ELSET,ELSET=TEST
111,
*ELSET,ELSET=EOUT
101,111
*ELSET,ELSET=EALL2
101,111
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL2
.25, 
******************
** TEST 3: MGAX1
******************
*NODE,NSET=NALL
201,  11.0, 1.
202,  11.5, 2.
203,  12.0, 3.
211,  11.0, 5.
212,  12.0, 6.
213,  13.0, 7.
*NSET,NSET=BOT
201,211,
*NSET,NSET=MID
212,
*NSET,NSET=TOP1
213,
*NSET,NSET=TOP
203,213
*ELEMENT,TYPE=MGAX1,ELSET=EALL
201, 201,202
202, 202,203
211, 211,212
212, 212,213
*ELSET,ELSET=REFER
201,202
*ELSET,ELSET=TEST
211,212
*ELSET,ELSET=EOUT
201,211
*ELSET,ELSET=EALL3
201,202,211,212
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL3
.25, 
******************
** TEST 4: MGAX2
******************
*NODE,NSET=NALL
301,  16.0, 1.
302,  16.5, 2.
303,  17.0, 3.
311,  16.0, 5.
312,  17.0, 6.
313,  18.0, 7.
*NSET,NSET=BOT
301,311,
*NSET,NSET=MID
312,
*NSET,NSET=TOP1
313,
*NSET,NSET=TOP
303,313
*ELEMENT,TYPE=MGAX2,ELSET=EALL
301, 301,302,303
311, 311,312,313
*ELSET,ELSET=REFER
301, 
*ELSET,ELSET=TEST
311, 
*ELSET,ELSET=EOUT
301,311
*ELSET,ELSET=EALL4
301,311
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL4
.25, 
****************************
** TEST 5: MAX1 WITH REBAR
****************************
*NODE,NSET=NALL
401,   1.0, 10.
402,   1.5, 11.
403,   2.0, 12.
411,   1.0, 14.
412,   2.0, 15.
413,   3.0, 16.
*NSET,NSET=BOT
401,411,
*NSET,NSET=MID
412,
*NSET,NSET=TOP1
413,
*NSET,NSET=TOP
403,413
*ELEMENT,TYPE=MAX1,ELSET=EALL
401,  401,402
402,  402,403
411,  411,412
412,  412,413
*ELSET,ELSET=REFER
401,402
*ELSET,ELSET=TEST
411,412
*ELSET,ELSET=EOUT
401,411
*ELSET,ELSET=EALL5
401,402,411,412
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL5
.25,
*REBAR LAYER,GEOMETRY=ANGULAR
REB5,0.1,180,0,ELAS,30,1  
****************************
** TEST 6: MAX2 WITH REBAR
****************************
*NODE,NSET=NALL
501,   6.0, 10.
502,   6.5, 11.
503,   7.0, 12.
511,   6.0, 14.
512,   7.0, 15.
513,   8.0, 16.
*NSET,NSET=BOT
501,511,
*NSET,NSET=MID
512,
*NSET,NSET=TOP1
513,
*NSET,NSET=TOP
503,513
*ELEMENT,TYPE=MAX2,ELSET=EALL
501, 501,502,503
511, 511,512,513
*ELSET,ELSET=REFER
501,
*ELSET,ELSET=TEST
511,
*ELSET,ELSET=EOUT
501,511
*ELSET,ELSET=EALL6
501,511
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL6
.25,
*REBAR LAYER,GEOMETRY=ANGULAR
REB6,0.1,19.2,0,ELAS,30,1   
*****************************
** TEST 7: MGAX1 WITH REBAR
*****************************
*NODE,NSET=NALL
601,  11.0, 10.
602,  11.5, 11.
603,  12.0, 12.
611,  11.0, 14.
612,  12.0, 15.
613,  13.0, 16.
*NSET,NSET=BOT
601,611,
*NSET,NSET=MID
612,
*NSET,NSET=TOP1
613,
*NSET,NSET=TOP
603,613
*ELEMENT,TYPE=MGAX1,ELSET=EALL
601, 601,602
602, 602,603
611, 611,612
612, 612,613
*ELSET,ELSET=REFER
601,602
*ELSET,ELSET=TEST
611,612
*ELSET,ELSET=EOUT
601,611
*ELSET,ELSET=EALL7
601,602,611,612
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL7
.25, 
*REBAR LAYER,GEOMETRY=ANGULAR
REB7,0.1,10.43,0,ELAS,30,1  
*****************************
** TEST 8: MGAX2 WITH REBAR
*****************************
*NODE,NSET=NALL
701,  16.0, 10.
702,  16.5, 11.
703,  17.0, 12.
711,  16.0, 14.
712,  17.0, 15.
713,  18.0, 16.
*NSET,NSET=BOT
701,711,
*NSET,NSET=MID
712,
*NSET,NSET=TOP1
713,
*NSET,NSET=TOP
703,713
*ELEMENT,TYPE=MGAX2,ELSET=EALL
701, 701,702,703
711, 711,712,713
*ELSET,ELSET=REFER
701, 
*ELSET,ELSET=TEST
711, 
*ELSET,ELSET=EOUT
701,711
*ELSET,ELSET=EALL8
701,711
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=EALL8
.25,
*REBAR LAYER,GEOMETRY=ANGULAR
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
1, 
S,E,THE,RBANG,RBROT
*EL FILE,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,THE,E,TEMP
*EL FILE,REBAR,FREQ=99
1, 
S,E,RBANG,RBROT
*END STEP