*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 8: CPS4
**
*NODE,NSET=NALL
1001,  15.0,  0.0
1002,  16.0,  0.0
1003,  16.0, -1.0
1004,  15.0, -1.0
1005,  17.0,  0.0
1006,  17.0, -1.0
1007,  15.0, -3.0
1008,  16.5, -3.1
1009,  16.5, -4.1
1010,  15.0, -4.0
1011,  18.0, -3.2
1012,  18.0, -4.2
*NSET,NSET=LEFT
1001,1004,1007,1010
*NSET,NSET=MID
1008,1009
*NSET,NSET=RIGHT
1011,1012
*NSET,NSET=TIP
1005,1006,1011,1012
*ELEMENT,TYPE=CPS4,ELSET=TEST8
1001, 1001,1004,1003,1002
1002, 1002,1003,1006,1005
1003, 1007,1010,1009,1008
1004, 1008,1009,1012,1011
*ELSET,ELSET=REFER
1001,1002
*ELSET,ELSET=TEST
1003,1004
**
** TEST 9: CPS4I
**
*NODE,NSET=NALL
2001,  15.0, -9.0
2002,  16.0, -9.0
2003,  16.0, -10.0
2004,  15.0, -10.0
2005,  17.0, -9.0
2006,  17.0, -10.0
2007,  15.0, -12.0
2008,  16.5, -12.1
2009,  16.5, -13.1
2010,  15.0, -13.0
2011,  18.0, -12.2
2012,  18.0, -13.2
*NSET,NSET=LEFT
2001,2004,2007,2010
*NSET,NSET=MID
2008,2009
*NSET,NSET=RIGHT
2011,2012
*NSET,NSET=TIP
2005,2006,2011,2012
*ELEMENT,TYPE=CPS4I,ELSET=TEST9
2001, 2001,2004,2003,2002
2002, 2002,2003,2006,2005
2003, 2007,2010,2009,2008
2004, 2008,2009,2012,2011
*ELSET,ELSET=REFER
2001,2002
*ELSET,ELSET=TEST
2003,2004
**
** TEST 10: CPS4R
**
*NODE,NSET=NALL
3001,  15.0, -18.0
3002,  16.0, -18.0
3003,  16.0, -19.0
3004,  15.0, -19.0
3005,  17.0, -18.0
3006,  17.0, -19.0
3007,  15.0, -21.0
3008,  16.5, -21.1
3009,  16.5, -22.1
3010,  15.0, -22.0
3011,  18.0, -21.2
3012,  18.0, -22.2
*NSET,NSET=LEFT
3001,3004,3007,3010
*NSET,NSET=MID
3008,3009
*NSET,NSET=RIGHT
3011,3012
*NSET,NSET=TIP
3005,3006,3011,3012
*ELEMENT,TYPE=CPS4R,ELSET=TEST10
3001, 3001,3004,3003,3002
3002, 3002,3003,3006,3005
3003, 3007,3010,3009,3008
3004, 3008,3009,3012,3011
*ELSET,ELSET=REFER
3001,3002
*ELSET,ELSET=TEST
3003,3004
**
** TEST 11: CPS8
**
*NODE,NSET=NALL
4001,  25.0,  0.0
4002,  26.0,  0.0
4003,  26.0, -1.0
4004,  25.0, -1.0
4005,  27.0,  0.0
4006,  27.0, -1.0
4007,  27.0, -0.5
4008,  25.0, -0.5
4009,  25.0, -3.0
4010,  26.5, -3.1
4011,  26.5, -4.1
4012,  25.0, -4.0
4013,  28.0, -3.2
4014,  28.0, -4.2
4015,  28.0, -3.7
4016,  25.0, -3.5
*NSET,NSET=LEFT
4001,4004,4008,4009,4012,4016
*NSET,NSET=MID
4010,4011
*NSET,NSET=RIGHT
4013,4014,4015
*NSET,NSET=TIP
4005,4006,4007,4013,4014,4015
*ELEMENT,TYPE=CPS8,ELSET=TEST11
4001, 4001,4004,4006,4005,4008,4003,4007,4002
4002, 4009,4012,4014,4013,4016,4011,4015,4010
*ELSET,ELSET=REFER
4001, 
*ELSET,ELSET=TEST
4002, 
**
** TEST 12: CPS8R
**
*NODE,NSET=NALL
5001,  25.0, -9.0
5002,  26.0, -9.0
5003,  26.0, -10.0
5004,  25.0, -10.0
5005,  27.0, -9.0
5006,  27.0, -10.0
5007,  27.0, -9.5
5008,  25.0, -9.5
5009,  25.0, -12.0
5010,  26.5, -12.1
5011,  26.5, -13.1
5012,  25.0, -13.0
5013,  28.0, -12.2
5014,  28.0, -13.2
5015,  28.0, -12.7
5016,  25.0, -12.5
*NSET,NSET=LEFT
5001,5004,5008,5009,5012,5016
*NSET,NSET=MID
5010,5011
*NSET,NSET=RIGHT
5013,5014,5015
*NSET,NSET=TIP
5005,5006,5007,5013,5014,5015
*ELEMENT,TYPE=CPS8R,ELSET=TEST12
5001, 5001,5004,5006,5005,5008,5003,5007,5002
5002, 5009,5012,5014,5013,5016,5011,5015,5010
*ELSET,ELSET=REFER
5001, 
*ELSET,ELSET=TEST
5002, 
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL
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
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*OUTPUT,HISTORY,VARIABLE=PRESELECT
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
E,S
*NODE OUTPUT
U,
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
*OUTPUT,HISTORY,VARIABLE=PRESELECT,FREQ=1
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
E,S,THE,TEMP
*NODE OUTPUT
U,
RF,
*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
*OUTPUT,FIELD,FREQ=0
*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
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF
*ELEMENT OUTPUT
S,E
*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
*OUTPUT,HISTORY,VARIABLE=PRESELECT
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT
E,S,
*NODE OUTPUT
U,
RF,
*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
*OUTPUT,FIELD,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: 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
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #12: 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
*END STEP
**
**
*STEP
 STEP #13: 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