*HEADING
 AXISYMMETRIC SHELL ELEMENTS: MODEL CHANGE VERIFICATION
*PREPRINT,MODEL=YES,HISTORY=NO
*****************
** TEST 1: SAX1
*****************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID1
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=SAX1,ELSET=EALL
*ELSET,ELSET=EALL1
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL1
.25,
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*****************
** TEST 2: SAX2
*****************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID1
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=SAX2,ELSET=EALL
*ELSET,ELSET=EALL2
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL2
.25,
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
**************************
** TEST 3: SAX1 WITH REBAR
**************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID1
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=SAX1,ELSET=EALL
*ELSET,ELSET=EALL3
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL3
.25,
*REBAR LAYER,ANGULAR SPACING
REB1,0.1,9.977,0,ELAS,30,1
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
**************************
** TEST 4: SAX2 WITH REBAR
**************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=MID1
*NSET,NSET=TOP1
*NSET,NSET=TOP
*ELEMENT,TYPE=SAX2,ELSET=EALL
*ELSET,ELSET=EALL4
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL4
.25,
*REBAR LAYER,ANGULAR SPACING
REB2,0.1,6.95,0,ELAS,30,1
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
************************************************
**  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
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
**
**
*STEP,NLGEOM
 STEP #10: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos
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
S,E
*EL FILE,ELSET=EOUT,REBAR,FREQ=99
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
S,E
*EL FILE,REBAR,FREQ=99
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
*elementoutput
inten,
*elementoutput,rebar
inten,
*NODE PRINT,FREQ=0
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
S,E
*EL FILE,REBAR,FREQ=99
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
*elementoutput
inten,
*elementoutput,rebar
inten,
*NODE PRINT,FREQ=0
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
S,E
*EL FILE,REBAR,FREQ=99
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
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
S,E,RBANG,RBROT
*EL FILE,REBAR,FREQ=99
S,E,RBANG,RBROT
*NODE PRINT,NSET=BOT,FREQ=999
RF,
*CLOAD,AMP=RAMP
TOP,1,150.
*END STEP