*HEADING
 SHELL ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*PREPRINT,MODEL=YES,HISTORY=NO
*RESTART,WRITE
*************************
** TEST 1: S4R
*************************
*NODE,NSET=NALL
101,  0.0, 1., 0.0
102,  1.0, 1., 0.0
103,  2.0, 1., 0.0
104,  0.0, 1., 0.5
105,  1.0, 1., 0.5
106,  2.0, 1., 0.5
107,  0.0, 1., 1.0
108,  1.0, 1., 1.0
109,  2.0, 1., 1.0
111, 0.0,  0., 0.0
112, 1.5, -.1, 0.0
113, 3.0, -.3, 0.0
114, 0.0,  .0, 0.5
115, 1.5, -.1, 0.5
116, 3.0, -.3, 0.5
117, 0.0,  .0, 1.0
118, 1.5, -.1, 1.0
119, 3.0, -.3, 1.0
*NSET,NSET=BOT
101,102,103,111,112,113
*NSET,NSET=LEFT
101,111,104,114,107,117
*NSET,NSET=MID
112,115,118
*NSET,NSET=RIGHT
113,116,119
*NSET,NSET=TIP
103,106,109,113,116,119
*ELEMENT,TYPE=S4R,ELSET=TEST1
101, 101,102,105,104
102, 102,103,106,105
103, 104,105,108,107
104, 105,106,109,108
111,111,112,115,114
112,112,113,116,115
113,114,115,118,117
114,115,116,119,118
*ELSET,ELSET=REFER
101,102,103,104
*ELSET,ELSET=TEST
111,112,113,114
*ELSET,ELSET=EOUT
101,111
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST1
.1, 
*************************
** TEST 2: S8R
*************************
*NODE,NSET=NALL
201, 5.0, 1.0, 0.0
202, 6.0, 1.0, 0.0
203, 7.0, 1.0, 0.0
204, 5.0, 1.0, 0.5
205, 6.0, 1.0, 0.5
206, 7.0, 1.0, 0.5 
207, 5.0, 1.0, 1.0
208, 6.0, 1.0, 1.0
209, 7.0, 1.0, 1.0
211, 5.0,  .0, 0.0
212, 6.5, -.1, 0.0
213, 8.0, -.3, 0.0
214, 5.0,  .0, 0.5
215, 6.5, -.1, 0.5
216, 8.0, -.3, 0.5 
217, 5.0,  .0, 1.0
218, 6.5, -.1, 1.0
219, 8.0, -.3, 1.0
*NSET,NSET=BOT
201,202,203,211,212,213
*NSET,NSET=LEFT
201,204,207,211,214,217
*NSET,NSET=MID
212,218
*NSET,NSET=RIGHT
213,216,219
*NSET,NSET=TIP
203,206,209,213,216,219
*ELEMENT,TYPE=S8R,ELSET=TEST2
201,201,203,209,207,202,206,208,204
211,211,213,219,217,212,216,218,214
*ELSET,ELSET=REFER
201, 
*ELSET,ELSET=TEST
211, 
*ELSET,ELSET=EOUT
201,211
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST2
.1, 
********************************
** TEST 3: S4R WITH ORIENTATION
********************************
*NODE,NSET=NALL
301, 10.0, 1.0, 0.0
302, 11.0, 1.0, 0.0
303, 12.0, 1.0, 0.0
304, 10.0, 1.0, 0.5
305, 11.0, 1.0, 0.5
306, 12.0, 1.0, 0.5 
307, 10.0, 1.0, 1.0
308, 11.0, 1.0, 1.0
309, 12.0, 1.0, 1.0
311, 10.0,  .0, 0.0
312, 11.5, -.1, 0.0
313, 13.0, -.3, 0.0
314, 10.0,  .0, 0.5
315, 11.5, -.1, 0.5
316, 13.0, -.3, 0.5 
317, 10.0,  .0, 1.0
318, 11.5, -.1, 1.0
319, 13.0, -.3, 1.0
*NSET,NSET=BOT
301,302,303,311,312,313
*NSET,NSET=LEFT
301,311,304,314,307,317
*NSET,NSET=MID
312,315,318
*NSET,NSET=RIGHT
313,316,319
*NSET,NSET=TIP
303,313,306,316,309,319
*ELEMENT,TYPE=S4R,ELSET=TEST3
301, 301,302,305,304
302, 302,303,306,305
303, 304,305,308,307
304, 305,306,309,308
311,311,312,315,314
312,312,313,316,315
313,314,315,318,317
314,315,316,319,318
*ELSET,ELSET=REFER
301,302,303,304
*ELSET,ELSET=TEST
311,312,313,314
*ELSET,ELSET=EOUT
301,311
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST3
.1, 
*ORIENTATION,NAME=ORIENT3,SYSTEM=C
0.,1.,0.,  0.,2.,0.
3,0.
**************************
** TEST 4: S4R WITH REBAR
**************************
*NODE,NSET=NALL
401, 0.0, 4.0, 0.0
402, 1.0, 4.0, 0.0
403, 2.0, 4.0, 0.0
404, 0.0, 4.0, 0.5
405, 1.0, 4.0, 0.5
406, 2.0, 4.0, 0.5
407, 0.0, 4.0, 1.0
408, 1.0, 4.0, 1.0
409, 2.0, 4.0, 1.0
411, 0.0, 3.0, 0.0
412, 1.5, 2.9, 0.0
413, 3.0, 2.7, 0.0
414, 0.0, 3.0, 0.5
415, 1.5, 2.9, 0.5
416, 3.0, 2.7, 0.5
417, 0.0, 3.0, 1.0
418, 1.5, 2.9, 1.0
419, 3.0, 2.7, 1.0
*NSET,NSET=BOT
401,402,403,411,412,413
*NSET,NSET=LEFT
401,411,404,414,407,417
*NSET,NSET=MID
412,415,418
*NSET,NSET=RIGHT
413,416,419
*NSET,NSET=TIP
403,413,406,416,409,419
*ELEMENT,TYPE=S4R,ELSET=TEST4
401, 401,402,405,404
402, 402,403,406,405
403, 404,405,408,407
404, 405,406,409,408
411,411,412,415,414
412,412,413,416,415
413,414,415,418,417
414,415,416,419,418
*ELSET,ELSET=REFER
401,402,403,404
*ELSET,ELSET=TEST
411,412,413,414
*ELSET,ELSET=EOUT
401,411
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST4
.1,
*REBAR LAYER
REBAR4,1.0,0.05,0.05,ELAS2,0,1
REBAR4_2,1.0,0.05,-0.05,ELAS2,0,1  
*****************************
** TEST 5: S8R WITH REBAR
*****************************
*NODE,NSET=NALL
501, 5.0,  4.0, 0.0
502, 6.0,  4.0, 0.0
503, 7.0,  4.0, 0.0
504, 5.0,  4.0, 0.5
505, 6.0,  4.0, 0.5
506, 7.0,  4.0, 0.5 
507, 5.0,  4.0, 1.0
508, 6.0,  4.0, 1.0
509, 7.0,  4.0, 1.0
511, 5.0,  3.0, 0.0
512, 6.5,  2.9, 0.0
513, 8.0,  2.7, 0.0
514, 5.0,  3.0, 0.5
515, 6.5,  2.9, 0.5
516, 8.0,  2.7, 0.5 
517, 5.0,  3.0, 1.0
518, 6.5,  2.9, 1.0
519, 8.0,  2.7, 1.0
*NSET,NSET=BOT
501,502,503,511,512,513
*NSET,NSET=LEFT
501,504,507,511,514,517
*NSET,NSET=MID
512,518
*NSET,NSET=RIGHT
513,516,519
*NSET,NSET=TIP
503,506,509,513,516,519
*ELEMENT,TYPE=S8R,ELSET=TEST5
501,501,503,509,507,502,506,508,504
511,511,513,519,517,512,516,518,514
*ELSET,ELSET=REFER
501, 
*ELSET,ELSET=TEST
511, 
*ELSET,ELSET=EOUT
501,511
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST5
.1,
*REBAR LAYER
REBAR5,1.0,0.05,0.05,ELAS2,0,1
REBAR5_2,1.0,0.05,-0.05,ELAS2,0,1
*************************************
** TEST 6: S8R WITH NODAL THICKNESS
*************************************
*NODE,NSET=NALL
601, 10.0,  4.0, 0.0
602, 11.0,  4.0, 0.0
603, 12.0,  4.0, 0.0
604, 10.0,  4.0, 0.5
605, 11.0,  4.0, 0.5
606, 12.0,  4.0, 0.5 
607, 10.0,  4.0, 1.0
608, 11.0,  4.0, 1.0
609, 12.0,  4.0, 1.0
611, 10.0,  3.0, 0.0
612, 11.5,  2.9, 0.0
613, 13.0,  2.7, 0.0
614, 10.0,  3.0, 0.5
615, 11.5,  2.9, 0.5
616, 13.0,  2.7, 0.5 
617, 10.0,  3.0, 1.0
618, 11.5,  2.9, 1.0
619, 13.0,  2.7, 1.0
*NSET,NSET=BOT
601,602,603,611,612,613
*NSET,NSET=LEFT
601,604,607,611,614,617
*NSET,NSET=MID
612,618
*NSET,NSET=RIGHT
613,616,619
*NSET,NSET=TIP
603,606,609,613,616,619
*NSET,NSET=L6
601,604,607,611,614,617
*NSET,NSET=M6
602,608,612,618
*NSET,NSET=R6
603,606,609,613,616,619
*NODAL THICKNESS
L6,.05
M6,.1
R6,.15
*ELEMENT,TYPE=S8R,ELSET=TEST6
601,601,603,609,607,602,606,608,604
611,611,613,619,617,612,616,618,614
*ELSET,ELSET=REFER
601, 
*ELSET,ELSET=TEST
611, 
*ELSET,ELSET=EOUT
601,611
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST6,NODAL THICKNESS
**************************
** TEST 7: S3R
**************************
*NODE,NSET=NALL
701, 0.0, 7.0, 0.0
702, 1.0, 7.0, 0.0
703, 2.0, 7.0, 0.0
704, 0.0, 7.0, 0.5
705, 1.0, 7.0, 0.5
706, 2.0, 7.0, 0.5
707, 0.0, 7.0, 1.0
708, 1.0, 7.0, 1.0
709, 2.0, 7.0, 1.0
711, 0.0, 6.0, 0.0
712, 1.5, 5.9, 0.0
713, 3.0, 5.7, 0.0
714, 0.0, 6.0, 0.5
715, 1.5, 5.9, 0.5
716, 3.0, 5.7, 0.5
717, 0.0, 6.0, 1.0
718, 1.5, 5.9, 1.0
719, 3.0, 5.7, 1.0
*NSET,NSET=BOT
701,702,703,711,712,713
*NSET,NSET=LEFT
701,711,704,714,707,717
*NSET,NSET=MID
712,715,718
*NSET,NSET=RIGHT
713,716,719
*NSET,NSET=TIP
703,713,706,716,709,719
*ELEMENT,TYPE=S3R,ELSET=TEST7
701, 701,702,705
702, 701,705,704
703, 702,703,706
704, 702,706,705
705, 704,705,708
706, 704,708,707
707, 705,706,709
708, 705,709,708
711, 711,712,715
712, 711,715,714
713, 712,713,716
714, 712,716,715
715, 714,715,718
716, 714,718,717
717, 715,716,719
718, 715,719,718
*ELSET,ELSET=REFER
701,702,703,704,705,706,707,708
*ELSET,ELSET=TEST
711,712,713,714,715,716,717,718
*ELSET,ELSET=EOUT
701,711
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST7
.1, 
**************************
** TEST 8: STRI3
**************************
*NODE,NSET=NALL
801, 5.0, 7.0, 0.0
802, 6.0, 7.0, 0.0
803, 7.0, 7.0, 0.0
804, 5.0, 7.0, 0.5
805, 6.0, 7.0, 0.5
806, 7.0, 7.0, 0.5
807, 5.0, 7.0, 1.0
808, 6.0, 7.0, 1.0
809, 7.0, 7.0, 1.0
811, 5.0, 6.0, 0.0
812, 6.5, 5.9, 0.0
813, 8.0, 5.7, 0.0
814, 5.0, 6.0, 0.5
815, 6.5, 5.9, 0.5
816, 8.0, 5.7, 0.5
817, 5.0, 6.0, 1.0
818, 6.5, 5.9, 1.0
819, 8.0, 5.7, 1.0
*NSET,NSET=BOT
801,802,803,811,812,813
*NSET,NSET=LEFT
801,811,804,814,807,817
*NSET,NSET=MID
812,815,818
*NSET,NSET=RIGHT
813,816,819
*NSET,NSET=TIP
803,813,806,816,809,819
*ELEMENT,TYPE=S3R,ELSET=TEST7
801, 801,802,805
802, 801,805,804
803, 802,803,806
804, 802,806,805
805, 804,805,808
806, 804,808,807
807, 805,806,809
808, 805,809,808
811, 811,812,815
812, 811,815,814
813, 812,813,816
814, 812,816,815
815, 814,815,818
816, 814,818,817
817, 815,816,819
818, 815,819,818
*ELSET,ELSET=REFER
801,802,803,804,805,806,807,808
*ELSET,ELSET=TEST
811,812,813,814,815,816,817,818
*ELSET,ELSET=EOUT
801,811
*************************
** TEST 9: S4
*************************
*NODE,NSET=NALL
901,  0.0, 1., 0.0
902,  1.0, 1., 0.0
903,  2.0, 1., 0.0
904,  0.0, 1., 0.5
905,  1.0, 1., 0.5
906,  2.0, 1., 0.5
907,  0.0, 1., 1.0
908,  1.0, 1., 1.0
909,  2.0, 1., 1.0
911, 0.0,  0., 0.0
912, 1.5, -.1, 0.0
913, 3.0, -.3, 0.0
914, 0.0,  .0, 0.5
915, 1.5, -.1, 0.5
916, 3.0, -.3, 0.5
917, 0.0,  .0, 1.0
918, 1.5, -.1, 1.0
919, 3.0, -.3, 1.0
*NSET,NSET=BOT
901,902,903,911,912,913
*NSET,NSET=LEFT
901,911,904,914,907,917
*NSET,NSET=MID
912,915,918
*NSET,NSET=RIGHT
913,916,919
*NSET,NSET=TIP
903,906,909,913,916,919
*ELEMENT,TYPE=S4,ELSET=TEST9
901, 901,902,905,904
902, 902,903,906,905
903, 904,905,908,907
904, 905,906,909,908
911,911,912,915,914
912,912,913,916,915
913,914,915,918,917
914,915,916,919,918
*ELSET,ELSET=REFER
901,902,903,904
*ELSET,ELSET=TEST
911,912,913,914
*ELSET,ELSET=EOUT
901,911
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST9
.1, 
************************************************
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
************************************************
*ELSET,ELSET=EALL
TEST,REFER
*NORMAL
TEST,LEFT,      -0.,    -1.,  0.
TEST,MID,     -.086, -.9963,  0.
TEST,RIGHT, -.11723, -.9931,  0.
*MATERIAL,NAME=ELAS
*ELASTIC
1000., .25
*EXPANSION
1.E-3,
*DENSITY
1.E-3,
*MATERIAL,NAME=ELAS2
*ELASTIC
1000., .25
*EXPANSION
1.E-3,
*DENSITY
1.E-5,
*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
NALL,3
LEFT,1,6
RIGHT,1,1,-1.0
RIGHT,2,2, 0.3
RIGHT,4,5, 0.0
RIGHT,6,6,.11750075
MID,1,1,-0.5
MID,2,2, 0.1
MID,4,5, 0.0
MID,6,6,.086105945
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*ELEMENT OUTPUT,ELSET=EOUT
E,S
*ELEMENT OUTPUT,ELSET=EOUT,REBAR
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
LEFT,1,6
*END STEP
**
**
*STEP,NLGEOM,INC=100
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.025,1.
*DLOAD
TEST, GRAV,150.,2.,2.,.1
REFER,GRAV,150.,2.,2.,.1
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.25,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM,INC=100
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.025,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY THERMAL LOAD
*STATIC
.25,1.0
*BOUNDARY,OP=NEW,FIXED
LEFT,1,6
*TEMPERATURE
NALL,100.
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
1,
S,THE,
E,TEMP
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=EOUT
S,THE
E,TEMP
*ELEMENT OUTPUT,ELSET=EOUT,REBAR
S,THE
E,TEMP
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
62, , ,
*BOUNDARY,OP=NEW
LEFT,1,6
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.0075,.075
*CLOAD
TIP,2,150.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*ELEMENT OUTPUT,ELSET=EOUT
E,S
*ELEMENT OUTPUT,ELSET=EOUT,REBAR
E,S
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
3.,30., 3
*CLOAD
TIP,2,150.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
3., 30., 3
*BOUNDARY,OP=NEW
LEFT,1,6
*CLOAD
TIP,2,150.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*output,field
*END STEP
**
**
*STEP,NLGEOM
 STEP #12: DYNAMICS
*DYNAMIC
.0075,.075
*BOUNDARY,OP=NEW
LEFT,1,6
*DLOAD,OP=NEW
*END STEP