*HEADING
C3D20 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
13,.5,0,0
14,1,.5,0
15,.5,1,0
16,0,.5,0
17,.5,0,1
18,1,.5,1
19,.5,1,1
20,0,.5,1
21,0,0,.5
22,1,0,.5
23,1,1,.5
24,0,1,.5
25,1.5,0,0
26,2,.5,0
27,1.5,1,0
28,1.5,0,1
29,2,.5,1
30,1.5,1,1
31,2,0,.5
32,2,1,.5
**
41,3,0,0
42,4,0,0
43,4,1.03,0
44,3,1.03,0
45,3,0,1
46,4,0,1
47,4,1.03,1
48,3,1.03,1
49,5,0,0
50,5,1.03,0
51,5,0,1
52,5,1.03,1
53,3.5,0,0
54,4,.515,0
55,3.5,1.03,0
56,3,.515,0
57,3.5,0,1
58,4,.515,1
59,3.5,1.03,1
60,3,.515,1
61,3,0,.5
62,4,0,.5
63,4,1.03,.5
64,3,1.03,.5
65,4.5,0,0
66,5,.515,0
67,4.5,1.03,0
68,4.5,0,1
69,5,.515,1
70,4.5,1.03,1
71,5,0,.5
72,5,1.03,.5
**
*NSET,NSET=BOT
1,2,9,5,6,11,21,22,13,17,25,28,31
41,42,49,45,46,51,61,62,53,57,65,68,71
*NSET,NSET=TIP
4,3,10,8,7,12,24,23,15,19,32,27,30
44,43,50,48,47,52,64,63,55,59,72,67,70
*NSET,NSET=DEFBC
44,43,50,48,47,52,64,63,55,59,72,67,70
*NSET,NSET=MIDSIDEW
14,16,18,20,29,26,54,56,58,60,69,66
*NSET,NSET=MIDSIDE
54,56,58,60,69,66
*ELEMENT,TYPE=C3D20,ELSET=TEST1
1,1,2,3,4,5,6,7,8,13,14,
15,16,17,18,19,20,21,22,23,24
**
2,2,9,10,3,6,11,12,7,25,26,27,14,
28,29,30,18,22,31,32,23
**
3,41,42,43,44,45,46,47,48,53,54,
55,56,57,58,59,60,61,62,63,64
**
4,42,49,50,43,46,51,52,47,65,66,67,54,
68,69,70,58,62,71,72,63
*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, 
MISESONLY,PRESSONLY
*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=9999
*ELEMENT OUTPUT
S, 
THE, 
E, 
TEMP, 
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
200, , 
*BOUNDARY,OP=NEW
BOT,1,3
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FREQ=0,FIELD
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.001,.01
*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, 
MISESONLY,PRESSONLY
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
90.,400., 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
90.,400., 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
90.,400., 5
*CLOAD
TIP,1,0.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
E, 
*END STEP
**
**
**
*STEP,NLGEOM
 STEP #13: DYNAMICS
*DYNAMIC
.001, .01
*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