*HEADING
 MEMBRANE ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*PREPRINT,MODEL=YES,HISTORY=NO
******************
** TEST 1: M3D4R
******************
*NODE,NSET=NALL
 1,  0.0, 1., 0.0
 2,  1.0, 1., 0.0
 3,  2.0, 1., 0.0
 4,  0.0, 1., 0.5
 5,  1.0, 1., 0.5
 6,  2.0, 1., 0.5
 7,  0.0, 1., 1.0
 8,  1.0, 1., 1.0
 9,  2.0, 1., 1.0
 11, 0.0, 0., 0.0
 12, 1.5, -.1, 0.0
 13, 3.0, -.3, 0.0
 14, 0.0,  .0, 0.5
 15, 1.5, -.1, 0.5
 16, 3.0, -.3, 0.5
 17, 0.0,  .0, 1.0
 18, 1.5, -.1, 1.0
 19, 3.0, -.3, 1.0
*NSET,NSET=NALLREF,GEN
 1,9
*NSET,NSET=BOT
 1,2,3,11,12,13
*NSET,NSET=LEFT
 1,11,4,14,7,17
*NSET,NSET=MID
 12,15,18
*NSET,NSET=RIGHT
 13,16,19
*NSET,NSET=TIP
 3,13,6,16,9,19
*ELEMENT,TYPE=M3D4R,ELSET=TEST1
 1, 1,2,5,4
 2, 2,3,6,5
 3, 4,5,8,7
 4, 5,6,9,8
 11,11,12,15,14
 12,12,13,16,15
 13,14,15,18,17
 14,15,16,19,18
*ELSET,ELSET=REFER
 1,2,3,4
*ELSET,ELSET=TEST
 11,12,13,14
*ELSET,ELSET=EOUT
 1,11
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST1
 2., 
*****************
** TEST 2: M3D4
*****************
*NODE,NSET=NALL
 101,  5.0, 1., 0.0
 102,  6.0, 1., 0.0
 103,  7.0, 1., 0.0
 104,  5.0, 1., 0.5
 105,  6.0, 1., 0.5
 106,  7.0, 1., 0.5
 107,  5.0, 1., 1.0
 108,  6.0, 1., 1.0
 109,  7.0, 1., 1.0
 111, 5.0,  .0, 0.0
 112, 6.5, -.1, 0.0
 113, 8.0, -.3, 0.0
 114, 5.0,  .0, 0.5
 115, 6.5, -.1, 0.5
 116, 8.0, -.3, 0.5
 117, 5.0,  .0, 1.0
 118, 6.5, -.1, 1.0
 119, 8.0, -.3, 1.0
*NSET,NSET=NALLREF,GEN
 101,109
*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,113,106,116,109,119
*ELEMENT,TYPE=M3D4,ELSET=TEST2
 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
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST2
 2., 
*****************
** TEST 3: M3D8
*****************
*NODE,NSET=NALL
 201, 10.0, 1., 0.0
 202, 11.0, 1., 0.0
 203, 12.0, 1., 0.0
 204, 10.0, 1., 0.5
 205, 11.0, 1., 0.5
 206, 12.0, 1., 0.5 
 207, 10.0, 1., 1.0
 208, 11.0, 1., 1.0
 209, 12.0, 1., 1.0
 211, 10.0,  .0, 0.0
 212, 11.5, -.1, 0.0
 213, 13.0, -.3, 0.0
 214, 10.0,  .0, 0.5
 215, 11.5, -.1, 0.5
 216, 13.0, -.3, 0.5 
 217, 10.0,  .0, 1.0
 218, 11.5, -.1, 1.0
 219, 13.0, -.3, 1.0
*NSET,NSET=NALLREF,GEN
 201,209
*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
*NSET,NSET=BOUND
 202,203,205,206,208,209,212,213,216,218,219
*ELEMENT,TYPE=M3D8,ELSET=TEST3
 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
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST3
 2., 
