*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 13: CAX4
**
*NODE,NSET=NALL
6001,  35.0,  0.0
6002,  36.0,  0.0
6003,  36.0, -1.0
6004,  35.0, -1.0
6005,  37.0,  0.0
6006,  37.0, -1.0
6007,  35.0, -3.0
6008,  36.5, -3.1
6009,  36.5, -4.1
6010,  35.0, -4.0
6011,  38.0, -3.2
6012,  38.0, -4.2
*NSET,NSET=LEFT
6001,6004,6007,6010
*NSET,NSET=MID
6008,6009
*NSET,NSET=RIGHT
6011,6012
*NSET,NSET=TIP
6005,6006,6011,6012
*ELEMENT,TYPE=CAX4,ELSET=TEST13
6001, 6001,6004,6003,6002
6002, 6002,6003,6006,6005
6003, 6007,6010,6009,6008
6004, 6008,6009,6012,6011
*ELSET,ELSET=REFER
6001,6002
*ELSET,ELSET=TEST
6003,6004
**
** TEST 14: CAX4H
**
*NODE,NSET=NALL
7001,  35.0,  -9.0
7002,  36.0,  -9.0
7003,  36.0,  -10.0
7004,  35.0,  -10.0
7005,  37.0,  -9.0
7006,  37.0,  -10.0
7007,  35.0,  -12.0
7008,  36.5,  -12.1
7009,  36.5,  -13.1
7010,  35.0,  -13.0
7011,  38.0,  -12.2
7012,  38.0,  -13.2
*NSET,NSET=LEFT
7001,7004,7007,7010
*NSET,NSET=MID
7008,7009
*NSET,NSET=RIGHT
7011,7012
*NSET,NSET=TIP
7005,7006,7011,7012
*ELEMENT,TYPE=CAX4H,ELSET=TEST14
7001, 7001,7004,7003,7002
7002, 7002,7003,7006,7005
7003, 7007,7010,7009,7008
7004, 7008,7009,7012,7011
*ELSET,ELSET=REFER
7001,7002
*ELSET,ELSET=TEST
7003,7004
**
** TEST 15: CAX4I
**
*NODE,NSET=NALL
8001,  35.0,  -18.0
8002,  36.0,  -18.0
8003,  36.0,  -19.0
8004,  35.0,  -19.0
8005,  37.0,  -18.0
8006,  37.0,  -19.0
8007,  35.0,  -21.0
8008,  36.5,  -21.1
8009,  36.5,  -22.1
8010,  35.0,  -22.0
8011,  38.0,  -21.2
8012,  38.0,  -22.2
*NSET,NSET=LEFT
8001,8004,8007,8010
*NSET,NSET=MID
8008,8009
*NSET,NSET=RIGHT
8011,8012
*NSET,NSET=TIP
8005,8006,8011,8012
*ELEMENT,TYPE=CAX4I,ELSET=TEST15
8001, 8001,8004,8003,8002
8002, 8002,8003,8006,8005
8003, 8007,8010,8009,8008
8004, 8008,8009,8012,8011
*ELSET,ELSET=REFER
8001,8002
*ELSET,ELSET=TEST
8003,8004
**
** TEST 16: CAX4R
**
*NODE,NSET=NALL
9001,  35.0,  -27.0
9002,  36.0,  -27.0
9003,  36.0,  -28.0
9004,  35.0,  -28.0
9005,  37.0,  -27.0
9006,  37.0,  -28.0
9007,  35.0,  -30.0
9008,  36.5,  -30.1
9009,  36.5,  -31.1
9010,  35.0,  -31.0
9011,  38.0,  -30.2
9012,  38.0,  -31.2
*NSET,NSET=LEFT
9001,9004,9007,9010
*NSET,NSET=MID
9008,9009
*NSET,NSET=RIGHT
9011,9012
*NSET,NSET=TIP
9005,9006,9011,9012
*ELEMENT,TYPE=CAX4R,ELSET=TEST16
9001, 9001,9004,9003,9002
9002, 9002,9003,9006,9005
9003, 9007,9010,9009,9008
9004, 9008,9009,9012,9011
*ELSET,ELSET=REFER
9001,9002
*ELSET,ELSET=TEST
9003,9004
**
** TEST 17: CAX8
**
*NODE,NSET=NALL
10001,  45.0,  0.0
10002,  46.0,  0.0
10003,  46.0, -1.0
10004,  45.0, -1.0
10005,  47.0,  0.0
10006,  47.0, -1.0
10007,  47.0, -0.5
10008,  45.0, -0.5
10009,  45.0, -3.0
10010,  46.5, -3.1
10011,  46.5, -4.1
10012,  45.0, -4.0
10013,  48.0, -3.2
10014,  48.0, -4.2
10015,  48.0, -3.7
10016,  45.0, -3.5
*NSET,NSET=LEFT
10001,10004,10008,10009,10012,10016
*NSET,NSET=MID
10010,10011
*NSET,NSET=RIGHT
10013,10014,10015
*NSET,NSET=TIP
10005,10006,10007,10013,10014,10015
*ELEMENT,TYPE=CAX8,ELSET=TEST17
10001, 10001,10004,10006,10005,10008,10003,10007,10002
10002, 10009,10012,10014,10013,10016,10011,10015,10010
*ELSET,ELSET=REFER
10001, 
*ELSET,ELSET=TEST
10002, 
**
** TEST 18: CAX8H
**
*NODE,NSET=NALL
20001,  45.0,  -9.0
20002,  46.0,  -9.0
20003,  46.0,  -10.0
20004,  45.0,  -10.0
20005,  47.0,  -9.0
20006,  47.0,  -10.0
20007,  47.0,  -9.5
20008,  45.0,  -9.5
20009,  45.0,  -12.0
20010,  46.5,  -12.1
20011,  46.5,  -13.1
20012,  45.0,  -13.0
20013,  48.0,  -12.2
20014,  48.0,  -13.2
20015,  48.0,  -12.7
20016,  45.0,  -12.5
*NSET,NSET=LEFT
20001,20004,20008,20009,20012,20016
*NSET,NSET=MID
20010,20011
*NSET,NSET=RIGHT
20013,20014,20015
*NSET,NSET=TIP
20005,20006,20007,20013,20014,20015
*ELEMENT,TYPE=CAX8H,ELSET=TEST18
20001, 20001,20004,20006,20005,20008,20003,20007,20002
20002, 20009,20012,20014,20013,20016,20011,20015,20010
*ELSET,ELSET=REFER
20001, 
*ELSET,ELSET=TEST
20002, 
**
** TEST 19: CAX8R
**
*NODE,NSET=NALL
30001,  45.0,  -18.0
30002,  46.0,  -18.0
30003,  46.0,  -19.0
30004,  45.0,  -19.0
30005,  47.0,  -18.0
30006,  47.0,  -19.0
30007,  47.0,  -18.5
30008,  45.0,  -18.5
30009,  45.0,  -21.0
30010,  46.5,  -21.1
30011,  46.5,  -22.1
30012,  45.0,  -22.0
30013,  48.0,  -21.2
30014,  48.0,  -22.2
30015,  48.0,  -21.7
30016,  45.0,  -21.5
*NSET,NSET=LEFT
30001,30004,30008,30009,30012,30016
*NSET,NSET=MID
30010,30011
*NSET,NSET=RIGHT
30013,30014,30015
*NSET,NSET=TIP
30005,30006,30007,30013,30014,30015
*ELEMENT,TYPE=CAX8R,ELSET=TEST19
30001, 30001,30004,30006,30005,30008,30003,30007,30002
30002, 30009,30012,30014,30013,30016,30011,30015,30010
*ELSET,ELSET=REFER
30001, 
*ELSET,ELSET=TEST
30002, 
**
**  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
*ELEMENT OUTPUT,ELSET=TEST13
S,E
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S,E
*NODE OUTPUT
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,FIELD,FREQ=999
*ELEMENT OUTPUT
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
*OUTPUT,FIELD
*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
*ELEMENT OUTPUT
S,E
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
E, 
S, 
*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,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
*OUTPUT,FIELD,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
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
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
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*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