*HEADING
C3D15V 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=C3D15V,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, 
*NODE FILE,FREQ=999
RF, 
*EL FILE,POSITION=CENTROID,FREQ=999
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: APPLY DLOAD TYPE BX
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, BX, 150.
REFER, BX, 150. 
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: APPLY DLOAD TYPE BY
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, BY, 150.
REFER,BY, 150.
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY DLOAD TYPE BZ
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, BZ, 150.
REFER,BZ, 150.
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: APPLY DLOAD TYPE HP1
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, HP1, 1.E3, 2., 0. 
REFER, HP1, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: APPLY DLOAD TYPE HP2
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, HP2, 1.E3, 2., 0. 
REFER, HP2, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #9: APPLY DLOAD TYPE HP3
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, HP3, 1.E3, 2., 0. 
REFER, HP3, 1.E3, 2., 0.
*END STEP
**
**
**
*STEP,NLGEOM
 STEP #10: APPLY DLOAD TYPE HP4
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, HP4, 1.E3, 2., 0. 
REFER, HP4, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #11: APPLY DLOAD TYPE P1
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, P1, 1.E3, 2., 0. 
REFER, P1, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #12: APPLY DLOAD TYPE P2
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, P2, 1.E3, 2., 0. 
REFER, P2, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #13: APPLY DLOAD TYPE P3
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, P3, 1.E3, 2., 0. 
REFER, P3, 1.E3, 2., 0.
*END STEP
**
**
**
*STEP,NLGEOM
 STEP #14: APPLY DLOAD TYPE P4
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, P4, 1.E3, 2., 0. 
REFER, P4, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #15: APPLY DLOAD TYPE CENT
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, CENT, 8.E3,0.,0.,0.,1.,0.,0. 
REFER, CENT, 8.E3,0.,0.,0.,1.,0.,0. 
*END STEP
**
**
*STEP,NLGEOM
 STEP #15: APPLY DLOAD TYPE CORIO
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, CORIO, 8.E3,0.,0.,0.,1.,0.,0. 
REFER, CORIO, 8.E3,0.,0.,0.,1.,0.,0. 
*END STEP