*HEADING
C3D27RH 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
101,.5,.5,.5
102,.5,.5,0
103,.5,.5,1
104,.5,0,.5
105,1,.5,.5
106,.5,1,.5
107,0,.5,.5
108,1.5,.5,.5
109,1.5,.5,0
110,1.5,.5,1
111,1.5,0,.5
112,2,.5,.5
113,1.5,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
201,3.5,.515,.5
202,3.5,.515,0
203,3.5,.515,1
204,3.5,0,.5
205,4,.515,.5
206,3.5,1.03,.5
207,3,.515,.5
208,4.5,.515,.5
209,4.5,.515,0
210,4.5,.515,1
211,4.5,0,.5
212,5,.515,.5
213,4.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
104,111,204,211
*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
106,113,206,213
*NSET,NSET=DEFBC
44,43,50,48,47,52,64,63,55,59,72,67,70
206,213
*NSET,NSET=MIDSIDEW
14,16,18,20,29,26,54,56,58,60,69,66
101,102,103,105,107,108,109,110,112
201,202,203,205,207,208,209,210,212
*NSET,NSET=MIDSIDE
54,56,58,60,69,66
201,202,203,205,207,208,209,210,212
*ELEMENT,TYPE=C3D27RH,ELSET=TEST1
1,1,2,3,4,5,6,7,8,13,14,
15,16,17,18,19,20,21,22,23,24,
101,102,103,104,105,106,107
**
2,2,9,10,3,6,11,12,7,25,26,27,14,
28,29,30,18,22,31,32,23,
108,109,110,111,112,113,105
**
3,41,42,43,44,45,46,47,48,53,54,
55,56,57,58,59,60,61,62,63,64,
201,202,203,204,205,206,207
**
4,42,49,50,43,46,51,52,47,65,66,67,54,
68,69,70,58,62,71,72,63,
208,209,210,211,212,213,205
*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, 
*OUTPUT,HIST,FREQ=1
*NODE OUTPUT,NSET=BOT
 RT
*OUTPUT,FIELD,FREQ=100
*NODE OUTPUT
 UT
*ELEMENT OUTPUT
 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.
*OUTPUT,FIELD
*ELEMENT OUTPUT
 GRAV
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: APPLY DLOAD TYPE BX
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST, BX, 150.
REFER, BX, 150. 
*OUTPUT,FIELD
*ELEMENT OUTPUT
 BF
*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.
*OUTPUT,FIELD
*ELEMENT OUTPUT
 HP
*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.
*OUTPUT,FIELD
*ELEMENT OUTPUT
 P,
*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. 
*OUTPUT,FIELD
*ELEMENT OUTPUT
 CENTMAG,
*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. 
*OUTPUT,FIELD
*ELEMENT OUTPUT
 CORIOMAG,
*END STEP