*HEADING
C3D8RH ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
**
*NODE,NSET=NALL
1,0,0,0
2,1,0,0
3,1,1,0
4,0,1,0
5,0,0,1
6,1,0,1
7,1,1,1
8,0,1,1
9,2,0,0
10,2,1,0
11,2,0,1
12,2,1,1
**
21,3,0,0
22,4,0,0
23,4,1.03,0
24,3,1.03,0
25,3,0,1
26,4,0,1
27,4,1.03,1
28,3,1.03,1
29,5,0,0
30,5,1.03,0
31,5,0,1
32,5,1.03,1
*NSET,NSET=BOT
1,2,9,5,6,11,21,22,29,25,26,31
*NSET,NSET=TIP
4,3,10,8,7,12,24,23,30,28,27,32
*NSET,NSET=DEFBC
24,23,30,28,27,32
*ELEMENT,TYPE=C3D8RH,ELSET=TEST1
1,1,2,3,4,5,6,7,8
2,2,9,10,3,6,11,12,7
3,21,22,23,24,25,26,27,28
4,22,29,30,23,26,31,32,27
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
3,4
**
**
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL
*MATERIAL,NAME=ELAS
*ELASTIC
1.E10, .25
*EXPANSION
1.E-3, 
*DENSITY
8.E3, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.,0.
*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
BOT,PINNED
DEFBC,2,2,-0.03
DEFBC,1,1,0.0
DEFBC,3,3,0.0
**MIDSIDE,2,2,-0.015
**MIDSIDE,1,1,0.0
**MIDSIDE,3,3,0.0
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
E, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S, 
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*ELEMENT OUTPUT,POSITION=CENTROID
E, 
S, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF, 
*END STEP
**
**
*STEP,NLGEOM
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
BOT,PINNED
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.25,1.
*DLOAD
TEST, GRAV,150.,0.,1.,0.
REFER,GRAV,150.,0.,1.,0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.25,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.25,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY THERMAL LOAD
*STATIC
.25,1.0
*BOUNDARY,OP=NEW,FIXED
BOT,1,3
TIP,1,3
**MIDSIDEW,1,3
*TEMPERATURE
NALL,100.
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
THE, 
E, 
TEMP, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S, 
THE, 
E, 
TEMP, 
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos
50, , ,-1000.
*BOUNDARY,OP=NEW
BOT,1,3
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.005,.05
*CLOAD
TIP,1,.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,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, 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
10.,350., 5
*CLOAD
TIP,1,0.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,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 #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
10.,350., 5
*BOUNDARY,OP=NEW
BOT,1,3
*CLOAD
TIP,1,0.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,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 #12: STEADY-STATE DYNAMICS, SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
10.,350., 5
*CLOAD
TIP,1,0.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,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,NLGEOM
 STEP #13: DYNAMICS
*DYNAMIC
.005, .05
*BOUNDARY,OP=NEW
BOT,PINNED
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
E, 
*NODE PRINT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TIP,1,.1
*END STEP