*HEADING
C3D TET ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
**RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 1: C3D4
**
*NODE,NSET=NALL
1,0.,0.,0.5
2,1.,0.,1.
3,1.,0.,0.
4,1.,1.,0.5
5,2.,0.,0.5
**
6,3.,0.,0.5
7,4.,0.,1.
8,4.,0.,0.
9,4.,1.03,0.5
10,5.,0.,0.5
*NSET,NSET=BOT
1,2,3,5,6,7,8,10
*NSET,NSET=TIP
4,9
*NSET,NSET=DEFBC
9, 
*ELEMENT,TYPE=C3D4,ELSET=TEST1
1,1,2,3,4
2,5,3,2,4
3,6,7,8,9
4,10,8,7,9
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
3,4
**
** TEST 2: C3D4H
**
*NODE,NSET=NALL
11,0.,2.5,0.5
12,1.,2.5,1.
13,1.,2.5,0.
14,1.,3.5,0.5
15,2.,2.5,0.5
**
16,3.,2.5,0.5
17,4.,2.5,1.
18,4.,2.5,0.
19,4.,3.53,0.5
20,5.,2.5,0.5
*NSET,NSET=BOT
11,12,13,15,16,17,18,20
*NSET,NSET=TIP
14,19
*NSET,NSET=DEFBC
19, 
*ELEMENT,TYPE=C3D4H,ELSET=TEST2
5,11,12,13,14
6,15,13,12,14
7,16,17,18,19
8,20,18,17,19
*ELSET,ELSET=REFER
5,6
*ELSET,ELSET=TEST
7,8
**
** TEST 3: C3D10
**
*NODE,NSET=NALL
101,0.,5.,0.5
102,1.,5.,1.
103,1.,5.,0.
104,1.,6.,0.5
111,.5,5.,0.75
112,1,5,0.5
113,.5,5,.25
114,.5,5.5,.5
115,1,5.5,.75
116,1,5.5,.25
105,2,5,.5
117,1.5,5,.75
118,1.5,5,.25
119,1.5,5.5,0.5
**
121,3.,5.,0.5
122,4.,5.,1.
123,4.,5.,0.
124,4.,6.03,0.5
131,3.5,5.,0.75
132,4,5,0.5
133,3.5,5,.25
134,3.5,5.515,.5
135,4,5.515,.75
136,4,5.515,.25
125,5,5,.5
137,4.5,5,.75
138,4.5,5.,.25
139,4.5,5.515,0.5
*NSET,NSET=BOT
101,102,103,105,111,112,113,117,118,121,122,123
125,131,132,133,137,138
*NSET,NSET=TIP
104,124
*NSET,NSET=MIDSIDE
134,135,136,139
*NSET,NSET=MIDSIDEW
134,135,136,139,114,115,116,119
*NSET,NSET=DEFBC
124, 
*ELEMENT,TYPE=C3D10,ELSET=TEST3
101,101,102,103,104,111,112,113,114,115,116
102,102,105,103,104,117,118,112,115,119,116
103,121,122,123,124,131,132,133,134,135,136
104,122,125,123,124,137,138,132,135,139,136
*ELSET,ELSET=REFER
101,102
*ELSET,ELSET=TEST
103,104
**
**
** TEST 4: C3D10H
**
*NODE,NSET=NALL
141,0.,7.5,0.5
142,1.,7.5,1.
143,1.,7.5,0.
144,1.,8.5,0.5
151,.5,7.5,0.75
152,1,7.5,0.5
153,.5,7.5,.25
154,.5,8,.5
155,1,8,.75
156,1,8,.25
145,2,7.5,.5
157,1.5,7.5,.75
158,1.5,7.5,.25
159,1.5,8,0.5
**
161,3.,7.5,0.5
162,4.,7.5,1.
163,4.,7.5,0.
164,4.,8.53,0.5
171,3.5,7.5,0.75
172,4,7.5,0.5
173,3.5,7.5,.25
174,3.5,8.015,.5
175,4,8.015,.75
176,4,8.015,.25
165,5,7.5,.5
177,4.5,7.5,.75
178,4.5,7.5,.25
179,4.5,8.015,0.5
*NSET,NSET=BOT
141,142,143,145,151,152,153,157,158,161,162
163,165,171,172,173,177,178
*NSET,NSET=TIP
144,164
*NSET,NSET=MIDSIDE
174,175,176,179
*NSET,NSET=MIDSIDEW
174,175,176,179,154,155,156,159
*NSET,NSET=DEFBC
164, 
*ELEMENT,TYPE=C3D10H,ELSET=TEST4
105,141,142,143,144,151,152,153,154,155,156
106,142,145,143,144,157,158,152,155,159,156
107,161,162,163,164,171,172,173,174,175,176
108,162,165,163,164,177,178,172,175,179,176
*ELSET,ELSET=REFER
105,106
*ELSET,ELSET=TEST
107,108
**
** TEST 5: C3D10M
**
*SYSTEM
0.,5.,0.
*NODE,NSET=NALL
201,0.,5.,0.5
202,1.,5.,1.
203,1.,5.,0.
204,1.,6.,0.5
211,.5,5.,0.75
212,1,5,0.5
213,.5,5,.25
214,.5,5.5,.5
215,1,5.5,.75
216,1,5.5,.25
205,2,5,.5
217,1.5,5,.75
218,1.5,5,.25
219,1.5,5.5,0.5
**
221,3.,5.,0.5
222,4.,5.,1.
223,4.,5.,0.
224,4.,6.03,0.5
231,3.5,5.,0.75
232,4,5,0.5
233,3.5,5,.25
234,3.5,5.515,.5
235,4,5.515,.75
236,4,5.515,.25
225,5,5,.5
237,4.5,5,.75
238,4.5,5.,.25
239,4.5,5.515,0.5
*NSET,NSET=BOT
201,202,203,205,211,212,213,217,218,221,222,223
225,231,232,233,237,238
*NSET,NSET=TIP
204,224
*NSET,NSET=MIDSIDE
234,235,236,239
*NSET,NSET=MIDSIDEW
234,235,236,239,214,215,216,219
*NSET,NSET=DEFBC
224, 
*ELEMENT,TYPE=C3D10M,ELSET=TEST5
201,201,202,203,204,211,212,213,214,215,216
202,202,205,203,204,217,218,212,215,219,216
203,221,222,223,224,231,232,233,234,235,236
204,222,225,223,224,237,238,232,235,239,236
*ELSET,ELSET=REFER
201,202
*ELSET,ELSET=TEST
203,204
**
**
** TEST 6: C3D10MH
**
*NODE,NSET=NALL
241,0.,7.5,0.5
242,1.,7.5,1.
243,1.,7.5,0.
244,1.,8.5,0.5
251,.5,7.5,0.75
252,1,7.5,0.5
253,.5,7.5,.25
254,.5,8,.5
255,1,8,.75
256,1,8,.25
245,2,7.5,.5
257,1.5,7.5,.75
258,1.5,7.5,.25
259,1.5,8,0.5
**
261,3.,7.5,0.5
262,4.,7.5,1.
263,4.,7.5,0.
264,4.,8.53,0.5
271,3.5,7.5,0.75
272,4,7.5,0.5
273,3.5,7.5,.25
274,3.5,8.015,.5
275,4,8.015,.75
276,4,8.015,.25
265,5,7.5,.5
277,4.5,7.5,.75
278,4.5,7.5,.25
279,4.5,8.015,0.5
*NSET,NSET=BOT
241,242,243,245,251,252,253,257,258,261,262
263,265,271,272,273,277,278
*NSET,NSET=TIP
244,264
*NSET,NSET=MIDSIDE
274,275,276,279
*NSET,NSET=MIDSIDEW
274,275,276,279,254,255,256,259
*NSET,NSET=DEFBC
264, 
*ELEMENT,TYPE=C3D10MH,ELSET=TEST6
205,241,242,243,244,251,252,253,254,255,256
206,242,245,243,244,257,258,252,255,259,256
207,261,262,263,264,271,272,273,274,275,276
208,262,265,263,264,277,278,272,275,279,276
*ELSET,ELSET=REFER
205,206
*ELSET,ELSET=TEST
207,208
**
**
**  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
E, 
S, 
*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
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
52, , 
*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
.0007,.007
*CLOAD
TIP,1,150.
*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
175.0,500.0, 5
*CLOAD
TIP,1,150.
*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
 STEP #11: FREQUENCY EXTRACTION FOR SUBSPACE PROJECTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
152, , 
*BOUNDARY,OP=NEW
BOT,1,3
*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
175.0,500.0, 5
*BOUNDARY,OP=NEW
BOT,1,3
*CLOAD
TIP,1,150.
*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 #13: STEADY-STATE DYNAMICS,SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
175.0,500.0, 5
*CLOAD
TIP,1,150.
*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 #14: DYNAMICS
*DYNAMIC
.0007, .007
*BOUNDARY,OP=NEW
BOT,PINNED
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
E, 
*NODE PRINT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TIP,1,150.
*END STEP