******************
** TEST 4: M3D8R
******************
*NODE,NSET=NALL
 301, 15.0, 1., 0.0
 302, 16.0, 1., 0.0
 303, 17.0, 1., 0.0
 304, 15.0, 1., 0.5
 305, 16.0, 1., 0.5
 306, 17.0, 1., 0.5 
 307, 15.0, 1., 1.0
 308, 16.0, 1., 1.0
 309, 17.0, 1., 1.0
 311, 15.0,  .0, 0.0
 312, 16.5, -.1, 0.0
 313, 18.0, -.3, 0.0
 314, 15.0,  .0, 0.5
 315, 16.5, -.1, 0.5
 316, 18.0, -.3, 0.5 
 317, 15.0,  .0, 1.0
 318, 16.5, -.1, 1.0
 319, 18.0, -.3, 1.0
*NSET,NSET=NALLREF,GEN
 301,309
*NSET,NSET=BOT
 301,302,303,311,312,313
*NSET,NSET=LEFT
 301,304,307,311,314,317
*NSET,NSET=MID
 312,318
*NSET,NSET=RIGHT
 313,316,319
*NSET,NSET=TIP
 303,306,309,313,316,319
*ELEMENT,TYPE=M3D8R,ELSET=TEST4
 301,301,303,309,307,302,306,308,304
 311,311,313,319,317,312,316,318,314
*ELSET,ELSET=REFER
 301, 
*ELSET,ELSET=TEST
 311, 
*ELSET,ELSET=EOUT
 301,311
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST4
 2., 
*****************
** TEST 5: M3D9R
*****************
*NODE,NSET=NALL
 401, 20.0, 1., 0.0
 402, 21.0, 1., 0.0
 403, 22.0, 1., 0.0
 404, 20.0, 1., 0.5
 405, 21.0, 1., 0.5
 406, 22.0, 1., 0.5 
 407, 20.0, 1., 1.0
 408, 21.0, 1., 1.0
 409, 22.0, 1., 1.0
 411, 20.0,  .0, 0.0
 412, 21.5, -.1, 0.0
 413, 23.0, -.3, 0.0
 414, 20.0,  .0, 0.5
 415, 21.5, -.1, 0.5
 416, 23.0, -.3, 0.5 
 417, 20.0,  .0, 1.0
 418, 21.5, -.1, 1.0
 419, 23.0, -.3, 1.0
*NSET,NSET=NALLREF,GEN
 401,409
*NSET,NSET=BOT
 401,402,403,411,412,413
*NSET,NSET=LEFT
 401,404,407,411,414,417
*NSET,NSET=MID
 412,415,418
*NSET,NSET=RIGHT
 413,416,419
*NSET,NSET=TIP
 403,406,409,413,416,419
*ELEMENT,TYPE=M3D9R,ELSET=TEST5
 401,401,403,409,407,402,406,408,404,405
 411,411,413,419,417,412,416,418,414,415
*ELSET,ELSET=REFER
 401, 
*ELSET,ELSET=TEST
 411, 
*ELSET,ELSET=EOUT
 401,411
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST5
 2., 
*****************************
** TEST 6: M3D4R WITH REBAR
*****************************
*NODE,NSET=NALL
 501,  0.0, 6., 0.0
 502,  1.0, 6., 0.0
 503,  2.0, 6., 0.0
 504,  0.0, 6., 0.5
 505,  1.0, 6., 0.5
 506,  2.0, 6., 0.5
 507,  0.0, 6., 1.0
 508,  1.0, 6., 1.0
 509,  2.0, 6., 1.0
 511, 0.0,  5.0, 0.0
 512, 1.5,  4.9, 0.0
 513, 3.0,  4.7, 0.0
 514, 0.0,  5.0, 0.5
 515, 1.5,  4.9, 0.5
 516, 3.0,  4.7, 0.5
 517, 0.0,  5.0, 1.0
 518, 1.5,  4.9, 1.0
 519, 3.0,  4.7, 1.0
*NSET,NSET=NALLREF,GEN
 501,509
*NSET,NSET=BOT
 501,502,503,511,512,513
*NSET,NSET=LEFT
 501,511,504,514,507,517
*NSET,NSET=MID
 512,515,518
*NSET,NSET=RIGHT
 513,516,519
*NSET,NSET=TIP
 503,513,506,516,509,519
