*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
 Coupled temperature-displacement elements
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
** TEST 30: CPEG4T
**
*NODE,NSET=NALL
40001,  55.0,  0.0
40002,  56.0,  0.0
40003,  56.0, -1.0
40004,  55.0, -1.0
40005,  57.0,  0.0
40006,  57.0, -1.0
40007,  55.0, -3.0
40008,  56.5, -3.1
40009,  56.5, -4.1
40010,  55.0, -4.0
40011,  58.0, -3.2
40012,  58.0, -4.2
*NODE,NSET=XTRA
99998,0.,0.
*NSET,NSET=LEFT
40001,40004,40007,40010
*NSET,NSET=MID
40008,40009
*NSET,NSET=RIGHT
40011,40012
*NSET,NSET=TIP
40005,40006,40011,40012
*ELEMENT,TYPE=CPEG4T,ELSET=TEST30
40001, 40001,40004,40003,40002
40002, 40002,40003,40006,40005
40003, 40007,40010,40009,40008
40004, 40008,40009,40012,40011
*ELSET,ELSET=REFER
40001,40002
*ELSET,ELSET=TEST
40003,40004
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL,REFNODE=XTRA
1.0, 
*MATERIAL,NAME=ELAS
*ELASTIC
1000., .25
*EXPANSION
1.E-3, 
*DENSITY
1.E-3, 
*CONDUCTIVITY
1.0, 
*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 
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
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
99998,3,5
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=999
E, S
*END STEP
**
**
*STEP
 STEP #2: REMOVE TEST ELEMENTS
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
1,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
LEFT,1,2
99998,3,5
*END STEP
**
**
*STEP
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
0.1,1.
*DLOAD
TEST,P4,-0.25
REFER,P4,-0.25
*END STEP
**
**
*STEP
 STEP #5: REMOVE TEST ELEMENT
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
.1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
.1,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #7: APPLY THERMAL LOAD
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
.1,1.0
*BOUNDARY,OP=NEW,FIXED
TIP,1
LEFT,1,2
*BOUNDARY,OP=NEW
NALL,11,11,100.
99998,3,5
*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
99998,3,5
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE 
*OUTPUT,FIELD,VAR=PRESELECT
*OUTPUT,HISTORY,VAR=PRESELECT
*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
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS
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
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
1000., 10000., 5
*BOUNDARY,OP=NEW
LEFT,1,2
99998,3,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
*END STEP
**
**
*STEP
 STEP #12: DYNAMICS
*DYNAMIC
.0004, .004
*BOUNDARY,OP=NEW
LEFT,1,2
99998,3,5
*EL FILE,FREQ=999
 S,E,THE,TEMP
*NODE FILE,FREQ=999
 U,
*EL PRINT,FREQ=999
S,E
*NODE PRINT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TIP,2,0.2
*OUTPUT,FIELD,VAR=PRESELECT,FREQ=1
*OUTPUT,HISTORY,VAR=PRESELECT,FREQ=1
*END STEP