*HEADING
 JOINTC ELEMENT: MODEL CHANGE VERIFICATION
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 1: JOINTC
**
*NODE,NSET=NALL
1,
2,
11,
12,.1,.1
*NSET,NSET=LEFT
1,11
*NSET,NSET=RIGHT
12, 
*NSET,NSET=TIP
2,12
*ELEMENT,TYPE=JOINTC,ELSET=J
1, 1,2
11,11,12
*JOINT,ELSET=J
*SPRING
1,
100.,
*SPRING
2,
200.,
*SPRING
3,
300.,
*SPRING
4,
400.
*SPRING
5,
500.,
*SPRING
6,
600.,
*ELEMENT,TYPE=MASS,ELSET=MASS
201, 2
211, 12
*MASS,ELSET=MASS
10., 
*ELSET,ELSET=TEST
11,211
*ELSET,ELSET=TEST1
11,
*ELSET,ELSET=TEST2
211,
*AMPLITUDE,NAME=RAMP
0.,0., .5,1.0
**
**
*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,AMP=RAMP
NALL,3
LEFT,1,6
RIGHT,1,1,-0.1
RIGHT,2,2,-0.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*ELEMENT OUTPUT,POSITION=CENTROID
E,
S,
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
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
LEFT,1,6
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.25,1.
*BOUNDARY,AMP=RAMP,TYPE=VELOCITY
2,1,1,.3
2,2,2,.3
12,1,1,.3
12,2,2,.3
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.25,1.0
*BOUNDARY,OP=NEW,FIXED
NALL,1,3
LEFT,6,6
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.25,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST1,
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST2,
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos
8, , 
*BOUNDARY,OP=NEW
LEFT,1,6
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #8: MODAL DYNAMIC
*MODAL DYNAMIC
.5,5.
*CLOAD
TIP,1,2.5
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
E, 
S, 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
.05,.5,5
*CLOAD
TIP,1,2.5
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF, 
*OUTPUT,FIELD,FREQ=0
*ELEMENT OUTPUT,POSITION=CENTROID
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
.05,.5,5
*BOUNDARY,OP=NEW
LEFT,1,6
*CLOAD
TIP,1,2.5
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=NALL
RF, 
*OUTPUT,FIELD, FREQ=0
*ELEMENT OUTPUT,POSITION=CENTROID
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
.05,.5,5
*CLOAD
TIP,1,2.5
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=NALL
RF, 
*OUTPUT,FIELD, FREQ=0
*ELEMENT OUTPUT,POSITION=CENTROID
*END STEP
**
**
*STEP,NLGEOM
 STEP #12: DYNAMICS
*DYNAMIC
.5,5.
*BOUNDARY,OP=NEW
LEFT,1,6
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*NODE PRINT,FREQ=999
RF, 
*CLOAD
TIP,1,2.5
*END STEP