*ELEMENT,TYPE=M3D4R,ELSET=TEST6
 501, 501,502,505,504
 502, 502,503,506,505
 503, 504,505,508,507
 504, 505,506,509,508
 511,511,512,515,514
 512,512,513,516,515
 513,514,515,518,517
 514,515,516,519,518
*ELSET,ELSET=REFER
 501,502,503,504
*ELSET,ELSET=TEST
 511,512,513,514
*ELSET,ELSET=EOUT
 501,511
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST6
 2., 
*REBAR LAYER
REB6,1.0,0.5,0,ELAS,30,1 
****************************
** TEST 7: M3D4 WITH REBAR
****************************
*NODE,NSET=NALL
 601,  5.0, 6., 0.0
 602,  6.0, 6., 0.0
 603,  7.0, 6., 0.0
 604,  5.0, 6., 0.5
 605,  6.0, 6., 0.5
 606,  7.0, 6., 0.5
 607,  5.0, 6., 1.0
 608,  6.0, 6., 1.0
 609,  7.0, 6., 1.0
 611, 5.0,  5.0, 0.0
 612, 6.5,  4.9, 0.0
 613, 8.0,  4.7, 0.0
 614, 5.0,  5.0, 0.5
 615, 6.5,  4.9, 0.5
 616, 8.0,  4.7, 0.5
 617, 5.0,  5.0, 1.0
 618, 6.5,  4.9, 1.0
 619, 8.0,  4.7, 1.0
*NSET,NSET=NALLREF,GEN
 601,609
*NSET,NSET=BOT
 601,602,603,611,612,613
*NSET,NSET=LEFT
 601,611,604,614,607,617
*NSET,NSET=MID
 612,615,618
*NSET,NSET=RIGHT
 613,616,619
*NSET,NSET=TIP
 603,613,606,616,609,619
*ELEMENT,TYPE=M3D4,ELSET=TEST7
 601, 601,602,605,604
 602, 602,603,606,605
 603, 604,605,608,607
 604, 605,606,609,608
 611,611,612,615,614
 612,612,613,616,615
 613,614,615,618,617
 614,615,616,619,618
*ELSET,ELSET=REFER
 601,602,603,604
*ELSET,ELSET=TEST
 611,612,613,614
*ELSET,ELSET=EOUT
 601,611
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST7
 2., 
*REBAR LAYER
REB7,1.0,0.5,0,ELAS,30,1 
****************************
** TEST 8: M3D8 WITH REBAR
****************************
*NODE,NSET=NALL
 701, 10.0, 6., 0.0
 702, 11.0, 6., 0.0
 703, 12.0, 6., 0.0
 704, 10.0, 6., 0.5
 705, 11.0, 6., 0.5
 706, 12.0, 6., 0.5 
 707, 10.0, 6., 1.0
 708, 11.0, 6., 1.0
 709, 12.0, 6., 1.0
 711, 10.0,  5.0, 0.0
 712, 11.5,  4.9, 0.0
 713, 13.0,  4.7, 0.0
 714, 10.0,  5.0, 0.5
 715, 11.5,  4.9, 0.5
 716, 13.0,  4.7, 0.5 
 717, 10.0,  5.0, 1.0
 718, 11.5,  4.9, 1.0
 719, 13.0,  4.7, 1.0
*NSET,NSET=NALLREF,GEN
 701,709
*NSET,NSET=BOT
 701,702,703,711,712,713
*NSET,NSET=LEFT
 701,704,707,711,714,717
*NSET,NSET=MID
 712,718
*NSET,NSET=RIGHT
 713,716,719
*NSET,NSET=TIP
 703,706,709,713,716,719
*ELEMENT,TYPE=M3D8,ELSET=TEST8
 701,701,703,709,707,702,706,708,704
 711,711,713,719,717,712,716,718,714
*ELSET,ELSET=REFER
 701, 
*ELSET,ELSET=TEST
 711, 
*ELSET,ELSET=EOUT
 701,711
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST8
 2.,
