*HEADING
 MEMBRANE ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*PREPRINT,MODEL=YES,HISTORY=NO
******************
** TEST 1: M3D4R
******************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D4R,ELSET=TEST1
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST1
2.,
*****************
** TEST 2: M3D4
*****************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D4,ELSET=TEST2
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST2
2.,
*****************
** TEST 3: M3D8
*****************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*NSET,NSET=BOUND
*ELEMENT,TYPE=M3D8,ELSET=TEST3
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST3
2.,
******************
** TEST 4: M3D8R
******************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D8R,ELSET=TEST4
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST4
2.,
*****************
** TEST 5: M3D9R
*****************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D9R,ELSET=TEST5
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST5
2.,
*****************************
** TEST 6: M3D4R WITH REBAR
*****************************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D4R,ELSET=TEST6
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*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
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D4,ELSET=TEST7
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*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
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D8,ELSET=TEST8
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST8
2.,
*REBAR LAYER
REB8,1.0,0.5,0,ELAS,30,1
******************
** TEST 9: M3D8R
******************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D8R,ELSET=TEST9
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST9
2.,
*REBAR LAYER
REB9,1.0,0.5,0,ELAS,30,1
*******************
** TEST 10: M3D9R
*******************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D9R,ELSET=TEST10
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*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
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D4R,ELSET=TEST11
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ORIENTATION,NAME=ORIENT11,SYSTEM=C
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST11,ORIENTATION=ORIENT11
2.,
***********************************
** TEST 12: M3D4 WITH ORIENTATION
***********************************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D4,ELSET=TEST12
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ORIENTATION,NAME=ORIENT12,SYSTEM=C
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST12,ORIENTATION=ORIENT12
2.,
*********************************************
** TEST 13: M3D8R WITH REBAR AND ORIENTATION
*********************************************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=M3D8R,ELSET=TEST13
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*ORIENTATION,NAME=ORIENT13,SYSTEM=C
*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
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*NSET,NSET=L14
*NSET,NSET=M14
*NSET,NSET=R14
*NODAL THICKNESS
L14,1.5
M14,2.0
R14,2.5
*ELEMENT,TYPE=M3D8R,ELSET=TEST14
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=TEST14,NODAL THICKNESS
*******************
** TEST 15: M3D9R
*******************
*NODE,NSET=NALL
*NSET,NSET=NALLREF,GEN
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*NSET,NSET=L15
*NSET,NSET=M15
*NSET,NSET=R15
*NODAL THICKNESS
L15,1.5
M15,2.0
R15,2.5
*ELEMENT,TYPE=M3D9R,ELSET=TEST15
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*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
1000., .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
*NODE PRINT,NSET=LEFT,FREQ=999
RF,
*EL PRINT,ELSET=EOUT,POSITION=CENTROID,FREQ=999
S,E
*EL PRINT,ELSET=EOUT,REBAR,FREQ=99
S,E,RBANG,RBROT
*EL FILE,ELSET=EOUT,REBAR,FREQ=99
S,E,RBANG,RBROT
*NODE FILE,NSET=LEFT,FREQ=999
RF,
*EL FILE,ELSET=EOUT,POSITION=CENTROID,FREQ=999
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
BOT,3
NALL,2
*END STEP
**
**
*STEP,NLGEOM,INC=100
 STEP #4: APPLY LOAD TO BOTH ELEMENTS
