*HEADING
C3D15VH PRISM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
*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
201,0.5,4,0.5
202,0.5,4.5,0.5
203,0,4.5,0.5
204,1,4.5,0.5
205,.5,5,0.5
**
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
301,3.5,4,0.5
302,3.5,4.515,0.5
303,3,4.515,0.5
304,4,4.515,0.5
305,3.5,5.03,0.5
*NSET,NSET=BOT
101,102,104,105,107,110,113,114
131,132,134,135,137,140,143,144
201,301
*NSET,NSET=TIP
103,106,117,116,115,118,119,120
133,136,147,146,145,148,149,150
205,305
*NSET,NSET=MIDSIDE
142,139,141,138,151,152
302,303,304
*NSET,NSET=MIDSIDEW
112,109,111,108,121,122
142,139,141,138,151,152
202,203,204
302,303,304
*NSET,NSET=DEFBC
133,136,147,146,145,148,149,150
305, 
*ELEMENT,TYPE=C3D15VH,ELSET=TEST3
101,101,102,103,104,105,106,107,108,
109,110,111,112,113,114,115,
201,202,203
**
102,103,102,117,106,105,116,108,122,
118,111,121,119,115,114,120,
202,204,205
**
103,131,132,133,134,135,136,137,138,
139,140,141,142,143,144,145,
301,302,303
**
104,133,132,147,136,135,146,138,152,
148,141,151,149,145,144,150,
302,304,305
*ELSET,ELSET=REFER
101,102
*ELSET,ELSET=TEST
103,104
**
**
*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,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
E, S,
THE,
*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
.001,.01
*CLOAD
TIP,1,-0.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
E, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
S,E
*NODE OUTPUT
RF, 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
50.0,600.0, 5
*CLOAD
TIP,1,-0.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
E, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
S,E
*NODE OUTPUT
RF, 
*END STEP
**
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
50.0,600.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
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
S,E
*NODE OUTPUT
RF, 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #12: STEADY-STATE DYNAMICS, SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
50.0,600.0, 5
*CLOAD
TIP,1,-0.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
E, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
S,E
*NODE OUTPUT
RF, 
*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