*HEADING
C3D20R 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=C3D20R,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, 
*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