*REBAR LAYER
REB8,1.0,0.5,0,ELAS,30,1  
******************
** TEST 9: M3D8R
******************
*NODE,NSET=NALL
 801, 15.0, 6., 0.0
 802, 16.0, 6., 0.0
 803, 17.0, 6., 0.0
 804, 15.0, 6., 0.5
 805, 16.0, 6., 0.5
 806, 17.0, 6., 0.5 
 807, 15.0, 6., 1.0
 808, 16.0, 6., 1.0
 809, 17.0, 6., 1.0
 811, 15.0,  5.0, 0.0
 812, 16.5,  4.9, 0.0
 813, 18.0,  4.7, 0.0
 814, 15.0,  5.0, 0.5
 815, 16.5,  4.9, 0.5
 816, 18.0,  4.7, 0.5 
 817, 15.0,  5.0, 1.0
 818, 16.5,  4.9, 1.0
 819, 18.0,  4.7, 1.0
*NSET,NSET=NALLREF,GEN
 801,809
*NSET,NSET=BOT
 801,802,803,811,812,813
*NSET,NSET=LEFT
 801,804,807,811,814,817
*NSET,NSET=MID
 812,818
*NSET,NSET=RIGHT
 813,816,819
*NSET,NSET=TIP
 803,806,809,813,816,819
*ELEMENT,TYPE=M3D8R,ELSET=TEST9
 801,801,803,809,807,802,806,808,804
 811,811,813,819,817,812,816,818,814
*ELSET,ELSET=REFER
 801, 
*ELSET,ELSET=TEST
 811, 
*ELSET,ELSET=EOUT
 801,811
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST9
 2.,
*REBAR LAYER
REB9,1.0,0.5,0,ELAS,30,1  
*******************
** TEST 10: M3D9R
*******************
*NODE,NSET=NALL
 901, 20.0, 6., 0.0
 902, 21.0, 6., 0.0
 903, 22.0, 6., 0.0
 904, 20.0, 6., 0.5
 905, 21.0, 6., 0.5
 906, 22.0, 6., 0.5 
 907, 20.0, 6., 1.0
 908, 21.0, 6., 1.0
 909, 22.0, 6., 1.0
 911, 20.0,  5.0, 0.0
 912, 21.5,  4.9, 0.0
 913, 23.0,  4.7, 0.0
 914, 20.0,  5.0, 0.5
 915, 21.5,  4.9, 0.5
 916, 23.0,  4.7, 0.5 
 917, 20.0,  5.0, 1.0
 918, 21.5,  4.9, 1.0
 919, 23.0,  4.7, 1.0
*NSET,NSET=NALLREF,GEN
 901,909
*NSET,NSET=BOT
 901,902,903,911,912,913
*NSET,NSET=LEFT
 901,904,907,911,914,917
*NSET,NSET=MID
 912,915,918
*NSET,NSET=RIGHT
 913,916,919
*NSET,NSET=TIP
 903,906,909,913,916,919
*ELEMENT,TYPE=M3D9R,ELSET=TEST10
 901,901,903,909,907,902,906,908,904,905
 911,911,913,919,917,912,916,918,914,915
*ELSET,ELSET=REFER
 901, 
*ELSET,ELSET=TEST
 911, 
*ELSET,ELSET=EOUT
 901,911
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST10
 2., 
*REBAR LAYER
REB10,1.0,0.5,0,ELAS,30,1 
************************************
** TEST 11: M3D4R WITH ORIENTATION
************************************
*NODE,NSET=NALL
 1501,  0.0,  11., 0.0
 1502,  1.0,  11., 0.0
 1503,  2.0,  11., 0.0
 1504,  0.0,  11., 0.5
 1505,  1.0,  11., 0.5
 1506,  2.0,  11., 0.5
 1507,  0.0,  11., 1.0
 1508,  1.0,  11., 1.0
 1509,  2.0,  11., 1.0
 1511, 0.0,  10.0, 0.0
 1512, 1.5,   9.9, 0.0
 1513, 3.0,   9.7, 0.0
 1514, 0.0,  10.0, 0.5
 1515, 1.5,   9.9, 0.5
 1516, 3.0,   9.7, 0.5
 1517, 0.0,  10.0, 1.0
 1518, 1.5,   9.9, 1.0
 1519, 3.0,   9.7, 1.0
