*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
 Axisymmetric twist elements 
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 33: CGAX4T
**
*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=CGAX4T,ELSET=TEST33
1, 1,4,3,2
2, 2,3,6,5
3, 7,10,9,8
4, 8,9,12,11
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
3,4
**
** TEST 34: CGAX4HT
**
*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=CGAX4HT,ELSET=TEST34
101, 101,104,103,102
102, 102,103,106,105
103, 107,110,109,108
104, 108,109,112,111
*ELSET,ELSET=REFER
101,102
*ELSET,ELSET=TEST
103,104
**
** TEST 35: CGAX8T
**
*NODE,NSET=NALL
201,  15.0,  0.0
202,  16.0,  0.0
203,  16.0, -1.0
204,  15.0, -1.0
205,  17.0,  0.0
206,  17.0, -1.0
207,  17.0, -0.5
208,  15.0, -0.5
209,  15.0, -3.0
210,  16.5, -3.1
211,  16.5, -4.1
212,  15.0, -4.0
213,  18.0, -3.2
214,  18.0, -4.2
215,  18.0, -3.7
216,  15.0, -3.5
*NSET,NSET=LEFT
201,204,208,209,212,216
*NSET,NSET=MID
210,211
*NSET,NSET=RIGHT
213,214,215
*NSET,NSET=TIP
205,206,207,213,214,215
*ELEMENT,TYPE=CGAX8T,ELSET=TEST35
201, 201,204,206,205,208,203,207,202
202, 209,212,214,213,216,211,215,210
*ELSET,ELSET=REFER
201, 
*ELSET,ELSET=TEST
202, 
**
**
** TEST 36: CGAX8HT
**
*NODE,NSET=NALL
401,  15.0,  0.0
402,  16.0,  0.0
403,  16.0, -1.0
404,  15.0, -1.0
405,  17.0,  0.0
406,  17.0, -1.0
407,  17.0, -0.5
408,  15.0, -0.5
409,  15.0, -3.0
410,  16.5, -3.1
411,  16.5, -4.1
412,  15.0, -4.0
413,  18.0, -3.2
414,  18.0, -4.2
415,  18.0, -3.7
416,  15.0, -3.5
*NSET,NSET=LEFT
401,404,408,409,412,416
*NSET,NSET=MID
410,411
*NSET,NSET=RIGHT
413,414,415
*NSET,NSET=TIP
405,406,407,413,414,415
*ELEMENT,TYPE=CGAX8HT,ELSET=TEST36
401, 401,404,406,405,408,403,407,402
402, 409,412,414,413,416,411,415,410
*ELSET,ELSET=REFER
401, 
*ELSET,ELSET=TEST
402, 
**
** TEST 37: CGAX8RHT
**
*NODE,NSET=NALL
501,  15.0, -9.0
502,  16.0, -9.0
503,  16.0, -10.0
504,  15.0, -10.0
505,  17.0, -9.0
506,  17.0, -10.0
507,  17.0, -9.5
508,  15.0, -9.5
509,  15.0, -12.0
510,  16.5, -12.1
511,  16.5, -13.1
512,  15.0, -13.0
513,  18.0, -12.2
514,  18.0, -13.2
515,  18.0, -12.7
516,  15.0, -12.5
*NSET,NSET=LEFT
501,504,508,509,512,516
*NSET,NSET=MID
510,511
*NSET,NSET=RIGHT
513,514,515
*NSET,NSET=TIP
505,506,507,513,514,515
*ELEMENT,TYPE=CGAX8RHT,ELSET=TEST37
501, 501,504,506,505,508,503,507,502
502, 509,512,514,513,516,511,515,510
*ELSET,ELSET=REFER
501, 
*ELSET,ELSET=TEST
502, 
**
** TEST 38: CGAX8RT
**
*NODE,NSET=NALL
601,  15.0, -18.0
602,  16.0, -18.0
603,  16.0, -19.0
604,  15.0, -19.0
605,  17.0, -18.0
606,  17.0, -19.0
607,  17.0, -18.5
608,  15.0, -18.5
609,  15.0, -21.0
610,  16.5, -21.1
611,  16.5, -22.1
612,  15.0, -22.0
613,  18.0, -21.2
614,  18.0, -22.2
615,  18.0, -21.7
616,  15.0, -21.5
*NSET,NSET=LEFT
601,604,608,609,612,616
*NSET,NSET=MID
610,611
*NSET,NSET=RIGHT
613,614,615
*NSET,NSET=TIP
605,606,607,613,614,615
*ELEMENT,TYPE=CGAX8RT,ELSET=TEST38
601, 601,604,606,605,608,603,607,602
602, 609,612,614,613,616,611,615,610
*ELSET,ELSET=REFER
601, 
*ELSET,ELSET=TEST
602, 
** TEST 38: CGAX4RT
**
*NODE,NSET=NALL
701,  5.0,  -30.0
702,  6.0,  -30.0
703,  6.0, -31.0
704,  5.0, -31.0
705,  7.0,  -30.0
706,  7.0, -31.0
707,  5.0, -33.0
708,  6.5, -33.1
709,  6.5, -34.1
710, 5.0, -34.0
711, 8.0, -33.2
712, 8.0, -34.2
*NSET,NSET=LEFT
701,704,707,710
*NSET,NSET=MID
708,709
*NSET,NSET=RIGHT
711,712
*NSET,NSET=TIP
705,706,711,712
*ELEMENT,TYPE=CGAX4RT,ELSET=TEST38
701, 701,704,703,702
702, 702,703,706,705
703, 707,710,709,708
704, 708,709,712,711
*ELSET,ELSET=REFER
701,702
*ELSET,ELSET=TEST
703,704
**
**
** TEST 39: CGAX4RHT
**
*NODE,NSET=NALL
801,  15.0, -30.0
802,  16.0, -30.0
803,  16.0, -31.0
804,  15.0, -31.0
805,  17.0, -30.0
806,  17.0, -31.0
807,  15.0, -33.0
808,  16.5, -33.1
809,  16.5, -34.1
810, 15.0, -34.0
811, 18.0, -33.2
812, 18.0, -34.2
*NSET,NSET=LEFT
801,804,807,810
*NSET,NSET=MID
808,809
*NSET,NSET=RIGHT
811,812
*NSET,NSET=TIP
805,806,811,812
*ELEMENT,TYPE=CGAX4RHT,ELSET=TEST39
801, 801,804,803,802
802, 802,803,806,805
803, 807,810,809,808
804, 808,809,812,811
*ELSET,ELSET=REFER
801,802
*ELSET,ELSET=TEST
803,804
**
**
**  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, 
LE, 
*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
LEFT,1,2
LEFT,5,5
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.1,1.
*DLOAD
TEST,P4,-0.25
REFER,P4,-0.25
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP, NLGEOM
 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,
LE,
*DLOAD, OP=NEW
TEST,  BR, 150.
REFER, BR, 150. 
TEST, BZ, 150.
REFER,BZ, 150.
*END STEP