*HEADING
 TEST OF REMOVING SUBSTRUCTURES: EXCAVATION USING SUBSTRUCTURES
******************
** GLOBAL MODEL **
******************
*NODE,NSET=NALL
21,    0.,   -66.66
25,  100.,   -66.66
41,    0.,   -33.33
45,  100.,   -33.33
61,    0.,     0.0
65,  100.,     0.0
*NGEN,NSET=NALL
21,25,1
41,45,1
61,65,1
*ELEMENT,TYPE=Z1,ELSET=LAYER1, FILE=psuprem1_top
1, 21,22,23,24,25 
*ELEMENT,TYPE=Z2,ELSET=LAYER2, FILE=psuprem1_middle
2, 21,22,23,24,25,41,42,43,44,45 
*ELEMENT,TYPE=Z3,ELSET=LAYER3, FILE=psuprem1_bottom
3, 41,42,43,44,45,61,62,63,64,65 
*SUBSTRUCTURE PROPERTY,ELSET=LAYER1
*SUBSTRUCTURE PROPERTY,ELSET=LAYER2
*SUBSTRUCTURE PROPERTY,ELSET=LAYER3
*RESTART,WRITE
**********************
** APPLY BODY LOADS **
**********************
*STEP
 APPLY BODY LOADS
*STATIC
1.0,1.0
*SLOAD
1,BOUND,1.0
2,BOUND,1.0
3,BOUND,1.0
*OUTPUT,FIELD,VAR=PRESELECT
*SUBSTRUCTURE PATH,ENTER ELEMENT=1
*OUTPUT,FIELD,VAR=PRESELECT
*EL FILE,FREQ=999
COORD,S
*EL PRINT
COORD,S
*NODE FILE,FREQ=999
U,
*NODE PRINT
U,
*SUBSTRUCTURE PATH,LEAVE
*SUBSTRUCTURE PATH,ENTER ELEMENT=2
*OUTPUT,FIELD,VAR=PRESELECT
*EL PRINT
COORD,S
*EL FILE,FREQ=999
COORD,S
*NODE FILE,FREQ=999
U,
*NODE PRINT
U,
*SUBSTRUCTURE PATH,LEAVE
*SUBSTRUCTURE PATH,ENTER ELEMENT=3
*OUTPUT,FIELD,VAR=PRESELECT
*EL PRINT
COORD,S
*EL FILE,FREQ=999
COORD,S
*NODE PRINT
U,
*NODE FILE,FREQ=999
U,
*SUBSTRUCTURE PATH,LEAVE
*EL PRINT,FREQ=0
*NODE PRINT
U, 
*END STEP
*STEP
 REMOVE TOP LAYER
*STATIC
1.0,1.0
*MODEL CHANGE,REMOVE
LAYER3,
*SUBSTRUCTURE PATH,ENTER ELEMENT=1
*OUTPUT,FIELD,VAR=PRESELECT
*EL PRINT
COORD,S
*EL FILE,FREQ=999
COORD,S
*NODE PRINT
U,
*NODE FILE,FREQ=999
U,
*SUBSTRUCTURE PATH,LEAVE
*SUBSTRUCTURE PATH,ENTER ELEMENT=2
*OUTPUT,FIELD,VAR=PRESELECT
*EL PRINT
COORD,S
*EL FILE,FREQ=999
COORD,S
*NODE PRINT
U,
*NODE FILE,FREQ=999
U,
*SUBSTRUCTURE PATH,LEAVE
*NODE PRINT
U,
*END STEP
*STEP
 REMOVE MIDDLE LAYER
*STATIC
1.0,1.0
*MODEL CHANGE,REMOVE
LAYER2,
*SUBSTRUCTURE PATH,ENTER ELEMENT=1
*OUTPUT,FIELD,VAR=PRESELECT
*EL PRINT
COORD,S
*EL FILE,FREQ=999
COORD,S
*NODE PRINT
U, 
*NODE FILE,FREQ=999
U, 
*SUBSTRUCTURE PATH,LEAVE
*EL PRINT,FREQ=0
*NODE PRINT
U, 
*END STEP