*NSET,NSET=NALLREF,GEN
 1501,1509
*NSET,NSET=BOT
 1501,1502,1503,1511,1512,1513
*NSET,NSET=LEFT
 1501,1511,1504,1514,1507,1517
*NSET,NSET=MID
 1512,1515,1518
*NSET,NSET=RIGHT
 1513,1516,1519
*NSET,NSET=TIP
 1503,1513,1506,1516,1509,1519
*ELEMENT,TYPE=M3D4R,ELSET=TEST11
 1501, 1501,1502,1505,1504
 1502, 1502,1503,1506,1505
 1503, 1504,1505,1508,1507
 1504, 1505,1506,1509,1508
 1511, 1511,1512,1515,1514
 1512, 1512,1513,1516,1515
 1513, 1514,1515,1518,1517
 1514, 1515,1516,1519,1518
*ELSET,ELSET=REFER
 1501,1502,1503,1504
*ELSET,ELSET=TEST
 1511,1512,1513,1514
*ELSET,ELSET=EOUT
 1501,1511
*ORIENTATION,NAME=ORIENT11,SYSTEM=C
 0.,1.,0.,  0.,2.,0.
 3,0.
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST11,ORIENTATION=ORIENT11
 2., 
***********************************
** TEST 12: M3D4 WITH ORIENTATION
***********************************
*NODE,NSET=NALL
 1601,  5.0,  11., 0.0
 1602,  6.0,  11., 0.0
 1603,  7.0,  11., 0.0
 1604,  5.0,  11., 0.5
 1605,  6.0,  11., 0.5
 1606,  7.0,  11., 0.5
 1607,  5.0,  11., 1.0
 1608,  6.0,  11., 1.0
 1609,  7.0,  11., 1.0
 1611, 5.0,  10.0, 0.0
 1612, 6.5,   9.9, 0.0
 1613, 8.0,   9.7, 0.0
 1614, 5.0,  10.0, 0.5
 1615, 6.5,   9.9, 0.5
 1616, 8.0,   9.7, 0.5
 1617, 5.0,  10.0, 1.0
 1618, 6.5,   9.9, 1.0
 1619, 8.0,   9.7, 1.0
*NSET,NSET=NALLREF,GEN
 1601,1609
*NSET,NSET=BOT
 1601,1602,1603,1611,1612,1613
*NSET,NSET=LEFT
 1601,1611,1604,1614,1607,1617
*NSET,NSET=MID
 1612,1615,1618
*NSET,NSET=RIGHT
 1613,1616,1619
*NSET,NSET=TIP
 1603,1613,1606,1616,1609,1619
*ELEMENT,TYPE=M3D4,ELSET=TEST12
 1601, 1601,1602,1605,1604
 1602, 1602,1603,1606,1605
 1603, 1604,1605,1608,1607
 1604, 1605,1606,1609,1608
 1611, 1611,1612,1615,1614
 1612, 1612,1613,1616,1615
 1613, 1614,1615,1618,1617
 1614, 1615,1616,1619,1618
*ELSET,ELSET=REFER
 1601,1602,1603,1604
*ELSET,ELSET=TEST
 1611,1612,1613,1614
*ELSET,ELSET=EOUT
 1601,1611
*ORIENTATION,NAME=ORIENT12,SYSTEM=C
 5.,1.,0.,  5.,2.,0.
 3,0.
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST12,ORIENTATION=ORIENT12
 2., 
*********************************************
** TEST 13: M3D8R WITH REBAR AND ORIENTATION
*********************************************
*NODE,NSET=NALL
 1701, 10.0,  11., 0.0
 1702, 11.0,  11., 0.0
 1703, 12.0,  11., 0.0
 1704, 10.0,  11., 0.5
 1705, 11.0,  11., 0.5
 1706, 12.0,  11., 0.5 
 1707, 10.0,  11., 1.0
 1708, 11.0,  11., 1.0
 1709, 12.0,  11., 1.0
 1711, 10.0, 10.0, 0.0
 1712, 11.5,  9.9, 0.0
 1713, 13.0,  9.7, 0.0
 1714, 10.0, 10.0, 0.5
 1715, 11.5,  9.9, 0.5
 1716, 13.0,  9.7, 0.5 
 1717, 10.0, 10.0, 1.0
 1718, 11.5,  9.9, 1.0
 1719, 13.0,  9.7, 1.0