*STATIC
0.025,1.
*DLOAD,OP=NEW
TEST, GRAV,35000.,.1,0.,.1
REFER,GRAV,35000.,.1,0.,.1
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: APPLY DLOAD TYPE BX
*STATIC
.25,1.0
*DLOAD, OP=NEW
TEST,  BX, 150.
REFER, BX, 150.
TEST, BZ, 150.
REFER,BZ, 150.
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: APPLY DLOAD TYPE CENT
*STATIC
.25,1.0
*DLOAD, OP=NEW
 TEST1, CENT, 6.,0.,0.5,0.,0.,0.,1.
 TEST2, CENT, 6.,0.,0.5,0.,0.,0.,1.
 TEST3, CENT, 6.,0.,0.5,0.,0.,0.,1.
 TEST4, CENT, 6.,0.,0.5,0.,0.,0.,1.
 TEST5, CENT, 6.,0.,0.5,0.,0.,0.,1.
 TEST6, CENT, 6.,0.,5.5,0.,0.,0.,1.
 TEST7, CENT, 6.,0.,5.5,0.,0.,0.,1.
 TEST8, CENT, 6.,0.,5.5,0.,0.,0.,1.
 TEST9, CENT, 6.,0.,5.5,0.,0.,0.,1.
 TEST10, CENT, 6.,0.,5.5,0.,0.,0.,1.
 TEST11, CENT, 6.,0.,10.5,0.,0.,0.,1.
 TEST12, CENT, 6.,0.,10.5,0.,0.,0.,1.
 TEST13, CENT, 6.,0.,10.5,0.,0.,0.,1.
 TEST14, CENT, 6.,0.,10.5,0.,0.,0.,1.
 TEST15, CENT, 6.,0.,10.5,0.,0.,0.,1.
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY DLOAD TYPE CORIO
*STATIC
.25,1.0
*DLOAD, OP=NEW
 TEST1, CORIO, 1.E5,0.,0.5,0.,0.,0.,1.
 TEST2, CORIO, 1.E5,0.,0.5,0.,0.,0.,1.
 TEST3, CORIO, 1.E5,0.,0.5,0.,0.,0.,1.
 TEST4, CORIO, 1.E5,0.,0.5,0.,0.,0.,1.
 TEST5, CORIO, 1.E5,0.,0.5,0.,0.,0.,1.
 TEST6, CORIO, 1.E5,0.,5.5,0.,0.,0.,1.
 TEST7, CORIO, 1.E5,0.,5.5,0.,0.,0.,1.
 TEST8, CORIO, 1.E5,0.,5.5,0.,0.,0.,1.
 TEST9, CORIO, 1.E5,0.,5.5,0.,0.,0.,1.
 TEST10, CORIO, 1.E5,0.,5.5,0.,0.,0.,1.
 TEST11, CORIO, 1.E5,0.,10.5,0.,0.,0.,1.
 TEST12, CORIO, 1.E5,0.,10.5,0.,0.,0.,1.
 TEST13, CORIO, 1.E5,0.,10.5,0.,0.,0.,1.
 TEST14, CORIO, 1.E5,0.,10.5,0.,0.,0.,1.
 TEST15, CORIO, 1.E5,0.,10.5,0.,0.,0.,1.
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: APPLY DLOAD TYPE CENTRIF
*STATIC
.25,1.0
*DLOAD, OP=NEW
 TEST1, CENTRIF, 100.,0.,0.5,0.,0.,0.,1.
 TEST2, CENTRIF, 100.,0.,0.5,0.,0.,0.,1.
 TEST3, CENTRIF, 100.,0.,0.5,0.,0.,0.,1.
 TEST4, CENTRIF, 100.,0.,0.5,0.,0.,0.,1.
 TEST5, CENTRIF, 100.,0.,0.5,0.,0.,0.,1.
 TEST6, CENTRIF, 100.,0.,5.5,0.,0.,0.,1.
 TEST7, CENTRIF, 100.,0.,5.5,0.,0.,0.,1.
 TEST8, CENTRIF, 100.,0.,5.5,0.,0.,0.,1.
 TEST9, CENTRIF, 100.,0.,5.5,0.,0.,0.,1.
 TEST10, CENTRIF, 100.,0.,5.5,0.,0.,0.,1.
 TEST11, CENTRIF, 100.,0.,10.5,0.,0.,0.,1.
 TEST12, CENTRIF, 100.,0.,10.5,0.,0.,0.,1.
 TEST13, CENTRIF, 100.,0.,10.5,0.,0.,0.,1.
 TEST14, CENTRIF, 100.,0.,10.5,0.,0.,0.,1.
 TEST15, CENTRIF, 100.,0.,10.5,0.,0.,0.,1.
*END STEP
**
**
*STEP,NLGEOM
 STEP #9: APPLY DLOAD TYPE ROTA
*STATIC
.25,1.0
*DLOAD, OP=NEW
 TEST1, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST2, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST3, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST4, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST5, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST6, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST7, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST8, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST9, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST10, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST11, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST12, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST13, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST14, ROTA, 100.,0.,0.,0.,0.,-1.,0.
 TEST15, ROTA, 100.,0.,0.,0.,0.,-1.,0.
*END STEP