*HEADING
 AXISYMMETRIC SHELL ELEMENTS: MODEL CHANGE VERIFICATION
*PREPRINT,MODEL=YES,HISTORY=NO
*****************
** TEST 1: SAX1
*****************
*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=MID1
12,
*NSET,NSET=TOP1
13,
*NSET,NSET=TOP
3,13
*ELEMENT,TYPE=SAX1,ELSET=EALL
1,  1,2
2,  2,3
11, 11,12
12, 12,13
*ELSET,ELSET=EALL1
1,2,11,12
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL1
.25,
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
11,12
*ELSET,ELSET=EOUT
1,11
*****************
** TEST 2: SAX2
*****************
*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=MID1
112,
*NSET,NSET=TOP1
113,
*NSET,NSET=TOP
103,113
*ELEMENT,TYPE=SAX2,ELSET=EALL
101, 101,102,103
111, 111,112,113
*ELSET,ELSET=EALL2
101,111
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL2
.25,
*ELSET,ELSET=REFER
101,
*ELSET,ELSET=TEST
111,
*ELSET,ELSET=EOUT
101,111
**************************
** TEST 3: SAX1 WITH REBAR
**************************
*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=MID1
212,
*NSET,NSET=TOP1
213,
*NSET,NSET=TOP
203,213
*ELEMENT,TYPE=SAX1,ELSET=EALL
201, 201,202
202, 202,203
211, 211,212
212, 212,213
*ELSET,ELSET=EALL3
201,202,211,212
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL3
.25,
*REBAR LAYER,GEOMETRY=ANGULAR
REB1,0.1,9.977,0,ELAS,30,1   
*ELSET,ELSET=REFER
201,202
*ELSET,ELSET=TEST
211,212
*ELSET,ELSET=EOUT
201,211
**************************
** TEST 4: SAX2 WITH REBAR
**************************
*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=MID1
312,
*NSET,NSET=TOP1
313,
*NSET,NSET=TOP
303,313
*ELEMENT,TYPE=SAX2,ELSET=EALL
301, 301,302,303
311, 311,312,313
*ELSET,ELSET=EALL4
301,311
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL4
.25,
*REBAR LAYER,GEOMETRY=ANGULAR
REB2,0.1,6.95,0,ELAS,30,1   
*ELSET,ELSET=REFER
301, 
*ELSET,ELSET=TEST
311, 
*ELSET,ELSET=EOUT
301,311
************************************************
**  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
MID1,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
*OUTPUT,HISTORY,VARIABLE=PRESELECT,FREQ=1
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
E,S
*ELEMENT OUTPUT,REBAR
S,E
*NODE OUTPUT
U,
RF,
*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: DISTRIBUTED LOAD
*STATIC
0.025,1.
*DLOAD
TEST, P, -1.
REFER,P, -1.
*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: APPLY DLOAD TYPE BX
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST,  BR, 10.
REFER, BR, 10. 
TEST, BZ, 10.
REFER,BZ, 10.
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: APPLY DLOAD TYPE BR
*STATIC
.25,1.0
*DLOAD, OP=NEW
**TEST,  GRAV, 15000., 1., 0.,
**REFER, GRAV, 15000., 1., 0.,
TEST,  BR, 15.
REFER, BR, 15.
*END STEP
**
**
*STEP,NLGEOM
 STEP #9: 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
**
**
*STEP,NLGEOM
 STEP #10: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
20, , 
*EL PRINT,ELSET=EOUT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*MODAL FILE 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: MODAL DYNAMIC
*MODAL DYNAMIC
.0025,.025
*CLOAD
TOP,1,150.
*NODE PRINT,FREQ=0
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
1, 
S,E
*EL FILE,ELSET=EOUT,REBAR,FREQ=99
1, 
S,E
*NODE FILE,NSET=BOT,FREQ=999
RF, 
*EL FILE,ELSET=EOUT,POSITION=CENTROID,FREQ=999
E,S
*OUTPUT,HISTORY,VARIABLE=PRESELECT,FREQ=1
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT
E,S
*ELEMENT OUTPUT,REBAR
S,E
*NODE OUTPUT
U,
RF,
*END STEP
**
**
*STEP,PERTURBATION
 STEP #12: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
5.,25., 5
*CLOAD
TOP,1,150.
*NODE PRINT,FREQ=0
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
1, 
S,E
*EL FILE,REBAR,FREQ=99
1, 
S,E
*NODE FILE,NSET=BOT,FREQ=999
RF, 
*EL FILE,ELSET=EOUT,POSITION=CENTROID,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #13: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
5., 25., 5
*CLOAD
TOP,1,150.
*output,field
*NODE PRINT,FREQ=0
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
1, 
S,E
*EL FILE,REBAR,FREQ=99
1, 
S,E
*NODE FILE,NSET=BOT,FREQ=999
RF, 
*EL FILE,ELSET=EOUT,POSITION=CENTROID,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #14: STEADY-STATE DYNAMICS, SUPSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
5., 25., 5
*CLOAD
TOP,1,150.
*output,field
*NODE PRINT,FREQ=0
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
1, 
S,E
*EL FILE,REBAR,FREQ=99
1, 
S,E
*NODE FILE,NSET=BOT,FREQ=999
RF, 
*EL FILE,ELSET=EOUT,POSITION=CENTROID,FREQ=0
*END STEP
**
**
*STEP,NLGEOM
 STEP #14: DYNAMICS
*DYNAMIC
.0025,.025
*NODE PRINT,FREQ=0
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
1, 
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
1, 
S,E,RBANG,RBROT
*EL FILE,REBAR,FREQ=99
1, 
S,E,RBANG,RBROT
*NODE PRINT,NSET=BOT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TOP,1,150.
*END STEP