*HEADING
 SUBSTRUCTURE ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=YES
**
**  GLOBAL USAGE
**
*NODE,NSET=NALL
2,  1.0, 1.
3,  2.0, 1.
12, 1.5, -.1
13, 3.0, -.3
*NSET,NSET=left
 2,12
*NSET,NSET=rightT
 13, 
*NSET,NSET=rightTR
 13,3
*ELEMENT,TYPE=Z01,ELSET=REFER, FILE=pmcp_substr_gen1
1, 2,3
*ELEMENT,TYPE=Z01,ELSET=TEST, FILE=pmcp_substr_gen1
2, 12,13
*SUBSTRUCTURE PROPERTY,ELSET=REFER
*SUBSTRUCTURE PROPERTY,ELSET=TEST
**
**
**
*STEP
 STEP #1: Arbitrarily DISPLACE TEST SPECIMEN  nodes
*STATIC
0.25,1.
*PRINT,MODEL CHANGE=YES
*BOUNDARY,op=new
left, encastre
rightT, 1,,-1.0
rightT, 2,, 2.
rightT, 6,,.1
*NODE PRINT,FREQ=999
U,RF
*NODE FILE,FREQ=999
RF, 
U, 
*OUTPUT,FIELD,VAR=PRESELECT
*END STEP
**
**
*STEP
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
left,encastre
**************************************** Substructure recovery
*SUBSTRUCTURE PATH,ENTER ELEMENT=1
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*EL FILE,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,VAR=PRESELECT
*SUBSTRUCTURE PATH,LEAVE
*SUBSTRUCTURE PATH,ENTER ELEMENT=2
*EL FILE,POSITION=CENTROID,FREQ=999
1,
S,E
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,VAR=PRESELECT
*SUBSTRUCTURE PATH,LEAVE
****************************************
*END STEP
**
**
*STEP
 STEP #4: APPLY INCREMENTAL DISPLACEMENT TO BOTH ELEMENTS 
*STATIC
0.25,1.
*BOUNDARY,OP=NEW,TYPE=VELOCITY
left, encastre
rightTR, 1,,-1.0
rightTR, 2,, 2.
rightTR, 6,,.1
*END STEP
**
**
*STEP
 STEP #5: REMOVE TEST ELEMENT, keeping nodes fixed
*STATIC
.25,1.0
*BOUNDARY,OP=NEW,FIX
left,encastre
rightTR,encastre
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.25,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
**
*STEP
 STEP #7:  RESTORE SIMPLE BC AND APPLY SLOAD CASE
*STATIC
.25,1.0
*BOUNDARY,op=new
left,encastre
*SLOAD
TEST, TEMPREF, 1.0
REFER, TEMPREF, 1.0
TEST, GRAVREF,1000.
REFER,GRAVREF,1000.
*END STEP
**