*HEADING
 SHELL ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*PREPRINT,MODEL=YES,HISTORY=NO
*RESTART,WRITE
*************************
** TEST 1: S4R
*************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=S4R,ELSET=TEST1
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*SHELL SECTION,TEMPERATURE,MATERIAL=ELAS,ELSET=TEST1
.1,
*************************
** TEST 2: S8R
*************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=S8R,ELSET=TEST2
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*SHELL SECTION,TEMPERATURE,MATERIAL=ELAS,ELSET=TEST2
.1,
********************************
** TEST 3: S4R WITH ORIENTATION
********************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=S4R,ELSET=TEST3
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST3
.1,
*ORIENTATION,NAME=ORIENT3,SYSTEM=C
**************************
** TEST 4: S4R WITH REBAR
**************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=S4R,ELSET=TEST4
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*SHELL SECTION,TEMPERATURE,MATERIAL=ELAS,ELSET=TEST4
.1,
*REBAR,ELEMENT=SHELL,NAME=REBAR4,MATERIAL=ELAS
TEST4,1.,.05, .05,1
*REBAR,ELEMENT=SHELL,NAME=REBAR4_2,MATERIAL=ELAS
TEST4,1.,.05,-.05,1
*****************************
** TEST 5: S8R WITH REBAR
*****************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=S8R,ELSET=TEST5
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*SHELL SECTION,TEMPERATURE,MATERIAL=ELAS,ELSET=TEST5
.1,
*REBAR,ELEMENT=SHELL,NAME=REBAR5,MATERIAL=ELAS
TEST5,1.,.05, .05,1
*REBAR,ELEMENT=SHELL,NAME=REBAR5_2,MATERIAL=ELAS
TEST5,1.,.05,-.05,1
*************************************
** TEST 6: S8R WITH NODAL THICKNESS
*************************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*NSET,NSET=L6
*NSET,NSET=M6
*NSET,NSET=R6
*NODAL THICKNESS
L6,.05
M6,.1
R6,.15
*ELEMENT,TYPE=S8R,ELSET=TEST6
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*SHELL SECTION,MATERIAL=ELAS,ELSET=TEST6,NODAL THICKNESS
**************************
** TEST 7: S3R
**************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=S3R,ELSET=TEST7
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*SHELL SECTION,TEMPERATURE,MATERIAL=ELAS,ELSET=TEST7
.1,
**************************
** TEST 8: STRI3
**************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=S3R,ELSET=TEST7
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*SHELL SECTION,TEMPERATURE,MATERIAL=ELAS,ELSET=TEST8
.1,
*************************
** TEST 9: S4
*************************
*NODE,NSET=NALL
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=MID
*NSET,NSET=RIGHT
*NSET,NSET=TIP
*ELEMENT,TYPE=S4,ELSET=TEST9
*ELSET,ELSET=REFER
*ELSET,ELSET=TEST
*ELSET,ELSET=EOUT
*SHELL SECTION,TEMPERATURE,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,
*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
S,E
*NODE FILE,FREQ=999
RF,
*EL FILE,ELSET=EOUT,POSITION=CENTROID,FREQ=999
E,S
*OUTPUT,HISTORY,VARIABLE=PRESELECTED,FREQ=999
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
U,
*ELEMENT OUTPUT
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 DLOAD TYPE P TO BOTH ELEMENTS
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST,P,.025
REFER,P,.025
*END STEP
**
**
*STEP,NLGEOM,INC=20
 STEP #5: APPLY DLOAD TYPE GRAV
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, GRAV,100.,2.,2.,.1
REFER,GRAV,100.,2.,2.,.1
*OUTPUT,HISTORY,VARIABLE=PRESELECTED,FREQ=999
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
U,RF
*ELEMENT OUTPUT
S,
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: APPLY DLOAD TYPE BX
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, BX, .25
REFER,BX, .25
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY DLOAD TYPE BY
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, BY, .25
REFER,BY, .25
*END STEP
**
**
*STEP,NLGEOM,INC=20
 STEP #8: APPLY DLOAD TYPE BZ
*STATIC
0.25,1.
*DLOAD,OP=NEW
TEST, BZ, .1
REFER,BZ, .1
*END STEP