*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 45: CPS3
**
*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=CPS3,ELSET=TEST45
1001, 1001,1004,1002
1002, 1002,1004,1003
1003, 1002,1003,1005
1004, 1005,1003,1006
1005, 1007,1010,1008
1006, 1008,1010,1009
1007, 1008,1009,1011
1008, 1011,1009,1012
*ELSET,ELSET=REFER
1001,1002,1003,1004
*ELSET,ELSET=TEST
1005,1006,1007,1008
**
** TEST 46: CPS6
**
*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
4017,  26.0, -0.5
4018,  26.5, -3.6
*NSET,NSET=LEFT
4001,4004,4008,4009,4012,4016
*NSET,NSET=MID
4010,4011,4018
*NSET,NSET=RIGHT
4013,4014,4015
*NSET,NSET=TIP
4005,4006,4007,4013,4014,4015
*ELEMENT,TYPE=CPS6,ELSET=TEST46
4001, 4001,4004,4005,4008,4017,4002
4002, 4005,4004,4006,4017,4003,4007
4003, 4009,4012,4013,4016,4018,4010
4004, 4013,4012,4014,4018,4011,4015
*ELSET,ELSET=REFER
4001,4002
*ELSET,ELSET=TEST
4003,4004
**
** TEST 47: CPS6M
**
*NODE,NSET=NALL
5001,  35.0,  0.0
5002,  36.0,  0.0
5003,  36.0, -1.0
5004,  35.0, -1.0
5005,  37.0,  0.0
5006,  37.0, -1.0
5007,  37.0, -0.5
5008,  35.0, -0.5
5009,  35.0, -3.0
5010,  36.5, -3.1
5011,  36.5, -4.1
5012,  35.0, -4.0
5013,  38.0, -3.2
5014,  38.0, -4.2
5015,  38.0, -3.7
5016,  35.0, -3.5
5017,  36.0, -0.5
5018,  36.5, -3.6
*NSET,NSET=LEFT
5001,5004,5008,5009,5012,5016
*NSET,NSET=MID
5010,5011,5018
*NSET,NSET=RIGHT
5013,5014,5015
*NSET,NSET=TIP
5005,5006,5007,5013,5014,5015
*ELEMENT,TYPE=CPS6M,ELSET=TEST47
5001, 5001,5004,5005,5008,5017,5002
5002, 5005,5004,5006,5017,5003,5007
5003, 5009,5012,5013,5016,5018,5010
5004, 5013,5012,5014,5018,5011,5015
*ELSET,ELSET=REFER
5001,5002
*ELSET,ELSET=TEST
5003,5004
**
**  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,FIELD,FREQ=999
*ELEMENT OUTPUT
S, 
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*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,P3,-0.25
REFER,P3,-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,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, 
*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
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
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 #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,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
*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