*NSET,NSET=NALLREF,GEN
 1701,1709
*NSET,NSET=BOT
 1701,1702,1703,1711,1712,1713
*NSET,NSET=LEFT
 1701,1704,1707,1711,1714,1717
*NSET,NSET=MID
 1712,1718
*NSET,NSET=RIGHT
 1713,1716,1719
*NSET,NSET=TIP
 1703,1706,1709,1713,1716,1719
*ELEMENT,TYPE=M3D8R,ELSET=TEST13
 1701,1701,1703,1709,1707,1702,1706,1708,1704
 1711,1711,1713,1719,1717,1712,1716,1718,1714
*ELSET,ELSET=REFER
 1701, 
*ELSET,ELSET=TEST
 1711, 
*ELSET,ELSET=EOUT
 1701,1711
*ORIENTATION,NAME=ORIENT13,SYSTEM=C
 10.,1.,0.,  10.,2.,0.
 3,0.
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST13,ORIENTATION=ORIENT13
 2., 
*REBAR LAYER
REB13,1.0,0.5,0,ELAS,30,1 
****************************************
** TEST 14: M3D8R WITH NODAL THICKNESS
****************************************
*NODE,NSET=NALL
 1801, 15.0,  11., 0.0
 1802, 16.0,  11., 0.0
 1803, 17.0,  11., 0.0
 1804, 15.0,  11., 0.5
 1805, 16.0,  11., 0.5
 1806, 17.0,  11., 0.5 
 1807, 15.0,  11., 1.0
 1808, 16.0,  11., 1.0
 1809, 17.0,  11., 1.0
 1811, 15.0, 10.0, 0.0
 1812, 16.5,  9.9, 0.0
 1813, 18.0,  9.7, 0.0
 1814, 15.0, 10.0, 0.5
 1815, 16.5,  9.9, 0.5
 1816, 18.0,  9.7, 0.5 
 1817, 15.0, 10.0, 1.0
 1818, 16.5,  9.9, 1.0
 1819, 18.0,  9.7, 1.0
*NSET,NSET=NALLREF,GEN
 1801,1809
*NSET,NSET=BOT
 1801,1802,1803,1811,1812,1813
*NSET,NSET=LEFT
 1801,1804,1807,1811,1814,1817
*NSET,NSET=MID
 1812,1818
*NSET,NSET=RIGHT
 1813,1816,1819
*NSET,NSET=TIP
 1803,1806,1809,1813,1816,1819
*NSET,NSET=L14
 1801,1804,1807,1811,1814,1817
*NSET,NSET=M14
 1802,1808,1812,1818
*NSET,NSET=R14
 1803,1806,1809,1813,1816,1819
*NODAL THICKNESS
 L14,1.5
 M14,2.0
 R14,2.5
*ELEMENT,TYPE=M3D8R,ELSET=TEST14
 1801,1801,1803,1809,1807,1802,1806,1808,1804
 1811,1811,1813,1819,1817,1812,1816,1818,1814
*ELSET,ELSET=REFER
 1801, 
*ELSET,ELSET=TEST
 1811, 
*ELSET,ELSET=EOUT
 1801,1811
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST14,NODAL THICKNESS
*******************
** TEST 15: M3D9R
*******************
*NODE,NSET=NALL
 1901, 20.0,  11., 0.0
 1902, 21.0,  11., 0.0
 1903, 22.0,  11., 0.0
 1904, 20.0,  11., 0.5
 1905, 21.0,  11., 0.5
 1906, 22.0,  11., 0.5 
 1907, 20.0,  11., 1.0
 1908, 21.0,  11., 1.0
 1909, 22.0,  11., 1.0
 1911, 20.0, 10.0, 0.0
 1912, 21.5,  9.9, 0.0
 1913, 23.0,  9.7, 0.0
 1914, 20.0, 10.0, 0.5
 1915, 21.5,  9.9, 0.5
 1916, 23.0,  9.7, 0.5 
 1917, 20.0, 10.0, 1.0
 1918, 21.5,  9.9, 1.0
 1919, 23.0,  9.7, 1.0
