*HEADING
 DATA DECK TO VERIFY BODY LOADS FOR DIFFERENT ELEMENTS WITH REBARS
**
** TEST OF REBAR: BODY LOADS WITH VARIOUS ELEMENTS WITH REBAR
**
** 1) CPS4R 
** 2) CPE4R
** 3) CAX4R
** 4) C3D8R
** 5) M3D4R
** 6) S4R
** 7) S4RS
** 8) S4RSW
** 9) SAX1
**
**  ANALYTIC RESULTS: 
**     For body load (BX)
**  NODE:       RF1: 
**   11       -2.5000
**   12       -2.5000
**   13       -2.5000
**   14       -2.5000
**   21       -2.5000
**   22       -2.5000
**   23       -2.5000
**   24       -2.5000
**   31       -52.360
**   32       -57.596
**   33       -57.596
**   34       -52.360
**   41       -1.2500
**   42       -1.2500
**   43       -1.2500
**   44       -1.2500
**   45       -1.2500
**   46       -1.2500
**   47       -1.2500
**   48       -1.2500
**   51       -.62500
**   52       -.62500
**   53       -.62500
**   54       -.62500
**   61       -.62500
**   62       -.62500
**   63       -.62500
**   64       -.62500
**   71       -.62500
**   72       -.62500
**   73       -.62500
**   74       -.62500
**   81       -.62500
**   82       -.62500
**   83       -.62500
**   84       -.62500
**   91       -49.742
**   92       -52.360
** 
**     For gravity load (GRAV)
**  NODE:       RF1:
**   11       -9.16667E-04
**   12       -9.16667E-04
**   13       -5.83333E-04
**   14       -5.83333E-04
**   21       -9.16667E-04
**   22       -9.16667E-04
**   23       -5.83333E-04
**   24       -5.83333E-04
**   31       -8.37758E-03
**   32       -8.90118E-03
**   33       -7.33038E-03
**   34       -6.80678E-03
**   41       -4.58333E-04
**   42       -4.58333E-04
**   43       -4.58333E-04
**   44       -4.58333E-04
**   45       -2.91667E-04
**   46       -2.91667E-04
**   47       -2.91667E-04
**   48       -2.91667E-04
**   51       -3.125E-04
**   52       -3.125E-04
**   53       -3.125E-04
**   54       -3.125E-04
**   61       -3.125E-04
**   62       -3.125E-04
**   63       -3.125E-04
**   64       -3.125E-04
**   71       -3.125E-04
**   72       -3.125E-04
**   73       -3.125E-04
**   74       -3.125E-04
**   81       -3.125E-04
**   82       -3.125E-04
**   83       -3.125E-04
**   84       -3.125E-04
**   91       -2.487E-02
**   92       -2.618E-02
**
** TEST CASE 1)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
**
*NODE,NSET=NALL
11,0.,0.
12,1.,0.
13,1.,1.
14,0.,1.
*ELEMENT,TYPE=CPS4R,ELSET=BASE1
11,11,12,13,14
*SOLID SECTION,ELSET=BASE1, MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB1
11, .5, .5, 0., .33333333, 1 
**
** TEST CASE 2)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPE4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
**
*NODE,NSET=NALL
21,1.5,0.
22,2.5,0.
23,2.5,1.
24,1.5,1.
*ELEMENT,TYPE=CPE4R,ELSET=BASE2
21,21,22,23,24
*SOLID SECTION,ELSET=BASE2, MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB2
21, .5, .5, 0., .33333333, 1 
**
** TEST CASE 3)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CAX4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
**
*NODE,NSET=NALL
31,3.,0.
32,4.,0.
33,4.,1.
34,3.,1.
*ELEMENT,TYPE=CAX4R,ELSET=BASE3
31,31,32,33,34
*SOLID SECTION,ELSET=BASE3, MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB3
31, .5, .5, 0., .33333333, 1, ,1.5
**
** TEST CASE 4)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: C3D8R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
**
*NODE,NSET=NALL
41,  4.5, 0.,  0.
42,  5.5, 0.,  0.
43,  5.5, 0., -1.
44,  4.5, 0., -1.
45,  4.5, 1.,  0.
46,  5.5, 1.,  0.
47,  5.5, 1., -1.
48,  4.5, 1., -1.
*ELEMENT,TYPE=C3D8R,ELSET=BASE4
41,41,42,43,44,45,46,47,48
*SOLID SECTION,ELSET=BASE4,MATERIAL=BASE
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB4
41, .5, .5, 90., .33333333, 1, 1
**
** TEST CASE 5)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDPLANE
** ELEMENT TYPE: M3D4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
**
*NODE,NSET=NALL
51,  0.0, 3.,  0.
52,  1.0, 3.,  0.
53,  1.0, 3., -1.
54,  0.0, 3., -1.
*ELEMENT,TYPE=M3D4R,ELSET=BASE5
51,51,52,53,54
*MEMBRANE SECTION,ELSET=BASE5,MATERIAL=BASE, POISSON=0.25
.25, 
*REBAR LAYER
REB5,.25,.5,,REBAR,0.,
**
**
** TEST CASE 6)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDPLANE
** ELEMENT TYPE: S4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
**
*NODE,NSET=NALL
61,  1.5, 3.,  0.
62,  2.5, 3.,  0.
63,  2.5, 3., -1.
64,  1.5, 3., -1.
*ELEMENT,TYPE=S4R,ELSET=BASE6
61,61,62,63,64
*SHELL SECTION,ELSET=BASE6,MATERIAL=BASE
.25,
*REBAR LAYER
REB6,.25,.5,,REBAR,0.,1 
**
**
** TEST CASE 7)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDPLANE
** ELEMENT TYPE: S4RS
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
**
*NODE,NSET=NALL
71,  3.0, 3.,  0.
72,  4.0, 3.,  0.
73,  4.0, 3., -1.
74,  3.0, 3., -1.
*ELEMENT,TYPE=S4RS,ELSET=BASE7
71,71,72,73,74
*SHELL SECTION,ELSET=BASE7,MATERIAL=BASE
.25, 
*REBAR LAYER
REB7,.25,.5,,REBAR,0.,1
**
**
** TEST CASE 8)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDPLANE
** ELEMENT TYPE: S4RSW
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
**
*NODE,NSET=NALL
81,  4.5, 3.,  0.
82,  5.5, 3.,  0.
83,  5.5, 3., -1.
84,  4.5, 3., -1.
*ELEMENT,TYPE=S4RSW,ELSET=BASE8
81,81,82,83,84
*SHELL SECTION,ELSET=BASE8,MATERIAL=BASE
.25, 
*REBAR LAYER
REB8,.25,.5,,REBAR,0.,1
**
**
** TEST CASE 9)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDPLANE
** ELEMENT TYPE: SAX1
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
**
*NODE,NSET=NALL
91,  6.0, 3.,  0.
92,  7.0, 3.,  0.
*ELEMENT,TYPE=SAX1,ELSET=BASE9
91,91,92
*SHELL SECTION,ELSET=BASE9,MATERIAL=BASE
0.25,
*REBAR LAYER,GEOMETRY=ANGULAR
REB9,.25,4.408,,REBAR,0.,1 
**
**
*MATERIAL,NAME=BASE
*ELASTIC
1.E4, 0.25
*DENSITY
1.E-4, 
*MATERIAL,NAME=REBAR
*ELASTIC
2.E4, 
*DENSITY
2.E-4, 
**
*RESTART,TIMEMARKS=YES,WRITE,NUMBER=1
**
*NSET, NSET=QA_TEST_ALLNODES
NALL, NALL, NALL, NALL
NALL, NALL, NALL, NALL
NALL, 
**
*STEP
*DYNAMIC,EXPLICIT
,.001
*DLOAD
BASE1,BX, 10.
BASE2,BX, 10.
BASE3,BR, 10.
BASE4,BX, 10.
BASE5,BX, 10.
BASE6,BX, 10.
BASE7,BX, 10.
BASE8,BX, 10.
BASE9,BR, 10.
*FILE OUTPUT,NUMBER INTERVAL=1,TIMEMARKS=YES
*NODE FILE
RF, 
*ENERGY FILE
***OUTPUT,HISTORY,TIME INTERVAL=1.0E-6
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLWK,ETOTAL,
**DT,
**
*BOUNDARY
NALL, 1,6
**
*OUTPUT, FIELD, NUMBER INTERVAL=1, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
RF, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.001
*END STEP
*STEP
*DYNAMIC,EXPLICIT
,.001
*DLOAD,OP=NEW
BASE1,GRAV, 10., 1.0, 0.0, 0.0
BASE2,GRAV, 10., 1.0, 0.0, 0.0
**BASE3,GRAV, 10., 1.0, 0.0, 0.0
BASE3,BR, 1.E-3
BASE4,GRAV, 10., 1.0, 0.0, 0.0
BASE5,GRAV, 10., 1.0, 0.0, 0.0
BASE6,GRAV, 10., 1.0, 0.0, 0.0
BASE7,GRAV, 10., 1.0, 0.0, 0.0
BASE8,GRAV, 10., 1.0, 0.0, 0.0
**BASE9,GRAV, 10., 1.0, 0.0, 0.0
BASE9,BR, 1.E-3
*FILE OUTPUT,NUMBER INTERVAL=1,TIMEMARKS=YES
*NODE FILE
RF, 
*ENERGY FILE
*OUTPUT, FIELD, NUMBER INTERVAL=1, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
RF, 
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.001
***************************
**REBAR OUTPUT AVAILABLE **
** for M3D4R ELEMENTS    **
***************************
***
*OUTPUT,FIELD,NUMBER INTERVAL=1
*ELEMENT OUTPUT,REBAR=REB5
LE,S
************
*OUTPUT,HISTORY,TIME INTERVAL=0.001
*ELEMENT OUTPUT,REBAR=REB5,ELSET=BASE5
S11,RBFOR,RBANG,RBROT
*END STEP