*HEADING
C3D PRISM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 1: C3D6
**
*NODE,NSET=NALL
1,0,0,0
2,1,0,0
3,0,1,0
4,0,0,1
5,1,0,1
6,0,1,1
7,1,1,0
8,1,1,1
**
11,3,0,0
12,4,0,0
13,3,1.03,0
14,3,0,1
15,4,0,1
16,3,1.03,1
17,4,1.03,0
18,4,1.03,1
*NSET,NSET=BOT
1,2,4,5,11,12,14,15
*NSET,NSET=TIP
6,3,8,7,16,13,18,17
*NSET,NSET=DEFBC
16,13,18,17
*ELEMENT,TYPE=C3D6,ELSET=TEST1
1,1,2,3,4,5,6
2,3,2,7,6,5,8
3,11,12,13,14,15,16
4,13,12,17,16,15,18
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
3,4
**
** TEST 2: C3D6H
**
*NODE,NSET=NALL
21,0,2,0
22,1,2,0
23,0,3,0
24,0,2,1
25,1,2,1
26,0,3,1
27,1,3,0
28,1,3,1
**
31,3,2,0
32,4,2,0
33,3,3.03,0
34,3,2,1
35,4,2,1
36,3,3.03,1
37,4,3.03,0
38,4,3.03,1
*NSET,NSET=BOT
21,22,24,25,31,32,34,35
*NSET,NSET=TIP
26,23,28,27,36,33,38,37
*NSET,NSET=DEFBC
36,33,38,37
*ELEMENT,TYPE=C3D6H,ELSET=TEST1
5,21,22,23,24,25,26
6,23,22,27,26,25,28
7,31,32,33,34,35,36
8,33,32,37,36,35,38
*ELSET,ELSET=REFER
5,6
*ELSET,ELSET=TEST
7,8
**
**
** TEST 3: C3D15
**
*NODE,NSET=NALL
101,0,4,0
102,1,4,0
103,0,5,0
104,0,4,1
105,1,4,1
106,0,5,1
107,.5,4,0
108,.5,4.5,0
109,0,4.5,0
110,.5,04,1
111,.5,4.5,1
112,0,4.5,1
113,0,04,.5
114,1,04,.5
115,0,5,.5
116,1,5,1
117,1,5,0
118,.5,5,0
119,.5,5,1
120,1,5,.5
121,1,4.5,1
122,1,4.5,0
**
131,3,4,0
132,4,4,0
133,3,5.03,0
134,3,4,1
135,4,4,1
136,3,5.03,1
137,3.5,4,0
138,3.5,4.515,0
139,3,4.515,0
140,3.5,4,1
141,3.5,4.515,1
142,03,4.515,1
143,03,04,.5
144,4,04,.5
145,03,5.03,.5
146,4,5.03,1
147,4,5.03,0
148,3.5,5.03,0
149,3.5,5.03,1
150,4,5.03,.5
151,4,4.515,1
152,4,4.515,0
*NSET,NSET=BOT
101,102,104,105,107,110,113,114
131,132,134,135,137,140,143,144
*NSET,NSET=TIP
103,106,117,116,115,118,119,120
133,136,147,146,145,148,149,150
*NSET,NSET=MIDSIDE
142,139,141,138,151,152
*NSET,NSET=MIDSIDEW
112,109,111,108,121,122
142,139,141,138,151,152
*NSET,NSET=DEFBC
133,136,147,146,145,148,149,150
*ELEMENT,TYPE=C3D15,ELSET=TEST3
101,101,102,103,104,105,106,107,108,
109,110,111,112,113,114,115
**
102,103,102,117,106,105,116,108,122,
118,111,121,119,115,114,120
**
103,131,132,133,134,135,136,137,138,
139,140,141,142,143,144,145
**
104,133,132,147,136,135,146,138,152,
148,141,151,149,145,144,150
*ELSET,ELSET=REFER
101,102
*ELSET,ELSET=TEST
103,104
**
**
** TEST 4: C3D15H
**
*NODE,NSET=NALL
201,0,6,0
202,1,6,0
203,0,7,0
204,0,6,1
205,1,6,1
206,0,7,1
207,.5,6,0
208,.5,6.5,0
209,0,6.5,0
210,.5,06,1
211,.5,6.5,1
212,0,6.5,1
213,0,06,.5
214,1,06,.5
215,0,7,.5
216,1,7,1
217,1,7,0
218,.5,7,0
219,.5,7,1
220,1,7,.5
221,1,6.5,1
222,1,6.5,0
**
231,3,6,0
232,4,6,0
233,3,7.03,0
234,3,6,1
235,4,6,1
236,3,7.03,1
237,3.5,6,0
238,3.5,6.515,0
239,3,6.515,0
240,3.5,6,1.0
241,3.5,6.515,1
242,03,6.515,1
243,03,06,.5
244,4,06,.5
245,03,7.03,.5
246,4,7.03,1
247,4,7.03,0
248,3.5,7.03,0
249,3.5,7.03,1
250,4,7.03,.5
251,4,6.515,1
252,4,6.515,0
*NSET,NSET=BOT
201,202,204,205,207,210,213,214
231,232,234,235,237,240,243,244
*NSET,NSET=TIP
203,206,217,216,215,218,219,220
233,236,247,246,245,248,249,250
*NSET,NSET=MIDSIDE
242,239,241,238,251,252
*NSET,NSET=MIDSIDEW
212,209,211,208,221,222
242,239,241,238,251,252
*NSET,NSET=DEFBC
233,236,247,246,245,248,249,250
*ELEMENT,TYPE=C3D15H,ELSET=TEST4
201,201,202,203,204,205,206,207,208,
209,210,211,212,213,214,215
**
202,203,202,217,206,205,216,208,222,
218,211,221,219,215,214,220
**
203,231,232,233,234,235,236,237,238,
239,240,241,242,243,244,245
**
204,233,232,247,236,235,246,238,252,
248,241,251,249,245,244,250
*ELSET,ELSET=REFER
201,202
*ELSET,ELSET=TEST
203,204
**
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
**
*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
*NODE OUTPUT
RF, 
*ELEMENT OUTPUT
S, 
*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=999
*ELEMENT OUTPUT
S, 
THE, 
E, 
TEMP, 
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
100, , 
*BOUNDARY,OP=NEW
BOT,1,3
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.002,.02
*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=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
50.0,150.0, 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
50.0,150.0, 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 #11: STEADY-STATE DYNAMICS, SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
50.0,150.0, 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,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