*NSET,NSET=NALLREF,GEN
 1901,1909
*NSET,NSET=BOT
 1901,1902,1903,1911,1912,1913
*NSET,NSET=LEFT
 1901,1904,1907,1911,1914,1917
*NSET,NSET=MID
 1912,1915,1918
*NSET,NSET=RIGHT
 1913,1916,1919
*NSET,NSET=TIP
 1903,1906,1909,1913,1916,1919
*NSET,NSET=L15
 1901,1904,1907,1911,1914,1917
*NSET,NSET=M15
 1902,1905,1908,1912,1915,1918
*NSET,NSET=R15
 1903,1906,1909,1913,1916,1919
*NODAL THICKNESS
 L15,1.5
 M15,2.0
 R15,2.5
*ELEMENT,TYPE=M3D9R,ELSET=TEST15
 1901,1901,1903,1909,1907,1902,1906,1908,1904,1905
 1911,1911,1913,1919,1917,1912,1916,1918,1914,1915
*ELSET,ELSET=REFER
 1901, 
*ELSET,ELSET=TEST
 1911, 
*ELSET,ELSET=EOUT
 1901,1911
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST15,NODAL THICKNESS
************************************************
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
************************************************
*ELSET,ELSET=EALL
 TEST,REFER
*MATERIAL,NAME=ELAS
*ELASTIC
 1.E3, .25
*EXPANSION
 1.E-3, 
*DENSITY
 1.E-3, 
*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
 NALLREF,2
 LEFT,1,2
 RIGHT,1,1,-1.0
 RIGHT,2,2, 0.3
 MID,1,1,-0.5
 MID,2,2, 0.1
*EL FILE,ELSET=EOUT,REBAR,FREQ=99
 S,E,RBANG,RBROT
*EL FILE,ELSET=EOUT,FREQ=999
 E, S
*NODE FILE,NSET=LEFT,FREQ=999
 RF, 
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT,NSET=LEFT
 RF, 
*ELEMENT OUTPUT,ELSET=EOUT
 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,FIXED
 LEFT,1
 NALL,2
 BOT,3
*END STEP
**
**
*STEP,NLGEOM,INC=100
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
 0.025,1.
*DLOAD
 TEST, GRAV,35000.,1.,0.,.1
 REFER,GRAV,35000.,1.,0.,.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
 NALL,2
 BOT,3
*TEMPERATURE
 NALL,100.
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos
 60, , ,
*BOUNDARY,OP=NEW
 LEFT,1
 NALL,2
 BOT,3
*EL PRINT,ELSET=EOUT,FREQ=0
*NODE PRINT,NSET=LEFT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*MODAL FILE 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
 .00025,.0025
*CLOAD
 TIP,1,150.
*EL FILE,ELSET=EOUT,REBAR,FREQ=99
 S,E,RBANG,RBROT
*EL FILE,ELSET=EOUT,FREQ=999
 E, S
*NODE FILE,NSET=LEFT,FREQ=999
 RF, 
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT,NSET=LEFT
 RF, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=EOUT
 E,S
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
 145.,150., 5
*CLOAD
 TIP,1,150.
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
 145., 150., 5
*BOUNDARY,OP=NEW
 LEFT,1
 NALL,2
 BOT,3
*CLOAD
 TIP,1,150.
*END STEP
**
**
*STEP,PERTURBATION
 STEP #12: STEADY-STATE DYNAMICS,SUBSPACE PROJECTION 
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
 145., 150., 5
*CLOAD
 TIP,1,150.
*END STEP
**
**
*STEP,NLGEOM
 STEP #13: DYNAMICS
*DYNAMIC
.00025,.0025
*CLOAD,AMP=RAMP
 TIP,1,150.
*output,field
*elementoutput
*elementoutput,rebar
*END STEP