** TESTS OF REBAR: UNIAXIAL ALONG X-DIRECTION, FOLLOWED
**                 BY ROTATION AT END
**                 ELASTIC BASE
**                 ELASTIC REBAR
** STEP 1)
** LOADING: DISPLACEMENT ALONG X-AXIS 
** ANALYTIC RESULTS:
**      SOLID STRESS    REBAR STRESS   REBAR STRAIN
**        ~ 0.0           1906.2         .09531
**      REACTION FORCE AT LEFT EDGE NODES
**        Fx = 866.46
**
** STEP 2) 
** LOADING: ROTATION ABOUT Y-AXIS AT RIGHT END
** ANALYTIC RESULTS:
**            REBAR STRESS   REBAR STRAIN
** REBAR 1       1185.0         .05925
** REBAR 2       2627.4         .13136
**      REACTION FORCE AT LEFT EDGE NODES
**        Fx = 866.46
**        My = 74.5
** TEST WITH SHELL OFFSET
**            REBAR STRESS   REBAR STRAIN
** REBAR 1       1906.2         .09531
** REBAR 2       2627.4         .1650
**
*HEADING
**
** TEST CASE 1)
** REBAR, 0 DEGREES
** ELEMENT TYPE: S4R
** 
**
*NODE,NSET=NALL
11,  0., 0.
12,  1., 0.
13,  1., 1.
14,  0., 1.
*NSET,NSET=BOT
11,12
*NSET,NSET=LEFT
11,14
*NSET,NSET=RIGHT
12,13
*NSET,NSET=TOP
13,14
*ELEMENT,TYPE=S4R,ELSET=SHELL1
11, 11,12,13,14
*SHELL SECTION,ELSET=SHELL1,MATERIAL=BASE
1.0,
*REBAR LAYER
REBAR1,0.5,1.0,0.250,REBAR,0.0,1
REBAR1_2,0.5,1.0,-0.250,REBAR,0.0,1 
*ELSET,ELSET=EALL
SHELL1,
**
**
** TEST CASE 2)
** REBAR, 0 DEGREES
** ELEMENT TYPE: S4RS
**
*NODE,NSET=NALL
21,  2., 0.
22,  3., 0.
23,  3., 1.
24,  2., 1.
*NSET,NSET=BOT
21,22
*NSET,NSET=TOP
23,24
*NSET,NSET=LEFT
21,24
*NSET,NSET=RIGHT
22,23
*ELEMENT,TYPE=S4RS,ELSET=SHELL2
21, 21,22,23,24
*SHELL SECTION,ELSET=SHELL2,MATERIAL=BASE
1.0,
*REBAR LAYER
REBAR2,0.5,1.0,0.250,REBAR,0.0,1
REBAR2_2,0.5,1.0,-0.250,REBAR,0.0,1  
*ELSET,ELSET=EALL
SHELL2,
**
** TEST CASE 3)
** REBAR, 0 DEGREES
** ELEMENT TYPE: S4RSW
**
*NODE,NSET=NALL
31,  4., 0.
32,  5., 0.
33,  5., 1.
34,  4., 1.
*NSET,NSET=BOT
31,32
*NSET,NSET=TOP
33,34
*NSET,NSET=LEFT
31,34
*NSET,NSET=RIGHT
32,33
*ELEMENT,TYPE=S4RSW,ELSET=SHELL3
31, 32,33,34,31
*SHELL SECTION,ELSET=SHELL3,MATERIAL=BASE
1.0,
*REBAR LAYER
REBAR3,0.5,1.0,0.250,REBAR,0.0,1
REBAR3_2,0.5,1.0,-0.250,REBAR,0.0,1  
*ELSET,ELSET=EALL
SHELL3,
**
** TEST CASE 4)
** REBAR, 0 DEGREES
** ELEMENT TYPE: S4R, SHELL OFFSET
**
*NODE,NSET=NALL
41,  6., 0.
42,  7., 0.
43,  7., 1.
44,  6., 1.
*NSET,NSET=BOT
41,42
*NSET,NSET=TOP
43,44
*NSET,NSET=LEFT
41,44
*NSET,NSET=RIGHT
42,43
*ELEMENT,TYPE=S4R,ELSET=SHELL4
41, 42,43,44,41
*SHELL SECTION,ELSET=SHELL4,MATERIAL=BASE,OFFSET=.25
1.0,
*REBAR LAYER
REBAR4,0.5,1.0,0.250,REBAR,0.0,1
REBAR4_2,0.5,1.0,-0.250,REBAR,0.0,1  
*ELSET,ELSET=EALL
SHELL4,
**
************************************************
** TEST CASE 5)
** REBAR, 0 DEGREES
** ELEMENT TYPE: S3R
** 
**
*NODE,NSET=NALL
101,  0., 2.
102,  1., 2.
103,  1., 3.
104,  0., 3.
*NSET,NSET=BOT
101,102
*NSET,NSET=LEFT
101,104
*NSET,NSET=RIGHT
102,103
*NSET,NSET=TOP
103,104
*ELEMENT,TYPE=S3R,ELSET=SHELL5
101, 101,102,104
102, 102,103,104
*SHELL SECTION,ELSET=SHELL5,MATERIAL=BASE
1.0,
*REBAR LAYER
REBAR5,0.5,1.0,0.250,REBAR,0.0,1
REBAR5_2,0.5,1.0,-0.250,REBAR,0.0,1 
*ELSET,ELSET=EALL
SHELL5,
**
**
** TEST CASE 6)
** REBAR, 0 DEGREES
** ELEMENT TYPE: S3RS
**
*NODE,NSET=NALL
121,  2., 2.
122,  3., 2.
123,  3., 3.
124,  2., 3.
*NSET,NSET=BOT
121,122
*NSET,NSET=TOP
123,124
*NSET,NSET=LEFT
121,124
*NSET,NSET=RIGHT
122,123
*ELEMENT,TYPE=S3RS,ELSET=SHELL6
121, 121,122,124
122, 123,124,122
*SHELL SECTION,ELSET=SHELL6,MATERIAL=BASE
1.0,
*REBAR LAYER
REBAR6,0.5,1.0,0.250,REBAR,0.0,1
REBAR6_2,0.5,1.0,-0.250,REBAR,0.0,1  
*ELSET,ELSET=EALL
SHELL6,
**
** TEST CASE 7)
** REBAR, 0 DEGREES
** ELEMENT TYPE: S3R, SHELL OFFSET
**
*NODE,NSET=NALL
131,  4., 2.
132,  5., 2.
133,  5., 3.
134,  4., 3.
*NSET,NSET=BOT
131,132
*NSET,NSET=TOP
133,134
*NSET,NSET=LEFT
131,134
*NSET,NSET=RIGHT
132,133
*ELEMENT,TYPE=S3R,ELSET=SHELL7
131, 131,132,134
132, 132,133,134
*SHELL SECTION,ELSET=SHELL7,MATERIAL=BASE,OFFSET=.25
1.0,
*REBAR LAYER
REBAR7,0.5,1.0,0.250,REBAR,0.0,1
REBAR7_2,0.5,1.0,-0.250,REBAR,0.0,1  
*ELSET,ELSET=EALL
SHELL7,
**
***********************************************
**
** TEST CASE 8)
** REBAR, 0 DEGREES
** ELEMENT TYPE: SAX1
**
*NODE,NSET=NALL
51,  0., 4.
52,  1., 4.
*NSET,NSET=LEFT
51,
*NSET,NSET=RIGHT-AXI
52, 
*NSET,NSET=BOT
51,52
*ELEMENT,TYPE=SAX1,ELSET=SHELL8
51, 51,52
*SHELL SECTION,ELSET=SHELL8,MATERIAL=BASE
1.0,
*REBAR LAYER,GEOMETRY=ANGULAR
REBAR8,0.5,4.4084,0.250,REBAR,0.0,1
REBAR8_2,0.5,4.4084,-0.250,REBAR,0.0,1  
*ELSET,ELSET=EALL
SHELL8,
**
**
*MATERIAL,NAME=BASE
*ELASTIC
1., 0.498
*DENSITY
1.E-4, 
*MATERIAL,NAME=REBAR
*ELASTIC
2.E4, 
*DENSITY
1.E-4, 
******************
** LOAD AXIALLY **
******************
*AMPLITUDE,NAME=SMOOTH,DEFINITION=SMOOTH STEP
0.0,0.0,.02,1.0
**
*ELSET,ELSET=QA_TEST_ALLELEMS
SHELL1, SHELL2, SHELL3, SHELL4
SHELL5,SHELL6,SHELL7,SHELL8,
**
*NSET, NSET=QA_TEST_ALLNODES
NALL, BOT, LEFT, RIGHT
TOP, NALL, BOT, TOP
LEFT, RIGHT, NALL, BOT
TOP, LEFT, RIGHT, NALL
BOT, TOP, LEFT, RIGHT
NALL, LEFT, RIGHT-AXI, BOT
**
*STEP
*DYNAMIC,EXPLICIT
,.02
*BOUNDARY,AMPLITUDE=SMOOTH
NALL,3
LEFT,1
TOP,2
BOT,2
LEFT,4,6
RIGHT-AXI,4,6
RIGHT,4,6
RIGHT,1,1, 0.1
RIGHT-AXI,1,1, 0.1
**
*RESTART,WRITE,NUMBER INTERVAL=1
*FILE OUTPUT,TIMEMARKS=YES,NUMBER INTERVAL=1
*EL FILE
S,LE
STH, 
*NODE FILE
RF,U
**
*OUTPUT, FIELD, TIMEMARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST_ALLELEMS
S,LE,
STH, 
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U,
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.02
************************
*OUTPUT, FIELD, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST_ALLELEMS,REBAR
S,LE,RBFOR
*END STEP
*********************
** LOAD IN BENDING **
*********************
*STEP
*DYNAMIC,EXPLICIT
,.02
*BOUNDARY,AMPLITUDE=SMOOTH
RIGHT,5,5, -.17453
RIGHT-AXI,6,6, -.17453
**
*END STEP