*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
 WITH GENERALIZED PLANE STRAIN ELEMENTS
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 10: CPEG8
**
*NODE,NSET=NALL
701,  5.0, -27.0
702,  6.0, -27.0
703,  6.0, -28.0
704,  5.0, -28.0
705,  7.0, -27.0
706,  7.0, -28.0
707,  7.0, -27.5
708,  5.0, -27.5
709,  5.0, -30.0
710,  6.5, -30.1
711,  6.5, -31.1
712,  5.0, -31.0
713,  8.0, -30.2
714,  8.0, -31.2
715,  8.0, -30.7
716,  5.0, -30.5
717,  5.5, -26.0
719,  5.6, -26.0
*NSET,NSET=LEFT
701,704,708,709,712,716
*NSET,NSET=MID
710,711
*NSET,NSET=RIGHT
713,714,715
*NSET,NSET=TIP
705,706,707,713,714,715
*ELEMENT,TYPE=CPEG8,ELSET=TEST201
701, 701,704,706,705,708,703,707,702
*ELEMENT,TYPE=CPEG8,ELSET=TEST202
702, 709,712,714,713,716,711,715,710
*ELSET,ELSET=REFER
701,
*ELSET,ELSET=TEST
702,
**
** TEST 21: CPEG8H
**
*NODE,NSET=NALL
801,  5.0, -36.0
802,  6.0, -36.0
803,  6.0, -37.0
804,  5.0, -37.0
805,  7.0, -36.0
806,  7.0, -37.0
807,  7.0, -36.5
808,  5.0, -36.5
809,  5.0, -39.0
810,  6.5, -39.1
811,  6.5, -40.1
812,  5.0, -40.0
813,  8.0, -39.2
814,  8.0, -40.2
815,  8.0, -39.7
816,  5.0, -39.5
*NSET,NSET=LEFT
801,804,808,809,812,816
*NSET,NSET=MID
810,811
*NSET,NSET=RIGHT
813,814,815
*NSET,NSET=TIP
805,806,807,813,814,815
*ELEMENT,TYPE=CPEG8H,ELSET=TEST211
801, 801,804,806,805,808,803,807,802
*ELEMENT,TYPE=CPEG8H,ELSET=TEST212
802, 809,812,814,813,816,811,815,810
*ELSET,ELSET=REFER
801, 
*ELSET,ELSET=TEST
802, 
**
** TEST 22: CPEG8R
**
*NODE,NSET=NALL
901,  5.0, -45.0
902,  6.0, -45.0
903,  6.0, -46.0
904,  5.0, -46.0
905,  7.0, -45.0
906,  7.0, -46.0
907,  7.0, -45.5
908,  5.0, -45.5
909,  5.0, -48.0
910,  6.5, -48.1
911,  6.5, -49.1
912,  5.0, -49.0
913,  8.0, -48.2
914,  8.0, -49.2
915,  8.0, -48.7
916,  5.0, -48.5
*NSET,NSET=LEFT
901,904,908,909,912,916
*NSET,NSET=MID
910,911
*NSET,NSET=RIGHT
913,914,915
*NSET,NSET=TIP
905,906,907,913,914,915
*ELEMENT,TYPE=CPEG8R,ELSET=TEST221
901, 901,904,906,905,908,903,907,902
*ELEMENT,TYPE=CPEG8R,ELSET=TEST222
902, 909,912,914,913,916,911,915,910
*ELSET,ELSET=REFER
901,920
*ELSET,ELSET=TEST
902, 
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=REFER,REFNODE=717
1.0, 
*SOLID SECTION,MATERIAL=ELAS,ELSET=TEST,REFNODE=719
1.0, 
*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
LEFT,1,2
RIGHT,1,1,-1.0
RIGHT,2,2, 0.2
MID,1,1,-0.5
MID,2,2, 0.1
719,3,3,-0.1
719,4,4,-0.01
719,5,5, 0.01
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=999
E, S
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*ELEMENT OUTPUT
S, 
*ELEMENT OUTPUT,POSITION=CENTROID
E, 
S, 
MISESONLY,PRESSONLY
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF, 
*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,1,2
*END STEP
**
**
*STEP
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.1,1.
*DLOAD
TEST,P4,-0.25
REFER,P4,-0.25
*END STEP
**
**
*STEP
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.1,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #7: APPLY THERMAL LOAD
*STATIC
.1,1.0
*BOUNDARY,OP=NEW,FIXED
TIP,1
LEFT,1,2
*TEMPERATURE
NALL,100.
*EL PRINT,FREQ=999
S,THE,TEMP
*EL FILE,FREQ=999
S,THE,TEMP
*END STEP
**
**
*STEP
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
32,
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.0004,.004
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=999
E,S
*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, 
MISESONLY,PRESSONLY
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
1000., 10000., 5
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=0
*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
 STEP #11: FREQUENCY EXTRACTION FOR SUBSPACE PROJECTION
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
124,
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*EL FILE ,FREQ=0
*NODE FILE ,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #12: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
1000., 10000., 5
*BOUNDARY,OP=NEW
LEFT,1,2
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=0
*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 #13: STEADY-STATE DYNAMICS,SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
1000., 10000., 5
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=0
*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
 STEP #14: DYNAMICS
*DYNAMIC
.0004, .004
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=999
S,E
*NODE PRINT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TIP,2,0.2
*END STEP