*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
 Axisymmetric thermally coupled twist elements
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 33: CGAX3T
**
*NODE,NSET=NALL
1,  5.0,  0.0
2,  6.0,  0.0
3,  6.0, -1.0
4,  5.0, -1.0
5,  7.0,  0.0
6,  7.0, -1.0
7,  5.0, -3.0
8,  6.5, -3.1
9,  6.5, -4.1
10, 5.0, -4.0
11, 8.0, -3.2
12, 8.0, -4.2
*NSET,NSET=LEFT
1,4,7,10
*NSET,NSET=MID
8,9
*NSET,NSET=RIGHT
11,12
*NSET,NSET=TIP
5,6,11,12
*ELEMENT,TYPE=CGAX3T,ELSET=TEST33
1, 1,4,2
2, 2,4,3
3, 2,3,5
4, 3,6,5
5, 7,10,8
6, 10,9,8
7, 8,9,11
8, 9,12,11
*ELSET,ELSET=REFER
1,2,3,4
*ELSET,ELSET=TEST
5,6,7,8
*ELSET,ELSET=REFER4
1,3
*ELSET,ELSET=TEST4
5,7
**
** TEST 34: CGAX3HT
**
*NODE,NSET=NALL
101, 5.0, -9.0
102, 6.0, -9.0
103, 6.0, -10.0
104, 5.0, -10.0
105, 7.0, -9.0
106, 7.0, -10.0
107, 5.0, -12.0
108, 6.5, -12.1
109, 6.5, -13.1
110, 5.0, -13.0
111, 8.0, -12.2
112, 8.0, -13.2
*NSET,NSET=LEFT
101,104,107,110
*NSET,NSET=MID
108,109
*NSET,NSET=RIGHT
111,112
*NSET,NSET=TIP
105,106,111,112
*ELEMENT,TYPE=CGAX3HT,ELSET=TEST34
101, 101,104,102
102, 104,103,102
103, 102,103,105
104, 103,106,105
105, 107,110,108
106, 110,109,108
107, 108,109,111
108, 109,112,111
*ELSET,ELSET=REFER
101,102,103,104
*ELSET,ELSET=TEST
105,106,107,108
*ELSET,ELSET=REFER4
101,103
*ELSET,ELSET=TEST4
105,107
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL
1.0, 
*MATERIAL,NAME=ELAS
*ELASTIC
1000., .25
*EXPANSION
1.E-3, 
*DENSITY
1.E-3, 
*CONDUCTIVITY
1.0, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.
*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
LEFT,1,2
NALL,5,5
RIGHT,1,1,-1.0
RIGHT,2,2, 0.2
MID,1,1,-0.5
MID,2,2, 0.1
*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
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
LEFT,1,2
LEFT,5,5
*END STEP
**
**
*STEP
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.1,1.
*DLOAD
TEST4,P3,-0.25
REFER4,P3,-0.25
*END STEP
**
**
*STEP
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.1,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY DLOAD TYPE BR & BZ
*STATIC
.25,1.0
*EL PRINT,POSITION=CENTROID,FREQ=999
S,
E,
*EL FILE,POSITION=CENTROID,FREQ=999
S,E
*DLOAD, OP=NEW
TEST,  BR, 150.
REFER, BR, 150. 
TEST, BZ, 150.
REFER,BZ, 150.
*END STEP